Bilinear products and the orthogonality of quasinormal modes on hyperboloidal foliations

This paper investigates the divergence issues in bilinear products for Schwarzschild black-hole quasinormal modes on hyperboloidal foliations, proposing regularization procedures and flux-based alternatives to define finite orthogonality relations and explicitly compute excitation factors.

The Gist

Imagine a black hole not as a cosmic vacuum cleaner, but as a giant, invisible bell. When two black holes crash into each other, they don't just disappear; they ring like a bell after being struck. This "ringing" is called a Quasinormal Mode (QNM). Just like a bell has a specific pitch and a specific way it fades away, a black hole has a unique "fingerprint" of sound waves that tell us about its mass, spin, and the laws of gravity itself.

Scientists want to listen to these sounds to test Einstein's theory of General Relativity. But there's a huge problem: The math is broken.

The Problem: The "Infinite" Bell

In standard physics, to analyze a sound, you need to measure it from start to finish. But black holes are tricky.

1. The Horizon: The edge of the black hole where nothing escapes.
2. Infinity: The far reaches of space.

When scientists try to do the math to separate these "ringing" sounds from the background noise, the equations blow up. The numbers get infinitely large at the horizon and at infinity. It's like trying to measure the volume of a bell, but your microphone explodes the moment you get too close to the bell or too far away. Because of this, the standard mathematical tools (called "inner products") that usually let us compare different sounds fail completely.

The Solution: A New Map (Hyperboloidal Foliations)

The authors of this paper say, "The bell isn't broken; our map is."

They introduce a new way of looking at spacetime called Hyperboloidal Foliations.

• The Old Way: Imagine trying to map a curved Earth using a flat piece of paper. The edges get stretched and distorted. That's what old coordinates do to black holes.
• The New Way: Imagine wrapping the Earth in a flexible, stretchy skin that hugs the curves perfectly. This new "skin" (the hyperboloidal framework) stretches the map so that the black hole's edge and the far reaches of space are brought to finite, manageable distances. Suddenly, the "infinite" numbers become finite and calculable.

The Twist: The Mirror Image (The J Operator)

Even with this new map, there's still a snag. To compare two different "ringing" sounds (modes) and prove they are distinct, you need a special mathematical operation. The authors use a tool called the J operator.

Think of the J operator as a Time-Reversing Mirror.

• If you have a sound wave traveling out from the black hole (a normal ring), the mirror flips it to show a wave traveling in from the future.
• The problem? When you look at this "mirror image" wave in the new map, it still behaves strangely at the edges. It's like looking at a reflection in a funhouse mirror that makes the edges of the room look like they are stretching into infinity again.

The Fix: Regularization (The "Mathematical Band-Aid")

The paper shows that while the math looks broken (divergent) at the edges, it's actually just a trick of perspective. The authors propose two clever ways to fix it, which they call Regularization:

1. The "Semi-Analytic" Approach: Instead of trying to calculate the whole thing at once, they break the problem into tiny, known pieces (like solving a puzzle using pre-made shapes) and then smoothly connect them. It's like calculating the area of a weird shape by filling it with perfect circles and triangles, then adjusting the edges.
2. The "Complex Contour" Approach: This is the most creative. Imagine the math is a landscape with a deep, impassable canyon (the infinity problem). Instead of trying to cross the canyon, they walk around it by stepping into a parallel dimension (the complex plane). They walk along a path where the numbers stay calm and finite, do their calculation, and then step back.

The Result: Listening to the Bell

By using these new maps and these "mathematical band-aids," the authors successfully:

• Proved Orthogonality: They showed that different black hole "ringing" modes are truly distinct and don't interfere with each other, just like different notes on a piano.
• Calculated Excitation: They figured out how to predict exactly how loud a specific ring will be based on the initial crash (the initial data).

Why This Matters

This isn't just abstract math. As our telescopes (like LIGO and the future Einstein Telescope) get better, we will hear the "ringing" of black holes with incredible clarity. This paper provides the instruction manual for decoding those sounds. It tells us how to translate the raw data from the universe into a clear understanding of what happened when those black holes collided, helping us test if Einstein was right about how gravity works in the most extreme environments in the universe.

In short: The authors fixed a broken mathematical tool by changing the map, using a mirror, and taking a detour through a parallel dimension, allowing us to finally hear the true song of the black hole.

Technical Summary

1. Problem Statement

Quasinormal modes (QNMs) are the characteristic damped oscillations of perturbed black holes, serving as "fingerprints" for black hole spectroscopy. A fundamental theoretical challenge in QNM analysis is the definition of a bilinear product (or inner product) under which QNMs with different frequencies are orthogonal.

• The Divergence Issue: In standard coordinate systems (e.g., Schwarzschild-Droste), QNM eigenfunctions diverge at the event horizon and spatial infinity. This prevents the definition of a standard $L^2$ inner product.
• The $J$-Operator Necessity: To define orthogonality, one must use a bilinear form involving a discrete symmetry operator $J$ (time and azimuthal reflection, $t \to -t, \phi \to -\phi$). This maps a retarded QNM solution to an "anti-QNM" (advanced solution).
• The Core Conflict: While the hyperboloidal framework (using coordinates that intersect future null infinity $\mathcal{I}^+$ and the horizon $\mathcal{H}^+$) renders QNM eigenfunctions finite and smooth, the anti-QNM solutions (required for the $J$-product) become singular at these same boundaries when expressed in future hyperboloidal coordinates. Consequently, the integrand of the bilinear product remains divergent, necessitating regularization.

2. Methodology

The authors extend the formalism of conserved currents (symplectic and energy) to the hyperboloidal framework to construct QNM orthogonality products.

A. Hyperboloidal Framework & The $J$ Operator

• Foliations: The paper introduces two distinct hyperboloidal foliations:
  1. Future Foliations ($\bar{x}^\mu$): Adapted to retarded (outgoing) boundary conditions. QNMs are regular here; anti-QNMs diverge.
  2. Past Foliations ($\check{x}^\mu$): Adapted to advanced (ingoing) boundary conditions. Anti-QNMs are regular here; QNMs diverge.
• The $J$ Operator: Defined as a coordinate pullback between future and past foliations. It maps a regular QNM in future coordinates to a singular anti-QNM in the same coordinates (and vice versa). This geometric interpretation clarifies that the divergence is a structural feature of the $J$-product, not a coordinate artifact.

B. Conserved Currents and Bilinear Forms

The authors construct bilinear products based on two conserved currents:

1. Symplectic Current: Derived from the wave operator structure.
2. Energy Current: Derived from the stress-energy tensor and the timelike Killing vector.
   They define an Extended Orthogonality Product that includes bulk integrals and boundary flux terms. They demonstrate that while the extended product makes the Hamiltonian formally self-adjoint, the flux terms vanish only under specific regularization schemes.

C. Regularization Strategies

To handle the divergences in the integrand (caused by the $e^{-2i\omega H(\sigma)}$ term from the $J$-mapping), two analytic continuation methods are employed:

1. Semi-Analytic Integration:
   • Expands the radial QNM functions using a Frobenius series (Leaver's method) around the horizon.
   • Expresses the bilinear integrals in terms of Tricomi hypergeometric functions ($U(a,b,z)$).
   • Computes the integral in a region where it converges ($\text{Re}(z) > 0$) and analytically continues the result to the physical QNM frequencies ($\text{Im}(\omega) < 0$).
2. Complex Contour Integration:
   • Deforms the radial integration path into the complex plane.
   • The contour runs parallel to the branch cut emanating from the horizon, ensuring the integrand decays exponentially at infinity.
   • This method avoids series expansions and relies on the damping provided by the complex path.

3. Key Contributions

• Geometric Interpretation of $J$: The paper provides a clear geometric definition of the $J$ operator within the hyperboloidal framework, showing it as a map between future and past foliations. This explains why the integrand diverges: it evaluates a regular mode against a singular anti-mode in the same coordinate system.
• Equivalence of Currents: The authors prove that the orthogonality relations derived from the symplectic current and the energy current yield consistent results when regularized, bridging a gap in the literature regarding the choice of current for QNM projections.
• Flux Cancellation: They demonstrate that the divergent bulk integrals and the boundary flux terms cancel each other out exactly when using regularization schemes, validating the use of bulk-only integrals for practical calculations.
• Excitation Coefficients: The paper establishes a rigorous link between orthogonality projections and Green's function techniques for calculating QNM excitation coefficients (amplitudes) directly from initial data.

4. Results

• Orthogonality Verification: Numerical tests on scalar perturbations of the Schwarzschild spacetime confirm that distinct QNMs are orthogonal ($\langle \Psi_I, \Psi_J \rangle = 0$ for $I \neq J$) to machine precision using both regularization methods.
• Convergence:
  • The semi-analytic method is robust for high overtones (high $n$) but computationally expensive due to series summation.
  • The complex contour method is highly efficient for low overtones but requires careful contour selection for modes with small real frequencies.
• Excitation Factors: The authors computed excitation factors ($b_n$) and amplitudes ($c_n$) for constant initial data. The results from both currents and both regularization methods agree within 6 significant digits, matching results from alternative hyperboloidal projection methods that do not require regularization.
• Limitations of Extended Products: While the "extended product" (including fluxes) is formally self-adjoint, it cannot be used to compute excitation coefficients directly from initial data alone because the flux terms require knowledge of the full time evolution of the system.

5. Significance

This work resolves a long-standing foundational issue in black hole perturbation theory: the definition of a well-behaved inner product for QNMs in a geometrically natural setting.

• For Gravitational Wave Astronomy: It provides a robust, geometrically informed framework for calculating QNM excitation amplitudes, which are crucial for interpreting gravitational wave signals from black hole mergers (ringdown phase).
• Theoretical Physics: It clarifies the role of boundary conditions and fluxes in non-Hermitian spectral problems, showing how conformal methods can tame singularities inherent in QNM problems.
• Future Applications: The framework is generalizable to rotating (Kerr) black holes and other perturbation types (vector/tensor), paving the way for precision black hole spectroscopy in next-generation detectors.