More Like this

1. Jets with a Twist: The Emergence of FR0 Jets in a 3D GRMHD Simulation of Zero-angular-momentum Black Hole Accretion

   https://doi.org/10.3847/1538-4357/ad0974  

   Lalakos, Aretaios; Tchekhovskoy, Alexander; Bromberg, Omer; Gottlieb, Ore; Jacquemin-Ide, Jonatan; Liska, Matthew; Zhang, Haocheng (March 2024, The Astrophysical Journal)

   Abstract Spinning supermassive black holes (BHs) in active galactic nuclei magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within 3 orders of magnitude in distance from the BH, well before reaching kiloparsec scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by a weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen to be much larger than the BH gravitational radius,R_B= 10³R_g. The BH develops dynamically important large-scale magnetic fields, forms a magnetically arrested disk (MAD), and launches relativistic jets that propagate well outsideR_Band suppress BH accretion to 1.5% of the Bondi rate, ${\dot{M}}_{\mathrm{B}}$. Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff–Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets 5%–10% of ${\dot{M}}_{\mathrm{B}}$reach the BH. Thereafter, the disk starts rocking back and forth by angles 90°–180°: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to ≲2% ${\dot{M}}_{\mathrm{B}}$. Because the BAD and RAD states tangle up the jets and destroy them well insideR_B, they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies. 

   more » « less
2. Bridging Scales in Black Hole Accretion and Feedback: Magnetized Bondi Accretion in 3D GRMHD

   https://doi.org/10.3847/2041-8213/ad1048  

   Cho 조, Hyerin 혜린; Prather, Ben S.; Narayan, Ramesh; Natarajan, Priyamvada; Su, Kung-Yi; Ricarte, Angelo; Chatterjee, Koushik (December 2023, The Astrophysical Journal Letters)

   Abstract Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magnetohydrodynamic (GRMHD) code KHARMA that allows us to span 7 orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with a Bondi radius ofR_B≈ 2 × 10⁵GM_•/c², whereM_•is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as ∼r⁻¹rather thanr^−3/2and the accretion rate $\dot{M}$is consequently suppressed by over 2 orders of magnitude below the Bondi rate ${\dot{M}}_{\mathrm{B}}$. We find continuous energy feedback from the accretion flow to the external medium at a level of $\sim {10}^{-2}\dot{M}{c}^2\sim 5\times {10}^{-5}{\dot{M}}_{\mathrm{B}}{c}^2$. Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback. 

   more » « less
3. A Unified Picture of Short and Long Gamma-Ray Bursts from Compact Binary Mergers

   https://doi.org/10.3847/2041-8213/ad096e  

   Gottlieb, Ore; Metzger, Brian D; Quataert, Eliot; Issa, Danat; Martineau, Tia; Foucart, Francois; Duez, Matthew D; Kidder, Lawrence E; Pfeiffer, Harald P; Scheel, Mark A (November 2023, The Astrophysical Journal Letters)

   Abstract The recent detections of the ∼10 s longγ-ray bursts (GRBs) 211211A and 230307A followed by softer temporally extended emission (EE) and kilonovae point to a new GRB class. Using state-of-the-art first-principles simulations, we introduce a unifying theoretical framework that connects binary neutron star (BNS) and black hole–NS (BH–NS) merger populations with the fundamental physics governing compact binary GRBs (cbGRBs). For binaries with large total masses,M_tot≳ 2.8M_⊙, the compact remnant created by the merger promptly collapses into a BH surrounded by an accretion disk. The duration of the pre-magnetically arrested disk (MAD) phase sets the duration of the roughly constant power cbGRB and could be influenced by the disk mass,M_d. We show that massive disks (M_d≳ 0.1M_⊙), which form for large binary mass ratiosq≳ 1.2 in BNS orq≲ 3 in BH–NS mergers, inevitably produce 211211A-like long cbGRBs. Once the disk becomes MAD, the jet power drops with the mass accretion rate as $\dot{M}\sim t^{-2}$, establishing the EE decay. Two scenarios are plausible for short cbGRBs. They can be powered by BHs with less massive disks, which form for otherqvalues. Alternatively, for binaries withM_tot≲ 2.8M_⊙, mergers should go through a hypermassive NS (HMNS) phase, as inferred for GW170817. Magnetized outflows from such HMNSs, which typically live for ≲1 s, offer an alternative progenitor for short cbGRBs. The first scenario is challenged by the bimodal GRB duration distribution and the fact that the Galactic BNS population peaks at sufficiently low masses that most mergers should go through an HMNS phase. 

   more » « less
4. Magnetically Driven Neutron-rich Ejecta Unleashed: Global 3D Neutrino–General Relativistic Magnetohydrodynamic Simulations of Collapsars Probe the Conditions for r -process Nucleosynthesis

   https://doi.org/10.3847/2041-8213/adc694  

   Issa, Danat; Gottlieb, Ore; Metzger, Brian D; Jacquemin-Ide, Jonatan; Liska, Matthew; Foucart, Francois; Halevi, Goni; Tchekhovskoy, Alexander (May 2025, The Astrophysical Journal Letters)

   Abstract Collapsars—rapidly rotating stellar cores that form black holes—can power gamma-ray bursts and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process (r-process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included sufficiently strong magnetic fields and the black hole (BH), both of which are essential for launching mass outflows. We presentνh-amr, a novel neutrino-transport general relativistic magnetohydrodynamic (νGRMHD) code, which we use to perform the first 3D globalνGRMHD collapsar simulations. We find a self-consistent formation of a weakly magnetized dense accretion disk, which has sufficient time to neutronize. Eventually, substantial magnetic flux accumulates near the BH, becomes dynamically important, leads to a magnetically arrested disk (MAD), and unbinds some of the neutron-rich material. However, the strong flux also hinders accretion, lowers density, and increases neutrino-cooling timescale, which prevents further disk neutronization. Typical collapsar progenitors with mass accretion rates, $\dot{M}\sim 0.1-1M_{\odot }{\mathrm{s}}^{-1}$, do not produce significant neutron-rich (Y_e < 0.25) ejecta. However, we find that MADs at higher mass accretion rates, $\dot{M}\gtrsim \text{few}{M}_{\odot }{\mathrm{s}}^{-1}$(e.g., for more centrally concentrated progenitors), can unbindM_ej ≲ M_⊙of neutron-rich ejecta. The outflows inflate a shocked cocoon that mixes with the infalling neutron-poor stellar gas and raises the final outflowY_e; however, the finalr-process yield may be determined earlier at the point of neutron capture freeze-out. Future work will explore under what conditions more typical collapsar engines becomer-process factories. 

   more » « less
5. Maximal Jet Energy of Gamma-Ray Bursts through the Blandford–Znajek Mechanism

   https://doi.org/10.3847/2041-8213/adaeb8  

   Wu, Zhao-Feng; Damoulakis, Michail; Beniamini, Paz; Giannios, Dimitrios (February 2025, The Astrophysical Journal Letters)

   Abstract Gamma-ray bursts (GRBs) are among the most energetic events in the Universe, driven by relativistic jets launched from black holes (BHs) formed during the collapse of massive stars or after the merger of two neutron stars. The jet power depends on the BH spin and the magnetic flux accreted onto it. In the standard thin disk model, jet power is limited by insufficient magnetic flux, even when the spin approaches maximum possible value. In contrast, the magnetically arrested disk (MAD) state limits jet energy by extracting significant angular momentum, braking BH rotation. We propose a unified model incorporating both standard thin disk and MAD states, identifying a universal curve for jet power per accretion rate as a function of the magnetic flux ratio, ${\Delta }_{\mathrm{eq}}={({\Phi }_{\mathrm{BH}}/{\Phi }_{\mathrm{MAD}})}_{\mathrm{eq}}$, at spin equilibrium. For long GRBs (lGRBs), the model predicts a maximum jet energy of ∼1.5% of the accretion energy, occurring at Δ_eq ∼ 0.4, where the BH equilibrium spin isa ∼ 0.5. Both long and short GRBs are unlikely to be produced by a MAD: for short GRBs, this requires an accreted mass orders of magnitude smaller than that available, while for lGRBs, the narrow progenitor mass distribution challenges the ability to produce the observed broad distribution of jet energies. This framework provides a consistent explanation for both standard and luminous GRBs, emphasizing the critical role of magnetic flux. Both long and short GRBs require magnetic flux distributions that peak around 10²⁷G cm². 

   more » « less