Abstract:The present paper studies the post-Newtonian dynamics of N bodies in general relativity. We derive covariant equations of translational and rotational motion of N extended bodies having arbitrary distribution of mass and velocity of matter by employing the set of global and local coordinate charts on curved spacetime manifold of N-body system along with the mathematical apparatus of the Cartesian STF tensors and Blanchet-Damour multipole formalism. We separate the self-field effects of the bodies from the external gravitational environment and construct the effective background spacetime manifold by making use of the asymptotic matching technique. We make worldline of the center of mass of each body identical with that of the origin of the body-adapted local coordinates by the appropriate choice of the dipole moments. The covariant equations of motion are obtained by applying the Einstein principle of equivalence and the Fermi-Walker law of transportation of the linear momentum and spin of each body. Our approach significantly extends the Mathisson-Papapetrou-Dixon covariant equations of motion beyond the spin-dipole particle approximation by accounting for the entire infinite set of the internal multipoles of the bodies which are gravitationally coupled with the curvature tensor of the background manifold and its covariant derivatives. The results of our study can be used for much more accurate prediction of orbital dynamics of extended bodies in inspiraling binary systems and construction of templates of gravitational waves at the merger stage when the strong gravitational interaction between the higher-order multipoles of the bodies play a dominant role. The covariant theory of the post-Newtonian equations of motion beyond the spin-dipole approximation is a solid foundation for future improvements in long-term accuracy of relativistic celestial ephemerides of the solar system bodies.

Equivariant optimisation for the gravitational $n$-body problem: a computational factory of symmetric orbits

Symmetries and periodic orbits for the $n$-body problem: about the computational approach

Symplectic methods in the numerical search of orbits in real-life planetary systems

Looking for Stable Celestial Systems Using Bayesian Optimisation

Computational symplectic topology and symmetric orbits in the restricted three-body problem

Unchained polygons and the N-body problem

Symplectic geometry and space mission design

orbitN: A symplectic integrator for planetary systems dominated by a central mass -- Insight into long-term solar system chaos

Periodic and Quasi-Periodic Orbits in the Collinear Four-Body Problem: A Variational Analysis

Numerical Methods and Closed Orbits in the Kepler-Heisenberg Problem

Covariant Equations of Motion Beyond the Spin-Dipole Particle Approximation

A "barbell" in a central force field: A case study in symmetry reduction

On the Stability of the Three Classes of Newtonian Three-Body Planar Periodic Orbits

Dynamical Configurations of Celestial Systems Comprised of Multiple Irregular Bodies

An efficient approach for searching three-body periodic orbits passing through Eulerian configuration

Gravitational spin-orbit coupling through the third-subleading post-Newtonian order: Exploring spin-gauge flexibility

A commend on "Three Classes of Newtonian Three-Body Planar Periodic Orbits" by Šuvakov and Dmitrašinović (PRL, 2013)

Symmetric Periodic Orbits in Three Sub-Problems of the N-body Problem

Parallelization of the Symplectic Massive Body Algorithm (SyMBA) $N$-body Code

A new code for orbit analysis and Schwarzschild modelling of triaxial stellar systems