Mathisson–Papapetrou–Dixon Equations

• The Mathisson–Papapetrou–Dixon equations describe the motion of extended test bodies with intrinsic spin in curved spacetime using a pole–dipole approximation and require a spin supplementary condition to define the centroid.
• They embody a robust framework that integrates covariant stress-energy conservation with relativity, enabling precise modeling of gravitational-wave generation and orbital dynamics in strong-field regimes.
• The formalism extends to include higher multipoles, various matter couplings, and alternative geometric structures, offering versatile analytical and numerical tools for astrophysical applications.

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 [(Costa et al., 2012); (Iosifidis, 18 Sep 2025); (Andersson et al., 8 Dec 2025); (Deriglazov et al., 2015)]. 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) (Han et al., 2016, Almonacid et al., 2017). Extensions to non-Riemannian geometry include explicit contributions from torsion and non-metricity, leading to additional coupling terms with the hypermomentum tensor (Iosifidis, 18 Sep 2025, Fabbri, 2024).

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 (Andersson et al., 8 Dec 2025, Witzany et al., 2018, Benavides-Gallego et al., 2023). 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$ (Velandia et al., 2016, Deriglazov et al., 2015, Benavides-Gallego et al., 2023).
• 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 (Costa et al., 2012).
• 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 (Witzany et al., 2018, Benavides-Gallego et al., 2023).

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” (Timogiannis et al., 2022, Benavides-Gallego et al., 2023). 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 (Witzany et al., 2018). 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) (Witzany et al., 2018).

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 (Deriglazov et al., 2015, Deriglazov et al., 2015). 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 (Costa et al., 2012).
• 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 (Costa et al., 2012).
• 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 (Velandia et al., 2016, Velandia-Heredia et al., 2017).
• 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 (Han et al., 2016).
• 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 [(Andersson et al., 8 Dec 2025); (Obukhov et al., 2014)].

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 (Andersson et al., 8 Dec 2025, Semerák, 2015).

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 [(Plyatsko et al., 2011); (Lukes-Gerakopoulos, 2017)]. 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 (Witzany et al., 2018).
• 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 (Andersson et al., 8 Dec 2025, Witzany et al., 2018, Velandia et al., 2016).

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) [(Fabbri, 2024); (Iosifidis, 18 Sep 2025); (Obukhov et al., 2014)]. 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 (Benavides-Gallego et al., 2023, Timogiannis et al., 2022).

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 (Semerák et al., 2015, Semerák, 2015). 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" (Costa et al., 2012)
• "Conserved quantities and integrability for massless spinning particles in general relativity" (Andersson et al., 8 Dec 2025)
• "Hamiltonians and canonical coordinates for spinning particles in curved space-time" (Witzany et al., 2018)
• "Lagrangian formulation for Mathisson-Papapetrou-Tulczyjew-Dixon (MPTD) equations" (Deriglazov et al., 2015)
• "Mathisson-Papapetrou-Tulczyjew-Dixon (MPTD) equations in ultra-relativistic regime and gravimagnetic moment" (Deriglazov et al., 2015)
• "Generalized Papapetrou's equations of motion for an extended test body within static and isotropic metrics" (Almonacid et al., 2017)
• "Spinning particles in vacuum spacetimes of different curvature types" (Semerák et al., 2015)
• "Spinning test body orbiting around a Kerr black hole: Comparing Spin Supplementary Conditions..." (Timogiannis et al., 2022)
• "Comparing spin supplementary conditions for particle motion around traversable wormholes" (Benavides-Gallego et al., 2023)
• "Numerical solution of Mathisson-Papapetrou-Dixon equations for spinning test particles in a Kerr metric" (Velandia-Heredia et al., 2017)
• Additional context throughout from (Velandia et al., 2016, Semerák, 2015, Plyatsko et al., 2011, Han et al., 2016, Loomis et al., 2017, Obukhov et al., 2014, Fabbri, 2024, Iosifidis, 18 Sep 2025, Lukes-Gerakopoulos, 2017).