Search

The most popular method to construct initial data for black-hole-binary simulations is the puncture method, in which compactified wormholes are given linear and angular momentum via the Bowen-York extrinsic curvature. When these data are evolved, they quickly approach a ``trumpet'' topology, suggesting that it would be preferable to use data that are in trumpet form from the outset. To achieve this, we extend the puncture method to allow the construction of Bowen-York trumpets, including an outline of an existence and uniqueness proof of the solutions. We construct boosted, spinning and binary Bowen-York puncture trumpets using a single-domain pseudospectral elliptic solver, and evolve the binary data and compare with standard wormhole-data results. We also show that for boosted trumpets the black-hole mass can be prescribed {\it a priori}, without recourse to the iterative procedure that is necessary for wormhole data., Comment: 15 pages, 14 figures. Published version

Binary black-hole systems with spins aligned or anti-aligned to the orbital angular momentum provide the natural ground to start detailed studies of the influence of strong-field spin effects on gravitational wave observations of coalescing binaries. Furthermore, such systems may be the preferred end-state of the inspiral of generic supermassive binary black-hole systems. In view of this, we have computed the inspiral and merger of a large set of binary systems of equal-mass black holes with spins parallel to the orbital angular momentum but otherwise arbitrary. Our attention is particularly focused on the gravitational-wave emission so as to quantify how much spin effects contribute to the signal-to-noise ratio, to the horizon distances, and to the relative event rates for the representative ranges in masses and detectors. As expected, the signal-to-noise ratio increases with the projection of the total black hole spin in the direction of the orbital momentum. We find that equal-spin binaries with maximum spin aligned with the orbital angular momentum are more than "three times as loud" as the corresponding binaries with anti-aligned spins, thus corresponding to event rates up to 30 times larger. We also consider the waveform mismatch between the different spinning configurations and find that, within our numerical accuracy, binaries with opposite spins S_1=-S_2 cannot be distinguished whereas binaries with spin S_1=S_2 have clearly distinct gravitational-wave emissions. Finally, we derive a simple expression for the energy radiated in gravitational waves and find that the binaries always have efficiencies E_rad/M > 3.6%, which can become as large as E_rad/M = 10% for maximally spinning binaries with spins aligned with the orbital angular momentum., Comment: 18 pages, 11 figures, matches published version

The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise., Comment: Proceedings of Numerical Relativity and Data Analysis NRDA 2008 (Syracuse NY, August 2008) - 14 pages, 1 figure

We analyze stationary slicings of the Schwarzschild spacetime defined by members of the Bona-Masso family of slicing conditions. Our main focus is on the influence of a non-vanishing offset to the extrinsic curvature, which forbids the existence of standard Cauchy foliations but at the same time allows gauge choices that are adapted to include null infinity (scri) in the evolution. These hyperboloidal slicings are especially interesting for observing outgoing gravitational waves. We show that the standard 1+log slicing condition admits no overall regular hyperboloidal slicing, but by appropriately combining with harmonic slicing, we construct a gauge condition that leads to a strongly singularity-avoiding hyperboloidal foliation that connects the black hole to scri., Comment: 18 pages, 11 figures, LaTeX; minor changes for the published version

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses., Comment: 56 pages, 25 figures; various clarifications; accepted to CQG

We quantify the consistency of numerical-relativity black-hole-binary waveforms for use in gravitational-wave (GW) searches with current and planned ground-based detectors. We compare previously published results for the $(\ell=2,| m | =2)$ mode of the gravitational waves from an equal-mass nonspinning binary, calculated by five numerical codes. We focus on the 1000M (about six orbits, or 12 GW cycles) before the peak of the GW amplitude and the subsequent ringdown. We find that the phase and amplitude agree within each code's uncertainty estimates. The mismatch between the $(\ell=2,| m| =2)$ modes is better than $10^{-3}$ for binary masses above $60 M_{\odot}$ with respect to the Enhanced LIGO detector noise curve, and for masses above $180 M_{\odot}$ with respect to Advanced LIGO, Virgo and Advanced Virgo. Between the waveforms with the best agreement, the mismatch is below $2 \times 10^{-4}$. We find that the waveforms would be indistinguishable in all ground-based detectors (and for the masses we consider) if detected with a signal-to-noise ratio of less than $\approx14$, or less than $\approx25$ in the best cases., Comment: 17 pages, 9 figures. Version accepted by PRD

These notes summarize basic concepts underlying numerical relativity and in particular the numerical modeling of black hole dynamics as a source of gravitational waves. Main topics are the 3+1 decomposition of general relativity, the concept of a well-posed initial value problem, the construction of initial data for general relativity, trapped surfaces and gravitational waves. Also, a brief summary is given of recent progress regarding the numerical evolution of black hole binary systems., Comment: 28 pages, lectures given at winter school 'Conceptual and Numerical Challenges in Femto- and Peta-Scale Physics' in Schladming, Austria, 2007

We study gravitational perturbations of Schwarzschild spacetime by solving a hyperboloidal initial value problem for the Bardeen-Press equation. Compactification along hyperboloidal surfaces in a scri-fixing gauge allows us to have access to the gravitational waveform at null infinity in a general setup. We argue that this hyperboloidal approach leads to a more accurate and efficient calculation of the radiation signal than the common approach where a timelike outer boundary is introduced. The method can be generalized to study perturbations of Kerr spacetime using the Teukolsky equation., Comment: 14 pages, 9 figures

We study the angular resolution of the gravitational wave detector LISA and show that numerical relativity can drastically improve the accuracy of position location for coalescing Super Massive Black Hole (SMBH) binaries. For systems with total redshifted mass above $10^7 M_{\odot}$, LISA will mainly see the merger and ring-down of the gravitational wave (GW) signal, which can now be computed numerically using the full Einstein equations. Using numerical waveforms that also include about ten GW cycles of inspiral, we improve inspiral-only position estimates by an order of magnitude. We show that LISA localizes half of all such systems at $z=1$ to better than 3 arcminutes and the best 20% to within one arcminute. This will give excellent prospects for identifying the host galaxy., Comment: 4 pages, 1 figure

We expand upon our previous analysis of numerical moving-puncture simulations of the Schwarzschild spacetime. We present a derivation of the family of analytic stationary 1+log foliations of the Schwarzschild solution, and outline a transformation to isotropic-like coordinates. We discuss in detail the numerical evolution of standard Schwarzschild puncture data, and the new time-independent 1+log data. Finally, we demonstrate that the moving-puncture method can locate the appropriate stationary geometry in a robust manner when a numerical code alternates between two forms of 1+log slicing during a simulation., Comment: 19 pages, 22 figures

The Effective-One-Body (EOB) formalism contains several flexibility parameters, notably $a_5$, $\vp$ and $\a$. We show here how to jointly constrain the values of these parameters by simultaneously best-fitting the EOB waveform to two, independent, numerical relativity (NR) simulations of inspiralling and/or coalescing binary black hole systems: published Caltech-Cornell {\it inspiral} data (considered for gravitational wave frequencies $M\omega\leq 0.1$) on one side, and newly computed {\it coalescence} data on the other side. The resulting, approximately unique, "best-fit" EOB waveform is then shown to exhibit excellent agreement with NR coalescence data for several mass ratios. The dephasing between this best-fit EOB waveform and published Caltech-Cornell inspiral data is found to vary between -0.0014 and +0.0008 radians over a time span of $\sim 2464M$ up to gravitational wave frequency $M\omega= 0.1$, and between +0.0013 and -0.0185 over a time span of 96M after $M\omega=0.1$ up to $M\omega=0.1565$. The dephasings between EOB and the new coalescence data are found to be smaller than: (i) $\pm 0.025$ radians over a time span of 730M (11 cycles) up to merger, in the equal mass case, and (ii) $\pm 0.05$ radians over a time span of about 950M (17 cycles) up to merger in the 2:1 mass-ratio case. These new results corroborate the aptitude of the EOB formalism to provide accurate representations of general relativistic waveforms, which are needed by currently operating gravitational wave detectors., Comment: 25 pages, 15 figures. Improved discussion about errors on numerical relativity data. Version published in Phys. Rev. D

We compare results from numerical simulations of spinning binaries in the "orbital hangup" case, where the binary completes at least nine orbits before merger, with post-Newtonian results using the approximants TaylorT1, T4 and Et. We find that, over the ten cycles before the gravitational-wave frequency reaches $M\omega = 0.1$, the accumulated phase disagreement between NR and 2.5PN results is less than three radians, and is less than 2.5 radians when using 3.5PN results. The amplitude disagreement between NR and restricted PN results increases with the black holes' spin, from about 6% in the equal-mass case to 12% when the black holes' spins are $S_i/M_i^2 = 0.85$. Finally, our results suggest that the merger waveform will play an important role in estimating the spin from such inspiral waveforms., Comment: 10 pages, 6 figures. Matches version published in PRD

We compare the phase evolution of equal-mass nonspinning black-hole binaries from numerical relativity (NR) simulations with post-Newtonian (PN) results obtained from three PN approximants: the TaylorT1 and T4 approximants, for which NR-PN comparisons have already been performed in the literature, and the recently proposed approximant TaylorEt. The accumulated phase disagreement between NR and PN results over the frequency range $M\omega = 0.0455$ to $M\omega = 0.1$ is greater for TaylorEt than either T1 or T4, but has the attractive property of decreasing monotonically as the PN order is increased., Comment: 6 pages, 4 figures

We provide a compact analytic formula to compute the spin of the black hole produced by the coalescence of two black holes. The expression, which uses an analytic fit of numerical-relativity data and relies on four assumptions, aims at modelling generic initial spin configurations and mass ratios. A comparison with numerical-relativity simulations already shows very accurate agreements with all of the numerical data available to date, but we also suggest a number of ways in which our predictions can be further improved., Comment: Extended discussion of physical assumptions to make them less cryptic; matches published version

Recent calculations of the recoil velocity in black-hole binary mergers have found kick velocities of $\approx2500 $km/s for equal-mass binaries with anti-aligned initial spins in the orbital plane. In general the dynamics of spinning black holes can be extremely complicated and are difficult to analyze and understand. In contrast, the ``superkick'' configuration is an example with a high degree of symmetry that also exhibits exciting physics. We exploit the simplicity of this ``test case'' to study more closely the role of spin in black-hole recoil and find that: the recoil is with good accuracy proportional to the difference between the $(l = 2, m = \pm 2)$ modes of $\Psi_4$, the major contribution to the recoil occurs within $30M$ before and after the merger, and that this is after the time at which a standard post-Newtonian treatment breaks down. We also discuss consequences of the $(l = 2, m = \pm 2)$ asymmetry in the gravitational wave signal for the angular dependence of the SNR and the mismatch of the gravitational wave signals corresponding to the north and south poles.

We analyze numerical-relativity (NR) waveforms that cover nine orbits (18 gravitational-wave cycles) before merger of an equal-mass system with low eccentricity, with numerical uncertainties of 0.25 radians in the phase and less than 2% in the amplitude; such accuracy allows a direct comparison with post-Newtonian (PN) waveforms. We focus on one of the PN approximants that has been proposed for use in gravitational-wave data analysis, the restricted 3.5PN ``TaylorT1'' waveforms, and compare these with a section of the numerical waveform from the second to the eighth orbit, which is about one and a half orbits before merger. This corresponds to a gravitational-wave frequency range of $M\omega = 0.0455$ to 0.1. Depending on the method of matching PN and NR waveforms, the accumulated phase disagreement over this frequency range can be within numerical uncertainty. Similar results are found in comparisons with an alternative PN approximant, 3PN ``TaylorT3''. The amplitude disagreement, on the other hand, is around 6%, but roughly constant for all 13 cycles that are compared, suggesting that only 4.5 orbits need be simulated to match PN and NR waves with the same accuracy as is possible with nine orbits. If, however, we model the amplitude up to 2.5PN order, the amplitude disagreement is roughly within numerical uncertainty up to about 11 cycles before merger., Comment: 14 pages, 18 figures. Modifications resulting from bug fixes in LAL, and extended analysis of numerical errors and phase agreement with PN, now including the 3PN TaylorT3 approximant. No change to main conclusions

Standard choices of quasi-circular orbit parameters for black-hole binary evolutions result in eccentric inspiral. We introduce a conceptually simple method, which is to integrate the post-Newtonian equations of motion through hundreds of orbits, and read off the values of the momenta at the separation at which we wish to start a fully general relativistic numerical evolution. For the particular case of non-spinning equal-mass inspiral with an initial coordinate separation of $D = 11M$ we show that this approach reduces the eccentricity by at least a factor of five to $e < 0.002$ as compared to using standard quasi-circular initial parameters., Comment: 6 pages, 3 figures, 1 table

We describe a modification of a fourth-order accurate ``moving puncture'' evolution code, where by replacing spatial fourth-order accurate differencing operators in the bulk of the grid by a specific choice of sixth-order accurate stencils we gain significant improvements in accuracy. We illustrate the performance of the modified algorithm with an equal-mass simulation covering nine orbits., Comment: 13 pages, 6 figures

We discuss possible origins for discrepancies observed in the radiated energies in head-on collisions of non-spinning binaries starting from Brill-Lindquist and superposed Kerr-Schild data. For this purpose, we discuss the impact of different choices of gauge parameters and a small initial boost of the black holes., Comment: Proceedings of the Eleventh Marcel Grossmann Meeting; 3 pages (limit)

We study the inspiral, merger and ringdown of unequal mass black hole binaries by analyzing a catalogue of numerical simulations for seven different values of the mass ratio (from q=M2/M1=1 to q=4). We compare numerical and Post-Newtonian results by projecting the waveforms onto spin-weighted spherical harmonics, characterized by angular indices (l,m). We find that the Post-Newtonian equations predict remarkably well the relation between the wave amplitude and the orbital frequency for each (l,m), and that the convergence of the Post-Newtonian series to the numerical results is non-monotonic. To leading order the total energy emitted in the merger phase scales like eta^2 and the spin of the final black hole scales like eta, where eta=q/(1+q)^2 is the symmetric mass ratio. We study the multipolar distribution of the radiation, finding that odd-l multipoles are suppressed in the equal mass limit. Higher multipoles carry a larger fraction of the total energy as q increases. We introduce and compare three different definitions for the ringdown starting time. Applying linear estimation methods (the so-called Prony methods) to the ringdown phase, we find resolution-dependent time variations in the fitted parameters of the final black hole. By cross-correlating information from different multipoles we show that ringdown fits can be used to obtain precise estimates of the mass and spin of the final black hole, which are in remarkable agreement with energy and angular momentum balance calculations., Comment: 51 pages, 28 figures, 16 tables. Many improvements throughout the text in response to the referee report. The calculation of multipolar components in Appendix A now uses slightly different conventions. Matches version in press in PRD

Binary black hole simulations have traditionally been computationally very expensive: current simulations are performed in supercomputers involving dozens if not hundreds of processors, thus systematic studies of the parameter space of binary black hole encounters still seem prohibitive with current technology. Here we show how the multi-layered refinement level code BAM can be used on dual processor workstations to simulate certain binary black hole systems. BAM, based on the moving punctures method, provides grid structures composed of boxes of increasing resolution near the center of the grid. In the case of binaries, the highest resolution boxes are placed around each black hole and they track them in their orbits until the final merger when a single set of levels surrounds the black hole remnant. This is particularly useful when simulating spinning black holes since the gravitational fields gradients are larger. We present simulations of binaries with equal mass black holes with spins parallel to the binary axis and intrinsic magnitude of S/m^2= 0.75. Our results compare favorably to those of previous simulations of this particular system. We show that the moving punctures method produces stable simulations at maximum spatial resolutions up to M/160 and for durations of up to the equivalent of 20 orbital periods., Comment: 20 pages, 8 figures. Final version, to appear in a special issue of Class. Quantum Grav. based on the New Frontiers in Numerical Relativity Conference, Golm, July 2006

We study in spherical symmetry the conformal compactification for hyperboloidal foliations with nonvanishing constant mean curvature. The conformal factor and the coordinates are chosen such that null infinity is at a fixed radial coordinate location., Comment: 3 pages, 1 figure; uses World Scientific style; submitted to the Proceedings of the Eleventh Marcel Grossmann Meeting, July 23-29, 2006, Berlin

Currently the most popular method to evolve black-hole binaries is the ``moving puncture'' method. It has recently been shown that when puncture initial data for a Schwarzschild black hole are evolved using this method, the numerical slices quickly lose contact with the second asymptotically flat end, and end instead on a cylinder of finite Schwarzschild coordinate radius. These slices are stationary, meaning that their geometry does not evolve further. We will describe these results in the context of maximal slices, and present time-independent puncture-like data for the Schwarzschild spacetime., Comment: Proceedings for 29th Spanish Relativity Meeting. Added more details about the time-independent solution, with reference to the analytic result of Baumgarte and Naculich

It is well-known that Bowen-York initial data contain spurious radiation. Although this ``junk'' radiation has been seen to be small for non-spinning black-hole binaries in circular orbit, its magnitude increases when the black holes are given spin. It is possible to reduce the spurious radiation by applying the puncture approach to multiple Kerr black holes, as we demonstrate for examples of head-on collisions of equal-mass black-hole binaries., Comment: 10 pages, 2 figures, submitted to special "New Frontiers in Numerical Relativity" issue of Classical and Quantum Gravity

When unequal-mass black holes merge, the final black hole receives a ``kick'' due to the asymmetric loss of linear momentum in the gravitational radiation emitted during the merger. The magnitude of this kick has important astrophysical consequences. Recent breakthroughs in numerical relativity allow us to perform the largest parameter study undertaken to date in numerical simulations of binary black hole inspirals. We study non-spinning black-hole binaries with mass ratios from $q=M_1/M_2=1$ to $q =0.25$ ($\eta = q/(1 + q)^2$ from 0.25 to 0.16). We accurately calculate the velocity of the kick to within 6%, and the final spin of the black holes to within 2%. A maximum kick of $175.2\pm11$ km s$^{-1}$ is achieved for $\eta = 0.195 \pm 0.005$., Comment: 4 pages, 4 figures. Version accepted by PRL

We present single and binary black hole simulations that follow the moving puncture paradigm of simulating black-hole spacetimes without excision, and use moving boxes mesh refinement. Focussing on binary black hole configurations where the simulations cover roughly two orbits, we address five major issues determining the quality of our results: numerical discretization error, finite extraction radius of the radiation signal, physical appropriateness of initial data, gauge choice and computational performance. We also compare results we have obtained with the BAM code described here with the independent LEAN code., Comment: 21 pages, 21 figures, three tables

Significant advances in numerical simulations of black-hole binaries have recently been achieved using the puncture method. We examine how and why this method works by evolving a single black hole. The coordinate singularity and hence the geometry at the puncture are found to change during evolution, from representing an asymptotically flat end to being a cylinder. We construct an analytic solution for the stationary state of a black hole in spherical symmetry that matches the numerical result and demonstrates that the evolution is not dominated by artefacts at the puncture but indeed finds the analytical result., Comment: 4 pages, 2 figures. Replaced with version that matches the one published in PRL: one extra figure, and modified abstract and introduction

We discuss the hyperboloidal evolution problem in general relativity from a numerical perspective, and present some new results. Families of initial data which are the hyperboloidal analogue of Brill waves are constructed numerically, and a systematic search for apparent horizons is performed. Schwarzschild-Kruskal spacetime is discussed as a first application of Friedrich's general conformal field equations in spherical symmetry, and the Maxwell equations are discussed on a nontrivial background as a toy model for continuum instabilities., Comment: 11 pages, 9 figures. To appear in the Proceedings of the Spanish Relativity Meeting (ERE 2005), Oviedo, Spain, 6-10 Sept 2005

We extend the notion of numerical stability of finite difference approximations to include hyperbolic systems that are first order in time and second order in space, such as those that appear in Numerical Relativity. By analyzing the symbol of the second order system, we obtain necessary and sufficient conditions for stability in a discrete norm containing one-sided difference operators. We prove stability for certain toy models and the linearized Nagy-Ortiz-Reula formulation of Einstein's equations. We also find that, unlike in the fully first order case, standard discretizations of some well-posed problems lead to unstable schemes and that the Courant limits are not always simply related to the characteristic speeds of the continuum problem. Finally, we propose methods for testing stability for second order in space hyperbolic systems., Comment: 18 pages, 9 figures

In this paper we present our recent work in developing a computer-algebra tool for systems of partial differential equations (PDEs), termed "Kranc". Our work is motivated by the problem of finding solutions of the Einstein equations through numerical simulations. Kranc consists of Mathematica based computer-algebra packages, that facilitate the task of dealing with symbolic tensorial calculations and realize the conversion of systems of partial differential evolution equations into parallelized C or Fortran code., Comment: LaTeX llncs style, 9 pages, 1 figure, to appaer in the proceedings of "SYNASC 2004 - 6th International Symposium on Symbolic and Numeric Algorithms for Scientific Computing", Timisoara, Romania, September 26-30 2004

We discuss critical gravitational collapse on the threshold of apparent horizon formation as a model both for the discussion of global aspects of critical collapse and for numerical studies in a compactified context. For our matter model we choose a self-gravitating massless scalar field in spherical symmetry, which has been studied extensively in the critical collapse literature. Our evolution system is based on Bondi coordinates, the mass function is used as an evolution variable to ensure regularity at null infinity. We compute radiation quantities like the Bondi mass and news function and find that they reflect the DSS behavior. Surprisingly, the period of radiation at null infinity is related to the formal result for the leading quasi-normal mode of a black hole with rapidly decreasing mass. Furthermore, our investigations shed some light on global versus local issues in critical collapse, and the validity and usefulness of the concept of null infinity when predicting detector signals., Comment: 17 pages, 13 figures, uses RevTeX4, changed order of figures and corrected references to match PRD version

We present a suite of Mathematica-based computer-algebra packages, termed "Kranc", which comprise a toolbox to convert (tensorial) systems of partial differential evolution equations to parallelized C or Fortran code. Kranc can be used as a "rapid prototyping" system for physicists or mathematicians handling very complicated systems of partial differential equations, but through integration into the Cactus computational toolkit we can also produce efficient parallelized production codes. Our work is motivated by the field of numerical relativity, where Kranc is used as a research tool by the authors. In this paper we describe the design and implementation of both the Mathematica packages and the resulting code, we discuss some example applications, and provide results on the performance of an example numerical code for the Einstein equations., Comment: 24 pages, 1 figure. Corresponds to journal version

We study the properties of the outgoing gravitational wave produced when a non-spinning black hole is excited by an ingoing gravitational wave. Simulations using a numerical code for solving Einstein's equations allow the study to be extended from the linearized approximation, where the system is treated as a perturbed Schwarzschild black hole, to the fully nonlinear regime. Several nonlinear features are found which bear importance to the data analysis of gravitational waves. When compared to the results obtained in the linearized approximation, we observe large phase shifts, a stronger than linear generation of gravitational wave output and considerable generation of radiation in polarization states which are not found in the linearized approximation. In terms of a spherical harmonic decomposition, the nonlinear properties of the harmonic amplitudes have simple scaling properties which offer an economical way to catalog the details of the waves produced in such black hole processes., Comment: 17 pages, 20 figures, abstract and introduction re-written

In recent years, many different numerical evolution schemes for Einstein's equations have been proposed to address stability and accuracy problems that have plagued the numerical relativity community for decades. Some of these approaches have been tested on different spacetimes, and conclusions have been drawn based on these tests. However, differences in results originate from many sources, including not only formulations of the equations, but also gauges, boundary conditions, numerical methods, and so on. We propose to build up a suite of standardized testbeds for comparing approaches to the numerical evolution of Einstein's equations that are designed to both probe their strengths and weaknesses and to separate out different effects, and their causes, seen in the results. We discuss general design principles of suitable testbeds, and we present an initial round of simple tests with periodic boundary conditions. This is a pivotal first step toward building a suite of testbeds to serve the numerical relativists and researchers from related fields who wish to assess the capabilities of numerical relativity codes. We present some examples of how these tests can be quite effective in revealing various limitations of different approaches, and illustrating their differences. The tests are presently limited to vacuum spacetimes, can be run on modest computational resources, and can be used with many different approaches used in the relativity community., Comment: REVTeX4, 19 pages

We discuss the application of computer algebra to problems commonly arising in numerical relativity, such as the derivation of 3+1-splits, manipulation of evolution equations and automatic code generation. Particular emphasis is put on working with abstract index tensor quantities as much as possible., Comment: 5 pages, no figures. To appear in the Proceedings of the Spanish Relativity Meeting (ERE 2002), Mao, Menorca, Spain, 22-24 Sept 2002

We present a method for computing the evolution of a spacetime containing a massive particle and a black hole. The essential idea is that the gravitational field is evolved using full numerical relativity, with the particle generating a non-zero source term in the Einstein equations. The matter fields are not evolved by hydrodynamic equations. Instead the particle is treated as a rigid body whose center follows a geodesic. The necessary theoretical framework is developed and then implemented in a computer code that uses the null-cone, or characteristic, formulation of numerical relativity. The performance of the code is illustrated in test runs, including a complete orbit (near r=9M) of a Schwarzschild black hole., Comment: 12 pages, 9 figures, RevTeX4, to appear in Phys. Rev. D

We present a fully nonlinear calculation of the waveform of the gravitational radiation emitted in the fission of a vacuum white hole. At early times, the waveforms agree with close-approximation perturbative calculations but they reveal dramatic time and angular dependence in the nonlinear regime. The results pave the way for a subsequent computation of the radiation emitted after a binary black hole merger., Comment: 11 pages, 6 figures, RevTeX4

I discuss the conformal approach to the numerical simulation of radiating isolated systems in general relativity. The method is based on conformal compactification and a reformulation of the Einstein equations in terms of rescaled variables, the so-called ``conformal field equations'' developed by Friedrich. These equations allow to include ``infinity'' on a finite grid, solving regular equations, whose solutions give rise to solutions of the Einstein equations of (vacuum) general relativity. The conformal approach promises certain advantages, in particular with respect to the treatment of radiation extraction and boundary conditions. I will discuss the essential features of the analytical approach to the problem, previous work on the problem - in particular a code for simulations in 3+1 dimensions, some new results, open problems and strategies for future work., Comment: 34 pages, submitted to the Proceedings of the 2001 Spanish Relativity meeting, eds. L. Fernandez and L. Gonzalez, to be published by Springer, Lecture Notes in Physics series

This talk reports on the status of an approach to the numerical study of isolated systems with the conformal field equations. We first describe the algorithms used in a code which has been developed at AEI in the last years, and discuss a milestone result obtained by Huebner. Then we present more recent results as examples to sketch the problems we face in the conformal approach to numerical relativity and outline a possible roadmap toward making this approach a practical tool., Comment: 20 pages, to appear in the proceedings of the conference ``The Conformal Structure of Spacetimes: Geometry, Analysis, Numerics'', ed. by J. Frauendiener and H. Friedrich, to be published by Springer Verlag, Lecture Notes in Physics series

We analyze a bifurcation phenomenon associated with critical gravitational collapse in a family of self-gravitating SU(2) $\sigma$-models. As the dimensionless coupling constant decreases, the critical solution changes from discretely self-similar (DSS) to continuously self-similar (CSS). Numerical results provide evidence for a bifurcation which is analogous to a heteroclinic loop bifurcation in dynamical systems, where two fixed points (CSS) collide with a limit cycle (DSS) in phase space as the coupling constant tends to a critical value., Comment: 4 pages, 5 figures, uses REVTeX

This talk reports on recent progress toward the semiglobal study of asymptotically flat spacetimes within numerical relativity. The development of a 3D solver for asymptotically Minkowski-like hyperboloidal initial data has rendered possible the application of Friedrich's conformal field equations to astrophysically interesting spacetimes. As a first application, the whole future of a hyperboloidal set of weak initial data has been studied, including future null and timelike infinity. Using this example we sketch the numerical techniques employed and highlight some of the unique capabilities of the numerical code. We conclude with implications for future work., Comment: 6 pages, published in "Relativistic Astrophysics: 20th Texas Symposium", ed. by J. Craig Wheeler and Hugo Martel, AIP Conference Proceedings 586 (Austin, Texas, 10-15 Dec. 2000)

We use null hypersurface techniques in a new approach to calculate the retarded waveform from a binary black hole merger in the close approximation. The process of removing ingoing radiation from the system leads to two notable features in the shape of the close approximation waveform for a head-on collision of black holes: (i) an initial quasinormal ringup and (ii) weak sensitivity to the parameter controlling the collision velocity. Feature (ii) is unexpected and has the potential importance of enabling the design of an efficient template for extracting the gravitational wave signal from the noise., Comment: 14 pages, 13 figures, revtex4

We discuss a sequence of numerically constructed geometries describing binary black hole event horizons -- providing the necessary input for characteristic evolution of the exterior spacetime. Our sequence approaches a single Schwarzschild horizon as one limiting case and also includes cases where the horizon's crossover surface is not hidden by a marginally anti-trapped surface (MATS)., Comment: 2 pages, uses World Scientific style file ws-p9-75x6-50.cls (enclosed), contribution to the 9th Marcel Grossmann meeting (MG9), Rome, July 2000

We present a characteristic algorithm for computing the perturbations of a Schwarzschild spacetime by means of solving the Teukolsky equations. Our methods and results are expected to have direct bearing on the study of binary black holes presently underway using a fully {\em nonlinear} characteristic code \cite{Gomez98a}., Comment: 2 pages, contribution to the 9th Marcel Grossmann meeting (MG9), Rome, July 2000

We present a characteristic algorithm for computing the perturbation of a Schwarzschild spacetime by means of solving the Teukolsky equation. We implement the algorithm as a characteristic evolution code and apply it to compute the advanced solution to a black hole collision in the close approximation. The code successfully tracks the initial burst and quasinormal decay of a black hole perturbation through 10 orders of magnitude and tracks the final power law decay through an additional 6 orders of magnitude. Determination of the advanced solution, in which ingoing radiation is absorbed by the black hole but no outgoing radiation is emitted, is the first stage of a two stage approach to determining the retarded solution, which provides the close approximation waveform with the physically appropriate boundary condition of no ingoing radiation., Comment: Revised version, published in Phys. Rev. D, 34 pages, 13 figures, RevTeX

We present an algorithm for calculating the complete data on an event horizon which constitute the necessary input for characteristic evolution of the exterior spacetime. We apply this algorithm to study the intrinsic and extrinsic geometry of a binary black hole event horizon, constructing a sequence of binary black hole event horizons which approaches a single Schwarzschild black hole horizon as a limiting case. The linear perturbation of the Schwarzschild horizon provides global insight into the close limit for binary black holes, in which the individual holes have joined in the infinite past. In general there is a division of the horizon into interior and exterior regions, analogous to the division of the Schwarzschild horizon by the r=2M bifurcation sphere. In passing from the perturbative to the strongly nonlinear regime there is a transition in which the individual black holes persist in the exterior portion of the horizon. The algorithm is intended to provide the data sets for production of a catalog of nonlinear post-merger wave forms using the PITT null code., Comment: Revised version, to appear in Phys. Rev. D. July 15 (2001), 41 pages, 11 figures, RevTeX/epsf/psfig

We report on the existence and phenomenology of type II critical collapse within the one-parameter family of SU(2) $\sigma$-models coupled to gravity. Numerical investigations in spherical symmetry show discretely self-similar (DSS) behavior at the threshold of black hole formation for values of the dimensionless coupling constant $\ccbeta$ ranging from 0.2 to 100; at 0.18 we see small deviations from DSS. While the echoing period $\Delta$ of the critical solution rises sharply towards the lower limit of this range, the characteristic mass scaling has a critical exponent $\gamma$ which is almost independent of $\ccbeta$, asymptoting to $0.1185 \pm 0.0005$ at large $\ccbeta$. We also find critical scaling of the scalar curvature for near-critical initial data. Our numerical results are based on an outgoing-null-cone formulation of the Einstein-matter equations, specialized to spherical symmetry. Our numerically computed initial-data critical parameters $p^*$ show 2nd order convergence with the grid resolution, and after compensating for this variation in $p^*$, our individual evolutions are uniformly 2nd order convergent even very close to criticality., Comment: 23 pages, includes 10 postscript figure files, uses REVTeX, epsf, psfrag, and AMS math fonts (amstex + amssymb); to appear in PRD15. Summary of revisions from v2: fix wrong formula in figure 6 caption and y-axis label, also minor wording changes and update publication status of refs 5-7

We present a class of numerical solutions to the SU(2) nonlinear $\sigma$-model coupled to the Einstein equations with cosmological constant $\Lambda\geq 0$ in spherical symmetry. These solutions are characterized by the presence of a regular static region which includes a center of symmetry. They are parameterized by a dimensionless ``coupling constant'' $\beta$, the sign of the cosmological constant, and an integer ``excitation number'' $n$. The phenomenology we find is compared to the corresponding solutions found for the Einstein-Yang-Mills (EYM) equations with positive $\Lambda$ (EYM$\Lambda$). If we choose $\Lambda$ positive and fix $n$, we find a family of static spacetimes with a Killing horizon for $0 \leq \beta < \beta_{max}$. As a limiting solution for $\beta = \beta_{max}$ we find a {\em globally} static spacetime with $\Lambda=0$, the lowest excitation being the Einstein static universe. To interpret the physical significance of the Killing horizon in the cosmological context, we apply the concept of a trapping horizon as formulated by Hayward. For small values of $\beta$ an asymptotically de Sitter dynamic region contains the static region within a Killing horizon of cosmological type. For strong coupling the static region contains an ``eternal cosmological black hole''., Comment: 20 pages, 6 figures, Revtex

We study event horizons of non-axisymmetric black holes and show how features found in axisymmetric studies of colliding black holes and of toroidal black holes are non-generic and how new features emerge. Most of the details of black hole formation and black hole merger are known only in the axisymmetric case, in which numerical evolution has successfully produced dynamical space-times. The work that is presented here uses a new approach to construct the geometry of the event horizon, not by locating it in a given spacetime, but by direct construction. In the axisymmetric case, our method produces the familiar pair-of-pants structure found in previous numerical simulations of black hole mergers, as well as event horizons that go through a toroidal epoch as discovered in the collapse of rotating matter. The main purpose of this paper is to show how new - substantially different - features emerge in the non-axisymmetric case. In particular, we show how black holes generically go through a toroidal phase before they become spherical, and how this fits together with the merger of black holes., Comment: 28 pages, 10 figures, uses REVTEX. Improved quality figures and additional color images are provided at http://www.phyast.pitt.edu/~shusa/EH/

We study the utilization of conformal compactification within the conformal approach to solving the constraints of general relativity for asymptotically flat initial data. After a general discussion of the framework, particular attention is paid to simplifications that arise when restricting to a class of initial data which have a certain $U(1)\times U(1)$ conformal symmetry., Comment: Contribution to the 8th Conference on Theoretical Physics, General Relativity and Gravitation in Bistrita-Nasaud, Romania; LaTeX, 13 pages