After its launch, a GRB jet propagates through a dense medium prior to its breakout. The interaction of the jet with the medium gives rise to the formation of a complex structured outflow, often referred to as a "structured jet". This structure is essential for our understanding of GRBs as it ultimately dictates their emission signatures. However, to date, the underlying physics which sets the post-breakout jet morphology remained unexplored, and its modeling has been mainly done by assuming ad-hoc functions. Here, using a set of 3D simulations we follow the evolution of hydrodynamic long and short GRB jets after breakout and provide a physically motivated post-breakout outflow structure of GRB jets (the structure of weakly magnetized jets is underway). Our simulations feature Rayleigh-Taylor fingers, which grow from the cocoon into the jet and destabilize it. The mixing of jet-cocoon material gives rise to a previously unidentified region between the two, which we denote the jet-cocoon interface (JCI). In lGRBs the mixing is strong, leading to most of the outflow energy to drift into the JCI. In sGRBs, where the medium is lighter, the mixing is weaker, and the JCI and the jet core hold a comparable amount of energy. Remarkably, the jet structure (core and JCI) in all systems can be characterized by simple angular power-law distributions of power and velocity, with power-law indices that depend solely on the level of mixing. This result supports the commonly used power-law angular distribution, and disfavors a Gaussian jet modeling. At larger angles, where the cocoon dominates, the structure is more complex including both an angular and a radial structure. The mixing shapes the prompt light curve and implies that typical afterglows of lGRBs are different from those of sGRBs. The predictions that we provide can be used to infer jet characteristics from prompt and afterglow observations.

Strong variability is a common characteristic of the prompt emission of gamma-ray bursts (GRB). This observed variability is widely attributed to an intermittency of the central engine, through formation of strong internal shocks in the GRB-emitting jet expelled by the engine. In this paper we study numerically the propagation of hydrodynamic jets, injected periodically by a variable engine, through the envelope of a collapsed star. By post-processing the output of 3D numerical simulations, we compute the net radiative efficiency of the outflow. We find that all intermittent jets are subject to heavy baryon contamination that inhibits the emission at and above the photosphere well below detection limits. This is in contrast to continuous jets that, as shown recently, produce a highly variable gamma-ray photospheric emission with high efficiency, owing to the interaction of the jet with the stellar envelope. Our results challenge the variable engine model for hydrodynamic jets, and may impose constraints on the duty cycle of GRB engines. If such systems exist in nature, they are not expected to produce bright gamma-ray emission, but should appear as X-ray, optical and radio transients that resemble a delayed GRB afterglow signal.

The first hours following a neutron star merger are considered to provide several UV/optical/NIR signals: $ \beta $-decay emission from free neutrons, radioactive decay of shocked heavy elements in the cocoon and cocoon's cooling emission. Here we consider two additional emission sources: $ \beta $-decay of free neutrons in the cocoon and synchrotron by the $ \beta $-decay electrons. We present 3D RHD simulations of jets that propagate in a multi-layer ejecta from the merger and calculate semi-analytically the resulting light curves. We find that the free neutrons emission at high latitudes is enhanced by the cocoon by a factor of a few to power a wide ($ \lesssim 60^\circ $) and brief ($ \sim 1 $ hour) UV signal that can reach an absolute magnitude of $\gtrsim$ -15, comparable with the cooling emission. If the ejected neutron matter mass is $ M_n \gtrsim 10^{-4}{\rm M_\odot} $, the synchrotron emission may yield a long ($ \sim 8 $ hours) quasi-isotropic UV/optical signal with an absolute magnitude between -12 and -15, depending on the magnetic field. Such a high mass of a mildly-relativistic component may partly obscure the cocoon's shocked r-process elements, thereby attenuating its radioactive decay emission. Future observations on these timescales, including null detections, may place constraints on the ejected neutron matter mass and shed light on the ejecta and jet-cocoon characteristics.

VLBI and JVLA observations revealed that GW170817 involved a narrow jet ($ \theta_j \approx 4^\circ $) that dominated the afterglow peak at our viewing angle, $ \theta_{\rm obs} \approx 20^\circ $. This implies that at the time of the afterglow peak, the observed signal behaved like an afterglow of a top-hat jet seen at $ \theta_{\rm obs} \gg \theta_j $, and it can be modeled by analytic expressions that describe such jets. We use a set of numerical simulations to calibrate these analytic relations and obtain generic equations for the peak time and flux of such an afterglow as seen from various observing angles. Using the calibrated equations and the estimated parameters of GW170817, we estimate the detectability of afterglows from future double neutron star mergers during the Advanced LIGO/Virgo observation run O3. GW170817 took place at a relatively low-density environment. Afterglows of similar events will be detectable only at small viewing angles, $ \theta_{\rm obs} \lesssim 20^\circ $, and only $\sim 20\% $ of the GW detections of these events will be accompanied by a detectable afterglow. At higher densities, more typical to sGRB sites, up to $ 70\% $ of the GW detections are expected to be followed by a detectable afterglow, typically at $ \theta_{\rm obs} \sim 30^\circ $. We also provide the latest time one should expect an afterglow detection. We find that for typical parameters, if the jet emission had not been detected within about a year after the merger, it is unlikely to be ever detected.

The primary dissipation mechanism in jets of gamma-ray bursts (GRBs), and the high efficiency of the prompt emission are long standing issues. One possibility is strong collimation of a weakly magnetized relativistic jet by the surrounding medium, which can considerably enhance the efficiency of the photospheric emission. We derive a simple analytic criterion for the radiative efficiency of a collimated jet showing that it depends most strongly on the baryon loading. We confirm this analytic result by 3D numerical simulations, and further find that mixing of jet and cocoon material at the collimation throat leads to a substantial stratification of the outflow as well as sporadic loading, even if the injected jet is uniform and continuous. One consequence of this mixing is a strong angular dependence of the radiative efficiency. Another is large differences in the Lorentz factor of different fluid elements that lead to formation of internal shocks. Our analysis indicates that in both long and short GRBs a prominent photospheric component cannot be avoided when observed within an angle of a few degrees to the axis, unless the asymptotic Lorentz factor is limited by baryon loading at the jet base to $\Gamma_\infty <100$ (with a weak dependence on outflow power). Photon generation by newly created pairs behind the collimation shock regulates the observed temperature at $\sim 50~\theta_0^{-1}$ keV, where $\theta_0$ is the initial jet opening angle, in remarkable agreement with the observed peak energies of prompt emission spectra. Further consequences for the properties of the prompt emission are discussed at the end.

The binary neutron star merger GW170817 was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 41+/-3 Mpc. The radio and X-ray afterglows of GW170817 exhibited delayed onset, a gradual rise in the emission with time as t^0.8, a peak at about 150 days post-merger, followed by a relatively rapid decline. To date, various models have been proposed to explain the afterglow emission, including a choked-jet cocoon and a successful-jet cocoon (a.k.a. structured jet). However, the observational data have remained inconclusive as to whether GW170817 launched a successful relativistic jet. Here we show, through Very Long Baseline Interferometry, that the compact radio source associated with GW170817 exhibits superluminal motion between two epochs at 75 and 230 days post-merger. This measurement breaks the degeneracy between the models and indicates that, while the early-time radio emission was powered by a wider-angle outflow (cocoon), the late-time emission was most likely dominated by an energetic and narrowly-collimated jet, with an opening angle of <5 degrees, and observed from a viewing angle of about 20 degrees. The imaging of a collimated relativistic outflow emerging from GW170817 adds substantial weight to the growing evidence linking binary neutron star mergers and short gamma-ray bursts.

The Hubble constant ($H_0$) measures the current expansion rate of the Universe, and plays a fundamental role in cosmology. Tremendous effort has been dedicated over the past decades to measure $H_0$. Notably, Planck cosmic microwave background (CMB) and the local Cepheid-supernovae distance ladder measurements determine $H_0$ with a precision of $\sim 1\%$ and $\sim 2\%$ respectively. A $3$-$\sigma$ level of discrepancy exists between the two measurements, for reasons that have yet to be understood. Gravitational wave (GW) sources accompanied by electromagnetic (EM) counterparts offer a completely independent standard siren (the GW analogue of an astronomical standard candle) measurement of $H_0$, as demonstrated following the discovery of the neutron star merger, GW170817. This measurement does not assume a cosmological model and is independent of a cosmic distance ladder. The first joint analysis of the GW signal from GW170817 and its EM localization led to a measurement of $H_0=74^{+16}_{-8}$ km/s/Mpc (median and symmetric $68\%$ credible interval). In this analysis, the degeneracy in the GW signal between the source distance and the weakly constrained viewing angle dominated the $H_0$ measurement uncertainty. Recently, Mooley et al. (2018) obtained tight constraints on the viewing angle using high angular resolution imaging of the radio counterpart of GW170817. Here we obtain a significantly improved measurement $H_0=68.9^{+4.7}_{-4.6}$ km/s/Mpc by using these new radio observations, combined with the previous GW and EM data. We estimate that 15 more localized GW170817-like events (comparable signal-to-noise ratio, favorable orientation), having radio images and light curve data, will potentially bring resolution to the tension between the Planck and Cepheid-supernova measurements, as compared to 50-100 GW events without such data.

The short Gamma-Ray Burst, GRB170817A, that followed the binary neutron star merger gravitational waves signal, GW170817, is not a usual sGRB. It is weaker by three orders of magnitude than the weakest sGRB seen before and its spectra, showing a hard early signal followed by a softer thermal spectrum, is unique. We show, first, that the $\gamma$-rays must have emerged from at least mildly relativistic outflow, implying that a relativistic jet was launched following the merger. We then show that the observations are consistent with the predictions of a mildly relativistic shock breakout: a minute $\gamma$-ray energy as compared with the total energy and a rather smooth light curve with a hard to soft evolution. We present here a novel analytic study and detailed numerical 2D and 3D relativistic hydrodynamic and radiation simulations that support the picture in which the observed $\gamma$-rays arose from a shock breakout of a cocoon from the merger's ejecta (Kasliwal 2017). The cocoon can be formed by either a choked jet which does not generate a sGRB (in any direction) or by a successful jet which generates an undetected regular sGRB along the system's axis pointing away from us. Remarkably, for the choked jet model, the macronova signal produced by the ejecta (which is partially boosted to high velocities by the cocoon's shock) and the radio that is produced by the interaction of the shocked cocoon material with the surrounding matter, agree with the observed UV/optical/IR emission and with current radio observations. Finally, we discuss the possibility that the jet propagation within the ejecta may photodissociate some of of the heavy elements and may affect the composition of a fraction of ejecta and, in turn, the opacity and the early macronova light.

Relativistic jets reside in high-energy astrophysical systems of all scales. Their interaction with the surrounding media is critical as it determines the jet evolution, observable signature, and feedback on the environment. During its motion the interaction of the jet with the ambient media inflates a highly pressurized cocoon, which under certain conditions collimates the jet and strongly affects its propagation. Recently, Bromberg et al. (2011) derived a general simplified (semi)analytic solution for the evolution of the jet and the cocoon in case of an unmagnetized jet that propagates in a medium with a range of density profiles. In this work we use a large suite of 2D and 3D relativistic hydrodynamic simulations in order to test the validity and accuracy of this model. We discuss the similarities and differences between the analytic model and numerical simulations and also, to some extent, between 2D and 3D simulations. Our main finding is that although the analytic model is highly simplified, it properly predicts the evolution of the main ingredients of the jet-cocoon system, including its temporal evolution and the transition between various regimes (e.g., collimated to uncollimated). The analytic solution predicts a jet head velocity that is faster by a factor of about 3 compared to the simulations, as long as the head velocity is Newtonian. We use the results of the simulations to calibrate the analytic model which significantly increases its accuracy. We provide an applet that calculates semi-analytically the propagation of a jet in an arbitrary density profile defined by the user at http://www.astro.tau.ac.il/~ore/propagation.html.

The gravitational waves from the first binary neutron star merger, GW170817, were accompanied by a multi-wavelength electromagnetic counterpart, from $\gamma$-rays to radio. The accompanying gamma-rays, seems at first to confirm the association of mergers with short gamma-ray bursts (sGRBs). The common interpretation was that we see an emission from an sGRB jet seen off-axis. However, a closer examination of the sub-luminous $\gamma$-rays and the peculiar radio afterglow were inconsistent with this simple interpretation. Here we present results of 3D and 2D numerical simulations that follow the hydrodynamics and emission of the outflow from a neutron star merger form its ejection and up to its deceleration by the circum-merger medium. Our results show that the entire set of $\gamma$-rays, X-rays and radio observations can be explained by the emission from a mildly relativistic cocoon material (Lorentz factor $\sim$2-5) that was formed while a jet propagated through the material ejected during the merger. The $\gamma$-rays are generated when the cocoon breaks out from the engulfing ejecta while the afterglow is produced by interaction of the cocoon matter with the interstellar medium. The strong early uv/optical signal may be a Lorentz boosted macronova/kilonova. The fate of the jet itself is currently unknown, but our full-EM models define a path to resolving between successful and choked jet scenarios, outputting coupled predictions for the image size, morphology, observed time-dependent polarization and light curve behavior from radio to X-ray. The predictive power of these models will prove key in interpreting the on-going multi-faceted observations of this unprecedented event.

GW170817 is the first gravitational wave detection of a binary neutron star merger. It was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 40 Mpc. It has been proposed that the observed gamma-ray, X-ray and radio emission is due to an ultra-relativistic jet launched during the merger, directed away from our line of sight. The presence of such a jet is predicted from models positing neutron star mergers as the central engines driving short-hard gamma-ray bursts (SGRBs). Here we show that the radio light curve of GW170817 has no direct signature of an off-axis jet afterglow. While we cannot rule out the existence of a jet pointing elsewhere, the observed gamma-rays could not have originated from such a jet. Instead, the radio data requires a mildly relativistic wide-angle outflow moving towards us. This outflow could be the high velocity tail of the neutron-rich material dynamically ejected during the merger or a cocoon of material that breaks out when a jet transfers its energy to the dynamical ejecta. The cocoon scenario can explain the radio light curve of GW170817 as well as the gamma-rays and X-rays (possibly also ultraviolet and optical emission), and hence is the model most consistent with the observational data. Cocoons may be a ubiquitous phenomenon produced in neutron star mergers, giving rise to a heretofore unidentified population of radio, ultraviolet, X-ray and gamma-ray transients in the local universe.

Short Gamma-Ray Bursts (SGRBs) are believed to arise from compact binary mergers (either neutron star-neutron star or black hole-neutron star). If so their jets must penetrate outflows that are ejected during the merger. As a jet crosses the ejecta it dissipates its energy, producing a hot cocoon which surrounds it. We present here 3D numerical simulations of jet propagation in mergers' outflows and we calculate the resulting emission. This emission consists of two components: the cooling emission, the leakage of the thermal energy of the hot cocoon, and the cocoon macronova that arises from the radioactive decay of the cocoon's material. This emission gives a brief (~ one hour) blue, wide angle signal. While the parameters of the outflow and jet are uncertain, for the configurations we have considered the signal is bright (~ -14 $-$ -15 absolute magnitude) and outshines all other predicted UV-optical signals. The signal is brighter when the jet breakout time is longer and its peak brightness does not depend strongly on the highly uncertain opacity. A rapid search for such a signal is a promising strategy to detect an electromagnetic merger counterpart. A detected candidate could be then followed by deep IR searches for the longer but weaker macronova arising from the rest of the ejecta.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

We present the first relativistic MHD numerical simulation of a magnetic jet that propagates through and emerges from the dynamical ejecta of a binary neutron star merger. Generated by the magnetized rotation of the merger remnant, the jet propagates through the ejecta and produces an energetic cocoon that expands at mildly relativistic velocities and breaks out of the ejecta. We show that if the ejecta has a low-mass ($\sim10^{-7} M_\odot$) high-velocity ($v\sim0.85$ c) tail, the cocoon shock breakout will generate $\gamma$-ray emission that is comparable to the observed short GRB170817A that accompanied the recent gravitational wave event GW170817. Thus, we propose that this GRB, which is quite different from all other short GRBs observed before, was produced by a different mechanism. We expect, however, that such events are numerous and many will be detected in coming LIGO-Virgo runs.