  We summarize a recent work done on the title's subject. First, we present the
asymptotic scheme of post-Newtonian (PN) approximation for general relativity
in the harmonic gauge. Then, we discuss the definition of the mass centers and
the derivation of equations for their motion, following that scheme. Finally,
we briefly analyze the reason why a new term has thus been found in the
equations of motion.
