A fast stroboscopic spectral method for rotating systems in numerical relativity

We present a numerical technique for solving evolution equations, as the wave equation, in the description of rotating astrophysical compact objects in comoving coordinates, which avoids the problems associated with the light cylinder. The technique implements a fast spectral matching between two domains in relative rotation: an inner spherical domain, comoving with the sources and lying strictly inside the light cylinder, and an outer inertial spherical shell. Even though the emphasis is placed on spectral techniques, the matching is independent of the specific manner in which equations are solved inside each domain, and can be adapted to different schemes. We illustrate the strategy with some simple but representative examples.Comment: 16 pages, 15 figure

Status of black-hole-binary simulations for gravitational-wave detection

It is now possible to theoretically calculate the gravitational-wave signal from the inspiral, merger and ringdown of a black-hole-binary system. The late inspiral, merger and ringdown can be calculated in full general relativity using numerical methods. The numerical waveforms can then be either stitched to inspiral waveforms predicted by approximation techniques (in particular post-Newtonian calculations) that start at an arbitrarily low frequency, or used to calibrate free parameters in analytic models of the full waveforms. In this review I summarize the status of numerical-relativity (NR) waveforms that include at least ten cycles of the dominant mode of the GW signal before merger, which should be long enough to produce accurate, complete waveforms for GW observations.Comment: 22 pages, 6 figures, Version to be published in CQG, NRDA 2008 Special Issu

Binary black hole merger in the extreme mass ratio limit

We discuss the transition from quasi-circular inspiral to plunge of a system of two nonrotating black holes of masses $m_1$ and $m_2$ in the extreme mass ratio limit $m_1m_2\ll (m_1+m_2)^2$. In the spirit of the Effective One Body (EOB) approach to the general relativistic dynamics of binary systems, the dynamics of the two black hole system is represented in terms of an effective particle of mass $\mu\equiv m_1m_2/(m_1+m_2)$ moving in a (quasi-)Schwarzschild background of mass $M\equiv m_1+m_2$ and submitted to an ${\cal O}(\mu)$ radiation reaction force defined by Pad\'e resumming high-order Post-Newtonian results. We then complete this approach by numerically computing, \`a la Regge-Wheeler-Zerilli, the gravitational radiation emitted by such a particle. Several tests of the numerical procedure are presented. We focus on gravitational waveforms and the related energy and angular momentum losses. We view this work as a contribution to the matching between analytical and numerical methods within an EOB-type framework.Comment: 14 pages, six figures. Revised version. To appear in the CQG special issue based around New Frontiers in Numerical Relativity conference, Golm (Germany), July 17-21 200

Characteristic extraction in numerical relativity: binary black hole merger waveforms at null infinity

The accurate modeling of gravitational radiation is a key issue for gravitational wave astronomy. As simulation codes reach higher accuracy, systematic errors inherent in current numerical relativity wave-extraction methods become evident, and may lead to a wrong astrophysical interpretation of the data. In this paper, we give a detailed description of the Cauchy-characteristic extraction technique applied to binary black hole inspiral and merger evolutions to obtain gravitational waveforms that are defined unambiguously, that is, at future null infinity. By this method we remove finite-radius approximations and the need to extrapolate data from the near zone. Further, we demonstrate that the method is free of gauge effects and thus is affected only by numerical error. Various consistency checks reveal that energy and angular momentum are conserved to high precision and agree very well with extrapolated data. In addition, we revisit the computation of the gravitational recoil and find that finite radius extrapolation very well approximates the result at \scri. However, the (non-convergent) systematic differences to extrapolated data are of the same order of magnitude as the (convergent) discretisation error of the Cauchy evolution hence highlighting the need for correct wave-extraction.Comment: 41 pages, 8 figures, 2 tables, added references, fixed typos. Version matches published version

Summary of sessions B1/B2 and B2: relativistic astrophysics and numerical relativity

The numerical relativity session at GR18 was dominated by physics results on binary black hole mergers. Several groups can now simulate these from a time when the post-Newtonian equations of motion are still applicable, through several orbits and the merger to the ringdown phase, obtaining plausible gravitational waves at infinity, and showing some evidence of convergence with resolution. The results of different groups roughly agree. This new-won confidence has been used by these groups to begin mapping out the (finite dimensional) initial data space of the problem, with a particular focus on the effect of black hole spins, and the acceleration by gravitational wave recoil to hundreds of km s?1 of the final merged black hole. Other work was presented on a variety of topics, such as evolutions with matter, extreme mass ratio inspirals and technical issues such as gauge choices