Construction of N-body time-symmetric initial data sets in general relativity

Given a collection of N asymptotically Euclidean ends with zero scalar curvature, we construct a Riemannian manifold with zero scalar curvature and one asymptotically Euclidean end, whose boundary has a neighborhood isometric to the disjoint union of a specified collection of sub-regions of the given ends. An application is the construction of time-symmetric solutions of the constraint equations which model N-body initial data.