⚛️ general relativity

Geodesic deviation to all orders via a tangent bundle formalism

This paper introduces a tangent bundle formalism using covariant differential forms to derive an exact, all-orders in-in Lagrangian and equations for finite geodesic deviation, providing explicit expressions up to the tenth order that confirm previous fourth-order results.

Original authors: Joon-Hwi Kim

Published 2026-07-02

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Joon-Hwi Kim

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

The Big Picture: Measuring the "Stretch" of Space

Imagine you are floating in space with a friend. You are both falling freely toward a massive object, like a black hole. In a perfect, flat universe, you would stay the same distance apart forever. But in our universe, space is curved.

As you fall, the curvature of space acts like a giant, invisible hand. It might pull your friend closer to you, or push them further away, or even twist your path. This change in distance between two falling objects is called geodesic deviation.

For over 90 years, physicists have had a formula to describe this, but it only works for tiny distances. It's like a ruler that is accurate for measuring a grain of sand, but if you try to measure a mountain with it, the math breaks down.

This paper introduces a new, super-powerful ruler. The author, Joon-Hwi Kim, has figured out how to calculate exactly how space stretches and twists between two objects, no matter how far apart they are. He didn't just stop at a few steps; he found the formula for infinite steps (all orders), providing a complete map of how space behaves.

The Old Way vs. The New Way

The Old Way (Synge Calculus):
Previously, to solve this problem, physicists used a method called "Synge calculus." Imagine trying to navigate a maze by looking at a map of the entire maze at once, trying to connect the start point to the finish point simultaneously. It's incredibly complex, involving a web of intricate rules and "bitensors" (mathematical objects that live at two different places at once). It's like trying to solve a puzzle where the pieces keep changing shape as you move them.

The New Way (Tangent Bundle Formalism):
Kim proposes a different approach. Instead of looking at the start and end points simultaneously, he treats the problem like a movie.

The Setup: Imagine a "Tangent Bundle" not as a scary math term, but as a control room for a spaceship. In this control room, you have two dials: one for where you are (position) and one for how fast you are moving (velocity).
The Engine: Kim defines a special "engine" (a vector field) in this control room. When you turn this engine on, it doesn't just move the ship; it simulates the flow of time and the curvature of space all at once.
The Magic Trick: He uses a mathematical tool called a "Lie derivative" (think of it as a shapeshifting brush) to paint the future path of the ship. By "brushing" the current state of the ship with this tool, he can instantly generate the exact position and velocity of the ship at any future moment, accounting for all the twists and turns of space.

The "Chemical" Recipe

One of the most creative parts of the paper is how the author organizes the results. The formulas for these infinite corrections are incredibly long and complex. To make sense of them, Kim uses an analogy from organic chemistry.

The Atoms: The building blocks of his formulas are the Riemann tensors. In our analogy, these are Carbon atoms. They represent the curvature of space.
The Chains: When you combine these "Carbon atoms" with the distance between the two falling objects, they form long chains.
The Molecules: The final equations look like complex molecules. Kim even provides a visual "chemical notation" where the math looks like a diagram of a molecule.
- A simple chain represents a small correction.
- A branching molecule represents a more complex interaction where the curvature twists in multiple directions.

The paper provides the "recipe" for these molecules up to a very high level of complexity (the 10th order), which is like having the blueprint for a skyscraper instead of just a house.

Why This Matters (According to the Paper)

The paper claims that this new method is better than the old one because:

It's Exact: It works for any distance, not just tiny ones.
It's Systematic: The "chemical" structure makes it easy to see the pattern and calculate higher levels of accuracy without getting lost in the math.
It's Verified: The author checked his results against previous work (which only went up to the 4th or 5th order) and found they matched perfectly. He then extended the calculation all the way to the 10th order.

A Surprising Connection

The paper ends with a fascinating twist. The author notes that the math he developed for falling particles in space looks very similar to the math needed to describe a spinning black hole.

He suggests that his "all-orders" formula for falling particles might be the key to unlocking the equations for how a spinning black hole moves. It's as if he was trying to fix a leaky faucet, but in doing so, he accidentally discovered the blueprint for a jet engine. He promises to explore this connection in a future paper.

Summary

In short, this paper replaces a confusing, static map of space with a dynamic, flowing engine. It translates complex gravitational physics into a "chemical" language of carbon chains, allowing scientists to calculate exactly how space stretches and twists between two objects, no matter how far apart they are.

Technical Summary: Geodesic Deviation to All Orders via a Tangent Bundle Formalism

Problem Statement
The geodesic deviation equation (GDE) is a cornerstone of general relativity, describing the evolution of the infinitesimal separation between two nearby free-falling test particles. While the standard linearized GDE is well-established, generalizing this equation to finite separations has remained a challenging problem. Previous efforts have successfully derived corrections up to fourth order in the separation vector $y^\mu$ and provided all-orders formulae in terms of complex tensor expressions involving Synge bitensors. However, obtaining explicit, fully expanded equations or Lagrangians at arbitrary orders requires solving intricate recursion relations derived from the Synge bitensor formalism, which relies on "in-out" boundary conditions (characterizing a geodesic by its two endpoints). The objective of this work is to derive the exact, all-orders GDE and its associated Lagrangian for finite geodesic deviation in a manifestly covariant form, expressed explicitly in terms of the Riemann tensor and its derivatives.

Methodology
The authors propose an alternative framework to the traditional Synge calculus by formulating geodesic deviation as an initial-value problem within the tangent bundle $TM$ of spacetime.

Tangent Bundle Formalism: The authors treat the tangent bundle $TM$ (with local coordinates $(x^\mu, y^\mu)$ ) as a $2d$ -dimensional manifold. They define a specific vector field $N$ on $TM$ that acts as the generator of geodesic deviation. This vector field is invariant under coordinate transformations and encodes the first-order formulation of the geodesic equation.
In-In Boundary Conditions: Unlike the Synge formalism which uses endpoints $x$ and $z$ (in-out), this approach specifies an initial point $x$ and a tangent vector $y$ at $x$ (in-in). The finite deviation is generated by the unit-time flow of the vector field $N$ , denoted as $e^N$ .
Covariant Lie Derivative and Dressing: To handle the calculations efficiently and maintain covariance at the original point $x$ , the authors introduce a covariant Lie derivative, $\mathcal{L}^D_N = D\iota_N + \iota_N D$ , where $D$ is the covariant exterior derivative. They utilize a "dressing identity," which states that the exponential of this operator, $e^{\mathcal{L}^D_N}$ , computes the value of a tensor at the deviated point and parallel-transports it back to the original point via Wilson lines (parallel propagators).
Recursive Structure: By applying the Cartan magic formula and the properties of the interior product $\iota_N$ , the authors derive a recursive structure for the expansion of the parallel propagator and the deviation vector. This leads to the identification of "Q-tensors," which are products of the Riemann tensor and its covariant derivatives contracted with the separation vector $y$ .

Key Contributions and Results

Exact All-Orders Formulae: The paper derives explicit, all-orders formulae for the Jacobi propagators (the tensors relating the deviation vector and its derivative at the observer's location to those at the test particle's location). These are expressed as infinite sums where the coefficients are products of binomial coefficients. The terms involve ordered partitions of integers, and the number of terms at each order relates to Fibonacci numbers.
Explicit Expressions up to Tenth Order: The authors provide the explicit expansion of the Jacobi propagators, the Lagrangian, and the geodesic deviation equation up to tenth order ( $O(y^{10})$ ). These results are provided in ancillary files (All.nb, Low.nb, Q.nb) and summarized in the appendices.
Lagrangian Formulation: A second-order Lagrangian for the all-orders GDE is derived. This Lagrangian describes the dynamics of the deviation vector $y^\mu$ relative to an observer's worldline. The authors verify that their second-order Lagrangian reproduces the known results up to $O(y^5)$ found by Vines (2014).
The All-Orders GDE: By varying the derived Lagrangian or directly manipulating the first-order equations of motion using the dressing identity, the authors obtain the explicit second-order geodesic deviation equation. The equation includes nonlinear tidal effects due to the finite size of the separation, coupling the Riemann tensor and its derivatives to powers of the separation vector $y$ and the relative velocity.
Zero-Torsion Identities: The formalism naturally yields a set of identities (Eqs. 33a–33c) relating the derived tensors ( $X, Y, \dot{X}, \dot{Y}$ ) to the zero-torsion condition of the Levi-Civita connection, providing consistency checks for the expansions.

Significance and Claims
The paper claims to establish a robust "in-in" formalism for geodesic deviation that serves as a computationally efficient alternative to the Synge bitensor calculus. By shifting the framework to the tangent bundle and utilizing covariant differential forms, the authors bypass the need for solving the complex recursion relations inherent in the Synge formalism.

The primary significance lies in the provision of explicit, fully expanded expressions for the GDE and its Lagrangian to high orders ( $O(y^{10})$ ), which were previously inaccessible in closed form. The authors note that their results agree with previous fourth-order calculations and Vines' $O(y^5)$ results, validating the new approach.

While the paper focuses on the mathematical derivation, it suggests potential applications in:

Gravitational Physics: Improving the accuracy of gravitational wave data analysis (referencing LIGO) and the analytical treatment of eccentric relativistic orbits.
Theoretical Extensions: The framework is described as versatile, with potential applications in manifestly covariant perturbation theories for gauge and gravitational interactions, worldline formalisms, and field space geometry.
Black Hole Physics: The authors note a surprising connection to the problem of finding all-orders-in-spin equations of motion for Kerr black holes, suggesting that the all-orders GDE for an "imaginary deviation" may deduce parts of the nonlinear completion of the Mathisson-Papapetrou-Dixon equations. This specific connection is noted as a subject for a follow-up article.

The paper concludes by presenting a novel "chemical" notation (visualizing tensor structures as carbon chains) to manage the complexity of the high-order terms, offering an intuitive way to represent the intricate algebraic structures of the GDE.