
kinematical evolution [13] of pa, Sab, it is natural to look for a Hamiltonian formulation1of
these equations. In the case of general relativity, two distinct approaches have been proposed
in the past. One of them is based on extending, to curved spacetime, the Lagrangian
formulation of a free spinning particle in flat spacetime [23]. Inspired by this work and
similar development in effective field theory [24], the authors of [25] provided a Hamiltonian
formulation of the dipolar MPTD system via a singular Legendre transformation of this
Lagrangian. Extensions to quadrupolar [26] and octupolar orders [27] have been derived,
allowing for the SSC to be understood as a gauge freedom [28]. One drawback of these
Lagrangian formulations, however, is their non-covariance, i.e., they require (i) a particular
3+1 split of the background spacetime and (ii) a fine-tuned choice of spin supplementary
condition with no clear physical interpretation.
The other approach can be traced back to Souriau’s symplectic approach [29] to the prob-
lem, cf. the recent note on Souriau’s pioneering ideas [30]. While our treatment does not
follow directly Souriau’s, it is solely based on symplectic geometry too, and borrow Souriau’s
general philosophy: the Hamiltonian system will be considered the prime ingredient, and
not a mere Legendre transform of some Lagrangian. The general symplectic structure that
we shall use is already present (but somewhat hidden) in Souriau’s [29] and Kunzle’s work
[31], but it seems that its first clear appearance is in2[32] (see also references [33]). Natu-
rally, our formulation will share some similarities with other covariant approaches proposed
in [32–35]. However, the works [32,34] only work at the level of the equations of motion for
simple orbits and do not exploit the Hamiltonian nature of the system. Similarly, the Hamil-
tonians proposed in [33,35] have a number of drawbacks that are unsatisfactory for practical
applications, in particular ad-hoc mass parameters too many phase space dimensions, both
issues related to unresolved constraints.Lastly, the necessary choice of spin supplementary
condition for the well-posedness of the MPTD system will lead us, inevitably, to algebraic
constraints. In the Lagrangian formalism, this is usually treated using the Dirac-Bergmann
algorithm [25,36]. Here, in the context of symplectic geometry, this will be handled within
the classical theory of “constrained Hamiltonian systems” following classical textbooks [37]
(or [38] for detailed proofs).
C. The question of integrability around black holes
At monopolar order, the spin tensor Sab is not involved in the description of the body.
The two MPTD equations reduce to (i) pais parallel-transported along the worldline and (ii)
pais tangent to the worldline. Combining the two implies that the worldline is a geodesic of
the background spacetime. The equations of geodesic motion are the most general setup to
describe “free systems”, as was first realized by Arnold in [39].Owing to their purely free (i.e.,
purely kinematical) evolution, such geodesic systems are usually simple enough that they
can sometimes be integrable, i.e., possess enough integrals of motion that make the solution
to their equations of motion solvable by quadrature (i.e., using algebraic manipulations
and standard calculus). In general relativity, the most remarkable example of this is the
integrability of geodesics in the Kerr spacetime, first demonstrated by Carter [40], five
years after Kerr’s discovery of the eponymous metric describing rotating black holes [41].
Integrability of Kerr geodesics is far from trivial: there are not enough spacetime symmetries
1By “Hamiltonian formulation”, we mean a system of ordinary differential equations generated by a scalar
function and some geometric structure on a phase space. We shall be more precise when necessary.
2However, in [32] it is claimed that these brackets are symplectic, although they are not.
3