We explore how a recently developed analytical black-hole binary spacetime can be extended using numerical simulations to go beyond the slow-inspiral phase. The analytic spacetime solves the Einstein field equations approximately, with the approximation error becoming progressively smaller the more separated the binary. To continue the spacetime beyond the slow-inspiral phase, we need to transition. Such a transition was previously explored at smaller separations. Here, we perform this transition at a separation of D=20M (large enough that the analytical metric is expected to be accurate), and evolve for six orbits. We find that small constraint violations can have large dynamical effects, but these can be removed by using a constraint-damping system like the conformal covariant formulation of the Z4 system. We find agreement between the subsequent numerical spacetime and the predictions of post-Newtonian theory for the waveform and inspiral rate that is within the post-Newtonian truncation error.

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.

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave astrophysics communities. The purpose of NINJA is to study the ability to detect gravitational waves emitted from merging binary black holes and recover their parameters with next-generation gravitational-wave observatories. We report here on the results of the second NINJA project, NINJA-2, which employs 60 complete binary black hole hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. In a "blind injection challenge" similar to that conducted in recent LIGO and Virgo science runs, we added 7 hybrid waveforms to two months of data recolored to predictions of Advanced LIGO and Advanced Virgo sensitivity curves during their first observing runs. The resulting data was analyzed by gravitational-wave detection algorithms and 6 of the waveforms were recovered with false alarm rates smaller than 1 in a thousand years. Parameter estimation algorithms were run on each of these waveforms to explore the ability to constrain the masses, component angular momenta and sky position of these waveforms. We also perform a large-scale monte-carlo study to assess the ability to recover each of the 60 hybrid waveforms with early Advanced LIGO and Advanced Virgo sensitivity curves. Our results predict that early Advanced LIGO and Advanced Virgo will have a volume-weighted average sensitive distance of 300Mpc (1Gpc) for $10M_{\odot}+10M_{\odot}$ ($50M_{\odot}+50M_{\odot}$) binary black hole coalescences. We demonstrate that neglecting the component angular momenta in the waveform models used in matched-filtering will result in a reduction in sensitivity for systems with large component angular momenta. [Abstract abridged for ArXiv, full version in PDF]

We study conformally-flat initial data for an arbitrary number of spinning black holes with exact analytic solutions to the momentum constraints constructed from a linear combination of the classical Bowen-York and conformal Kerr extrinsic curvatures. The solution leading to the largest intrinsic spin, relative to the ADM mass of the spacetime epsilon_S=S/M^2_{ADM}, is a superposition with relative weights of Lambda=0.783 for conformal Kerr and (1-Lambda)=0.217 for Bowen-York. In addition, we measure the spin relative to the initial horizon mass M_{H_0}, and find that the quantity chi=S/M_{H_0}^2 reaches a maximum of \chi^{max}=0.9856 for Lambda=0.753. After equilibration, the final black-hole spin should lie in the interval 0.9324<chi_{final}<0.9856. We perform full numerical evolutions to compute the energy radiated and the final horizon mass and spin. We find that the black hole settles to a final spin of chi_{final}^{max}=0.935 when Lambda=0.783. We also study the evolution of the apparent horizon structure of this "maximal" black hole in detail.

As 2 black holes bound to each other in a close binary approach merger their inspiral time becomes shorter than the characteristic inflow time of surrounding orbiting matter. Using an innovative technique in which we represent the changing spacetime in the region occupied by the orbiting matter with a 2.5PN approximation and the binary orbital evolution with 3.5PN, we have simulated the MHD evolution of a circumbinary disk surrounding an equal-mass non-spinning binary. Prior to the beginning of the inspiral, the structure of the circumbinary disk is predicted well by extrapolation from Newtonian results. The binary opens a low-density gap whose radius is roughly two binary separations, and matter piles up at the outer edge of this gap as inflow is retarded by torques exerted by the binary; nonetheless, the accretion rate is diminished relative to its value at larger radius by only about a factor of 2. During inspiral, the inner edge of the disk at first moves inward in coordination with the shrinking binary, but as the orbital evolution accelerates, the rate at which the inner edge moves toward smaller radii falls behind the rate of binary compression. In this stage, the rate of angular momentum transfer from the binary to the disk slows substantially, but the net accretion rate decreases by only 10-20%. When the binary separation is tens of gravitational radii, the rest-mass efficiency of disk radiation is a few percent, suggesting that supermassive binary black holes in galactic nuclei could be very luminous at this stage of their evolution. If the luminosity were optically thin, it would be modulated at a frequency that is a beat between the orbital frequency of the disk's surface density maximum and the binary orbital frequency. However, a disk with sufficient surface density to be luminous should also be optically thick; as a result, the periodic modulation may be suppressed.

We describe the Einstein Toolkit, a community-driven, freely accessible computational infrastructure intended for use in numerical relativity, relativistic astrophysics, and other applications. The Toolkit, developed by a collaboration involving researchers from multiple institutions around the world, combines a core set of components needed to simulate astrophysical objects such as black holes, compact objects, and collapsing stars, as well as a full suite of analysis tools. The Einstein Toolkit is currently based on the Cactus Framework for high-performance computing and the Carpet adaptive mesh refinement driver. It implements spacetime evolution via the BSSN evolution system and general-relativistic hydrodynamics in a finite-volume discretization. The toolkit is under continuous development and contains many new code components that have been publicly released for the first time and are described in this article. We discuss the motivation behind the release of the toolkit, the philosophy underlying its development, and the goals of the project. A summary of the implemented numerical techniques is included, as are results of numerical test covering a variety of sample astrophysical problems.