Einstein–Infeld–Hoffmann equations

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 \not = A} \frac{G m_B \vec{n}_{BA}}{r_{AB}^2} \\ & {} \quad{} + \frac{1}{c^2} \sum_{B \not = A}   \frac{G m_B \vec{n}_{BA}}{r_{AB}^2} \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{1}{c^2} \sum_{B \not = A} \frac{G m_B}{r_{AB}^2}\left[\vec{n}_{AB}\cdot(4\vec{v}_A-3\vec{v}_B)\right](\vec{v}_A-\vec{v}_B) \\ & {} \quad {} + \frac{7}{2c^2} \sum_{B \not = A}{ \frac{G m_B \vec{a}_B }{r_{AB}}} + O (c^{-4}) \end{align} $$

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}_A=d^2\vec{x}_A/dt^2$$ is the barycentric acceleration vector of body A
 * $$r_{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 A
 * $$m_A$$ is the mass of body A.
 * $$c$$ is the speed of light
 * $$G$$ is the gravitational constant
 * and the big O notation is used to indicate that terms of order c&minus;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.