Gravitational Wave Scattering in Spinless WQFT

1. Problem Statement

The paper addresses the scattering of gravitational waves (GWs) off a black hole (BH), a process fundamental to understanding wave propagation in strong-field gravity and the tidal response of compact objects.

• Context: While Black Hole Perturbation Theory (BHPT) has long provided exact results for this scattering via the Teukolsky/Regge-Wheeler equations, there is a growing need to connect these results to Effective Field Theory (EFT) approaches, specifically Worldline Quantum Field Theory (WQFT).
• Motivation: WQFT models black holes as point particles. To accurately model binary black hole scattering at high Post-Minkowskian (PM) orders (beyond the current 5PM state-of-the-art), one must determine "non-minimal" operators in the WQFT action that encode finite-size effects (tidal Love numbers).
• The Gap: A unique matching calculation between WQFT amplitudes and General Relativity (GR) results (derived from BHPT) is required to fix the Wilson coefficients of these non-minimal operators. Previous work had computed WQFT amplitudes to Next-to-Leading Order (NLO) but lacked a direct comparison to BHPT phase shifts.
• Goal: This paper aims to compute the gravitational wave scattering amplitude in spinless WQFT up to Next-to-Next-to-Leading Order (NNLO, or $O(G^3)$) and demonstrate its equivalence to the BHPT phase shift, thereby validating the framework for future higher-precision matching.

2. Methodology

A. Theoretical Framework: WQFT vs. BHPT

• WQFT Approach: The black hole is modeled as a spinless point particle with a worldline $x^\mu(\tau)$. The interaction is described by the Einstein-Hilbert action coupled to the worldline. The scattering is computed using Feynman diagrams with retarded propagators (Schwinger-Keldysh formalism) to ensure causality.
• BHPT Approach: The scattering is analyzed in the rest frame of the black hole using the Teukolsky equation for curvature perturbations. The solution is expanded in spin-weighted spherical harmonics, yielding a partial-wave scattering matrix characterized by a phase shift $\delta$.
• The Matching Strategy: The authors do not match the raw scattering amplitudes directly. Instead, they utilize the exponential representation of the S-matrix:
  $$\hat{S} = e^{i\hat{N}}$$
  They prove that the matrix elements of the operator $\hat{N}$ (the Magnus operator) map directly onto the BHPT phase shift $\delta$ via a partial-wave transformation. This is crucial because $\hat{N}$ is infrared (IR) finite and exhibits better behavior in the forward limit compared to the standard T-matrix.

B. Computational Techniques

1. Diagram Generation:
   • The authors employ a recursive Berends-Giele-like relation derived from the classical limit of the WQFT action to generate the graviton two-point function (the source of the scattering amplitude).
   • They identify 6 one-loop and 20 two-loop diagrams.
   • A distinction is made between "potential" gravitons (off-shell, static) and "active" gravitons (on-shell, carrying the wave frequency $\omega$).
2. Integration:
   • Dimensional Regularization: Calculations are performed in $D = 4 - 2\epsilon$ dimensions.
   • Loop Reduction: Tensor integrals are reduced to scalar master integrals using Passarino-Veltman reduction.
   • Differential Equations: The master integrals are solved using the method of differential equations with respect to the scattering variable $x = \sin(\theta/2)$.
   • Boundary Conditions: The authors analyze the boundary behavior ($x \to 0$, forward limit) to determine integration constants. They identify two dominant kinematic regions (forward and general wave) that simplify the boundary conditions, reducing the number of independent constants required.
3. Extraction of $\hat{N}$:
   • The T-matrix elements are computed first.
   • The $\hat{N}$-matrix elements are extracted perturbatively using the relation $\hat{S} = 1 + i\hat{T} = e^{i\hat{N}}$, expanding order by order (e.g., $\hat{N}^{(2)} = \hat{T}^{(2)} - \frac{1}{2}(\hat{T}^{(1)})^2$).

3. Key Contributions

1. Proof of Exponentiation in Partial-Wave Space:
   The paper provides a rigorous proof that, in the absence of dissipation, the exponential representation of the S-matrix ($\hat{S} = e^{i\hat{N}}$) maps directly to the BHPT phase shift. This establishes that the WQFT amplitude itself exponentiates in partial-wave space, a non-trivial result linking quantum-inspired EFT techniques to classical GR perturbation theory.
2. NNLO ($O(G^3)$) Calculation:
   The authors explicitly compute the scattering amplitude for a spinless black hole up to two loops (NNLO). This is the first time such a high-order calculation has been matched directly to BHPT results for wave scattering.
3. Efficient Diagrammatic and Integration Techniques:
   • Development of a recursive generator for the graviton two-point function.
   • Demonstration that the forward-limit divergence of the T-matrix is regularized in the $\hat{N}$-matrix, allowing for a clean projection onto 4D spin-weighted spherical harmonics without needing higher-dimensional regularization for the matching step.
4. Validation of the Minimal WQFT Action:
   The calculation confirms that the minimal point-particle WQFT action (without non-minimal tidal operators) correctly reproduces the GR scattering amplitude up to $O(G^3)$. This implies that non-minimal operators (tidal effects) only enter at higher orders ($O(G^4)$ in the action, $O(G^6)$ in scattering), consistent with theoretical expectations.

4. Key Results

• Amplitude Expressions: The paper provides explicit analytic expressions for the $\hat{N}$-matrix elements for both helicity-preserving ($\sigma=1$) and helicity-reversing ($\sigma=-1$) scattering up to $O(\epsilon_{PM}^3)$.
  • The results involve polylogarithms ($Li_2$), logarithms, and rational functions of the scattering variable $x$.
  • For helicity-reversing scattering, the $O(G^2)$ term vanishes, and the $O(G^3)$ term is non-zero.
• Matching to BHPT:
  • By projecting the $\hat{N}$-matrix elements onto spin-weighted spherical harmonics, the authors recover the known BHPT phase shifts $\delta_{\ell}^{(n)}$ for $n=1, 2, 3$.
  • The matching is exact for the leading, NLO, and NNLO terms.
  • The regularization of the forward divergence (helicity-preserving LO) is shown to reproduce the BHPT result when the infrared scale is fixed appropriately.
• Gauge Invariance: The computed T-matrix element is verified to be gauge invariant, satisfying Weinberg's soft graviton theorem (exponentiation of IR divergences).

5. Significance and Outlook

• Foundation for Tidal Matching: This work lays the essential groundwork for determining the Wilson coefficients of non-minimal operators in the WQFT action. Since the minimal action reproduces BHPT up to $O(G^3)$, any deviation at higher orders (or for spinning black holes) can be uniquely attributed to finite-size effects (Love numbers).
• Spin Extension: The methodology is readily extensible to spinning black holes. The authors note that while non-minimal operators for spinless scattering appear at high PM orders, spin effects can be incorporated at lower loop orders, making the matching of spinning WQFT to the Teukolsky equation a viable next step.
• Computational Efficiency: The demonstration that the $\hat{N}$-matrix element is IR-finite and well-behaved in the forward limit suggests it is the superior observable for matching EFT to GR, avoiding the complexities of IR divergences inherent in the T-matrix.
• Future Directions: The authors plan to extend this matching to higher PM orders and to include spin, which will be crucial for the precision modeling of gravitational waveforms from binary black hole mergers detected by observatories like LIGO/Virgo and future space-based detectors.

In summary, this paper successfully bridges the gap between Worldline EFT and Black Hole Perturbation Theory, validating the WQFT approach for wave scattering and providing the necessary technical tools to extract tidal properties of black holes from gravitational wave data.