A Semi-Analytical Loss Cone Theory for Tidal Disruption Event Rates Around Kerr Black Holes

This paper presents the first semi-analytical framework for tidal disruption event rates around spinning Kerr black holes, demonstrating that while spin induces a first-order bias favoring the disruption of retrograde stars, the global event rate remains remarkably insensitive to black hole spin.

The Gist

Imagine a galaxy as a giant, bustling city, and at its very center sits a supermassive black hole—a cosmic vacuum cleaner so powerful it can tear apart anything that gets too close. Sometimes, a wandering star gets unlucky, gets pulled in, and is ripped to shreds by the black hole's gravity. This spectacular explosion of light is called a Tidal Disruption Event (TDE).

For a long time, scientists have tried to predict how often these events happen. Think of it like trying to guess how many people will fall into a specific hole in the ground in a crowded park. The "hole" is the danger zone around the black hole.

This paper, written by Wenkang Xin, tackles a tricky problem: What happens when the black hole is spinning?

The Spinning Black Hole Problem

Most previous theories treated black holes like stationary, non-spinning spheres. But in reality, these giants spin incredibly fast. Imagine a spinning top. When a top spins, it doesn't just pull things in straight down; it drags space itself around with it.

Because of this spin, the "danger zone" (where a star gets ripped apart) isn't a perfect circle anymore. It becomes lopsided.

• The Analogy: Imagine a spinning merry-go-round with a trapdoor in the middle. If you run toward the trapdoor with the spin (prograde), the door might feel smaller or harder to reach. If you run against the spin (retrograde), the door might feel wider and easier to fall into.
• The Paper's Finding: The author calculated that stars orbiting against the spin of the black hole are actually more likely to get torn apart than stars orbiting with it. It's like the spinning black hole is "scooping up" the counter-spinning stars more efficiently.

The "Loss Cone" Game

To understand how stars get into this danger zone, the author uses a concept called the "Loss Cone."

• The Analogy: Imagine a giant, invisible funnel pointing at the black hole. Stars are like marbles rolling on a table. Most marbles roll safely around the funnel. But if a marble gets hit by another marble (a gravitational nudge from a neighbor), it might get knocked into the funnel. Once it's in the funnel, it's doomed.
• The 1D vs. 2D Problem: Old theories treated the funnel as a simple circle (1D). You only had to worry about how fast the marble was moving. But because the black hole is spinning, the funnel is now a 3D shape that depends on the angle of the marble's approach (2D). The author built a new mathematical framework to handle this complex, angled funnel.

The Big Surprise: The Total Count Doesn't Change Much

Here is the most counter-intuitive part of the paper.

Even though the spinning black hole makes it much easier for some stars (the counter-spinning ones) to get destroyed, and harder for others (the co-spinning ones), the total number of stars getting destroyed stays almost exactly the same.

• The Analogy: Imagine a busy highway with a toll booth. If the booth is slightly tilted, cars might get through the left side faster and the right side slower. However, if the total number of cars on the highway is huge and they are spread out evenly, the total number of cars passing through the booth per hour doesn't change much. The tilt just shifts which cars get through, not how many.
• The Paper's Conclusion: The "global rate" (the total number of TDEs) is remarkably robust. It doesn't care much if the black hole is spinning or not. The spin just changes the mix of stars that get eaten (more counter-spinning ones), but the total count remains stable.

Why Does This Matter?

While the total number of events doesn't change, the type of events does.

1. Star Preferences: The black hole is essentially "eating" the counter-spinning stars first. Over time, the population of stars near the black hole will change, leaving behind mostly stars that spin in the same direction as the black hole.
2. The Aftermath: When a star is torn apart, it creates a disk of debris. If the star was counter-spinning, that debris disk will be tilted at a weird angle relative to the black hole's spin. This affects how the resulting flash of light looks to us.

Summary

Wenkang Xin's paper is like building a new, more accurate map for a spinning city.

• Old Map: Assumed the city was round and still.
• New Map: Accounts for the city spinning, showing that the "danger zones" are tilted.
• Key Insight: Even though the danger zones are tilted, the total number of accidents (TDEs) remains roughly the same. The spin just changes who gets in trouble (the counter-spinning stars) and how the wreckage looks, but it doesn't change the overall traffic flow.

This work provides a simpler, semi-analytical way to calculate these rates without needing massive, time-consuming computer simulations, helping astronomers better understand the population of black holes and stars in our universe.

Technical Summary

Technical Summary: A Semi-Analytical Loss Cone Theory for Tidal Disruption Event Rates Around Kerr Black Holes

Problem Statement
Tidal Disruption Events (TDEs) occur when stars are scattered onto near-radial orbits and torn apart by a supermassive black hole (BH). While the supply of stars is governed by Newtonian stellar dynamics, the disruption threshold is set by general relativistic (GR) tidal dynamics. Existing theoretical frameworks for TDE rates face two primary limitations:

1. Non-spinning Assumption: Classical loss cone theory assumes a non-spinning (Schwarzschild) BH, ignoring the fact that astrophysical BHs are likely rapidly rotating (Kerr).
2. Dimensionality: Classical treatments model diffusion in the magnitude of angular momentum ($L$) only (1D), assuming spherical symmetry. However, a spinning BH breaks spherical symmetry, making the disruption and capture thresholds dependent on the orbital inclination (the orientation of the angular momentum vector).

Consequently, there is a lack of a self-consistent framework that incorporates spin-dependent, inclination-dependent boundaries into the stellar dynamical calculation of TDE rates.

Methodology
The dissertation constructs a semi-analytical framework that bridges GR tidal dynamics and 2D stellar diffusion.

• Relativistic Boundaries (Chapter 2): The author employs a novel tidal tensor formalism based on successive frame transformations (ZAMO to Local Freely Falling frame) rather than parallel-transported tetrads. This allows for the direct computation of the tidal tensor at the pericentre of a geodesic in Kerr spacetime.

  • The study defines two critical boundaries as functions of angular momentum magnitude ($L$) and inclination ($x = L_z/L$): the tidal disruption boundary ($L_{tde}$) and the direct capture boundary ($L_{cap}$, defined by Innermost Bound Spherical Orbits).
  • It is demonstrated that $L_{tde}(x)$ is approximately linear in $x$ with a small negative slope, while $L_{cap}(x)$ exhibits strong non-linearity for high spins.

• 1D Loss Cone Theory with Direct Capture (Chapter 3): Before extending to 2D, the author rigorously revisits the classical 1D loss cone problem.

  • A boundary layer analysis is performed to solve the Fokker-Planck equation across all regimes (empty, intermediate, and full).
  • A closed-form analytical expression for the direct capture fraction ($C_{cap}$) is derived for nested loss cones (where $L < L_{cap} < L_{tde}$). This fraction depends on the loss cone fullness parameter $q$ and the ratio $\phi_c = L_{cap}/L_{tde}$.

• 2D Loss Cone Theory (Chapter 4): The core contribution is the formulation of a 2D Fokker-Planck equation describing simultaneous diffusion in angular momentum magnitude ($R = L^2/L_c^2$) and orientation ($x$).

  • Diffusion Coefficients: Using vector perturbation and coordinate transformation methods, the author derives diffusion coefficients showing that for highly eccentric orbits in an isotropic background, dynamical friction vanishes. The collision operator reduces to a pure divergence form (flux-conserving) with no cross-terms between $R$ and $x$.
  • Perturbative Solution: Assuming the loss cone boundary is linear in $x$ with a small slope $\epsilon$, the author performs a perturbative expansion.
  • Limits: Analytical solutions are derived for the Empty Loss Cone (diffusion-limited) and Full Loss Cone (phase-space density limited) regimes.
  • Interpolation: Heuristic interpolation formulae are proposed for the intermediate regime based on the structural similarity between the limits and the 1D reciprocal-sum approximation.

Key Results

1. Preferential Disruption: The analysis reveals a strong first-order ($O(\epsilon)$) local bias favoring the disruption of retrograde stars ($x < 0$) over prograde stars ($x > 0$). This is because the loss cone boundary is larger for retrograde orbits, allowing a higher flux of stars to enter the disruption zone at those inclinations.
2. Global Rate Robustness: Despite the strong local inclination bias, the global TDE rate (integrated over all inclinations) is remarkably insensitive to BH spin.
   • In both the empty and full limits, the global flux retains the same functional form as the classical 1D isotropic case.
   • Spin-induced anisotropy only introduces a second-order ($O(\epsilon^2)$) geometric correction to the flux normalization, which is negligible for typical galactic parameters.
   • This robustness arises from the conservation of flux and the symmetry of the isotropic background distribution, which enforces that leading-order corrections must be even functions of the spin direction.
3. Direct Capture: The derived capture fraction $C_{cap}$ shows that direct capture is negligible in the empty regime but constitutes a fixed fraction of events in the full regime. Accurate modeling of the intermediate regime is essential for massive BHs where $L_{cap} \approx L_{tde}$.

Significance and Claims
The paper claims to provide the first semi-analytical framework capable of incorporating spin-dependent loss cone boundaries into TDE rate theory. Its primary significance lies in:

• Decoupling Local vs. Global Effects: It demonstrates that while BH spin fundamentally alters the inclination distribution of disrupted stars (creating a retrograde bias), it does not significantly alter the total event rate in isotropic galaxies. This challenges the notion that spin is a primary driver of global rate variations.
• Tractability: By deriving closed-form expressions for capture fractions and perturbative solutions for 2D diffusion, the work offers a computationally tractable route to include spin and inclination in population-level TDE studies, avoiding the need for purely numerical simulations for rate estimation.
• Astrophysical Implications: The retrograde bias implies that the surviving stellar population will develop a prograde bias over time, and the preferential consumption of retrograde stars acts as a secular spin-down mechanism for the BH. Furthermore, the inclination bias affects the misalignment of debris discs, influencing the observable flare dynamics (e.g., Lense-Thirring precession).

The author concludes that while the total number of TDEs is robust against spin, future population studies must account for the inherent retrograde bias to accurately decode observed TDE demographics and constrain BH spin distributions.