Black Hole Initial Data Research Articles

Marginally Trapped Tubes Generated from Nonlinear Scalar Field Initial Data

We show that the maximal future development of asymptotically flat spherically symmetric black hole initial data for a self-gravitating nonlinear scalar field, also called a Higgs field, contains a connected, achronal marginally trapped tube which is asymptotic to the event horizon of the black hole, provided the initial data is sufficiently small and decays like O(r^{-1/2}), and the potential function V is nonnegative with bounded second derivative. This result can be loosely interpreted as a statement about the stability of `nice' asymptotic behavior of marginally trapped tubes under certain small perturbations of Schwarzschild.

Equilibrium initial data for moving puncture simulations: the stationary 1 + log slicing

We discuss a ‘stationary 1 + log’ slicing condition for the construction of solutions to Einstein's constraint equations. For stationary spacetimes, these initial data give a stationary foliation when evolved with ‘moving puncture’ gauge conditions that are often used in black hole evolutions. The resulting slicing is time independent and agrees with the slicing generated by being dragged along a timelike Killing vector of the spacetime. When these initial data are evolved with moving puncture gauge conditions, numerical errors arising from coordinate evolution should be minimized. While these properties appear very promising, suggesting that this slicing condition should be an attractive alternative to, for example, maximal slicing, we demonstrate in this paper that solutions can be constructed only for a small class of problems. For binary black hole initial data, in particular, it is often assumed that there exists an approximate helical Killing vector that generates the binary's orbit. We show that 1 + log slices that are stationary with respect to such a helical Killing vector cannot be asymptotically flat, unless the spacetime possesses an additional axial Killing vector.

Extreme Bowen–York initial data

The Bowen–York family of spinning black hole initial data depends essentially on one, positive, free parameter. The extreme limit corresponds to making this parameter equal to zero. This choice represents a singular limit for the constraint equations. We prove that in this limit a new solution of the constraint equations is obtained. These initial data have similar properties to the extreme Kerr and Reissner–Nordström black hole initial data. In particular, in this limit one of the asymptotic ends changes from asymptotically flat to cylindrical. The existence proof is constructive, we actually show that a sequence of Bowen–York data converges to the extreme solution.

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.

Black hole initial data with a horizon of prescribed geometry

The purpose of this work is to construct asymptotically flat, time symmetric initial data with an apparent horizon of prescribed intrinsic geometry. To do this, we use the parabolic partial differential equation for prescribing scalar curvature. In this equation the horizon geometry is contained within the freely specifiable part of the metric. This contrasts with the conformal method in which the geometry of the horizon can only be specified up to a conformal factor.

Measuring eccentricity in binary black-hole initial data

Initial data for evolving black-hole binaries can be constructed via many techniques, and can represent a wide range of physical scenarios. However, because of the way that different schemes parametrize the physical aspects of a configuration, it is not always clear what a given set of initial data actually represents. This is especially important for quasiequilibrium data constructed using the conformal thin-sandwich approach. Most initial-data studies have focused on identifying data sets that represent binaries in quasicircular orbits. In this paper, we consider initial-data sets representing equal-mass black-hole binaries in eccentric orbits. We will show that effective-potential techniques can be used to calibrate initial data for black-hole binaries in eccentric orbits. We will also examine several different approaches, including post-Newtonian diagnostics, for measuring the eccentricity of an orbit. Finally, we propose the use of the ``Komar-mass difference'' as a useful, invariant means of parametrizing the eccentricity of relativistic orbits.

Filling the holes: Evolving excised binary black hole initial data with puncture techniques

We follow the inspiral and merger of equal-mass black holes (BHs) by the moving puncture technique and demonstrate that both the exterior solution and the asymptotic gravitational waveforms are unchanged when the initial interior solution is replaced by constraint-violating junk initial data. We apply this result to evolve conformal thin-sandwich (CTS) binary BH initial data by filling their excised interiors with arbitrary, but smooth, initial data and evolving with standard puncture gauge choices. The waveforms generated for both puncture and filled-CTS initial data are remarkably similar, and there are only minor differences between irrotational and corotational CTS BH binaries. Even the interior solutions appear to evolve to the same constraint-satisfying solution at late times, independent of the initial data.

Energy-minimizing two black holes initial data

An attempt to construct the ‘‘ground state’’ vacuum initial data for the gravitational field surrounding two black holes is presented. The ground state is defined as the gravitational initial data minimizing the ADM mass within the class of data for which the masses of the holes and their distance are fixed. To parametrize different geometric arrangements of the two holes (and, therefore, their distance), we use an appropriately chosen scale factor. A method for analyzing the variations of the ADM mass and the masses (areas) of the horizons in terms of gravitational degrees of freedom is proposed. The Misner initial data are analyzed in this context: it is shown that they do not minimize the ADM mass.

Reducing orbital eccentricity in binary black hole simulations

Binary black hole simulations starting from quasi-circular (i.e., zero radial velocity) initial data have orbits with small but nonzero orbital eccentricities. In this paper, the quasi-equilibrium initial-data method is extended to allow nonzero radial velocities to be specified in binary black hole initial data. New low-eccentricity initial data are obtained by adjusting the orbital frequency and radial velocities to minimize the orbital eccentricity, and the resulting (∼5 orbit) evolutions are compared with those of quasi-circular initial data. Evolutions of the quasi-circular data clearly show eccentric orbits, with eccentricity that decays over time. The precise decay rate depends on the definition of eccentricity; if defined in terms of variations in the orbital frequency, the decay rate agrees well with the prediction of Peters (1964 Phys. Rev. 136 1224–32). The gravitational waveforms, which contain ∼8 cycles in the dominant l = m = 2 mode, are largely unaffected by the eccentricity of the quasi-circular initial data. The overlap between the dominant mode in the quasi-circular evolution and the same mode in the low-eccentricity evolution is about 0.99.

Binary black hole initial data from matched asymptotic expansions

We present an approximate metric for a binary black-hole spacetime to construct initial data for numerical relativity. This metric is obtained by asymptotically matching a post-Newtonian metric for a binary system to a perturbed Schwarzschild metric for each hole. In the inner zone near each hole, the metric is given by the Schwarzschild solution plus a quadrupolar perturbation corresponding to an external tidal gravitational field. In the near zone, well outside each black hole but less than a reduced wavelength from the center of mass of the binary, the metric is given by a post-Newtonian expansion including the lowest-order deviations from flat spacetime. When the near zone overlaps each inner zone in a buffer zone, the post-Newtonian and perturbed Schwarzschild metrics can be asymptotically matched to each other. By demanding matching (over a 4-volume in the buffer zone) rather than patching (choosing a particular 2-surface in the buffer zone), we guarantee that the errors are small in all zones. The resulting piecewise metric is made formally C{sup {infinity}} with smooth transition functions so as to obtain the finite extrinsic curvature of a 3-slice. In addition to the metric and extrinsic curvature, we present explicit results for the lapse and themore » shift, which can be used as initial data for numerical simulations. This initial data is not accurate all the way to the asymptotically flat ends inside each hole, and therefore must be used with evolution codes which employ black hole excision rather than puncture methods. This paper lays the foundations of a method that can be straightforwardly iterated to obtain initial data to higher perturbative order.« less

Relativistic three-body effects in black hole coalescence

Three-body interactions are expected to be common in globular clusters and in galactic cores hosting supermassive black holes. We consider an equal-mass binary black hole system in the presence of a third black hole. Using numerically generated binary black hole initial data sets, and first and second-order post-Newtonian (1PN and 2PN) techniques, we find that the presence of the third black hole has non-negligible relativistic effects on the location of the binary’s innermost stable circular orbit (ISCO), and that these effects arise at 2PN order. For a stellar-mass black hole binary in orbit about a supermassive black hole, the massive black hole has stabilizing effects on the orbiting binary, leading to an increase in merger time and a decrease of the terminal orbital frequency, and an amplification of the gravitational radiation emitted from the binary system by up to 6%.

Approximate initial data for binary black holes

We construct approximate analytical solutions to the constraint equations of general relativity for binary black holes of arbitrary mass ratio in quasicircular orbit. We adopt the puncture method to solve the constraint equations in the transverse-traceless decomposition and consider perturbations of Schwarzschild black holes caused by boosts and the presence of a binary companion. A superposition of these two perturbations then yields approximate, but fully analytic binary black hole initial data that are accurate to first order in the inverse of the binary separation and the square of the black holes' momenta.

Circular orbits and spin in black-hole initial data

The construction of initial data for black-hole binaries usually involves the choice of free parameters that define the spins of the black holes and essentially the eccentricity of the orbit. Such parameters must be chosen carefully to yield initial data with the desired physical properties. In this paper, we examine these choices in detail for the quasiequilibrium method coupled to apparent-horizon/quasiequilibrium boundary conditions. First, we compare two independent criteria for choosing the orbital frequency, the "Komar-mass condition" and the "effective-potential method," and find excellent agreement. Second, we implement quasi-local measures of the spin of the individual holes, calibrate these with corotating binaries, and revisit the construction of non-spinning black hole binaries. Higher-order effects, beyond those considered in earlier work, turn out to be important. Without those, supposedly non-spinning black holes have appreciable quasi-local spin; furthermore, the Komar-mass condition and effective potential method agree only when these higher-order effects are taken into account. We compute a new sequence of quasi-circular orbits for non-spinning black-hole binaries, and determine the innermost stable circular orbit of this sequence.

Conformal thin-sandwich puncture initial data for boosted black holes

We apply the puncture approach to conformal thin-sandwich black-hole initial data. We solve numerically the conformal thin-sandwich puncture (CTSP) equations for a single black hole with nonzero linear momentum. We show that conformally-flat solutions for a boosted black hole have the same maximum gravitational radiation content as the corresponding Bowen-York solution in the conformal transverse-traceless decomposition. We find that the physical properties of these data are independent of the free slicing parameter.

Inner boundary conditions for black hole initial data derived from isolated horizons

We present a set of boundary conditions for solving the elliptic equations in the Initial Data Problem for space-times containing a black hole, together with a number of constraints to be satisfied by the otherwise freely specifiable standard parameters of the Conformal Thin Sandwich formulation. These conditions altogether are sufficient for the construction of a horizon that is instantaneously in equilibrium in the sense of the Isolated Horizons formalism. We then investigate the application of these conditions to the Initial Data Problem of binary black holes and discuss the relation of our analysis with other proposals that exist in the literature.

Distorted black hole initial data using the puncture method

We solve for single distorted black hole initial data using the puncture method, where the Hamiltonian constraint is written as an elliptic equation in R^3 for the nonsingular part of the metric conformal factor. With this approach we can generate isometric and non--isometric black hole data. For the isometric case, our data are directly comparable to those obtained by Bernstein et al., who impose isometry boundary conditions at the black hole throat. Our numerical simulations are performed using a parallel multigrid elliptic equation solver with adaptive mesh refinement. Mesh refinement allows us to use high resolution around the black hole while keeping the grid boundaries far away in the asymptotic region.

Excision boundary conditions for black-hole initial data

We define and extensively test a set of boundary conditions that can be applied at black hole excision surfaces when the Hamiltonian and momentum constraints of general relativity are solved within the conformal thin-sandwich formalism. These boundary conditions have been designed to result in black holes that are in quasiequilibrium and are completely general in the sense that they can be applied with any conformal three-geometry and slicing condition. Furthermore, we show that they retain precisely the freedom to specify an arbitrary spin on each black hole. Interestingly, we have been unable to find a boundary condition on the lapse that can be derived from a quasiequilibrium condition. Rather, we find evidence that the lapse boundary condition is part of the initial temporal gauge choice. To test these boundary conditions, we have extensively explored the case of a single black hole and the case of a binary system of equal-mass black holes, including the computation of quasi-circular orbits and the determination of the inner-most stable circular orbit. Our tests show that the boundary conditions work well.

Publisher's Note: Quasi-equilibrium binary black hole initial data for dynamical evolutions [Phys. Rev. D70, 084033 (2004)

We present a formalism for constructing quasi-equilibrium binary black hole initial data suitable for numerical evolution. We construct quasi-equilibrium models by imposing an approximate helical Killing symmetry appropriate for quasi-circular orbits. We use the sum of two Kerr-Schild metrics as our background metric, thereby improving on conformally flat backgrounds that do not accommodate rotating black holes and providing a horizon-penetrating lapse, convenient for implementing black hole excision. We set inner boundary conditions at an excision radius well inside the apparent horizon and construct these boundary conditions to incorporate the quasi-equilibrium condition and recover the solution for isolated black holes in the limit of large separation. We use our formalism both to generate initial data for binary black hole evolutions and to construct a crude quasi-equilibrium, inspiral sequence for binary black holes of fixed irreducible mass.

Black ring formation in particle systems

It is known that the formation of apparent horizons with nonspherical topology is possible in higher-dimensional spacetimes. One of these is the black ring horizon with ${S}^{1}\ifmmode\times\else\texttimes\fi{}{S}^{D\ensuremath{-}3}$ topology where $D$ is the spacetime dimension number. In this paper, we investigate the black ring horizon formation in systems with $n$-particles. We analyze two kinds of system: the high-energy $n$-particle system and the momentarily-static $n$-black hole initial data. In the high-energy particle system, we prove that the black ring horizon does not exist at the instant of collision for any $n$. But there remains a possibility that the black ring forms after the collision and this result is not sufficient. Because calculating the metric of this system after the collision is difficult, we consider the momentarily-static $n$-black hole initial data that can be regarded as a simplified $n$-particle model and numerically solve the black ring horizon that surrounds all the particles. Our results show that there is the minimum particle number that is necessary for the black ring formation and this number depends on $D$. Although many particle number is required in five-dimensions, $n=4$ is sufficient for the black ring formation in the $D\ensuremath{\ge}7$ cases. The black ring formation becomes easier for larger $D$. We provide a plausible physical interpretation of our results and discuss the validity of Ida and Nakao's conjecture for the horizon formation in higher-dimensions. Finally we briefly discuss the probable methods of producing the black rings in accelerators.

Quasi-equilibrium binary black hole initial data for dynamical evolutions

We present a formalism for constructing quasiequilibrium binary black hole initial data suitable for numerical evolution. We construct quasiequilibrium models by imposing an approximate helical Killing symmetry appropriate for quasicircular orbits. We use the sum of two Kerr-Schild metrics as our background metric, thereby improving on conformal flat backgrounds that do not accommodate rotating black holes and providing a horizon-penetrating lapse, convenient for implementing black hole excision. We set inner boundary conditions at an excision radius well inside the apparent horizon and construct these boundary conditions to incorporate the quasiequilibrium condition and recover the solution for isolated black holes in the limit of large separation. We use our formalism both to generate initial data for binary black hole evolutions and to construct a crude quasiequilibrium, inspiral sequence for binary black holes of fixed irreducible mass.