Nonlocal-in-time tail effects in gravitational scattering to fifth Post-Minkowskian and tenth self-force orders

Using worldline effective field theory and the novel Sparse Integral Reducer (SpideR) integration algorithm, this paper derives the nonlocal-in-time conservative tail effects in gravitational scattering up to fifth Post-Minkowskian and tenth self-force orders, expressing the results in terms of multiple polylogarithms and confirming consistency with existing literature through sixth post-Newtonian order.

The Gist

Imagine two massive objects, like black holes or neutron stars, zooming past each other in the vast emptiness of space. They don't crash; they just swing around one another, like two skaters passing on ice, before flying off in different directions. This is called "gravitational scattering."

For decades, physicists have been trying to write the perfect "rulebook" for how these objects move. This rulebook needs to be incredibly precise because modern telescopes (like LISA) are about to listen to the faint whispers of these cosmic dances. To do this, scientists use a method called Post-Minkowskian (PM) expansion, which is like building a model of gravity layer by layer, adding more and more detail with each step.

This paper, written by a team of physicists, tackles a very specific, tricky layer of that rulebook: the 5th layer (5PM).

The Problem: The "Echo" Effect

When these two heavy objects move, they don't just move through empty space; they ripple the fabric of spacetime itself, sending out gravitational waves (ripples).

Here's the catch: These ripples don't just fly away forever. Some of them bounce off the "static" gravity field created by the objects themselves. Think of it like shouting in a canyon. You shout, the sound hits the walls, and an echo comes back to you.

In physics, this echo is called a "tail effect." It's a "nonlocal-in-time" effect, which is a fancy way of saying: What happens right now depends on what happened in the past. The objects are reacting to their own echoes.

The problem is that these echoes make the math incredibly messy. If you try to use the rules for flying objects (scattering) to predict how objects orbit each other (bound systems), these echoes cause the math to break down unless you carefully separate the "instant" effects from the "echo" effects.

The Solution: A New Mathematical Tool

To solve this, the authors had to calculate the exact size and shape of these "echoes" up to a very high level of precision (the 10th order of a specific mass-ratio expansion).

The math involved was so complex that standard computer tools couldn't handle it. It was like trying to solve a Sudoku puzzle where the grid suddenly grew from 9x9 to a size that would fill a stadium.

So, the team built a new digital tool called SPI
D
E
R (Sparse Integral Reducer).

• The Analogy: Imagine you have a massive pile of tangled headphones. A standard tool tries to untangle them one by one, which takes forever. SPI
  D
  E
  R is like a smart robot that looks at the whole pile, figures out the pattern of the tangles, and creates a set of instructions to untangle any knot in that pile instantly. It uses a clever trick called "finite-field arithmetic" (doing math with remainders of prime numbers) to keep the numbers small and manageable before reconstructing the final answer.

What They Found

Using this new tool, the team successfully calculated the "tail" contributions to the gravitational scattering angle.

• The Result: They found that the math describing these echoes involves complex numbers called "multiple polylogarithms" (think of them as advanced, multi-layered logarithmic functions).
• The Check: They compared their results with existing, highly accurate calculations from a different approach (Post-Newtonian expansion) and found perfect agreement. This confirms their new tool and method are working correctly.

Why It Matters

The ultimate goal isn't just to calculate scattering; it's to understand binary inspirals (two objects spiraling into each other).

Currently, physicists have the "scattering" data (how they fly past) and the "bound" data (how they orbit), but they can't perfectly translate one into the other because of these "echo" effects. This paper provides the missing piece of the puzzle. By isolating and calculating the "echo" part, they have cleared the path to finally extract the "instant" part of the gravity rulebook.

Once they have that, they can use the "Boundary-to-Bound" dictionary (a mathematical translation tool) to turn the scattering data into a perfect model for orbiting black holes. This will help future gravitational wave detectors listen to the universe with unprecedented clarity.

In short: The authors built a super-smart calculator (SPI
D
E
R) to solve a massive math problem about gravitational echoes. They proved their method works by matching it with known results, and now they have the key ingredient needed to build a perfect map of how black holes dance together.

Technical Summary

1. Problem Statement

The paper addresses a fundamental conceptual obstruction in General Relativity: the difficulty of converting scattering data (hyperbolic encounters) into observables for gravitationally bound systems (elliptic orbits) using the Post-Minkowskian (PM) expansion.

• The Core Issue: While the "local-in-time" conservative dynamics of binary systems can be derived from scattering data via the "boundary-to-bound" (B2B) dictionary, the full scattering amplitude contains nonlocal-in-time hereditary effects. Specifically, tail effects arise when emitted gravitational waves scatter off the static curvature of the binary's own field.
• The Gap: Previous calculations at 5PM order (fifth order in Newton's constant $G$) were limited to the leading Self-Force (SF) order (1SF). To fully isolate the local-in-time conservative dynamics required for high-precision waveform modeling (especially for future detectors like LISA), one must compute these nonlocal tail contributions to higher orders in the mass ratio ($\nu = m_1 m_2 / (m_1+m_2)^2$).
• Objective: The authors aim to derive the nonlocal-in-time tail contributions to the gravitational scattering angle at 5PM order and through the 10th order in the Self-Force expansion (10SF).

2. Methodology

The authors utilize the Worldline Effective Field Theory (WEFT) formalism combined with advanced multiloop integration techniques.

A. Theoretical Framework

• Radial Action: The nonlocal contribution is governed by an integral over the source energy spectrum multiplied by a logarithm of the center-of-mass gravitational wave frequency.
• Expansion Strategy: The integrand is expanded in powers of the small mass ratio ($m_2/m_1 \ll 1$). This splits the calculation into three parts:
  1. A term independent of loop momenta (related to total radiated energy).
  2. A term treated as a standard linear propagator with non-integer power (previously solved at 5PM/1SF).
  3. The New Contribution: The remaining terms in the expansion, which involve increasingly high powers of linear propagators (up to $\pm 12$ at 10SF).

B. Computational Innovation: SPI

DE
R
The primary technical challenge is the "Integration-by-Parts" (IBP) reduction of scalar integrals. Standard algorithms (like Laporta's) become computationally intractable due to the massive index ranges of the propagators at high SF orders.

• Solution: The authors developed a novel algorithm called **SPI
  DE
  R (Sparse Integral Reducer).
• Mechanism: Unlike traditional methods that solve for a specific list of target integrals, SPI
  DE
  R constructs symbolic reduction rules at the level of "sectors."
  • It uses Finite-Field (FF) arithmetic and Rational Reconstruction (RR) to handle large expressions efficiently.
  • It avoids the "forward solving" bottleneck by applying pre-derived symbolic rules to reduce integrals on the fly.
  • This allows the computation of a reduction graph containing $\sim 10^8$ integrals entirely in main memory.

C. Integration and Boundary Conditions

• Master Integrals: The reduced integrals are solved using the method of differential equations, transformed into a canonical form.
• Special Functions: The solutions are expressed in terms of Multiple Polylogarithms (MPLs) up to transcendental weight three.
• Boundary Conditions: The integration constants are determined by taking the limit of small relative velocity ($v \to 0$), where the three-loop integrals factorize into a two-loop integral and a tadpole.

3. Key Contributions

1. Derivation of 5PM/10SF Tail Effects: The paper provides the first complete analytic derivation of nonlocal-in-time tail corrections to the scattering angle up to the 10th order in the mass-ratio expansion.
2. **Development of SPI
   DE
   R: A new, highly efficient IBP reduction tool capable of handling the extreme complexity of high-order SF expansions, overcoming the limitations of existing public codes (like Kira, Fire, LiteRed).
3. Analytic Results: The authors provide explicit expressions for the scattering angle coefficients involving rational functions and MPLs, valid for arbitrary mass ratios (via the $\nu$ expansion).
4. Validation: The results show perfect agreement with existing Post-Newtonian (PN) literature up to 6PN order in the overlapping regime.

4. Key Results

• Scattering Angle Structure: The nonlocal contribution to the scattering angle $\chi$ is expressed as a sum over powers of the impact parameter $b$ and the mass ratio $\nu$. The coefficients involve a basis of 26 special functions ($G_i(x)$) constructed from MPLs with an alphabet $\{0, 1, 2, 1 \pm i\}$.
• Explicit Expansion: The authors provide the PN-expanded form of the scattering angle (Eq. 9) up to 10PN order, including logarithmic terms and rational coefficients dependent on $\nu$.
• Boundary Integrals: The paper lists the analytic values of the required boundary integrals (Appendix A), which involve generalized hypergeometric functions (${}_4F_3$) and their $\epsilon$-expansions.
• Agreement: The computed nonlocal tail terms match the "W1-only" results from previous PN studies (Bini, Damour, Geralico) through 6PN, confirming the correctness of the 5PM/10SF derivation.

5. Significance and Implications

• Isolating Local Dynamics: The primary utility of this work is enabling the subtraction of nonlocal tail effects from the total 5PM scattering amplitude. This leaves behind the "local-in-time" conservative dynamics, which is the necessary input for constructing the Hamiltonian for bound binary inspirals.
• Waveform Modeling: These results are critical for third-generation gravitational wave detectors (LISA, Cosmic Explorer, Einstein Telescope), which require extreme precision (up to 5PM and high SF orders) to model the inspiral phase of massive black hole binaries.
• Resolution of Singularities: The paper discusses the "spurious poles" at $\gamma=3$ (where $\gamma$ is the Lorentz factor) that appear in potential-region calculations. By isolating the tail sector, the authors clarify the separation between conservative and dissipative sectors, offering a path to resolve ambiguities in the B2B dictionary (e.g., the "memory" contribution vs. potential contributions).
• Future Directions: The authors note that while the tail sector is now under control, the full 5PM conservative dynamics still requires a careful treatment of "memory-like" effects and the handling of singularities in the potential sector, which they plan to address in future work.

In summary, this paper represents a major leap in the precision of gravitational scattering calculations, bridging the gap between high-order PM scattering and high-order PN bound dynamics through the development of novel computational tools and the rigorous calculation of hereditary tail effects.