More Like this

1. 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)

   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
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)

   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. Cooling Envelope Model for Tidal Disruption Events

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

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

   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
4. The Peak of the Fallback Rate from Tidal Disruption Events: Dependence on Stellar Type

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

   Bandopadhyay, Ananya ; Fancher, Julia ; Athian, Aluel ; Indelicato, Valentino ; Kapalanga, Sarah ; Kumah, Angela ; Paradiso, Daniel A. ; Todd, Matthew ; Coughlin, Eric R. ; Nixon, C. J. ( January 2024 , The Astrophysical Journal Letters)

   A star completely destroyed in a tidal disruption event (TDE) ignites a luminous flare that is powered by the fallback of tidally stripped debris to a supermassive black hole (SMBH) of massM_•. We analyze two estimates for the peak fallback rate in a TDE, one being the “frozen-in” model, which predicts a strong dependence of the time to peak fallback rate,t_peak, on both stellar mass and age, with 15 days ≲t_peak≲ 10 yr for main sequence stars with masses 0.2 ≤M_⋆/M_⊙≤ 5 andM_•= 10⁶M_⊙. The second estimate, which postulates that the star is completely destroyed when tides dominate the maximum stellar self-gravity, predicts thatt_peakis very weakly dependent on stellar type, with$t_{\mathrm{peak}}=\left(23.2\pm 4.0\mathrm{days}\right){\left(M_{•}/{10}^6{M}_{\odot }\right)}^{1/2}$for 0.2 ≤M_⋆/M_⊙≤ 5, while$t_{\mathrm{peak}}=\left(29.8\pm 3.6\mathrm{days}\right){\left(M_{•}/{10}^6{M}_{\odot }\right)}^{1/2}$for a Kroupa initial mass function truncated at 1.5M_⊙. This second estimate also agrees closely with hydrodynamical simulations, while the frozen-in model is discrepant by orders of magnitude. We conclude that (1) the time to peak luminosity in complete TDEs is almost exclusively determined by SMBH mass, and (2) massive-star TDEs power the largest accretion luminosities. Consequently, (a) decades-long extra-galactic outbursts cannot be powered by complete TDEs, including massive-star disruptions, and (b) the most highly super-Eddington TDEs are powered by the complete disruption of massive stars, which—if responsible for producing jetted TDEs—would explain the rarity of jetted TDEs and their preference for young and star-forming host galaxies.

    

   more » « less
5. The Luminosity Function of Tidal Disruption Events from Fallback-powered Emission: Implications for the Black Hole Mass Function

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

   Coughlin, Eric R. ; Nicholl, Matt ( May 2023 , The Astrophysical Journal Letters)

   Tidal disruption events (TDEs), in which a star is destroyed by the gravitational field of a supermassive black hole (SMBH), are being observed at a high rate owing to the advanced state of survey science. One of the properties of TDEs that is measured with increasing statistical reliability is the TDE luminosity function,$d{\dot{N}}_{\mathrm{TDE}}/\mathit{dL}$, which is the TDE rate per luminosity (i.e., how many TDEs are within a given luminosity range). Here we show that if the luminous emission from a TDE is directly coupled to the rate of return of tidally destroyed debris to the SMBH, then the TDE luminosity function is in good agreement with observations and scales as ∝L^−2.5for high luminosities, provided that the SMBH mass function${\mathit{dN}}_{•}/{\mathit{dM}}_{•}$—the number of SMBHs (N_•) per SMBH mass (M_•)—is approximately flat in the mass range over which we observe TDEs. We also show that there is a cutoff in the luminosity function at low luminosities that is a result of direct captures, and this cutoff has been tentatively observed. If${\mathit{dN}}_{•}/{\mathit{dM}}_{•}$is flat, which is in agreement with some observational campaigns, these results suggest that the fallback rate feeds the accretion rate in TDEs. Contrarily, if${\mathit{dN}}_{•}/d\log M_{•}$is flat, which has been found theoretically and is suggested by other observational investigations, then the emission from TDEs is likely powered by another mechanism. Future observations and more TDE statistics, provided by the Rubin Observatory/LSST, will provide additional evidence as to the reality of this tension.

    

   more » « less