More Like this

1. 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
2. Stars Crushed by Black Holes. II. A Physical Model of Adiabatic Compression and Shock Formation in Tidal Disruption Events

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

   Coughlin, Eric R.; Nixon, C. J. (February 2022, The Astrophysical Journal)

   Abstract We develop a Newtonian model of a deep tidal disruption event (TDE), for which the pericenter distance of the star,r_p, is well within the tidal radius of the black hole,r_t, i.e., whenβ≡r_t/r_p≫ 1. We find that shocks form forβ≳ 3, but they are weak (with Mach numbers ∼1) for allβ, and that they reach the center of the star prior to the time of maximum adiabatic compression forβ≳ 10. The maximum density and temperature reached during the TDE follow much shallower relations withβthan the previously predicted ${\rho }_{\max }\propto {\beta }^3$and $T_{\max }\propto {\beta }^2$scalings. Belowβ≃ 10, this shallower dependence occurs because the pressure gradient is dynamically significant before the pressure is comparable to the ram pressure of the free-falling gas, while aboveβ≃ 10, we find that shocks prematurely halt the compression and yield the scalings ${\rho }_{\max }\propto {\beta }^{1.62}$and $T_{\max }\propto {\beta }^{1.12}$. We find excellent agreement between our results and high-resolution simulations. Our results demonstrate that, in the Newtonian limit, the compression experienced by the star is completely independent of the mass of the black hole. We discuss our results in the context of existing (affine) models, polytropic versus non-polytropic stars, and general relativistic effects, which become important when the pericenter of the star nears the direct capture radius, atβ∼ 12.5 (2.7) for a solar-like star disrupted by a 10⁶M_⊙(10⁷M_⊙) supermassive black hole. 

   more » « less
3. Scaling of Magnetic Reconnection Electron Bulk Heating in the High-Alfvén-speed and Low-β Regime of Earth’s Magnetotail

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

   Øieroset, M; Phan, T D; Oka, M; Drake, J F; Fuselier, S A; Gershman, D J; Maheshwari, K; Giles, B L; Zhang, Q; Guo, F; et al (August 2023, The Astrophysical Journal)

   Abstract We have surveyed 21 reconnection exhaust events observed by Magnetospheric MultiScale in the low-plasma-βand high-Alfvén-speed regime of the Earth’s magnetotail to investigate the scaling of electron bulk heating produced by reconnection. The ranges of inflow Alfvén speed and inflow electronβ_ecovered by this study are 800–4000 km s⁻¹and 0.001–0.1, respectively, and the observed heating ranges from a few hundred electronvolts to several kiloelectronvolts. We find that the temperature change in the reconnection exhaust relative to the inflow, ΔT_e, is correlated with the inflow Alfvén speed,V_Ax,in, based on the reconnecting magnetic field and the inflow plasma density. Furthermore, ΔT_eis linearly proportional to the inflowing magnetic energy per particle, ${m_{\mathrm{i}}{V}_{\mathrm{Ax},\mathrm{in}}}^2$, and the best fit to the data produces the empirical relation ΔT_e= 0.020 ${m_{\mathrm{i}}{V}_{\mathrm{Ax},\mathrm{in}}}^2$, i.e., the electron temperature increase is on average ∼2% of the inflowing magnetic energy per particle. This magnetotail study extends a previous magnetopause reconnection study by two orders of magnitude in both magnetic energy and electronβ, to a regime that is comparable to the solar corona. The validity of the empirical relation over such a large combined magnetopause–magnetotail plasma parameter range ofV_A∼ 10–4000 km s⁻¹andβ_e∼ 0.001–10 suggests that one can predict the magnitude of the bulk electron heating by reconnection in a variety of contexts from the simple knowledge of a single parameter: the Alfvén speed of the ambient plasma. 

   more » « less
4. Cooling Envelope Model for Tidal Disruption Events

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

   Metzger, Brian D. (September 2022, The Astrophysical Journal Letters)

   Abstract We present a toy model for the thermal optical/UV/X-ray emission from tidal disruption events (TDEs). Motivated by recent hydrodynamical simulations, we assume that the debris streams promptly and rapidly circularize (on the orbital period of the most tightly bound debris), generating a hot quasi-spherical pressure-supported envelope of radiusR_v∼ 10¹⁴cm (photosphere radius ∼10¹⁵cm) surrounding the supermassive black hole (SMBH). As the envelope cools radiatively, it undergoes Kelvin–Helmholtz contractionR_v∝t⁻¹, its temperature risingT_eff∝t^1/2while its total luminosity remains roughly constant; the optical luminosity decays as $\nu L_{\nu }\propto R_v^2{T}_{\mathrm{eff}}\propto t^{-3/2}$. Despite this similarity to the mass fallback rate ${\dot{M}}_{\mathrm{fb}}\propto t^{-5/3}$, envelope heating from fallback accretion is subdominant compared to the envelope cooling luminosity except near optical peak (where they are comparable). Envelope contraction can be delayed by energy injection from accretion from the inner envelope onto the SMBH in a regulated manner, leading to a late-time flattening of the optical/X-ray light curves, similar to those observed in some TDEs. Eventually, as the envelope contracts to near the circularization radius, the SMBH accretion rate rises to its maximum, in tandem with the decreasing optical luminosity. This cooling-induced (rather than circularization-induced) delay of up to several hundred days may account for the delayed onset of thermal X-rays, late-time radio flares, and high-energy neutrino generation, observed in some TDEs. We compare the model predictions to recent TDE light-curve correlation studies, finding both agreement and points of tension. 

   more » « less
5. Suppression of Collisionless Magnetic Reconnection in the High Ion β, Strong Guide Field Limit

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

   Giai, Carlos_A; Haggerty, Colby_C; Shay, Michael_A; Cassak, Paul_A; Davis, Zachary_K (December 2024, The Astrophysical Journal)

   Abstract In magnetic reconnection, the ion bulk outflow speed and ion heating have been shown to be set by the available reconnecting magnetic energy, i.e., the energy stored in the reconnecting magnetic field (B_r). However, recent simulations, observations, and theoretical works have shown that the released magnetic energy is inhibited by upstream ion plasma betaβ_i—the relative ion thermal pressure normalized to magnetic pressure based on the reconnecting field—for antiparallel magnetic field configurations. Using kinetic theory and hybrid particle-in-cell simulations, we investigate the effects ofβ_ion guide field reconnection. While previous works have suggested that guide field reconnection is uninfluenced byβ_i, we demonstrate that the reconnection process is modified and the outflow is reduced for sufficiently large ${\beta }_i>(B_r^2+B_g^2)/B_r^2$. We develop a theoretical framework that shows that this reduction is consistent with an enhanced exhaust pressure gradient, which reduces the outflow speed as $v_{\mathrm{out}}\propto 1/\sqrt{{\beta }_i}$. These results apply to systems in which guide field reconnection is embedded in hot plasmas, such as reconnection at the boundary of eddies in fully developed turbulence like the solar wind or the magnetosheath as well as downstream of shocks such as the heliosheath or the mergers of galaxy clusters. 

   more » « less