🔢 mathematics

Gravity between two masses beyond Newton's law

This paper resolves the long-standing challenge of calculating the static gravitational interaction between two arbitrary masses beyond Newton's law by deriving high-order corrections and a closed-form expression for the first backreaction, utilizing a stationary interaction energy framework that bridges the gap between test-body limits and comparable mass systems.

Original authors: Zhengwen Liu, Xing-Bo Yuan

Published 2026-10-30

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Zhengwen Liu, Xing-Bo Yuan

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

Gravity is the force that keeps our feet on the ground and the planets in their orbits. For centuries, we understood this force through the lens of Isaac Newton, who described it as a simple pull between two objects that grows stronger as they get closer. This view worked perfectly for most of human history. However, in the early twentieth century, Albert Einstein revealed a deeper truth: gravity is not just a force, but a curvature of space and time itself. Massive objects do not merely pull on one another; they bend the fabric of the universe, and that bending tells other objects how to move.

This shift from a simple pull to a curved geometry introduced a profound complication. In Newton's world, gravity was additive; you could simply add up the effects of different masses. In Einstein's universe, gravity creates its own gravity. The energy of the gravitational field itself acts as a source of further curvature, creating a complex, self-interacting web. This nonlinearity means that calculating the interaction between two massive objects is incredibly difficult. Even when the objects are sitting still, not moving or radiating energy, the mathematics of their interaction involves an infinite hierarchy of feedback loops. For over a century, physicists have struggled to untangle these loops to predict exactly how two masses influence each other in this curved environment, a task that has become urgently important as we now listen to the gravitational waves of colliding black holes and neutron stars.

A team of researchers has now taken a massive leap forward in solving this static puzzle. They have calculated the precise gravitational interaction between two point masses, determining the strength of their pull with a level of detail that reaches the twelfth order in Newton's constant beyond the basic Newtonian description. While previous work had managed to map out the interaction up to the eighth order of complexity, this new study pushes the boundary to the twelfth, covering a vast range of mass ratios from a tiny probe orbiting a giant star to two objects of equal size. The researchers achieved this by focusing on the stationary energy of the system, a method that allowed them to bypass the need to solve for every single detail of the gravitational field. Instead, they used a mathematical principle where the stability of the system's energy reveals the necessary information, effectively filtering out the noise and leaving only the essential interactions.

One of the most significant achievements of this work is the derivation of a single, closed formula that describes the first major correction to the simple "probe" scenario. In the past, if one mass was tiny compared to the other, physicists could use a known solution called the Schwarzschild metric to describe the interaction perfectly. However, when the smaller mass begins to contribute its own gravity, it alters the geometry, creating a "backreaction." The team found that this first backreaction, which occurs when the smaller mass starts to warp the space around the larger one, can be expressed as a neat, exact function. This formula captures the entire infinite series of corrections for this specific effect, rather than just a few terms of an approximation. It reveals that as the objects get closer and the gravitational field gets stronger, the interaction changes in a predictable, smooth way, eventually crossing a point where the correction flips from positive to negative.

The researchers also mapped out the complete interaction for objects of any mass ratio, not just when one is tiny. They discovered that as the masses become more comparable, the interaction becomes dominated by higher-order effects where both objects are warping space simultaneously. In the case of two equal masses, the contribution from these complex, mutual warping effects is enormous, dwarfing the simple pull and the first backreaction. The team found that the coefficients of these terms grow rapidly, suggesting that the interaction becomes increasingly complex as the objects approach each other. Despite this complexity, the underlying mathematics remained surprisingly orderly, with the results coming out as precise rational numbers rather than messy, irrational constants. This suggests a hidden simplicity within the chaotic nonlinear structure of Einstein's theory.

The study relied on a combination of analytical proofs and high-precision numerical calculations to verify their results. They cross-checked their findings using three different independent methods: expanding the solution around one object, solving the equations for two objects on a grid, and evaluating the interactions through a series of mathematical pairings. These methods agreed with each other to an extraordinary degree of precision, confirming the validity of their calculations up to the twelfth order. The work bridges the gap between the simple case of a tiny probe and the complex reality of two comparable masses, providing a complete map of the static gravitational potential. This map serves as a fundamental building block for understanding the dynamics of binary systems, offering the precise theoretical predictions needed to interpret the signals from gravitational wave detectors and to model the orbits of black holes and neutron stars as they spiral toward each other.

Technical Summary: Gravity Between Two Masses Beyond Newton's Law

Problem Statement
Newton's law of universal gravitation describes the interaction between two masses via a simple inverse-distance potential. In Einstein's General Relativity, gravity is non-linear: mass curves spacetime, and the geometry itself acts as a source for further curvature. Even in the static limit (where bodies are stationary and non-radiating), the interaction between two masses involves an infinite hierarchy of field couplings. While the interaction for a test mass (a probe) in the background of a single mass is known exactly (the Schwarzschild solution), extracting the complete static interaction for two comparable masses remains a fundamental analytical challenge. Previous high-order calculations have been limited by the computational cost of advancing one power of Newton's constant ( $G$ ) at a time, typically reaching only up to $G^8$ with full mass dependence.

Methodology
The authors formulate the problem around the stationary energy of the gravitational field. Their approach relies on three core theoretical mechanisms:

Auxiliary Linear Sources: Each point mass is replaced by an auxiliary linear source. The stationary energy of the field is computed as a double power series in these sources, with coefficients given by dimensionless connected correlators ( $\bar{\Gamma}_{a,b}$ ). The physical couplings are restored algebraically at the end via a source-to-mass transformation, avoiding new spatial integrations.
Variational Stationarity: The authors exploit the principle that a variation of the fields about their classical solution alters the energy only at second order. Consequently, field responses known through order $p$ in a source determine the interaction energy through order $2p+1$ . This "2n+1 principle" allows the authors to eliminate the need to solve for the highest-order field responses explicitly; instead, the energy is determined by pairings of lower-order field responses.
Complementary Expansions: The calculation utilizes two overlapping expansion strategies:
- Exact Background Expansion: Expanding about the exact Schwarzschild field of one body sums all powers of that body's source. This yields the probe term (0SF) and the first backreaction (1SF) to all orders in $G$ .
- Flat Space Expansion: Expanding about flat space while retaining full dependence on both masses resolves the complete mass dependence at finite orders in $G$ .
- Dimensional Regularization: The calculation is performed in $d = 3 + \epsilon$ spatial dimensions to handle point-mass singularities. The authors carefully track the cancellation of poles and the emergence of finite local terms, utilizing exact $d$ -dimensional data for regular parts of the fields where necessary.

Key Contributions and Results

Exact First Backreaction (1SF):
The authors derive a closed-form expression for the first backreaction (the term linear in the mass ratio $q = m_2/m_1$ ) to all orders in $G$ .
- The potential is expressed as $V/(m_2 c^2) = \Phi_0(\alpha) + q\Phi_1(\alpha) + O(q^2)$ , where $\alpha = Gm_1/(rc^2)$ .
- $\Phi_1(\alpha)$ is given explicitly in Equation (3) of the paper.
- This result proves that all static 1SF coefficients are purely rational numbers, free of transcendental constants, at every order in $G$ .
- The derivation relies on a "collective cancellation" where the sum of subtracted generic angular modes (multipoles) vanishes, leaving only contributions from the monopole and dipole.
Static Potential to Twelfth Order ( $G^{12}$ ):
The authors determine the complete static interaction between two arbitrary point masses through the 12th order in Newton's constant ( $G^{12}$ ).
- This extends the previous frontier (which was $G^8$ with full mass dependence) by four orders.
- The result is presented as a polynomial in the masses and $G$ , with coefficients that are rational numbers.
- The calculation required evaluating 21 independent connected correlators. Eleven were fixed by the exact 0SF and 1SF functions; three were known from prior work; and seven were new.
- Of the seven new correlators, four were evaluated analytically. The remaining three ( $\bar{\Gamma}_{4,8}, \bar{\Gamma}_{5,7}, \bar{\Gamma}_{6,6}$ ) were identified numerically to approximately 22 significant digits and matched to rational numbers with relative differences below $4 \times 10^{-23}$ .
Diagrammatic Reduction:
The methodology demonstrates a massive reduction in computational complexity. For instance, the 74,092 Feynman diagrams contributing to the $G^{11}$ term collapse into just 18 pairings of field responses. This bypasses the need to solve for individual high-order field components, relying instead on integrated pairings.

Significance and Claims
The paper claims to bridge the gap between the infinitesimal probe limit and the regime of comparable masses. By providing the first exact closed-form expression for the first backreaction and the complete static potential up to $G^{12}$ , the work:

Exposes a "hidden analytical simplicity" within Einstein's theory, specifically the rational nature of the coefficients and the collective cancellations in the mode sums.
Provides essential building blocks for modeling orbital dynamics and gravitational-wave generation, particularly for the conservative sector of the interaction.
Offers a framework (stationary energy, field-response pairings, exact backgrounds) that the authors suggest can be extended to moving sources and velocity-dependent sectors, potentially aiding post-Newtonian and post-Minkowskian calculations.

The authors remain modest regarding future implications, noting that whether the collective cancellations observed in the static sector survive in time-dependent scenarios (moving sources) remains an open question. They also note that the rationality of coefficients at orders beyond $G^{12}$ is currently an open hypothesis.

More like this