The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $gamma$-
ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.
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 collapse of a massive star with low angular momentum content is commonly argued to result in the formation of a black hole without an accompanying bright transient. Our goal in this Letter is to understand the flow in and around a newly-formed bl
ack hole, involving accretion and rotation, via general relativistic hydrodynamics simulations aimed at studying the conditions under which infalling material can accrete without forming a centrifugally supported structure and, as a result, generate no effective feedback. If the feedback from the black hole is, on the other hand, significant, the collapse would be halted and we suggest that the event is likely to be followed by a bright transient. We find that feedback is only efficient if the specific angular momentum of the infalling material at the innermost stable circular orbit exceeds that of geodesic circular flow at that radius by at least $approx 20%$. We use the results of our simulations to constrain the maximal stellar rotation rates of the disappearing massive progenitors PHL293B-LBV and N6946-BH1, and to provide an estimate of the overall rate of disappearing massive stars. We find that about a few percent of single O-type stars with measured rotational velocities are expected to spin below the critical value before collapse and are thus predicted to vanish without a trace.
The most popular model for short gamma-ray bursts (sGRBs) involves the coalescence of binary neutron stars. Because the progenitor is actually hidden from view, we must consider under which circumstances such merging systems are capable of producing
a successful sGRB. Soon after coalescence, winds are launched from the merger remnant. In this paper, we use realistic wind profiles derived from global merger simulations in order to investigate the interaction of sGRB jets with these winds using numerical simulations. We analyze the conditions for which these axisymmetric winds permit relativistic jets to breakout and produce a sGRB. We find that jets with luminosities comparable to those observed in sGRBs are only successful when their half-opening angles are below ~20{deg}. This jet collimation mechanism leads to a simple physical interpretation of the luminosities and opening angles inferred for sGRBs. If wide, low luminosity jets are observed, they might be indicative of a different progenitor avenue such as the merger of a neutron star with a black hole. We also use the observed durations of sGRB to place constraints on the lifetime of the wind phase, which is determined by the time it takes the jet to breakout. In all cases we find that the derived limits argue against completely stable remnants for binary neutron star mergers that produce sGRBs.
The central engine of short gamma-ray bursts (sGRBs) is hidden from direct view, operating at a scale much smaller than that probed by the emitted radiation. Thus we must infer its origin not only with respect to the formation of the {it trigger} - t
he actual astrophysical configuration that is capable of powering a sGRB - but also from the consequences that follow from the various evolutionary pathways that may be involved in producing it. Considering binary neutron star mergers we critically evaluate, analytically and through numerical simulations, whether the neutrino-driven wind produced by the newly formed hyper-massive neutron star can allow the collimated relativistic outflow that follows its collapse to actually produce a sGRB or not. Upon comparison with the observed sGRB duration distribution, we find that collapse cannot be significantly delayed (<=100 ms) before the outflow is choked, thus limiting the possibility that long-lived hyper-massive remnants can account for these events. In the case of successful breakthrough of the jet through the neutrino-driven wind, the energy stored in the cocoon could contribute to the precursor and extended emission observed in sGRBs.
