Higher-order geodesic deviations applied to the Kerr metric

Starting with an exact and simple geodesic, we generate approximate geodesics by summing up higher-order geodesic deviations within a general relativistic setting, without using Newtonian and post-Newtonian approximations. We apply this method to the problem of closed orbital motion of test particles in the Kerr metric spacetime. With a simple circular orbit in the equatorial plane taken as the initial geodesic, we obtain finite eccentricity orbits in the form of Taylor series with the eccentricity playing the role of a small parameter. The explicit expressions of these higher-order geodesic deviations are derived using successive systems of linear equations with constant coefficients, whose solutions are of harmonic oscillator type. This scheme gives best results when applied to orbits with low eccentricities, but with arbitrary possible values of (GM/Rc²).