skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Cooling Envelope Model for Tidal Disruption Events
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 radiusRv∼ 1014cm (photosphere radius ∼1015cm) surrounding the supermassive black hole (SMBH). As the envelope cools radiatively, it undergoes Kelvin–Helmholtz contractionRv∝t−1, its temperature risingTeff∝t1/2while its total luminosity remains roughly constant; the optical luminosity decays as ν L ν R v 2 T eff t 3 / 2 . Despite this similarity to the mass fallback rate M ̇ fb 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
Award ID(s):
2009255
PAR ID:
10371935
Author(s) / Creator(s):
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal Letters
Volume:
937
Issue:
1
ISSN:
2041-8205
Format(s):
Medium: X Size: Article No. L12
Size(s):
Article No. L12
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract 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,tpeak, on both stellar mass and age, with 15 days ≲tpeak≲ 10 yr for main sequence stars with masses 0.2 ≤M/M≤ 5 andM= 106M. The second estimate, which postulates that the star is completely destroyed when tides dominate the maximum stellar self-gravity, predicts thattpeakis very weakly dependent on stellar type, with t peak = 23.2 ± 4.0 days M / 10 6 M 1 / 2 for 0.2 ≤M/M≤ 5, while t peak = 29.8 ± 3.6 days M / 10 6 M 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
  2. 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 ofRB≈ 2 × 105GM/c2, whereMis 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−1rather thanr−3/2and the accretion rate M ̇ is consequently suppressed by over 2 orders of magnitude below the Bondi rate M ̇ B . We find continuous energy feedback from the accretion flow to the external medium at a level of 10 2 M ̇ c 2 5 × 10 5 M ̇ 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. Abstract 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 N ̇ TDE / 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 dN / 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 dN / 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 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
  4. Abstract State transitions in black hole X-ray binaries are likely caused by gas evaporation from a thin accretion disk into a hot corona. We present a height-integrated version of this process, which is suitable for analytical and numerical studies. With radiusrscaled to Schwarzschild units and coronal mass accretion rate m ̇ c to Eddington units, the results of the model are independent of black hole mass. State transitions should thus be similar in X-ray binaries and an active galactic nucleus. The corona solution consists of two power-law segments separated at a break radiusrb∼ 103(α/0.3)−2, whereαis the viscosity parameter. Gas evaporates from the disk to the corona forr>rb, and condenses back forr<rb. Atrb, m ̇ c reaches its maximum, m ̇ c , max 0.02 ( α / 0.3 ) 3 . If atr≫rbthe thin disk accretes with m ̇ d < m ̇ c , max , then the disk evaporates fully before reachingrb, giving the hard state. Otherwise, the disk survives at all radii, giving the thermal state. While the basic model considers only bremsstrahlung cooling and viscous heating, we also discuss a more realistic model that includes Compton cooling and direct coronal heating by energy transport from the disk. Solutions are again independent of black hole mass, andrbremains unchanged. This model predicts strong coronal winds forr>rb, and aT∼ 5 × 108K Compton-cooled corona forr<rb. Two-temperature effects are ignored, but may be important at small radii. 
    more » « less
  5. Abstract We present13CO(J= 1 → 0) observations for the EDGE-CALIFA survey, which is a mapping survey of 126 nearby galaxies at a typical spatial resolution of 1.5 kpc. Using detected12CO emission as a prior, we detect13CO in 41 galaxies via integrated line flux over the entire galaxy and in 30 galaxies via integrated line intensity in resolved synthesized beams. Incorporating our CO observations and optical IFU spectroscopy, we perform a systematic comparison between the line ratio 12 / 13 I [ 12 CO ( J = 1 0 ) ] / I [ 13 CO ( J = 1 0 ) ] and the properties of the stars and ionized gas. Higher 12 / 13 values are found in interacting galaxies compared to those in noninteracting galaxies. The global 12 / 13 slightly increases with infrared colorF60/F100but appears insensitive to other host-galaxy properties such as morphology, stellar mass, or galaxy size. We also present azimuthally averaged 12 / 13 profiles for our sample up to a galactocentric radius of 0.4r25(∼6 kpc), taking into account the13CO nondetections by spectral stacking. The radial profiles of 12 / 13 are quite flat across our sample. Within galactocentric distances of 0.2r25, the azimuthally averaged 12 / 13 increases with the star formation rate. However, Spearman rank correlation tests show the azimuthally averaged 12 / 13 does not strongly correlate with any other gas or stellar properties in general, especially beyond 0.2r25from the galaxy centers. Our findings suggest that in the complex environments in galaxy disks, 12 / 13 is not a sensitive tracer for ISM properties. Dynamical disturbances, like galaxy interactions or the presence of a bar, also have an overall impact on 12 / 13 , which further complicates the interpretations of 12 / 13 variations. 
    more » « less