Given a smooth globally hyperbolic (3+1)-dimensional spacetime satisfying the Einstein vacuum equations (possibly with cosmological constant) and an inextendible timelike geodesic, we construct a family of metrics depending on a small parameter \epsilon>0 with the following properties. (1) They solve the Einstein vacuum equations modulo \mathcal{O}(\epsilon^\infty). (2) Away from the geodesic they tend to the original metric as \epsilon\to 0. (3) Their \epsilon^{-1}-rescalings near every point of the geodesic tend to a fixed subextremal Kerr metric (assuming a condition on the mode stability at zero frequency which we verify in the very slowly rotating case). Our results apply on spacetimes which do not admit nontrivial Killing vector fields in a neighborhood of a point on the geodesic. They also apply in a neighborhood of the domain of outer communications of subextremal Kerr and Kerr-de Sitter spacetimes, in which case our metrics model extreme mass ratio inspirals if we choose the timelike geodesic to cross the event horizon. The metrics which we construct here depend on ϵ and the (rescaled) coordinates on the original spacetime in a log-smooth fashion. This in particular justifies the formal perturbation theoretic setup in work of Gralla-Wald on gravitational self-force in the case of small black holes.