From geodesics of the multipole solutions to the perturbed Kepler problem

A static and axisymmetric solution of the Einstein vacuum equations with a finite number of relativistic multipole moments (RMM) is written in multipole symmetry adapted (MSA) coordinates up to certain order of approximation, and the structure of its metric components is explicitly shown. From the equation of equatorial geodesics, we obtain the Binet equation for the orbits and it allows us to determine the gravitational potential that leads to the equivalent classical orbital equations of the perturbed Kepler problem. The relativistic corrections to Keplerian motion are provided by the different contributions of the RMM of the source starting from the monopole (Schwarzschild correction). In particular, the perihelion precession of the orbit is calculated in terms of the quadrupole and 2⁴-pole moments. Since the MSA coordinates generalize the Schwarzschild coordinates, the result obtained allows measurement of the relevance of the quadrupole moment in the first order correction to the perihelion frequency-shift.