Periastron advance in black-hole binaries

The general relativistic (Mercury-type) periastron advance is calculated here for the first time with exquisite precision in full general relativity. We use accurate numerical relativity simulations of spinless black-hole binaries with mass ratios 1/8≤m(1)/m(2)≤1 and compare with the predictions of several analytic approximation schemes. We find the effective-one-body model to be remarkably accurate and, surprisingly, so also the predictions of self-force theory [replacing m(1)/m(2)→m(1)m(2)/(m(1)+m(2))(2)]. Our results can inform a universal analytic model of the two-body dynamics, crucial for ongoing and future gravitational-wave searches.     

Accurate evolutions of orbiting black-hole binaries without excision

We present a new algorithm for evolving orbiting black-hole binaries that does not require excision or a corotating shift. Our algorithm is based on a novel technique to handle the singular puncture conformal factor. This system, based on the Baumgarte-Shapiro-Shibata-Nakamura formulation of Einstein's equations, when used with a "precollapsed" initial lapse, is nonsingular at the start of the evolution and remains nonsingular and stable provided that a good choice is made for the gauge. As a test case, we use this technique to fully evolve orbiting black-hole binaries from near the innermost stable circular orbit regime. We show fourth-order convergence of waveforms and compute the radiated gravitational energy and angular momentum from the plunge. These results are in good agreement with those predicted by the Lazarus approach.     

Exploring intermediate and massive black-hole binaries with the Einstein Telescope

We discuss the capability of a third-generation ground-based detector such as the Einstein Telescope (ET) to enhance our astrophysical knowledge through detections of gravitational waves emitted by binaries including intermediate-mass and massive black holes. The design target for such instruments calls for improved sensitivity at low frequencies, specifically in the ~ 1-10{\sim 1-10} Hz range. This will allow the detection of gravitational waves generated in binary systems containing black holes of intermediate mass, ~ 100-10000  M_\odot{\sim100-10000\,\,M_{\odot}} . We primarily discuss two different source types—mergers between two intermediate mass black holes (IMBHs) of comparable mass, and intermediate-mass-ratio inspirals (IMRIs) of smaller compact objects with mass ~ 1-10  M_\odot{\sim1-10\,\,M_{\odot}} into IMBHs. IMBHs may form via two channels: (i) in dark matter halos at high redshift through direct collapse or the collapse of very massive metal-poor Population III stars, or (ii) via runaway stellar collisions in globular clusters. In this paper, we will discuss both formation channels, and both classes of merger in each case. We review existing rate estimates where these exist in the literature, and provide some new calculations for the approximate numbers of events that will be seen by a detector like the Einstein Telescope. These results indicate that the ET may see a few to a few thousand comparable-mass IMBH mergers and as many as several hundred IMRI events per year. These observations will significantly enhance our understanding of galactic black-hole growth, of the existence and properties of IMBHs and of the astrophysics of globular clusters. We finish our review with a discussion of some more speculative sources of gravitational waves for the ET, including hypermassive white dwarfs and eccentric stellar-mass compact-object binaries.     

Numerical simulations of black-hole binaries and gravitational wave emission

We review recent progress in numerical relativity simulations of black-hole (BH) spacetimes. Following a brief summary of the methods employed in the modeling, we summarize the key results in two major areas of BH physics: (i) BHs as sources of gravitational waves (GWs) and (ii) astrophysical systems involving BHs. We conclude with a list of the most urgent tasks for numerical relativity in these areas.     

Regularization of spherical and axisymmetric evolution codes in numerical relativity

Several interesting astrophysical phenomena are symmetric with respect to the rotation axis, like the head-on collision of compact bodies, the collapse and/or accretion of fields with a large variety of geometries, or some forms of gravitational waves. Most current numerical relativity codes, however, cannot take advantage of these symmetries due to the fact that singularities in the adapted coordinates, either at the origin or at the axis of symmetry, rapidly cause the simulation to crash. Because of this regularity problem it has become common practice to use full-blown Cartesian three-dimensional codes to simulate axi-symmetric systems. In this work we follow a recent idea of Rinne and Stewart and present a simple procedure to regularize the equations both in spherical and axi-symmetric spaces. We explicitly show the regularity of the evolution equations, describe the corresponding numerical code, and present several examples clearly showing the regularity of our evolutions.     

Consistency of post-Newtonian waveforms with numerical relativity

General relativity predicts the gravitational wave signatures of coalescing binary black holes. Explicit waveform predictions for such systems, required for optimal analysis of observational data, have so far been achieved primarily using the post-Newtonian (PN) approximation. The quality of this treatment is unclear, however, for the important late-inspiral portion. We derive late-inspiral waveforms via a complementary approach, direct numerical simulation of Einstein's equations. We compare waveform phasing from simulations of the last approximately 14 cycles of gravitational radiation from equal-mass, nonspinning black holes with the corresponding 2.5PN, 3PN, and 3.5PN orbital phasing. We find phasing agreement consistent with internal error estimates for either approach, suggesting that PN waveforms for this system are effective until the last orbit prior to final merger.     

Status of numerical relativity

I describe the current status of numerical relativity from my personal point of view. Here, I focus mainly on explaining the numerical implementations necessary for simulating general relativistic phenomena such as the merger of compact binaries and stellar collapse, emphasizing the well-developed current status of such implementations that enable simulations for several astrophysical phenomena. Some of our latest results for simulation of binary neutron star mergers are briefly presented.     

Hangup kicks: still larger recoils by partial spin-orbit alignment of black-hole binaries

We revisit the scenario of the gravitational radiation recoil acquired by the final remnant of a black-hole-binary merger by studying a set of configurations that have components of the spin both aligned with the orbital angular momentum and in the orbital plane. We perform a series of 42 new full numerical simulations for equal-mass and equal-spin-magnitude binaries. We extend previous recoil fitting formulas to include nonlinear terms in the spins and successfully include both the new and known results. The new predicted maximum velocity approaches 5000 km/s for spins partially aligned with the orbital angular momentum, which leads to an important increase of the probabilities of large recoils in generic astrophysical mergers. We find non-negligible probabilities for recoils of several thousand km/s from accretion-aligned binaries.     

Gravitational waves from sub-lunar-mass primordial black-hole binaries: a new probe of extradimensions

In many brane world models, gravity is largely modified at the electroweak scale approximately 1 TeV. In such models, primordial black holes (PBHs) with a lunar mass M approximately 10(-7)M([circle dot]) might have been produced when the temperature of the Universe was at approximately 1 TeV. If a significant fraction of the dark halo of our galaxy consists of these lunar mass PBHs, a huge number of BH binaries will exist in our neighborhood. Third generation detectors such as EURO can detect gravitational waves from these binaries, and can also determine their chirp mass. With a new detector designed to be sensitive at high frequency bands greater, similar 1 kHz, the existence of extradimensions could be confirmed.     

Formation and evolution of X-ray binaries

We review recent progress in theoretical understanding of X-ray binaries, which has largely been driven by new observations. We select several topics including formation of compact low-mass X-ray binaries, the evolutionary connection between low-mass X-ray binaries and binary and millisecond radio pulsars, and ultraluminous X-ray sources, to illustrate the interplay between theories and observations.     

A radiation scalar for numerical relativity

This Letter describes a scalar curvature invariant for general relativity with a certain, distinctive feature. While many such invariants exist, this one vanishes in regions of space-time which can be said unambiguously to contain no gravitational radiation. In more general regions which incontrovertibly support nontrivial radiation fields, it can be used to extract local, coordinate-independent information partially characterizing that radiation. While a clear, physical interpretation is possible only in such radiation zones, a simple algorithm can be given to extend the definition smoothly to generic regions of space-time.