Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the final merger of two black holes in a binary system releases more power than the combined light from all the stars in the visible universe . this vast energy comes in the form of gravitational waves , which travel across the universe at the speed of light , bearing the waveform signature of the merger . today , ground - based gravitational - wave detectors stand poised to detect the mergers of stellar black - hole binaries , the corpses of massive stars . in addition , planning is underway for a space - based detector that will observe the mergers of massive black holes , awesome behemoths at the centers of galaxies , with masses of @xmath0 , where @xmath1 is the mass of the sun . since template matching forms the basis of most gravitational - wave data analysis , knowledge of the merger waveforms is crucial . calculating these waveforms requires solving the full einstein equations of general relativity on a computer in three spatial dimensions plus time . as you might imagine , this is a formidable task . in fact , numerical relativists have attempted to solve this problem for many years , only to encounter a host of puzzling instabilities causing the computer codes to crash before they could compute any sizable portion of a binary orbit . remarkably , in the past few years a series of dramatic breakthroughs has occurred in numerical relativity ( nr ) , yielding robust and accurate simulations of black - hole mergers for the first time .    in this article , we review these breakthroughs and the wealth of new knowledge about black - hole mergers that is emerging , highlighting key applications to astrophysics and gravitational - wave data analysis . we focus on comparable - mass black - hole binaries , with component mass ratios @xmath2 , where @xmath3 and @xmath4 , @xmath5 are the individual black - hole masses . we will frequently also refer to the _ symmetric mass ratio _ @xmath6 for simplicity , we choose to set @xmath7 and @xmath8 ; with this , we can scale the dynamics and waveforms for black - hole binaries with the total system mass @xmath9 . in particular , we can express both length and time scales in terms of the mass , giving @xmath10 . we begin by setting both the scientific and historical contexts . in sec . [ sec : astro_grav ] we provide a brief overview of astrophysical black - hole binaries as sources for gravitational - wave detectors . we next turn to a historical overview in sec . [ sec : history ] , surveying efforts to evolve black - hole mergers on computers , spanning more than four decades and culminating with the recent triumphs . having thus set the stage , we focus on more in - depth discussions of the key components underlying successful black - hole merger simulations , discussing computational methodologies in sec . [ sec : comp ] , including numerical - relativity techniques and black - hole binary initial data . section  [ sec : gws ] is the heart of this review . here we discuss the key results from numerical relativity simulations of black - hole mergers , following a historical development and concentrating on the merger dynamics and the resulting gravitational waveforms . these results have opened up a variety of exciting applications in general relativity , gravitational waves , and astrophysics . we discuss synergistic interactions between numerical relativity and analytic approaches to modeling gravitational dynamics and waveforms in sec . [ sec : nrpn ] , and applications of the results to gravitational - wave data analysis in sec . [ sec : gwda ] . the impact of merger simulations on astrophysics is presented in sec . [ sec : astro ] , which includes discussions of recoiling black holes and potential electromagnetic signatures of the final merger . we conclude with a look at the frontiers and future directions of this field in sec . [ sec : frontiers ] . before we begin , we mention several other resources that may interest our readers . the review articles by @xcite and @xcite provide interesting surveys of numerical relativity several years before the breakthroughs in black - hole merger simulations . the article by @xcite is an early review of the recent successes , covering some of the same topics that we discuss here . @xcite reviewed the status of black - hole simulations producing long waveforms ( including at least ten cycles of the dominant gravitational - wave mode ) and their application to gravitational - wave data analysis . finally , the text books by @xcite and @xcite provide many more details on the mathematical and computational aspects of numerical relativity than we can include here , and serve as useful supplements to our discussions . black holes and gravitational waves are surely among the most exotic and amazing predictions in all of physics . these two offspring of einstein s general relativity are brought together in black - hole binaries , expected to be among the strongest emitters of gravitational radiation . we begin by presenting some basic properties of black holes and gravitational waves . for fuller discussions and more details , see @xcite and @xcite . a black hole forms when matter collapses to infinite density , producing a singularity of infinite curvature in the fabric of spacetime . each black hole is surrounded by an event horizon , at which the escape velocity is the speed of light . the event horizon is a global property of the spacetime , since it is defined by the paths of  outgoing " photons that are the boundary between photon trajectories that must fall inward , and those that can escape to infinity . the photons defining the event horizon hover at finite radius at the surface of the black hole . since , in principle , mass ( energy ) can fall into the event horizon at late times  which will move the location at which photon paths can hover  we must know the entire future development of the system to locate the event horizon .    when black holes merge , a single event horizon forms whose area is at least as large as the sum of the individual horizons . since numerical relativists want to know when this occurs during the course of a calculation , they rely on a related concept known as an _ apparent horizon _ , whose location depends only on the properties of the spacetime at any given time @xcite . for quiescent black holes , the apparent and event horizons coincide ; for more general holes , the apparent horizon is always inside the event horizon ( with restrictions on the behavior of the matter involved ) . so , in terms of causality in a numerically - generated spacetime , any physical phenomenon found inside an apparent horizon should not leak out and affect the spacetime outside . the simplest black hole is nonrotating and is described by the spherically symmetric schwarzschild solution to the einstein equations of general relativity in vacuum ( i.e. , with no `` matter '' sources in the spacetime ) . a schwarzschild black hole is fully specified by one quantity , its mass @xmath9 . the horizon is located at coordinate @xmath11 ( in schwarzschild coordinates ) ; its area is @xmath12 . more general black holes can have both charge and spin . since a charged black hole in astrophysics will generally be neutralized rapidly by any surrounding plasma , we can consider only rotating , uncharged black holes . stationary ( i.e. , time independent ) black holes are described by the axisymmetric kerr solution . a kerr black hole is fully specified by two quantities , its mass @xmath9 and its angular momentum per unit mass @xmath13 . the event horizon is located at the boyer - lindquist @xcite radius @xmath14 , where @xmath15 the area of the event horizon is @xmath16 . equation   requires @xmath17 ; when @xmath18 the black hole is said to be maximally rotating or `` extremal '' . notice that @xmath19 when @xmath20 , and that @xmath14 decreases as @xmath13 increases , thus bringing the location of the horizon deeper into the potential well as the black - hole spin increases . photons and test particles in the vicinity of a single black hole can experience either stable or unstable orbits . for a schwarzschild black hole of mass @xmath9 , the _ innermost stable circular orbit _ ( isco ) occurs at @xmath21 for massive test particles . in the case of a kerr black hole , the isco is closer in for co - rotating test particles , and farther out for counter - rotating particles . while the concept of an isco is strictly defined only for massive test particles , it has proven useful for studies of the spacetime around two black holes spiralling together on quasicircular orbits . imagine that you put the two black holes on an _ instantaneously _ circular orbit around each other ; at that moment they have neither nonzero radial velocity nor nonzero radial acceleration . at any given separation , the black holes have some angular momentum . the isco is the separation where that angular momentum is a minimum , in analogy to the test particle definition . black holes at closer separations would be expected to fall inward , toward the center , even without radiating angular momentum via gravitational radiation . gravitational waves are ripples in the curvature of spacetime itself . they carry energy and momentum and travel at the speed of light , bearing the message of disturbances in the gravitational field .    as with electromagnetic waves , gravitational waves can be decomposed into multipolar contributions that reflect the nature of the source that generates them . recall that electromagnetic radiation has no monopole contribution due to the conservation of total charge . by analogy , conservation of total mass - energy guarantees that there can be no monopole gravitational radiation . since dipolar variations of charge and currents are possible , electromagnetic waves can have a dipole character . however , conservation of linear and angular momenta removes any possibility of dipolar gravitational waves , so the leading - order contribution to gravitational radiation is quadrupolar . gravitational waves are thus generated by systems with time - varying mass quadrupole moments @xcite . in the wave zone , a gravitational wave is described as a perturbation @xmath22 to a smooth underlying spacetime . the wave amplitude is @xmath23 where @xmath24 is the quadrupole moment of the source , @xmath25 is the distance from the source , and @xmath26 is the mass in the source that is undergoing quadrupolar changes . this shows that the strongest gravitational waves will be produced by large masses moving at high velocities , such as binaries of compact stars and black holes . the wave on the left is purely in the @xmath27 polarization state , and the one on the right is purely in the @xmath28 polarization state . the gravitational waves produce tidal forces in planes transverse to the propagation direction . from @xcite . reprinted with permission from aaas.,width=336 ]    a gravitational wave is purely transverse , acting tidally in directions perpendicular to its propagation direction . when a gravitational wave impinges on a detector of length scale @xmath29 , it produces a length change in that detector @xmath30 . by substituting in typical values for compact objects in the universe into eq .  , one can see that astrophysical sources typically yield wave amplitudes of @xmath31 at the earth . consequently , precision measurements are needed to make detections . gravitational waves have two polarization components , known as @xmath32 and @xmath33 for a linearly polarized wave . figure  [ fig : forcelines ] shows the corresponding lines of force for sinusoidal gravitational waves propagating along the @xmath34 axis . in the purely @xmath27 polarization on the left , the wave stretches along one axis and squeezes along the other , alternating sinusoidally as the wave passes . the @xmath28 polarization wave on the right acts similarly , stretching and squeezing along axes rotated by @xmath35 . in general , a gravitational wave is a superposition of these two states , conveniently written as a complex waveform strain @xmath22 , where @xmath36      astronomers have found evidence for black holes throughout the universe on a remarkable range of scales . the smallest of these , stellar black holes , have masses in the range @xmath37 and form as the end - products of massive star evolution . there is good observational evidence for the existence of stellar black holes , based on dynamical measurements of the masses of compact objects in transient systems that undergo x - ray outbursts . since neutron stars can not have masses @xmath38 , compact objects more massive than this must be black holes @xcite . intermediate - mass black holes ( imbhs ) have masses in the range @xmath39 . imbhs may form as the result of multiple mergers of smaller objects in the centers of dense stellar clusters in the present universe @xcite , assuming mass loss from stellar winds is not significant @xcite . they may also arise from the evolution of very massive stars early in the history of the universe , forming black - hole `` seeds '' in the centers of massive halos ( the precursors of the galaxies we see today ) early in the history of the universe , to redshifts @xmath40 @xcite . currently the best observational evidence for imbhs comes from models of ultraluminous x - ray sources @xcite . finally , massive black holes ( mbhs ) have masses in the range @xmath41 and are found at the centers of galaxies , including our own milky way galaxy . the observational case for the existence of mbhs is quite strong , based on dynamical models of stars and gas believed to be moving in the potential well of the central mbh @xcite . black - hole binaries are binary systems in which each component is a black hole . as mentioned above , we focus here on comparable - mass binaries , which are expected to produce the strongest gravitational - wave signals . stellar black - hole binaries may form as the result of binaries composed of two massive stars ; see @xcite and references therein . stellar black - hole binaries may also arise from dynamical processes in which a black hole is captured into an orbit around another black hole in dense stellar environments @xcite . imbh binaries can also form through dynamical processes in stellar clusters @xcite and from mergers of massive halos at high redshifts . since both stellar and imbh binaries are `` dark ''  that is , they are generally not surrounded by gas which might produce electromagnetic radiation  we have few observational constraints on these types of black - hole binaries . this situation will change dramatically , however , with the detection of gravitational radiation from these systems @xcite , as gravitational waves bring direct information about the dynamical behavior of the orbiting masses and do not rely on electromagnetic emissions from nearby matter .    since essentially all galaxies are believed to contain an mbh at the center and to undergo a merger with another galaxy at least once during the history of the universe , mbh binaries can arise when their host galaxies merge @xcite ; see also @xcite and references therein . however , due to the vast cosmic distances involved , and the small angular separations on the sky expected for mbh binaries , only a few candidates are currently known through electromagnetic observations @xcite . when they form , mbh binaries typically have relatively wide separations , and the gravitational radiation they emit is very weak and insufficient to cause the binary to coalesce within the age of the universe . however , various processes such as gaseous dissipation and @xmath42-body interactions with stars can remove orbital energy from the binary and cause the black holes to spiral together @xcite ; see also @xcite and @xcite and references therein . eventually , the black holes reach separations at which gravitational radiation reaction becomes the dominant energy - loss mechanism , leading to the final coalescence of the black holes and the emission of strong gravitational waves @xcite . mergers of comparable - mass black - hole binaries are expected to be among the strongest sources of gravitational waves . this final death spiral of a black - hole binary encompasses three stages : inspiral , merger , and ringdown @xcite . in the early stages of the inspiral , the orbits of most astrophysical black - hole binaries will circularize due to the emission of gravitational radiation @xcite . during the inspiral , the orbital time scale is much shorter than the time scale on which the orbital parameters change ; consequently , the black holes spiral together on quasicircular orbits . since the black holes have wide separations , they can be treated as point particles . the inspiral dynamics and waveforms can be calculated using post - newtonian ( pn ) equations , which result from a systematic expansion of the full einstein equations in powers of @xmath43 , where @xmath44 is the binary separation @xcite . the inspiral phase produces gravitational waves in the characteristic form of a _ chirp _ , which is a sinusoid with both frequency and amplitude increasing with time . as the black holes spiral inward , they eventually reach the strong - field , dynamical regime of general relativity . in this merger stage , the orbital evolution is no longer quasi - adiabatic ; rather , the black holes plunge together and coalesce into a single , highly distorted remnant black hole , surrounded by a common horizon . since the point - particle and pn approximations break down , numerical relativity simulations of the einstein equations in three dimensions are needed to calculate the merger . due to the difficulty of these simulations , the resulting gravitational waveforms were completely unknown until recently . finally , the highly distorted remnant black hole settles down into a quiescent rotating kerr black hole by shedding its nonaxisymmetric modes through gravitational - wave emission . we call this process the `` ringdown , '' in analogy to how a bell that has been struck sheds its distortions as sound waves . various analytic techniques of black - hole perturbation theory @xcite form the basis of ringdown calculations , producing gravitational waveforms in the shape of exponentially damped sinusoids @xcite . the characteristic gravitational - wave frequency of a quasicircular black - hole binary , produced by the dominant ( highest - order ) quadrupole component , is @xmath45 where @xmath46 is the orbital frequency . astrophysical black - hole binaries produce gravitational waves that span three frequency regimes , depending on the black - hole masses @xcite . stellar black - hole binaries and the lower mass end of the imbh binaries radiate in the high frequency band , @xmath47 , which is already being observed by ground - based laser interferometer detectors such as ligo @xcite , and will be observed by the advanced detectors by @xmath48 @xcite . low frequency gravitational waves cover the band @xmath49 and will be observed by the space - based laser interferometer lisa , currently under development @xcite . mbh binaries with masses @xmath50 will be very strong sources for lisa , with the lower mass systems visible out to redshifts @xmath40 @xcite ; the inspirals of imbh binaries will also be detectable @xcite . finally , the very low frequency band @xmath51 will be observed by pulsar timing arrays @xcite . this band is expected to dominated by gravitational waves from a very large population of unresolved mbh binaries @xcite with possibly a few discrete sources @xcite . the quest to calculate the gravitational - wave signals from the merger of two black holes spans more than four decades and encompasses key developments in theoretical and experimental general relativity , astrophysics , and computational science . in this section , we begin by delineating these threads in general terms , and then turn to a more detailed account of select milestones along the path toward successful simulations of black - hole mergers .      at the end of the 18th century @xcite and @xcite first speculated , using newtonian gravity , that a star could become so compact that the escape velocity from its surface would exceed the speed of light . in the 20th century , scientists realized that such black holes could form as the final state of total gravitational collapse in general relativity @xcite ; john wheeler would later popularize the term `` black hole '' to describe such an object @xcite . beginning in the 1960s , many highly energetic astrophysical phenomena were discovered with physical properties pointing to extremely strong gravitational fields as underlying mechanisms ; among these are quasars and x - ray binaries such as cygnus x-1 @xcite , the first credible black - hole candidate . as discussed in sec . [ sec : astro_grav : astro_bhs ] , today astrophysical black holes are believed to exist on a vast range of scales throughout the universe , and black - hole binaries are considered to be strong sources of gravitational waves . einstein s equations of general relativity form a coupled set of nonlinear partial differential equations , in which dynamic curved spacetime takes the role of newton s gravitational field and interacts nonlinearly with massive bodies . gravitational waves were first recognized as solutions to the linearized , weak - field einstein equations early in the past century . by mid - century , gravitational waves were recognized as real physical phenomena , carrying energy and being capable of producing a response when impinging on a detector . this development spawned a major branch of experimental general relativity , with concepts for the first gravitational - wave detectors appearing in the 1960s ; see @xcite for a review . the discovery of two neutron stars in a binary system by @xcite provided an astrophysical laboratory for the first _ indirect _ detection of gravitational radiation . decades of observation revealed the binary orbit to be shrinking by precisely the amount expected if the system were emitting gravitational waves according to general relativity @xcite ; hulse and taylor were awarded the nobel prize in 1993 . today , the prognosis for _ direct _ detection of gravitational waves is excellent , with the first events expected from the advanced ground - based interferometers around the middle of this coming decade . as you might expect , einstein s equations pose formidable obstacles to anyone who would dare to probe the physics within . throughout most of the 20th century , relativists uncovered a fairly small number of exact solutions by exploiting symmetries , and made progress toward more general problems using various perturbative expansions . by the 1960s , computers had become powerful enough to encourage attempts at solving einstein s equations numerically , to uncover physics beyond the realm of perturbation theory . the subsequent development of numerical relativity has been made possible in part by continued increases in computer power and advances in algorithms and computational methods . most numerical - relativity simulations start with the idea of decomposing four - dimensional spacetime into a stack of curved three - dimensional spacelike slices threaded by a congruence of timelike curves @xcite . arnowit , deser , and misner ( adm ) pioneered this `` 3 + 1 '' approach as the basis for a canonical formulation of the dynamics of general relativity @xcite . in this hamiltonian formulation , the three - metric @xmath52 on the spatial slices takes the role of the `` configuration variables . '' quantities based on the extrinsic curvature @xmath53 , which is roughly the time derivative of @xmath52 , play the role of `` conjugate momenta . '' variation of an action with respect to @xmath52 produces a set of six first - order evolution equations for the conjugate momenta ; varying the momenta gives six first - order evolution equations for @xmath52 . adm also introduced four lagrange multipliers as freely specifiable gauge or coordinate conditions , representing the four coordinate degrees of freedom in general relativity . variation of these lagrange multipliers yields four equations that must hold on each slice : a hamiltonian constraint , and three momentum constraints . originally intended as a tool for quantizing gravity , the adm formalism later became the basis for most work in numerical relativity . as we discuss in sec . [ sec : comp ] below , key elements in this approach are solving the cauchy problem , beginning with the initial data on a spacelike slice , and then evolving that data forward in time . the constraint equations form the basis for this initial value problem . appropriate choices for the gauge conditions are crucial ingredients for today s successful black - hole merger simulations . @xcite made the first known attempt to simulate the head - on collision of two equal - mass black holes on a computer in 1964 , using a two - dimensional axisymmetric approach . their simulation ran for 50 timesteps to a duration of @xmath54 ; at this point , they decided that the simulation was no longer accurate enough to warrant continued evolution , and stopped the code . smarr and eppley @xcite returned to this problem in the mid-1970s , again employing two - dimensional axisymmetry but now using the adm formalism , specialized coordinates , and improved coordinate conditions . although they encountered problems with instabilities and large numerical errors , they managed to evolve the collision and extract information about the emitted gravitational waves . smarr and eppley had used the most powerful computers of their day . going to the next step , orbiting black holes in three dimensions , was deemed to be not feasible at the time , due to unresolved questions about the instabilities and insufficient computer power . the black - hole merger problem lay dormant for over a decade .    in the 1990s , attention returned to black - hole mergers as the ligo project began to move forward , and black - hole mergers were recognized as the strongest sources for this detector . since the signal - to - noise ratio for such ground - based detectors is fairly modest , having a template for the merger waveform is a key part of the data analysis strategy . numerical relativists revisited the problem of colliding two black holes head - on with modern techniques and more powerful computers @xcite . in the mid-1990s , the national science foundation funded a computational grand challenge grant for a large multi - institution collaboration aimed at evolving black - hole mergers in three dimensions and calculating the resulting gravitational waveforms . around the same time , a large and very active numerical relativity group arose at the newly - formed albert - einstein institut ( aei ) in potsdam , germany . during the late 1990s and into the early 2000s , two developments on the experimental side further increased the desire for black - hole merger simulations : the ground - based gravitational - wave detectors started taking data , and interest grew in lisa and its potential for observing gravitational waves from massive black - hole binary mergers . while no one expected the task at hand  developing computer codes to solve the full einstein equations in three dimensions for the final few orbits and merger of two black holes  to be simple , numerical relativists found that the problem was far more difficult than anticipated . producing a waveform useful for gravitational - wave detection purposes typically would require running a simulation for a duration of several hundred @xmath9 . however , a variety of instabilities plagued the codes , causing them to crash well before any significant portion of an orbit could be achieved . nevertheless , during a period encompassing a little over a decade , much important work was accomplished that laid the foundations for later success . key milestones include these developments :    * initial data for binary black holes near the isco @xcite ; * new methods for representing black holes on computational grids such as punctures @xcite and excision @xcite ; * recognition of the importance of hyperbolicity in formulating the einstein equations for numerical solution @xcite ; * improved formulations of the einstein equations @xcite ; * fully three - dimensional evolution codes and their use in evolving distorted black holes @xcite , boosted black holes @xcite , head - on collisions @xcite , and grazing collisions @xcite ; * coordinate conditions that keep the slices from crashing into singularities and the spatial coordinates from falling into the black holes as the evolution proceeds @xcite ; * the cactus computational toolkit , which provided a framework for developing numerical - relativity codes and analysis tools used by many groups ; * modern adaptive mesh refinement finite - difference ( including carpet with cactus , bam @xcite , and paramesh @xcite ) , and multi - domain spectral @xcite infrastructures for numerical relativity .    throughout this period , the length of black - hole evolutions gradually increased . improvements in the formalisms allowed simulations of single black holes and , later , two black holes to increase in duration to @xmath55 . the addition of new slicing and shift conditions again increased the evolution times to @xmath56 . the lazarus project took a novel approach to combine these relatively short - duration binary simulations with perturbation techniques for the late - time behavior to produce a hybrid model for a black - hole merger @xcite . by late 2003 , @xcite carried out the first complete orbit of two equal - mass , nonspinning black holes . while this simulation lasted @xmath57 , the code crashed shortly after the orbit was completed and the gravitational waves were not extracted . overall , progress was slow , difficult , and incremental . however , the situation was about to change dramatically .      in early 2005 , frans pretorius electrified the relativity community when he achieved the first evolution of an equal - mass black - hole binary through its final orbit , merger , and ringdown using techniques very different from those employed by the rest of the community @xcite . later in 2005 the groups at the university of texas at brownsville ( utb ) and nasa s goddard space flight center independently and simultaneously discovered a new method , called `` moving punctures , '' that also produced successful black - hole mergers @xcite . their presentations were given back - to - back at a workshop on numerical relativity , to the amazement of each other and the assembled participants .    since the moving - puncture approach was based on underlying techniques used by several other groups , it was rapidly and readily adopted by most in the community , producing a growing avalanche of results . 2006 was a year of many firsts : the first simulations of unequal - mass black holes and the accompanying recoil of the remnant hole @xcite , the first mergers of spinning black holes @xcite , the first long waveforms ( @xmath58 orbits ) @xcite , the first comparisons with pn results @xcite , and the first systematic parameter study in numerical relativity @xcite . the year 2007 opened with the discovery of `` superkicks ''  recoil velocities exceeding @xmath59 @xcite , and in 2008 a black - hole binary merger with mass ratio @xmath60 was accomplished @xcite , while @xcite carried out the first long - term evolution of generic spinning and precessing black - hole binaries . as this article was being written in 2009 the state of the art continues to advance , with the first simulations of black - hole mergers using spectral numerical techniques @xcite , and the first steps towards modeling the flows of gas around the merging black holes @xcite . applications of the merger results in areas such as comparisons with pn expressions for the waveforms , astrophysical computations of black - hole merger rates , and the development of templates for gravitational - wave data analysis have accompanied these technical developments in the simulations . the study of black - hole mergers using numerical relativity is thriving indeed . in this section , we discuss the mathematical and numerical foundations underlying current black - hole merger simulations , highlighting the key issues involved in achieving successful evolutions . for more detailed treatments , we direct the interested reader to @xcite or @xcite . the central task of numerical relativity is solving einstein s field equations @xmath61 where einstein s tensor @xmath62 represents the curvature of spacetime , the energy - momentum tensor @xmath63 contains the matter sources , and @xmath64 . by convention , an index @xmath65 selects a `` time '' component , and @xmath66 selects a `` space '' component . @xmath62 depends on the first and second derivatives of the metric tensor @xmath67 . for vacuum black - hole spacetimes , @xmath68 . note that all of the tensor fields discussed here are symmetric in the indices , e.g. @xmath69 . the metric characterizes the geometry of spacetime by giving the infinitesimal spacetime interval @xmath70 through the following definition : @xmath71 where we use the einstein summation convention which implies summation over all values of a given index that appears twice in an expression . many physical implications of the metric are immediately apparent from eq . ( [ eq : interval ] ) . for example , when @xmath72 , the resulting metric - determined relationship between the time and space coordinates yields the paths that light rays must follow in this spacetime . the dependence of einstein s tensor on the metric can be simply illustrated in coordinates known as  harmonic " ; in vacuum eq .   takes the form : @xmath73 where @xmath74 is the flat - space wave - operator and @xmath75 represents all terms nonlinear in the metric . if @xmath75 is interpreted as an effective source term , this is a simple wave equation . the familiar form of this equation suggests that its cauchy problem can be solved by specifying the metric and its first derivative on an initial , spatial surface and then integrating in time , as for an ordinary wave equation . einstein s equations admit various formulations and coordinate conditions , which should be tailored to the problem at hand  in this case , numerical simulation of black - hole spacetimes . regardless of these choices , current numerical practices universally involve an initial three - dimensional slice of spacetime that is evolved forward in time . here , we review the history and current most common choices of initial data , black - hole representations , formulation , coordinate conditions , and some details of the numerics . the impetus , and most important outcome , of all these developments is the ability to generate gravitational waveforms from black - hole binary sources over many cycles . the cauchy problem concerns solution of the field equations given initial data specified on an initial ( typically spatial ) hypersurface . in practice , the cauchy problem is more conveniently investigated in a `` 3 + 1 '' formulation , explicitly based on a foliation of the spacetime into three - dimensional spatial slices parametrized by a time coordinate . a common 3 + 1 formulation inspired by adm @xcite divides up the components of the metric according to their relationships with space and time , such that the line element takes the form : @xmath76 where @xmath77 is called the lapse function , @xmath78 the shift vector , and @xmath79 is the spatial three - metric . we write the time coordinate @xmath80 and the spatial coordinate indices @xmath81 . note that contraction with @xmath67 or its inverse @xmath82 is used to lower or raise indices of four - dimensional tensors , respectively , while @xmath52 or its inverse @xmath83 is used to lower or raise indices of three - dimensional tensors . the lapse and shift represent coordinate freedom in the metric ; we can choose these quantities arbitrarily . however , since the three - metric @xmath52 ( and its first and second spatial derivatives ) determines the intrinsic curvature of the slice , it carries the information about the gravitational field and thus is constrained by the physics . the meaning of the lapse and shift can be understood by considering two successive spatial slices separated by an infinitesimal time interval @xmath84 ( fig . [ fig : smarr_3plus1 ] ) . an observer along a vector normal to the first slice will measure an elapsed proper time of @xmath85 in evolving to the second slice , and a change in spatial coordinate of @xmath86 . since einstein s equations are second order in time , we must also specify the initial time derivative of the three - metric . rather than specifying this derivative directly , we define a new quantity : @xmath87 where @xmath88 is an ordinary partial derivative and @xmath89 is the spatial covariant derivative . note that the spacetime covariant derivative @xmath90 is a partial derivative with a correction such that it transforms as a vector and satisfies @xmath91 . @xmath89 is the projection of @xmath90 onto the spatial slice and is equivalent to a three - dimensional covariant derivative formed from @xmath52 . for the case of euclidean normal coordinates , @xmath92 and @xmath93 , eq .  ( [ eq : kij ] ) reduces to the simple expression @xmath94 . if we define a unit vector @xmath95 normal to the spatial slice , we can show that eq . ( [ eq : kij ] ) is equivalent to @xmath96 . as suggested by this expression , @xmath53 is a measure of the change of the normal vector as it is transported along the slice . in this way @xmath53 gives an extrinsic measure of the curvature of a three - dimensional spatial slice with respect to its embedding in four - dimensional spacetime . it is therefore called the extrinsic curvature . depending on the formulation , the extrinsic curvature might or might not come into the evolution equations ; however , @xmath53 is almost universally utilized when calculating initial data .     and @xmath97 are depicted . @xmath77 is the lapse , and @xmath98 represents the proper time lapse between slices . @xmath99 is the shift , and @xmath100 represents the amount by which the spatial coordinates shift between slices . @xmath101 is normal to the slice at @xmath102 . if a ray parallel to @xmath101 intersects the @xmath102 slice at a point @xmath103 , then it will intersect the @xmath97 slice at @xmath104 . @xmath105 is a coordinate time vector . if a ray parallel to @xmath106 intersects the @xmath102 slice at point @xmath103 , then it will also intersect the @xmath97 slice at @xmath103 . , width=336 ]    of the ten component equations of eq . , six determine the time - evolution of the metric , while four must be satisfied on a spatial slice at any given time , and are thus constraint equations . with an appropriate choice of time coordinate , and assuming vacuum spacetime , these four constraint equations are equivalent to the condition @xmath107 . the time - time component @xmath108 is called the _ hamiltonian constraint _ , and the time - space components @xmath109 are called the _ momentum constraint_. these take the form of conditions on the extrinsic curvature : @xmath110 where @xmath111 is the three - dimensional ricci scalar associated with the three - metric @xmath52 , and @xmath112 , sometimes called the _ mean curvature_. these constraint equations must be solved in order to obtain an initial spatial slice consistent with einstein s equations . how does one represent an exotic object such as a black hole in a numerical simulation ? in particular , how can one use finite computational methods repeatedly to model an object which , analytically , contains physical and/or coordinate singularities ? fortunately , two successful strategies have emerged to meet this challenge . the unusual topology of black holes offers one way out . as @xcite originally showed , a black hole can be considered a  bridge " or  wormhole " connecting one universe or  worldsheet " , to a second worldsheet ( see fig .  [ fig : wormhole_flatspace ] ) . exploiting this topology , a continuous spatial slice which avoids the physical singularity contained within the event horizon of each black hole can be constructed as follows . starting with a spatial slice of schwarzschild spacetime , remove the interior of the event horizon . identify the resulting spherical boundary with the spherical boundary of an identical copy of this space . the two - dimensional analog would be to take two sheets of paper , cut out a disk from each , and then glue the resulting circular edges together . each sheet , or copy of schwarzschild spacetime in this example , is called a worldsheet . the identified spherical boundaries connecting the worldsheets form what is referred to as the  throat " of the wormhole .    to complete this construction , we require an appropriate coordinate system to continuously cover both worldsheets . brill and lindquist discovered coordinates that will prove convenient , in which the three - metric of a particular spatial slice of schwarzschild spacetime is given by @xcite @xmath113 where @xmath25 is a radial coordinate . in these coordinates , the event horizon is at @xmath114 . we can consider each of the worldsheets described above as being separately labeled by such coordinates . designate one worldsheet as @xmath115 and the other as @xmath116 , and call their radial coordinates @xmath117 and @xmath118 respectively . since the interior of the event horizon has been removed from each worldsheet , assume that @xmath119 and @xmath120 . the metric on each worldsheet has the form of eq . ( [ eq : bls ] ) . as noted by brill and lindquist , this form of the metric is unchanged by the transformation @xmath121 , and @xmath122 when @xmath114 . so if we define a new coordinate @xmath123 by @xmath124 mapping spatial infinity on worldsheet @xmath116 to @xmath125 , we obtain a single continuous coordinate system that applies to worldsheet @xmath115 for @xmath126 and worldsheet @xmath116 for @xmath127 . this avoidance of the physical singularity comes at the expense of a coordinate singularity in the metric at @xmath125 , called a `` puncture '' . however , it turns out this coordinate singularity can be confined to a single scalar variable , as suggested in eq . ( [ eq : bls ] ) . with a suitable change of variables and other means , numerical simulations haven proven capable of handling this irregular scalar field . thus the  puncture " method is one way to represent a black hole that is amenable to computation . another strategy is excision , first proposed by unruh @xcite . given that no physical information can escape an event horizon to influence the exterior , the interior of a black hole can in principle be excised from the computational domain . this relies on the fact that all physical information propagates inwards from the event horizon towards the physical singularity , i.e. lightcones tilt inwards , and nothing physical propagates outward from the horizon . extrapolation can be used for the boundary condition of the excised region , and any nonphysical numerical error that escapes the horizon should be negligible . this approach is not constrained to the particular coordinates required by the puncture method , but it can be more difficult due to the need for precise positioning of the excision boundary and accurate extrapolation . excision was used successfully for orbiting binary simulations by @xcite and continues to be used by the caltech - cornell dual - coordinate spectral code @xcite .    either of these methods can be useful in representing black holes on the initial data slice . surprisingly , we will also see that either of these representations can be made robust enough to persist as the black holes evolve . the starting point for successful simulations of black - hole mergers is finding initial data for astrophysically realistic inspiralling black holes . if we were simulating the orbits of stars in newtonian gravity , this would be a simple procedure . for example , we could simply specify the masses and spins of the stars , along with their positions and velocities on orbits derived from the dynamics of point particles , and then evolve the system numerically to allow it to `` relax '' into orbits appropriate for bodies of finite size . in general relativity , however , the initial data must satisfy the constraint equations  ( [ eq : ham ] ) -  ( [ eq : mom ] ) , which are cast in terms of the @xmath128-metric @xmath52 and the extrinsic curvature @xmath53 . since there is no obvious , or natural , connection between these field variables and the astrophysical properties of inspiralling black holes , obtaining suitable initial data is a major challenge . building on the earlier work of @xcite , york developed a general procedure for solving the constraint equations to produce initial data for the cauchy problem in the 1970s ( see @xcite for a review ) . this approach generally requires solving a coupled elliptic system of four nonlinear field equations . we can break this problem down into more manageable pieces with some simplifying assumptions . while these choices do come at the cost of some loss of generality and astrophysical realism in the initial data for two black holes , it has been seen that for sufficiently long evolutions the final orbits and waveform signatures from the black - hole evolution are largely insensitive to this level of detail in the initial data , at least in the case of equal - mass nonspinning black holes . the first simplification is to choose traceless extrinsic curvature , @xmath129 . with this , the hamiltonian constraint is decoupled from the momentum constraints , and can be solved separately . to find solutions corresponding to multiple black holes , we generally further assume that the initial slice is conformally flat . that is , the three - metric is the product of a scalar conformal factor with a flat metric , @xmath130 with this , the problem reduces to solving a ( typically nonlinear ) equation for the scalar field @xmath131 . @xcite found a simple solution , representing @xmath42 black holes momentarily at rest , which gives @xmath132 where @xmath133 is the mass associated with the @xmath134th black hole and @xmath135 represents its coordinate center . each @xmath135 corresponds to the location of a puncture , as described in the last section . this is a valuable solution , but not generally useful , because these black holes lack momentum and spin . an explicit solution by bowen for the momentum constraint was the last crucial step in defining a procedure for calculating initial data for multiple black holes with specified linear and angular momenta @xcite . the bowen - york prescription employed a two - sheeted topology found by misner for the black - hole interiors @xcite . later brandt and brgmann generalized the procedure for the brill - lindquist topology @xcite . the brandt - brgmann puncture data for arbitrary momenta and spin is now widely used because of its ease of implementation . more recently york developed another modeling ansatz , known as the conformal - thin - sandwich ( cts ) approach @xcite , which has certain additional advantages . for example , it turns out the conventional brandt - brgmann puncture can not yield a spin parameter @xmath136 greater than @xmath137 @xcite , while cts data can go higher @xcite . as discussed above , the spatial metric is chosen to be conformally flat , then instead of providing an ansatz for the extrinsic curvature @xmath53 , the initial time - derivative of the conformal metric is specified ( generally to vanish ) . in addition , a condition can be imposed on the slicing ( see below ) . the result is a coupled system of elliptic equations which is solved to enforce the constraints and optionally the choice of slicing condition . boundary conditions may additionally be supplied to enforce rotational motion . either approach provides an ansatz for constructing three - dimensional binary black - hole initial _ field _ data for a specified choice of particle - like parameters , masses , positions , momenta , and spins . inevitably there are differences from the nearly quiescent evolving systems that we seek to represent . generally there is some level of spurious radiation generated from a period of initial transient dynamics through which the system relaxes to become quiescent on sub - orbital time scales . in particular , extraneous radiation content is an unavoidable consequence of conformal flatness , which post - newtonian analysis has shown must deviate from the physically relevant inspiralling binary solution @xcite . this spurious radiation is seen in plots of the gravitational waveforms produced by mergers ; see fig .  [ fig : universal_wf ] and other waveform plots in sec . [ sec : gws : universal ] . often the simulations will also undergo a period of initial gauge - evolution which , though physically inconsequential , may affect the quality of the simulation numerically . for puncture initial data @xcite analyzed some of these gauge dynamics , and developed a promising approach (  trumpet " data ) to mitigate it ( see sec . [ sec : movepunc ] for related gauge issues ) . however , for simulations lasting for several orbits these modest transient effects are generally negligible . most black - hole binary simulation studies are designed to represent the astrophysical population of systems , which have circularized before the gravitational radiation becomes observationally significant . these simulations begin with circularly inspiralling initial data configurations . even before reliable numerical relativity simulations were possible , considerable attention had been given to prescribing initial data for these near - circular configurations ( see @xcite for a review ) . within the cts approach , it is particularly natural to impose initially circular motion upon the system , resulting in an initial data prescription for quasicircular orbits @xcite . for either the cts or brandt - brgmann data , quasicircular parameters may be chosen by constraints on either an effective gravitational potential or total energy of the system @xcite . for more realistic inspiralling trajectories , the pn approximation may be used @xcite , or , for higher accuracy , an iterative procedure involving short numerical evolutions @xcite . the einstein evolution equations , which determine the time - development of the initial data , form a set of at least six coupled , nonlinear propagation equations . the exact formulation of these equations depends on the choice of evolved variables as well as how constraint identities are incorporated . although many formulations are possible , not all are equivalent from a numerical point of view . the choice of formulation can affect how numerical errors grow in time , and whether these errors can converge to zero as the resolution of the computational grid is increased . an ideal formulation is well - posed , meaning that : ( i ) for a given set of initial data , a unique solution exists , and ( ii ) solutions depend continuously on perturbations of the initial data . to determine whether a given formulation has the desirable property of well - posedness , we first consider it as a matrix equation : @xmath138 where @xmath139 is a vector , the components of which represent the evolved variables such as components of the metric ; the linear operator @xmath140 is the principal part , containing the highest - order spatial derivatives ; and @xmath141 contains additional terms which may be nonlinear . certain conditions on the eigenvalues and eigenvectors of the operator @xmath140 in eq . are used to characterize the equations , typically as either `` strongly hyperbolic '' or `` weakly hyperbolic '' @xcite . if the eigenvalues are real and the eigenvectors form a complete set ( so any solution can be written as a linear combination of them ) , the equations are strongly hyperbolic . then , assuming adequate boundary conditions , the equations are well - posed and all solutions are bounded by a function that grows ( at most ) exponentially at a rate independent of the initial data . strong hyperbolicity is a key ingredient for successful black - hole merger simulations .    on the other hand , if the eigenvalues are real but the eigenvectors are not complete , the equations are weakly hyperbolic . in this case the equations are not well - posed , and they permit solutions which grow at rates that depend on the initial data . weak hyperbolicity implies that small numerical errors in the initial data may grow at a rate which depends on the resolution . it then becomes difficult to show that the numerical solution converges to a well - defined continuum solution , and at high resolutions the simulation may become unstable @xcite .      the quest for workable formulations of einstein s equations has proceeded along two parallel lines of development . one originated with consideration of  harmonic coordinates " , so called because the coordinates satisfy the wave equation @xmath142 , where @xmath143 is the curved - space wave operator . in these coordinates , einstein s equations can be written such that the principal part resembles a wave equation in terms of the metric as in eq .  . in this form , einstein s equations are manifestly hyperbolic @xcite . however , the harmonic coordinate condition is too restrictive for numerical purposes , so generalized harmonic coordinates were eventually developed by introducing a source term into the coordinate condition , i.e. @xmath144 @xcite , a suitable choice for which preserves strong hyperbolicity . the subsequent introduction of _ constraint - damping _ terms , which tend to drive the constraints towards zero , further ensured stability @xcite . this formulation is manifestly second - order in both time and space , and has been implemented numerically as such @xcite , but for more efficient numerical integration a first - order - in - time formulation was also developed @xcite , and is currently being used by some groups .      a second line of development originated with the invention of the adm formulation of einstein s equations @xcite . this formulation was refined by @xcite , who suggested evolving the specific quantity of the extrinsic curvature ( eq .  ) . the associated evolution equations specify the first - order time - derivatives for the three - metric @xmath52 and the extrinsic curvature @xmath53 , as well as for the gauge variables @xmath77 and @xmath78 . however , the adm formulation is weakly hyperbolic ( see , e.g. , chapter 5 of @xcite ) , and attempts at stable numerical evolutions with it were not successful . there are essentially two ways to modify the adm equations that may affect their hyperbolicity ( while keeping them first - order in time and maintaining their 3 + 1 character ) . one way is to add the constraints , which may change the principal part of the equations without changing the associated physics . the other is to define new independent variables . simple addition of constraints failed to give a strongly hyperbolic formulation however . eventually it was found that strongly hyperbolic versions could be constructed from the adm equations by the promotion of certain derived quantities to the status of independently evolved variables @xcite . although strong hyperbolicity is important for a stable formulation , it is not the only property critical for accurate evolution . of the various strongly hyperbolic 3 + 1 formulations , one eventually emerged as more successful than the rest . its development began with the observation that the numerical accuracy of @xmath145 , the determinant of the three - metric , and @xmath146 , the trace of the extrinsic curvature , could best be preserved by evolving these two quantities independently @xcite . subsequently , to remove the redundancy in evolving the full three - metric and extrinsic curvature , the conformal three - metric @xmath147 and conformal traceless extrinsic curvature @xmath148 were substituted @xcite . to eliminate certain terms with second derivatives ( which contribute to the principal part of this system ) , a new independently evolved variable @xmath149 was also introduced @xcite . further improvements were made shortly afterwards , such as the addition of constraints to eliminate more second - derivative terms @xcite . the resulting formulation , now commonly known as bssn , proved to be robustly stable and accurate for numerical purposes . note that the above choices which comprise its form were empirically motivated , by identifying and eliminating terms which tended to compromise numerical accuracy . some analytic justification for its success was later found by @xcite , who showed that the bssn formulation avoids exponentially growing modes and most zero - speed modes , which accumulate numerical error . finally @xcite and @xcite proved that bssn is also strongly hyperbolic , given an appropriate choice of gauge .    a few noteworthy refinements of the bssn formulation followed . @xcite and @xcite suggested adding specific terms proportional to the  gamma constraint " ( @xmath150 minus its definition ) and the hamiltonian constraint into the evolution equations . these terms have the effect of damping the constraints and thereby improving stability . later , for evolution of punctures specifically , @xcite evolved the conformal factor in the form @xmath151 , which vanishes at the puncture ( an improvement in accuracy over the previous standard of @xmath152 , which is singular at the puncture ) . this change of variables , however , introduces the potentially singular factor @xmath153 into the evolution equations , necessitating an arbitrary restriction on the minimum value of @xmath154 . an alternative variable @xmath155 was subsequently suggested because resulting appearances of @xmath156 mostly cancel with factors of @xmath157 @xcite .      the choice of gauge or coordinate conditions , like the choice of formulation , has important consequences on the numerics , especially the stability of the simulation . important considerations include how to deal with the extreme conditions of black holes such as the physical singularities , the possible coordinate singularities , the strong - field gradients , and the dynamical , surrounding spacetime . the coordinates must accommodate these features in a way that is numerically tractable . an early consideration was how to avoid the physical singularity of a black hole through an appropriate slicing condition . one such choice , _ maximal slicing _ , required solution of a numerically expensive elliptic equation @xcite . later , an interest in constructing a hyperbolic evolution system led to a generalization of harmonic slicing known as the bona - mass family @xcite . one particular member of this family , called 1+log slicing , was found to also be singularity - avoiding , like maximal slicing , but at less numerical cost @xcite . meanwhile , it was recognized that a shift vector was required to counter the large field gradients , or `` slice - stretching '' , incurred in the presence of a black hole @xcite . a hyperbolic shift condition called the `` gamma driver '' was introduced that fulfilled this requirement @xcite : @xmath158 where @xmath159 depends on a conformal three - metric @xmath160 of the evolving spatial slice , @xmath161 is an auxiliary variable , @xmath78 is the shift , and @xmath162 is some scalar field . @xmath163 is a damping parameter that fine - tunes the growth of the shift , which affects the coordinate size of the black - hole horizons , which in turn has bearing on the required numerical resolution @xcite . this shift condition also has the desirable tendency to drive the coordinates to quiescence in synchrony with the physics , for example after a binary merges into a stationary black hole . clearly if @xmath150 vanishes , @xmath161 is damped to zero and @xmath78 approaches a ( typically small ) constant . initially these conditions were intended for use with punctures that remained at fixed coordinate positions ( i.e. comoving coordinates ) , by choosing @xmath162 to vanish at the punctures . however , this led to large gradients between the merging black holes since there the metric had to vanish as the physical distance contracted . pathologies also arose from the twisting of the coordinates as the black holes orbited each other . these issues , which resulted in large numerical errors and instabilities , were eventually resolved with the breakthrough discovery that slight modifications of the 1+log and gamma - driver conditions allowed arbitrary motion of the punctures with robust stability @xcite . in particular , in the shift condition , a critical alteration was  unpinning " the puncture by no longer requiring the factor @xmath162 to vanish at the puncture but rather to remain constant . these modifications were subsequently refined to eliminate zero - speed modes , thus preventing the possibility of error build - up @xcite . this was accomplished in part by the addition of advection terms : @xmath164 it is noteworthy that these particular gauge choices , coupled with the bssn equations , also result in strong hyperbolicity @xcite . a further refinement resulted from the observation that , assuming that @xmath165 initially , eqs .   and can be integrated to yield @xmath166 ; this substitution allows the removal of the equation for @xmath161 , leaving only a single shift equation @xcite : @xmath167 use of this or similar gauge conditions became known as the  moving puncture " method , and proved to be very successful as it became increasingly widespread among the numerical - relativity community . analysis has verified that moving punctures are valid black - hole solutions , although the initial character of the punctures changes significantly during the course of evolution . the second worldsheet shown in fig . [ fig : wormhole_flatspace ] , for example , invariably becomes disconnected in numerical simulations , at a point within the horizon near the throat , due to the action of the shift vector effectively shifting computational grid points onto the first worldsheet . meanwhile the spatial coordinates evolve such that the @xmath168 singularities in the conformal factor @xmath131 become @xmath169 singularities . for a single nonspinning black hole , the numerical result rapidly asymptotes to an exact form of the schwarzschild solution called a  trumpet " , recently investigated by @xcite and @xcite , which turns out to be a type of solution first considered by @xcite . development of generalized harmonic coordinates initially proceeded independently of the above 3 + 1-formulated conditions . as mentioned , in harmonic coordinates the dalembertian of each coordinate vanishes . in generalized harmonic coordinates , the wave equation for each coordinate is allowed a source term , i.e. @xmath170 these  gauge driving " source terms @xmath171 can be either algebraically specified or evolved such that hyperbolicity is preserved @xcite . the first successful numerical orbit of black holes involved a source term for the time coordinate that effectively kept the lapse close to its minkowski value of unity , while the spatial coordinates remained harmonic @xcite . this was accomplished by evolving the source term itself , according to @xmath172\alpha^{-1}\ ] ] where @xmath173 and @xmath174 are constants . more recently , to dampen extraneous gauge dynamics during the inspiral and merger of generic binaries , @xcite found the following gauge driver to be successful : @xmath175 ^ 3 \\ h_i & = & -\mu_0\left[\log\left(\frac{\sqrt{g}}{\alpha}\right)\right]^2\frac{\beta_i}{\alpha}\end{aligned}\ ] ] where @xmath176 is a specified function of time that starts at zero and eventually increases monotonically to unity . additional noteworthy advances in coordinate conditions were developed for the generalized harmonic formulation to facilitate spectral methods ( see sec . [ sec : comp : computation ] below ) but are in principle generalizable to other frameworks . one of these is the use of multiple coordinate patches , where the coordinates are chosen according to the local physical geometry for optimal accuracy . in the generalized harmonic formulation , characteristic speeds are readily available for use in constructing physical boundary conditions between patches @xcite . the other advance is that of  dual coordinates " . because moving punctures involve irregular fields , they are not easily made compatible with spectral methods , which are sensitive to irregularities . an alternative , that of a moving , excised region , can leave a large trail of interpolation error in its wake . the remaining option is to use comoving coordinates @xcite , but these can result in instabilities due to the steep field gradients and other pathologies mentioned in sec . [ sec : comp : gauge ] . dual coordinates were invented to exploit the advantages of both comoving and non - comoving coordinates while avoiding their disadvantages . this is accomplished by computing all tensor field components in the non - comoving coordinate system ( thus avoiding steep gradients ) , yet evolving them as functions of the comoving coordinates ( thus allowing stationary excision boundaries ) @xcite . the initial data , formulation , and gauge are all , in principle , analytic choices applicable to an infinite , continuous manifold . one must finally make choices pertaining explicitly to the finite , discrete mechanics of the computations that will numerically approximate the above analytic specifications . an immediate consideration is the fact that the simulated domain must have finite extent . various conformal compactification schemes which map spatial or null infinity to a finite boundary have been tested @xcite , but currently it is more common to impose artificial boundaries at finite spatial coordinates and apply some form of either radiative boundary condition @xcite or constraint - preserving boundary condition @xcite .    to mitigate the effect of inward - propagating errors from the outer boundaries , some form of mesh refinement is typically employed to push the outer boundaries as far away as possible from regions of interest ( e.g. wave sources and extraction regions ) . with a limited number of grid points available , their density is judiciously chosen to be highest near the strong field gradients of the black holes , moderate throughout regions where wave propagation is studied , and coarser beyond that . if the simulated black holes are very dynamic then some algorithm for automatically adapting the mesh refinement is necessary . meanwhile the interfaces between refinement regions require interpolation .    on this grid , spatial derivatives are computed in one of two ways . they may be computed with finite differencing stencils across uniformly spaced grid points , derived from taylor expansions , currently up to eighth - order accurate @xcite . alternatively , the spectral approach may be used , in which coefficients of an expansion in basis functions are computed to some order , on a number of collocation points comparable to the number of basis modes , from which the derivatives are then obtained analytically @xcite . in the former case , dissipative terms are often added to the evolution equations to reduce noise @xcite . in the latter case , spectral methods often include a smoothing step for a similar end . last , time must be advanced in discrete steps . various explicit time - integration algorithms have been tested , for example the iterated crank - nicholson scheme , but the fourth - order runge - kutta algorithm has become most widely used , due to its superior accuracy and efficiency . in some codes the timestep size is made to vary with spatial resolution , for even greater efficiency , albeit at the expense of the complication of time - interpolation . implicit timestepping schemes are also being investigated as a means to greatly increase the timestep size @xcite . one of the most important end - results of a simulation of merging black holes is a computation of the emitted gravitational radiation . for this purpose , it is useful to calculate the  weyl tensor " , @xmath177 . the weyl tensor , like the einstein tensor , is constructed from derivatives of the metric . unlike the einstein tensor , the weyl tensor has degrees of freedom that do not necessarily depend on a massive source ; it can be nonzero while the einstein tensor vanishes . certain components of the weyl tensor form a complex quantity called @xmath178 , one of five  weyl scalars " used to classify spacetimes @xcite . in a special ,  transverse - traceless " , spherical coordinate system , @xmath178 can be expressed as follows : @xmath179 where the subscripts @xmath180 denote orthonormal tetrad components . in a spacetime with gravitational radiation this quantity typically falls off as @xmath181 @xcite , and can be associated with outgoing radiation at spatial infinity ( @xmath182 ) in asymptotically flat spacetimes @xcite . in terms of the strain introduced in the last section , this can be written @xmath183 , \ ] ] where , in terms of a metric perturbation @xmath184 , @xmath185 what makes @xmath178 a particularly useful measure of the radiation is that to first order in @xmath186 it is coordinate - invariant . although the linearized radiation interpretation for @xmath178 is only strictly valid in the limit as @xmath182 , in numerical simulations it is often extracted on a sphere of large but finite radius . this will approximate the expected radiative behavior if @xmath187 ; if spheres of multiple radii are used then it is also possible to extrapolate the results to infinity . typically @xmath178 is computed on the computational grid points from the metric variables and then interpolated onto the spherical extraction surface . recently , cauchy - characteristic extraction methods have been applied to allow direct evaluation of radiation at future null infinity @xcite . the product of @xmath178 with a spherical harmonic function is then integrated over the sphere , as it is useful to extract specific modes of the field . the structure of the gravitational field is such that it is convenient to expand in spin - weight -2 spherical harmonics @xmath188 @xcite ( just as , analogously , it proves convenient to use the spin - weighted -1 harmonics in expanding the electromagnetic field ) . thus : @xmath189 unlike electrodynamics , there is no dipole moment of the gravitational field in conventional general relativity , so the dominant contribution is the @xmath190 quadrupole moment . other information that is important in analyzing spacetimes relates to the properties of the black holes themselves . this is commonly taken from apparent horizons , surfaces from which all light rays must go inward . it is useful to compute this during a simulation , for example by finding surfaces on which the expansion of light rays are minimized @xcite . when a black hole is sufficiently isolated , it is possible to define a mass and spin magnitude from data evaluated on the horizon @xcite , as well as linear momentum @xcite . these quantities can be used to characterize the physical properties of the black holes , and to label the spacetimes when comparing with pn theory . the final merger of a black - hole binary takes place in the dynamical , strong - field regime of general relativity . for many years , numerical relativists wondered what they would find when they probed this arena using black - hole merger simulations . what would the merger waveforms look like ? how strongly might artifacts in the initial data affect the mergers ? how might the effects of unequal masses and spins change the picture obtained by studying the simplest case of nonspinning , equal - mass mergers ? would they uncover any unexpected phenomena ? recent breakthroughs in black - hole merger simulations are providing important tools for addressing such questions .    in this section we explore the dynamics of black - hole mergers and the resulting gravitational waveforms . we start with a look at the lazarus approach , which provides hybrid waveforms by combining analytic methods with brief numerical simulations . we then begin a more comprehensive discussion of black - hole merger physics , starting with mergers of nonspinning , equal - mass black holes . next we consider mergers of unequal - mass holes , followed by mergers with spin . throughout this discussion we aim to provide a historical context while highlighting the key physical results .      in the late 1990s , numerical relativity had advanced sufficiently to allow brief simulations of two black holes in three dimensions . by `` brief '' we mean here durations of @xmath191 , which is a small fraction of the estimated binary orbital period of @xmath192 near the isco . while most of the community focused on technical developments aimed at extending the duration of these simulations , a small band of collaborators took a different approach @xcite . their novel idea was to use full numerical - relativity simulations to calculate the strong - field approach to merger during @xmath191 of simulation time , and then to calculate the remaining evolution using perturbative techniques @xcite . they began with traditional puncture black holes ( which remain fixed in the grid ) on quasicircular orbits near the isco , and evolved them towards merger . then , just before the simulation became unstable , they stopped the calculation and extracted the physical data for the merging black holes and the emerging gravitational waves . this data was then interpreted as initial data for a highly distorted single black hole , and evolved using techniques of black - hole perturbation theory . they called this method of reviving a ( nearly ) dead calculation the _ lazarus _ approach . the hybrid waveforms produced by the lazarus collaboration gave the first indications of what might be expected from the final merger of black - hole binaries @xcite . they ran a suite of simulations , varying the initial black - hole separations and the time at which they stopped the calculation and made the transition from the fully numerical evolution to the perturbed black hole evolution . figure  [ fig : laz_wf ] shows the dominant quadrupole @xmath193 component of re(@xmath178 ) , which corresponds to the @xmath27 polarization , for the case of equal - mass , nonspinning black holes . note the simple shape of the waveform smoothly tying together what might be the end of an ( inspiral ) chirp with a damped ( ringdown ) sinusoid .    by tuning the mass and spin of the background hole for the perturbative evolution , @xcite also determined that the merger remnant was a kerr ( spinning ) hole of mass @xmath194 , and spin @xmath195 . the lazarus method was also applied to mergers of black holes with spins either aligned or anti - aligned with the orbital angular momentum @xcite ; in all cases , a similar waveform shape was seen . how generic was this simple shape and , in particular , would it also arise in situations where the black holes complete one or more orbits before merging ? part of @xmath196 , corresponding to the @xmath27 polarization ( measured on the system s axis ) , is shown for ten simulations : five different initial black - hole separations ( designated qc-0 , etc . ) , each with two transition times to perturbative evolution . ] early in 2005 , frans pretorius stunned the community by achieving the first robust simulation of two equal - mass black holes through a single plunge orbit , merger , and ringdown @xcite . the resulting gravitational waveforms are shown in fig . [ fig : frans_wf ] , where @xmath25re(@xmath178 ) is plotted versus @xmath197 . here @xmath25 is the coordinate distance from the center of the grid to the sphere on which the waveform is extracted ( see sec .  [ sec : extract ] ) , and @xmath198 is the mass of a single black hole . ( note that the time axes for all other waveform plots in this paper are scaled by either the total system mass @xmath9 or the mass of the final , remnant black hole @xmath199 ) . , due to lower resolution in the outer regions ]    here curves are plotted for waveforms extracted on four spheres with successively larger radii . the curves have been shifted in time so that the waveforms overlap . note that the overall waveform shape is simple . the merger yields a single black hole that is spinning , with spin parameter @xmath200 , where @xmath199 is the mass of the final black hole that forms .    to carry out these mergers , pretorius @xcite used techniques that are very different from the more traditional approach based on bssn and punctures , used by nearly all other numerical relativists . as discussed in sec . [ sec : harm ] , he used a generalized harmonic formulation in which the einstein equations are written with second - order time and space derivatives . the spatial domain was compactified , so that the outer boundary of each slice was mapped to spatial infinity . moreover , he excised the black holes ( each formed from the collapse of scalar field `` blobs '' ) and evolved their motion across the grid using adaptive mesh refinement ( see secs . [ sec : bhongrid ] and  [ sec : comp : computation ] ) . pretorius had clearly developed a robust method for evolving black - hole binary mergers . in the wake of his remarkable success , many numerical relativists began studying his methods . however , another surprise was soon to emerge .    in the autumn of 2005 two research groups , one at utb and the other at goddard , independently developed a powerful technique for evolving mergers of puncture black holes within the traditional bssn approach @xcite . as discussed in sec . [ sec : movepunc ] , the standard puncture method requires the black holes to remain fixed on the numerical grid . in this new `` moving puncture '' approach , simple but novel gauge conditions allow the punctures to move across the grid , producing accurate and stable merger evolutions @xcite .    although the utb and goddard codes were both based on the bssn approach , they were independently written and featured somewhat different implementations of moving punctures . in addition , the goddard code used second - order finite differencing ( standard in the community at the time ) . they employed a box - in - box fixed mesh refinement to produce the high resolution in the region around the black holes needed to compute the dynamics accurately , while maintaining adequate resolution in more distant regions and a large enough computational domain to allow accurate extraction of the gravitational waves . the utb code used fourth - order finite differencing , and the `` fish - eye '' coordinate transformation to produce higher resolution around the black - hole orbits and lower resolution in the wave extraction regions . both groups successfully evolved an equal - mass , nonspinning binary through the final plunge orbit , merger , and ringdown , and extracted the gravitational waveforms . in these simulations , the black holes completed @xmath201 orbit before merger , defined as the time at which a single apparent horizon was formed . figure  [ fig : utb_tracks ] shows the tracks traced by the black - hole punctures , with the apparent horizons superimposed , for the utb run @xcite . these simulations produced plunge waveforms with a simple shape similar to that found by pretorius ; compare fig . [ fig : frans_wf ] from pretorius s simulations with fig . [ fig : gsfc_wf ] , which shows waveforms from the goddard simulation plus a corresponding lazarus waveform .     and @xmath202 . the solid circles denote the centroids of the apparent horizons every @xmath203 during the evolution . the first common horizon , marking the time of merger , forms at @xmath202 , just before the punctures complete @xmath204 orbit . ]    , @xmath205 mode of @xmath206 was extracted from the numerical simulation on spheres of radii @xmath207 , and @xmath208 for the medium- and high - resolution runs . the waveforms extracted at different radii have been rescaled by @xmath25 , and shifted in time to account for the wave propagation time between the two extraction spheres . curves are shown for medium ( @xmath209 in the finest grid ) and high ( @xmath210 ) resolutions . a waveform from a lazarus calculation starting with the same initial data is shown for comparison @xcite . the lazarus waveform shown in this figure is scaled differently from those in fig . [ fig : laz_wf ] . ]      in most cases , astrophysical black - hole binaries spiral together over many orbits , radiating away any orbital eccentricity in the form of gravitational waves . by the time the black holes reach the final inspiral , their orbits will be quasicircular . all mergers of such equal - mass nonspinning binaries should thus produce the same waveform signature , subject only to rescaling with the total system mass @xmath9 .    over the years , concerns had been raised within the relativity community about the effects of deviations from astrophysical initial conditions on the waveforms . for example , spurious eccentricity in the orbits could arise from starting conditions that did not approximate quasicircular inspiral accurately . moreover , the commonly - used conformal flatness prescription ( see sec .  [ sec : comp : initialdata ] ) for the initial data is different from the conditions experienced by an astrophysical binary ; these differences would result in spurious gravitational radiation being present in the initial conditions for the binary simulations . how would factors such as these influence the merger waveforms ?    having developed robust techniques for evolving black - hole mergers , numerical relativists eagerly pursued longer simulations with the black holes starting at wider separations , completing more orbits , and producing longer wavetrains . the utb group ran a simulation in which the black holes completed nearly 1.5 orbits before the formation of a single apparent horizon @xcite . they observed that the resulting waveform was quite similar to those produced by their earlier simulation that started near the plunge @xcite , with a brief oscillatory signal at the beginning . @xcite anticipated `` that the plunge waveform , when starting from quasicircular orbits , has a generic shape that is essentially independent of the initial separation of the binary . '' the goddard group , simultaneously pursuing this same goal , produced the first `` universal '' waveform for the merger of equal - mass , nonspinning black holes @xcite . using approximately quasicircular initial conditions , they ran a series of four simulations , starting the black holes at increasingly wider separations . in all , the binaries completed @xmath211 , 2.5 , 3.6 , and 4.2 orbits before forming a common apparent horizon ; the systems emitted just under 4% of their energy as gravitational radiation , and the final black holes were spinning with @xmath212 . to compare the results of these models , they chose the fiducial time @xmath102 to be the moment of peak amplitude in the gravitational radiation ; this typically occurs within a few @xmath9 of merger . it is worth noting how close these results are to those obtained from the lazarus approach to the same system ( see sec . [ sec : gws : laz ] ) . orbit before merger ( denoted by * ) . ]    the binary orbital dynamics of these runs is clearly seen in fig . [ fig : tracks ] , which shows the tracks of the black - hole punctures . here only one puncture track is shown for each binary , with the trajectories oriented to superpose at the fiducial @xmath102 . the tracks show the effects of eccentricity in the early stages of each run ; plots of the separation @xmath25 versus time show that the initial eccentricity decreases as the initial separation increases . as the holes spiral together into the strong - field regime , the eccentricity diminishes . the puncture tracks lock on to a single universal trajectory , independent of the initial conditions , through the final orbit , plunge , and merger . figure  [ fig : universal_wf ] shows the corresponding universal gravitational waveform for equal - mass , nonspinning black holes . specifically , it shows the real part of the @xmath190 , @xmath205 mode of @xmath206 versus time ; for this case , the quadrupole @xmath213 modes strongly dominate all other modes . here the signals produced by each run have been shifted in time so that @xmath102 marks the moment of peak radiation amplitude . starting from @xmath214 , the waveforms for the final burst of radiation show nearly perfect agreement , with differences at the level of @xmath215 . for the preceding few orbits , the waveforms agree in amplitude and phase to @xmath216 , except for the brief initial bursts of spurious radiation ( see sec .  [ sec : comp : initialdata ] ) . overall , the merger stage lasts @xmath57 , converting @xmath217 of the initial total mass @xmath9 into gravitational - wave energy . the gravitational wave released during this burst has a luminosity @xmath218 , which is greater than the total luminosity of all the stars in the visible universe . for stellar black - hole binary mergers , this luminosity will last for a few milliseconds , while for mbh binaries it will last for several minutes . , @xmath205 mode ) produced by the four simulations whose puncture tracks are shown in fig . [ fig : tracks ] @xcite . the signals have been shifted in time so that the peak radiation amplitude occurs at @xmath219 . the agreement among the waveforms is excellent , with differences @xmath215 for the final burst of radiation starting at @xmath220 . ] note that this universal waveform has a simple shape . the signal starts with a short inspiral chirp , increases smoothly in amplitude during the plunge and merger , and then transitions to the damped ringdown sinusoid . overall , the frequency increases monotonically , reaching a maximum value that stays constant during the ringdown @xcite . since black - hole - merger waveforms will be used to help find signals in data from gravitational - wave detectors , it was important to verify that various groups consistently achieve the same results . figure  [ fig : compare_wf ] provides a comparison of waveforms computed by pretorius and the groups at utb and goddard for equal - mass , nonspinning black - hole mergers @xcite ; here the real part of @xmath206 s dominant @xmath221 , @xmath205 mode is shown , with all three waveforms aligned in time so that the moment of peak gravitational radiation amplitude occurs at @xmath222 . all three simulations evolved through the last @xmath223 orbits before merger ; the goddard run is r1 from @xcite and shown in figs . [ fig : tracks ] and  [ fig : universal_wf ] . note that all evolutions produce the same overall waveform shape . the goddard and utb groups evolved black holes with zero spin starting on similar quasicircular initial orbits , and produced nearly identical waveforms . pretorius used co - rotating initial conditions which impart a small spin @xmath224 to each black hole . this difference is consistent with the slightly higher frequencies seen in pretorius results during the ringdown , @xmath225 .    . the differences between pretorius waveform and the others during the ringdown , @xmath226 , are consistent with his simulation producing a slightly faster rotating final black hole due to these initial spins . reproduced by permission of the institute of physics . ] how robust is this universal waveform to larger amounts of orbital eccentricity and gravitational radiation in the initial data ? @xcite studied a series of equal - mass , nonspinning binaries starting from roughly the same initial orbital period and having varying amounts of initial orbital eccentricity @xmath227 . for @xmath228 , the initial eccentricity is radiated away , and the binaries circularize and begin a universal plunge at @xmath229 before the time of peak radiation amplitude , producing a final black hole with @xmath230 . for @xmath231 , the black holes do not complete any orbits , but rather plunge together and merge ; the final black hole spin reaches a maximum value @xmath232 around @xmath233 . figure  [ fig : psu - ecc - tracks ] shows the puncture tracks , and fig .  [ fig : psu - ecc - wfs ] shows the real part of the strain s @xmath221 , @xmath205 mode , for two nonspinning , equal - mass black - hole mergers with nonzero eccentricity from this study ; in both cases , @xcite show quantitatively that the waveform for the final burst of radiation is the same as the universal waveform .    in a complementary study , @xcite found that the final spin of the remnant black hole is essentially insensitive to the eccentricity for binaries that do not plunge immediately . @xcite examined the effects of spurious gravitational waves by superposing a tunable packet of gravitational radiation on an equal - mass , nonspinning binary . they found the binary evolution and the spin of the remnant black hole to be robust to modest amounts of added gravitational radiation . ( left ) and @xmath234 ( right ) from @xcite . in these cases , the initial eccentricity is radiated away and the binaries circularize before merger . ] strain mode ( labeled @xmath32 here ) for black - hole mergers shown in fig . [ fig : psu - ecc - tracks ] from @xcite . these runs have initial orbital eccentricity @xmath235 ( top ) and @xmath234 ( bottom ) and enter a universal plunge by @xmath229 before the gravitational radiation reaches its peak value . ]      the black - hole simulations discussed so far cover the last @xmath236 four orbits before merger . these achievements marked the triumph of solving a long - standing problem of fundamental importance to general relativity : the two - body problem for the final merger of equal - mass schwarzschild black holes driven by gravitational radiation reaction . these simulations provided the first look at the dynamics and waveforms of black holes evolving and merging in the nonlinear , strong - field regime . they also pointed the way towards applications in detecting the final gravitational - wave burst from black - hole mergers @xcite . however , for more advanced applications to gravitational - wave data analysis and comparisons with pn analytical treatments , longer waveforms that start many cycles before merger are essential . the goddard group @xcite produced the first such long waveforms for equal - mass , nonspinning black holes starting @xmath58 orbits or @xmath237 gravitational - wave cycles before merger . they used improved initial conditions with low eccentricity @xmath238 , and focused on improving accuracy while the black holes traversed the relatively long inspiral . they investigated the observability of black - hole mergers with ground- and space - based gravitational - wave detectors @xcite ( see sec.[sec : gwda : direct ] ) and successfully applied their long waveforms to comparisons with pn results , focusing on the waveform phases @xcite . and applies when the black holes are far enough apart that the black hole speeds remain well below the speed of light . we refer to pn results by the order at which the series is truncated . for example , `` 2  pn '' means that terms of order @xmath239 are retained . see sec . [ sec : nrpn ] for a deeper discussion of the pn approximation in the context of numerical relativity . ]    shortly thereafter , the jena group simulated a binary inspiralling for nine orbits ( @xmath240 gravitational - wave cycles ) before merger @xcite . with this , they made the first quantitative comparisons with both the pn phase _ and _ amplitude , and quantified the level of error in the quadrupole approximation . they used higher - order finite differencing @xcite and initial binary parameters calculated using the pn approximation to reduce the initial eccentricity significantly @xcite , enabling a precise measurement of the waveform phase . the caltech - cornell group currently holds the record for the longest and most accurate black - hole binary evolution , starting 16 orbits and 32 gravitational - wave cycles before merger @xcite . using their spectral code ( see sec .  [ sec : comp : computation ] ) , they begin with a very small initial orbital eccentricity @xmath241 and evolve with very high accuracy through a relatively long inspiral , then merger and ringdown . the impressive trajectories of their black holes are shown in fig . [ fig : cc_traj ] and the accompanying gravitational waveforms in fig . [ fig : cc_wf ] .         showing the @xmath190 , @xmath205 component of re(@xmath206 ) the left panel shows a zoom of the inspiral waveform , and the right panel a zoom of the merger and ringdown . ] comparison of waveforms from different groups remains important in the push for higher accuracy and use in gravitational - wave data analysis . the samurai project @xcite sets the current state of the art for studying the consistency of black - hole binary waveforms . this effort focuses on comparing waveforms from equal - mass , nonspinning binaries , starting with at least six orbits ( or 12 gravitational - wave cycles ) before merger and continuing through the ringdown . they focus on comparing the amplitude @xmath242 and phase @xmath243 for the @xmath190 , @xmath205 mode of @xmath206 , defined as @xmath244 the gravitational - wave frequency of this mode is then @xmath245 . they compare the results from five independent numerical codes : the moving - puncture codes from the aei , goddard , jena , and penn state groups , and the caltech - cornell spectral code . figure  [ fig : samurai - phase ] compares the gravitational - wave amplitudes and fig . [ fig : samurai - amps ] the gravitational - wave phases as a function of frequency for the five waveforms . qualitatively , the results appear to be quite consistent . quantitatively , they concluded that these waveforms have sufficient accuracy to be used for detection with all current and planned ground - based detectors . as a function of the dimensionless frequency @xmath246 for the five codes from the samurai comparison @xcite . the scale along the top of the panel labels the frequency in khz scaled with respect to the total binary mass in solar units . ]    , @xmath205 component of the gravitational waves produced by the five codes in the samurai comparison @xcite are shown as a function of frequency . ] early in 2006 , numerical relativists took the next step in opening up the parameter space of binary black - hole mergers by simulating nonspinning binaries with unequal masses . this added a new parameter , the mass ratio @xmath247 , to the problem and brought an additional need for adaptive mesh refinement to achieve adequate resolution around the smaller black hole . as shown in fig . [ fig : unequal_tracks ] , the smaller black hole also moves faster , completing an orbit around the center of mass in the same time as the larger hole and thus requiring smaller timesteps for its evolution . these factors combine to make simulations of unequal - mass binaries more technically challenging ; currently , numerical relativists are able to simulate mass ratios up to @xmath60 @xcite .      decomposition into spin - weighted spherical harmonic modes ( see sec . [ sec : extract ] ) provides the basis for an in - depth study of black - hole mergers . for the equal - mass case @xmath248 , the @xmath221 , @xmath249 quadrupole mode is dominant and the odd-@xmath250 modes are suppressed by symmetry . as @xmath247 increases , the sub - dominant modes become more important , affecting both the evolution of the source and the emitted radiation . @xcite carried out a set of unequal - mass mergers with mass ratios in the range @xmath251 . they found that , to leading order , the total energy emitted during merger scales @xmath252 while the spin of the final black hole scales @xmath253 , where @xmath254 is the symmetric mass ratio . they also studied the multipolar structure of the gravitational waves . figure  [ fig : berti - amps ] shows several @xmath255 modes of the radiation produced in their simulations for the case @xmath256 . they demonstrated that the higher modes carry larger fractions of the total energy as @xmath247 increases ; in particular , the @xmath257 mode generally carries @xmath216 of the emitted energy for @xmath258 . the goddard group @xcite performed a complementary study of nonspinning unequal - mass mergers for mass ratios @xmath259 . they found that the overall simple waveform shape first discovered for equal - mass mergers extends to the cases with @xmath260 ; this is easily seen when the gravitational waves are scaled by @xmath254 . figure  [ fig : wfs_eta_scaled ] shows the strain @xmath261 , including the @xmath221 and @xmath262 modes , for an observer located at distance @xmath25 on the system s orbital axis . the waveforms are aligned so that the maximum amplitude of the @xmath190 , @xmath205 mode occurs at @xmath102 . the differences in phase during the final ringdown portion of the waveforms are consistent with the differing spins of the final black holes ; for example , with nonspinning black holes , a larger mass ratio @xmath247 results in a smaller final spin for the remnant black hole . when viewed off - axis , the waveforms show modest amplitude and phase modulations , while still preserving the overall simplicity in shape .     for different mass ratios , from @xcite . @xmath32 is calculated for different mass ratios using the @xmath221 and @xmath262 modes . the observer is located at distance @xmath25 along the axis of the system . ] @xcite also found that each of the individual spherical harmonic components is circularly polarized during the inspiral , merger , and ringdown , as seen by distant observers on the system s rotational axis . during the inspiral , the phase and frequency of the different @xmath255 components are the same for each mass ratio ; this is expected since the waveform phase is directly related to the orbital phase . more interestingly , for the @xmath263 modes , this relationship continues to hold during merger and ringdown .    based on these observations , they developed a simple conceptual picture in which each @xmath255 mode of the gravitational radiation is produced separately by the @xmath255 mode of some _ implicit rotating source_. the fixed relationship between the phase and frequency of the @xmath263 modes shows that these components of the implicit source rotate rigidly during the entire coalescence from inspiral through ringdown . in comparison , the @xmath264 source components peel away from the main trend of the @xmath263 parts during merger , indicating less rigid rotation . @xcite have carried out the highest mass ratio merger simulation to date , with @xmath60 . this binary radiates @xmath265 of its mass as gravitational radiation as it undergoes @xmath266 orbits before merging to form a black hole with @xmath267 ; these values fit the scaling relations found by @xcite and @xcite . figure  [ fig:10to1_wfs1 ] , from @xcite shows that the @xmath268 and @xmath269 modes of re(@xmath270 ) exhibit the simple waveform shape individually ; this trend continues through @xmath271 , the highest mode they studied . unequal - mass black - hole binaries can form as a result of galaxy mergers , dynamical processes in star clusters , and the final stages of binary stellar evolution . mergers of unequal - mass binaries bring a qualitatively new phenomenon of great importance to astrophysical scenarios of black - hole growth and retention : recoils or kicks .    as discussed in sec . [ sec : astro_grav : basic ] , the gravitational waves emitted by a merging black - hole binary carry away linear momentum . if the pattern of gravitational - wave emission is not symmetrical ( _ i.e. , _ if there is more radiation emitted in some direction than in others ) , then global conservation of momentum requires the center of mass , and thus the remnant black hole that forms , to recoil in the opposite direction . the situation for an unequal - mass , nonspinning binary is shown schematically in fig . [ fig : schematic_kick ] , from @xcite , who attributed this argument to alan wiseman ; we summarize their discussion here . the two holes are orbiting in a plane about their center of mass . since the smaller black hole moves faster , its wave emission undergoes more `` forward beaming '' than that of the larger hole . instantaneously this gives a net flux of momentum parallel to the smaller hole s velocity , and an opposing recoil or kick at the center of mass . the direction of this kick changes continually as the black holes orbit . if the orbit was circular , the center of mass would move in a circle and suffer no net recoil . however , since the black holes are spiralling together ( due to the energy and angular momentum carried away by the gravitational waves ) , the instantaneous kicks do not cancel exactly but rather accumulate . this causes the final merged black hole to have a non - zero net recoil in the orbital plane . is emitted parallel to the smaller hole s velocity . momentum conservation requires that the system center of mass recoil in the opposite direction , @xmath272 . reproduced with kind permission from springer science+business media : fig .  1 of @xcite . ]    this recoil has been studied by several authors , including @xcite . more recent analytic treatments using pn approximations were carried out by @xcite , while some comparable - mass estimates were also made with the lazarus approach @xcite . however , since the dominant part of the kick depends sensitively on the strong - field regime close to merger , where the orbits are less circular , an accurate calculation of the kick requires full numerical relativity simulations of the final stages of binary inspiral , merger and ringdown @xcite . the goddard group @xcite carried out the first accurate calculation of black - hole recoil for the merger of a nonspinning binary with @xmath273 . figure  [ fig : gsfc - kick.eps ] shows the kick velocity as a function of time for three simulations in which the black holes start out on orbits with increasingly larger initial separations @xmath274 . the kick velocity grows rapidly during merger , reaching a peak value @xmath191 after the time of peak radiation amplitude in @xmath178 , and then dropping to a lower , final value . it is clear from this figure that the black holes must start far enough apart in order to get a consistent and reliable value for the final recoil velocity . for this @xmath273 case , their results spanned the range @xmath275 , with a best estimate of @xmath276 .    following on from this , and other unequal - mass simulations by @xcite , @xcite performed the first systematic parameter study of nonspinning black - hole binary mergers with mass ratios in the range @xmath277 . they found that the maximum recoil velocity @xmath278 occurs for the mass ratio @xmath279 . for the case @xmath60 , @xcite found a recoil velocity @xmath280 . astrophysical black holes are generally expected to be spinning . including the spin vector of both black holes in the binary adds six more dimensions to the parameter space , giving seven in all when the mass ratio @xmath247 is included . exploration of this large parameter space began in 2006 , with the simplest cases of equal masses and spins , and is gradually growing to incorporate more generic binaries . the utb group carried out the first black - hole binary mergers with spin @xcite . they simulated three binaries , each having equal masses and starting on quasicircular orbits with initial orbital frequency @xmath281 , giving an initial orbital period @xmath282 . in the aligned case , both holes have spin parameters @xmath283 and vectors parallel to the orbital angular momentum @xmath284 ; in the anti - aligned case , the black holes have the same spin parameters but the vectors are anti - parallel to @xmath284 . a nonspinning equal - mass binary was run for comparison . the @xmath221 , @xmath205 components of the waveforms are plotted in fig .  [ fig : utb - spin ] , which shows all three cases starting out at the same time @xmath102 at the same orbital period and gravitational - wave frequency @xmath285 . notice that the aligned system takes longer to merge , completing more orbits and producing a longer wavetrain . this behavior is caused by a spin - orbit interaction that produces an effective repulsive force between the black holes . in the anti - aligned case , this interaction yields an effective attraction , causing the binary to merge more quickly and with fewer orbits . the nonspinning case is intermediate between the other two . all three cases produce similar waveforms with a clean , simple shape . mode is shown for three cases : aligned spins ( solid / red ) , anti - aligned spins ( dotted / black ) , and no spins ( dashed / blue ) . all three simulations start with the same initial orbital period at @xmath102 . figure kindly provided by m. campanelli . ] it is also interesting to compare the total energy radiated as gravitational waves @xmath286 and the spin parameter of the final black hole @xmath287 for these three cases . the aligned case yields @xmath288 and @xmath289 while the anti - aligned case gives @xmath290 and @xmath291 . the nonspinning case is again between these two cases , with @xmath292 and @xmath212 @xcite . the black - hole mergers discussed so far have resulted from simple binary dynamics with no precession and have produced simple waveforms with a similar shape , at least when examined mode - by - mode in a spin - weighted spherical harmonic decomposition . might more `` generic '' precessing binaries generate waveforms with more complex patterns ? this important question pertains not only to basic orbital dynamics in general relativity , but also to strategies for developing templates to search for signals in data from gravitational - wave detectors ( see sec . [ sec : gwda ] ) . the parameter space for such binaries is very large , and explorations of this space have only recently started . the utb group , now at the rochester institute of technology ( rit ) , pioneered the study of precessing black - hole binary mergers . they evolved a generic black - hole binary with unequal masses and unequal , non - aligned precessing spins that undergoes @xmath293 orbits before merger and produces a relatively long wavetrain @xcite . this system has a mass ratio @xmath294 and arbitrarily oriented spins with magnitudes @xmath295 and @xmath296 . the initial conditions were determined from point - mass evolutions using 2.5  pn and 3.5  pn parameters . figure  [ fig : rit-3d - tracks ] shows the difference in the black - hole trajectories @xmath297 . for a non - precessing binary , @xmath297 would have no component in the @xmath34 direction . in this generic case , the precession of the total system spin drives a precession of the orbital plane , producing evolution of the trajectory in the @xmath34 direction .     for the generic binary evolution from @xcite . the precession of the system spin induces precession of the orbital plane . ] the resulting waveforms demonstrate the effects of precession on the amplitudes of the sub - dominant modes . figure  [ fig : rit - precess ] shows the @xmath221,@xmath298 mode of the strain @xmath32 . the numerical - relativity evolution is shown by the solid curve ; note the amplitude modulations induced by the precession . the dotted curve shows a solution of the pn equations of motion at 3.5  pn order . , @xmath298 mode of the strain ( @xmath32 ) is shown for the numerical - relativity simulation ( solid ) and a 3.5  pn calculation ( dotted ) . the amplitude oscillations are induced by precession . ]      the merger of a binary consisting of nonspinning black holes produces a single spinning black hole with spin direction aligned with the orbital angular momentum @xmath284 of the binary . in this case , the spin of the final hole arises purely from @xmath284 . when the individual black holes have spins that are not aligned with @xmath284 , the spin of the final hole will , in general , not be aligned with the initial spin but rather will have a different direction . this phenomenon , in which the spin direction of the final black hole differs significantly from the spin directions of the individual holes prior to merger , is known as a `` spin flip . '' the simplest situation in which a spin flip can occur is a binary with both spins anti - aligned with @xmath284 . the case of generic binaries , with unequal masses and arbitrarily - oriented spins has a much larger parameter space . the rit group @xcite has examined spin flips for several such cases . figure  [ fig : rit-3d - flip - tracks ] shows the puncture trajectories for an equal - mass binary with equal spins having magnitude @xmath299 and initially pointing @xmath35 above the initial orbital plane . the black holes execute @xmath300 orbits before merger , during which time the spins , shown as arrows along each trajectory , precess by @xmath301 . the evolution of the spin direction for this case is shown in fig .  [ fig : rit - spin - flip ] . the final black hole has @xmath302 with a spin direction flipped by @xmath303 with respect to the component spins just before merger .    in a related phenomenon , the direction of the _ total _ angular momentum ( @xmath304 ) may change . this case was studied using general principles and extrapolations from point - particle results by @xcite , who predicted that binaries with @xmath305 will experience a flip in total angular momentum direction provided that the initial spins are equal and anti - aligned with @xmath284 , and that the individual black hole spins obey @xmath306 . in particular , they predicted that , for @xmath307 , a schwarzschild black hole should form if @xmath308 .    motivated by this work , @xcite simulated a series of binaries with @xmath309 and having equal anti - aligned spins in the range @xmath310 $ ] . they found that a schwarzschild hole is formed when the initial spin is @xmath311 .     above the initial orbital plane . the spin directions are shown as arrows along each trajectory every @xmath312 until merger . ]    , from @xcite . the spin direction of the component holes , shown by the red arrows plotted every @xmath312 until merger , varies continuously due to precession . the blue arrow shows the final black - hole spin , with direction flipped discontinuously from that of the component holes just prior to merger . ] the mergers of asymmetrical spinning binaries will produce recoiling remnant black holes . in the simplest case , the binary components have equal masses and spin magnitudes , and the asymmetry arises from the spin vectors pointing in opposite directions . in early 2007 , several new results appeared in quick succession , generating considerable excitement among numerical relativists as well as astrophysicists eager to apply these results . in this section we discuss the basic results and defer to sec . [ sec : astro ] their impact on astrophysics . @xcite simulated the merger of equal - mass black holes with equal spin magnitudes @xmath313 and having one spin vector aligned and the other anti - aligned with the orbital angular momentum . these mergers produced kicks directed purely in the orbital plane with magnitudes @xmath314 . shortly thereafter @xcite showed that , for the specific case @xmath315 , @xmath316 ; linear extrapolation of results to the maximally spinning case @xmath317 yields @xmath318 . @xcite carried out a systematic study of equal - mass mergers with spin magnitude @xmath319 and direction aligned with the orbital angular momentum . the other black hole has spin magnitude @xmath320 , and direction both aligned and anti - aligned . they found a maximum kick velocity @xmath321 .    these results all demonstrated that mergers of spinning black holes can produce significantly larger kick velocities than nonspinning mergers however , new results would soon show that kicks from the mergers of spinning black holes could get much larger indeed . the idea of `` superkicks '' was first discussed by the rit group @xcite . they observed that a pn treatment @xcite predicts that the recoil due to spin is maximized when @xmath322 , @xmath323 , and the spin directions are anti - aligned with each other and lying in the orbital plane . the jena group first simulated binaries in this configuration and demonstrated the phenomenon of superkicks . using @xmath324 , they found a resulting kick velocity @xmath325 @xcite . figure  [ fig : jena - kick ] shows the coordinate positions of the black - hole punctures from one of their simulations . notice that the trajectories move out of the initial orbital plane and that the final black hole is kicked in the @xmath326 direction . the rit group @xcite carried out similar simulations with @xmath327 . their results show kicks perpendicular to the orbital plane with magnitudes up to @xmath328 ; using the expression from @xcite they predict a maximum recoil @xmath329 for the case of maximally spinning black holes , @xmath317 . @xcite performed a multipolar analysis of recoil from black - hole mergers , for both unequal masses and non - zero , non - precessing spins . they found that multipole moments up to and including @xmath330 are sufficient to accurately reproduce the final recoil velocity ( within @xmath331 ) , and that only a few dominant modes contribute significantly to it ( within @xmath332 ) . @xcite studied the role of spin in producing superkicks . they showed that the recoil velocity is almost entirely due to the asymmetry between the @xmath333 and @xmath334 modes of @xmath178 . in addition , the major contribution to the recoil occurs within a period @xmath335 before and after merger , after the time at which a standard pn treatment of the evolution breaks down .    , initially oppositely directed in the orbital plane . during the evolution , the black holes move out of the original plane , and the final black hole recoils with velocity @xmath336 in the negative @xmath34-direction . ] einstein s theory of general relativity implies predictions for the dynamics and gravitational - wave generation for generic black - hole binaries at all stages of the merger process . as we have seen , numerical relativity now provides an excellent tool for concretely deriving these predictions beginning in the last orbits and continuing through the conclusion of the coalescence in the ringdown of the remnant black hole . there is no hard limit to how early numerical relativity can be applied in these simulations , but the considerable costs of these simulations tends to limit the duration to tens of orbital periods . though the pn expansion provides only an approximation to general relativity , it succeeds in making accurate and efficient predictions over very long time - scales , while the black holes remain well separated . late in the merger , however , as the black - hole velocities increase , these errors grow , and pn theory is no longer reliable    a full understanding of general relativity s predictions for the merger process must draw on both approaches . the maturation of numerical relativity has fueled efforts toward synthesizing a complete understanding of black - hole mergers built from the results and techniques of both these fields . even a coarse general review of the pn formalism and results is outside the scope of this presentation . decades of work in this area have already been covered @xcite . here we focus exclusively on key synergistic areas where fruitful research interactions drawing from both numerical relativity and pn results and formalisms are yielding a more complete general - relativistic understanding of black - hole binary systems . the broad range of research that interfaces these two approaches can be grouped into four categories , based on the different ways that pn theory impacts numerical - relativity research . pn theory : 1 ) provides independent results for cross - checking and comparing with nr ; 2 ) provides models for the initial values needed to begin astrophysically realistic nr simulations ; 3 ) provides insight for interpreting nr results ; and 4 ) may provide the basis for empirical models representing the combined knowledge drawn from pn and nr investigations . as mentioned , pn theory is based on an approximate expansion of einstein s equations in powers of velocity @xmath337 providing general relativistic corrections to the newtonian , small - velocity limit . this approach has a long history of success as the primary framework for deriving explicit predictions from general relativity for physical and astrophysical observations . the relatively small velocities involved in most observations have made low - order pn theory an ideal tool for testing general relativity as the standard model for gravitational physics in many contexts including earth - orbit experiments , solar - system dynamical observations , and in even precision pulsar timing observations of binary neutron star systems @xcite . gravitational - wave researchers have applied pn theory to represent general relativity - based signal expectations for the vast majority of gravitational - wave searches for anticipated observations of black - hole and neutron - star binary inspirals , and likewise in the development and planning for current and future gravitational - wave instruments . the early epoch of black - hole binary inspiral is also well - described by low - order pn theory . for near - circular inspirals , @xmath338 , where the binary separation @xmath44 is scaled by the total system mass @xmath9 . at large @xmath44 , the pn expansion provides excellent predictions , even at low order . as the black holes lose energy and sink closer together , the velocity grows , requiring higher - order pn terms for sufficiently accurate predictions . currently , pn predictions are available up to 3.5  pn order ( 2.5  pn for spin - terms ) . these higher - order expansions seem to provide pn predictions sufficiently accurate for the analysis of data from current instruments for all but the last several orbits . once the separation approaches , say , @xmath339 , the accuracy of the pn expansion diminishes . even at ( hypothetical ) arbitrarily high orders , the expansion may fail , if the expansion parameter @xmath340 exceeds the series radius of convergence . precisely when the pn approximation effectively fails depends on the details of the system being studied and the requisite accuracy , but generally , for the last orbits and merger pn theory is no longer reliable .    the strongest gravitational radiation is generated in the late stages of inspiral or merger where the internal consistency of the pn approach has been , at best , difficult to assess @xcite . numerical simulation can treat the late portions of the mergers , with practical resource limitations on the duration of the simulations , and consequently how far apart the black holes can be at the start of the simulation . would it be practical , however , to run numerical simulations long enough to `` overlap '' with the part of the problem treated successfully by pn methods , or would there be an intermediate region requiring yet another approach @xcite ?    the first numerical inspiral results were quickly compared against pn calculations @xcite . these first comparisons yielded promising indications that the gap between nr and pn could be bridged , but also made clear that pn results were not uniquely determined for the final part of the inspiral . numerical results can be verified by internal consistency studies , examining , for instance , the convergence of the results toward consistency with einstein s equations as the resolution is increased . external verification , however , would strengthen the case that these new results are indeed correct . even more importantly , these comparisons would allow an independent check of the late inspiral pn predictions which are not easily confirmed by self - consistency studies . roughly one year after the first robust numerical - relativity results , simulations lasting @xmath341 were conducted , and quantitatively compared with various pn approximations . quantitatively cross - checking waveform comparisons requires some care . how do the differences between pn and numerical results compare with intrinsic numerical error estimate how do different variants of pn predictions compare with the numerical results ? because each waveform comes with no meaningful absolute reference in time and phase , how can the freedom to offset these parameters be controlled for the comparison ? effective comparisons also require longer waveforms than those produced in the earliest simulations . the first comparison addressing these issues came at the end of 2006 @xcite . the study considered the case of equal - mass non - spinning black - hole mergers , comparing a 3.5  pn waveform with numerical results covering the last seven orbits . the compared waveforms are shown in fig . [ fig : pnprlwf ] . in the comparison the waveform phases agreed within 1 rad of phase drift , for a little over ten gravitational - wave cycles preceding the last orbit before merger , comparable to numerical error estimates . this result gave clear indication that pn waveforms could be accurate in the last orbits approaching merger , and that pn and nr , combined , could treat the complete waveform signal . not all equally valid variants ( approximants ) of the approximate pn waveforms agree this well with numerical results . to understand this , it is worthwhile to consider pn results in a little more detail . the pn approximation is formally understood as an expansion in powers of the speed of the merging black holes , @xmath340 . a concrete result of the theory will typically be a taylor expansion for a specific dependent variable , in terms of a chosen independent variable . the obvious choice , to express the waveform itself ( gravitational - wave strain ) as a function of time , would give a poor result since the sinusoidal shape of the waveform is difficult to approximate by a polynomial . pn waveform results are typically expressed as separate expansions for the orbital phase , and polarization component amplitudes . these are given as expansions in the time to merger @xmath342 or orbital frequency @xmath343 . the orbital phase information may also be expressed by an expansion for @xmath344 , referred to as the chirp - rate ; see @xcite for a review of the approach and a particular explicit waveform expansion . the pn information may also be encoded in an hamiltonian formulation description of the dynamics , which may be integrated numerically . researchers also choose between taylor series and other `` resummed '' expressions of the results based , for instance , on pad approximants @xcite .    in the comparisons described above , the pn waveform phasing was found to agree with numerical results derived from a taylor series expansion for chirprate @xmath345 expanded in powers of frequency . this is known as the taylort4 approximant in the language of the ground - based gravitational - wave community . later studies with longer and more accurate numerical simulations provided more precise tests of phases and amplitudes for a variety of pn approximants @xcite . @xcite found even closer agreement ( than found in @xcite ) with the taylort4 approximant , to within 0.05 rad over nearly 30 cycles before the orbital frequency reaches @xmath346 at roughly an orbit before merger .     up to @xmath347 , while the comparison in right panel extends this range to @xmath348 . ]    figure  [ fig : boylepncomp ] shows the results of four pn approximants with terms kept to various pn orders . the excellent agreement of the taylort4 approximant with numerical results is not matched by other equally consistent pn waveform approximants in the late portions of the waveforms . however , all four approximants agree to within about a radian if the comparison is cut off at waveform frequency @xmath349 , about five orbits before merger [ note that this is the frequency of the dominant @xmath350 mode , twice the orbital frequency @xmath351 . while most attention has so far been given to the nonspinning cases , some results exist for other regions of parameter space . for equal - mass configurations of black holes with spins aligned with the orbital angular momentum , the phase agreement of the taylort4 approximant is not generally best , and phase differences ( over ten cycles preceding @xmath352 ) are much larger than those seen in the nonspinning case , with differences near 1 rad or more for spins @xmath306 @xcite . pn spin effects are not yet derived beyond 2.5  pn order , possibly limiting the accuracy of the pn waveforms in these comparisons . some pn comparisons have also been examined for other configurations , including a precessing system @xcite and for eccentric configurations of equal - mass , non - spinning mergers @xcite .    in comparing fourier - domain pn approximants with numerical results for a set of non - spinning mergers with various mass ratios , @xcite discovered that the abrupt truncation of the fourier - domain waveforms , inducing gibbs oscillations in the time - domain , can approximately emulate the physical quasinormal ringdown that terminates the numerical - simulation waveforms . even better agreement is possible if @xmath254 , the symmetric mass ratio , is allowed to extend beyond its physical range ( to @xmath353 ) or if an extra pseudo-4pn order term is added . more recent studies have confirmed that these simple models are directly useful in some gravitational - wave observation applications @xcite . these results have motivated further development of phenomenological full - waveform models .      for the early part of the waveforms various pn approximants agree at 3.5  pn order , the best currently available , providing excellent approximations . at late times , the results of different approximants diverge , and the numerical - simulation results are the only way to accurately derive the predictions of general relativity . can an analytic waveform description be developed that simultaneously encodes the results of both complementary theoretical approaches ? the generally simple features of the numerically simulated merger waveforms raise hopes that fairly accurate approximate analytic waveform descriptions may be produced with simple dependence on the binary system parameters , providing an efficient means also of interpolating from a sparse sampling of the full parameter space for which numerical - simulation studies are completed . these waveform models would be computationally efficient ( compared to numerical simulations ) and may have applications in a broad class of observational data analysis applications . ajith and collaborators have developed a fourier - domain full - waveform model phenomenological approach resembling the treatment of @xcite , but with greater emphasis on matching the numerical waveforms @xcite . they begin combining information from the pn and nr - based waveforms by first stitching together time - series data for the dominant ( @xmath190,@xmath205 ) component waveforms , including a long pn precursor joined to a numerical merger waveform . fourier transforms of waveforms are then fit to a parametrized model of the waveform similar to waveform families previously applied in phenomenological treatments of purely pn waveforms @xcite . their waveform model for mergers of non - spinning binaries has the form @xmath354 where @xmath355 is a constant related to the masses , and @xmath356 is the distance to the source . the effective fourier amplitude @xmath357 is modeled piecewise , @xmath358 with distinct power - law segments before and after a merger frequency @xmath359 , and a lorentzian @xmath360 decay beyond the transition frequency @xmath361 demarking the ringdown of the final black hole . the fourier phase is represented by a single power - law expansion @xmath362 for each waveform , the free amplitude and phase coefficients are determined by a fit to the stitched nr - pn hybrid waveforms . it was found that the mass - ratio parameter dependence could be modeled by simple fits to quadratic functions of the symmetric mass ratio . these waveforms have been applied in gravitational - wave data analysis studies ( see sec .  [ sec : gwda : direct];@xcite ) . recently @xcite have developed a generalization of this model for non - precessing system of spinning black holes . stitching together pn and nr waveforms can be avoided , and closer contact with the basic physics maintained , by requiring direct pn consistency in the full waveform model . some pn approximants can resemble the numerical late - merger results . if such approximants can be tuned to agree with numerical results right up to merger by careful choice of adjustable parameters which only affect pn consistency beyond known pn order , the result would simultaneously encode pn and numerical - relativity results . the effective - one - body ( eob ) family of pn hamiltonian models has promise as a tunable model for encoding both pn and nr results @xcite . in the time domain , techniques have also been developed for extending these waveforms into the ringdown of the final black hole . analysis of the early numerical results in comparison with an untuned eob model gave promising indications that the waveforms could be closely approximated this way @xcite . with tuning , it appears that this construction can provide an analytic but potentially very accurate , approximation to the complete coalescence waveform .    in the eob model , the binary motion is recast as the motion of a single _ effective _ body of mass @xmath363 moving about a central potential , as is familiar from newtonian mechanics . in the general relativistic version of this framework , the effective body s motion follows a geodesic ( to 2  pn order ) around a modified version of a schwarzschild metric . the motion of the effective body is described by an effective hamiltonian , which , for systems of nonspinning black holes , may take the form    @xmath364 } \,.\end{aligned}\ ] ]    the expressions for @xmath365 and @xmath366 are chosen so that the pn - expansion of the hamiltonian is consistent with the results from pn theory to known order ( 3  pn for nonspinning black holes ) @xcite . the nonconservative contribution to the motion , arising from the loss of angular momentum to gravitational radiation , is encoded in an additional flux term entering as an external force in hamilton s equations , which is also constrained to be consistent with pn theory ( typically to 3.5pn order ) . hamilton s equations are integrated to derive the effective body s motion , and hence the black - hole trajectories . waveforms are constructed using pn extension of the quadrupole formula relating the motion of the black holes to the amplitude and phase of the multipolar radiation components . the last part of the radiation , arising after merger , is completed by continuously matching a sum of quasinormal ringdown modes to the waveforms truncated near the point of merger .    while consistent with pn theory , the formalism can be adjusted to also match the wave forms derived from numerical - relativity simulations . the flux function , as well as @xmath367 and @xmath368 , can be adjusted by the addition of higher - order pn terms to modify the strong - field dynamics to agree with nr without violating pn consistency . an early implementation of this approach , encoding the combined results pn and nr for nonspinning black - hole mergers was presented by @xcite . figure  [ fig : faithful ] shows their result comparing numerical and adjusted eob - model waveforms for a 4:1 mass - ratio merger . component of the gravitational - wave strain from the merger of binary with mass ratio 4:1 . the waveform is for an observer located an at inclination angle @xmath369 from the axis of rotation . [ fig : faithful],scaledwidth=45.0% ]    subsequent work involving more accurate numerical simulations , and more careful tuning of the eob model has improved eob model waveforms @xcite . recent comparisons of improved eob models with high - accuracy numerical results from @xcite yield differences comparable to numerical errors showing phase agreement within about @xmath370 rad through @xmath371 30 gravitational - wave cycles @xcite . recently the eob approach has been extended to match cases of non - precessing spinning black - hole mergers@xcite . there remain open questions in extending these waveform models to generic spin configurations , more extreme mass ratios , and full multipolar content . before a numerical simulation can commence , a numerical relativist requires a model for the initial configuration of the two black holes and the geometric fields which represent them . as discussed in sec . [ sec : comp : initialdata ] this involves not only solving the general relativistic constraint equations , but also providing some ansatz for the free data . in some cases , some information from pn theory may be applied to produce an initial data model that more precisely represents the desired physical configuration .    to understand how pn information can be applied , it will be useful to revisit the discussion of numerical initial data . typical simulations begin with a modeling ansatz , the brandt - brgmann model , for instance , which produces initial _ field _ data from a given set of particle - like parameters related to the black holes masses , spins , initial positions , and momenta . usually , the initial data models make no attempt to represent the gravitational radiation fields which would have been previously generated by the motion of the black holes , and do not utilize pn information in this reduction of the field degrees of freedom to particle parameters . information from pn theory is frequently applied in choosing the specific particle parameters which correspond to the sought - after simulation , particularly for circular inspiral configurations . work has also begun toward making richer use of pn information for improving the field ansatz . most of the simulations discussed have been designed to represent black holes in circularized orbits . before long - lasting numerical - simulation results were available , comparative studies of numerical initial data sets with pn - derived information provided a gauge of the results [ see for instance @xcite ] . without evolutions , these studies focused on theoretical properties of the black - hole configuration space such as the isco ( see sec . [ sssec : astro_grav : basic : bh ] ) . still subtle imbalances in the initial data , such as excess angular momentum for the chosen separation , can lead to eccentricity in the simulation which impacts the simulated radiation waveforms @xcite . the residual eccentricity can be reduced via a straightforward iterative procedure @xcite , but this requires repeated simulations lasting several orbits ; expensive in both time and computational resources . to minimize this eccentricity without resorting to iteration , several groups use trajectory information from pn theory in setting up the parameters for the numerical initial data . this approach enables simulations with eccentricities @xmath372 for equal - mass non - spinning mergers with initial separations @xmath373 @xcite . the technique is helpful in simulations with unequal masses and nonvanishing spins as well , though the residual eccentricity is generally larger @xcite . work is underway on techniques allowing the richer application of pn information from numerical initial data modeling . the common assumption of a conformally flat spatial metric in numerical initial data , disagrees with pn results at the 2  pn order @xcite , making this a likely leading source of modeling error which produces spurious initial transients in the simulations . though not yet as well - developed as the widely applied models , an alternative approach applies metric information from pn theory for a non - trivial initial conformal metric @xcite . with these techniques it is also possible to encode in the initial data information about the prior radiation generated by the system before the `` initial '' time of the simulation @xcite . we have noted several specific areas where research in numerical simulations makes contact with pn theory . these alone do not provide a full picture of the interplay between the two theoretical approaches . most broadly , pn theory provides a foundation for interpreting numerical results . numerical relativists draw widely from pn - based background knowledge of the black - hole binary . many important phenomena , such as the isco , spin - orbit coupling , spin precession , and ` gravitational rocket ' kicks were first studied in the pn approximation , providing a foundation for subsequent numerical studies . small examples can be noted throughout this review . in particular we point to the value of pn results in interpreting mergers of spinning black holes ( sec . [ sec : gws : spinning ] ) and kicks ( sec . [ sec : gws : uneq ] ) . analytic formulas for approximately expressing the final state of the black hole , mass , spin , and momentum have also drawn heavily from insights based on the pn treatment ( see sec .  [ sec : astro : recoil_formulae ] . ) several detectors are active or in their planning stages to detect the gravitational - wave signals from astrophysical processes . prominent among these are the ground - based interferometric detectors  ligo ( u.s.a . ) @xcite , geo ( germany ) , virgo ( italy ) , tama ( japan ) , and aigo ( australia )  sensitive to frequencies in the range @xmath374 , as well as next - generation instruments such as the einstein telescope @xcite . also planned to launch at the end of the next decade is the laser interferometer space antenna ( lisa ) , with complementary frequency sensitivities between @xmath375 . all of these instruments are subject to a variety of noise sources ; in the ground - based detectors , these sources will completely overwhelm the signal , unless filtered intelligently . figure  [ fig : ligo_virgo_noise ] shows the design sensitivity curve resulting from these disparate noise sources , for the ligo and virgo detectors .    , while virgo does a little better ( note that achieved sensitivities may not perfectly mirror the design curves ) . taken from @xcite . reproduced by permission of the institute of physics.,width=288 ]    burst analysis can extract many signals , but optimal results in gravitational - wave data analysis require _ matched filtering _ of the noisy input signal @xcite . fundamental to this approach is some measure of overlap between a physical signal @xmath376 and filtering template waveform @xmath377 . for this it is convenient to use a frequency - space inner product @xmath378 between signals , defined as @xcite @xmath379 , \label{eq : inner_product_da}\end{aligned}\ ] ] where @xmath380 and @xmath381 are the fourier transforms of the signals , and @xmath382 is the ( one - sided ) noise power spectral density of the detector we are interested in . this can be normalized to produce the _ match _ ( or _ overlap _ ) @xcite between two waveforms .    to use requires the availability of a set of waveform _ templates _  simple , few - parameter analytic model waveforms to filter against . even for a relatively simple system such as a black - hole binary , the combination of possible source and detector configurations leads to a 17-dimensional parameter space . some of these parameters are intrinsic to the source : total binary mass @xmath383 , symmetric mass ratio @xmath254 , spins @xmath384 , time to coalescence @xmath385 , eccentricity @xmath227 , and eccentric phase @xmath386 . the remaining parameters have to do with the relation between source and detector : ( luminosity ) distance to source @xmath387 , inclination @xmath388 , orbital phase @xmath389 , waveform polarization @xmath131 , and position of the binary on the detector s sky @xmath390 . restricting consideration to binaries that have circularized by the time they enter the detector s window @xcite allows us to neglect the eccentricity parameters @xmath391 . crucially also all observables have a simple dependence on the total ( redshifted ) mass , @xmath392 . similarly , the observed waveforms have a trivial dependence on the time to coalescence @xmath385 . as this means only one theoretical waveform has to be generated to cover all astrophysical masses and coalescence times , @xmath392 and @xmath385 are sometimes treated as extrinsic parameters instead . nevertheless , numerical relativists are left with a seven - dimensional parameter space to cover ( @xmath254 , @xmath393 , @xmath394 ) , with no obvious shortcuts to lighten the workload further . moreover , numerical methods will always struggle with finite accuracy . equipped with a parametrized set of template waveforms @xmath395 , we can test an incoming detector data stream @xmath396 , consisting of a ( possible ) `` exact '' gravitational wave @xmath397 and detector noise @xmath398 by calculating the _ signal - to - noise ratio _ ( snr ) @xcite : @xmath399 in the case where the signal @xmath397 in the data stream corresponds perfectly ( up to overall scaling ) with the model waveform @xmath400 for some value of the model parameters @xmath401 , we obtain the optimal snr : @xmath402 . we refer to this simply as the `` snr '' @xmath403 for the remainder of this section . expected snrs depend strongly on the detector in question . generally speaking , the ground - based detectors in operation ( ligo , geo , virgo , etc ) are expected to observe with only modest optimal snrs , while lisa is expected to observe mbh mergers with snrs of hundreds @xcite .    of course , real model waveforms will be inaccurate , due to incomplete knowledge of the underlying physics , or perhaps a desire for simplicity . in this case , the achieved snr can not be optimal . @xcite defined the _ fitting factor _ ( ff ) as the reduction in signal from such imperfect templates : @xmath404    a fitting factor of 1 means that the exact signal lies somewhere in the model space . however , for a particular detector , it is impossible to generate templates arbitrarily close together in any one parameter dimension . we can develop a sense of how closely spaced a set of templates must be so that any physical signal will be detected by at least one of them with high likelihood . @xcite developed a geometric picture of a template - space metric defining the distance between neighboring closely - spaced templates . with these tools , he was able to link the number of templates @xmath405 to the desired _ minimal match _ ( mm )  the worst value of the match between the signal and any template . the connection between these concepts is nicely laid out by @xcite . figure  [ fig : lindblom_fig3 ] , from @xcite , demonstrates how @xmath406 and @xmath407 combine to produce an effective mismatch @xmath408 . , model waveform with the same physical parameters @xmath409 , actual `` best fit '' model waveform @xmath410 , and template bank waveforms @xmath411 , @xmath412 . [ figure taken from @xcite],width=288 ]    the percentage of detected signals scales as the cube of the fitting factor . for this reason , detector analysts require rather stringent fitting factors from their template banks , e.g. , 97% ( implying detection of 90% of all signals present ) . for initial ligo , for instance , @xcite estimate that @xmath371 100 zero - spin templates ( produced from @xmath371 10 different numerical simulations ) would suffice to detect all nonspinning merging binaries . obviously we expect that the addition of spins of significant magnitude and arbitrary direction will greatly increase this number . detection is only a first step , although an important one . when signal strengths are high , as is expected with lisa for massive binaries , we might expect to be able to read off some of the system s intrinsic and extrinsic parameters . achieving this requires templates that accurately cover the parameter space of the binary , with a one - to - one relationship with the parameters of the binary  that is , that the best - fit model waveform @xmath410 and closest - parameters model waveform @xmath409 in fig . [ fig : lindblom_fig3 ] coincide closely ( within statistical errors ) . it also requires that we have an understanding of the probabilistic correlations between parameters .    for high-@xmath403 signals , the simplest way to estimate parameter uncertainties is using the _ fisher information matrix _ @xcite , which uses simple derivatives of the template waveform shape with respect to the physical parameters : @xmath413 in the high - snr limit , the parameters are assumed to have a multivariate gaussian distribution , and the fisher matrix is the inverse of the covariance matrix between parameters . more sophisticated approaches that do not assume high snr include markov - chain monte carlo methods @xcite . different requirements can be made on the quality of the developed templates , as described by @xcite . the first is that they be _ effectual_. roughly speaking , this means that the templates must be capable in the bulk of detecting gravitational - wave signals . in the language of matched filtering above , we demand that the fitting factor @xmath414 be extremely close to unity for any signal .    more stringently , we might demand that the templates be _ faithful _ : that is , that the template waveform corresponds to underlying source parameters that correspond closely to the parameters of the detected waveform source . in the language of @xcite , the model waveforms @xmath409 and @xmath410 coincide closely ( see fig . [ fig : lindblom_fig3 ] ) . this requirement will be crucial for identifying the astrophysical source of the radiation . traditionally , templates used in the search for comparable - mass black - hole binaries have depended wholly on pn theory . as detection strategies are most sensitive to phase discrepancies between template and observed signal ( see below ) , the determination of waveform phase to high accuracy has been crucial . modern pn templates include phase terms up to 3.5pn order for nonspinning binaries , with 2.5  pn order spin corrections . many templates use this high - accuracy phase with the leading - order quadrupole amplitude to produce _ restricted _ templates ; however , amplitude corrections are known to 2.5  pn order ( higher for certain low-@xmath415 modes ) @xcite . while such pn - driven templates offer excellent phase accuracy during the long , slow inspiral of the binary , their usefulness becomes questionable as we approach the last pre - merger orbits of the binary . we may ask , then , what does our new numerical insight into the final moments of merger bring us ?      once full numerical waveforms became available , several groups used the results to test older predictions of the effect of the merger segment on observability . for the equal - mass nonspinning case , the numerical merger gave results consistent with @xcite for initial ligo , though with a smaller merger snr @xcite . the boost in snr from merger is significantly greater for advanced ligo and lisa , as shown in fig . [ fig : ligosnr ] , from @xcite .     before merger ) is used ; the dotted curve uses the signal between this time and when final plunge commences ( @xmath416 to @xmath417 before merger ; the thin solid curve uses just the merger waveform ( starting @xmath417 before merger ) . the thick solid curve is the combined result of using the entire waveform . the snr , and hence distance reach , of the detector is clearly greatly enhanced by the final merger portion.,title="fig:",width=336 ]   before merger ) is used ; the dotted curve uses the signal between this time and when final plunge commences ( @xmath416 to @xmath417 before merger ; the thin solid curve uses just the merger waveform ( starting @xmath417 before merger ) . the thick solid curve is the combined result of using the entire waveform . the snr , and hence distance reach , of the detector is clearly greatly enhanced by the final merger portion.,title="fig:",width=336 ]    the relevance of the final plunge and merger to the overall snr of the binary depends on the total mass of the binary and where it falls in the window of the detector . figure  [ fig : lisahchar ] , adapted from @xcite , shows this for lisa , plotting the `` characteristic signal strain '' @xmath418 , the late - merger and ringdown signal will fall outside lisa s sensitivity band . for such low masses , inspiral - only waveforms should be adequate for data - analysis purposes . curve , we mark times before the peak amplitude ( circle )  one hour ( square ) , one day ( diamond ) and one month ( triangle).,width=336 ]    the ligo and virgo detectors are sensitive to frequencies above @xmath419 ( see fig . [ fig : ligo_virgo_noise ] ) . the merger of binary systems is marked by a chirp gravitational - wave signal , whose monotonically increasing frequency saturates at the dominant quasinormal mode , which depends on the post - merger hole s mass and spin . for a nonspinning binary of total mass @xmath9 , this @xmath420 . this means that a merging binary of mass greater than @xmath421 will never be seen by ligo . in contrast , most post - newtonian templates stop at @xmath422 , the frequency at the innermost stable circular orbit ( well defined only for test particles ; see sec . [ sssec : astro_grav : basic : bh ] ) . this is @xmath423 for an incomplete pn waveform to be useful , we want the missing merger - ringdown section ( i.e. the part of the signal with @xmath424 ) to contribute as little as possible to the total snr . for instance , if we take the high end of ligo s sensitivity window to be @xmath425 , this is @xmath422 for a binary of total mass @xmath426 . then , for binaries with @xmath427 , @xmath428 will fall within the visible window ; that is , ligo will see the final stages of the binary merger , where no acceptable pn waveform is available . recent work by @xcite supplied a more precise answer : inspiral - only pn templates can be trusted for all ligo configurations for systems with @xmath429 . above this critical mass , there is a gap in reliable information , which must be filled by numerical results . figure  [ fig : panetal_snrdist ] , adapted from @xcite , shows the importance of the last stages of merger in initial ligo as a function of mass range . the thin blue curves are hybrid nr / eob waveforms , while the thick red waveforms are generated from these by `` whitening '' : that is , they have been fourier - transformed to the frequency domain , then rescaled by @xmath430 , and finally re - transformed to the time domain . the whitened amplitude in a segment is proportional to the contribution of that segment to the total snr ; each marked segment accounts for 10% of the total . as the total system mass @xmath9 increases , the whitened amplitude becomes more bulked toward the merger time . for the largest total mass shown , @xmath431 , more than 90% of the signal power comes from the last cycle + merger + ringdown . increases , the final merger and ringdown account for a greater fraction of the total power . adapted from @xcite.,width=336 ]    in sec . [ sec : gws : longer_wfs ] , we referenced the recent `` samurai '' project @xcite , which establishes the consistency of all of the `` long '' ( pre - merger duration @xmath432 ) equal - mass nonspinning waveforms within their stated numerical accuracy . as well as direct comparisons of phase and amplitude errors over time , they conducted mismatch tests in the regime of the ligo and virgo ground - based detectors . they demonstrated that quadrupole waveforms ( @xmath433 ) from the five numerical codes have tiny mismatches ( @xmath434 ) for binary masses above @xmath435 with `` enhanced '' ligo , and masses above @xmath436 with advanced ligo , virgo and advanced virgo . all waveforms would be indistinguishable for snrs below 14 , entirely reasonable for the current generation of ground - based detectors ; this is shown for the advanced ligo detector in fig .  [ fig : samurai_measuredaligo ] . below which different `` long '' numerical waveforms are indistinguishable for the advanced ligo detector . even the earliest , lowest - accuracy published long waveforms are indistinguishable for snr levels below 14 , as demonstrated in this figure from the samurai paper @xcite.,width=336 ]    studies using incomplete ( non - merger ) pn waveforms have shown that introducing fuller harmonics of spinning and precessing binary systems greatly enhances parameter estimation , removing some parameter correlations and hence reducing certain errors by several orders of magnitude ( though other parameters are significantly less improved ) @xcite . it is expected that the introduction of full numerical waveforms for the merger portion will lead to further improvements , at least in regions of the detector s window where the binary merger is visible . a recent study by @xcite used phenomenological merger - inspiral - ringdown templates ( discussed in sec . [ sec : gwda : templates ] ) to assess these improvements for advanced ligo and advanced virgo . they showed that parameter - estimate accuracy is improved substantially compared to inspiral - only waveforms for system masses @xmath437 ; in particular , the average sky - position error for an @xmath438 equal - mass , nonspinning binary at a luminosity distance of 1 gpc drops to about one - tenth of a square degree . since they only uses the dominant mode for this investigation , there will be systematic errors as well . numerical verification of the extent of parameter estimation improvements for lisa has led to differing estimates @xcite . in particular , @xcite find that 50% of injected signals could be located on the detector sky within an error box of 3 arc min , while @xmath385 could be measured to less than a second and @xmath387 within 1.5% ( ignoring issues of weak lensing ) . they attributed these impressive results to a combination of very high snr ( between 900 and 9000 ) and the use of higher multipoles of the radiation . in contrast , @xcite found smaller ( by about one order of magnitude ) improvements over inspiral - only waveforms . this may in part be attributable to the use of different mass ratios and different - length waveforms by the two groups , but has not been fully resolved . one of the eventual aims of numerical - relativity simulations of binary mergers is the production of a set of gravitational - wave `` templates '' that will cover the many - dimensional space of astrophysical parameters . numerical merger simulations are still orders of magnitude too slow to run on the fly to filter incoming detector data streams ; thus there is a need to develop analytic waveform expressions that encode the numerics . it has been demonstrated in the physical regime of nonspinning binary mergers @xcite that inspiral - only pn - based templates become inconsistent with each other ( and with full waveforms ) at total binary masses @xmath439 for initial and advanced ligo configurations ; above this mass , full templates including the numerical - simulation - based understanding of merger / ringdown become necessary . one such set of templates incorporating numerical - relativity data is due to @xcite . these were extensions of the stationary - phase approximation ( spa ) templates used in ligo / virgo data analysis , with input from equal - mass numerical waveforms supplied by frans pretorius @xcite , as well as equal- and unequal - mass numerical waveforms from the goddard group . these templates are effectual  they achieved fitting factors of @xmath440  but they did this by extending source parameters into unphysical space .    to date a handful of useful faithful full - waveform template banks has been produced , aimed at covering all mass ratios for nonspinning binaries . the `` phenomenological '' templates of @xcite are simple three - segment curves in frequency space , with matching parameters tuned by numerical data generated by the ` bam ` @xcite and ` ccatie ` @xcite codes . the templates of @xcite , coded in the lsc algorithm library as `` eobnr '' , are an extension of effective - one - body waveforms , with a single `` pseudo-4pn '' parameter tuned by numerical data ( initially from the goddard hahndol code , and later using waveforms from the caltech - cornell group . both the phenomenological and eobnr template banks are discussed in sec . [ sec : nrpn : analytic ] . both template banks are faithful in the sense described above , in the restricted @xmath441 parameter space , and both have been used in data - analysis injection tests @xcite . currently , eobnr templates are being used both for injection and filtering in the high - mass region of ligo test analysis , and also for injection for the ringdown band . they were also both used in the matched - filter analyses of the ninja project ( see sec . [ sec : gwda : inject ] ) , performing as well as the inspiral - only templates in detection , and considerably better in the limited parameter estimation attempted . figure  [ fig : ninja_eobnr_parest ] , from @xcite , shows the accuracy with which the total mass @xmath9 and time of merger were extracted from all detected injections in ninja , when the eobnr templates were used . results for the phenomenological templates were similar for @xmath9 , but performed slightly less well for time of merger , as it is not an explicit part of this model ( given the relatively small number of samples in the ninja tests , such minor differences may not be significant ) . , for all detected signals , while the lower shows the error in merger time . adapted from @xcite . reproduced by permission of the institute of physics.,width=288 ]    , for all detected signals , while the lower shows the error in merger time . adapted from @xcite . reproduced by permission of the institute of physics.,width=288 ]    as significant spin is expected in astrophysical black holes at all scales , it is important to develop templates that encode spin also . @xcite demonstrated that matched filtering with nonspinning merger waveforms does a good job in detecting waveforms from binaries with significant spins . as their spinning waveforms come from aligned - spin mergers where the total spin is zero ( one hole s spin is aligned with the orbital angular momentum , while the other s is anti - aligned ) , their configurations will lack some prominent spin - related effects , including spin - orbit effects ( pulling in or pushing out the isco of the orbit ) and precession of the orbital plane . their conclusions , then , are likely to be optimistic , and more generic configurations can be expected to do much worse with nonspinning templates . the phenomenological templates of ajith _ have recently been extended to include non - precessing spinning systems @xcite . these have been shown to be both effectual and faithful in searches over non - precessing binary hybrid waveforms with ligo , where earlier zero - spin templates failed badly . for maximum simplicity , these new phenomenological template families are ultimately parametrized by three numbers : the total mass @xmath9 , the symmetric mass ratio @xmath254 , and a single effective spin @xmath442 where @xmath443 . an obvious application of numerical waveforms for data analysis is in directly testing analysis algorithms . as current templates are overwhelmingly based on pn information  either extended heuristically to include merger and ringdown , or terminating before merger  so current detection and parameter estimation algorithms are also based on incomplete information . recently , a multi - group collaboration ( the ninja project ) has performed the first direct insertion of explicit numerical waveforms into a simulated ligo datastream , and investigated detection efficiency and systematic errors @xcite . essentially all active numerical - relativity groups contributed waveforms ; these were scaled randomly to represent different - mass sources , resulting in 126 different injections . nine data - analysis groups analyzed the post - injection data stream , with a variety of methods . these fell into three main groups : matched filtering , burst analysis , and bayesian techniques . for the matched filtering , some groups used standard `` taylorf2 '' inspiral - only templates generated by pn theory , some used ringdown - only templates , and some used full inspiral - merger - ringdown templates , such as the eobnr and phenomenological waveforms discussed below and in sec . [ sec : nrpn : analytic ] . both the taylorf2 and the full template waveforms yielded quite high detection rates  around 80 of the 126 injected signals in triple coincidence between the ligo detectors  while the ringdown templates performed more poorly , detecting only 45 signals in triple coincidence . however for parameter estimation of the total mass @xmath9 , most of the taylorf2 estimates have an error of 40% or more , significantly more than for the eobnr and phenomenological estimates . figure  [ fig : ninja_eobnr_parest ] shows errors in detected mass ( upper panel ) and coalescence time ( lower panel ) when using the eobnr templates . while mass errors are still quite large , between @xmath444 and @xmath445 for the bulk of detected signals , they are significantly less biased than those from the taylorf2 templates ( compare with fig .  8 of @xcite ) . estimation of the coalescence time @xmath385 was better for the eobnr and phenomenological template banks as well , as might be expected ( since taylorf2 templates must necessarily cut off before coalescence ) . however , the eobnr and phenomenological templates typically covered a larger mass range as well ( up to @xmath446 compared to @xmath447 for taylorf2 ) , which will bias the results when the injected signals had a large mass . burst analysis is a more abstract approach to the problem , using partial information about signals , such as dominant frequency and approximate duration . in ninja , two burst techniques were used  the `` q pipeline '' and the hilbert - huang transform ( hht ) . the q pipeline fits signals to a sine - gaussian model parametrized by a central frequency and a @xmath24 factor ( @xmath448 , the number of oscillations under the gaussian ) . at the single - detector level , the detection rate is comparable to that of the eobnr and phenomenological matched - filter searches . parameter estimation is limited to estimates of the peak power time ( close to @xmath385 ) and peak frequency ; for these parameters , performance was comparable to the ringdown matched - filter search . the hilbert - huang transform decomposes input data into `` intrinsic mode functions '' , each characterized by a single frequency scale . applied to the signal , the hht generates a high - resolution time - frequency map of the data @xcite . the hht applied to the ninja data found @xmath449 coincident events , competitive with the best matched - filter results ; however , no parameter estimation is currently possible with this method . bayesian analysis attempts to reconstruct the posterior probability density function of a parametrized signal based on the observed signal . two variant bayesian approaches were taken in ninja analysis , involving different monte carlo methods  markov chain monte carlo ( mcmc ) and nested sampling . in both cases , a parametrized waveform model is needed . for the mcmc study , a restricted pn waveform was used , only 1.5  pn in phase , but including leading spin effects . detection of injected signals performed well , but the high mass of the injections meant that the bulk of the snr occurred in the merger - ringdown phase ; as a result the inspiral - only model masses merger time was significantly biased . the nested sampling effort instead uses the taylorf2 templates ( restricted amplitude , 2.0  pn in phase ) as well as phenomenological templates . in detection , the phenomenological templates significantly out - perform the taylorf2 model ; in parameter estimation , it achieves results consistent with those of the matched - filter analysis . overall , current results are encouraging , with high detection rates of injected waveforms , and good , if limited , accuracy in estimated parameters when attempted with faithful template banks . however , the project was hampered by the limited length of the contributed waveforms and unrealistically stationary ligo and virgo noise profiles . follow - on work in this area will demand longer numerical waveforms that span the ligo sensitivity window , or more likely , a systematic blending of long pn inspiral waveforms to late - inspiral / merger numerical waveforms , but also may be extended to other gravitational - wave sources . the recent successes in modeling of binary black - hole mergers have captured the attention of astronomers and astrophysicists . while black - hole binaries have been studied both observationally and theoretically for many years , most efforts have been primarily within the framework of newtonian gravity . the new developments in numerical relativity now supply the missing piece : the effects of the final merger in the strong - field general relativistic regime . in this section we highlight three key areas of current interest : recoiling black holes , the spin of the merged black hole , and electromagnetic signatures from the final merger . one of the more dramatic implications of gravitational - wave computations for astrophysical observations is the possibility that gravitational radiation - induced recoil will eject black holes from their host galaxies . the largest recoil velocities predicted by numerical relativity exceed the escape velocities of many galaxies ( @xmath450 ) . calculating the probability of this kind of event requires a detailed understanding of the recoil on mass ratios , spin magnitudes , and spin orientations , together with some expectations about the distributions of these parameters . it would be useful to have in hand a simple , analytic formula expressing the relevant dependencies of the recoil velocity . the highly nonlinear , strong - field interactions of merger dynamics determine the final recoil velocity of the merged remnant . numerical simulations are required to accurately compute such results . however , in attempting to construct an ansatz for a phenomenological formula , one might initially assume a form consistent with pn predictions . this will suggest the functional dependence on mass and spin , and their leading powers . for example , pn analysis @xcite suggests that the recoil due to unequal masses may be proportional to @xmath451 ( the  fitchett formula " ) , where @xmath254 is the symmetric mass ratio . meanwhile , to leading pn orders , the spins should contribute through components of @xmath452 and @xmath453 @xcite . such an ansatz can be further supported and constrained by symmetry arguments . an economical approach is to consider the recoil velocity as a taylor expansion in all six spin components ( three for each black hole ) , where the coefficients are functions of mass and initial separation . rotation , parity , and exchange symmetries of the binary system then impose relationships between these coefficients @xcite . by such arguments , for example , the component of the recoil velocity parallel to the orbital angular momentum can be shown to be proportional , to leading order in spin , to the dot product of the coordinate separation vector between the black holes and some linear combination of the spins @xcite . unknown coefficients in the ansatz can then be fit to numerical results , and correction terms can be added as needed for numerical agreement . with these considerations , a tentative formula has been taking shape through the combined efforts of the numerical - relativity community , of the following form @xcite : @xmath454 \nonumber\\                & + & \frac{k_s(q-1)\eta^2}{(q+1)^3 } \nonumber\\                & \times&\left[q^2\alpha^{\perp}_1 \cos(\phi_1-\phi_1)+\alpha^{\perp}_2\cos(\phi_2-\phi_2)\right ] \label{eq : v_parallel}\end{aligned}\ ] ] where @xmath455 is the contribution due to mass asymmetry , @xmath456 is the contribution due to spin that yields a kick perpendicular to the orbital angular momentum , @xmath457 is the contribution due to spin that yields a kick parallel to the orbital angular momentum , @xmath458 is the projection of the dimensionless spin vector @xmath459 of black hole @xmath134 along the orbital angular momentum , @xmath460 is the magnitude of its projection @xmath461 into the orbital plane , @xmath462 refers to the angle made by @xmath461 with respect to some reference angle in the orbital plane , and @xmath463 and @xmath464 are constants for a given mass ratio and initial separation . here the spins are considered to have been measured at some point before merger or ideally , arbitrarily close to merger . @xmath463 and @xmath464 encode the amount of precession of each spin before merger . one can then evaluate the formula over ranges of expected mass ratios and spin magnitudes , and for all spin orientations , to compute the probability of a given recoil speed . for example , for mass ratios between 1 and 10 , and spin magnitudes of 0.9 , the probability of exceeding @xmath59 is predicted to be @xmath216 . studies of various model - dependent speed probabilities are given by @xcite and @xcite .      among the significant astrophysical consequences of gravitational - radiation recoil , growth rates of black holes can be affected . for example , the recoil of massive black holes imparting less than the escape velocity of their host galaxies tends to have regulatory effects on the amount of mass accreted @xcite . recoil velocities that exceed galactic escape velocities can impact massive black - hole growth in a different way , by reducing the chances of subsequent mergers @xcite . this in turn may modestly reduce the rate of coalescence events observable by lisa @xcite . the recoil might also have more directly observable consequences . some quasars are thought to originate from the coalescence of massive black holes during galactic mergers . quasars ejected from their host galaxies might therefore be expected . a study by @xcite found no evidence for such events , indicating that they are rare . however , recent observations of two rapidly moving extragalactic quasars @xcite are strong candidates for such recoiled black holes @xcite . the final spin of the merged remnant of a binary is also of astrophysical interest . it would be useful to know the probability of various spin magnitudes in constructing matched filtering templates and estimating gravitational signal detectability for ligo and lisa @xcite . the ability to predict a final spin given initial binary parameters also implies that observation of the spin of a black hole may help understand its origin . in particular , an understanding of the relationship between the final spin of a black hole and the orbital angular momentum of its binary precursor may help explain the formation of x - shaped jets @xcite .    as with the kick , it would be advantageous to construct a simple analytic formula with which to estimate the spin . observations from numerical simulations suggest this might be possible . for example , in the case of nonspinning , equal - mass binaries undergoing circular inspiral , the final spin is insensitive to the initial separation , being determined by the terminal dynamics of the merger . such evolution results in a universal final spin of @xmath465 @xcite . additionally , for unequal - mass , nonspinning binaries , the final spin scales roughly linearly with symmetric mass ratio @xcite . various approaches to constructing an ansatz for the final spin have been proposed . for aligned - spin , equal - mass mergers , @xcite produced a simple formula quadratic in the ( common ) initial spin ; this satisfied the `` cosmic censorship '' hypothesis : @xmath466 . for nonspinning mergers , @xcite simply assumed @xmath467 , for fitting parameters @xmath13 and @xmath468 , which agreed with data from nonspinning binaries reasonably well . @xcite arrived at a more generic ansatz , for initial black holes of arbitrary spins , using the pn approximation . @xcite proposed a different expression for generic binaries , using a set of assumptions about the approximate conservation of the magnitudes and relative angles of certain angular momentum vectors . several other formulas have been suggested @xcite , each in good agreement with some subset of available numerical data . however , broad agreement on an analytic formula for the final spin from generically precessing binaries has yet to be achieved . we have seen above that black - hole mergers are expected to be `` loud , '' since they produce strong gravitational - wave signals . but will they also be `` bright ? '' that is , will there be an accompanying display of photons , detectable by telescopes in any frequency range of the full electromagnetic spectrum ? the answer to this question depends critically on the amount and distribution of gas and magnetic fields surrounding the merging binary . for stellar black - hole binaries , and intermediate - mass black - hole binaries in stellar clusters , any matter in an accretion disk would be consumed by one or both black holes and disappear relatively quickly , making an electromagnetic counterpart highly unlikely . massive black - hole binaries formed from galaxy mergers , however , present a very different situation . in a gas - rich or `` wet '' merger , there is likely enough gas available to feed accretion disks around each black hole that eventually evolve into a circumbinary disk ; this provides a source of gas and magnetic fields that could generate detectable electromagnetic emission . even in the case of gas - poor or `` dry '' mergers , the thin hot gas present is such ( elliptical ) galaxies might be sufficient to produce an electromagnetic signature . such electromagnetic signatures would be valuable to astrophysics . identification of the source on the sky would help confirm and characterize the merger , and probe accretion physics . a measurement of its redshift using electromagnetic radiation , taken together with the determination of luminosity distance using gravitational waves observed by lisa @xcite , would provide an independent calibration of the distance scale and a precise probe of cosmology including the nature of the mysterious dark energy @xcite . differences in arrival times of the electromagnetic and gravitational - wave signatures could also test fundamental principles such as the relative propagation speed of photons and gravitational waves .    with the recent successes in numerical - relativity modeling of black - hole mergers , we find considerable interest in understanding possible electromagnetic signals from these events . most work to date has focused on mechanisms that could produce emission from a surrounding accretion disk , including signals induced by a recoiling merged black hole encountering the disk @xcite . moreover , explorations of possible electromagnetic signals that could arise in the dynamic spacetime near the binary during its last few orbits and merger are just now beginning . @xcite recently studied the effects of a merging black - hole binary on a surrounding magnetic field . the equal - mass , nonspinning holes start out with separation @xmath469 on quasicircular orbits somewhat outside the isco . the magnetic field is initially poloidal and assumed to be generated by currents in a distant circumbinary disk located at @xmath470 . the electric field around the binary is initially zero . there is no matter near the black holes , as predicted when the circumbinary disk is thin @xcite . they solved the coupled einstein - maxwell equations for a binary evolving in presence of externally sourced electromagnetic fields . as the system evolves , the inspiralling black holes stir up the fields .     before the merger , from @xcite . the electric field lines are twisted around the black holes , while the magnetic field lines are mainly aligned with the @xmath34 axis . ]    , except at a later time , @xmath471 before merger . ] figures  [ fig : em_early ] and  [ fig : em_later ] show the electric and magnetic field lines at @xmath472 and @xmath471 before the merger , respectively . the magnetic fields are mostly aligned with the @xmath34 axis , while the electric fields are twisted around the black holes . the binary dynamics induce oscillations in the electromagnetic energy flux , with a period approximately half of the dominant ( quadrupole ) gravitational - wave signal . the energy in the electromagnetic field is also enhanced gradually . @xcite have certainly made an interesting start on this important problem . can this scenario generate detectable electromagnetic emission ? the answer awaits more realistic simulations , including spinning black holes . other astrophysical scenarios allow both matter and magnetic fields in the vicinity of the binary close to merger . in this case , calculation of possible electromagnetic signatures requires solving the equations of general relativistic magnetohydrodynamics in a dynamically evolving spacetime governed by the einstein equations . @xcite took a step towards solving this problem by mapping the flow of pressureless matter , modeled as non - interacting point particles , in the dynamical spacetime around the merging black holes . they started with a distribution of particles around the black - hole binary , and then evolved the binary using numerical relativity while tracing the motion of the particles along geodesics as the binary evolves . to estimate the energetics of the flow , they detected `` collisions '' by looking for particles within a small distance ( typically , @xmath473 ) of each other and then calculating the lorentz factors @xmath474 in the center - of - mass frame of each collision . they began with equal - mass inspiralling black - hole binaries @xmath475 orbits before merger , and considered both black holes to have either zero spin , or @xmath476 with spin vectors aligned with the binary orbital angular momentum . a single black hole of mass @xmath9 is also evolved as a control case . approximately 75,000 geodesic particles are then initially distributed uniformly throughout a solid annulus centered on the binary and having inner radius @xmath477 , outer radius @xmath478 , and vertical full thickness @xmath479 . particles within the inner radius are excluded , to avoid transient signatures from particles initially near the horizons . note that these circumbinary disks are geometrically thick , and thus potentially have high enough inward radial speeds to keep up with the shrinking separation of the inspiralling binary , providing a source of gas near the black holes during their merger @xcite . @xcite studied two initial velocity configurations for the particles . in the `` orbital '' configuration , the initial velocities are randomly distributed around a tangential velocity @xmath480 that would give a circular orbit in a schwarzschild spacetime of mass @xmath9 , resulting in a scale height of @xmath481 . in contrast , the `` isotropic '' configuration is an extreme ( that is , very hot ) case in which the particles only have random velocities , with each component sampled from a gaussian distribution of standard deviation @xmath482 .    for the orbital initial velocities , the particles remain mostly in a rotational configuration ; some particles do enter the region around the black holes , but are soon ejected by a gravitational slingshot . the collisions typically occur at relatively low velocities with lorentz factors @xmath483 . by the time of merger , there are essentially no particles near the black holes , and thus there is no energetic signature of the merger . the situation is quite different with the isotropic initial conditions , which provide a continual influx of particles throughout the evolution . moreover , there is a clear signature of the merger visible in the maximum lorentz factor , which takes on values @xmath484 during the inspiral and spikes up to @xmath485 just before merger for the nonspinning case . for rotating black holes , @xmath486 during the inspiral , spiking to @xmath487 just before merger . in addition , gravitational torques from the merging binary effectively stir the particle distribution , leading to high - velocity outflows with particles reaching @xmath488 or more . @xcite point out that , in realistic astrophysical disks , viscosity would cause angular momentum transport , bringing material inward towards the merging black holes . this could produce a scenario between the two extremes they studied , with a clear merger signal . more recent work by @xcite , using hydrodynamical simulations of gas around black - hole binaries , provided strong support in this regard . further confirmation of this interesting possibility , along with other mechanisms involving general relativistic magnetohydrodynamics , awaits more detailed simulations with gas and magnetic fields , currently in development . five years ago , we did not know what obstacles might conspire to prevent the effective application of numerical - relativity simulations in deriving predictions of general relativity to address questions of astrophysical black - hole binaries . up to that time , the numerical relativity researchers addressed their research primarily to questions of theoretical and computational physics . then , quite suddenly , the last obstacles vanished . as reported , numerical relativity can now be applied to understand general relativity s predictions of strong - field gravitational physics and to begin addressing questions of astrophysics and gravitational - wave data analysis . we expect future work in this field to be guided not by internal problems in gravitational theory , but by questions of astrophysics , and other areas where strong - field gravitational theory applies . these questions will motivate both continued , more detailed , investigation of some of the phenomena which have already been revealed , and the development of new capabilities to bring simulation results to bear on new questions . though the work has begun in applying numerical - relativity results to gravitational - wave observations , much more work is still needed . it is now clear that the signal - to - noise ( snr ) ratio from the last moments of the merger will dominate over the inspiral signal in many potential observations . though numerical - relativity simulations have shed considerable light on the basic features of the burst of gravitational radiation that completes the predicted signals , applications in data analysis will require _ comprehensive _ quantitative knowledge of the merger signals generated over the full parameter space of mergers . while we discussed examples of mergers that sample the binary black - hole parameter space along what may be its principal axes , a full quantification of the signal space will require a long systematic program of investigation . even with the aid of empirical models for encoding numerical results , many simulations must be conducted to qualify possible models . initial models suitable for detection and ( perhaps ) parameter estimation of low - snr signals , will require much further development for application to the high - snr observations expected with future instrument upgrades , such as advanced ligo , and future instruments such as lisa . considerably more simulations of high quality will be required to achieve these goals . in addition to understanding the signals generated by mergers through the bulk of parameter space in greater depth , relativists must also extend the region of parameter space covered by simulations . typical current simulations are limited in the range of parameters that can be practically covered . currently popular initial data models , particularly those assuming conformal flatness , are limited to @xmath489 @xcite . alternatives more suited to rapidly spinning kerr black holes are less well developed in other ways . while astrophysical limits on spin magnitudes remain unclear , it is likely that larger spins will have to be treated with possibly novel initial data models . other , related challenges in the simulations may also exist . current simulations are also restricted to comparable mass ratios . simulations beyond @xmath490 require somewhat heroic investments of computational resources with present techniques . astrophysical mergers may occur , however , over a broad range of mass ratios . in the limit of extremely large mass ratios , approximation techniques are possible that treat the motion of the smaller object as a perturbed geodesic in the spacetime of the larger black hole . such analytic methods look promising in this limit , where numerical simulations are less practical . however , there is a rather large middle ground , say @xmath491 , where new techniques may be required . one problem is that the small spatial scale of the smaller black hole requires extremely small timesteps for stable explicit time integrations , even though there is nothing interesting happening on these timescales . alternatives such as implicit schemes @xcite may open the door to a broader range of applications . most of the simulations discussed also focus on circular , or nearly circular orbits . while this is clearly an important portion of the black - hole - binary population , scenarios have been proposed in which one black hole captures another from nearly parabolic initial encounters . understanding the potential signals from such systems may require new techniques appropriate for long periods of slow , effectively newtonian evolution , punctuated by brief periods of strong gravitational interaction . as discussed , the numerical discovery of strong gravitational - wave recoils in asymmetric mergers of spinning black holes has had significant impact in areas of astrophysics beyond direct gravitational - wave observations . these results have spurred excitement in the black - hole astronomy community about other possible applications of numerical relativity . since other areas of astronomy are based on ( primarily ) electromagnetic observations , these interactions naturally lead to applications involving the interactions of black holes and electromagnetically visible matter . there are many relevant phenomena of interest : accretion disks around black holes , possibly disturbed by gravitational recoil , black - hole  neutron - star binaries , neutron - star  neutron - star binaries , jets from active galactic nuclei , quasars , and the mysterious origins of gamma ray bursts .    a key direction for future work has been adding physics beyond the purely gravitational simulations we have reviewed here . as discussed , efforts are currently underway to use and develop magnetohydrodynamics in a dynamical spacetime . as such simulations are applied with increasing realism to more astrophysical scenarios , refinement of the techniques will be required . adequately resolving shocks and turbulence interacting with magnetic fields , achieving robust accuracy in the presence of very strong magnetic fields , and enforcing the vanishing divergence of the magnetic fields , all in curved spacetime , are among the challenges that need to be met . astrophysical mergers typically involve nearly circular orbits , or in extreme scenarios , initially parabolic encounters . numerical simulations can be applied to other sorts of black - hole interactions as well .    for high - velocity black - hole encounters , one expects a phase transition dividing the spacetimes in which the black holes merge from those in which they pass each other hyperbolically , never to approach each other again : see fig . [ fig : crit_traj ] . numerical simulation studies are beginning to elucidate critical behavior near this phase transition @xcite . these fascinating phenomena at the extreme limit of strong gravitational dynamics have been discussed in the context of upcoming experimental measurements . in the trans - planckian energy limit , some particle collisions may be describable by classical black - hole dynamics @xcite . in highly speculative theories involving large extra spatial dimensions , it is conceivable that tev scale physics may be sufficient to produce such black - hole - like collisions @xcite . these possibilities motivate the first numerical - simulation studies of ultra - relativistic black - hole collisions @xcite . einstein gave us general relativity nearly 100 years ago . over the past century this theory , our standard model of gravitational physics , has passed experimental and observational tests over ranging from laboratory - scale physics experiments and solar - system tests to observations of compact astronomical objects and cosmology @xcite . these measurements have , so far , provided no need for a refinement of einstein s theory of gravity [ though some models to explain cosmic acceleration do involve alternative gravity @xcite ] . numerical simulations are now revealing the detailed predictions of einstein s theory for the final merger dynamics of a black - hole binary , and its record in the emitted gravitational radiation . gravitational - wave observations of these events will expose such gravitational phenomena to measurements in a strong - gravity regime far beyond anything which has previously been tested . these include the first tests of higher - order terms in the pn expansion and , indeed , of the physical predictions which numerical relativity is just now revealing . the possibility that general relativity will find its limits in these observations motivates a better understanding of what waveforms might be predicted by alternative theories of gravity @xcite . alternative - gravity models of black - hole mergers would also likely require numerical simulation to derive waveform predictions , but little work has been done in this direction so far [ see one related example from @xcite ] . if einstein s theory proves correct in predicting the detailed physics revealed by numerical relativity , it would stand as a truly incredible achievement of scientific induction , divining the details of phenomena vastly removed from physical observations on which the theory was founded . if not , then these observations , together with a confident understanding of einstein s predictions founded in numerical relativity simulations , may indeed lay the foundation for the next theory of gravity . this review draws on the work of a broad research community and would not have been possible without the many contributions and support of our colleagues . we especially want to thank a few individuals who have made particularly valuable contributions . william d. boggs , bernd brgmann , alessandra buonanno , mark hannam , richard matzner , cole miller , and harald pfeiffer gave insightful and helpful comments on our article . manuela campanelli and harald pfeiffer supplied us with figures from their simulations that were not in the published literature . we also benefited from many useful discussions with sean mcwilliams . we acknowledge support from nasa grant no . 06-befs06 - 19 . b.j.k . was supported in part by an appointment to the nasa postdoctoral program at the goddard space flight center , administered by oak ridge associated universities through a contract with nasa .
Abstract:

Answer: understanding the predictions of general relativity for the dynamical interactions of two black holes has been a long - standing unsolved problem in theoretical physics . black - hole mergers are monumental astrophysical events , releasing tremendous amounts of energy in the form of gravitational radiation , and are key sources for both ground- and space - based gravitational - wave detectors . the black - hole merger dynamics and the resulting gravitational waveforms can only be calculated through numerical simulations of einstein s equations of general relativity . for many years , numerical relativists attempting to model these mergers encountered a host of problems , causing their codes to crash after just a fraction of a binary orbit could be simulated . recently , however , a series of dramatic advances in numerical relativity has allowed stable , robust black - hole merger simulations . this remarkable progress in the rapidly maturing field of numerical relativity , and the new understanding of black - hole binary dynamics that is emerging is chronicled . important applications of these fundamental physics results to astrophysics , to gravitational - wave astronomy , and in other areas are also discussed .