Spin the black circle II: tidal heating and torquing of a rotating black hole by a test mass on generic orbits

This paper presents a numerical study of energy and angular momentum fluxes for a test particle on generic orbits around a Kerr black hole, revealing complex flux behaviors in eccentric and hyperbolic trajectories and proposing a new factorized, resummed analytical model that significantly improves the prediction of superradiant regimes and flux accuracy in the strong-field limit compared to existing post-Newtonian expressions.

The Gist

Imagine a black hole not as a static, silent vacuum cleaner, but as a spinning, cosmic whirlpool. Now, imagine throwing a small pebble (a test mass) into this whirlpool. Depending on how you throw it, the pebble might circle gently, swing wildly in an oval, or zoom past in a straight line before being flung away.

This paper, titled "Spin the black circle II," is a detailed study of what happens to the black hole's "skin" (its event horizon) when these different kinds of pebbles interact with it. The authors are essentially asking: How much energy and spin does the black hole steal from the pebble, or how much does it give back?

Here is the breakdown of their discovery using everyday analogies:

1. The Setting: The Spinning Whirlpool

Think of the black hole as a giant, spinning merry-go-round.

• The Horizon: This is the edge of the merry-go-round. Once you cross it, you can't get off.
• Tidal Heating and Torquing: As the pebble orbits, it creates ripples in the fabric of space-time (like a boat creating waves in a lake). These ripples hit the black hole's edge. Sometimes, the black hole absorbs the energy of these ripples (heating up, gaining mass). Sometimes, if the black hole is spinning fast enough, it actually pushes back against the ripples, stealing energy from the black hole's own spin and giving it to the pebble. This is called Superradiance.

2. The Experiment: Three Types of Throws

The researchers simulated three different ways the pebble could move around the black hole:

• Circular Orbits: The pebble spins in a perfect circle, like a satellite.
• Eccentric Orbits: The pebble swings in a stretched-out oval, getting very close to the black hole and then swinging far away, like a comet.
• Hyperbolic Orbits: The pebble zooms in, swings around the black hole in a quick U-turn, and shoots back out into space, never to return.

The Surprise:
For the circular orbits, scientists already knew the rules. But for the eccentric and hyperbolic orbits, things got messy and exciting.

• The "Flashlight" Effect: In a perfect circle, the energy flow is steady. But in an oval or a fly-by, the energy flow is like a strobe light. It flashes on and off, sometimes even reversing direction!
• The "Tug-of-War": The paper found that the black hole doesn't just steadily gain or lose energy. It engages in a complex dance. At one moment, it's stealing spin from the pebble; a split second later, it might be giving spin back. The "instantaneous" flow of energy can change signs multiple times during a single orbit.

3. The Problem: The Math is Too Hard

Scientists have mathematical formulas (equations) to predict how much energy is exchanged.

• The Old Maps: The existing formulas work great for the "perfect circle" scenario. But when the pebble swings in a wild oval or zooms past, the old maps get lost. They are like trying to use a flat map of a city to navigate a mountain; they work okay on the flat parts but fail miserably in the steep, complex areas (the "strong-field" regime near the black hole).
• The Error: When the pebble gets very close to the black hole, the old formulas were off by as much as 100%. That's like predicting a car will go 60 mph, but it actually goes 120 mph or stops dead.

4. The Solution: A New, Smarter Formula

The authors created a new, "resummed" formula. Think of this as upgrading from a flat paper map to a 3D GPS with real-time traffic updates.

• The "Factorized" Approach: They broke the problem down. Instead of one giant, confusing equation, they separated the "circular" part (which they know well) from the "non-circular" part (the wobbles and swings).
• The "Superradiance" Switch: They identified a specific "switch" in the math that determines when the black hole starts stealing energy versus giving it back. They found that for wild orbits, this switch depends on the exact speed and direction of the pebble, not just its average speed.
• The Result: Their new formula is much better. It predicts the "switch" point (when the energy flow flips) with about 90% accuracy for most scenarios. It works for both the energy flow and the spin flow.

5. Why Does This Matter?

You might ask, "Who cares about a pebble orbiting a black hole?"

• Listening to the Universe: We have detectors (like LISA in the future) that will "hear" the gravitational waves from black holes. To understand what we hear, we need perfect models of how black holes behave.
• The "Fingerprint": If a black hole is spinning and eating a star, the way it absorbs or spits out energy leaves a specific "fingerprint" on the gravitational waves. If our math is wrong, we might misinterpret the signal and think we're seeing a different kind of object.
• Testing Einstein: This research helps us test Einstein's theory of General Relativity in the most extreme conditions possible. If the black hole behaves exactly as the new formula predicts, Einstein wins again. If not, we might find new physics.

Summary Analogy

Imagine you are trying to predict how much water a spinning bucket will splash out when you throw a ball at it.

• Old Method: You only practiced throwing the ball in a perfect circle. Your prediction works great for circles but fails when you throw the ball in a zig-zag or a straight line.
• This Paper: The authors went out and actually threw the ball in zig-zags and straight lines. They watched the water splash in real-time, saw it splash forward, then backward, then forward again. They then wrote a new rulebook that accounts for these wild splashes, allowing us to predict the water flow accurately even when the ball is thrown wildly.

This work is a crucial step in building a "user manual" for the most extreme objects in our universe, ensuring that when we finally listen to the symphony of colliding black holes, we understand every note.

Technical Summary

1. Problem Statement

The paper addresses the dynamics of energy and angular momentum fluxes absorbed by the event horizon of a rotating (Kerr) black hole due to the tidal interaction with an orbiting test mass. While the phenomenon of tidal heating and torquing is well-understood for circular orbits, its behavior on generic orbits (eccentric and hyperbolic) remains less explored, particularly in the strong-field regime.

Key challenges identified include:

• Superradiance Complexity: For circular orbits, the onset of superradiance (energy extraction from the black hole) is determined simply by the orbital frequency relative to the horizon frequency ($\Omega_H$). For generic orbits, the interplay between superradiant and non-superradiant regimes is more complex, potentially leading to multiple sign changes in instantaneous fluxes.
• Analytical Limitations: Existing Post-Newtonian (PN) analytical expressions for horizon fluxes often fail to accurately predict flux magnitudes and sign changes in the strong-field regime (close encounters) or for high eccentricities.
• Lack of Numerical Benchmarks: There was a need for a comprehensive numerical dataset covering a wide range of orbital parameters and black hole spins to validate and improve analytical models.

2. Methodology

The authors employed a hybrid approach combining high-precision numerical relativity simulations with the development and testing of improved analytical models.

A. Numerical Framework

• Setup: The team simulated a test particle of mass $\mu$ orbiting a Kerr black hole of mass $M$ and dimensionless spin $\hat{a}$. They neglected radiation-reaction effects on the orbit (geodesic approximation), which is valid for the timescales of interest in the test-mass limit.
• Solver: They utilized Teukode, a time-domain solver for the 2+1 Teukolsky equation. This solver uses a horizon-penetrating, hyperboloidal foliation to compute gravitational perturbations ($\psi_0$) and extract waveforms without extrapolation.
• Flux Calculation: Horizon fluxes of energy ($\dot{M}$) and angular momentum ($\dot{S}$) were computed using the Newman-Penrose scalar $\psi_0$ and the formalism developed by Poisson. The calculation involved summing over azimuthal modes ($m=0, 1, 2$).
• Parameter Space: The study covered 257 configurations, including:
  • Circular orbits: 8 simulations.
  • Eccentric orbits: 140 simulations (varying eccentricity $e_0 \in [0.1, 0.96]$ and semilatus rectum).
  • Hyperbolic orbits: 109 simulations (varying initial energy $E_0$ and angular momentum).
  • Spin: Black hole spins $\hat{a} \in [-0.8, 1.0]$.

B. Analytical Framework & Resummation Strategies

The authors tested and refined analytical expressions derived in their previous work (Paper I) and standard PN expansions. Key strategies included:

• Factorization: Isolating a "superradiance prefactor" that determines the sign of the flux based on the difference between a "tidal frequency" ($\Omega_T$) and the horizon frequency ($\Omega_H$).
• Multipolar Decomposition: Analyzing fluxes mode-by-mode ($m=0, 1, 2$) to understand how different modes contribute to sign changes.
• Resummation:
  • CNC Factorization: Splitting fluxes into a circular term and a non-circular correction factor.
  • K-Reparameterization: Introducing a variable $K$ to ensure non-circular corrections reduce exactly to unity for circular orbits, improving convergence.
  • DIN Resummation: Using Damour-Iyer-Nagar strategies to resum the circular part of the fluxes using 11PN test-mass results (Fujita).
  • Exact Frequency Replacement: Replacing PN approximations of orbital frequency ($p_\phi/r^2$) with exact time-domain derivatives ($\dot{\phi}$) in the prefactors.

3. Key Contributions

1. Extensive Numerical Catalog: The paper presents the most extensive numerical survey of horizon fluxes for generic orbits in Kerr spacetime to date, providing a reference dataset for future modeling.
2. Discovery of Complex Phenomenology: The authors uncovered that for eccentric and hyperbolic orbits, instantaneous fluxes can exhibit multiple peaks and sign changes. This indicates a complex interplay where superradiant and non-superradiant regimes coexist or alternate during a single encounter.
3. Mode-Dependent Superradiance: They demonstrated that different azimuthal modes ($m$) can change signs at different times, and the condition for superradiance is not uniform across all modes for generic orbits.
4. Improved Analytical Model (Model 2): They proposed a factorized and resummed representation that significantly outperforms standard PN expressions. This model:
   • Uses exact orbital frequencies in the superradiance prefactor.
   • Includes non-spinning corrections to the tidal frequency.
   • Utilizes K-reparameterized non-circular corrections.
   • Incorporates high-order (11PN) circular flux data.

4. Key Results

• Sign Changes: The numerical data confirms that for positive spins ($\hat{a} > 0$), fluxes can switch from negative (absorption) to positive (extraction) and vice versa during an orbit. The critical parameters for these switches differ for energy ($\dot{M}$) and angular momentum ($\dot{S}$), and for different $m$-modes.
• Predictivity of Superradiance Onset: The proposed analytical model predicts the onset frequency of the superradiant regime (where flux sign changes) to within 10% for $\gtrsim 73\%$ of configurations for both energy and angular momentum fluxes.
• Accuracy Metrics:
  • Circular Orbits: The model reproduces known results with high accuracy, correctly capturing the single sign change near the Last Stable Orbit (LSO).
  • Eccentric Orbits: Median relative errors are $\sim 13\%$ for $\dot{M}$ and $\sim 7\%$ for $\dot{S}$. About 43–57% of configurations are accurate within 10%.
  • Hyperbolic Orbits: Median relative errors are $\sim 22\%$ for $\dot{M}$ and $\sim 21\%$ for $\dot{S}$. About 24–28% of configurations are accurate within 10%.
• Error Correlation: The peak orbital frequency ($\Omega_{max}$) is identified as the strongest predictor of model error. Errors increase significantly when $\Omega_{max}$ pushes the PN expansion beyond its validity domain (strong-field regime).
• Quasi-Normal Modes (QNMs): The authors observed that for highly eccentric orbits and close encounters, QNM excitations dominate the fluxes during the outgoing leg, causing deviations that current analytical models cannot capture.

5. Significance

• Gravitational Wave Modeling: The results provide crucial improvements for the Effective-One-Body (EOB) formalism, which is used to generate waveform templates for gravitational wave detectors like LISA (Laser Interferometer Space Antenna). Accurate modeling of horizon fluxes is essential for extreme mass-ratio inspirals (EMRIs) and high-eccentricity binaries.
• Testing General Relativity: By refining the understanding of tidal heating and torquing, this work aids in testing the "no-hair" theorem and the nature of black hole horizons using future GW observations.
• Analytical Framework: The proposed "factorized and resummed" approach offers a robust template for extending analytical descriptions to other strong-field phenomena where standard PN expansions fail.
• Future Directions: The paper highlights the need for further work on comparable-mass binaries, spin-induced precession, and the connection between horizon dynamics (shear) and the waveform at infinity, particularly for plunging orbits where event horizon teleology poses challenges.

In summary, this work bridges the gap between numerical simulations and analytical approximations for black hole horizon fluxes, providing a validated, high-accuracy model for generic orbits that is essential for the next generation of gravitational wave astronomy.