Einstein–Infeld–Hoffmann equations explained

The Einstein–Infeld–Hoffmann equations of motion, jointly derived by Albert Einstein, Leopold Infeld and Banesh Hoffmann, are the differential equations describing the approximate dynamics of a system of point-like masses due to their mutual gravitational interactions, including general relativistic effects. It uses a first-order post-Newtonian expansion and thus is valid in the limit where the velocities of the bodies are small compared to the speed of light and where the gravitational fields affecting them are correspondingly weak.

Given a system of N bodies, labelled by indices A = 1, ..., N, the barycentric acceleration vector of body A is given by:

\begin{align} \vec{a} _{A} & = \sum _{B} (G m _{B} \vec{n)/(_{BA}) } \\& \quad + \frac \sum_ \frac \left[v_A^2+2v_B^2 - 4(\vec{v}_A \cdot \vec{v}_B) - \frac{3}{2} (\vec{n}_{AB} \cdot \vec{v}_B)^2 \right. \\ & {} \qquad {} \left. {} - 4 \sum_{C \not = A} \frac{G m_C}{r_{AC}} - \sum_{C \not = B} \frac{G m_C}{r_{BC}} + \frac{1}{2}((\vec{x}_B-\vec{x}_A) \cdot \vec{a}_B) \right] \\& \quad + \frac \sum_ \frac\left[\vec{n}_{AB}\cdot(4\vec{v}_A-3\vec{v}_B)\right](\vec_A-\vec_B) \\& \quad + \frac \sum_ + O (c^)\end

where:

\vec{x} _{A} is the barycentric position vector of body A\vec{v} _{A=d\vec{x}} _{A/dt} is the barycentric velocity vector of body A\vec{a} ^{2} _{A/dt} is the barycentric acceleration vector of body Ar _{AB} = |\vec{x} _{A-\vec{x}} _{B|} is the coordinate distance between bodies A and B\vec{n} _{AB} = (\vec{x} _{A-\vec{x}} _{B)/r} _{AB} is the unit vector pointing from body B to body Am _{A} is the mass of body A.c is the speed of lightG is the gravitational constant

and the big O notation is used to indicate that terms of order c^-4 or beyond have been omitted.

The coordinates used here are harmonic. The first term on the right hand side is the Newtonian gravitational acceleration at A; in the limit as c → ∞, one recovers Newton's law of motion.

The acceleration of a particular body depends on the accelerations of all the other bodies. Since the quantity on the left hand side also appears in the right hand side, this system of equations must be solved iteratively. In practice, using the Newtonian acceleration instead of the true acceleration provides sufficient accuracy.^[1]

Further reading

• Einstein . A. . Infeld . L. . Hoffmann . B. . 1938 . The Gravitational Equations and the Problem of Motion . . Second series . 39 . 1 . 65–100 . 1968714 . 1938AnMat..39...65E . 10.2307/1968714 .
• Book: Kovalevsky, Jean . Seidelmann . P. Kenneth . 2004 . Fundamentals of Astrometry . limited . New York . . 173 . 0521642167 .
• Book: Landau, Lev . Lifshitz . Evgeny . 1971 . The classical theory of fields . Oxford . . 337 .

Notes and References

1. Standish, Williams. Orbital Ephemerides of the Sun, Moon, and Planets, Pg 4. Web site: Archived copy . 2010-04-03 . dead . https://web.archive.org/web/20110205225629/http://iau-comm4.jpl.nasa.gov/XSChap8.pdf . 2011-02-05 .