We describe early success in the evolution of binary black-hole spacetimes with a numerical code based on a generalization of harmonic coordinates. Indications are that with sufficient resolution this scheme is capable of evolving binary systems for enough time to extract information about the orbit, merger, and gravitational waves emitted during the event. As an example we show results from the evolution of a binary composed of two equal mass, nonspinning black holes, through a single plunge orbit, merger, and ringdown. The resultant black hole is estimated to be a Kerr black hole with angular momentum parameter a 0:70. At present, lack of resolution far from the binary prevents an accurate estimate of the energy emitted, though a rough calculation suggests on the order of 5% of the initial rest mass of the system is radiated as gravitational waves during the final orbit and ringdown.

The gravitational wave observations GW150914 and GW151226 by Advanced LIGO provide the first opportunity to learn about physics in the extreme gravity environment of coalescing binary black holes. The LIGO Scientific Collaboration and the Virgo Collaboration have verified that this observation is consistent with Einstein's theory of General Relativity, constraining the presence of certain parametric anomalies in the signal. This paper expands their analysis to a larger class of anomalies, highlighting the inferences that can be drawn on non-standard theoretical physics mechanisms that could otherwise have affected the observed signals. We find that these gravitational wave events constrain a plethora of mechanisms associated with the generation and propagation of gravitational waves, including the activation of scalar fields, gravitational leakage into large extra dimensions, the variability of Newton's constant, the speed of gravity, a modified dispersion relation, gravitational Lorentz violation and the strong equivalence principle. Though other observations limit many of these mechanisms already, GW150914 and GW151226 are unique in that they are direct probes of dynamical strong-field gravity and of gravitational wave propagation. We also show that GW150914 constrains inferred properties of exotic compact object alternatives to Kerr black holes. We argue, however, that the true potential for GW150914 to both rule out exotic objects and constrain physics beyond General Relativity is severely limited by the lack of understanding of the coalescence regime in almost all relevant modified gravity theories. This event thus significantly raises the bar that these theories have to pass, both in terms of having a sound theoretical underpinning, and reaching the minimal level of being able to solve the equations of motion for binary merger events. We conclude with a discussion of the additional inferences that can be drawn if the lower-confidence observation of an electromagnetic counterpart to GW150914 holds true, or such a coincidence is observed with future events; this would provide dramatic constraints on the speed of gravity and gravitational Lorentz violation.

We investigate the dynamics and gravitational-wave (GW) emission in the binary merger of equal-mass black holes as obtained from numerical relativity simulations. The simulations were performed with an evolution code based on generalized harmonic coordinates developed by Pretorius, and used quasiequilibrium initial-data sets constructed by Cook and Pfeiffer. Results from the evolution of three sets of initial data are explored in detail, corresponding to different initial separations of the black holes, and exhibit between 2 -8 GW cycles before coalescence. We find that to a good approximation the inspiral phase of the evolution is quasicircular, followed by a ''blurred, quasicircular plunge'' lasting for about 1-1.5 GW cycles. After this plunge the GW frequency decouples from the orbital frequency, and we define this time to be the start of the merger phase. Roughly 10 -15 M separates the time between the beginning of the merger phase and when we are able to extract quasinormal ring-down modes from gravitational waves emitted by the newly formed black hole. This suggests that the merger lasts for a correspondingly short amount of time, approximately 0.5-0.75 of a full GW cycle. We present first-order comparisons between analytical models of the various stages of the merger and the numerical results-more detailed and accurate comparisons will need to await numerical simulations with higher accuracy, better control of systemic errors (including coordinate artifacts), and initial configurations where the binaries are further separated. During the inspiral, we find that if the orbital phase is well modeled, the leading order Newtonian quadrupole formula is able to match both the amplitude and phase of the numerical GW quite accurately until close to the point of merger. We provide comparisons between the numerical results and analytical predictions based on the adiabatic post-Newtonian (PN) and nonadiabatic resummed-PN models (effective-one-body and Padé models). For all models considered, 3PN and 3.5PN orders match the inspiral numerical data the best. From the ring-down portion of the GW, we extract the fundamental quasinormal mode and several of the overtones. Finally, we estimate the optimal signal-to-noise ratio (SNR) for typical binaries detectable by GW experiments. We find that, when the merger and ring-down phases are included, binaries with total mass larger than 40M (sources for ground-based detectors) are brought in band and can be detected with signal-to-noise up to 15 at 100 Mpc, whereas for binaries with total mass larger than 2 10 6 M (sources for space-based detectors) the SNR can be 10 4 at 1 Gpc.