⚛️ high-energy theory

Recursion Relations for Classical Gravity

This paper derives and applies a new set of recursion relations to calculate the perturbative scattering of two black holes in Einstein gravity, successfully recovering results up to the third post-Minkowskian order that include gravitational radiation back-reaction.

Original authors: Poul H. Damgaard, Kwangeon Kim, Kanghoon Lee

Published 2026-07-15

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Poul H. Damgaard, Kwangeon Kim, Kanghoon Lee

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, invisible trampoline made of spacetime. Usually, we think of this trampoline as perfectly flat and calm. But when you drop two heavy bowling balls (black holes) onto it, they don't just sit there; they roll toward each other, warping the fabric as they go. This paper is about figuring out exactly how those two balls move and how the trampoline ripples back at them, using a brand-new set of mathematical "recursive" instructions.

The Main Discovery: A Step-by-Step Recipe
The authors, Poul H. Damgaard, Kwangeon Kim, and Kanghoon Lee, have cooked up a new way to calculate the dance of these two massive objects. Instead of trying to solve the entire messy equation of gravity all at once (which is like trying to eat a whole pizza in one bite), they created a recursive recipe.

Think of it like building a tower of blocks. You start with a flat floor (flat space). Then you add the first layer of blocks (the first approximation of gravity). Once that layer is done, you use it to build the second layer, then the third. The paper shows that you can keep stacking these layers higher and higher, and each new layer is built entirely from the pieces you already calculated for the layers below it.

They tested this recipe up to the third post-Minkowskian order. In plain English, this means they calculated the interaction with extreme precision, going beyond the simple "first guess" and even the "second guess," to include the third level of detail. At this specific level, something fascinating happens: the gravity doesn't just pull the balls together; it also creates ripples (gravitational waves) that travel away and then push back on the balls. The authors found that their recursive method perfectly captures this "back-reaction," matching the known results from other complex methods.

What They Are NOT Doing (and What They Rule Out)
It is crucial to understand what this paper is not about, because the authors are very specific.

No Black Holes (Yet): Even though the title mentions "scattering of two black holes," the authors explicitly state that at any finite step of their calculation, there are no black holes in the problem. There are no event horizons, no singularities. They are scattering two massive point-sources. The "black hole" only emerges if you could theoretically add up infinite layers of their recipe. For now, they are just calculating the path of two heavy points moving through a flat background that gets increasingly bumpy.
No Effective Field Theory: Many modern physicists use a tool called "Effective Field Theory" to simplify gravity problems. This paper explicitly rejects that approach. They argue that you don't need to throw away parts of the theory to get the classical answer. Instead, they go straight back to the original, raw equations of motion (Einstein's equations) and solve them directly. They believe this "old-school" approach, when combined with their new recursive tricks, is actually cleaner and more direct.
Not a Simulation: This isn't a computer simulation where they guessed the answer and checked if it looked right. They derived the equations mathematically and solved them analytically. They proved that if you follow their recursion steps, you get the exact same answer as the most trusted methods in the field.

The "Magic" Ingredient: Causality
One of the most playful and clever parts of their method is how they handle time. Gravity travels at the speed of light, so the effect of one ball on the other must happen after the cause. In math, this is called "causality."

The authors found that if you simply use retarded Green functions (a fancy math term for "solutions that respect the arrow of time"), the universe automatically does the hard work for you. You don't need to manually subtract out "super-classical" nonsense or worry about energy loss. By strictly enforcing that effects happen after causes, the math naturally includes the energy lost to gravitational waves (the dissipative part) and the energy that stays in the system (the conservative part). It's as if they built a time-machine into the equations, and the universe politely handed them the correct answer for the energy loss.

The Verdict
The paper is a proof of concept. The authors have successfully demonstrated that this recursive method works up to the third order of precision. They have recovered the correct scattering angle and momentum kick, including the subtle effects of radiation.

They are confident that this method works for these specific orders. However, they are honest about the future: they admit that going to the fourth order (and beyond) will be significantly harder because the integrals (the math sums) become much more complicated. They haven't solved the fourth order yet, but they have shown that the door is open and the path is clear.

In short, they took a problem that was thought to be too messy for direct calculation, gave it a new set of algebraic tools, and showed that with a little patience and a strict respect for the flow of time, you can build the solution block by block, exactly as nature intended.

Technical Summary: Recursion Relations for Classical Gravity

Problem Statement
The paper addresses the perturbative evaluation of the scattering of two massive bodies in Einstein gravity, specifically within the post-Minkowskian (PM) expansion. The goal is to solve Einstein's field equations order-by-order in Newton's constant $G$ to determine physical observables such as the scattering angle and momentum kick. While recent amplitude-based and worldline approaches have advanced high-order calculations, this work returns to the classical equations of motion to establish a systematic, algebraic recursive formalism. A primary challenge addressed is the inclusion of dissipative effects (gravitational radiation back-reaction) at the third post-Minkowskian (3PM) order, a regime where the choice of Green function prescriptions becomes critical for capturing causality and energy loss.

Methodology
The authors employ a perturbative expansion of the Einstein-Hilbert action coupled to point-like matter sources, utilizing a specific set of variables and algebraic structures to manage complexity:

Gothic (Landau-Lifshitz) Variables and Doubled Formalism: The metric is expanded using "gothic" variables ( $\mathfrak{g}^{\mu\nu} = \sqrt{-g}g^{\mu\nu}$ ) and an auxiliary inverse field. This "doubled" set of variables avoids the infinite geometric series expansion typically required when inverting the metric, simplifying the interaction vertices.
Recursive Currents: The solution is formulated in terms of currents in Fourier space. The authors introduce a "doubled" Fourier transform structure, assigning separate momentum variables ( $\ell_1, \ell_2$ ) to each of the two scattering bodies. This allows the metric perturbations and matter trajectories to be expressed as convolutions of these currents.
Algebraic Recursion Relations: By defining a specific convolution bracket notation, the authors derive algebraic recursion relations for:
- Graviton Currents: Derived from the Einstein equations, relating the $n$ -th order metric perturbation to lower-order currents and the energy-momentum pseudotensor.
- Matter Currents: Derived from the geodesic equation, determining the trajectory corrections at each order based on the metric computed at previous orders.
Iterative Loop Integrals: The momentum integrals arising from the convolution brackets are evaluated systematically. The authors classify these into three families (Family-I, Family-II, and Radiation integrals) based on the structure of delta functions (velocity cuts) in the integrands. They utilize the method of regions (potential and radiation) and integration-by-parts (IBP) reductions to master integrals, often evaluating them as iterations of one-loop integrals.
Causal Prescription: Crucially, the formalism employs retarded Green functions. The authors argue that this choice automatically enforces causality and correctly captures all dissipative effects (radiation back-reaction) without needing to manually subtract super-classical terms, a feature often required in amplitude-based approaches.

Key Contributions and Results

Derivation of Recursion Relations: The paper provides a complete set of algebraic recursion relations for both graviton and matter currents up to arbitrary order in $G$ . This formalism allows the solution of the coupled Einstein-geodesic system without recourse to Feynman diagrams.
1PM and 2PM Verification: The formalism is tested at the first (1PM) and second (2PM) post-Minkowskian orders. The results for the momentum kick and scattering angle are shown to reproduce the known Westpfahl result (1985) and subsequent modern calculations, validating the recursive approach.
3PM Calculation with Radiation Back-Reaction: The primary technical achievement is the explicit calculation of the 3PM momentum kick. This order is significant as it is the first where gravitational radiation reaction becomes non-negligible.
- The authors successfully separate the conservative and dissipative parts of the momentum kick.
- The calculation involves evaluating complex two-loop integrals (Family-II and "mushroom" integrals) arising from the radiation sector.
- The final result for the 3PM momentum kick, including the $m_1^2 m_2^2$ terms responsible for radiation reaction, is found to be in full agreement with existing literature (e.g., references [66], [62], [63]).
Role of Retarded Green Functions: The paper demonstrates that using retarded Green functions is sufficient to capture the full dissipative dynamics, including the subtle radiation back-reaction, directly from the equations of motion.

Significance and Claims
The authors claim that this work revitalizes the original program of solving the classical equations of motion for the two-body problem, which had been largely superseded by amplitude-based methods due to computational complexity.

Alternative Framework: The paper presents a "third way" for computing classical gravity observables, distinct from both the standard amplitude-based KMOC formalism and the worldline effective field theory approaches.
Simplicity and Automation: By relying on algebraic recursion and avoiding effective field theory subtractions, the method offers a conceptually direct path to high-order results. The authors note that the recursive nature makes the equations amenable to symbolic manipulation programs.
Convergence and Validity: The paper acknowledges that while the perturbative expansion for a single body (Schwarzschild) is known to converge to the non-perturbative solution, the convergence for the two-body problem is not guaranteed. However, the successful recovery of the 3PM result, including dissipative effects, suggests that the expansion provides a valid approximation for scattering far from the horizon scales.
Modest Scope: The authors explicitly state that the current work is limited to the 3PM order. They note that higher orders (4PM and beyond) will require significantly more involved integration techniques, similar to the challenges faced in other approaches. The paper does not claim to solve the full non-perturbative two-body problem but rather demonstrates the viability of the recursive classical approach up to the order where radiation reaction first appears.

More like this