Tidal Disruption Events Rate Research Articles

The rate of extreme coronal line emitting galaxies in the Sloan Digital Sky Survey and their relation to tidal disruption events

ABSTRACT High-ionization iron coronal lines (CLs) are a rare phenomenon observed in galaxy and quasi-stellar object spectra that are thought to be created by high-energy emission from active galactic nuclei and certain types of transients. In cases known as extreme coronal line emitting galaxies (ECLEs), these CLs are strong and fade away on a time-scale of years. The most likely progenitors of these variable CLs are tidal disruption events (TDEs), which produce sufficient high-energy emission to create and sustain the CLs over these time-scales. To test the possible connection between ECLEs and TDEs, we present the most complete variable ECLE rate calculation to date and compare the results to TDE rates from the literature. To achieve this, we search for ECLEs in the Sloan Digital Sky Survey (SDSS). We detect sufficiently strong CLs in 16 galaxies, more than doubling the number previously found in SDSS. Using follow-up spectra from the Dark Energy Spectroscopic Instrument and Gemini Multi-Object Spectrograph, Wide-field Infrared Survey Explorer mid-infrared observations, and Liverpool Telescope optical photometry, we find that none of the nine new ECLEs evolve in a manner consistent with that of the five previously discovered variable ECLEs. Using this sample of five variable ECLEs, we calculate the galaxy-normalized rate of variable ECLEs in SDSS to be $R_\mathrm{G}=3.6~^{+2.6}_{-1.8}~(\mathrm{statistical})~^{+5.1}_{-0.0}~(\mathrm{systematic})\times 10^{-6}~\mathrm{galaxy}^{-1}~\mathrm{yr}^{-1}$. The mass-normalized rate is $R_\mathrm{M}=3.1~^{+2.3}_{-1.5}~(\mathrm{statistical})~^{+4.4}_{-0.0}~(\mathrm{systematic})\times 10^{-17}~\mathrm{M_\odot ^{-1}}~\mathrm{yr}^{-1}$ and the volumetric rate is $R_\mathrm{V}=7~^{+20}_{-5}~(\mathrm{statistical})~^{+10}_{-0.0}~(\mathrm{systematic})\times 10^{-9}~\mathrm{Mpc}^{-3}~\mathrm{yr}^{-1}$. Our rates are one to two orders of magnitude lower than TDE rates from the literature, which suggests that only 10–40 per cent of all TDEs produce variable ECLEs. Additional uncertainties in the rates arising from the structure of the interstellar medium have yet to be included.

Tidal Disruption of a Star on a Nearly Circular Orbit

We consider Roche lobe overflow (RLO) from a low-mass star on a nearly circular orbit, onto a supermassive black hole (SMBH). If mass transfer is unstable, its rate accelerates in a runaway process, resulting in highly super-Eddington mass accretion rates, accompanied by an optically thick outflow emanating from the SMBH vicinity. This produces a 1–4 week long, bright optical/ultraviolet flare, accompanied by a 1–10 year long X-ray precursor and post-cursor emitted from the accretion flow onto the SMBH. Such “Circular Tidal Disruption Events” (TDEs) represent a new class of nuclear transients, occurring at up to 1%–10% of the canonical parabolic tidal disruption event rate. Near-breakup rotation and strong tidal deformation of the star prior to disruption could lead to strong magnetic fields, making circular TDEs possible progenitors of jetted TDEs. Outflows prior to the final stellar disruption produce a circumnuclear environment (CNM) with ∼10−2 M⊙ at distances of ∼0.01–0.1 pc, likely leading to bright radio emission, and also similar to the CNM inferred for jetted TDEs. We discuss broader connections between circular TDEs and other recently identified classes of transients associated with galactic nuclei, such as repeating TDEs and quasi-periodic X-ray eruptions, as well as possible connections to luminous fast blue optical transients such as AT2018cow. We also discuss observational signatures of the analogous RLO of a white dwarf around an intermediate-mass BH, which may be a multimessenger source in the LISA era.

Counting the Unseen. I. Nuclear Density Scaling Relations for Nucleated Galaxies

The volumetric rate of tidal disruption events (TDEs) encodes information on the still-unknown demographics of central massive black holes (MBHs) in low-mass galaxies (≲109 M ⊙). Theoretical TDE rates from model galaxy samples can extract this information, but this requires accurately defining the nuclear stellar density structures. This region is typically dominated by nuclear star clusters (NSCs), which have been shown to increase TDE rates by orders of magnitude. Thus, we assemble the largest available sample of parsec-scale 3D density profiles that include NSC components. We deproject the point-spread-function-deconvolved surface-brightness profiles of 91 nearby galaxies of varying morphology and combine these with nuclear mass-to-light ratios estimated from measured colors or spectral synthesis to create 3D mass density profiles. We fit the inner 3D density profile to find the best-fit power-law density profile in each galaxy. We compile this information as a function of galaxy stellar mass to fit new empirical density scaling relations. These fits reveal positive correlations between galaxy stellar mass and central stellar density in both early- and late-type galaxies. We find that early-type galaxies have somewhat higher densities and shallower profiles relative to late-type galaxies at the same mass. We also use the density profiles to estimate the influence radius of each galaxy’s MBH and find that the sphere of influence was likely resolved in most cases. These new relations will be used in future works to build mock galaxy samples for dynamical TDE rate calculations, with the aim of constraining MBH demographics in low-mass galaxies.

Extended Emission-line Regions in Poststarburst Galaxies Hosting Tidal Disruption Events

We report the discovery of an extended emission-line region (EELR) in MUSE observations of Markarian 950, a nearby (z = 0.01628) poststarburst (PSB) galaxy that hosted the tidal disruption event (TDE) iPTF 16fnl. The EELR requires a nonstellar ionizing continuum with a luminosity Lion,min≳1043 erg s−1, inconsistent with the current weak state (L IR,AGN < 2.5 × 1042 erg s−1) of the galactic nucleus. The ionized gas has low velocity (∼–50 km s−1) and low turbulence (σ gas ≲ 50 km s−1) and is kinematically decoupled from the stellar motions, indicating that the gas kinematics is not active galactic nucleus (AGN) driven. Markarian 950 is the third PSB galaxy to host a weak nuclear ionizing source as well as an EELR and a TDE. The overall properties of these three galaxies, including the kinematics and accretion history, are unusual but strikingly similar. We estimate that the incidence of EELRs in PSB-TDE hosts is a factor of ∼10 × higher than in other PSB galaxies. This suggests that a gas-rich postmerger environment is a key ingredient in driving elevated TDE rates. Based on the current observations, we cannot rule out that the EELRs may be powered through an elevated TDE rate in these galaxies. If the EELRs are not TDE powered, the presence of intermittent AGN activity, and in particular the fading of the AGN, may be associated with an increased TDE rate and/or an increased rate of detecting TDEs.

Strong lensing of tidal disruption events: Detection rates in imaging surveys

Tidal disruption events (TDEs) are multi-messenger transients in which a star is tidally destroyed by a supermassive black hole at the center of galaxies. The Rubin Observatory Legacy Survey of Space and Time (LSST) is anticipated to detect hundreds to thousands of TDEs annually, such that the first gravitationally lensed TDE may be observed in the coming years. Using Monte-Carlo simulations, we quantify the rate of both unlensed and lensed TDEs as a function of limiting magnitudes in four different optical bands ($u$, $g$, $r$, and $i$) for a range of TDE temperatures that match observations. Dependent on the temperature and luminosity model, we find that $g$ and $r$ bands are the most promising bands with unlensed TDE detections that can be as high as $ sim $ annually. By populating a cosmic volume with realistic distributions of TDEs and galaxies that can act as gravitational lenses, we estimate that a few lensed TDEs (depending on the TDE luminosity model) can be detected annually in $g$ or $r$ bands in the LSST survey, with TDE redshifts in the range of $ sim 0.5$ to $ sim 2$. The ratio of lensed to unlensed detections indicates that we may detect $ sim 1$ lensed event for every $10^ $ unlensed events, which is independent of the luminosity model. The number of lensed TDEs decreases as a function of the image separations and time delays, and most of the lensed TDE systems are expected to have image separations below $ sim and time delays within sim 30$ days. At fainter limiting magnitudes, the $i$ band becomes notably more successful. These results suggest that strongly lensed TDEs are likely to be observed within the coming years and such detections will enable us to study the demographics of black holes at higher redshifts through the lensing magnifications. Our simulated catalogs of lensed TDEs are publicly available.

Repeating partial disruptions and two-body relaxation

Two-body relaxation may drive stars onto near-radial orbits around a massive black hole, resulting in a tidal disruption event (TDE). In some circumstances, stars are unlikely to undergo a single terminal disruption, but rather to have a sequence of many grazing encounters with the black hole. It has long been unclear what is the physical outcome of this sequence: each of these encounters can only liberate a small amount of stellar mass, but may significantly alter the orbit of the star. We study the phenomenon of repeating partial tidal disruptions (pTDEs) by building a semi-analytical model that accounts for mass loss and tidal excitations. In the empty loss cone regime, where two-body relaxation is weak, we estimate the number of consecutive partial disruption that a star can undergo, on average, before being significantly affected by two-body encounters. We find that in this empty loss cone regime, a star will be destroyed in a sequence of weak pTDEs, possibly explaining the tension between the low observed TDE rate and its higher theoretical estimates.

Exploring Low-mass Black Holes through Tidal Disruption Events in the Early Universe: Perspectives in the Era of the JWST, Roman Space Telescope, and LSST Surveys

The James Webb Space Telescope (JWST) has recently uncovered the presence of low-luminosity active galactic nuclei (AGNs) at z = 4–11. Spectroscopic observations have provided estimates of the nuclear black hole (BH) masses for these sources, extending the low-mass boundary down to M • ∼ 106–7 M ⊙. Despite this breakthrough, the observed lowest mass of BHs is still ≳1–2 orders of magnitude heavier than the predicted mass range of their seed population, thereby leaving the initial mass distribution of massive BHs poorly constrained. In this paper, we focus on UV-to-optical (in the rest frame) flares of stellar tidal disruption events (TDEs) embedded in low-luminosity AGNs as a tool for exploring low-mass BH populations with ≲104–6 M ⊙. We provide an estimate of the TDE rate over z = 4–11, associated with the properties of JWST-detected AGN host galaxies, and we find that deep and wide survey programs with JWST and the Roman Space Telescope (RST) can detect and identify TDEs up to z ≃ 4–7. The predicted detection numbers of TDEs at z > 4 in 1 yr are NTDE∼2–10(0.2–2) for the JADES-Medium (and COSMOS-Web) survey with JWST and NTDE∼2–10(8–50) for the deep (and wide) tiers of the High Latitude Time Domain Survey with RST. We further discuss survey strategies for hunting for transient high-redshift TDEs in wide-field surveys with RST, as well as a joint observation campaign with the Vera C. Rubin Observatory for enhancing the detection number. The high-redshift TDE search will give us a unique opportunity to probe the mass distribution of early BH populations.

A search for intermediate-mass black holes in compact stellar systems through optical emissions from tidal disruption events

ABSTRACT Intermediate-mass black holes (IMBHs) are expected to exist in globular clusters (GCs) and compact stellar systems (CSSs) in general, but none have been conclusively detected. Tidal disruption events (TDEs), where a star is tidally disrupted by the gravitational field of a black hole, have been observed to occur around the supermassive black holes (SMBHs) found at the centres of galaxies, and should also arise around IMBHs, especially in the dense stellar cores of CSSs. However, to date none have been observed in such environments. Using data from the Zwicky Transient Facility, we search for TDEs associated with CSSs, but none are found. This non-detection allows us to set an upper limit on the TDE rate in CSSs of $n_\text{TDE,Total}\lessapprox 10^{-7} \, \mathrm{ CSS}^{-1}\, \text{yr}^{-1}$, which is 2 dex below the observed TDE rate involving SMBHs interacting with 1 M⊙ main-sequence stars in the nuclei of massive galaxies. We also consider ultracompact dwarfs (UCDs) formed through a tidal stripping process in the surveyed volume. On the assumption that these CSSs contain SMBHs and TDE rates are comparable to current observed optical rates in galactic nuclei ($\approx 3.2 \times 10^{-5}\, \text{gal}^{-1}\, \text{yr}^{-1}$), we determine an upper limit for the number of UCDs formed through a tidal stripping process in the surveyed volume to be NGC,Strip &amp;lt; 1.4 × 104, which we estimate represents $\lt 6~{{\ \rm per\ cent}}$ of the population of GCs &amp;gt;106 M⊙.

Cosmic hide and seek: The volumetric rate of X-ray quasi-periodic eruptions

Multiwavelength extragalactic nuclear transients, particularly those detectable as multi-messengers, are among the primary drivers for the next-generation observatories. X-ray quasi-periodic eruptions (QPEs) are the most recent and perhaps most peculiar addition to this group. Here, we report a first estimate of the volumetric rate of QPEs based on the first four discoveries with the eROSITA X-ray telescope onboard the Spectrum Roentgen Gamma observatory. Under the assumption, supported by a suite of simulated light curves, that these four sources sample the intrinsic population somewhat homogeneously, we correct for their detection efficiency and compute a QPE abundance of Rvol = 0.60−0.43+4.73 × 10−6 Mpc−3 above an intrinsic average log L0.5−2.0 keVpeak &gt; 41.7. Since the exact lifetime of QPEs (τlife) is currently not better defined than between a few years or few decades, we convert this to a formation rate of ℛvol/τlife ≈ 0.6 × 10−7(τlife/10 yr)−1 Mpc−3 yr−1. As a comparison, this value is a factor ∼10 τlife times smaller than the formation rate of tidal disruption events. The origin of QPEs is still debated, although lately most models suggest that they are the electromagnetic counterpart of extreme mass ratio inspirals (EMRIs). In this scenario, the QPE rate would thus be the first ever constraint (i.e. a lower limit) to the EMRI rate from observations alone. Future discoveries of QPEs and advances in their theoretical modeling will consolidate or rule out their use for constraining the number of EMRIs detectable by the LISA mission.

Changing-look Active Galactic Nuclei Behavior Induced by Disk-captured Tidal Disruption Events

Recent observations of changing-look active galactic nuclei (AGNs) hint at a frequency of accretion activity not fully explained by tidal disruption events (TDEs) stemming from relaxation processes in nuclear star clusters (NSCs), traditionally estimated to occur at rates of 10−4–10−5 yr−1 per galaxy. In this Letter, we propose an enhanced TDE rate through the AGN disk capture process, presenting a viable explanation for the frequent transitions observed in changing-look AGNs. Specifically, we investigate the interaction between the accretion disk and retrograde stars within NSCs, resulting in the rapid occurrence of TDEs within a condensed time frame. Through detailed calculations, we derive the time-dependent TDE rates for both relaxation-induced TDE and disk-captured TDE. Our analysis reveals that TDEs triggered by the disk capture process can notably amplify the TDE rate by several orders of magnitude during the AGN phase. This mechanism offers a potential explanation for the enhanced high-energy variability characteristic of changing-look AGNs.

A New Population of Mid-infrared-selected Tidal Disruption Events: Implications for Tidal Disruption Event Rates and Host Galaxy Properties

Most tidal disruption events (TDEs) are currently found in time-domain optical and soft X-ray surveys, both of which are prone to significant obscuration. The infrared (IR), however, is a powerful probe of dust-enshrouded environments; hence, we recently performed a systematic search of NEOWISE mid-IR data for nearby, obscured TDEs within roughly 200 Mpc. We identified 18 TDE candidates in galactic nuclei, using difference imaging to uncover nuclear variability among significant host galaxy emission. These candidates were selected based on the following IR light-curve properties: (1) L W2 ≳ 1042 erg s−1 at peak; (2) fast rise, followed by a slow, monotonic decline; (3) no significant prior variability; and (4) no evidence for active galactic nucleus (AGN) activity in Wide-field Infrared Survey Explorer (WISE) colors. The majority of these sources showed no variable optical counterpart, suggesting that optical surveys indeed miss numerous obscured TDEs. Using narrow-line ionization levels and variability arguments, we identified six sources as possible underlying AGN, yielding a total of 12 TDEs in our gold sample. This gold sample yields a lower limit on the IR-selected TDE rate of (2.0 ± 0.3) × 10−5 galaxy−1 yr−1 ((1.3 ± 0.2) × 10−7 Mpc−3 yr−1), which is comparable to optical and X-ray TDE rates. The IR-selected TDE host galaxies do not show a green valley overdensity nor as a preference for quiescent, Balmer strong galaxies, which are both overrepresented in optical and X-ray TDE samples. This IR-selected sample represents a new population of dusty TDEs that have historically been missed by optical and X-ray surveys and helps alleviate tensions between observed and theoretical TDE rates and the so-called missing energy problem.

An Explanation for the Overrepresentation of Tidal Disruption Events in Post-starburst Galaxies

Tidal disruption events (TDEs) provide a valuable probe for studying the dynamics of stars in the nuclear environments of galaxies. Recent observations show that TDEs are strongly overrepresented in post-starburst or “green valley” galaxies, although the underlying physical mechanism remains unclear. Considering the possible interaction between stars and active galactic nucleus (AGN) disks, the TDE rates can be greatly changed compared to those in quiescent galactic nuclei. In this work, we revisit TDE rates by incorporating an evolving AGN disk within the framework of “loss cone” theory. We numerically evolve the Fokker–Planck equations by considering star–disk interactions, in situ star formation in the unstable region of the outer AGN disk, and evolution of the accretion process of supermassive black holes. We find that the TDE rates are enhanced by about 2 orders of magnitude shortly after the AGN transitions into an inactive stage. During this phase, the accumulated stars rapidly scatter into the loss cone due to the disappearance of the inner standard thin disk. Our results provide an explanation for the overrepresentation of TDEs in post-starburst galaxies.

Tidal Disruption Events from the Combined Effects of Two-body Relaxation and the Eccentric Kozai–Lidov Mechanism

Tidal disruption events (TDEs) take place when a star ventures too close to a supermassive black hole (SMBH) and becomes ruptured. One of the leading proposed physical mechanisms often invoked in the literature involves weak two-body interactions experienced by the population of stars within the host SMBH’s sphere of influence, commonly referred to as two-body relaxation. This process can alter the angular momentum of stars at large distances and place them into nearly radial orbits, thus driving them to disruption. On the other hand, gravitational perturbations from an SMBH companion via the eccentric Kozai–Lidov (EKL) mechanism have also been proposed as a promising stellar disruption channel. Here we demonstrate that the combination of EKL and two-body relaxation in SMBH binaries is imperative for building a comprehensive picture of the rates of TDEs. Here we explore how the density profile of the surrounding stellar distribution and the binary orbital parameters of an SMBH companion influence the rate of TDEs. We show that this combined channel naturally produces disruptions at a rate that is consistent with observations and also naturally forms repeated TDEs, where a bound star is partially disrupted over multiple orbits. Recent observations show stars being disrupted in short-period orbits, which is challenging to explain when these mechanisms are considered independently. However, the diffusive effect of two-body relaxation, combined with the secular nature of the eccentricity excitations from EKL, is found to drive stars on short eccentric orbits at a much higher rate. Finally, we predict that rTDEs are more likely to take place in the presence of a steep stellar density distribution.

Uncovering Hidden Massive Black Hole Companions with Tidal Disruption Events

Dynamical perturbations from supermassive black hole (SMBH) binaries can increase the rates of tidal disruption events (TDEs). However, most previous work focuses on TDEs from the heavier black hole in the SMBH binary (SMBHB) system. In this work, we focus on the lighter black holes in SMBHB systems and show that they can experience a similarly dramatic increase in their TDE rate due to perturbations from a more massive companion. While the increase in TDEs around the more massive black hole is mostly due to chaotic orbital perturbations, we find that, around the smaller black hole, the eccentric Kozai–Lidov mechanism is dominant and capable of producing a comparably large number of TDEs. In this scenario, the mass derived from the light curve and spectra of TDEs caused by the lighter SMBH companion is expected to be significantly smaller than the SMBH mass estimated from galaxy scaling relations, which are dominated by the more massive companion. This apparent inconsistency can help find SMBHB candidates that are not currently accreting as active galactic nuclei and that are at separations too small for them to be resolved as two distinct sources. In the most extreme cases, these TDEs provide us with the exciting opportunity to study SMBHBs in galaxies where the primary SMBH is too massive to disrupt Sun-like stars.

Tidal disruption rate suppression by the event horizon of spinning black holes

ABSTRACT The rate of observable tidal disruption events (TDEs) by the most massive (≳ few × 107 M⊙) black holes (BHs) is suppressed due to direct capture of stars. This suppression effect depends on the space–time geometry and holds the promise of probing the spin distribution of dormant BHs at the centres of galaxies. By extending the frozen-in approximation commonly used in the Newtonian limit, we propose a general relativistic criterion for the tidal disruption of a star of given interior structure. The rate suppression factor is then calculated for different BH masses, spins, and realistic stellar populations. We find that either a high-BH spin (≳ 0.5) or a young stellar population (≲1 Gyr) allows TDEs to be observed from BHs significantly more massive than $10^8\, \mathrm{M}_{\odot }$. We call this spin-age degeneracy (SAD). This limits our utility of the TDE rate to constrain the BH spin distribution, unless additional constraints on the age of the stellar population or the mass of the disrupted star can be obtained by modelling the TDE radiation or the stellar spectral energy distribution near the galactic nuclei.

Integral Field Spectroscopy of 13 Tidal Disruption Event Hosts from the Zwicky Transient Facility Survey

The host galaxies of tidal disruption events (TDEs) have been shown to possess peculiar properties, including high central light concentrations, unusual star formation histories, and “green” colors. The ubiquity of these large-scale galaxy characteristics among TDE host populations suggests that they may serve to boost the TDE rate in such galaxies by influencing the nuclear stellar dynamics. We present the first population study of integral field spectroscopy for 13 TDE host galaxies across all spectral classes and X-ray brightnesses with the purpose of investigating their large-scale properties. We derive the black hole masses via stellar kinematics (i.e., the M–σ relation) and find masses in the range , with a distribution dominated by black holes with M BH ∼ 106 M ⊙. We find one object with M BH ≳ 108 M ⊙, above the “Hills mass”, which if the disrupted star was of solar type, allows a lower limit of a ≳ 0.16 to be placed on its spin, lending further support to the proposed connection between featureless TDEs and jetted TDEs. We also explore the level of rotational support in the TDE hosts, quantified by (V/σ) e , a parameter that has been shown to correlate with the stellar age and may explain the peculiar host-galaxy preferences of TDEs. We find that the TDE hosts exhibit a broad range in (V/σ) e following a similar distribution as E + A galaxies, which have been shown to be overrepresented among TDE host populations.

Loss cone shielding

ABSTRACT A star wandering close enough to a massive black hole can be ripped apart by the tidal forces of the black hole. The advent of wide-field surveys at many wavelengths has quickly increased the number of tidal disruption events (TDEs) observed, and has revealed that (i) observed TDE rates are lower than theoretical predictions and (ii) E+A galaxies are significantly overrepresented. This overrepresentation further worsens the tension between observed and theoretically predicted TDEs for non-E+A galaxies. Classical loss cone theory focuses on the cumulative effect of many weak scatterings. However, a strong scattering can remove a star from the distribution before it can get tidally disrupted. Most stars undergoing TDEs come from within the radius of influence, the densest environments of the Universe. In such environments, close encounters rare elsewhere become non-negligible. We revise the standard loss cone theory to take into account classical two-body interactions as well as strong scattering, collisions, tidal captures, and study under which conditions close encounters can shield the loss cone. We (i) analytically derive the impact of strong scattering and other close encounters, (ii) compute time-dependent loss cone dynamics including both weak and strong encounters, and (iii) derive analytical solutions to the Fokker–Planck equation with strong scattering. We find that (i) TDE rates can be reduced to up to an order of magnitude and (ii) strong shielding preferentially reduces deeply plunging stars. We also show that stellar overdensities, one possible explanation for the E + A preference, can fail to increase TDE rates when taking into account strong scattering.

Tidal Disruption Event Demographics with the Zwicky Transient Facility: Volumetric Rates, Luminosity Function, and Implications for the Local Black Hole Mass Function

We conduct a systematic tidal disruption event (TDE) demographics analysis using the largest sample of optically selected TDEs. A flux-limited, spectroscopically complete sample of 33 TDEs is constructed using the Zwicky Transient Facility over 3 yr (from 2018 October to 2021 September). We infer the black hole (BH) mass (M BH) with host galaxy scaling relations, showing that the sample M BH ranges from 105.1 M ⊙ to 108.2 M ⊙. We developed a survey efficiency corrected maximum volume method to infer the rates. The rest-frame g-band luminosity function can be well described by a broken power law of , with L bk = 1043.1 erg s−1. In the BH mass regime of 105.3 ≲ (M BH/M ⊙) ≲ 107.3, the TDE mass function follows , which favors a flat local BH mass function (). We confirm the significant rate suppression at the high-mass end (M BH ≳ 107.5 M ⊙), which is consistent with theoretical predictions considering direct capture of hydrogen-burning stars by the event horizon. At a host galaxy mass of M gal ∼ 1010 M ⊙, the average optical TDE rate is ≈3.2 × 10−5 galaxy−1 yr−1. We constrain the optical TDE rate to be [3.7, 7.4, and 1.6] × 10−5 galaxy−1 yr−1 in galaxies with red, green, and blue colors.

The Correlations of Stellar Tidal Disruption Rates with Properties of Massive Black Holes and Their Host Galaxies

Stars can be either disrupted as tidal disruption events (TDEs) or swallowed whole by massive black holes (MBHs) at galactic centers when they approach sufficiently close to these MBHs. In this work, we investigate the correlations of such stellar consumption rates with both the MBH mass M BH and the inner slope of the host-galaxy mass density distribution α. We introduce a simplified analytical power-law model with a power-law stellar-mass density distribution surrounding MBHs and separate the contributions of two-body relaxation and stellar orbital precession for the stellar orbital angular momentum evolution in nonspherical galaxy potentials. The stellar consumption rates derived from this simplified model can be well consistent with the numerical results obtained with a more realistic treatment of stellar distributions and dynamics around MBHs, providing an efficient way to estimate TDE rates. The origin of the correlations of stellar consumption rates with M BH and α is explained by the dependence of this analytical model on those MBH/host-galaxy properties and by the separation of the stellar angular momentum evolution mechanisms. We propose that the strong positive correlation between the rates of stellar consumption due to two-body relaxation and α provides one interpretation for the overrepresentation of TDEs found in some rare E+A/poststarburst galaxies. We find high TDE rates for giant stars, up to those for solar-type stars. Understanding the origin of the correlations of the stellar consumption rates will be necessary for obtaining the demographics of MBHs and their host galaxies via TDEs.

The Luminosity Function of Tidal Disruption Events from Fallback-powered Emission: Implications for the Black Hole Mass Function

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, , 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.5 for high luminosities, provided that the SMBH mass function —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 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 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.