skip to main content


Title: On the Jet–Ejecta Interaction in 3D GRMHD Simulations of a Binary Neutron Star Merger Aftermath
Abstract Short γ -ray burst (sGRB) jets form in the aftermath of a neutron star merger, drill through disk winds and dynamical ejecta, and extend over four to five orders of magnitude in distance before breaking out of the ejecta. We present the first 3D general-relativistic magnetohydrodynamic sGRB simulations to span this enormous scale separation. They feature three possible outcomes: jet+cocoon, cocoon, and neither. Typical sGRB jets break out of the dynamical ejecta if (i) the bound ejecta’s isotropic equivalent mass along the pole at the time of the BH formation is ≲10 −4 M ⊙ , setting a limit on the delay time between the merger and BH formation, otherwise, the jets perish inside the ejecta and leave the jet-inflated cocoon to power a low-luminosity sGRB; (ii) the postmerger remnant disk contains a strong large-scale vertical magnetic field, ≳10 15 G; and (iii) if the jets are weak (≲10 50 erg), the ejecta’s isotropic equivalent mass along the pole must be small (≲10 −2 M ⊙ ). Generally, the jet structure is shaped by the early interaction with disk winds rather than the dynamical ejecta. As long as our jets break out of the ejecta, they retain a significant magnetization (≲1), suggesting that magnetic reconnection is a fundamental property of sGRB emission. The angular structure of the outflow isotropic equivalent energy after breakout consistently features a flat core followed by a steep power-law distribution (slope ≳3), similar to hydrodynamic jets. In the cocoon-only outcome, the dynamical ejecta broadens the outflow angular distribution and flattens it (slope ∼1.5).  more » « less
Award ID(s):
2107839 1815304 2031997
NSF-PAR ID:
10353390
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
The Astrophysical Journal Letters
Volume:
933
Issue:
1
ISSN:
2041-8205
Page Range / eLocation ID:
L2
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present the first numerical simulations that track the evolution of a black hole–neutron star (BH–NS) merger from premerger tor≳ 1011cm. The disk that forms after a merger of mass ratioq= 2 ejects massive disk winds (3–5 × 10−2M). We introduce various postmerger magnetic configurations and find that initial poloidal fields lead to jet launching shortly after the merger. The jet maintains a constant power due to the constancy of the large-scale BH magnetic flux until the disk becomes magnetically arrested (MAD), where the jet power falls off asLjt−2. All jets inevitably exhibit either excessive luminosity due to rapid MAD activation when the accretion rate is high or excessive duration due to delayed MAD activation compared to typical short gamma-ray bursts (sGRBs). This provides a natural explanation for long sGRBs such as GRB 211211A but also raises a fundamental challenge to our understanding of jet formation in binary mergers. One possible implication is the necessity of higher binary mass ratios or moderate BH spins to launch typical sGRB jets. For postmerger disks with a toroidal magnetic field, dynamo processes delay jet launching such that the jets break out of the disk winds after several seconds. We show for the first time that sGRB jets with initial magnetizationσ0> 100 retain significant magnetization (σ≫ 1) atr> 1010cm, emphasizing the importance of magnetic processes in the prompt emission. The jet–wind interaction leads to a power-law angular energy distribution by inflating an energetic cocoon whose emission is studied in a companion paper.

     
    more » « less
  2. ABSTRACT

    We consider conditions for jet breakout through ejecta following mergers of neutron stars and provide simple relations for the breakout conditions. We demonstrate that: (i) break-out requires that the isotropic-equivalent jet energy Ej exceeds the ejecta energy Eej by Ej ≥ Eej/βej, where βej = Vej/c, Vej is the maximum velocity of the ejecta. If the central engine terminates before the breakout, the shock approaches the edge of the ejecta slowly ∝ 1/t; late breakout occurs only if at the termination moment the head of the jet was relatively close to the edge. (ii) If there is a substantial delay between the ejecta’s and the jet’s launching, the requirement on the jet power increases. (iii) The forward shock driven by the jet is mildly strong, with Mach number M ≈ 5/4 (increasing with time delay td); (iii) the delay time td between the ejecta and the jet’s launching is important for $t_\mathrm{ d} \gt t_0= ({3 }/{16}) {c M_{\mathrm{ ej}} V_{\mathrm{ ej}}}/{L_\mathrm{ j}} = 1.01 {\rm \mathrm{ s}} M_{\mathrm{ ej}, -2} L_{\mathrm{ j}, 51} ^{-1} \left({\beta _{\mathrm{ ej}}} /{0.3} \right)$, where Mej is ejecta mass, Lj is the jet luminosity (isotropic equivalent). For small delays, t0 is also an estimate of the break-out time.

     
    more » « less
  3. Abstract

    The spin of a newly formed black hole (BH) at the center of a massive star evolves from its natal value due to two competing processes: accretion of gas angular momentum that increases the spin and extraction of BH angular momentum by outflows that decreases the spin. Ultimately, the final, equilibrium spin is set by a balance between both processes. In order for the BH to launch relativistic jets and power aγ-ray burst (GRB), the BH magnetic field needs to be dynamically important. Thus, we consider the case of a magnetically arrested disk (MAD) driving the spin evolution of the BH. By applying the semianalytic MAD BH spin evolution model of Lowell et al. to collapsars, we show that if the BH accretes ∼20% of its initial mass, its dimensionless spin inevitably reaches small values,a≲ 0.2. For such spins, and for mass accretion rates inferred from collapsar simulations, we show that our semianalytic model reproduces the energetics of typical GRB jets,Ljet∼ 1050erg s−1. We show that our semianalytic model reproduces the nearly constant power of typical GRB jets. If the MAD onset is delayed, this allows powerful jets at the high end of the GRB luminosity distribution,Ljet∼ 1052erg s−1, but the final spin remains low,a≲ 0.3. These results are consistent with the low spins inferred from gravitational wave detections of binary BH mergers. In a companion paper by Gottlieb et al., we use GRB observations to constrain the natal BH spin to bea≃ 0.2.

     
    more » « less
  4. Context. Accreting black holes (BHs) may be surrounded by a highly magnetized plasma threaded by an organized poloidal magnetic field. Nonthermal flares and power-law spectral components at high energy could originate from a hot, collisionless, and nearly force-free corona. The jets we often observe from these systems are believed to be rotation-powered and magnetically driven. Aims. We study axisymmetric BH magnetospheres, where a fraction of the magnetic field lines anchored in a surrounding disk are connected to the event horizon of a rotating BH. For different BH spins, we identify the conditions and sites of magnetic reconnection within 30 gravitational radii. Methods. With the fully general relativistic particle-in-cell code GRZeltron , we solve the time-dependent dynamics of the electron–positron pair plasma and of the electromagnetic fields around the BH. The aligned disk is represented by a steady and perfectly conducting plasma in Keplerian rotation, threaded by a dipolar magnetic field. Results. For prograde disks around Kerr BHs, the topology of the magnetosphere is hybrid. Twisted open magnetic field lines crossing the horizon power a Blandford-Znajek jet, while open field lines with their footpoint beyond a critical distance on the disk could launch a magneto-centrifugal wind. In the innermost regions, coupling magnetic field lines ensure the transfer of significant amounts of angular momentum and energy between the BH and the disk. From the Y point at the intersection of these three regions, a current sheet forms where vivid particle acceleration via magnetic reconnection takes place. We compute the synchrotron images of the current sheet emission. Conclusions. Our estimates for jet power and BH–disk exchanges match those derived from purely force-free models. Particles are accelerated at the Y point, which acts as a heat source for the so-called corona. It provides a physically motivated ring-shaped source of hard X-rays above the disk for reflection models. Episodic plasmoid ejection might explain millisecond flares observed in Cygnus X-1 in the high-soft state, but are too fast to account for daily nonthermal flares from Sgr A * . Particles flowing from the Y point down to the disk could produce a hot spot at the footpoint of the outermost closed magnetic field line. 
    more » « less
  5. ABSTRACT

    Evidence is mounting that recent multiwavelength detections of fast blue optical transients (FBOTs) in star-forming galaxies comprise a new class of transients, whose origin is yet to be understood. We show that hydrogen-rich collapsing stars that launch relativistic jets near the central engine can naturally explain the entire set of FBOT observables. The jet–star interaction forms a mildly relativistic shocked jet (inner cocoon) component, which powers cooling emission that dominates the high velocity optical signal during the first few weeks, with a typical energy of ∼1050–1051 erg. During this time, the cocoon radial energy distribution implies that the optical light curve exhibits a fast decay of $L \,\, \buildrel\propto \over \sim \,\,t^{-2.4}$. After a few weeks, when the velocity of the emitting shell is ∼0.01 c, the cocoon becomes transparent, and the cooling envelope governs the emission. The interaction between the cocoon and the dense circumstellar winds generates synchrotron self-absorbed emission in the radio bands, featuring a steady rise on a month time-scale. After a few months the relativistic outflow decelerates, enters the observer’s line of sight, and powers the peak of the radio light curve, which rapidly decays thereafter. The jet (and the inner cocoon) becomes optically thin to X-rays ∼day after the collapse, allowing X-ray photons to diffuse from the central engine that launched the jet to the observer. Cocoon cooling emission is expected at higher volumetric rates than gamma-ray bursts (GRBs) by a factor of a few, similar to FBOTs. We rule out uncollimated outflows, however, both GRB jets and failed collimated jets are compatible with all observables.

     
    more » « less