⚛️ general relativity

Ringdown in Vaidya spacetimes: time-dependent frequencies, Penrose limit and time-domain analyses

This paper investigates the characterization of ringdown waves in dynamical Vaidya spacetimes by extending the correspondence between quasinormal frequencies and unstable null geodesics (via the Penrose limit) from static black holes to time-dependent scenarios, validating the approach through comparisons with numerically calculated waveforms.

Original authors: Chul-Moon Yoo, Masashi Kimura, Akihiro Ishibashi, Rikuto Ohashi

Published 2026-02-23

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Chul-Moon Yoo, Masashi Kimura, Akihiro Ishibashi, Rikuto Ohashi

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 a black hole not as a static, frozen monster, but as a living, breathing entity that is constantly eating or spitting out matter. When two black holes crash into each other, they don't just stop; they "ring" like a struck bell, sending out ripples in space-time called gravitational waves. This final ringing phase is called the ringdown.

For decades, physicists have known a secret code to decode these rings. In a calm, unchanging universe (a static black hole), the pitch and speed of the ring are directly tied to the behavior of light trapped in a circular orbit right around the black hole's edge. It's like knowing the sound of a bell by studying the vibration of the metal right where the hammer hits.

This paper asks a big question: What happens when the black hole is changing? What if it's actively swallowing a cloud of gas (accreting) while it rings? Does the old code still work?

Here is the breakdown of their discovery, using some everyday analogies:

1. The "Bell" and the "Light Trap"

In a normal, static black hole, there is a specific zone called the photon sphere. Imagine this as a racetrack right above the black hole's surface where light can run in circles. However, it's a very unstable track; if a photon (a particle of light) gets nudged even slightly, it either falls into the black hole or escapes to infinity.

Physicists found that the "ringing" of the black hole (the gravitational waves) is essentially the sound of these unstable light orbits vibrating.

The Pitch (Frequency): Determined by how fast the light runs around the track.
The Fade (Damping): Determined by how quickly the light falls off the track.

2. The "Penrose Limit": A Zoom Lens

To study this, the authors use a mathematical tool called the Penrose Limit. Think of this as a super-powered microscope or a zoom lens.

Instead of looking at the whole messy universe, the Penrose Limit zooms in so tightly on that unstable light track that the rest of the universe disappears.
In this zoomed-in view, the complex, curved space-time of the black hole simplifies into a flat, easy-to-solve "plane wave."
In a static universe, this zoomed-in view perfectly predicts the sound of the bell.

3. The Problem: The Black Hole is Eating

The real universe isn't static. Black holes eat matter. When a black hole eats, its mass changes, and the "racetrack" for the light moves and stretches.

The Question: If we zoom in on a moving, changing racetrack (a dynamical photon sphere) using our Penrose microscope, can we still predict the sound of the bell?
The Analogy: Imagine trying to predict the sound of a bell while someone is stretching the metal of the bell while you are hitting it. The rules get messy.

4. The Experiment: Simulation vs. Theory

The authors did two things:

The Theory (The Zoom Lens): They used the Penrose Limit on a "Vaidya spacetime" (a mathematical model of a black hole eating or spitting out "null dust"—a fancy term for streams of light-like particles). They calculated what the ringdown should sound like if the zoom lens was perfect.
The Reality (The Full Simulation): They ran massive computer simulations of the entire black hole and the waves traveling through it, without using the zoom lens. This is the "real" sound.

5. The Findings: It's Complicated, But Close

Here is what they discovered:

The "Redshift" Effect: When the black hole is changing, the waves traveling from the light track to the observer get stretched (redshifted), like a siren changing pitch as an ambulance drives away. This effect messes up the simple prediction from the zoom lens.
The "Scattering" Effect: As the waves travel through the changing space-time, they bounce and scatter off the changing geometry, like a ball bouncing through a room where the walls are moving. The zoom lens (Penrose Limit) only sees the ball right at the start; it doesn't see the bouncing on the way out.
The Ratio Saves the Day: While the absolute pitch and volume were hard to predict with the simple zoom lens, the ratio between the pitch and the fade (how fast it dies out) was much more stable.
- Analogy: Even if the bell is being stretched and the sound is traveling through a windy room, the relationship between the tone and the decay remains surprisingly consistent with the simple model, provided you wait long enough and look from far enough away.

The Bottom Line

The paper concludes that the "Zoom Lens" (Penrose Limit) is a powerful tool, but it has limits in a dynamic universe.

It works well for understanding the local physics right next to the black hole.
It needs a correction to account for the journey the waves take to reach us (the scattering and redshift).

Why does this matter?
As we listen to the universe with gravitational wave detectors (like LIGO), we are hearing black holes that are often in the middle of eating or merging. If we want to use these "rings" to test Einstein's theory of gravity or measure black hole properties, we need to understand how the "eating" changes the sound. This paper provides a new map for translating the complex, messy sounds of a changing black hole back into the simple, clean physics of light orbits.

In short: They proved that even when a black hole is busy eating, the "ringing" still carries the fingerprint of the light trapped around it, but you have to account for the "noise" of the journey to hear it clearly.

1. Problem Statement

The paper addresses the characterization of gravitational wave ringdown signals (Quasinormal Modes, or QNMs) in dynamical spacetimes, specifically the Vaidya spacetime, which models a spherically symmetric black hole accreting or radiating null dust.

Context: In static, stationary spacetimes (like Schwarzschild or Kerr), a deep correspondence exists between the high-frequency (eikonal) QNMs and the properties of unstable circular photon orbits (UCOP). Specifically, the real part of the QNM frequency ( $\omega_{Re}$ ) corresponds to the orbital angular velocity ( $\Omega$ ), and the imaginary part ( $\omega_{Im}$ , the damping rate) corresponds to the Lyapunov exponent ( $\lambda$ ) of the orbit. This is often derived using the Penrose Limit (PL), which approximates the spacetime near a null geodesic as a plane wave.
The Gap: It is unclear if this QNM-geodesic correspondence and the utility of the Penrose Limit hold in dynamical backgrounds where the mass function $M(V)$ changes with time. Astrophysical black holes are rarely perfectly isolated or stationary; they interact with matter (accretion) and radiation.
Objective: To extend the PL analysis to the dynamical photon sphere in Vaidya spacetimes and determine the extent to which the PL geometry can predict the time-dependent ringdown frequencies and damping rates observed in full numerical simulations.

2. Methodology

The authors employ a multi-pronged approach combining analytical approximations, frequency-domain calculations, and time-domain numerical simulations.

A. Analytical Framework: The Penrose Limit (PL)

Static Case Review: They first review the derivation of the PL metric around a UCOP in a static spacetime. The metric takes the form of a plane wave (Brinkmann coordinates): $ds^2 = 2dudv + (\alpha^2 x^2 - \beta^2 y^2)du^2 + dx^2 + dy^2$ $d s^{2} = 2 d u d v + (α^{2} x^{2} - β^{2} y^{2}) d u^{2} + d x^{2} + d y^{2}$ .
- $\beta$ relates to the angular velocity $\Omega$ .
- $\alpha$ relates to the Lyapunov exponent $\lambda$ .
- Solving the wave equation in this metric yields the QNM spectrum: $\omega \approx \Omega(\ell + 1/2) - i\lambda(n + 1/2)$ .
Dynamical Extension (Vaidya):
- They consider the Vaidya metric with a time-dependent mass function $M(V)$ .
- They identify the trajectory of the dynamical photon sphere $r_{PS}(V)$ , which evolves as the mass changes.
- Adiabatic Approximation: They calculate the PL geometry along this evolving trajectory. Assuming the change is slow enough, they define time-dependent effective parameters $\Omega_{ad}(V)$ and $\lambda_{ad}(V)$ derived from the parallelly transported zweibeins and the Riemann tensor components along the null geodesic.
- They predict a time-dependent complex frequency $\omega_{PL, ad}(V)$ .

B. Frequency-Domain Analysis (Constant Accretion)

For a linear mass increase $M(V) = M_1 + M'V$ , the spacetime is conformally static.
The authors transform the master equations for scalar, electromagnetic, and gravitational perturbations into a Schrödinger-like form with an effective potential.
They solve for QNMs using Leaver's method (continued fractions) to obtain exact frequencies in the conformal static frame.

C. Time-Domain Numerical Simulations

They perform full numerical simulations of tensor perturbations in the full Vaidya spacetime (both constant and time-dependent accretion rates) using the double-null formalism.
Extraction Method: Instead of a global Fourier transform (which fails for time-dependent frequencies), they extract local frequencies by analyzing the waveform peaks.
- For peaks at $V_n$ and $V_{n-1}$ , they calculate the instantaneous real frequency $\omega_{Re, n}$ and damping rate $\omega_{Im, n}$ based on the phase difference and amplitude decay.
- They specifically analyze the ratio $R_n = |\omega_{Im, n} / \omega_{Re, n}|$ to minimize redshift effects.

3. Key Contributions and Results

A. Constant Accretion Rate ( $M' = \text{const}$ )

Conformal Staticity: The Vaidya spacetime with constant accretion is conformally equivalent to a static spacetime.
Validation: The authors found that the QNMs calculated via the PL analysis, the frequency-domain Leaver method, and the time-domain numerical simulations agree well.
Correction Factor: A small deviation exists between the PL prediction and the exact numerical results due to the specific field spin (scalar vs. gravitational). This is quantified by a correction factor $C \approx 1$ (derived from the Schwarzschild limit), which allows the PL to accurately predict the ringdown behavior even in the presence of constant accretion.

B. Time-Dependent Accretion Rate (General Case)

Setup: They simulated a scenario where the accretion rate varies (e.g., a smooth transition from $M_1$ to $M_2$ ).
Frequency Evolution: The instantaneous frequency $\omega_{Re}(t)$ and damping rate $\omega_{Im}(t)$ evolve monotonically from the values associated with mass $M_1$ to those of $M_2$ .
The Ratio $R$ : The ratio $R = \omega_{Im}/\omega_{Re}$ $R = ω_{I m} / ω_{R e}$ is the most robust observable.
- PL Prediction: The adiabatic PL analysis predicts a specific time-dependent trajectory for $R$ .
- Numerical Results: The numerical simulation shows a qualitatively similar trend to the PL prediction (a temporary reduction in $R$ during the accretion phase).
- Quantitative Discrepancy: There is a significant quantitative deviation (approx. 3 times the difference in the $M_2=3M_1$ case) between the PL prediction and the numerical $R_n$ .
Cause of Deviation: The authors attribute this discrepancy to scattering effects during the wave propagation from the dynamical photon sphere to the distant observer.
- The PL captures only the local geometry at the photon sphere.
- In a dynamical spacetime, the wave propagates through a time-varying background, accumulating phase shifts and scattering that the local PL approximation cannot capture.
- The redshift effect is dominant, shifting both real and imaginary parts equally, which is why the ratio $R$ is the preferred metric for comparison.

C. Parameter Dependence

Accretion Rate: As the accretion period becomes longer (adiabatic limit) for a fixed maximum rate, the numerical results converge closer to the PL prediction.
Angular Momentum ( $\ell$ ): The agreement does not significantly improve for higher $\ell$ (tested up to $\ell=10$ ). The $\ell=4$ case is already sufficient to be in the eikonal regime; the remaining error is due to the global scattering, not the local approximation.
Observer Distance: The results converge as the observer radius $r_{obs}$ increases, confirming that the deviation is a propagation effect.

4. Significance and Conclusions

Validity of PL in Dynamics: The study confirms that the Penrose Limit provides a valid local description of the ringdown origin in dynamical spacetimes. The time-dependent parameters $\Omega_{ad}$ and $\lambda_{ad}$ correctly track the evolution of the photon sphere.
Limitations: The PL is insufficient for predicting the exact waveform observed at infinity in highly dynamical scenarios because it neglects global scattering effects and the history of the wave propagation through the time-dependent background.
Physical Insight: The ratio of the damping rate to the frequency ( $R$ ) is shown to be a robust observable that encodes the geometry of the dynamical photon sphere, largely independent of the redshift effects that plague the individual frequency components.
Astrophysical Implication: Even in non-stationary environments (like accreting black holes), the ringdown signal carries information about the instantaneous geometry of the photon sphere. However, extracting this information requires accounting for the scattering history, which the simple PL cannot provide alone.

In summary, the paper successfully bridges the gap between static QNM-geodesic correspondence and dynamical spacetimes, demonstrating that while the local physics (PL) dictates the emission properties, the global propagation dynamics introduce significant corrections that must be modeled for precise gravitational wave spectroscopy in evolving systems.