Calculation of a regularized Teukolsky Green function in Schwarzschild spacetime

This paper derives exact expressions for the factors in the Hadamard form of the retarded Green function for the Teukolsky equation on Schwarzschild spacetime by utilizing a conformal direct-product metric to obtain a separable direct part and its multipolar modes, ultimately providing an improved representation of the Green function near coincidence for gravitational perturbations.

The Gist

Imagine you are standing in a vast, dark ocean (the universe) near a massive whirlpool (a black hole). If you drop a pebble, it creates ripples. In physics, we call these ripples "fields" (like light or gravity). To understand how these ripples behave, scientists use a mathematical tool called a Green's Function. Think of this function as a "ripple map" that tells you exactly how a disturbance at one point in space and time affects every other point.

However, calculating this map near a black hole is incredibly difficult. It's like trying to draw a map of a storm where the center is a singularity—a point of infinite chaos. The math breaks down right at the source of the ripple because the numbers go to infinity.

This paper, by David Aruquipa, Marc Casals, and Brien Nolan, is about creating a better, more accurate "ripple map" for black holes, specifically for different types of waves (light, gravity, and scalar fields).

Here is the breakdown of their work using simple analogies:

1. The Problem: The "Infinite Spike"

When you try to calculate the ripple map, you run into a "spike" of infinite energy right where the wave starts.

• The Old Way: Previously, scientists tried to fix this by breaking the map into tiny pieces (like pixels) and adding them up. But this is like trying to draw a sharp spike using only smooth, round dots. The result is a blurry blob instead of a sharp spike. This "blurring" ruins the accuracy of the calculation, especially when you need to know exactly what happens right next to the source.
• The Goal: They wanted to separate the "infinite spike" (which is predictable and messy) from the "smooth part" (which is the interesting physics we actually want to study).

2. The Solution: The "Hadamard Recipe"

The authors use a mathematical recipe called the Hadamard form. Imagine the ripple map is a cake.

• The Direct Part (The Crust): This is the "spiky" part. It only exists exactly where the wave travels in a straight line (or a light-like path). It's the messy, infinite part.
• The Tail Part (The Filling): This is the smooth part. It represents the wave bouncing around the black hole, lingering in the past. This is the part that carries the complex information about the black hole's gravity.

The authors' main achievement is figuring out the exact recipe for the "Crust" (the Direct Part) for the first time for spinning fields (like light and gravity) on a black hole.

3. The Trick: Unfolding the Black Hole

To solve this, they used a clever geometric trick.

• The Analogy: Imagine the black hole's spacetime is a complex, crumpled piece of paper. It's hard to draw on.
• The Trick: They realized this paper could be "unfolded" into two simpler sheets stuck together:
  1. A flat, 2D sheet representing time and distance from the center ($M_2$).
  2. A sphere representing the angles around the black hole ($S_2$).
• By separating the problem into these two sheets, they could solve the math for each sheet individually and then stitch them back together.

4. The "Euler Angle" Dance

One of the most beautiful parts of their work involves the sphere ($S_2$).

• To describe how a wave moves across the sphere, they had to track how the wave "rotates" as it travels.
• They discovered a deep connection between the path the wave takes (a geodesic) and Euler angles (a way to describe 3D rotation, like how a gymnast twists in the air).
• The Metaphor: Imagine a dancer spinning on a stage. The authors found a formula that tells you exactly how much the dancer's "spin" changes based on how far they walked and how they turned. This allowed them to write down the "Direct Part" of the ripple map in a clean, exact formula.

5. The Result: A Sharper Map

Once they had the exact formula for the "Crust" (the Direct Part), they did something very smart:

1. They calculated the "Crust" exactly.
2. They subtracted it from the total "Ripple Map."
3. What was left was the "Filling" (the Non-Direct Part).

Why is this better?
Because they removed the messy, infinite "Crust," the remaining "Filling" is much smoother and easier to calculate.

• Before: If you tried to calculate the ripple map 1 second after the pebble dropped, the old method was blurry and inaccurate.
• Now: Their new method gives a sharp, accurate picture much closer to the moment the pebble was dropped (down to about 1 second instead of 4 or 6).

6. Why Does This Matter?

This isn't just about math for math's sake.

• Self-Force: If a small object (like a star or a planet) orbits a black hole, it creates its own ripples. These ripples push back on the object, changing its orbit. This is called "self-force." To calculate this, you need a perfect map of the ripples right next to the object.
• Gravitational Waves: As we detect more black hole collisions with telescopes, we need incredibly precise models to understand what we are hearing. This paper provides a better "decoder ring" for those signals.

Summary

The authors took a messy, impossible-to-calculate problem (ripples near a black hole), unfolded the universe into two simpler shapes, found a dance-like formula for how the waves rotate, and used that to surgically remove the "infinite noise." The result is a much clearer, more accurate view of how gravity and light behave near the edge of a black hole.

Technical Summary

1. Problem Statement

Calculating the retarded Green function (GF) for field perturbations (scalar, electromagnetic, and gravitational) on black hole spacetimes is a fundamental challenge in general relativity, particularly for self-force calculations.

• The Difficulty: The GF is a distribution with singular support at the coincidence point ($x=x'$) and along null geodesics connecting the points. Standard numerical approaches using multipolar ($\ell$-mode) decompositions suffer from "smearing" the Dirac-$\delta$ singularity into Gaussian peaks. This contamination prevents accurate calculation of physical quantities (like self-force) near the singularity.
• The Goal: To improve the representation of the GF near coincidence by separating the singular "direct" part (governed by the Hadamard form) from the regular "tail" (non-direct) part. The authors aim to derive exact expressions for the direct part of the Bardeen-Press-Teukolsky (BPT) equation for generic spin $s$ fields, enabling the subtraction of the direct modes from the full GF modes to yield a more accurate non-direct part.

2. Methodology

The authors employ a 2+2 conformal decomposition of the Schwarzschild spacetime, leveraging the fact that Schwarzschild is conformal to the direct product of a 2-dimensional Lorentzian spacetime ($M_2$) and a 2-sphere ($S^2$).

• Conformal Rescaling: They define a conformally rescaled metric $\hat{g}_{\mu\nu} = r^{-2}g_{\mu\nu}$ on $\hat{M} = M_2 \times S^2$. This allows the Teukolsky operator to be separated into a radial-temporal part ($M_2$) and an angular part ($S^2$).
• Hadamard Form Decomposition: The direct part of the Green function, $sU(x, x')$, satisfies a transport equation. The authors factorize this direct biscalar into two independent factors:
  $$s\hat{U}(x, x') = s\bar{U}(x_i, x_i') \cdot s\breve{U}(x_A, x_A')$$
  where $s\bar{U}$ is the factor on $M_2$ and $s\breve{U}$ is the factor on $S^2$.
• Analytic Derivation:
  • $S^2$ Factor: They solve the transport equation on the sphere explicitly. This reveals a deep connection between the transport equation, geodesics on $S^2$, and Euler angles ($\bar{\alpha}, \bar{\beta}$). They derive a closed-form geometric expression involving these angles.
  • $M_2$ Factor: They express the $M_2$ factor in terms of the square root of the van Vleck determinant (known for the scalar case) multiplied by a spin-dependent coefficient. This coefficient is calculated as an integral along geodesics in $M_2$.
• Mode Decomposition: Using the derived analytic forms, they compute the multipolar $\ell$-modes of the direct part. They utilize spin-weighted spherical harmonics (SWSH) and addition theorems to project the geometric direct term onto $\ell$-modes.
• Regularization: They subtract these calculated direct $\ell$-modes from the full GF $\ell$-modes (previously calculated in [13]) to obtain the non-direct part ($sG_{nd}$).

3. Key Contributions

1. Exact Geometric Expression for Direct Term: The paper provides the first fully geometrical, analytic expression for the direct Hadamard term of the BPT equation for arbitrary spin $s$ on Schwarzschild spacetime.
   • The $S^2$ factor is given in closed form as $e^{-is(\bar{\alpha}+\bar{\beta})}(\gamma/\sin\gamma)^{1/2}$, where $\gamma$ is the geodesic distance and $\bar{\alpha}, \bar{\beta}$ are specific Euler angles.
   • The $M_2$ factor is expressed via the van Vleck determinant and an integral along geodesics involving the spin parameter.
2. Closed-Form Solutions for Specific Worldlines: Exact solutions for the direct term are derived for:
   • Circular timelike geodesics (constant radius).
   • Radial null geodesics.
   • These solutions show exponential growth/decay depending on the spin $s$ and the radial position relative to the photon ring ($r=3M$).
3. Multipolar Modes of the Direct Part: The authors derive an analytic formula (Eq. 5.20) for the $\ell$-modes of the direct part, $sG^d_\ell$, expressed in terms of Jacobi polynomials and the $M_2$ transport coefficients.
4. Regularized Green Function Calculation: By subtracting the direct modes from the full modes, they obtain a "non-direct" Green function that is significantly more accurate near coincidence than the raw mode-sum.

4. Key Results

• Angular Factor: The angular dependence of the direct term is explicitly linked to Euler angles, showing that the spin-weighted nature of the field introduces a phase factor $e^{-is(\bar{\alpha}+\bar{\beta})}$ dependent on the geodesic path on the sphere.
• Spin-Dependent Behavior: For non-zero spin ($s \neq 0$), the direct term includes a spin-dependent prefactor that grows or decays exponentially with time, unlike the scalar case ($s=0$).
• Numerical Validation:
  • The authors validated their exact elliptic integral representations against small coordinate expansions and near-coincidence expansions.
  • They demonstrated that the small coordinate expansion is valid up to $\Delta t \approx 20M$ for a circular orbit at $r=6M$.
• Improved Non-Direct Part:
  • For gravitational perturbations ($s=-2$), the subtraction of the direct part extends the region of validity of the mode-sum calculation closer to coincidence.
  • Circular Orbit ($r=6M$): Validity improves from $\Delta t \approx 4M$ to $\Delta t \approx 1.6M$.
  • Static Worldline ($r=6M$): Validity improves from $\Delta t \approx 2.9M$ to $\Delta t \approx 1.0M$.
• Limitations: The non-direct part still does not perfectly represent the GF exactly at coincidence ($\Delta t = 0$) because the subtraction leaves a "smeared" step-function behavior due to the finite mode sum. However, it is a superior approximation for the "intermediate" time regime crucial for self-force integrals.

5. Significance

• Self-Force Calculations: This work is critical for the "mode-sum" regularization scheme used in calculating the gravitational self-force on particles orbiting black holes. By providing an exact analytic form for the singular direct part, it allows for a cleaner subtraction, reducing numerical errors in the self-force calculation which is highly sensitive to the behavior of the Green function near the particle.
• Generalization: While previous work (e.g., Casals et al. [8]) handled the scalar case ($s=0$), this paper generalizes the method to generic spin fields (electromagnetic and gravitational), which are physically more relevant for astrophysical black hole perturbations.
• Geometric Insight: The derivation highlights a novel interplay between the geometry of geodesics on $S^2$, Euler angles, and the transport equations of the Hadamard coefficients, offering a deeper geometric understanding of the Green function's structure.
• Future Applications: The exact expressions for the direct term can be used to model the early stages of gravitational wave ringdowns and serve as a benchmark for numerical relativity codes. The authors note that while the non-direct part is improved, a separate calculation of the Hadamard tail ($sV$) near coincidence is still required for a complete solution, which they leave for future work.