We briefly discuss a method to construct a global, analytic, approximate spacetime for precessing, spinning binary black holes. The spacetime construction is broken into three parts: the inner zones are the spacetimes close to each black hole, and are approximated by perturbed Kerr solutions; the near zone is far from the two black holes, and described by the post-Newtonian metric; and finally the wave (far) zone, where retardation effects need to be taken into account, is well modeled by the post-Minkowskian metric. These individual spacetimes are then stitched together using asymptotic matching techniques to obtain a global solution that approximately satisfies the Einstein field equations. Precession effects are introduced by rotating the black hole spin direction according to the precessing equations of motion, in a way that is consistent with the global spacetime construction.

We construct a new global, fully analytic, approximate spacetime which accurately describes the dynamics of non-precessing, spinning black hole binaries during the inspiral phase of the relativistic merger process. This approximate solution of the vacuum Einstein's equations can be obtained by asymptotically matching perturbed Kerr solutions near the two black holes to a post-Newtonian metric valid far from the two black holes. This metric is then matched to a post-Minkowskian metric even farther out in the wave zone. The procedure of asymptotic matching is generalized to be valid on all spatial hypersurfaces, instead of a small group of initial hypersurfaces discussed in previous works. This metric is well suited for long term dynamical simulations of spinning black hole binary spacetimes prior to merger, such as studies of circumbinary gas accretion which requires hundreds of binary orbits.