🔭 astrophysics

TDEs on FIRE: Illuminating the Cosmic Evolution of Tidal Disruption Rates

Using FIRE-2 cosmological simulations, this study reveals that the per-galaxy tidal disruption rate peaks at redshift $z \sim 2.5$ and correlates strongly with star formation and central stellar density, while highlighting the significant role of satellite galaxies in probing intermediate-mass black holes across cosmic time.

Original authors: Rudrani Kar Chowdhury, Lixin Dai, Janet N. Y. Chang, Tsang Keung Chan

Published 2026-06-04

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Rudrani Kar Chowdhury, Lixin Dai, Janet N. Y. Chang, Tsang Keung Chan

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a vast, bustling city where massive black holes are the "mayors" sitting in the center of every galaxy. Usually, these mayors are quiet, but sometimes, a stray star wanders too close and gets ripped apart by the mayor's immense gravity. This dramatic event is called a Tidal Disruption Event (TDE). It's like a cosmic fireworks display that flashes brightly for a few weeks before fading away.

For a long time, astronomers have studied these fireworks in our "local neighborhood" of the universe (nearby galaxies). But they didn't know what was happening in the distant past, when the universe was young and the galaxies were still being built.

This paper, titled "TDEs on FIRE," uses a super-powerful computer simulation called FIRE-2 to travel back in time and count how often these fireworks happen between redshifts $z=1$ and $z=10$ (a time range covering the early universe up to about 8 billion years ago).

Here is the story of what they found, explained simply:

1. The Cosmic Fireworks Calendar

The researchers created a "calendar" of these events. They found that the frequency of these stellar fireworks wasn't constant; it had a schedule.

The Rise: In the very early universe, the rate of these events started low.
The Peak: The action got most intense around 2.5 billion years after the Big Bang (a time astronomers call redshift $z \approx 2.5$ ). This was the "golden hour" for TDEs.
The Decline: After that peak, the rate started to slow down as the universe got older and closer to today.

2. The Connection to Star Birth

Why did the fireworks peak at that specific time? The paper found a strong link between TDEs and Star Formation Rate (SFR).

The Analogy: Think of a galaxy as a factory. The "Star Formation Rate" is how fast the factory is churning out new cars (stars). The "TDE Rate" is how many of those cars crash into the mayor's office.
The Finding: The more stars the galaxy factory was producing, the more likely it was that a star would get lost and crash into the black hole. The peak in TDEs happened right when the universe was having its biggest "star-birth party."

3. The Size of the Mayor Matters

The paper looked at how the size of the black hole (the "Mayor") affects the crash rate.

The Sweet Spot: They found that the most frequent crashes happen with medium-sized black holes (called Intermediate-Mass Black Holes, or IMBHs).
Too Small or Too Big: If the black hole is too small, it doesn't have enough gravity to catch stars often. If it's too huge (a Supermassive Black Hole), it's actually less likely to rip a star apart because its gravity is so smooth and wide that it might just swallow the star whole without making a flashy explosion.
The Trend: This "sweet spot" for medium black holes exists in the early universe just as it does in our local neighborhood today. It suggests that the rules of how black holes and galaxies grow together have stayed the same for billions of years.

4. The "Crowded Room" Effect

The researchers also looked at how crowded the center of the galaxy is.

The Analogy: Imagine a dance floor. If the room is packed tight with dancers (stars) and the floor is steep (a steep density slope), it's much easier for someone to get pushed into the center.
The Finding: Galaxies with very dense centers and steep star clusters produced the most TDEs. The simulation showed that in the early universe, these "crowded dance floors" were very common, contributing to the high rate of events.

5. The Hidden Guests (Satellite Galaxies)

Not all galaxies are the big, main ones. Some are smaller "satellite" galaxies orbiting the big ones.

The Surprise: The paper found that these smaller, satellite galaxies are actually fireworks factories too. In the early universe, a huge chunk of all the TDEs happened in these small satellites, not just the big main galaxies.
Why it matters: This is exciting because these small satellites often host the medium-sized black holes that are hard to find. If we look for these "off-center" fireworks, we might finally find the missing middle-sized black holes that have been hiding in the dark.

6. Can We See Them?

The paper concludes by checking if our future telescopes (like the Nancy Grace Roman Space Telescope and the Vera Rubin Observatory) can actually spot these events.

The Verdict: Yes! The simulation suggests that most of these events in satellite galaxies are far enough away from the main galaxy to be seen clearly by these new telescopes. They won't get lost in the glare of the main galaxy.

Summary

In short, this paper used a high-definition cosmic movie (the FIRE-2 simulation) to show us that:

Tidal disruption events were most common when the universe was young and busy making stars.
The rate of these events is tightly linked to how fast a galaxy is making new stars.
Medium-sized black holes are the "stars" of the show, causing the most disruptions.
Small satellite galaxies are major contributors to these events, offering a new way to find hidden black holes.

This work gives astronomers a roadmap for what to look for when the next generation of telescopes starts scanning the sky, helping us understand how black holes and galaxies have grown together since the dawn of time.

Technical Summary: TDEs on FIRE

Problem Statement
While tidal disruption events (TDEs) have been extensively characterized in the local universe, their prevalence and cosmic evolution at high redshifts ( $z > 1$ ) remain largely unconstrained both observationally and theoretically. Current observational samples are flux-limited and overwhelmingly local, with only a handful of jetted candidates detected at $z \sim 1$ . Theoretical predictions face significant uncertainties due to the lack of resolved stellar profiles in high-redshift galaxies and poor constraints on the black hole mass function (BHMF) for low-mass black holes, which dominate the TDE population. Consequently, there is no robust prediction for the redshift evolution of TDE rates (TDR) or their dependence on host galaxy properties across cosmic time.

Methodology
This study utilizes the FIRE-2 (Feedback In Realistic Environments) cosmological zoom-in simulations, specifically the Massive Halo suite (simulations A1, A2, A4, and A8), to compute intrinsic TDE rates over the redshift range $z = 1–10$ . The methodology integrates three key components:

Simulation Framework: The FIRE-2 simulations employ the GIZMO code with high spatial resolution (reaching parsec scales) and comprehensive modeling of star formation, evolution, and stellar feedback. The Massive Halo suite includes a sophisticated model for the formation and growth of massive black holes (MBHs), seeding them at $M_{seed} = 1.4 \times 10^4 M_\odot$ when dark matter halos exceed a specific mass threshold. Black holes are treated as collisionless sink particles with artificial dynamical mass to ensure they remain at the galaxy center.
Stellar Profile Analysis: The authors extract radial stellar density profiles within 1 kpc of the galaxy centers (rather than the BH position, to account for potential dynamical friction offsets). These profiles are fitted with a double power-law function. Systems with unphysical inner slopes (steeper than outer slopes) or numerical noise are filtered out. The sample is restricted to galaxies with $M_{gal} > 10^7 M_\odot$ and $M_{BH} < 10^8 M_\odot$ (the threshold above which stars are swallowed whole rather than tidally disrupted).
Rate Calculation: TDE rates are calculated using loss cone theory and the publicly available code phaseflow. The code solves for the stellar distribution function $f(E)$ via Eddington inversion from the fitted density profiles. The total consumption rate is derived by integrating the scattering of stars into the loss cone due to two-body relaxation, assuming a monochromatic Sun-like stellar population ( $m_\star = 1 M_\odot, r_\star = 1 R_\odot$ ).

Key Contributions and Results

Redshift Evolution of TDR: The study provides the first prediction of TDR evolution across $z = 1–10$ . The averaged TDR rises from the early universe, peaks at approximately $4 \times 10^{-4} \text{ yr}^{-1}$ near $z \sim 2.5$ , and declines to $\sim 10^{-5} \text{ yr}^{-1}$ at $z = 1$ . A power-law fit for $z \ge 2.5$ is provided: $\log_{10}(\text{TDR}) \approx z^{-0.2} - 4$ .
Correlation with Star Formation: A strong correlation is found between TDR and the host galaxy's star formation rate (SFR). The TDR tracks the cosmic SFR history, peaking at $z \sim 2.5$ . The scaling relation is derived as $\text{TDR} \propto \text{SFR}^{0.69}$ .
Dependence on Masses ( $M_{BH}$ and $M_{gal}$ ):
- Black Hole Mass: The TDR increases with $M_{BH}$ in the intermediate-mass black hole (IMBH) regime, peaks at $M_{BH} \sim 1.8 \times 10^5 M_\odot$ , and then declines for supermassive black holes (SMBHs). This "turnover" trend qualitatively matches findings in the local universe, despite differences in normalization and slopes.
- Galaxy Mass: Similarly, the TDR rises with $M_{gal}$ for small galaxies, peaks at $M_{gal} \sim 5 \times 10^9 M_\odot$ , and declines for more massive systems.
Stellar Structural Drivers: The TDR is strongly driven by central stellar properties. Higher inner stellar densities ( $\rho_0$ ) and steeper inner density slopes ( $\gamma$ ) correlate with higher rates. The study highlights that the highest TDRs occur in systems with high $\rho_0$ , steep $\gamma$ , and small black hole radii of influence ( $r_{inf}$ ).
Satellite Galaxies: Satellite galaxies (non-primary hosts) exhibit TDRs comparable to primary galaxies. Their fractional contribution to the total TDR increases significantly at high redshifts, suggesting they are crucial for probing IMBHs and early galaxy assembly. Most satellite TDEs in the sample fall within the angular resolution limits of next-generation facilities like the Nancy Grace Roman Space Telescope, the Vera Rubin Observatory, and ULTRASAT.

Significance and Claims
The paper claims that cosmological simulations, specifically FIRE-2, offer a promising avenue for constraining the cosmic evolution of TDE rates. The work demonstrates that the qualitative trends of TDR with respect to black hole and galaxy mass, observed in the local universe, persist from high redshift to the present day, suggesting consistent black hole-galaxy scaling relations across cosmic time.

The authors modestly acknowledge limitations, including the absence of AGN feedback in the Massive Halo suite (which may artificially enhance SFRs and central densities at $z \lesssim 2$ ), the inability to fully resolve nuclear star clusters due to finite resolution, and the restriction of the simulation to $z=1$ (preventing constraints at $z=0$ ). Despite these caveats, the study establishes a framework for comparing theoretical predictions with upcoming observations from wide-area, high-cadence surveys, aiming to validate models of black hole growth, stellar dynamics, and galaxy co-evolution.