Now we impose the harmonic coordinate condition on the metric,. The explicit form of in the harmonic gauge is
Now we rewrite the relaxed Einstein equations into an integral form,)
We solve the Einstein equations as follows. First we split the integral region into two zones: the near zone and the far zone.
The near zone is the region containing the gravitational wave source where the wave character of the gravitational radiation is not manifest. In other words, in the near zone the retardation effect on the field is negligible. The near zone covers the whole source system. The size of the near zone is about a little larger than one wave length of the gravitational wave emitted by the source. In this paper we take the near zone as a sphere centered at some fixed point and enclosing the binary system. The radius of this sphere is set to be , where is arbitrary but larger than the size of the binary and the wave length of the gravitational radiation. The scaling of the near zone radius is derived from the dependence of the wavelength of the gravitational radiation emitted due to the orbital motion of the binary. Note that roughly speaking the frequency of such a wave is about twice the Keplerian frequency of the binary. The center of the near zone sphere would be determined, if necessary, for example, to be the center of mass of the near zone. The outside of the near zone is the far zone where the retardation effect of the field is crucial.
For the near zone field point , we write the field as
The far zone contribution can be evaluated with the DIRE method developed in . An explicit calculation shows that apparently there are the far zone contributions to the near zone field at 3 PN order. However, these 3 PN contributions are merely a gauge. Pati and Will showed that the far zone contribution does not affect the equations of motion up to 3 PN order inclusively and that the far zone contribution first appears at 4 PN order. This result is consistent with the earlier result of Blanchet and Damour  who used the multipolar-post-Minkowskian formalism. We follow the DIRE method and checke that the far zone contribution does not affect the equations of motion up to 3 PN order in Appendix A. Henceforth we shall focus our attention on the near zone contribution and do not write down the far zone contribution in the following calculation of the field.
As for the homogeneous solution, we shall ignore it for simplicity. If we take random initial data for the field  supposed to be of 1 PN order , they are irrelevant to the dynamics of the binary system up to the radiation reaction order . As we have assumed in the previous Section 3.3 that the magnitude of the free data of the gravitational field on the initial hypersurface is 2.5 PN order, we expect that the homogeneous solution does not affect the equations of motion up to 3 PN order. We leave a full implementation of the initial value formulation on the field as future work.
It is worth noticing that when we let become sufficiently large, then the condition corresponds to a no-incoming radiation condition at (Minkowskian) past null infinity (see e.g. ). Equation (69) can be written down as Kirchhoff’s formula,
Now we shall devote ourselves to the evaluation of the near zone contribution to the near zone field,
This work is licensed under a Creative Commons Attribution-Noncommercial-No Derivative Works 2.0 Germany License.