Collisions of boosted black holes: perturbation theory prediction of gravitational radiation

We consider general relativistic Cauchy data representing two nonspinning, equal-mass black holes boosted toward each other. When the black holes are close enough to each other and their momentum is sufficiently high, an encompassing apparent horizon is present so the system can be viewed as a single, perturbed black hole. We employ gauge-invariant perturbation theory, and integrate the Zerilli equation to analyze these time-asymmetric data sets and compute gravitational wave forms and emitted energies. When coupled with a simple Newtonian analysis of the infall trajectory, we find striking agreement between the perturbation calculation of emitted energies and the results of fully general relativistic numerical simulations of time-symmetric initial data.Comment: 5 pages (RevTex 3.0 with 3 uuencoded figures), CRSR-107

Cauchy-perturbative matching and outer boundary conditions I: Methods and tests

We present a new method of extracting gravitational radiation from three-dimensional numerical relativity codes and providing outer boundary conditions. Our approach matches the solution of a Cauchy evolution of Einstein's equations to a set of one-dimensional linear wave equations on a curved background. We illustrate the mathematical properties of our approach and discuss a numerical module we have constructed for this purpose. This module implements the perturbative matching approach in connection with a generic three-dimensional numerical relativity simulation. Tests of its accuracy and second-order convergence are presented with analytic linear wave data.Comment: 13 pages, 6 figures, RevTe

Cauchy-perturbative matching and outer boundary conditions: computational studies

We present results from a new technique which allows extraction of gravitational radiation information from a generic three-dimensional numerical relativity code and provides stable outer boundary conditions. In our approach we match the solution of a Cauchy evolution of the nonlinear Einstein field equations to a set of one-dimensional linear equations obtained through perturbation techniques over a curved background. We discuss the validity of this approach in the case of linear and mildly nonlinear gravitational waves and show how a numerical module developed for this purpose is able to provide an accurate and numerically convergent description of the gravitational wave propagation and a stable numerical evolution.Comment: 20 pages, RevTe

A multiple scales approach to crack front waves

Perturbation of a propagating crack with a straight edge is solved using the method of matched asymptotic expansions (MAE). This provides a simplified analysis in which the inner and outer solutions are governed by distinct mechanics. The inner solution contains the explicit perturbation and is governed by a quasi-static equation. The outer solution determines the radiation of energy away from the tip, and requires solving dynamic equations in the unperturbed configuration. The outer and inner expansions are matched via the small parameter L/l defined by the disparate length scales: the crack perturbation length L and the outer length scale l associated with the loading. The method is first illustrated for a scalar crack model and then applied to the elastodynamic mode I problem. The dispersion relation for crack front waves is found by requiring that the energy release rate is unaltered under perturbation. The wave speed is calculated as a function of the nondimensional parameter kl where k is the crack front wavenumber, and dispersive properties of the crack front wave speed are described for the first time. The example problems considered here demonstrate that the potential of using MAE for moving boundary value problems with multiple scales.Comment: 25 pages, 5 figure

Understanding initial data for black hole collisions

Numerical relativity, applied to collisions of black holes, starts with initial data for black holes already in each other's strong field. The initial hypersurface data typically used for computation is based on mathematical simplifying prescriptions, such as conformal flatness of the 3-geometry and longitudinality of the extrinsic curvature. In the case of head on collisions of equal mass holes, there is evidence that such prescriptions work reasonably well, but it is not clear why, or whether this success is more generally valid. Here we study these questions by considering the ``particle limit'' for head on collisions of nonspinning holes. Einstein's equations are linearized in the mass of the small hole, and described by a single gauge invariant spacetime function psi, for each multipole. The resulting equations have been solved by numerical evolution for collisions starting from various initial separations, and the evolution is studied on a sequence of hypersurfaces. In particular, we extract hypersurface data, that is psi and its time derivative, on surfaces of constant background Schwarzschild time. These evolved data can then be compared with ``prescribed'' data, evolved data can be replaced by prescribed data on any hypersurface, and evolved further forward in time, a gauge invariant measure of deviation from conformal flatness can be evaluated, etc. The main findings of this study are: (i) For holes of unequal mass the use of prescribed data on late hypersurfaces is not successful. (ii) The failure is likely due to the inability of the prescribed data to represent the near field of the smaller hole. (iii) The discrepancy in the extrinsic curvature is more important than in the 3-geometry. (iv) The use of the more general conformally flat longitudinal data does not notably improve this picture.Comment: 20 pages, REVTEX, 26 PS figures include

Towards a Singularity-Proof Scheme in Numerical Relativity

Progress in numerical relativity has been hindered for 30 years because of the difficulties of avoiding spacetime singularities in numerical evolution. We propose a scheme which excises a region inside an apparent horizon containing the singularity. Two major ingredients of the scheme are the use of a horizon-locking coordinate and a finite differencing which respects the causal structure of the spacetime. Encouraging results of the scheme in the spherical collapse case are given.Comment: 9 page

Robust evolution system for Numerical Relativity

The paper combines theoretical and applied ideas which have been previously considered separately into a single set of evolution equations for Numerical Relativity. New numerical ingredients are presented which avoid gauge pathologies and allow one to perform robust 3D calculations. The potential of the resulting numerical code is demonstrated by using the Schwarzschild black hole as a test-bed. Its evolution can be followed up to times greater than one hundred black hole masses.Comment: 11 pages, 4 figures; figure correcte

Einstein and Yang-Mills theories in hyperbolic form without gauge-fixing

The evolution of physical and gauge degrees of freedom in the Einstein and Yang-Mills theories are separated in a gauge-invariant manner. We show that the equations of motion of these theories can always be written in flux-conservative first-order symmetric hyperbolic form. This dynamical form is ideal for global analysis, analytic approximation methods such as gauge-invariant perturbation theory, and numerical solution.Comment: 12 pages, revtex3.0, no figure

Initial data for Einstein's equations with superposed gravitational waves

A method is presented to construct initial data for Einstein's equations as a superposition of a gravitational wave perturbation on an arbitrary stationary background spacetime. The method combines the conformal thin sandwich formalism with linear gravitational waves, and allows detailed control over characteristics of the superposed gravitational wave like shape, location and propagation direction. It is furthermore fully covariant with respect to spatial coordinate changes and allows for very large amplitude of the gravitational wave.Comment: Version accepted by PRD; added convergence plots, expanded discussion. 9 pages, 9 figure