On the equation of motion of compact binaries in Post-Newtonian   approximation
  A third post-Newtonian (3 PN) equation of motion for two spherical compact stars in a harmonic coordinate has been derived based on the surface integral approach and the strong field point particle limit. The strong field point particle limit enables us to incorporate a notion of a self-gravitating regular star into general relativity. The resulting 3 PN equation of motion is Lorentz invariant, unambiguous, and conserves an energy of the binary orbital motion. 
Yousuke Itoh