---
title: Mathisson–Papapetrou–Dixon Equations
url: https://www.emergentmind.com/topics/mathisson-papapetrou-dixon-equations
type: topic
---

# Mathisson–Papapetrou–Dixon Equations

The Mathisson–Papapetrou–Dixon (MPD) equations govern the motion of extended test bodies with intrinsic spin (pole–dipole approximation) in curved spacetime and are central to relativistic astrophysics, gravitational-wave modeling, and classical/quantum field correspondence. They comprise a closed dynamical system for the evolution of the worldline, linear momentum, and antisymmetric spin tensor of a compact object, truncated at dipole order. The system is intrinsically underdetermined and requires a spin supplementary condition (SSC) to select the physical centroid, with the choice of SSC encoding a gauge-like freedom in describing the center of mass. The MPD formalism admits diverse generalizations—different matter couplings (tensors, scalar fields), inclusion of higher multipoles, extensions to non-Riemannian backgrounds with torsion and non-metricity, and Hamiltonian or Lagrangian structures crucial for numerical integration in post-Newtonian and wave-generation applications.

## 1. Tensorial Structure and Physical Content

The standard MPD system arises from integrating the covariant conservation of the stress-energy tensor truncated at pole–dipole order. The essential dynamical variables are the worldline $z^\mu(\tau)$, normalized tangent $U^\mu = dz^\mu/d\tau$, linear momentum $P^\mu$, and antisymmetric spin tensor $S^{\mu\nu}=-S^{\nu\mu}$. The coupled equations read
\[
\frac{D P^\mu}{D\tau} = -\frac{1}{2}R^\mu{}_{\nu\alpha\beta}U^\nu S^{\alpha\beta}\,, \qquad 
\frac{D S^{\alpha\beta}}{D\tau} = 2P^{[\alpha}U^{\beta]} \,,
\]
where $\frac{D}{D\tau}$ is the covariant derivative along $z^\mu$, and $R^\mu{}_{\nu\alpha\beta}$ is the Riemann tensor of the background spacetime [1206.7093; 2509.14757; 2512.07677; 1509.04926]. This system contains more unknowns than equations (12 for $P^\mu$, $S^{\mu\nu}$, but only 8 equations), reflecting the necessity of physically defining the representative worldline within the extended body.

The pole–dipole formalism neglects quadrupole and higher corrections, which can be systematically included for more accurate modeling of deformable bodies (e.g., via $J^{\alpha\beta\gamma\delta}$ for quadrupoles) [1611.07602; 1703.07874]. Extensions to non-Riemannian geometry include explicit contributions from torsion and non-metricity, leading to additional coupling terms with the hypermomentum tensor [2509.14757; 2411.09651].

## 2. Spin Supplementary Conditions: Gauge Structure and Centroids

To close the MPD system, a spin supplementary condition (SSC) of the form $S^{\mu\nu}V_\nu=0$ is imposed, where $V^\mu$ is a timelike vector encoding the choice of observer used to define the body's centroid [2512.07677; 1808.06582; 2306.17394]. The principal SSCs include:

- **Tulczyjew–Dixon (TD):** $V^\mu=P^\mu/m$, $S^{\mu\nu}P_\nu=0$. This selects the centroid such that the mass-dipole moment vanishes in the momentum frame; the solution for the worldline is unique; the hidden momentum $P^\mu$ is generally not parallel to $U^\mu$ [1609.03630; 1509.04926; 2306.17394].

- **Mathisson–Pirani (MP):** $V^\mu=U^\mu$, $S^{\mu\nu}U_\nu=0$. This choice aligns the spin to be purely spatial in the rest frame of the reference worldline. It admits a gauge-like freedom, yielding a family of helical solutions (zitterbewegung) that encode different centroids for the same physical body [1206.7093].

- **Ohashi–Kyrian–Semerák (OKS):** $V^\mu=w^\mu$ is a parallel-transported timelike unit vector. Here $S^{\mu\nu}w_\nu=0$ enforces a fixed inertial centroid; the spin is Fermi-transported and hidden momentum vanishes [1808.06582; 2306.17394].

Each SSC corresponds to a different slicing of the stress-energy multipole expansion and is associated with a different but physically equivalent (at the dipole level) centroid within the Møller disk of possible worldlines. The freedom in choosing SSC is a benign gauge ambiguity, resolvable by appropriate “centroid shifts” [2206.11149; 2306.17394]. Differences between SSCs enter at higher order in the spin expansion for observables such as orbit frequencies and ISCO properties but can be bridged analytically by explicit shifts in coordinates and spin measures.

## 3. Hamiltonian and Lagrangian Structure

Hamiltonian formulations for the MPD system are vital for structure-preserving numerical integration, construction of effective-one-body models, and involvement in post-Newtonian expansions [1808.06582]. For each SSC, there exists a corresponding constrained Hamiltonian:

- **TD-Hamiltonian:** $H_{\text{TD}}$ involves a spin-dependent deformation of the metric and generates the equations on the $S^{\mu\nu}P_\nu=0$ surface.

- **MP-Hamiltonian:** $H_{\text{MP}}$ involves projecting out the spin directions orthogonal to $U^\mu$.

- **KS-Hamiltonian:** For the inertial (KS) condition, spin is a constant of motion and the dynamics simplifies.

Canonical spin coordinates can be constructed via tetrad decomposition of $S^{\mu\nu}$, allowing the use of high-order symplectic integrators—critical for long-term stability in modeling extreme mass ratio inspirals (EMRIs) [1808.06582].

The Lagrangian formalism for spinning bodies without auxiliary variables reproduces the MPD system, with explicit constraints automatically yielding an effective metric $G_{\mu\nu}(S)$, inherited from the spin structure [1509.04926; 1509.05357]. In the ultra-relativistic regime, the choice of effective metric becomes crucial for physical consistency.

## 4. Physical Effects and Observables

The inclusion of spin in the motion of test bodies induces rich physics:

- **Helical Motions and Zitterbewegung:** In flat spacetime under the Mathisson–Pirani SSC, the presence of helical solutions with frequency $\omega = M/S_\star$ precisely matching the Dirac zitterbewegung frequency for the electron ($2M_e/\hbar$) indicates a classical-quantum correspondence [1206.7093].

- **Hidden Momentum:** For SSCs with $P^\mu \not\parallel U^\mu$, the total momentum splits into kinetic and hidden pieces. Hidden momentum exchanges between spin and orbital degrees of freedom ensure total $P^\mu$ conservation—even when the CM follows a non-geodesic trajectory [1206.7093].

- **Spin-Orbit and Spin-Curvature Couplings:** In BH spacetimes (Kerr, Schwarzschild), the MPD equations—when truncated at the pole-dipole level and supplemented by the TD or MP condition—yield observable consequences such as the Lense–Thirring precession, gravitomagnetic clock effect (time difference up to $\sim 10^{-8}$ s due to spin coupling), and spin-corrected ISCO positions [1609.03630; 1712.01404].

- **Higher Multipole Extensions:** Quadrupole corrections contribute to the variation of rotational velocity and can induce chaos for nonphysical spin magnitudes, but for realistic spins, EMRI orbits remain regular [1611.07602].

- **Conserved Quantities:** In spacetimes with symmetries (Killing vectors, Killing–Yano tensors), the MPD equations admit linear and quadratic first integrals—generalized Carter constants and orbital angular momentum—regardless of SSC for leading-order spin [2512.07677; 1404.6977].

The generalization to massless particles (photon and graviton wave packets) under adapted SSCs reproduces the gravitational spin Hall effect and provides a complete set of first integrals in type D backgrounds [2512.07677; 1512.06253].

## 5. Applications and Numerical Integration

For numerical integration, the MPD system, together with appropriate SSC, is recast as a closed set of first-order ODEs. Reparameterization invariance under proper time or other affine parameters is preserved, but physical invariants (mass, spin, normalization) depend on both SSC and time parameterization [1110.1967; 1709.08942]. Efficient schemes involve:

- Dimensionless variables for coordinate, velocity, and spin components.
- Use of conserved quantities (energy, angular momentum) as algebraic constraints.
- High-order symplectic integrators, exploiting canonical Poisson bracket structure in phase space [1808.06582].
- Symmetrization and constraint-preserving methods to ensure fidelity to normalization and SSC at each step.

Astrophysical applications include precise modeling of gravitational waves for space-based detectors (LISA), probe of deviations from geodesicity in measurements of ISCOs and precession, and, for wave packets, incorporation of spin Hall deviations in photon and graviton propagation [2512.07677; 1808.06582; 1609.03630].

## 6. Generalizations and Theoretical Impact

The MPD formalism generalizes to scalar–tensor theories (modifications in the momentum and spin evolution due to direct scalar coupling), metric–affine backgrounds (nonzero torsion and non-metricity), and coupling to classical or quantum spinor fields (where the MPD system follows from field-theoretic conservation laws) [2411.09651; 2509.14757; 1404.6977]. The latter provides insight into the hydrodynamic and multipolar reconstruction of quantum field behavior at macroscopic scales.

In higher-order expansions, care is required as the pole–dipole truncation ceases to be accurate for large spins or objects with significant internal structure. The observed convergence across SSCs for physical spins, breakdown for large spins or orbits close to ISCO, and restoration by inclusion of centroid and spin corrections corroborate the gauge-like nature of the SSC ambiguity in physical predictions [2306.17394; 2206.11149].

## 7. Mathematical Structure and Geometric Interpretation

The geometric structure of the MPD equations illuminates the deep links between internal body structure, spacetime curvature, and the observer dependence of multipole expansions. The algebraic type (Petrov classification) of the background determines which curvature invariants (Weyl scalars $\Psi_2$, etc.) enter the force, and alignment of spin-eigenplanes with principal null directions can dramatically simplify the form of the evolution equations [1505.01069; 1512.06253]. In the massless limit, the reduction to a minimal set of Weyl scalars ($\Psi_1,\Psi_2$) further showcases the interplay between null propagation and spin-orbit interactions.

---

**References:**

- "Mathisson's helical motions demystified" [1206.7093]
- "Conserved quantities and integrability for massless spinning particles in general relativity" [2512.07677]
- "Hamiltonians and canonical coordinates for spinning particles in curved space-time" [1808.06582]
- "Lagrangian formulation for Mathisson-Papapetrou-Tulczyjew-Dixon (MPTD) equations" [1509.04926]
- "Mathisson-Papapetrou-Tulczyjew-Dixon (MPTD) equations in ultra-relativistic regime and gravimagnetic moment" [1509.05357]
- "Generalized Papapetrou's equations of motion for an extended test body within static and isotropic metrics" [1703.07874]
- "Spinning particles in vacuum spacetimes of different curvature types" [1505.01069]
- "Spinning test body orbiting around a Kerr black hole: Comparing Spin Supplementary Conditions..." [2206.11149]
- "Comparing spin supplementary conditions for particle motion around traversable wormholes" [2306.17394]
- "Numerical solution of Mathisson-Papapetrou-Dixon equations for spinning test particles in a Kerr metric" [1712.01404]
- Additional context throughout from [1609.03630], [1512.06253], [1110.1967], [1611.07602], [1701.01545], [1404.6977], [2411.09651], [2509.14757], [1709.08942].

Source: https://www.emergentmind.com/topics/mathisson-papapetrou-dixon-equations