Gravitational wave scattering at $\mathcal{O}(G^4)$: Murua construction and elliptics

Technical Summary: Gravitational Wave Scattering at O(G⁴): Murua Construction and Elliptics

Problem and Context
This work addresses the theoretical modeling of gravitational wave scattering off a Schwarzschild black hole (BH) background, a process critical for validating Effective Field Theory (EFT) descriptions of compact objects in General Relativity (GR). While EFT techniques model compact bodies as point particles corrected by higher-dimensional operators, these models must be matched against exact GR solutions (specifically Black Hole Perturbation Theory, or BHPT) to fix Wilson coefficients and validate the point-particle approximation. Previous work established the equivalence between spinless point particles and Schwarzschild BHs up to the third post-Minkowskian (3PM, $O(G^3)$) order. This paper extends that validation to the fourth order, $O(G^4)$ (or 4PM), requiring a three-loop calculation. A key challenge at this order is the extraction of the physical scattering phase shift from the scattering amplitude, a process that traditionally involves complex cut-subtraction procedures to handle infrared divergences and causality prescriptions.

Methodology
The authors employ Worldline Quantum Field Theory (WQFT) to compute the scattering amplitude of a gravitational wave off a spinless point particle. The calculation proceeds through the following steps:

1. Framework: The scattering is modeled using the WQFT action for a massive point particle coupled to linearized gravity. The amplitude is computed as an S-matrix element expanded in powers of the perturbative parameter $\epsilon_{PM} \equiv 2Gm\omega$.
2. Murua Construction: A central methodological innovation is the direct extraction of the "Magnusian" (the matrix element of the Magnus operator $\hat{N}$, where $\hat{S} = e^{i\hat{N}}$) from the master integrals. Instead of calculating the scattering amplitude $\hat{T}$ and performing explicit cut subtractions to convert it to $\hat{N}$, the authors apply Murua coefficients [42–44] directly to the master integrals.
   • This approach assigns specific weights to integrals based on the number of "active" propagators that can go on-shell.
   • For an integral with $n$ active propagators, the Murua value is a weighted sum over all possible causality prescriptions (retarded/advanced), effectively averaging over them to produce a real, infrared-finite result directly.
3. Integration and Differential Equations: The three-loop calculation involves 70 diagrams reduced to a basis of 20 master integrals. The authors utilize the method of differential equations, bringing the system to a canonical form.
   • A significant technical feature at this order is the emergence of an elliptic sector. The canonical differential equation includes a Picard-Fuchs equation solved by complete elliptic integrals of the first kind, $K(x^2)$.
   • Boundary conditions are determined by expanding near $x=0$ (forward limit) and $x=1$ (backward limit), utilizing techniques such as loop-by-loop integration, parametric representations, and high-precision numerical reconstruction via AMFlow.

Key Results

1. Validation of Point-Particle Equivalence: The authors explicitly project the computed Magnusian onto a basis of spin-weighted spherical harmonics and match the result to the phase shift predicted by BHPT. Numerical checks for angular momentum modes $\ell \in [2, 20]$ confirm equality with a relative precision better than $10^{-34}$. This proves that the point-particle WQFT description accurately reproduces Schwarzschild BH physics up to $O(G^4)$.
2. Efficiency of Murua Coefficients: The Murua construction successfully bypasses the need for explicit cut subtractions. The resulting Magnusian is manifestly real and infrared-safe, avoiding the complex pole structures and regularization scheme dependencies inherent in the raw scattering amplitude.
3. Mathematical Structure:
   • Elliptic Functions: This is the first instance where elliptic functions appear in the momentum-space description of gravitational wave scattering at this order.
   • Function Space: While the scattering amplitude contains iterated integrals of weight three, the Magnusian is restricted to weight two and below. This suggests a simplification in the function space of the physical observable compared to the amplitude.
   • Infrared Behavior: The Magnusian remains finite in the forward limit ($x \to 0$) and backward limit ($x \to 1$), whereas the amplitude exhibits quadratic divergences and infrared poles that must be resummed (e.g., via the Newtonian phase).

Significance and Claims
The paper claims to advance the program of matching WQFT to BHPT by one order in Newton's constant, solidifying the indistinguishability of Schwarzschild black holes and spinless point particles up to $O(G^4)$. The authors highlight two primary contributions:

• Methodological: The demonstration that Murua coefficients provide a "drastically improved" and "immensely economic" method for extracting the Magnusian, effectively replacing the cumbersome cut-subtraction machinery used in previous work.
• Physical: The confirmation that no new non-minimal operators (which would describe finite-size effects or tidal deformations) are required to describe Schwarzschild BHs at this precision, as the point-particle limit suffices.

The authors note that while elliptic functions appear here, the function space remains relatively simple compared to the two-body scattering problem at the same order (where elliptic functions appear quadratically). They suggest that the simplicity of the function space in WQFT may allow for all-loop statements regarding the structure of gravitational wave scattering. The work sets the stage for future calculations at $O(G^5)$, where non-minimal operators will become relevant, and for the inclusion of spin effects to model Kerr black holes.