  Accurate calculation of the gradual inspiral motion in an extreme mass-ratio
binary system, in which a compact-object inspirals towards a supermassive
black-hole requires calculation of the interaction between the compact-object
and the gravitational perturbations that it induces. These metric perturbations
satisfy linear partial differential equations on a curved background spacetime
induced by the supermassive black-hole. At the point particle limit the
second-order perturbations equations have source terms that diverge as
$r^{-4}$, where $r$ is the distance from the particle. This singular behavior
renders the standard retarded solutions of these equations ill-defined. Here we
resolve this problem and construct well-defined and physically meaningful
solutions to these equations. We recently presented an outline of this
resolution [E. Rosenthal, Phys. Rev. D 72, 121503 (2005)]. Here we provide the
full details of this analysis. These second-order solutions are important for
practical calculations: the planned gravitational-wave detector LISA requires
preparation of waveform templates for the expected gravitational-waves.
Construction of templates with desired accuracy for extreme mass-ratio binaries
requires accurate calculation of the inspiral motion including the interaction
with the second-order gravitational perturbations.
