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.


INTRODUCTION
A tidal disruption event (TDE) occurs when a star passes close enough to a supermassive black hole (SMBH) with mass M BH ≲ 10 8 M ⊙ .If the star is disrupted by the tidal force, half of the stellar debris will become gravitationally bound and fall back towards the SMBH, which will trigger a short-lived flare of emission (e.g., Rees 1988;Evans & Kochanek 1989;Phinney 1989;Ulmer 1999).The matter fallback rate and the light curve theoretically follow t −5/3 over the timescale of months to years, which is a typical feature for TDEs (Rees 1988).In the past few years, more than 100 TDE candidates have been reported (e.g., French et al. 2020b;Hammerstein et al. 2023;Yao et al. 2023) and the number will dramatically increase due to current and upcoming multi-band time-domain surveys (Rau et al. 2009;Cham-bers et al. 2016;Shappee et al. 2014;Bellm et al. 2019;Ivezić et al. 2019;Merloni et al. 2012;Wang et al. 2023b).
With the increasing number of observed TDEs, it is possible to estimate the TDE rates, which strongly depend on the environment surrounding the central SMBH (typically ∼ 0.1 − 10 pc for galaxies with M BH ∼ 10 6−8 M ⊙ , see Stone et al. 2020;French et al. 2020b, for reviews), including the number density distribution and the orbital dynamics of stars.Based on the different TDE samples, the TDE rates are typically around 10 −5 − 10 −4 yr −1 gal −1 (e.g., Magorrian & Tremaine 1999;Wang & Merritt 2004;Arcavi et al. 2014;French et al. 2016;van Velzen et al. 2018;Stone et al. 2020;Yao et al. 2023), which is in agreement with that estimated from two-body relaxation (see Merritt 2013;Stone et al. 2020, for reviews).However, several recent studies show that TDE hosts are significantly overrepresented in the post-starburst or E+A galaxies, whose spectra are characterized by a lack of strong emission lines but with strong Balmer absorption features (French et al. 2017;Law-Smith et al. 2017;French et al. 2020b;Hammerstein et al. 2021Hammerstein et al. , 2023;;Yao et al. 2023).Such post-starburst spectra indicate a low current star formation rate and a recent (≲ 1 Gyr) star formation burst (but just ended, see Law-Smith et al. 2017;French et al. 2020b).These Balmer-strong galaxies (including post-starburst or E+A galaxies) are commonly found in the green valley, which represents a transitional population between blue star-forming galaxies and red-quenched galaxies (see Salim 2014;Schawinski et al. 2014).Accounting for possible selection effects, the TDE rate is estimated to be ∼ 25 − 50 times overrepresented in this subclass of galaxies (Law-Smith et al. 2017;Hammerstein et al. 2021), suggesting a close correlation between TDE occurrence and galactic evolution.
Several explanations were proposed to account for the overenhancement in post-starburst galaxies.The gas involved in the galaxy merger loses a substantial fraction of its specific angular momentum due to tidal torques and dynamical friction processes, falling quickly into the center.Therefore, the galaxies may have an overdense central star cluster during this stage, which leads to a shorter timescale of two-body relaxation (Stone & Metzger 2016;Stone & van Velzen 2016;French et al. 2020a;Chen et al. 2023).Indeed, high-resolution SDSS and HST imaging showed that the TDE host galaxies have higher central concentrations on pc-kpc scales (Law-Smith et al. 2017;Graur et al. 2018;French et al. 2020a).The second possible channel involves the presence of an SMBH binary.Gravitational perturbations from an SMBH companion can significantly modify the orbits of surrounding stars (Ivanov et al. 2005;Chen et al. 2011;Melchor et al. 2023).The TDE rates could be 2-4 orders of magnitude larger than that estimated for the case of a single SMBH.However, the duration of the phase of enhanced tidal disruption is quite short, which is only ∼ 10 5 yr (Ivanov et al. 2005).It was also proposed that after the merger of the SMBH binary, the gravitational recoil imparted to the merged SMBH fills its loss cone and results in a TDE rate as high as ∼ 0.1 yr −1 (Stone & Loeb 2011).There are also several other dynamical effects for the post-starburst preference, which include radial anisotropies (Stone et al. 2018), nuclear triaxiality (Merritt & Poon 2004) and eccentric nuclear stellar disk (Madigan et al. 2018), where the secular effects can dramatically increase TDE rates.However, the nuclear cluster mass should be relatively small in these scenarios, so as not to give rise to mass precession that quenches coherent secular evolution (see Stone et al. 2020, for a review).
In this work, we propose an explanation for the overrepresentation of TDEs in post-starburst galaxies by considering the possible interaction between the nuclear stars and AGN disk, where the variation in nuclear stellar density distribution and orbital dynamics within the parsec scale will strongly affect the TDE rates.To achieve high TDE rates, both a high number density and the presence of high eccentric orbits of stars (around the loss cone boundary) are nec-essary.Interestingly, observations indicate that active galactic nuclei (AGNs) prefer to reside in hosts that have experienced a recent starburst (∼ a few×10 8 yr ago, see Davies et al. 2007;Schawinski et al. 2009;Wild et al. 2010;Kaviraj et al. 2011;Hopkins 2012).In these starburst galaxies, which could be triggered by galactic mergers (Hopkins & Quataert 2010;Volonteri et al. 2015;Liao et al. 2023), both the number density distribution and orbital dynamics of stars could be affected by the starburst and the subsequent AGN activity.Thus, the TDE rates are naturally connected to the history of galactic evolution.We calculate TDE rates by solving the Fokker-Planck (FP) equations taking into consideration the influence of the accretion disk, which includes a stable inner and a star-forming unstable outer disk surrounding SMBHs (Pan & Yang 2021;Wang et al. 2023a).The insitu star formation in the disk (Wang et al. 2021;Fan & Wu 2023) and the inward migration resulting from star-disk interactions lead to the accumulation of stars with low eccentricity in the nuclear region of galaxies, due to orbital circularization caused by density wave (e.g., Tanaka & Ward 2004).Once the disk dissipates, the stars are re-scattered towards the loss cone boundary, resulting in a rapid increase in TDE rates.
This paper is organized as follows: we present our main assumptions and numerical method in Section 2, and show the evolution of TDE rates in Section 3. We summarize and discuss our results in Section 4.

MODEL
The TDE rates can be calculated by integrating the flux in the direction of angular momentum at the boundary of the loss cone, which strongly depends on the environment surrounding the central SMBH and varies significantly during different stages of galactic evolution.The presence or absence of the accretion disk will significantly affect star formation, stellar evolution, and stellar orbits in the nuclear region.However, the precise evolution of the accretion process in the galactic center remains unclear.Following the positive correlation between the accretion rate of the central SMBH and the star formation rate (SFR) of galaxies given by simulations and observations, we simply assume a fast rise and slow decay for the evolution of accretion rate (Behroozi et al. 2013;Carnall et al. 2018), which is described by a doublepower-law form: where a and b are the falling and rising slopes, respectively, and τ is related to the time at which the accretion rate peaks.As shown in Figure 1, we divide the galaxy into three stages according to the accretion rate: the active stage, the transitional stage and the quiescent stage.
2.1.The AGN disk and the star-disk interactions The accretion disk is normally optically thick and geometrically thin in bright AGNs (Shakura & Sunyaev 1973).The standard thin disk is generally gravitationally unstable, leading to the collapse of self-gravitating clouds and forming stars at distances greater than several thousand Schwarzschild radii (Toomre 1964).In this work, we adopt the TQM disk in the active stage (see Thompson et al. 2005), which consists of an inner standard thin disk and an outer disk undergoing star formation.The star formation rate in the unstable region of the TQM disk, Ṁ⋆ , can be estimated from the well-known Kennicutt-Schmidt law ( Σ Ṁ⋆ ∼ 0.017Σ gas Ω gas , see Kennicutt 1998).In this work, we follow Dittmann & Miller (2020) to obtain the Ṁ⋆ , which is determined by the viscosity parameter α and the efficiency parameter ε (i.e., the rest mass energy from star formation is converted into radiation, see Dittmann & Miller 2020, for in detail).
In low-luminosity AGNs, the standard thin disk will switch into an optically thin, geometrically thick advectiondominated accretion flow (ADAF) when the accretion rate ṁ (≡ Ṁ/ ṀEdd ) falls below the critical value ṁcrit,1 (typically ṁcrit,1 ∼ 0.01, see Narayan & Yi 1994).The physical mechanism driving the disk transition is still unclear, but the evaporation/condensation processes may play a key role (e.g., Qiao & Liu 2009).In this work, we adopt a truncation radius of r tr = 17.3 ṁ−0.886 α 0.07 (in the unit of Schwarzschild radius r s , Taam et al. 2012) to represent the transition between these two accretion modes.Once the standard thin disk transitions to an ADAF, we neglect star-disk interactions due to the low gas density in the hot ADAF region.We assume that the standard disk is completely converted to the ADAF when the accretion rate ṁ < ṁcrit,2 ∼ 0.001.Since the truncation radius r tr ∼ 10 4 r s at this critical accretion rate.
In the AGN disk, the migration timescale and eccentricity damping timescale caused by the density wave excitation are well described by the following expression (Cresswell & Nelson 2008): , τ e = τ wave 0.78 1 − 0.14(e/h) 2 + 0.06(e/h) 3 (2) with where m s is the mass of stars, M = M(< a) is the total mass consisting of the SMBH, stars and disk gas within the radius a, e, a and Ω are the orbital eccentricity, semi-major axis and angular velocity of stars, Σ and h = H/a are the disk surface density and the disk aspect ratio, respectively.

FP equations and TDE rates
In this section, we introduce the orbit-averaged FP equations and the numerical set-up, which is mainly based on the work of Cohn & Kulsrud (1978); Stone & Metzger (2016); Pan & Yang (2021); Broggi et al. (2022); Wang et al. (2023a), where the numerical method for the FP equations can be found in Pan & Yang (2021); Broggi et al. (2022); Wang et al. (2023a).We consider a central SMBH with mass M BH surrounded by a spherically stellar cluster with total mass M s and a star-forming disk.The stars in the stellar cluster are assumed to have a single mass component m s .The initial stars' distribution follows Tremaine's cluster model (Tremaine et al. 1994), which is given by where r a = 4GM BH /σ 2 * ≡ 4r h is the radius of density transition (r h is the influential radius of the SMBH with mass M BH ), σ * is the stellar velocity dispersion given by M BH − σ ⋆ relaxation (Tremaine et al. 2002;Gültekin et al. 2009), (5) and γ ∼ 1.2 − 1.8 is the density scaling power index (Tremaine et al. 1994;Binney & Tremaine 2008).With the above density profiles, the collective gravitational potential background can be derived as: Following Cohn & Kulsrud (1978), we define the specific orbit energy and dimensionless specific orbital angular momentum: where J is the specific orbital angular momentum and J c (E) is the specific orbital angular momentum of a star with specific energy E on a circular orbit.In this case, the initial distribution function in the (E, R)-phase space is given by where n s is the number density of stars (Tremaine et al. 1994;Binney & Tremaine 2008).
Given initial distributions of stars, f i (t = 0, E, R), their evolution is governed by the orbit-averaged Fokker-Planck equation (Cohn & Kulsrud 1978) where with the diffusion coefficients {D EE , D RR , D ER } and the advection coefficients {D E , D R } as functions of f i (t, E, R) (which is defined in Cohn & Kulsrud 1978).We evolve the f i (t = 0, E, R) to a stable state f i (t = t f , E, R) according to the above FP equations and set the stable state as the initial state for subsequent calculations.
To calculate the TDE rates, we first obtain the evolution of the stars' distribution function by solving the FP equations above the three stages of galaxies: In the active stage, the stars can be divided into two populations: the disk component (stars embedded within the TQM disk) and the cluster component (stars stay outside of the TQM disk).Initially, we set the distribution function of cluster-component stars ).The evolution of cluster-and disk-component distribution functions is governed by the orbit-averaged FP equations (Pan & Yang 2021;Wang et al. 2023a) where with X = { f , g}, the diffusion coefficients {D EE , D RR , D ER }, and the advection coefficients {D X E , D X R } incorporating twobody scatterings and star-disk interactions.The advection coefficients in the disk are modified as where D E and D R are the advection terms in the quiescent stage and τ mig , τ e are given by equation ( 2) .It should be noted that the star-disk interactions will dominate the evolution of distribution function in the cluster only for the low inclination orbiters (θ ≲ π/6) since the migration torque by density wave decays quickly as the orbital inclination increase (e.g., Cresswell & Nelson 2008).Thus, we ignore the effect of the star-disk interactions on the clustercomponent stars, which means that the rate of TDEs in the active stage could be close to that in the quiescent stage (see also, MacLeod & Lin 2020; Teboul et al. 2022).
The source term S cap arises from stars captured by the disk and S dSF arises from in-situ star formation in the outer unsta-ble region of the disk, with where µ cap = 0.01 is a parameter quantifying the disk capture efficiency (see Pan & Yang 2021), τ mig is the migration timescale by star-disk interactions, Ṁ⋆ (E, R) = Ṁ⋆ (E) δ (R − 1) and the Ṁ⋆ (E) is obtained from the TQM disk model (Dittmann & Miller 2020).Here, we simply assume that the stars formed in the TQM disk are all in circular orbits, i.e., R = 1 (Derdzinski & Mayer 2022).For numerical simplification, we adopt the phenomenological treatment (S cap ) for the disk capture process rather than calculate the change in orbital inclination (θ ) of the stars' distribution function by solving the FP equations with the diffusion (D θ θ ) and advection (D θ ) terms (see Pan et al. 2022).
In the transitional stage, the inner thin disk will be gradually evaporated within a truncation radius (r tr ), and we neglect the star-disk interactions for the region within r < r tr .When the standard accretion disk disappears completely, the galaxy will enter into the quiescent stage, where the advection term D g E and D g R are substituted by D E and D R in the equation ( 12) (see Cohn & Kulsrud 1978;Pan & Yang 2021;Broggi et al. 2022;Wang et al. 2023a).Here, we do not include the process of the scattering of the disk-component stars into the cluster after the disk disappears.This process can be described by the diffusion (D θ θ ) and advection coefficients (D θ ) in the FP equations.However, both the diffusion coefficients D Eθ and D Rθ are equal to zero (see the appendix B in Pan et al. 2022, and reference therein), which means that the values of energy and angular momentum do not change as the stars' inclination changes in this process.Therefore, the D Eθ , D Rθ will not affect the TDE rates.
In the active stage, we set F R = 0 in the loss cone boundary due to the short eccentricity decay timescale of star-disk interactions (τ wave /τ disk ∼ 10 −5 , where τ disk is the duty cycle of AGN, see Pan & Yang 2021).Furthermore, the distribution function of stars in the quiescent stage is set as the initial conditions of the active stage in our simulation.In the transitional stage with intermediate accretion rate, we adopt F R = 0 for r > r tr and the standard loss cone boundary conditions for r ≤ r tr (see Cohn & Kulsrud 1978).In the quiescent stage, we still adopt the standard loss cone boundary conditions.The TDE rates are determined by the flux in the R-direction: 3. RESULTS Following the above procedure, we consider a model consisting of stars orbiting around an SMBH, where part of the stars are embedded in the accretion disk while others still stay in a spherical star cluster.For simplicity, we set a monochromatic star mass distribution m s = M ⊙ .We assume that the density power-law index of γ = 1.5 (ρ ∝ r −γ , see Pan & Yang 2021;Wang et al. 2023a, for details).In this work, we assume Ṁ0 = ṀEdd , a = 1.7, b = 4.0, τ = 2 × 10 6 yr, which will not affect the main result of overrepresentation of TDEs in the transitional stage.

The simulation results of the FP equations
In Figure 2, we show the stars' distribution function f (t, E, R) with of the different evolution stage in (E, R) space in the unit of pc −3 • σ −3 with the SMBH mass M BH = 10 6.75 M ⊙ .In Figure 2(a), we show the steady state of the evolution in the quiescent stage as the initial conditions of the active stage, where we assume that the accretion disk will present in some stage of galaxy evolution or duty-cycle effect.The distribution functions of active, transitional and quiescent stages are shown in Figure 2(b), (c) and (d), respectively.In the active stage, star-disk interactions (e.g., migration and damping of eccentricity) and star formation in the disk will dominate the evolution of the stars' distribution.As shown in Figure 2(b) with the accretion rate ṁ = 0.06, stars will gather at the boundary of R = 1 in the (E, R) space due to the circularizing induced by the density wave excitation.The stars in the circular orbits are far from the loss cone boundary, which hardly contributes to the TDE rates.In the transitional stage shown in Figure 2(c) with accretion rate ṁ = 0.009, the stars will be excited to the high eccentricity under the twobody scattering and enter the loss cone boundary as the stardisk interactions and star formation gradually vanish.In the quiescent stage, the distribution of stars will gradually return to the initial condition (see Figure 2d).

Evolution of the TDE rates
We present the evolution of TDE rates with different disk parameters (α, ε) in Figure 3.The green, orange, blue and red lines show the disk parameters α = 0.1, 0.3 and ε = 0.1, 0.4 with M BH = 10 6.75 M ⊙ , respectively.The two red dotted lines divide the evolutionary stage into the active stage, transitional stage and quiescent stage.The active, transitional and quiescent stages are divided by the two red dotted lines.In the active stage, the low-inclination stars would fall into accretion disk and be in circular orbits due to star-disk interactions within a timescale of 10 2−4 yr (e.g., Tremaine et al. 2002;Pan & Yang 2021) , while the high-inclination stars remain in the cluster.Therefore, the TDE rates in the active stage are mainly caused by the scattering of stars in high-inclination orbits, which are close to that in the quiescent galaxies as shown in Figure 3.During the transitional stage, when the inner standard thin disk disappearances, many previously embedded stars within the accretion disk are scattered, leading to a rapid increase in TDE rates.We find that the TDE rates during this stage are about 10 −2 − 10 −1 yr −1 gal −1 , which are enhanced by 2-3 orders of magnitude compared to that in the normal quiescent galaxies.
After the transitional stage, the TDE rates gradually decrease and return to levels similar to those in quiescent galaxies after about 1 Gyr.For the four different sets of disk parameters, the average TDE rates over 1 Gyr are about 320 (65, 22 and 5) times higher than those in normal quiescent galaxies for the green (orange, blue and red) lines as shown in Figure 3, which is consistent with the overrepresentation of TDEs in post-starburst galaxies (e.g., French et al. 2020a).

DISCUSSIONS
In this work, we explore the TDE rates based on a modified loss cone theory, where both star-disk interactions and star formation in the outer disk are considered.We find that the star formation in the outer unstable region of the AGN disk and the star-disk interactions will lead to the accumulation of stars in the nuclear region of galaxies.As the orbital circularization induced by the density wave disappears in the low-luminosity AGNs and quiescent galaxies, the stars will be scattered to the loss cone boundary and lead to a rapid increase in the TDE rates, which offers a possible explanation for the overrepresentation of TDE rates observed in green valley (or post-starburst) galaxies (French et al. 2016(French et al. , 2017;;Law-Smith et al. 2017;French et al. 2020b;Hammerstein et al. 2021Hammerstein et al. , 2023;;Yao et al. 2023).
Considering the in-situ formation and capture effect, it is not surprising for the presence of overdense nuclear star clusters and compact remnants in the AGN disk.The disk stars experience much different interactions than the stars in a vacuum, where the disk stars could suffer star-disk interactions, accretion, etc.The possible gravitational wave event associated with EM counterparts of GW190521 from the merger of two BHs is suggested to be the possible candidate that formed in AGN disk (Ashton et al. 2021).Furthermore, it is suggested that several AGNs may be super-solar metallicity based on the line ratios from the broad lines, which may be enriched by the supernovae due to the disk stars (e.g., Toyouchi et al. 2022;Fan & Wu 2023).The TDE rates are closely correlated with the number density of stars in the nuclear region, where the nuclear star clusters may play a key role in connecting the SMBH growth and properties of host galaxies.Many observations and theories showed a strong correlation between starburst and AGN activity (e.g., Brotherton et al. 1999;Satyapal et al. 2005;Wu & Cao 2006;Kaviraj et al. 2011;Dai et al. 2018), which could be triggered by the galactic mergers (Hopkins & Quataert 2010;Volonteri et al. 2015;Liao et al. 2023).The starburst, nuclear star cluster and SMBH activity in the host galaxy should be closely correlated, even though there is a possible delay (about a few 10 8 yr) between the peak of star formation and AGN activity (Davies et al. 2007;Schawinski et al. 2009;Wild et al. 2010;Hopkins 2012;Blank & Duschl 2016).Future highresolution optical/UV observations on the nuclear star clusters with different nuclear SMBH activities can further shed light on this issue.
However, to explain the TDE rates overenhancement in green valley (or post-starburst) galaxies, we have adopted several assumptions, which are worth further discussion: • Due to the complexity of black hole accretion and the limited observational data available, the evolution of the accretion rate and the disk model are quite uncertain.We assume that the accretion rate evolves as a double-pow-law form and divide AGN evolution into three stages: active stage ( ṁ > ṁcrit,1 ), transitional stage ( ṁcrit,2 ≤ ṁ ≤ ṁcrit,1 ) and quiescent stage ( ṁ < ṁcrit,2 ).The transitional accretion rate could greatly affect the time at which the peak of the TDE rates occur, while the disk model parameters (e.g., α and ε, see Dittmann & Miller 2020) could affect the magnitude of the peak of the TDE rates.However, the uncertainties of the disk model and the accretion rate's evolution will not affect the main conclusion that the enhancement of TDE rates will happen just after the disk transition.
• Though the correlation between starburst and AGN activity has been suggested by many authors, the quantitative relationship between accretion rate and SFR is still unknown.It is difficult to confirm whether the transitional stage of the AGN corresponds to the stage of green valley (or poststarburst) galaxies.
• The spherical symmetry of the stars' distribution will be broken by the star-forming disk, which could affect the background gravitational potential φ (r).However, the peak of star formation rate in the disk is about 0.01 M ⊙ • yr −1 for SMBH mass 10 6.75 M ⊙ , which is located at 10 5 -10 6 r g .In this region, the background gravitational potential is still dominated by the central SMBH and the original spherical symmetry star cluster where the gravitational potential asymmetry caused by the accretion disk and disk-component stars is less than ten percent.Therefore, we still assume that the background gravitational potential φ (r) is a constant function as shown in equation ( 6), which is spherically symmetric and stable.
• The Kozai-Lidov effect (Lidov 1962;Kozai 1962) could be triggered by the accretion disk (Karas & Šubr 2007), which can cause a secular oscillation of the clustercomponent stars' inclination and eccentricity.The analytical and simulated studies have shown that the Kozai-Lidov effect could increase TDE rates by several times for the disk with mass M d ∼ 0.01M BH (Karas & Šubr 2007;Terquem & Ajmia 2010;Kennedy et al. 2016), which means that the TDE rate could be higher than that in the active stage as shown in Figure 3.Note that the origin of large amplitude AGN variability is poorly understood, which could make it difficult to distinguish it from TDEs (see van Velzen et al. 2020, for a review).Up to now, only a few TDE candidates in AGNs have been found, most of which were discovered in lowluminosity galaxies (e.g., Blanchard et al. 2017;Yan & Xie 2018;Liu et al. 2020;Homan et al. 2023).In this work, we neglect the Kozai-Lidov effect since we do not aim at the TDE rate in the active stage.In the transitional stage, the number of the disk-component stars is about 2-3 orders of magnitude higher than that of the cluster-component stars within the region where the peak of star formation rate in the disk is located (∼ 10 5 − 10 6 r g ).However, for stars with a low (or zero) inclination angle with respect to the disk, the Kozai-Lidov process can not be efficiently excited (e.g., Naoz et al. 2013).Compared to the TDE rate enhancement by in-situ star formation, TDE rate enhancement by the Kozai-Lidov effect should be less important.Thus, we do not include this effect in our calculation, which needs further study in our future work.Note that we also neglect the influence of resonant relaxation since it has no significant influence on the time average TDE rate (Rauch & Tremaine 1996).
For the different TDE overenhancement mechanisms, more observational tests are preferred to distinguish them.For instance, in the cases that the TDE host galaxies exhibit signs of SMBH binary activity and/or evidence of kicked SMBH, it may indicate a significant influence of SMBH binary on TDEs.For other secular dynamical effects, such as radial anisotropies, nuclear triaxiality and eccentric nuclear stellar disk, the TDE host galaxies should be distinguishable based on their morphology.Our model is similar to the model of overdense central star clusters, which could be proved by future higher-resolution imaging.Indeed, high-resolution HST images have shown that the TDE host galaxies have higher central concentrations on pc-kpc scale (Law-Smith et al. 2017;Graur et al. 2018;French et al. 2020a).In addition, we calculated the TDE rates by considering possible star-disk interactions in both the active and non-active phases of galaxies.Our model suggests a significant enhancement in TDE rates during the turn-off of AGNs, coinciding with the disappearance of the standard thin disk as the accretion rate decreases.In this case, we connect the TDE rates with the evolution of galaxies and predict that a fraction of SMBHs will exhibit weak nuclear activities in TDE galaxies.It should be noted that this fraction may also suffer uncertainties in the timescale for transitioning from an active to a non-active state.In the currently observed TDEs, several candidates are situated within low-luminosity AGNs (e.g., Blanchard et al. 2017;Yan & Xie 2018;Liu et al. 2020;Shi et al. 2022;Homan et al. 2023), which support our predictions.Further observations of TDEs and their host galaxies in the time-domain era can shed light on the relationship between SMBH activities and the evolution of the galactic nuclear environment.Furthermore, future high-sensitivity X-ray surveys targeting TDE host galaxies can reveal the nuclear activities of SMBHs directly.
We thank the referee for the thorough and helpful suggestions that greatly improved the manuscript.M.W. thanks Xiao Fan, Jiancheng Wu and Xiangli Lei for helpful discussions.Y.M. is supported by the university start-up fund provided by Huazhong University of Science and Technology.Q.W. is supported by the National SKA Program of China (2022SKA0120101), the National Natural Science Foundation of China (grants U1931203 and 12233007) and the science research grants from the China Manned Space Project (No. CMS-CSST-2021-A06).The authors acknowledge Beijing PARATERA Tech CO., Ltd. for providing HPC resources that have contributed to the research results reported within this paper.

Figure 1 .
Figure 1.Schematic diagram of accretion disk evolution: the supermassive black hole (SMBH) resides at the galaxy's center, accompanied by a surrounding accretion disk and a population of stars.From left to right, Figures (a), (b), and (c) represent the SMBH stay in the active stage, transitional stage and quiescent stage, respectively.The red dashed lines in Figure (a) and (b) show the star in the circular orbit and eccentric orbit, respectively.The active stage (a): apart from the stars staying outside of the AGN disk, many stars are formed in the outer star-forming disk or captured by the disk.The stars embedded in the AGN disk normally stay in circular orbit as the disk gas, which are far from the loss cone boundary.The transitional stage (b): the standard thin disk (Shakura-Sunyaev disk, or SSD, see Shakura & Sunyaev 1973) gradually switches into an optically thin, geometrically thick advection-dominated accretion flow (ADAF) within a truncation radius(r tr ) after the accretion rate is lower than a critical value.The disk-component stars can be quickly scattered into the loss cone boundary through two-body relaxation and lead to a rapid increase in the TDE rates.The quiescent stage(c): the galaxy enters the quiescent stage at a very low accretion rate, where the SSD would fully disappear.

Figure 2 .
Figure 2. The distribution function of the disk-component stars in the (E, R) space, g(t, E, R), in the unit of pc −3 • σ −3 with SMBH mass M BH = 10 6.75 M ⊙ .(a), (b), (c) and (d) show the distribution function in the initial stage, active stage ( ṁ = 0.06), transitional stage ( ṁ = 0.009) and the quiescent stage according to the different accretion rate, respectively.The stars will gather at the R = 1 boundary in the (E, R) space in the active stage, while they will be scattered to the loss cone boundary in the transitional stage.

Figure 3 .
Figure 3. Evolution of TDE rates over time with M BH = 10 6.75 M ⊙ .The blue, orange, green and red show the evolution of TDE rates with the different disk parameters ( α = 0.1, 0.3 and ε = 0.1, 0.4, which is adapted from Dittmann & Miller 2020), respectively.The two red dotted lines divide the evolutionary stage into three parts: the active stage, the transitional stage and the quiescent stage.