Evolution of three-dimensional gravitational waves: Harmonic slicing case

Nonlinear metric perturbation enhancement of primordial gravitational waves

We present the evolution of the full set of Einstein equations during preheating after inflation. We study a generic supersymmetric model of hybrid inflation, integrating fields and metric fluctuations in a 3-dimensional lattice. We take initial conditions consistent with Einstein's constraint equations. The induced preheating of the metric fluctuations is not large enough to backreact onto the fields, but preheating of the scalar modes does affect the evolution of vector and tensor modes. In particular, they do enhance the induced stochastic background of gravitational waves during preheating, giving an energy density in general an order of magnitude larger than that obtained by evolving the tensor fluctuations in an homogeneous background metric. This enhancement can improve the expectations for detection by planned gravitational wave observatories.

Characteristic Evolution and Matching

I review the development of numerical evolution codes for general relativity based upon the characteristic initial value problem. Progress is traced from the early stage of 1D feasibility studies to 2D axisymmetric codes that accurately simulate the oscillations and gravitational collapse of relativistic stars and to current 3D codes that provide pieces of a binary black hole spacetime. Cauchy codes have now been successful at simulating all aspects of the binary black hole problem inside an artificially constructed outer boundary. A prime application of characteristic evolution is to eliminate the role of this artificial outer boundary via Cauchy-characteristic matching, by which the radiated waveform can be computed at null infinity. Progress in this direction is discussed.

Spectral Methods for Numerical Relativity

Equations arising in general relativity are usually too complicated to be solved analytically and one must rely on numerical methods to solve sets of coupled partial differential equations. Among the possible choices, this paper focuses on a class called spectral methods in which, typically, the various functions are expanded in sets of orthogonal polynomials or functions. First, a theoretical introduction of spectral expansion is given with a particular emphasis on the fast convergence of the spectral approximation. We then present different approaches to solving partial differential equations, first limiting ourselves to the one-dimensional case, with one or more domains. Generalization to more dimensions is then discussed. In particular, the case of time evolutions is carefully studied and the stability of such evolutions investigated. We then present results obtained by various groups in the field of general relativity by means of spectral methods. Work, which does not involve explicit time-evolutions, is discussed, going from rapidly-rotating strange stars to the computation of black-hole-binary initial data. Finally, the evolution of various systems of astrophysical interest are presented, from supernovae core collapse to black-hole-binary mergers.

Numerical Hydrodynamics and Magnetohydrodynamics in General Relativity

This article presents a comprehensive overview of numerical hydrodynamics and magneto-hydrodynamics (MHD) in general relativity. Some significant additions have been incorporated with respect to the previous two versions of this review (2000, 2003), most notably the coverage of general-relativistic MHD, a field in which remarkable activity and progress has occurred in the last few years. Correspondingly, the discussion of astrophysical simulations in general-relativistic hydrodynamics is enlarged to account for recent relevant advances, while those dealing with general-relativistic MHD are amply covered in this review for the first time. The basic outline of this article is nevertheless similar to its earlier versions, save for the addition of MHD-related issues throughout. Hence, different formulations of both the hydrodynamics and MHD equations are presented, with special mention of conservative and hyperbolic formulations well adapted to advanced numerical methods. A large sample of numerical approaches for solving such hyperbolic systems of equations is discussed, paying particular attention to solution procedures based on schemes exploiting the characteristic structure of the equations through linearized Riemann solvers. As previously stated, a comprehensive summary of astrophysical simulations in strong gravitational fields is also presented. These are detailed in three basic sections, namely gravitational collapse, black-hole accretion, and neutron-star evolutions; despite the boundaries, these sections may (and in fact do) overlap throughout the discussion. The material contained in these sections highlights the numerical challenges of various representative simulations. It also follows, to some extent, the chronological development of the field, concerning advances in the formulation of the gravitational field, hydrodynamics and MHD equations and the numerical methodology designed to solve them. To keep the length of this article reasonable, an effort has been made to focus on multidimensional studies, directing the interested reader to earlier versions of the review for discussions on one-dimensional works.

ELECTRONIC SUPPLEMENTARY MATERIAL

Supplementary material is available for this article at 10.12942/lrr-2008-7.

Geometry and regularity of moving punctures

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.

Recoil velocities from equal-mass binary-black-hole mergers

The final evolution of a binary-black-hole system gives rise to a recoil velocity if an asymmetry is present in the emitted gravitational radiation. Measurements of this effect for nonspinning binaries with unequal masses have pointed out that kick velocities approximately 175 km/s can be reached for a mass ratio approximately 0.36. However, a larger recoil can be obtained for equal-mass binaries if the asymmetry is provided by the spins. Using two independent methods we show that the merger of such binaries yields velocities as large as approximately 440 km/s for black holes having unequal spins that are antialigned and parallel to the orbital angular momentum.

Maximum gravitational recoil

Recent calculations of gravitational radiation recoil generated during black-hole binary mergers have reopened the possibility that a merged binary can be ejected even from the nucleus of a massive host galaxy. Here we report the first systematic study of gravitational recoil of equal-mass binaries with equal, but counteraligned, spins parallel to the orbital plane. Such an orientation of the spins is expected to maximize the recoil. We find that recoil velocity (which is perpendicular to the orbital plane) varies sinusoidally with the angle that the initial spin directions make with the initial linear momenta of each hole and scales up to a maximum of approximately 4000 km s-1 for maximally rotating holes. Our results show that the amplitude of the recoil velocity can depend sensitively on spin orientations of the black holes prior to merger.

Maximum kick from nonspinning black-hole binary inspiral

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 nonspinning black-hole binaries with mass ratios from q=M1/M2=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+/-11 km s(-1) is achieved for eta=0.195+/-0.005.

Formation of black hole and emission of gravitational waves

Numerical simulations were performed for the formation process of rotating black holes. It is suggested that Kerr black holes are formed for wide ranges of initial parameters. The nature of gravitational waves from a test particle falling into a Kerr black hole as well as the development of 3D numerical relativity for the coalescing binary neutron stars are discussed.

Formation of supermassive black holes through fragmentation of torodial supermassive stars

We investigate new paths to supermassive black hole formation by considering the general relativistic evolution of a differentially rotating polytrope with a toroidal shape. We find that this polytrope is unstable to nonaxisymmetric modes, which leads to a fragmentation into self-gravitating, collapsing components. In the case of one such fragment, we apply a simplified adaptive mesh refinement technique to follow the evolution to the formation of an apparent horizon centered on the fragment. This is the first study of the onset of nonaxisymmetric dynamical instabilities of supermassive stars in full general relativity.

Accurate evolution of orbiting binary black holes

We present a detailed analysis of binary black hole evolutions in the last orbit and demonstrate consistent and convergent results for the trajectories of the individual bodies. The gauge choice can significantly affect the overall accuracy of the evolution. It is possible to reconcile certain gauge-dependent discrepancies by examining the convergence limit. We illustrate these results using an initial data set recently evolved by Brügmann et al. [Phys. Rev. Lett. 92, 211101 (2004)10.1103/PhysRevLett.92.211101]. For our highest resolution and most accurate gauge, we estimate the duration of this data set's last orbit to be approximately 59MADM.

Accurate evolutions of orbiting black-hole binaries without excision

We present a new algorithm for evolving orbiting black-hole binaries that does not require excision or a corotating shift. Our algorithm is based on a novel technique to handle the singular puncture conformal factor. This system, based on the Baumgarte-Shapiro-Shibata-Nakamura formulation of Einstein's equations, when used with a "precollapsed" initial lapse, is nonsingular at the start of the evolution and remains nonsingular and stable provided that a good choice is made for the gauge. As a test case, we use this technique to fully evolve orbiting black-hole binaries from near the innermost stable circular orbit regime. We show fourth-order convergence of waveforms and compute the radiated gravitational energy and angular momentum from the plunge. These results are in good agreement with those predicted by the Lazarus approach.

Collapse of magnetized hypermassive neutron stars in general relativity

Hypermassive neutron stars (HMNSs)--equilibrium configurations supported against collapse by rapid differential rotation--are possible transient remnants of binary neutron-star mergers. Using newly developed codes for magnetohydrodynamic simulations in dynamical spacetimes, we are able to track the evolution of a magnetized HMNS in full general relativity for the first time. We find that secular angular momentum transport due to magnetic braking and the magnetorotational instability results in the collapse of an HMNS to a rotating black hole, accompanied by a gravitational wave burst. The nascent black hole is surrounded by a hot, massive torus undergoing quasistationary accretion and a collimated magnetic field. This scenario suggests that HMNS collapse is a possible candidate for the central engine of short gamma-ray bursts.

Collapse of magnetized hypermassive neutron stars in general relativity

Hypermassive neutron stars (HMNSs)--equilibrium configurations supported against collapse by rapid differential rotation--are possible transient remnants of binary neutron-star mergers. Using newly developed codes for magnetohydrodynamic simulations in dynamical spacetimes, we are able to track the evolution of a magnetized HMNS in full general relativity for the first time. We find that secular angular momentum transport due to magnetic braking and the magnetorotational instability results in the collapse of an HMNS to a rotating black hole, accompanied by a gravitational wave burst. The nascent black hole is surrounded by a hot, massive torus undergoing quasistationary accretion and a collimated magnetic field. This scenario suggests that HMNS collapse is a possible candidate for the central engine of short gamma-ray bursts.

Collapse of magnetized hypermassive neutron stars in general relativity

Hypermassive neutron stars (HMNSs)--equilibrium configurations supported against collapse by rapid differential rotation--are possible transient remnants of binary neutron-star mergers. Using newly developed codes for magnetohydrodynamic simulations in dynamical spacetimes, we are able to track the evolution of a magnetized HMNS in full general relativity for the first time. We find that secular angular momentum transport due to magnetic braking and the magnetorotational instability results in the collapse of an HMNS to a rotating black hole, accompanied by a gravitational wave burst. The nascent black hole is surrounded by a hot, massive torus undergoing quasistationary accretion and a collimated magnetic field. This scenario suggests that HMNS collapse is a possible candidate for the central engine of short gamma-ray bursts.

Dynamical determination of the innermost stable circular orbit of binary neutron stars

We determine the innermost stable circular orbit (ISCO) of binary neutron stars (BNSs) by performing dynamical simulations in full general relativity. Evolving quasiequilibrium (QE) binaries that begin at different separations, we bracket the location of the ISCO by distinguishing stable circular orbits from unstable plunges. We study Gamma=2 polytropes of varying compactions in both corotational and irrotational equal-mass binaries. For corotational binaries, we find an ISCO orbital angular frequency somewhat smaller than that determined by applying turning-point methods to QE initial data. For the irrotational binaries, the initial data sequences terminate before reaching a turning point, but we find that the ISCO frequency is reached prior to the termination point. Our findings suggest that the ISCO frequency varies with compaction but does not depend strongly on the stellar spin.

Numerical Hydrodynamics in General Relativity

The current status of numerical solutions for the equations of ideal general relativistic hydrodynamics is reviewed. With respect to an earlier version of the article, the present update provides additional information on numerical schemes, and extends the discussion of astrophysical simulations in general relativistic hydrodynamics. Different formulations of the equations are presented, with special mention of conservative and hyperbolic formulations well-adapted to advanced numerical methods. A large sample of available numerical schemes is discussed, paying particular attention to solution procedures based on schemes exploiting the characteristic structure of the equations through linearized Riemann solvers. A comprehensive summary of astrophysical simulations in strong gravitational fields is presented. These include gravitational collapse, accretion onto black holes, and hydrodynamical evolutions of neutron stars. The material contained in these sections highlights the numerical challenges of various representative simulations. It also follows, to some extent, the chronological development of the field, concerning advances on the formulation of the gravitational field and hydrodynamic equations and the numerical methodology designed to solve them.

ELECTRONIC SUPPLEMENTARY MATERIAL

Supplementary material is available for this article at 10.12942/lrr-2003-4.

Rotating Stars in Relativity

Rotating relativistic stars have been studied extensively in recent years, both theoretically and observationally, because of the information they might yield about the equation of state of matter at extremely high densities and because they are considered to be promising sources of gravitational waves. The latest theoretical understanding of rotating stars in relativity is reviewed in this updated article. The sections on the equilibrium properties and on the nonaxisymmetric instabilities in f-modes and r-modes have been updated and several new sections have been added on analytic solutions for the exterior spacetime, rotating stars in LMXBs, rotating strange stars, and on rotating stars in numerical relativity.

ELECTRONIC SUPPLEMENTARY MATERIAL

Supplementary material is available for this article at 10.12942/lrr-2003-3.

3D grazing collision of two black holes

We present results for two colliding black holes (BHs), with angular momentum, spin, and unequal mass. For the first time, gravitational waveforms are computed for a grazing collision from a full 3D numerical evolution. The collision can be followed through the merger to form a single BH, and through part of the ringdown period of the final BH. The apparent horizon is tracked and studied, and physical parameters, such as the mass of the final BH, are computed. The total energy radiated in gravitational waves is shown to be consistent with the total initial mass of the spacetime and the apparent horizon mass of the final BH.

Computational Cosmology: From the Early Universe to the Large Scale Structure

In order to account for the observable Universe, any comprehensive theory or model of cosmology must draw from many disciplines of physics, including gauge theories of strong and weak interactions, the hydrodynamics and microphysics of baryonic matter, electromagnetic fields, and spacetime curvature, for example. Although it is difficult to incorporate all these physical elements into a single complete model of our Universe, advances in computing methods and technologies have contributed significantly towards our understanding of cosmological models, the Universe, and astrophysical processes within them. A sample of numerical calculations (and numerical methods applied to specific issues in cosmology are reviewed in this article: from the Big Bang singularity dynamics to the fundamental interactions of gravitational waves; from the quark-hadron phase transition to the large scale structure of the Universe. The emphasis, although not exclusively, is on those calculations designed to test different models of cosmology against the observed Universe.

Numerical Hydrodynamics in General Relativity

The current status of numerical solutions for the equations of ideal general relativistic hydrodynamics is reviewed. Different formulations of the equations are presented, with special mention of conservative and hyperbolic formulations well-adapted to advanced numerical methods. A representative sample of available numerical schemes is discussed and particular emphasis is paid to solution procedures based on schemes exploiting the characteristic structure of the equations through linearized Riemann solvers. A comprehensive summary of relevant astrophysical simulations in strong gravitational fields, including gravitational collapse, accretion onto black holes and evolution of neutron stars, is also presented.

ELECTRONIC SUPPLEMENTARY MATERIAL

Supplementary material is available for this article at 10.12942/lrr-2000-2.