An equation of motion for relativistic compact binaries is derived through
the third post-Newtonian (3 PN) approximation of general relativity. The strong
field point particle limit and multipole expansion of the stars are used to
solve iteratively the harmonically relaxed Einstein equations. We take into
account the Lorentz contraction on the multipole moments defined in our
previous works. We then derive a 3 PN acceleration of the binary orbital motion
of the two spherical compact stars based on a surface integral approach which
is a direct consequence of local energy momentum conservation. Our resulting
equation of motion admits a conserved energy (neglecting the 2.5 PN radiation
reaction effect), is Lorentz invariant and is unambiguous: there exist no
undetermined parameter reported in the previous works. We shall show that our 3
PN equation of motion agrees physically with the Blanchet and Faye 3 PN
equation of motion if λ=−1987/3080, where λ is the parameter
which is undetermined within their framework. This value of λ is
consistent with the result of Damour, Jaranowski, and Sch\"afer who first
completed a 3 PN iteration of the ADM Hamiltonian in the ADMTT gauge using the
dimensional regularization.Comment: 52 pages, no figure, Appendices B and D added. Phys. Rev. D in pres