skip to main content


Title: Black Hole to Photosphere: 3D GRMHD Simulations of Collapsars Reveal Wobbling and Hybrid Composition Jets
Abstract Long-duration γ -ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth via black hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over six orders of magnitude in space and time. The collapsing stellar envelope forms an accretion disk, which drags inwardly the magnetic flux that accumulates around the BH, becomes dynamically important, and launches bipolar jets. The jets reach the photosphere at ∼10 12 cm with an opening angle θ j ∼ 6° and a Lorentz factor Γ j ≲ 30, unbinding ≳90% of the star. We find that (i) the disk–jet system spontaneously develops misalignment relative to the BH rotational axis. As a result, the jet wobbles with an angle θ t ∼ 12°, which can naturally explain quiescent times in GRB lightcurves. The effective opening angle for detection θ j + θ t suggests that the intrinsic GRB rate is lower by an order of magnitude than standard estimates. This suggests that successful GRBs are rarer than currently thought and emerge in only ∼0.1% of supernovae Ib/c, implying that jets are either not launched or choked inside most supernova Ib/c progenitors. (ii) The magnetic energy in the jet decreases due to mixing with the star, resulting in jets with a hybrid composition of magnetic and thermal components at the photosphere, where ∼10% of the gas maintains magnetization σ ≳ 0.1. This indicates that both a photospheric component and reconnection may play a role in the prompt emission.  more » « less
Award ID(s):
2107802 2107839 1815304 2031997
NSF-PAR ID:
10349747
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
The Astrophysical Journal Letters
Volume:
933
Issue:
1
ISSN:
2041-8205
Page Range / eLocation ID:
L9
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. ABSTRACT

    We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, l, of the collapsing material and the magnetic field, B, on the BH horizon: (i) subrelativistic outflow (low l and high B), (ii) stationary accretion shock instability (SASI; high l and low B), (iii) relativistic jets (high l and high B). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices (≳ 2), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the light curve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova + GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disc-jet flows, even in the absence of any tilt in the progenitor.

     
    more » « less
  4. Context. There has been significant technological and scientific progress in our ability to detect, monitor, and model the physics of γ -ray bursts (GRBs) over the 50 years since their first discovery. However, the dissipation process thought to be responsible for their defining prompt emission is still unknown. Recent efforts have focused on investigating how the ultrarelativistic jet of the GRB propagates through the progenitor’s stellar envelope for different initial composition shapes, jet structures, magnetisation, and, consequently, possible energy dissipation processes. Study of the temporal variability – in particular the shortest duration of an independent emission episode within a GRB – may provide a unique way to distinguish the imprint of the inner engine activity from geometry and propagation related effects. The advent of new high-energy detectors with exquisite time resolution now makes this possible. Aims. We aim to characterise the minimum variability timescale (MVT) defined as the shortest duration of individual pulses that shape a light curve for a sample of GRBs in the keV–MeV energy range and test correlations with other key observables such as the peak luminosity, the Lorentz factor, and the jet opening angle. We compare these correlations with predictions from recent numerical simulations for a relativistic structured – possibly wobbling – jet and assess the value of temporal variability studies as probes of prompt-emission dissipation physics. Methods. We used the peak detection algorithm MEPSA to identify the shortest pulse within a GRB time history and preliminarily calibrated MEPSA to estimate the full width at half maximum duration. We then applied this framework to two sets of GRBs: Swift GRBs (from 2005 to July 2022) and Insight Hard Modulation X-ray Telescope (Insight-HXMT) GRBs (from June 2017 to July 2021, including the exceptional 221009A). We then selected 401 GRBs with measured redshift to test for correlations. Results. We confirm that, on average, short GRBs have significantly shorter MVTs than long GRBs. The MVT distribution of short GRBs with extended emission such as 060614 and 211211A is compatible only with that of short GRBs. This is important because it provides a new clue concerning the progenitor’s nature. The MVT for long GRBs with measured redshift anti-correlates with peak luminosity; our analysis includes careful evaluation of selection effects. We confirm the anti-correlation with the Lorentz factor and find a correlation with the jet opening angle as estimated from the afterglow light curve, along with an inverse correlation with the number of pulses. Conclusions. The MVT can identify the emerging putative new class of long GRBs that are suggested to be produced by compact binary mergers. For otherwise typical long GRBs, the different correlations between MVT and peak luminosity, Lorentz factor, jet opening angle, and number of pulses can be explained within the context of structured, possibly wobbling, weakly magnetised relativistic jets. 
    more » « less
  5. Abstract

    We present a comprehensive study of 29 short gamma-ray bursts (SGRBs) observed ≈0.8−60 days postburst using Chandra and XMM-Newton. We provide the inferred distributions of the SGRB jet opening angles and true event rates to compare against neutron star merger rates. We perform a uniform analysis and modeling of their afterglows, obtaining 10 opening angle measurements and 19 lower limits. We report on two new opening angle measurements (SGRBs 050724A and 200411A) and eight updated values, obtaining a median value of 〈θj〉 ≈ 6.°1 [−3.°2, +9.°3] (68% confidence on the full distribution) from jet measurements alone. For the remaining events, we inferθj≳ 0.°5–26°. We uncover a population of SGRBs with wider jets ofθj≳ 10° (including two measurements ofθj≳ 15°), representing ∼28% of our sample. Coupled with multiwavelength afterglow information, we derive a total true energy of 〈Etrue,tot〉 ≈ 1049–1050erg, which is consistent with magnetohydrodynamic jet launching mechanisms. Furthermore, we determine a range for the beaming-corrected event rate ofRtrue3601800Gpc−3yr−1, set by the inclusion of a population of wide jets on the low end, and the jet measurements alone on the high end. From a comparison with the latest merger rates, our results are consistent with the majority of SGRBs originating from binary neutron star mergers. However, our inferred rates are well above the latest neutron star–black hole merger rates, consistent with at most a small fraction of SGRBs originating from such mergers.

     
    more » « less