Using the relativistic MHD code MPI-AMRVAC and a radiative transfer code in post-processing, we explore the influence of the magnetic-field configuration and transverse stratification of an over-pressured jet on its morphology, on the moving shock dy
namics, and on the emitted radio light curve. First, we investigate different large-scale magnetic fields with their effects on the standing shocks and on the stratified jet morphology. Secondly, we study the interaction of a moving shock wave with the standing shocks. We calculate the synthetic synchrotron maps and radio light curves and analyse the variability at two frequencies 1 and 15.3 GHz and for several observation angles. Finally, we compare the characteristics of our simulated light curves with radio flares observed from the blazar 3C 273 with OVRO and VLBA in the MOJAVE survey between 2008 and 2019. We find that, in a structured, over-pressured relativistic jet, the presence of the large-scale magnetic field structure changes the properties of the standing shock waves and leads to an opening of the jet. When crossing such standing shocks, moving shock waves accompanying overdensities injected in the base of the jet are causing very luminous radio flares. The observation of the temporal structure of these flares under different viewing angles probes the jet at different optical depths. At 1 GHz and for small angles, the self-absorption caused by the moving shock wave becomes more important and leads to a drop in the observed flux after it interacts with the brightest standing knot. A weak asymmetry is seen in the shape of the simulated flares, resulting from the remnant emission of the shocked standing shocks. The characteristics of the simulated flares and the correlation of peaks in the light curve with the crossing of moving and standing shocks favor this scenario as an explanation of the observed radio flares of 3C 273.
Several observations of astrophysical jets show evidence of a structure in the direction perpendicular to the jet axis, leading to the development of spine & sheath models of jets. Most studies focus on a two-component jet consisting of a highly rela
tivistic inner jet and a slower - but still relativistic - outer jet surrounded by an unmagnetized environment. These jets are believed to be susceptible to a relativistic Rayleigh-Taylor-type instability, depending on the effective inertia ratio of the two components. We extend previous studies by taking into account the presence of a non-zero toroidal magnetic field. Different values of magnetization are examined, to detect possible differences in the evolution and stability of the jet. We find that the toroidal field, above a certain level of magnetization $sigma$, roughly equal to 0.01, can stabilize the jet against the previously mentioned instabilities and that there is a clear trend in the behaviour of the average Lorentz factor and the effective radius of the jet when we continuously increase the magnetization. The simulations are performed using the relativistic MHD module from the open source, parallel, grid adaptive, MPI-AMRVAC code.
The winds of massive stars create large (>10 pc) bubbles around their progenitors. As these bubbles expand they encounter the interstellar coherent magnetic field which, depending on its strength, can influence the shape of the bubble. We wish to inv
estigate if, and how much, the interstellar magnetic field can contribute to the shape of an expanding circumstellar bubble around a massive star. We use the MPI-AMRVAC code to make magneto-hydrodynamical simulations of bubbles, using a single star model, combined with several different field strengths: B = 5, 10, and 20 muG for the interstellar magnetic field. This covers the typical field strengths of the interstellar magnetic fields found in the galactic disk and bulge. Furthermore, we present two simulations that include both a 5 muG interstellar magnetic field and a 10,000 K interstellar medium and two different ISM densities to demonstrate how the magnetic field can combine with other external factors to influence the morphology of the circumstellar bubbles. Our results show that low magnetic fields, as found in the galactic disk, inhibit the growth of the circumstellar bubbles in the direction perpendicular to the field. As a result, the bubbles become ovoid, rather than spherical. Strong interstellar fields, such as observed for the galactic bulge, can completely stop the expansion of the bubble in the direction perpendicular to the field, leading to the formation of a tube-like bubble. When combined with a warm, high-density ISM the bubble is greatly reduced in size, causing a dramatic change in the evolution of temporary features inside the bubble. The magnetic field of the interstellar medium can affect the shape of circumstellar bubbles. This effect may have consequences for the shape and evolution of circumstellar nebulae and supernova remnants, which are formed within the main wind-blown bubble.
Many evolved stars travel through space at supersonic velocities, which leads to the formation of bow shocks ahead of the star where the stellar wind collides with the interstellar medium (ISM). Herschel observations of the bow shock of $alpha$-Orion
is show that the shock is almost free of instabilities, despite being, at least in theory, subject to both Kelvin-Helmholtz and Rayleigh-Taylor instabilities. A possible explanation for the lack of instabilities lies in the presence of an interstellar magnetic field. We wish to investigate whether the magnetic field of the interstellar medium (ISM) in the Orion arm can inhibit the growth of instabilities in the bow shock of $alpha$-Orionis. We used the code MPI-AMRVAC to make magneto-hydrodynamic simulations of a circumstellar bow shock, using the wind parameters derived for $alpha$-Orionis and interstellar magnetic field strengths of $B,=,1.4,, 3.0$, and $5.0, mu$G, which fall within the boundaries of the observed magnetic field strength in the Orion arm of the Milky Way. Our results show that even a relatively weak magnetic field in the interstellar medium can suppress the growth of Rayleigh-Taylor and Kelvin-Helmholtz instabilities, which occur along the contact discontinuity between the shocked wind and the shocked ISM. The presence of even a weak magnetic field in the ISM effectively inhibits the growth of instabilities in the bow shock. This may explain the absence of such instabilities in the Herschel observations of $alpha$-Orionis.
We present analytic flux prescriptions for broadband spectra of self-absorbed and optically thin synchrotron radiation from gamma-ray burst afterglows, based on one-dimensional relativistic hydrodynamic simulations. By treating the evolution of criti
cal spectrum parameters as a power-law break between the ultrarelativistic and non-relativistic asymptotic solutions, we generalize the prescriptions to any observer time. Our aim is to provide a set of formulas that constitutes a useful tool for accurate fitting of model-parameters to observational data, regardless of the dynamical phase of the outflow. The applicability range is not confined to gamma-ray burst afterglows, but includes all spherical outflows (also jets before the jet-break) that produce synchrotron radiation as they adiabatically decelerate in a cold, power-law medium. We test the accuracy of the prescriptions and show that numerical evidence suggests that typical relative errors in the derivation of physical quantities are about 10 per cent. A software implementation of the presented flux prescriptions combined with a fitting code is freely available on request and on-line. Together they can be used in order to directly fit model parameters to data.
Numerical models of the wind-blown bubble of massive stars usually only account for the wind of a single star. However, since massive stars are usually formed in clusters, it would be more realistic to follow the evolution of a bubble created by seve
ral stars. We develope a two-dimensional (2D) model of the circumstellar bubble created by two massive stars, a 40 solar mass star and a 25 solar mass star, and follow its evolution. The stars are separated by approximately 16 pc and surrounded by a cold medium with a density of 20 particles per cubic cm. We use the MPI-AMRVAC hydrodynamics code to solve the conservation equations of hydrodynamics on a 2D cylindrical grid using time-dependent models for the wind parameters of the two stars. At the end of the stellar evolution (4.5 and 7.0 million years for the 40 and 25 solar mass stars, respectively), we simulate the supernova explosion of each star. Each star initially creates its own bubble. However, as the bubbles expand they merge, creating a combined, aspherical bubble. The combined bubble evolves over time, influenced by the stellar winds and supernova explosions. The evolution of a wind-blown bubble created by two stars deviates from that of the bubbles around single stars. In particular, once one of the stars has exploded, the bubble is too large for the wind of the remaining star to maintain and the outer shell starts to disintegrate. The lack of thermal pressure inside the bubble also changes the behavior of circumstellar features close to the remaining star. The supernovae are contained inside the bubble, which reflects part of the energy back into the circumstellar medium.
We present the first-ever simulations of non-ideal magnetohydrodynamical (MHD) stellar winds coupled with disc-driven jets where the resistive and viscous accretion disc is self-consistently described. The transmagnetosonic, collimated MHD outflows a
re investigated numerically using the VAC code. Our simulations show that the inner outflow is accelerated from the central object hot corona thanks to both the thermal pressure and the Lorentz force. In our framework, the thermal acceleration is sustained by the heating produced by the dissipated magnetic energy due to the turbulence. Conversely, the outflow launched from the resistive accretion disc is mainly accelerated by the magneto-centrifugal force. We also show that when a dense inner stellar wind occurs, the resulting disc-driven jet have a different structure, namely a magnetic structure where poloidal magnetic field lines are more inclined because of the pressure caused by the stellar wind. This modification leads to both an enhanced mass ejection rate in the disc-driven jet and a larger radial extension which is in better agreement with the observations besides being more consistent.
