for almost a century now after its discovery , the schwarzschild solution describing a static , spherically symmetric vacuum spacetime in einstein s theory of general relativity ( gr ) has been a valuable tool for a wealth of theoretical and experimental studies , including in particular weak - field tests of the theory . for more than half a century ,
the validity of the solution beyond the schwarzschild radius @xmath0 and the consequent existence of bhs was thought to be a mathematical artifact of the field equations .
this picture has changed dramatically in more recent decades .
crucial events in this sense were the mathematical discovery of the kerr solution describing rotating bhs @xcite , and the astrophysical x - ray and radio observations in the 1960s and 1970s .
cosmic x - ray sources , first detected in 1963 @xcite , were soon associated with compact stellar objects and , in particular , with mass exchange from an ordinary star onto a compact object in binary systems @xcite .
this interpretation received further support when the x - ray source cyg x-1 was identified with the supergiant boib star hd  226868 @xcite , which established cyg x-1 as the first observed x - ray binary . over the next few years ,
the explanation of these systems in terms of accretion onto compact stellar objects in binary systems became widely accepted ( see @xcite for a review ) .
measurements of the binary period and radial velocity provide lower limits on the mass of the compact companion .
the values obtained for some x - ray binaries comfortably exceed the upper mass limit of about @xmath1 for neutron stars .
recent investigations of cyg x-1 , for example , predict @xmath2 for the compact member @xcite . also starting in the early 1960s , optical and radio observations identified `` star like '' objects with surprisingly large redshifts @xcite .
if interpreted as cosmological in origin , the redshifts implied enormous distances and , accordingly , luminosities . over the next decade
the cosmological origin of these `` quasi - stellar '' sources ( or quasars ) became clear , and accretion onto supermassive bhs ( smbhs ) @xcite was accepted as the most plausible explanation of their energetics @xcite . by
now , observations of stellar dynamics near the centers of galaxies and iron k@xmath3 emission line profiles provide strong evidence that many galaxies harbor bhs at their centers @xcite .
bh masses are generally closely correlated to the bulge velocity dispersion and luminosity @xcite . for the milky way
, astrometric and radial velocity observations of short - period stars orbiting the galactic center @xcite predict an unseen gravitational mass of @xmath4 .
electromagnetic observations provide strong but _ indirect _ evidence for the existence of astrophysical bhs .
a more direct way of probing the bh nature of astrophysical compact objects is offered by the newly emerging field of gw astronomy .
first - generation ground - based laser - interferometric detectors ( such as ligo , virgo and geo600 ) have reached design sensitivity , and advanced versions of these detectors will be operational within a few years @xcite .
one of the most prominent sources of detectable gws for these detectors is the inspiral and coalescence of stellar mass and , possibly , intermediate - mass bh ( imbh ) binaries . at lower frequencies , space - based laser interferometric observations
@xcite have the potential to complement ground - based efforts with high signal - to - noise ratio observations of massive bh binaries .
space - based detectors are guaranteed to observe galactic binaries containing white dwarfs and/or neutron stars , and they may also detect the inspiral of stellar compact objects into massive bhs . due to the weak interaction of gws with any type of matter , including the detectors , digging physical signals from the noisy data stream represents a daunting task and heavily relies on so - called _ matched filtering _
techniques @xcite .
matched filtering is a common choice to search for signals of known form in noisy data , and works by cross - correlating the actual `` signal '' ( i.e. , the detector s output ) with a set of theoretical templates .
accurate modeling of the gw sources thus plays a vital role in maximizing the scientific output from these experimental efforts .    event rate estimates for imbhs and stellar - mass bh binaries are very uncertain , but advanced earth - based detectors are expected to observe several events per year ( see @xcite for a review ) .
_ noted that recent metallicity measurements require revisions in the population synthesis models used for event rate estimates , suggesting that stellar - mass bh - bh binaries may be the first gw sources to be detected @xcite .
the observations of two possible x - ray binary precursors of bh - bh binaries provide a much needed observational constraint on compact binaries containing at least one bh , suggesting again that bh - bh binary mergers may be more common than anticipated @xcite .
population synthesis models have large uncertainties , related e.g. to the poorly known common envelope phase , and some models would necessarily be ruled out if advanced ligo does _ not _ detect bh binaries . in this sense , gw detectors are _ guaranteed _ to put constraints on compact binary formation in the near future .
schutz @xcite also pointed out that coherent data analysis using three or more detectors may further increase event rates by factors of a few , making prospects of doing astrophysics with earth - based detectors even brighter .
estimates for event rates of smbh binaries detectable by space - based interferometers range from @xmath5 to @xmath6s per year , with a most likely estimate of @xmath7 @xcite .    over the past decade , bh modeling has also attracted a great deal of interest in the context of high - energy physics . according to the gauge / gravity duality , gravitational physics in anti
- de sitter ( ads ) backgrounds describes field theories on the boundary of that spacetime @xcite . according to this duality , bhs in ads
are dual to an equilibrium field theory , and therefore interacting bhs or bhs off - equilibrium may represent interesting physics from the dual point of view .
in fact , there is some evidence that some of the physics behind quark - gluon plasma formation in heavy - ion collisions is captured by this duality , which on the gravitational side corresponds to shock wave collisions and bh formation in ads @xcite .
understanding bh spacetimes is also of interest for high - energy physics in connection with extra - dimensional scenarios .
solutions to the hierarchy problem in physics have been developed where all standard - model interactions are confined to a 3 + 1 dimensional brane embedded in a spacetime with large extra dimensions or an extra dimension with a warp factor @xcite . in these models gravity
propagates in the entire higher - dimensional spacetime ( the `` bulk '' ) , and the fundamental planck scale could be as low as @xmath8 .
these energy scales are accessible to experiments such as the large hadron collider ( lhc ) .
if the fundamental planck scale is so low , thorne s hoop conjecture suggests the exciting possibility of bh production in high - energy parton - parton collisions @xcite , i.e.  a direct observable signature of the existence of extra - dimensions .
accurate modeling of energy and momentum losses in the form of gws during the collisions and determination of the bh formation cross section will be crucial for the analysis of data from such experiments ( see e.g.  @xcite ) , and they should be amenable to classical calculations in @xmath9-dimensional gr @xcite .    finally , the understanding of bh dynamics under extreme conditions is fascinating from a mathematical - physics point of view .
questions such as bh stability , cosmic censorship , etcetera can only be addressed by solving the full nonlinear set of the field equations .    for many of these scenarios ,
it is convenient to divide the coalescence of a bh binary in the framework of gr into three stages which mirror the different tools used for their modeling : ( i ) the inspiral phase , where the interaction between the individual holes is still sufficiently weak to justify the use of post - newtonian ( pn ) techniques @xcite ; ( ii ) the final orbits , plunge and merger , which are governed by the strong - field regime of einstein s equations and can be described only by fully numerical simulations ; ( iii ) the ringdown of the remnant bh , which is amenable to a perturbative treatment @xcite .    in this article
we will focus on the second , strong - field stage , on the numerical tools dedicated to its study , and results obtained following the 2005 breakthroughs @xcite . as we shall see , however , a comprehensive understanding of bh dynamics requires a close interplay of numerical and analytical studies , and we will discuss in various cases the interface between numerical relativity and approximation techniques .    following a brief summary of the framework employed in numerical relativity ( sec .
[ sec : numerics ] ) , in secs .
[ sec : gws]-[sec : highenergy ] we will present the main results obtained from numerical studies of bhs in the context of gw detection , astrophysics and high - energy physics , respectively .
we conclude in sec .
[ sec : conclusions ] with a discussion of future applications of numerical relativity . at the end of each section
we will provide references to the literature for further reading .
* _ notation : _ * we shall be using geometrical units , such that the gravitational constant @xmath10 and the speed of light @xmath11 .
we denote spacetime indices @xmath12 by greek letters , and spatial indices @xmath13 by latin letters .
in einstein s gr , spacetime is modelled as a @xmath9-dimensional manifold with a metric @xmath14 of lorentzian signature @xmath15 . with the exception of sec .  [
sec : highenergy ] we will restrict our discussion to @xmath16 , i.e. three spatial and one time dimension .
the metric is determined by the einstein equations @xmath17 with @xmath18 in vacuum . while analytic solutions ( as , for example , the schwarzschild and kerr solutions and the friedmann - lematre - robertson - walker solution ) have been found for systems with high symmetry , solutions for dynamical configurations without special symmetries generally require the use of numerical tools . for this purpose
the einstein equations need to be cast as an _ initial value formulation _ such that the cauchy - kowalewski theorem guarantees a uniquely determined evolution of appropriately chosen initial data . in the case of einstein s equations
, it is not obvious that such a formulation can be obtained : a straightforward calculation shows that eq .
( [ eq : einstein ] ) contains second time derivatives of the metric components @xmath19 , but not of @xmath20 .
this fact mirrors the coordinate or gauge freedom characteristic of gr , and implies the existence of so - called _ constraint _ equations which impose conditions on the metric components and their first derivatives on each spatial slice with constant coordinate time @xmath21 . in the following
we will discuss two approaches which provide an initial value formulation of the einstein equations and have led to successful , long - term stable numerical evolutions : the generalized harmonic ( gh ) formulation employed in pretorius breakthrough bh binary simulations @xcite and the canonical arnowitt - deser - misner ( adm ) split @xcite , reformulated by york @xcite , which forms the starting point for the _ moving puncture _
breakthroughs @xcite .    * _ adm 3 + 1 split : _ * spacetime is decomposed into a one - parameter family of three - dimensional spatial slices .
we choose coordinates adapted to this formulation , such that the coordinate time @xmath21 labels each slice and @xmath22 denotes points on the slice . on each hypersurface
, there exists a unique timelike , future pointing unit normal field @xmath23 which defines a projection operator @xmath24 onto the hypersurface .
the geometry of the hypersurface is completely determined by the _ three - metric _ or _ first fundamental form _
the coordinate freedom is represented by the shift vector @xmath26 and the lapse function @xmath27 ( with @xmath28 ) , which relate spatial coordinates on ( and measure separation in proper time between ) different hypersurfaces .
the projections @xmath29 , @xmath30 , @xmath31 of the einstein equations then lead to six evolution equations for the three - metric @xmath32 , three momentum constraints and the hamiltonian constraint : @xmath33 where @xmath34 , @xmath35 and @xmath36 denote the ricci tensor , ricci scalar and covariant derivative associated with the three - metric , and the extrinsic curvature @xmath37 has been introduced to obtain a first - order system in time .
it turns out that this formulation of the einstein equations is only _ weakly hyperbolic _
@xcite , and therefore not suitable for long - term stable numerical evolutions . motivated by these stability problems , several modifications of the adm system have been investigated .
the most popular is the baumgarte - shapiro - shibata - nakamura ( bssn ) formulation @xcite , which rearranges the degrees of freedom via a split of the extrinsic curvature into trace and trace - free parts , a conformal transformation and promotion of the contracted conformal christoffel symbols to the status of independent variables : @xmath38 the resulting set of evolution and constraint equations is given in eqs .
( 15)-(30 ) of @xcite .
subject to minor modifications , such as replacing the conformal factor by @xmath39 or adding the auxiliary constraint arising from the definition of @xmath40 in eq .
( [ eq : bssn ] ) , this formulation is employed in the present generation of moving puncture evolution codes @xcite which fix the gauge via `` 1+log '' slicing and the @xmath41-driver condition @xcite .
* _ gh formulation : _ * harmonic coordinates are defined by the condition @xmath42 , where @xmath43 represents the scalar wave operator .
these coordinates have played an important role in the analysis of the cauchy problem in gr @xcite . in harmonic coordinates ,
the einstein equations take on a manifestly hyperbolic form , which allows for a generalization to arbitrary coordinate or gauge choices @xcite .
the first step is to introduce four arbitrary source functions @xmath44 such that the coordinates obey @xmath45 and treat these functions as independent evolution variables .
then one considers the gh system @xmath46 where @xmath47 .
( [ eq : gh ] ) is equivalent to the einstein equations , subject to the validity of the constraint ( [ eq : h ] ) . in expanded form ,
the gh system is given by @xmath48 and the constraints @xmath49 are preserved by virtue of the bianchi identities provided that @xmath50 and @xmath51 vanish on the initial hypersurface @xmath52 . a key ingredient in the numerical evolution of the gh system is the addition of a constraint damping term @xmath53 c_{\mu}$ ] to the right hand side of eq .
( [ eq : ghexpanded ] ) @xcite . here
, @xmath54 is a non - vanishing timelike vector field . finally , it is interesting to note that the hamiltonian and momentum constraints are automatically satisfied on the initial hypersurface if @xmath55 @xcite .
the coordinates were determined in pretorius initial simulations by @xmath56 and the evolution of the metric proceeds according to eq .
( [ eq : ghexpanded ] ) with the constraint damping term of gundlach _ et al .
_ @xcite . for further discussion of gauge choices in the ghg system , see also @xcite .
the spectral code originally developed by the caltech - cornell group ( see @xcite and references therein ) employs a first - order version of the gh system @xcite with dual - coordinate frames and bh excision .
furthermore , the first - order gh formulation facilitates the specification of constraint - preserving boundary conditions @xcite . during the inspiral phase ,
they evolve the gauge functions @xmath57 such that they remain constant in a frame comoving with the bhs , but switch to a dynamic evolution during the plunge and ringdown ; see sec .
iii in @xcite .
while more complex in structure , this framework enables their code to generate bh evolutions with exceptional accuracy @xcite .
* _ initial data : _ * the construction of initial data requires solving the hamiltonian and momentum constraints . here
we only summarize the key concepts ; we refer the reader to cook s review article @xcite for details .
most work on the construction of initial data is based on the york - lichnerowicz split @xcite , which involves a conformal transformation of metric and extrinsic curvature and separates the latter into trace and trace - free part .
more recently , the _ thin - sandwich _
construction @xcite , which replaces the extrinsic curvature by the time derivative of the metric , has become a popular alternative . in either case , the resulting elliptic equations simplify substantially under the assumption of conformal flatness and a spatially constant @xmath58 .
while this formalism provides a convenient method to solve the constraint equations , we still need to ensure that the initial data represent a physically realistic system , typically two bhs with specific spins and momenta .
this can be achieved by generalizing the schwarzschild solution , which is obtained in the above framework in conformally flat form with conformal factor @xmath59 .
a generalization to initial data of @xmath60 bhs starting from rest is directly obtained by the construction of misner @xcite or brill and lindquist @xcite .
remarkably , an analytic solution for the extrinsic curvature can still be obtained for bhs with initial linear momenta @xmath61 and spins @xmath62 @xcite . by applying a compactification to the internal asymptotically flat region ,
brandt & brgmann @xcite arrived at the so - called _ puncture data _ construction , which is the method of choice for most numerical evolutions employing the bssn formulation . in order to overcome an upper limit of @xmath63 for the dimensionless spin of bhs in
puncture type initial data @xcite , lovelace _ et al . _
@xcite generated initial data based on a generalization of the kerr - schild solution @xcite and were thus able to evolve a bh binary with spin magnitude 0.95 .
in alternative to generalizing analytically known bh solutions , the presence of horizons in the initial data can be encoded in the form of boundary conditions for the metric and extrinsic curvature @xcite as determined by the isolated horizon framework @xcite .
initial data along these lines have been constructed in @xcite , and form the starting point of most of the numerical evolutions using the gh system ( see e.g. @xcite ) .    * _ diagnostics : _ * extracting physical information from numerical simulations of the einstein equations is nontrivial for two reasons . first , the evolved variables are dependent on the coordinate conditions ; second , it is often not possible to define local quantities familiar from newtonian physics .
the first difficulty requires the calculation of gauge - independent variables .
the second difficulty is alleviated by the _
isolated horizon _
framework @xcite , which facilitates the calculation of bh mass and spin in the limit of isolated bhs ; in practice , this framework is often applied when the bhs are farther apart and their interaction is considered negligible .
local properties of the bhs are encoded in their apparent horizon @xcite . in particular , the irreducible mass can be expressed in terms of the apparent horizon area : @xmath64 . in the limit of an isolated bh , the spin can be derived from the integration of the rotational killing vectors over the horizon according to eq .
( 8) of @xcite .
in practice it has been found that using _ flat - space _ rotational killing vectors yields reasonable results @xcite , but we also note the improved method to compute approximations to killing vectors in @xcite and an alternative method @xcite based on weyl scalars to compute quasi - local bh quantities without solving the killing equation .
it is important to bear in mind , however , the approximate nature of any spin calculation in bh binary simulations .    in comparison , it is more straightforward to extract global quantities of the spacetime , most notably the total mass - energy @xmath65 and the linear and angular momentum @xmath66 , @xmath67 @xcite .
these are given by surface integrals at spatial infinity , see e.g.  eqs .
( 7.15 ) , ( 7.56 ) and ( 7.63 ) in @xcite . arguably the most important diagnostic quantity in bh simulations is the gw signal . the most common method to extract gws
is based on the newman - penrose formalism @xcite and derives the complex newman - penrose scalar @xmath68 from contraction of the weyl tensor with suitably chosen tetrad vectors @xcite , sec .
iii a in @xcite .
it is often convenient to decompose @xmath68 into multipoles @xmath69 using spherical harmonics of spin - weight @xmath70 : @xmath71 .
contemporary numerical codes typically evaluate @xmath68 at finite coordinate radius .
this results in systematic errors , due to ambiguities in the tetrad choice and the neglect of gw backscattering .
while these errors can be reduced by extracting gws at a sequence of radii and extrapolating to infinite radius ( @xmath68 asymptotically falls off @xmath72 @xcite ) , a cleaner method is to calculate @xmath68 at infinity as provided by _
cauchy - characteristic extraction _
methods @xcite . from @xmath68 ,
one straightforwardly obtains the energy , linear and angular momentum radiated in the form of gws ; see e.g.  eqs .
( 49)-(51 ) in @xcite . in gw data analysis
it is more common to work with the wave strain @xmath73 decomposed into multipoles : @xmath74 . in practice , after integrating @xmath68 twice in time the signal can be severely affected by nonlinear drifts , which can be controlled to a significant extent by performing the integration in the frequency domain @xcite .
alternatively , gws can be extracted by viewing the metric in the far field regime as a perturbation of the schwarzschild spacetime and employing the formalism of regge , wheeler @xcite and zerilli @xcite .
one thus obtains the regge - wheeler - moncrief and zerilli - moncrief master functions @xmath75 , @xmath76 , which can be converted into the multipolar components of the gw strain @xmath77 according to eq .
( 49 ) of @xcite .
for a comprehensive summary of the numerical framework for evolving einstein s field equations , we refer the reader to the books by alcubierre @xcite and baumgarte and shapiro @xcite .
further details , including mathematical aspects of the equations and characteristic techniques , can be found in the review articles @xcite .
bh binary systems represent one of the most promising sources of detectable gws .
the parameters of a bh binary are commonly divided into _ intrinsic _ and _ extrinsic _ @xcite .
intrinsic parameters characterize physical properties of the system , such as the total mass @xmath78 , the mass ratio @xmath79 , the individual bh spins @xmath80 , @xmath81 and the orbital eccentricity .
extrinsic parameters such as sky position , distance , orbital inclination angle , arrival time and initial phase of the wave , in contrast , depend on the source location relative to the observer , and do not directly enter the gw source modeling process .
the majority of numerical studies of bh binary spacetimes performed to date have focussed on comparable mass - ratios @xmath82 and moderate spin magnitudes @xmath83 , with particular emphasis on spins ( anti-)aligned with the orbital angular momentum .
we note , however , the following explorations outside this `` best charted '' subset of the parameter space .
circular binaries with mass ratios up to @xmath84 have been studied in @xcite ; comparisons for this value of @xmath85 with fully perturbative calculations have been obtained in the limit of head - on collisions in @xcite , who report a discrepancy of @xmath86 in gw energy and momentum , probably in large part due to the discretization error of the numerical simulations .
the first direct comparison of numerical results with leading - order self - force prediction was studied in the context of periastron advance in bh binaries in @xcite .
bh binaries with spins @xmath87 have been evolved in @xcite , but exceeding an apparent barrier of @xmath88 @xcite requires departure from the conformal flatness assumption @xcite .
the present maximum is @xmath89 for aligned spins evolved for 12.5 orbits by lovelace _
finally , we note that emission of gws efficiently circularizes bh binary orbits @xcite .
bh binaries are therefore expected to have vanishing eccentricity by the time they enter the frequency window of ground - based detectors , and for this reason most work on gw source modeling has focused on the quasi - circular limit .
the derivation of bh momenta , including a small radial component , for quasi - circular initial data in numerical relativity is based either on integrating the pn equations of motion @xcite or iterative procedures using several numerical simulations @xcite . the decay of eccentricity during the inspiral phase as well as periastron advance has been measured in numerical simulations in @xcite . for numerical studies of bh binaries with significant eccentricity as well as arguments why eccentric binaries may represent relevant sources of gws after all , especially in the more extreme mass ratio regime and for space - based detectors ,
we refer the reader to @xcite .    * _ gravitational waveforms from bh binaries : _ * for illustration , we display in fig .  [ fig : hlm ] the real parts of the @xmath90 and @xmath91 multipoles of the gw strain @xmath77 obtained for the inspiral and coalescence of a binary     and @xmath91 obtained from the inspiral and merger of a nonspinning binary with mass - ratio @xmath92 .
the dashed red curve represents pn predictions matched to the numerical ( solid black ) signal .
the vertical dotted lines mark the time when the frequency of the @xmath93 multipole has reached @xmath94.,height=238 ]    of nonspinning bhs with mass ratio @xmath92 @xcite .
in the course of the inspiral , both amplitude and frequency of the gw signal gradually increase .
close to merger , defined as formation of a common apparent horizon , the gw amplitude reaches a maximum .
eventually it drops exponentially as the merged hole rings down to a stationary kerr state : see e.g.  fig .
18 in @xcite .
the qualitative features of the gw signal emitted by different types of binaries can be summarized as follows .
( i ) the inspiral of two nonspinning , equal - mass bhs is the case most intensively studied by numerical relativity .
the wave signal is dominated by the @xmath93 multipolar component , which carries @xmath95 of the total radiated energy . the longest numerical simulation by scheel _
@xcite covers 16 orbits and reports a ratio of final to initial bh mass @xmath96 and a final spin @xmath97 .
( ii ) as the mass - ratio decreases , higher - order modes become more prominent @xcite and the fractional energy in @xmath93 drops to about @xmath98 in the limit @xmath99 @xcite . motivated by the strong relation of phase and frequency of different @xmath100 multipoles , baker _
_ @xcite view the @xmath100 radiation modes as being generated by the corresponding momenta of an _ implicit rotating source _ ; see also @xcite .
( iii ) gw multipoles with odd @xmath101 are suppressed in the equal - mass limit due to symmetry .
( iv ) bh binaries with spins aligned with the orbital angular momentum emit stronger gw signals due to an effect sometimes referred to as _ orbital hang - up _ : equal - mass binaries with @xmath102 radiate about twice as much energy and angular momentum when compared to the nonspinning case @xcite . for anti - aligned spins the gw energy decreases by a similar factor . ( v ) orbital eccentricity introduces a nonmonotonic behavior of the gw frequency ; cf .
fig .  6 in @xcite .
( vi ) spin - spin and spin - orbit couplings cause a precession of the orbital plane in the case of bh binaries with spins that are not aligned or anti - aligned with the orbital angular momentum .
this precession manifests itself in a modulation of the gw amplitude emitted in a fixed angular direction : cf .
fig .  1 in @xcite .
* _ tests of general relativity : _ * numerical simulations of binary bh mergers provide an opportunity to study the nonlinearities of the theory and to test the kerr nature of astrophysical bhs . if gr is the correct theory of gravitation , all bhs in the universe should be described by the kerr solution @xcite , which depends on two parameters : the mass @xmath78 and the dimensionless spin @xmath103 .
the remnant of a binary bh merger settles down to the kerr solution by emitting gravitational radiation at characteristic ( complex ) quasinormal frequencies @xmath104 that depend only on @xmath103 . here
@xmath100 are the usual angular indices , and @xmath105 is the `` overtone number '' : modes with small @xmath105 have a longer damping time and should dominate the radiation @xcite .
the measurement of the real and imaginary part of a single quasinormal mode contains , in principle , enough information to determine the mass and spin of the bh . together with the accurate determination of the individual bh masses from the inspiral phase this already allows for tests of gr , by testing the gr prediction for mass loss during inspiral and merger and at the same time the strong field dynamics of gr .
the measurement of _ any _ other frequency or damping time then provides a further interesting test of the kerr nature of the final bh @xcite .
the feasibility of these tests depends on the characteristics of a given gw detector and on the relative excitation of the modes ( see @xcite for detailed studies ) .
quantifying the excitation of quasinormal modes for generic initial data is a formidable technical problem @xcite , but numerical merger simulations allow sensible estimates of the relative quasinormal mode amplitudes @xcite .
the relative mode amplitude depends on the binary parameters , and therefore , in principle , the measurement of a multi - mode signal with single or ( better ) multiple gw detectors can be used to measure the binary mass ratio or the inclination of the final bh spin with respect to the line of sight .
we also remark that the increased signal - to - noise ratio coming from the merger has been shown to _ improve _ the bounds on alternative theories of gravity that could come from observations of the inspiral only .
keppel and ajith , for example , discussed this possibility in the context of massive graviton theories , where the graviton mass would modify the dispersion relation of gws @xcite .
formulations of the evolution equations in alternative theories of gravity are in their infancy @xcite and the numerical exploration of bh binaries in alternative theories of gravity has just started @xcite .
this is an interesting line of research , and in the near future we may have more concrete ways to quantify deviations from gr in strong - field mergers .
* _ gw template banks : _ * the key challenge in gw data analysis is to accurately predict and isolate the features of gravitational waveforms in the data stream from gw detectors .
this is necessary in order to ( i ) detect the presence of a gw signal of astrophysical origin , and ( ii ) determine the parameters of the emitting source .
these goals can be achieved by matched filtering , i.e.  by cross - correlating the data stream @xmath106 ( composed of detector noise @xmath105 plus a potential gw signal @xmath77 ) with a bank of theoretical waveform templates @xmath107 , where the indices @xmath108 @xmath109 denote the @xmath110 intrinsic and extrinsic source parameters . for this purpose
it is imperative to construct a large set of accurate numerical waveforms which cover the relevant parameter space .
this task far exceeds the capacity of purely numerical methods . if we assume a seven - dimensional parameter space
one mass ratio and three spin parameters per hole
a mere ten waveforms per parameter would correspond to @xmath111 waveforms .
furthermore , using binary waveform models including the long inspiral phase plus merger optimizes the signal - to - noise ratio in gw observations , and the inclusion of merger and ringdown has been shown to provide significant improvements in source localization and distance calculation @xcite . because current numerical simulations cover only about 30 cycles of the inspiral , the generation of complete waveforms requires the combined use of pn and numerical methods .
the construction of such gw template banks currently proceeds along either of the following two paths .
\(i ) the _ effective - one - body _ ( eob )
approach @xcite maps the dynamics of the two - body problem in gr into the motion of a particle in an effective metric .
the components of the effective metric are currently determined to 3pn order .
the eob method improves upon this model by using additional _ pseudo - pn _ terms of higher order between @xmath112 and @xmath113 @xcite . ] which are not derived from pn expressions , but calibrated via comparison with numerical results @xcite .
further improvements come from using a resummed version of the pn expanded results and from modeling nonadiabatic effects in the inspiral : see e.g.  sec .
iv in @xcite .
the inspiral - plunge waveform resulting from this construction is then matched to a merger - ringdown signal composed of a superposition of quasinormal oscillation modes of a kerr bh .
the total number of free parameters varies between the individual eob models currently investigated , but in all cases , the general strategy is to calibrate these parameters by comparison with a finite number of numerical simulations .
( ii ) _ phenomenological waveform templates _ are based on _ hybrid waveforms _ : the early inspiral is modeled by pn techniques , and the resulting signal is matched ( within a specified window ) to a numerical waveform describing the last orbits , merger and ringdown .
the resulting set of hybrid waveforms are then approximated by a model containing a number of phenomenological parameters which are mapped to the physical parameters of the binaries , such as the mass ratio @xmath85 .
the phenomenological models initially developed for nonspinning binaries in @xcite , and then extended to the case of spins aligned with the orbital angular momentum @xcite , are based on hybrid waveforms matched in the _ time domain _ , whereas the more recent study by santamaria _ et al . _
@xcite performs the matching in the frequency domain . in spite of such differences in their construction ,
the final result of all these phenomenological models are closed - form analytic expressions for the waveforms in the frequency domain .
an exploration of phenomenological models for equal - mass binaries with spin precession is given in @xcite .
both types of template banks are employed in the analysis of gw data .
the recent search for gws from binary bh inspiral , merger and ringdown @xcite used phenomenological models for injection , and an eob model for injection and for the search templates .
this search constrained the rate of mergers of binaries with individual masses in the range @xmath114~m_{\odot}$ ] to be no more than @xmath115 at 90% confidence .
numerical waveforms have also been used inside the gw community - wide _ ninja _ project @xcite to study the sensitivity of existing gw search algorithms used in the analysis of observational data . for this purpose ,
numerical relativity waveforms were injected into a simulated data stream designed to mimic ground - based detector noise .
the algorithms detected a significant fraction of the injections , but likely require further development , in particular for the purpose of measuring source parameters @xcite .
a further community wide effort , the _ nrar _ project @xcite , is dedicated to a systematic exploration of the complete bh binary parameter space to develop optimally - calibrated template families to be used in gw data analysis .
one of the most exciting prospects of gw detection is the idea of looking for electromagnetic counterparts to the gw signal .
while a _ network _ of earth - based detectors is required to localize a gw source by triangulation ( see e.g.  @xcite ) , a single space - based detectors like lisa may suffice to localize low - redshift sources and determine their luminosity distance .
the inclusion of merger and/or of higher multipolar components of the radiation in the waveform models has been shown to provide significant improvements in source localization and distance determination @xcite .
for this reason it is very important to construct phenomenological models including higher multipoles of the full inspiral / merger / ringdown waveform .
the eob model has recently been extended in this direction @xcite . * _ accuracy requirements : _ * clearly , the two approaches for template construction mentioned above require accurate numerical waveforms . quantifying these accuracy requirements
is a nontrivial task .
numerical uncertainties due to finite resolution are directly tested by _ convergence analysis _ , and the error incurred by extracting gws at finite radius can be estimated using extrapolation of results from various radii to infinity @xcite .
the uncertainties in phase and amplitude thus derived , however , depend on the alignment of the waveforms in phase and time ; see e.g.  fig .  2 in @xcite . because offsets in phase and time are free parameters in gw data analysis , it is more convenient to measure accuracy and agreement of different waveforms in terms of quantities which take such alignment into account by construction .
such a measure is obtained from the inner product between two waveforms @xmath116 and @xmath117 used in the theory of parameter estimation @xcite : @xmath118 where the tilde and asterisk denote fourier transform and complex conjugate , respectively , and @xmath119 is the one - sided power spectral density of the detector strain noise @xcite .
in practice , this inner product is to be understood as maximized over constant offsets in time and phase ( @xmath120 and @xmath121 ) between the two waveforms . the signal - to - noise ratio that would be obtained for a physical signal @xmath122 and a model waveform @xmath123
is then given by @xmath124 .
it is related to the optimal signal - to - noise ratio @xmath125 for a perfect model waveform by the _ mismatch _
@xcite : @xmath127 to leading order in @xmath128 , this definition of the mismatch implies @xcite @xmath129 where @xmath130 is that part of the waveform error parallel to @xmath122 .
based on these definitions , accuracy requirements are obtained as follows . two waveforms differing by @xmath131 will be indistinguishable for a detector if @xmath132 @xcite , such that acceptable waveform models to be used in parameter estimation must lie within a sphere of unit radius of the exact waveform . for the purpose of event detection , we note that the fractional loss in detection of gw signals due to imperfect templates is proportional to the third power of the reduction in the snr , i.e.  @xmath133 .
a value of @xmath134 is typically deemed acceptable , corresponding to a mismatch @xmath135 . for detection efficiency , there arises , however , a complication due to the discrete nature of the template bank . as discussed in detail in sec .
ii b of @xcite , the effective mismatch in a real gw search is given by a contribution due to the actual modeling procedure , plus a term which accounts for the discrete spacing in the parameter space : see their fig .  3 and eq .
template banks adopted in ligo searches use 0.3 for the latter @xcite and , in consequence , lindblom _ at al .
_ @xcite recommend a maximum mismatch @xmath136 . in terms of the right - hand side of eq .
( [ eq : mm2 ] ) , dubbed _ inaccuracy function _ in ref .
@xcite , we summarize the accuracy requirements as ( cf .
( 10 ) in ref .
@xcite ) @xmath137 accuracy measurements employing the mismatch or @xmath128 have been used to study numerical , semi - analytic and hybrid waveforms . according to the samurai project @xcite ,
relative mismatches of the @xmath138 mode between a variety of purely numerical waveforms for binaries with total mass @xmath139 are better than @xmath140 .
_ @xcite found errors in the construction of hybrid waveforms to be dominated by pn contributions .
numerical waveforms should be long enough to enable matching at lower frequencies , where pn approximations are more accurate : @xmath141 orbits before merger for the equal - mass , nonspinning case , and @xmath142 orbits for binaries with spins @xmath143 .
@xcite report significantly more demanding requirements of @xmath144 orbits for nonspinning , equal - mass binaries .
similarly , boyle @xcite concludes that longer nr waveforms or more accurate pn predictions will be required .
the discrepancy largely arises from the use of more stringent accuracy thresholds : @xmath136 in @xcite ( instead of @xmath145 in @xcite ) for detection , and @xmath146 for parameter estimation .
most recently , ohme _ et al .
_ @xcite have argued that maximising the match over intrinsic source parameters as well as phase and time shifts may reduce length requirements on numerical waveforms .
a comparison of phenomenological and eob template banks based on the inaccuracy function has been performed by damour _ et al .
_ @xcite . for advanced ground - based detectors ,
they conclude that current phenomenological models deviate from the eob ( considered as a target model ) beyond acceptable thresholds over a wide range of the parameter space , for both , detection and parameter estimation .
further reading on bh simulations in the context of gw physics is available in @xcite .
astrophysical observations have provided the strongest evidence to date that bhs exist and play an important role in many physical processes in our universe .
astrophysical observations of x - ray binaries or quasars , as mentioned in sec .
[ sec : intro ] , have a good deal to tell us about bhs . in recent years , numerical relativity simulations of bhs have also deepened our insight into a variety of astrophysical systems .
we will first summarize what we have learned from `` traditional '' electromagnetic observations about the populations of bhs we can expect to exist in the universe . then we will discuss some astrophysical implications of numerical relativity .    * _ expected bh populations : _ * we have already mentioned the two main classes of bhs identified by astrophysical observations .
( i ) solar - mass bhs with @xmath147 are usually found in x - ray binaries @xcite and are formed as the end - product of the evolution of massive stars .
( ii ) smbhs with @xmath148 are believed to reside in most active galactic nuclei ( agn ) @xcite .
the assembly of these smbhs is likely to result from a combination of bh mergers and accretion of surrounding matter over cosmological timescales .
evidence for a third class of bhs with @xmath149 ( imbhs ) is tentative at present @xcite .
gw frequencies are inversely proportional to the system s total mass , and therefore bh masses play a crucial role in binary bh detectability .
earth - based detectors , such as ligo and virgo , have an optimal sensitivity band corresponding to stellar - mass bhs and imbhs , while the planned space - based detector lisa is most sensitive to high - mass imbhs and smbhs .
as discussed in the previous section , the expected gw pattern from bh binaries depends on the masses and spins of the binary members .
we shall further see below that the gravitational recoil imparted upon the remnants of coalescing binaries strongly depends on spins and mass ratio .
so , what do `` traditional '' electromagnetic observations tell us about bh masses and spins ?
current spin estimates for accreting , stellar mass bhs are usually obtained by modeling the shape of their x - ray spectrum or by analyzing the skew in fe  k@xmath3 emission lines .
the resulting estimates are all , to some extent , model - dependent , but frequently yield moderate values @xmath150 and , for some cases , values close to the kerr limit @xmath151 ; see e.g.  @xcite and references therein .
theoretical arguments and observations suggest that stellar - mass bhs in binaries _ retain the spin they had at birth _ : neither accretion nor angular momentum extraction are likely to change significantly their mass or spin @xcite . for smbhs ,
spin estimates relying on the shape of the fe  k@xmath3 line yield values ranging from moderate @xmath152 ( e.g.  @xcite to near - critical spins ( e.g.  @xcite ) .
smbhs are expected to grow by a combination of mergers and accretion .
their spin depends sensitively on the details of these processes and of their growth @xcite .
smbh assembly via bh mergers does not appear to be able to account for the large observed spin values ; this may indicate that accretion dominates over mergers @xcite . in any case ,
bh spins encode the history of their formation , and it would be extremely useful to have detailed knowledge of the bh spin distribution function . ultimately , such a detailed census of bh parameters is one of the key targets of future gw observations @xcite .
for the purpose of numerical modeling of bhs , present measurements indicate that all spin magnitudes ( @xmath153 ) should be covered .
we have already noted that stellar - mass bhs typically have masses in the range @xmath154~m_{\odot}$ ] .
unfortunately , there are no confirmed observations of stellar - mass _ binary _ bh candidates , so estimates of mass ratios must rely on theoretical models .
for stellar - mass binaries , population synthesis codes suggest that @xmath85 should always be quite close to unity @xcite .
measurements of smbh masses are obtained by observations of stellar motion near galactic centers , as mentioned in sec .
[ sec : intro ] , as well as motion of gas discs @xcite , applications of the virial theorem to the velocity dispersion of stars @xcite and reverberation mappings applied to more distant agns @xcite .
an exhaustive list of galaxies with smbh mass measurements in the range @xmath155 is presented by graham _
the impact of different smbh assembly models on the mass and mass ratio distribution of detectable binaries has been discussed by various authors .
the general consensus is that mass ratios @xmath156 ( and down to @xmath157 ) should be common @xcite .
in contrast with the case of stellar - mass bhs , there is by now some observational evidence for smbh binaries @xcite .
these observations are not sufficient to tightly constrain smbh merger rates , but they are at least broadly consistent with merger - tree models predicting tens to hundreds of events during the typical lifetime of a space - based gw detector @xcite .
a more speculative kind of source for earth - based detectors consists of the intermediate mass ratio inspiral of a compact object ( neutron star or bh ) into an imbh .
another promising source for advanced earth - based interferometers ( albeit with highly uncertain event rates ) are imbh - imbh mergers .
the initial inspiral of these binaries could be detected via space - based interferometers , while the ringdown phase is in the optimal bandwidth for second- and third - generation detectors such as advanced ligo and et , that could therefore be used for `` follow - up '' ringdown searches @xcite .
as in the case of bh spins , the observations imply that numerical relativity needs to cover a wide range in @xmath85 .
this can be done most efficiently by bridging the gap between numerical studies and the perturbative modeling of extreme - mass - ratio binaries ( e.g.  @xcite ) .    from this discussion , it is clear that a detailed knowledge of the spin evolution as well as the generation of gravitational recoil in bh binary mergers is important for understanding the cosmological evolution of smbhs over cosmological times .
we will discuss these effects in turn .
* _ black - hole spins resulting from mergers : _ * prior to the 2005 numerical relativity breakthroughs , it was known that `` minor mergers '' ( @xmath158 ) of a large , rotating bh with an isotropic distribution of small objects would tend to spin down the hole @xcite .
numerical merger simulations showed that the merger of comparable - mass , nonspinning bhs leads to a final kerr bh with spin parameter @xmath159 .
the simulations were followed by the development of several models to predict the spin of merger remnants as a function of the binary parameters for generic mass ratios and spins .
the first studies focussed on performing numerical evolutions of the last few orbits of bh binaries , and used the results to calibrate formulae that map the initial binary parameters to values for the spin of the merged hole @xcite .
it became clear , however , that the binary inspiral up to the last orbits has the potential to significantly affect the spin distribution ( e.g.  @xcite ) and therefore should be included , for example via pn modeling , in the derivation of maps from initial parameters to the merger remnant properties @xcite .
predictions for the final spin based on the extrapolation of test - particle calculations to finite mass ratios have been developed in @xcite , and show remarkably good agreement with numerical results in the comparable - mass regime .
a comprehensive review of all spin formulae is beyond the scope of this article , but we refer the reader to the review in @xcite and the discussion in sec .  v of @xcite .
if we assume an ensemble of bh binaries with initially randomly oriented spins , the final spin distribution is peaked around @xmath160 @xcite .
@xcite reported spin flips by @xmath161 with respect to the initial spin direction of the larger hole ; spin flips of this magnitude could provide an explanation for x - shaped radio sources @xcite .
_ @xcite demonstrated that spin precession over the course of a long inspiral of thousands of orbits tends to align ( antialign ) the binary bh spins with each other if the spin of the more massive bh is initially partially aligned ( antialigned ) with the orbital angular momentum , thus increasing ( decreasing ) the average final spin .
spin precession is stronger for comparable masses , and it could produce significant spin alignment before merger for both smbhs and stellar - mass bh binaries . * _ gravitational recoil : _ * one of the most spectacular results obtained from numerical bh binary simulations is the quantitative prediction of the magnitude of _ gravitational recoils _ ( or _ kicks _ ) . in the 1960s
it was realized that the emission of linear momentum via gws must impart a kick on the source due to momentum conservation @xcite . in the 1980s ,
the kick generated in compact binary inspirals and plunges was studied in the framework of pn theory and bh perturbation theory @xcite .
however , the astrophysical relevance of the gravitational recoil following bh binary mergers remained an open question until the recent numerical relativity breakthroughs .    for the case of an equal - mass , nonspinning bh binary , the net linear momentum emitted in gws vanishes due to symmetry .
nonzero recoils are therefore only generated in systems where this symmetry is broken through ( i ) a mass ratio @xmath162 ( equivalently , a symmetric mass - ratio parameter @xmath163 ) or ( ii ) nonvanishing spins .
for zero spins , early numerical studies in the range @xmath164 $ ] found that the kick velocity is well approximated by @xcite @xmath165 this result ( whose functional form is inspired by fitchett s analytical work @xcite ) was confirmed by simulations for smaller mass ratios ( @xmath166 , @xcite ) and by analytical work @xcite .
the maximal recoil obtained from eq .
( [ eq : nonspinningkick ] ) is @xmath167 for @xmath168 .
small orbital eccentricity @xmath169 should introduce corrections proportional to @xmath170 @xcite .
spinning binaries are characterized by seven free parameters ( the mass ratio plus three components for each bh spin ) , and numerical studies inevitably focussed on subsets of the parameter space .
it soon became clear that the spin interaction dominates over mass ratio effects .
the first studies for equal - mass binaries with spins parallel to the orbital angular momentum revealed kicks of several hundreds km / s , with an extrapolated maximum of @xmath171 for extremal spin magnitudes @xcite .
shortly thereafter , the discovery of the so - called _ superkicks _ marked one of the most surprising outcomes of numerical relativity : binaries for which the spins are perpendicular to the orbital angular momentum and anti - aligned with each other can generate kicks of _ thousands _ of km / s , with an extrapolated maximum @xmath172 for near - extremal spin magnitudes @xcite .
most recently , bh binaries with spins partially aligned with the orbital angular momentum but whose spin projections into the orbital plane still correspond to the superkick configurations have been found to result in even larger maximum recoils of up to @xmath173  km / s @xcite , so - called `` hang - up kicks '' .
many galaxies harbor smbhs at their centers , and galaxy mergers should generally lead to the merger of their central bhs @xcite .
typical escape velocities range from @xmath174 for dwarf galaxies up to @xmath175 for giant elliptic galaxies @xcite .
large kicks would displace or eject the merged hole from its host , with possibly observable consequences : a softening of the stellar density gradient in the galactic nucleus , off - center radio - loud active galactic nuclei , off - nuclear x - ray sources in nearby galaxies and the generation of electromagnetic signals via interaction of the bh with its gaseous environment @xcite .
bh ejection represents a potential obstacle for bh growth via merger , and thus puts constraints on merger - history models , which must be able to explain the assembly of smbhs by redshifts @xmath176 @xcite .
observed redshifts of broad - line relative to narrow - line regions in quasar spectra may be interpreted as a smoking gun of bh ejection due to gravitational recoil @xcite , but alternative interpretations ( such as a bh binary , or superposed emission regions from two interacting galaxies ) are possible @xcite .
the interest in astrophysical consequences of large recoils led to numerical studies of the bh parameter space @xcite .
phenomenological kick formulas inspired by pn studies @xcite and similar to eq .
( [ eq : nonspinningkick ] ) were proposed to map the input parameters of a given binary configuration to the final kick velocity .
available numerical results span mass ratios in the range @xmath177 and spin magnitudes @xmath178 , and are well described by eq .
( 2 ) of ref .  @xcite ( see also @xcite ) .
the apparent incompatibility between large recoils and the existence of bhs at galactic centers can be resolved by mechanisms that would _
align _ the individual bh spins with the orbital angular momentum of the binary .
one such mechanism are gas torques in the so - called `` wet '' ( gas - rich ) mergers , that would produce partial alignment `` early on '' in the inspiral phase @xcite .
if partial alignment in gas - rich mergers is the norm , as suggested also by the spin measurements discussed above , pn spin effects will lead to _ further _ alignment of the spins with the orbital angular momentum , significantly reducing the typical values of recoil velocities @xcite .
this mechanism has even been shown capable of efficient suppression of large recoil velocities in the case of the above mentioned hang - up kicks @xcite .    * _ merger simulations with matter and the blandford - znajek effect : _ * many numerical relativity groups are presently working on binary bh simulations in the presence of matter .
most of this work is trying to understand the signature of electromagnetic counterparts to binary bh mergers : for example , such a signature could be produced by merger events @xcite or recoiling bhs @xcite `` shocking '' the surrounding accretion disks .
the georgia tech group presented the first simulations in full gr of equal - mass , spinning bhs merging in a gas cloud .
they found that shocks , accretion and relativistic beaming can produce electromagnetic signatures correlated with gws , especially when spins are aligned with the orbital axis @xcite .
later work by the group focussed on hot , radiatively inefficient accretion flows . in this case , bh binaries exhibit a flare followed by a sudden drop in luminosity associated with the merger , and quasi - periodic oscillations correlated with the gws during inspiral @xcite .
the urbana group simulated equal - mass , nonspinning bh binaries embedded in gas clouds under different assumptions for the motion of the binary with respect to the gas .
they found evidence that the accretion rate and luminosity due to bremsstrahlung and synchrotron emission would be enhanced with respect to a single bh of the same mass as the binary , possibly being detectable by lsst @xcite .
recent work suggest that the circumbinary disk surrounding bh binaries should not produce detectable electromagnetic counterparts @xcite
. however , the _ magnetic fields _ produced by the circumbinary disk could affect the binary dynamics @xcite , produce stronger electromagnetic counterparts @xcite and extract energy from the orbiting bhs , which ultimately merge within the standard blandford - znajek scenario , generating electromagnetic emission along dual or single jets that could be observable to large distance @xcite ( but see also @xcite ) .
for further reading , we recommend the following reviews .
an excellent summary of bh mass and spin measurements and ( possible ) evidence for event horizons based on `` traditional '' electromagnetic astronomy is given by narayan @xcite .
the use of electromagnetic observations of bhs and neutron stars to probe strong - field gravity is reviewed by psaltis @xcite .
a more extended summary of numerical results on gravitational recoil can be found in zlochower _
_ @xcite . for a more general review of numerical bh simulations we recommend centrella _
the breakthroughs in numerical relativity have opened the door to tackle many fundamental problems in bh physics and to address questions of wider interest .
a new golden era in bh physics is just starting .
we summarize below what we think have been the main developments in this relatively short period of time .
* _ mathematical physics and fundamental issues . _ * a decades - old problem in gr concerns the high - energy collision of two bhs . due to the dominance of the gravitational interaction at large energies
, this process is thought to describe general trans - planckian scattering of particles : at very large center of mass ( cm ) energies all interactions are `` frozen '' while gravity is boosted , thus all the details about internal structure of the colliding particles should be washed out .
this is further supported by thorne s hoop conjecture ( recently tested in highly dynamical situations @xcite ) , which predicts bh formation from generic high - energy collisions of particles @xcite .
because the final bh horizon cloaks all details about the structure of the colliding objects , _ matter does not matter _ and one can for simplicity study bhs as representing a wide class of colliding objects @xcite .    for large cm energies ,
bh collisions are arguably the most violent and nonlinear process one could conceive of . evolving einstein s equations in such extreme regimes
poses new problems , such as the need to deal with all the different scales involved , and raises several questions . a fundamental issue concerns cosmic censorship .
in four spacetime dimensions , bhs have an upper bound on their angular momentum , given by , while that of a spinning top @xmath179 , so ( even though the event horizon generators move at the speed of light ) bhs rotate `` slowly '' , as measured by their kerr parameters . ] @xmath180 is it possible to tune the impact parameter in such a way as to produce a final object spinning above the kerr bound , or does nature somehow conspire to always radiate enough angular momentum as to produce a bh ?
this problem has been addressed in recent years @xcite .
> b_{\rm scat}$ ] ) , a prompt merger ( @xmath181<b^*$ ] ) and a nonprompt merger ( @xmath182<b_{\rm scat}$ ] ) . taken from ref .
, width=151 ]    we show the trajectory of two equal - mass , nonspinning , colliding bhs in fig .
[ fig : crit_traj ] , for cm velocities of around @xmath183 .
we notice three distinct regimes . for impact parameters above the _ scattering threshold _
( @xmath184 , the gravitational interaction between the two colliding bhs is relatively weak and they scatter to infinity . for impact parameters below the _ threshold of immediate merger _ ( @xmath185 ) ,
the bhs promptly collide and merge into a single bh .
the transition between these two states for @xmath186 is a nonprompt merger , which puts the bhs in a so - called zoom - whirl orbit , i.e.  the number of orbits @xmath187 .
the results show that delayed mergers are necessary in order for the system to radiate excess angular momentum , so that the end product obeys ( [ kerrbound ] ) .
the radiated energy and angular momentum increase strongly as the impact parameter approaches the threshold of immediate merger , as can be seen in fig .
[ fig : e_b ] .
for _ head - on _ collisions , numerical results suggest that a fraction @xmath188 of the cm energy is radiated in the limit of very large cm energy @xcite . for finely tuned impact parameters , the
total radiated energy ( in gws ) can be of order @xmath189 for a collision with @xmath190 in the cm frame , while the spin of the final hole can reach @xmath191 of the bound ( [ kerrbound ] ) .
an open issue is how much cm energy can be radiated for larger boosts .
finally , peak luminosities of around @xmath192 are attained in these simulations .
this number is orders of magnitude above any electromagnetic phenomena and close to the maximum conjectured possible value @xcite .
the scattering threshold as a function of boost was studied by shibata _
et al . _ and sperhake _
_ @xcite . for the range of boosts studied it is well approximated by @xmath193 here @xmath78 is the total adm mass and @xmath194 the velocity of each hole as measured in the cm frame .    * _ higher - dimensional black holes . _ * although bh solutions and their properties should in principle fall under the `` mathematical physics '' category , the activity in this specific field has been so intense that it deserves a subsection of its own .
asymptotically flat higher - dimensional black objects have a much richer structure than their four - dimensional counterparts .
for instance , spherical topology is not the only allowed topology for objects with a horizon .
one can also have , _
e.g. _ , black rings , with a donut - like topology .
remarkably , these two different horizon topologies coexist for certain regions in phase - space @xcite .
explorations of the stability of general higher - dimensional bhs are in their infancy .
generically it has been conjectured that , for @xmath195 , ultra - spinning myers - perry bhs will be unstable @xcite .
this instability has been confirmed by an analysis of linearized axisymmetric perturbations in @xmath196 @xcite .
clearly , the study of the nonlinear development of these instabilities requires numerical methods , such as the ones reviewed here .
a study of this type was very recently presented for nonaxisymmetric perturbations in @xmath197 @xcite , where it was found that a single spinning five dimensional myers - perry bh is unstable , for sufficiently large rotation parameter ( confirming previous conjectures @xcite ) .
general equilibrium states in anti - de sitter backgrounds only recently started to be explored : nonaxisymmetric solutions have finally been built @xcite , confirming previous conjectures about their existence @xcite , and braneworld bhs are now also starting to be explored @xcite .
fully dynamical situations are still uncharted territory .    * _ tev - scale gravity scenarios .
_ * the above studies are of direct relevance to some high - energy scenarios .
an outstanding problem in high - energy physics is the extremely large ratio between the four dimensional planck scale , @xmath198 gev , and the electroweak scale , @xmath199 gev .
it has been proposed that this _ hierarchy problem _ can be resolved by confining the standard model to a brane in a higher dimensional space @xcite .    in such models ,
the fundamental planck scale  the energy at which gravitational interactions become strong  could be as low as 1 tev .
thus , high - energy colliders , such as the large hadron collider ( lhc ) , may directly probe strongly coupled gravitational physics @xcite .
in fact , such tests may even be routinely available in the collisions of ultra high - energy cosmic rays with the earth s atmosphere @xcite , or in astrophysical bh environments @xcite ( for reviews see @xcite ) .
the production of bhs at trans - planckian collision energies ( compared to the fundamental planck scale ) should be well described by using classical gr extended to @xmath9  dimensions @xcite .
the challenge is then to use the classical framework to determine the cross section for production and , for each initial setup , the fractions of the collision energy and angular momentum that are lost in the higher - dimensional space by gw emission .
this information will be of paramount importance to improve the modelling of microscopic bh production in event generators .
the first models for bh production in parton - parton collisions used a simple black disk approach to estimate the cross section for production @xcite .
as we already described , only recently exact results for highly relativistic collisions where obtained in four dimensions , using numerical relativity techniques @xcite .
the extension to higher - dimensional spacetimes is a topic of current investigations @xcite .
the first numerical results concerning low - energy collisions have been reported last year @xcite , initial data for boosted bh binaries have been constructed using an extended puncture approach in @xcite and some results for high - energy collisions have been presented in @xcite .    * _ ads / cft and holography . _ * in 199798 , a powerful new technique known as the ads / cft correspondence or , more generally , the gauge - gravity duality , was introduced and rapidly developed @xcite .
this holographic correspondence provides an effective description of a nonperturbative , strongly coupled regime of certain gauge theories in terms of higher - dimensional classical gravity in ads backgrounds .
in particular , equilibrium and nonequilibrium properties of strongly coupled thermal gauge theories are related to the physics of higher - dimensional bhs , black branes and their fluctuations .
these studies revealed intriguing connections between the dynamics of bh horizons and hydrodynamics , and offer new perspectives on notoriously difficult problems , such as the bh information loss paradox , the nature of bh singularities and quantum gravity .
numerical relativity in anti - de sitter backgrounds is bound to contribute enormously to our understanding of the gauge - gravity duality and is likely to have important applications in the interpretation of observations @xcite .
for instance , in the context of the gauge - gravity duality , high - energy collisions of bhs have a dual description in terms of _ a ) _  high - energy collisions with balls of deconfined plasma surrounded by a confining phase , and _
b ) _  the rapid localized heating of a deconfined plasma .
these are the type of events that may have direct observational consequences for the experiments at brookhaven s relativistic heavy ion collider ( rhic ) @xcite .
numerical relativity in anti - de sitter is particularly difficult , and so far only very special situations have been handled @xcite . among these
, we note the recent numerical proof that pure ads space is unstable against collapse to bhs @xcite .
bhs in the context of high - energy physics have also been discussed by pretorius @xcite and in the white paper @xcite .
a recent review of numerical simulations of black strings in five spacetime dimensions is given by lehner and pretorius @xcite .
an overview about numerical results on bhs in @xmath200 dimensions is given in @xcite .
finally we note kanti s review on bhs at the lhc @xcite .
we conclude this review with a brief summary of the most urgent problems to be tackled by numerical relativity in the three main areas in the near future .
* _ gravitational wave physics : _ * in order to meet tight accuracy requirements in gw template generation , especially for parameter estimation , we either need longer numerical waveforms or a denser spacing of templates in the parameter space . both approaches will require computational resources that grow non - linearly as a function of the accuracy , either because of the slow nature of the inspiral at larger bh separations ( see sec .  2 in @xcite for a quantitative discussion ) or because of the dimensionality of the parameter space .
a second major task is the extension of existing template models to generic binaries with precessing spins and/or smaller mass ratios , bridging the gap to perturbative modeling of extreme mass ratio binaries .    * _ astrophysics : _ * while existing formulae for the kick and final spin resulting from the coalescence of bhs have been helpful for astrophysical studies , the calibration of generic models as proposed by boyle _
@xcite requires a more comprehensive exploration of the parameter space .
the relatively young field of numerical relativity simulations of bhs surrounded by accretion disks will likely improve our insight into expected optical counterparts to bh binary mergers , and thus provide a vital tool for the area of _ multi - messanger astrophysics_.    * _ high - energy physics : _ * arguably the most urgent task is the determination of the scattering cross - section and the loss of energy and momentum in gws in trans - planckian scattering of bhs in _ arbitrary _ dimensions . in view of the stability problems encountered in present studies
, it also appears desirable to obtain a better understanding of the well - posedness of the formulations currently employed for such studies .
second , the modeling of bhs in generic spacetimes such as de sitter and anti - de sitter is still in its infancy , even if preliminary studies of bhs in non - asymptotically flat spacetimes are encouraging . a better understanding of boundary issues
, well - posedness of the formulations and diagnostics such as wave extraction tools and horizon finders will be required to open up this uncharted and fertile ground of research .
* acknowledgements : * u.s .
acknowledges support from the ramn y cajal programme and grant no .
fis2011 - 30145-c03 - 03 of the ministry of education and science of spain , the fp7-people-2011-cig grant _ cbheo 293412 _ , nsf grants phy-0601459 , phy-0652995 and the sherman fairchild foundation to caltech .
e.b.s research is supported by nsf grant phy-0900735 and by nsf career grant phy-1055103 .
this work was supported by the dybho-256667 erc starting grant , and by fct - portugal through projects ptdc / fis/098025/2008 , ptdc / fis/098032/2008 , cte - ast/098034/2008 , cern / fp/116341/ 2010 and cern / fp/123593/2011 .
this research was supported by the dybho256667 erc starting grant , the nrhep295189 fp7people2011irses grant by allocations through the teragrid advanced support program at the san diego supercomputing center and the national institute for computational sciences under grant phy-090003 , the centro de supercomputacin de galicia ( cesga ) under project numbers icts - cesga-175 , icts - cesga-200 and icts - cesga-221 , the barcelona supercomputing center ( bsc ) under project aect-2011 - 2 - 0006 , aect-2011 - 3 - 0007 , aect-2012 - 1 - 0008 , the deisa extreme computing initiatives deci-6 and deci-7 prace-2ip project dybho .