Dynamics and stability of relativistic GRB blast waves

In gamma-ray-bursts (GRB), ultra-relativistic blast waves are ejected into the circumburst medium. We analyse in unprecedented detail the deceleration of a self-similar Blandford-McKee blast wave from a Lorentz factor 25 to the nonrelativistic Sedov phase. Our goal is to determine the stability properties of its frontal shock. We carried out a grid-adaptive relativistic 2D hydro-simulation at extreme resolving power, following the GRB jet during the entire afterglow phase. We investigate the effect of the finite initial jet opening angle on the deceleration of the blast wave, and identify the growth of various instabilities throughout the coasting shock front. We find that during the relativistic phase, the blast wave is subject to pressure-ram pressure instabilities that ripple and fragment the frontal shock. These instabilities manifest themselves in the ultra-relativistic phase alone, remain in full agreement with causality arguments, and decay slowly to finally disappear in the near-Newtonian phase as the shell Lorentz factor drops below 3. From then on, the compression rate decreases to levels predicted to be stable by a linear analysis of the Sedov phase. Our simulations confirm previous findings that the shell also spreads laterally because a rarefaction wave slowly propagates to the jet axis, inducing a clear shell deformation from its initial spherical shape. The blast front becomes meridionally stratified, with decreasing speed from axis to jet edge.Comment: accepted for publication in A&

Transverse stability of relativistic two-component jets

Context: Astrophysical jets from various sources seem to be stratified, with a fast inner jet and a slower outer jet. As it is likely that the launching mechanism for each component is different, their interface will develop differential rotation, while the outer jet radius represents a second interface where disruptions may occur. Aims: We explore the stability of stratified, rotating, relativistic two-component jets, in turn embedded in static interstellar medium. Methods: In a grid-adaptive relativistic hydrodynamic simulation with the AMRVAC code, the non-linear azimuthal stability of two-component relativistic jets is investigated. We simulate until multiple inner jet rotations have been completed. Results: We find evidence for the development of an extended shear flow layer between the two jet components, resulting from the growth of a body mode in the inner jet, Kelvin-Helmholtz surface modes at their original interface, and their nonlinear interaction. Both wave modes are excited by acoustic waves which are reflected between the symmetry axis and the interface of the two jet components. Their interaction induces the growth of near stationary, counterrotating vortices at the outer edge of the shear flow layer. The presence of a heavy external jet allows to slow down their further development, and maintain a collimated flow. At the outer jet boundary, small-scale Rayleigh-Taylor instabilities develop, without disrupting the jet configuration. Conclusion: We demonstrate that the cross-section of two-component relativistic jets, with a heavy, cold outer jet, is non-linearly stable.Comment: Accepted in A&A 24/09/200

AMRVAC and Relativistic Hydrodynamic simulations for GRB afterglow phases

We apply a novel adaptive mesh refinement code, AMRVAC, to numerically investigate the various evolutionary phases in the interaction of a relativistic shell with its surrounding cold Interstellar Medium (ISM). We do this for both 1D isotropic as well as full 2D jetlike fireball models. This is relevant for Gamma Ray Bursts, and we demonstrate that, thanks to the AMR strategy, we resolve the internal structure of the shocked shell-ISM matter, which will leave its imprint on the GRB afterglow. We determine the deceleration from an initial Lorentz factor $\gamma=100$ up to the almost Newtonian $\gamma\sim{\cal O}(2)$ phase of the flow. We present axisymmetric 2D shell evolutions, with the 2D extent characterized by their initial opening angle. In such jetlike GRB models, we discuss the differences with the 1D isotropic GRB equivalents. These are mainly due to thermally induced sideways expansions of both the shocked shell and shocked ISM regions. We found that the propagating 2D ultrarelativistic shell does not accrete all the surrounding medium located within its initial opening angle. Part of this ISM matter gets pushed away laterally and forms a wide bow-shock configuration with swirling flow patterns trailing the thin shell. The resulting shell deceleration is quite different from that found in isotropic GRB models. As long as the lateral shell expansion is merely due to ballistic spreading of the shell, isotropic and 2D models agree perfectly. As thermally induced expansions eventually lead to significantly higher lateral speeds, the 2D shell interacts with comparably more ISM matter and decelerates earlier than its isotropic counterpart.Comment: 12 pages, accepted in MNRAS, 12/01/200

Synchrotron radiation of self-collimating relativistic MHD jets

The goal of this paper is to derive signatures of synchrotron radiation from state-of-the-art simulation models of collimating relativistic magnetohydrodynamic (MHD) jets featuring a large-scale helical magnetic field. We perform axisymmetric special relativistic MHD simulations of the jet acceleration region using the PLUTO code. The computational domain extends from the slow magnetosonic launching surface of the disk up to 6000^2 Schwarzschild radii allowing to reach highly relativistic Lorentz factors. The Poynting dominated disk wind develops into a jet with Lorentz factors of 8 and is collimated to 1 degree. In addition to the disk jet, we evolve a thermally driven spine jet, emanating from a hypothetical black hole corona. Solving the linearly polarized synchrotron radiation transport within the jet, we derive VLBI radio and (sub-) mm diagnostics such as core shift, polarization structure, intensity maps, spectra and Faraday rotation measure (RM), directly from the Stokes parameters. We also investigate depolarization and the detectability of a lambda^2-law RM depending on beam resolution and observing frequency. We find non-monotonic intrinsic RM profiles which could be detected at a resolution of 100 Schwarzschild radii. In our collimating jet geometry, the strict bi-modality in polarization direction (as predicted by Pariev et al.) can be circumvented. Due to relativistic aberration, asymmetries in the polarization vectors across the jet can hint to the spin direction of the central engine.Comment: Submitted to Ap

Can the magnetic field in the Orion arm inhibit the growth of instabilities in the bow shock of Betelgeuse?

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$-Orionis 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.Comment: 5 pages, including 7 figures. The published version will include 4 animations. Accepted for publication in A&

Computing the dust distribution in the bowshock of a fast moving, evolved star

We study the hydrodynamical behavior occurring in the turbulent interaction zone of a fast moving red supergiant star, where the circumstellar and interstellar material collide. In this wind-interstellar medium collision, the familiar bow shock, contact discontinuity, and wind termination shock morphology forms, with localized instability development. Our model includes a detailed treatment of dust grains in the stellar wind, and takes into account the drag forces between dust and gas. The dust is treated as pressureless gas components binned per grainsize, for which we use ten representative grainsize bins. Our simulations allow to deduce how dust grains of varying sizes become distributed throughout the circumstellar medium. We show that smaller dust grains (radius <0.045 micro-meters) tend to be strongly bound to the gas and therefore follow the gas density distribution closely, with intricate finestructure due to essentially hydrodynamical instabilities at the wind-related contact discontinuity. Larger grains which are more resistant to drag forces are shown to have their own unique dust distribution, with progressive deviations from the gas morphology. Specifically, small dust grains stay entirely within the zone bound by shocked wind material. The large grains are capable of leaving the shocked wind layer, and can penetrate into the shocked or even unshocked interstellar medium. Depending on how the number of dust grains varies with grainsize, this should leave a clear imprint in infrared observations of bowshocks of red supergiants and other evolved stars.Comment: Accepted for publication in ApJL, 4 figure

Using numerical models of bow shocks to investigate the circumstellar medium of massive stars

Many massive stars travel through the interstellar medium at supersonic speeds. As a result they form bow shocks at the interface between the stellar wind. We use numerical hydrodynamics to reproduce such bow shocks numerically, creating models that can be compared to observations. In this paper we discuss the influence of two physical phenomena, interstellar magnetic fields and the presence of interstellar dust grains on the observable shape of the bow shocks of massive stars. We find that the interstellar magnetic field, though too weak to restrict the general shape of the bow shock, reduces the size of the instabilities that would otherwise be observed in the bow shock of a red supergiant. The interstellar dust grains, due to their inertia can penetrate deep into the bow shock structure of a main sequence O-supergiant, crossing over from the ISM into the stellar wind. Therefore, the dust distribution may not always reflect the morphology of the gas. This is an important consideration for infrared observations, which are dominated by dust emission. Our models clearly show, that the bow shocks of massive stars are useful diagnostic tools that can used to investigate the properties of both the stellar wind as well as the interstellar medium.Comment: 7 pages, 4 figures, to be published in the Journal of Physics: Conference Series (JPCS) as part of the proceedings of the 13th Annual International Astrophysics Conferenc

No visible optical variability from a relativistic blast wave encountering a wind-termination shock

Gamma-ray burst afterglow flares and rebrightenings of the optical and X-ray light curve have been attributed to both late time inner engine activity and density changes in the medium surrounding the burster. To test the latter, we study the encounter between the relativistic blast wave from a gamma-ray burster and a stellar wind termination shock. The blast wave is simulated using a high performance adaptive mesh relativistic hydrodynamics code, AMRVAC, and the synchrotron emission is analyzed in detail with a separate radiation code. We find no bump in the resulting light curve, not even for very high density jumps. Furthermore, by analyzing the contributions from the different shock wave regions we are able to establish that it is essential to resolve the blast wave structure in order to make qualitatively correct predictions on the observed output and that the contribution from the reverse shock region will not stand out, even when the magnetic field is increased in this region by repeated shocks. This study resolves a controversy in recent literature.Comment: 4 figures, submitted to MNRAS letter

Effect of angular opening on the dynamics of relativistic hydro jets

Context. Relativistic jets emerging from AGN cores transfer energy from the core to their surrounding ISM/IGM. Because jets are observed to have finite opening angles, one needs to quantify the role of conical versus cylindrical jet propagation in this energy transfer. Aims. We use FR-II AGN jets parameter with finite opening angles. We study the effect of the variation of the opening angle on the dynamics and energy transfer of the jet. We also point out how the characteristics of this external medium, such as its density profile, play a role in the dynamics. Methods. This study exploits our parallel AMR code MPI-AMRVAC with its special relativistic hydrodynamic model, incorporating an equation of state with varying effective polytropic index. We studied mildly under-dense jets up to opening angles of 10 degrees, at Lorentz factors of about 10, inspired by observations. Instantaneous quantification of the various ISM volumes and their energy content allows one to quantify the role of mixing versus shock-heated cocoon regions over the time intervals. Results. We show that a wider opening angle jet results in a faster deceleration of the jet and leads to a wider cocoon dominated by Kelvin-Helmholtz and Rayleigh-Taylor instabilities. The energy transfer mainly occurs in the shocked ISM region by both the frontal bow shock and cocoon-traversing shock waves, in a roughly 3 to 1 ratio to the energy transfer of the mixing zone, for a 5 degree opening angle jet. A rarefaction wave induces a dynamically formed layered structure of the jet beam. Conclusions. Finite opening angle jets can efficiently transfer significant fractions (25 % up to 70 %) of their injected energy over a growing region of shocked ISM matter. The role of the ISM stratification is prominent for determining the overall volume that is affected by relativistic jet injection

Le choix de l’accompagnement dans la petite entreprise en démarrage (The small enterprise in launching phase: choice of accompanimen)

L’intérêt du dirigeant de l’entreprise en phase de démarrage est de recourir à des structures d’accompagnements afin de combler les insuffisances techniques, managériales et organisationnelles. La décision du choix de l’accompagnateur se fonde sur des critères de perceptions et de crédibilité que le dirigeant peut discerner à travers certains signes. The interest of a company’s leader in the launching phase is to resort to structures of accompaniment in order to fill technical, managerial and organizational inefficiencies. The choice of the coach is based on perception and credibility criteria which the leader can distinguish through certain signs.small enterprise, launching phase, PME