  Accurate calculation of the motion of a compact object in a background
spacetime induced by a supermassive black hole is required for the future
detection of such binary systems by the gravitational-wave detector LISA.
Reaching the desired accuracy requires calculation of the second-order
gravitational perturbations produced by the compact object. At the point
particle limit the second-order gravitational perturbation equations turn out
to have highly singular source terms, for which the standard retarded solutions
diverge. Here we study a simplified scalar toy-model in which a point particle
induces a nonlinear scalar field in a given curved spacetime. The corresponding
second-order scalar perturbation equation in this model is found to have a
similar singular source term, and therefore its standard retarded solutions
diverge. We develop a regularization method for constructing well-defined
causal solutions for this equation. Notably these solutions differ from the
standard retarded solutions, which are ill-defined in this case.
