The value of that we found, , is perfectly consistent with the relation (1) and the result of  ().
Finally, let us discuss the ambiguity in the 3 PN equations of motion previously derived by Blanchet and Faye in [25, 27]. In their formalism, Dirac delta distributions are used to achieve the point particle limit. The (Lorentz invariant generalized) Hadamard partie finie regularization has been extensively employed to regularize divergences caused by their use of a singular source. In fact, unless regularization is employed, divergences occur both in the evaluation of the 3 PN field where (Poisson) integrals diverge at the location of the stars and in their derivation of the equations of motion where a substitution of the metric into a geodesic equation causes divergences.
When we regularize some terms at a point, say, , where the terms are singular, using the Hadamard partie finie regularization, roughly speaking we take an angular average of the finite part of the terms in the neighborhood of the singular point. Then if there are logarithmic terms such as , we should take an angular average over some sphere centered on with a finite radius. The radius of the sphere is arbitrary but we do not ignore it because we should ensure the argument of the logarithms to be dimensionless.
The problem here is that there is a priori no reason to expect that the radius for each star introduced to regularize the field and another radius for that star introduced to regularize the geodesic equation coincide with each other. Thus the Blanchet and Faye 3 PN equations of motion have four arbitrary constants instead of two in our equations of motion. In our framework, we can see the origin of the number if we assume that we have defined a different body zone in the derivation of the equations of motion from used in the derivation of the 3 PN field. However, in reality, we have only one body zone for each star. In our formalism the field is expressed in terms of the four-momentum (and multipole moments) which are defined as volume integrals over the body zone. On the other hand our general form of the equations of motion has been derived based on the conservation law of the four-momentum, and thus we evaluate the surface integrals in the general form of the equations of motion over the boundary of the body zone.
In fact, [25, 27] have shown that two of the four arbitrary constants can be removed by using a gauge freedom remaining in the harmonic gauge condition; the two places where the singular points exist are in some sense ambiguous. The remaining two turn out to appear as the ratios () where and are the four regularization parameters (roughly speaking, the radii of and in the terminology in the previous paragraph). Blanchet and Faye [25, 27] then proved that assuming the equations of motion are polynomials of the two masses of the stars, those two ratios should satisfy where and are pure numbers. Then they showed that in order for their equations of motion to admit conserved energy, then , while no argument was found to fix .
The above argument in turn means that the Blanchet and Faye 3 PN equations of motion do not give a conserved energy unless is different from . Damour, Jaranowski, and Schäfer  pointed out that there is an unsatisfactory feature in the generalized partie finie regularization which contradicts with the mathematical structure of general relativity. Indeed, by using dimensional regularization which is pointed out by them to be more satisfactory in this regard,  derived an unambigous ADM Hamiltonian in the ADM transeverse traceless gauge. Later, Blanchet et al.  used dimensional regularization and found that their new equations of motion physically agree with ours and admit a conserved energy.
This work is licensed under a Creative Commons Attribution-Noncommercial-No Derivative Works 2.0 Germany License.