Long gamma ray bursts (LGRBs) are associated to the collapse of a massive star and the formation of a relativistic jet. As the jet propagates through the star, it forms an extended, hot cocoon. The dynamical evolution of the jet/cocoon system and its
interaction with the environment has been studied extensively both analytically and numerically. On the other hand, the role played by the supernova (SN) explosion associated with LGRBs in determining the outcome of the system has been barely considered. In this paper, we discuss the large landscape of outcomes resulting from the interaction of the SN, jet and cocoon. We show that the outcome depends mainly on three timescales: the times for the cocoon and supernova shock wave to break through the surface of the progenitor star, and the time needed for the cocoon to engulf completely the progenitor star. The delay between the launch of the SN shock moving through the progenitor star and the jet can be related to these three timescales. Depending on the ordering of these time scales, the jet-cocoon might propagate inside the SN ejecta or the other way around, and the outcome for the properties of the explosion would be different. We discuss the imprint of the complex interaction between the jet-cocoon and the supernova shock on the emergent thermal and non-thermal radiation.
The discovery of GRB 170817A, the first unambiguous off-axis short gamma-ray burst arising from a neutron star merger, has challenged our understanding of the angular structure of relativistic jets. Studies of the jet propagation usually assume that
the jet is ejected from the central engine with a top-hat structure and its final structure, which determines the observed light curve and spectra, is primarily regulated by the interaction with the nearby environment. However, jets are expected to be produced with a structure that is more complex than a simple top-hat, as shown by global accretion simulations. We present numerical simulations of short GRBs launched with a wide range of initial structures, durations and luminosities. We follow the jet interaction with the merger remnant wind and compute its final structure at distances $gtrsim 10^{11}$~cm from the central engine. We show that the final jet structure, as well as the resulting afterglow emission, depend strongly on the initial structure of the jet, its luminosity and duration. While the initial structure at the jet is preserved for long-lasting SGRBs, it is strongly modified for jets barely making their way through the wind. This illustrates the importance of combining the results of global simulations with propagation studies in order to better predict the expected afterglow signatures from neutron star mergers. Structured jets provide a reasonable description of the GRB 170817A afterglow emission with an off-axis angle $theta_{rm obs} approx 22.5^circ$.
The first neutron star binary merger detected in gravitational waves, GW170817 and the subsequent detection of its emission across the electromagnetic spectrum showed that these systems are viable progenitors of short $gamma$-ray bursts (sGRB). The a
fterglow signal of GW170817 has been found to be consistent with a structured GRB jet seen off-axis, requiring significant amounts of relativistic material at large angles. This trait can be attributed to the interaction of the relativistic jet with the external wind medium. Here we perform numerical simulations of relativistic jets interacting with realistic wind environments in order to explore how the properties of the wind and central engine affect the structure of successful jets. We find that the angular energy distribution of the jet depends primarily on the ratio between the lifetime of the jet and the time it takes the merger remnant to collapse. We make use of these simulations to constrain the time it took for the merger remnant in GW170817 to collapse into a black hole based on the angular structure of the jet as inferred from afterglow observations. We conclude that the lifetime of the merger remnant in GW170817 was $approx 1-1.7$s, which, after collapse, triggered the formation of the jet.
The discovery of jets from tidal disruption events (TDEs) rejuvenated the old field of relativistic jets powered by accretion onto supermassive black holes. In this Chapter, we first review the extensive multi-wavelength observations of jetted TDEs.
Then, we show that these events provide valuable information on many aspects of jet physics from a new prospective, including the on-and-off switch of jet launching, jet propagation through the ambient medium, $gamma/$X-ray radiation mechanism, jet composition, and the multi-messenger picture. Finally, open questions and future prospects in this field are summarized.
The early optical emission of gamma-ray bursts gives an opportunity to understand the central engine and first stages of these events. About 30% of GRBs present flares whose origin is still a subject of discussion. We present optical photometry of GR
B 180620A with the COATLI telescope and RATIR instrument. COATLI started to observe from the end of prompt emission at $T+39.3$~s and RATIR from $T+121.4$~s. We supplement the optical data with the X-ray light curve from emph{Swift}/XRT. %The optical and X-ray light curves show very unusual behavior with features clearly beyond the standard fireball model. We observe an optical flare from $T+110$ to $T+550$~s, with a temporal index decay $alpha_mathrm{O,decay}=1.32pm 0.01$, and a $Delta t/t=1.63$, which we interpret as the signature of a reverse shock component. After the initial normal decay the light curves show a long plateau from $T+500$ to $T+7800$~s both in X-rays and the optical before decaying again after an achromatic jet break at $T+7800$~s. Fluctuations are seen during the plateau phase in the optical. Adding to the complexity of GRB afterglows, the plateau phase (typically associated with the coasting phase of the jet) is seen in this object after the ``normal decay phase (emitted during the deceleration phase of the jet) and the jet break phase occurs directly after the plateau. We suggest that this sequence of events can be explained by a rapid deceleration of the jet with $t_dlesssim 40$ s due to the high density of the environment ($approx 100$ cm$^{-3}$) followed by reactivation of the central engine which causes the flare and powers the plateau phase.
Relativistic supernovae constitute a sub-class of type Ic supernovae (SNe). Their non-thermal, radio emission differs notably from that of regular type Ic supernovae as they have a fast expansion speed (with velocities $sim$ 0.6-0.8 c) which can not
be explained by a standard, spherical SN explosion but advocates for a quickly evolving, mildly relativistic ejecta associated with the SN. In this paper, we compute the synchrotron radiation emitted by the cocoon of a long gamma-ray burst jet (GRB). We show that the energy and velocity of the expanding cocoon, and the radio non-thermal light curves and spectra are consistent with those observed in relativistic SNe. Thus, the radio emission from this events is not coming from the SN shock front, but from the mildly relativistic cocoon produced by the passage of a GRB jet through the progenitor star. We also show that the cocoon radio emission dominates the GRB emission at early times for GRBs seen off-axis, and the flux can be larger at late times compared with on-axis GRBs if the cocoon energy is at least comparable with respect to the GRB energy.
We present hydrodynamic simulations of the hot cocoon produced when a relativistic jet passes through the gamma-ray burst (GRB) progenitor star and its environment, and we compute the lightcurve and spectrum of the radiation emitted by the cocoon. Th
e radiation from the cocoon has a nearly thermal spectrum with a peak in the X-ray band, and it lasts for a few minutes in the observer frame; the cocoon radiation starts at roughly the same time as when $gamma$-rays from a burst trigger detectors aboard GRB satellites. The isotropic cocoon luminosity ($sim 10^{47}$ erg s$^{-1}$) is of the same order of magnitude as the X-ray luminosity of a typical long-GRB afterglow during the plateau phase. This radiation should be identifiable in the Swift data because of its nearly thermal spectrum which is distinct from the somewhat brighter power-law component. The detection of this thermal component would provide information regarding the size and density stratification of the GRB progenitor star. Photons from the cocoon are also inverse-Compton (IC) scattered by electrons in the relativistic jet. We present the IC lightcurve and spectrum, by post-processing the results of the numerical simulations. The IC spectrum lies in 10 keV--MeV band for typical GRB parameters. The detection of this IC component would provide an independent measurement of GRB jet Lorentz factor and it would also help to determine the jet magnetisation parameter.
Observations of several protostellar jets show systematic differences in radial velocity transverse to the jet propagation direction, which have been interpreted as evidence of rotation in the jets. In this paper we discuss the origin of these veloci
ty shifts, and show that they could be originated by rotation in the flow, or by side to side asymmetries in the shock velocity, which could be due to asymmetries in the jet ejection velocity/density or in the ambient medium. For typical poloidal jet velocities (~ 100-200 km/s), an asymmetry >~ 10% can produce velocity shifts comparable to those observed. We also present three dimensional numerical simulations of rotating, precessing and asymmetric jets, and show that, even though for a given jet there is a clear degeneracy between these effects, a statistical analysis of jets with different inclination angles can help to distinguish between the alternative origins of transverse velocity shifts. Our analysis indicate that side to side velocities asymmetries could represent an important contribution to transverse velocity shifts, being the most important contributor for large jet inclination angles (with respect the the plane of the sky), and can not be neglected when interpreting the observations.
We present 3D, adaptive mesh refinement simulations of G2, a cloud of gas moving in a highly eccentric orbit towards the galactic center. We assume that G2 originates from a stellar wind interacting with the environment of the Sgr A* black hole. The
stellar wind forms a cometary bubble which becomes increasingly elongated as the star approaches periastron. A few months after periastron passage, streams of material begin to accrete on the central black hole with accretion rates $dot{M} sim 10^{-8}$ M$_odot$ yr$^{-1}$. Predicted Br$gamma$ emission maps and position-velocity diagrams show an elongated emission resembling recent observations of G2. A large increase in luminosity is predicted by the emission coming from the shocked wind region during periastron passage. The observations, showing a constant Br$gamma$ luminosity, remain puzzling, and are explained here assuming that the emission is dominated by the free-wind region. The observed Br$gamma$ luminosity ($sim 8 times 10^{30}$ erg s$^{-1}$) is reproduced by a model with a $v_w=50$ km s$^{-1}$ wind velocity and a $10^{-7}$ M$_odot$ yr$^{-1}$ mass loss rate if the emission comes from the shocked wind. A faster and less dense wind reproduces the Br$gamma$ luminosity if the emission comes from the inner, free wind region. The extended cometary wind bubble, largely destroyed by the tidal interaction with the black hole, reforms a few years after periastron passage. As a result, the Br$gamma$ emission is more compact after periastron passage.
The circumstellar medium (CSM) of a massive star is modified by its winds before a supernova (SN) explosion occurs, and thus the evolution of the resulting supernova remnant (SNR) is influenced by both the geometry of the explosion as well as the com
plex structure of the CSM. Motivated by recent work suggesting the SNR W49B was a jet-driven SN expanding in a complex CSM, we explore how the dynamics and the metal distributions in a jet-driven explosion are modified by the interaction with the surrounding environment. In particular, we perform hydrodynamical calculations to study the dynamics and explosive nucleosynthesis of a jet-driven SN triggered by the collapse of a 25 solar masses Wolf-Rayet star and its subsequent interaction with the CSM up to several hundred years following the explosion. We find that although the CSM has small-scale effects on the structure of the SNR, the overall morphology and abundance patterns are reflective of the initial asymmetry of the SN explosion. Thus, we predict that jet-driven SNRs, such as W49B, should be identifiable based on morphology and abundance patterns at ages up to several hundred years, even if they expand into a complex CSM environment.
