Data For Binary Black Holes Research Articles

Constructing a boosted, spinning black hole in the damped harmonic gauge

The damped harmonic gauge is important for numerical relativity computations based on the generalized harmonic formulation of Einstein's equations, and is used to reduce coordinate distortions near binary black hole mergers. However, currently there is no prescription to construct quasiequilibrium binary black hole initial data in this gauge. Instead, initial data are typically constructed using a superposition of two boosted analytic single black hole solutions as free data in the solution of the constraint equations. Then, a smooth time-dependent gauge transformation is done early in the evolution to move into the damped harmonic gauge. Using this strategy to produce initial data in damped harmonic gauge would require the solution of a single black hole in this gauge, which is not known analytically. In this work we construct a single boosted, spinning, equilibrium BH in damped harmonic coordinates as a regular time-independent coordinate transformation from Kerr-Schild coordinates. To do this, we derive and solve a set of 4 coupled, nonlinear, elliptic equations for this transformation, with appropriate boundary conditions. This solution can now be used in the construction of damped harmonic initial data for binary black holes.

RIO: a new computational framework for accurate initial data of binary black holes

We present a computational framework (Rio) in the ADM 3+1 approach for numerical relativity. This work enables us to carry out high resolution calculations for initial data of two arbitrary black holes. We use the transverse conformal treatment, the Bowen-York and the puncture methods. For the numerical solution of the Hamiltonian constraint we use the domain decomposition and the spectral decomposition of Galerkin-Collocation. The nonlinear numerical code solves the set of equations for the spectral modes using the standard Newton-Raphson method, LU decomposition and Gaussian quadratures. We show the convergence of the Rio code. This code allows for easy deployment of large calculations. We show how the spin of one of the black holes is manifest in the conformal factor.

New code for equilibriums and quasiequilibrium initial data of compact objects. II. Convergence tests and comparisons of binary black hole initial data

COCAL is a code for computing equilibriums or quasiequilibrium initial data of single or binary compact objects based on finite difference methods. We present the results of supplementary convergence tests of COCAL code using time symmetric binary black hole data (Brill-Lindquist solution). Then, we compare the initial data of binary black holes on the conformally flat spatial slice obtained from COCAL and KADATH, where KADATH is a library for solving a wide class of problems in theoretical physics including relativistic compact objects with spectral methods. Data calculated from the two codes converge nicely towards each other, for close as well as largely separated circular orbits of binary black holes. Finally, as an example, a sequence of equal mass binary black hole initial data with corotating spins is calculated and compared with data in the literature.

New code for equilibriums and quasiequilibrium initial data of compact objects

We present a new code, named COCAL - Compact Object CALculator, for the computation of equilibriums and quasiequilibrium initial data sets of single or binary compact objects of all kinds. In the cocal code, those solutions are calculated on one or multiple spherical coordinate patches covering the initial hypersurface up to the asymptotic region. The numerical method used to solve field equations written in elliptic form is an adaptation of self-consistent field iterations in which Green's integral formula is computed using multipole expansions and standard finite difference schemes. We extended the method so that it can be used on a computational domain with excised regions for a black hole and a binary companion. Green's functions are constructed for various types of boundary conditions imposed at the surface of the excised regions for black holes. The numerical methods used in cocal are chosen to make the code simpler than any other recent initial data codes, accepting the second order accuracy for the finite difference schemes. We perform convergence tests for time symmetric single black hole data on a single coordinate patch, and binary black hole data on multiple patches. Then, we apply the code to obtain spatially conformally flat binary black hole initial data using boundary conditions including the one based on the existence of equilibrium apparent horizons.

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.

Time-symmetric initial data for binary black holes in numerical relativity

We look for physically realistic initial data in numerical relativity which are in agreement with post-Newtonian approximations. We propose a particular solution of the time-symmetric constraint equation, appropriate to two momentarily static black holes, in the form of a conformal decomposition of the spatial metric. This solution is isometric to the post-Newtonian metric up to the 2PN order. It represents a non-linear deformation of the solution of Brill and Lindquist, i.e. an asymptotically flat region is connected to two asymptotically flat (in a certain weak sense) sheets, that are the images of the two singularities through appropriate inversion transformations. The total ADM mass M as well as the individual masses m_1 and m_2 (when they exist) are computed by surface integrals performed at infinity. Using second order perturbation theory on the Brill-Lindquist background, we prove that the binary's interacting mass-energy M-m_1-m_2 is well-defined at the 2PN order and in agreement with the known post-Newtonian result.

Can a combination of the conformal thin-sandwich and puncture methods yield binary black hole solutions in quasiequilibrium?

We consider combining two important methods for constructing quasiequilibrium initial data for binary black holes: the conformal thin-sandwich formalism and the puncture method. The former seeks to enforce stationarity in the conformal three-metric and the latter attempts to avoid internal boundaries, like minimal surfaces or apparent horizons. We show that these two methods make partially conflicting requirements on the boundary conditions that determine the time slices. In particular, it does not seem possible to construct slices that are quasistationary and that avoid physical singularities while simultaneously are connected by an everywhere positive lapse function, a condition which must be obtained if internal boundaries are to be avoided. Some relaxation of these conflicting requirements may yield a soluble system, but some of the advantages that were sought in combining these approaches will be lost.

Modeling gravitational radiation from coalescing binary black holes

With the goal of bringing theory, particularly numerical relativity, to bear on an astrophysical problem of critical interest to gravitational wave observers we introduce a model for coalescence radiation from binary black hole systems. We build our model using the Lazarus approach, a technique that bridges far and close limit approaches with full numerical relativity to solve Einstein equations applied in the truly nonlinear dynamical regime. We specifically study the post-orbital radiation from a system of equal-mass non-spinning black holes, deriving waveforms which indicate strongly circularly polarized radiation of roughly 3% of the system's total energy and 12% of its total angular momentum in just a few cycles. To support this result we first establish the reliability of the late-time part of our model, including the numerical relativity and close-limit components, with a thorough study of waveforms from a sequence of black hole configurations that varies from previously treated head-on collisions to a representative target for ``ISCO'' data corresponding to the end of the inspiral period. We then complete our model with a simple treatment for the early part of the spacetime based on a standard family of initial data for binary black holes in circular orbit. A detailed analysis shows strong robustness in the results as the initial separation of the black holes is increased from $5.0$ to $7.8M$ supporting our waveforms as a suitable basic description of the astrophysical radiation from this system. Finally, a simple fitting of the plunge waveforms is introduced as a first attempt to facilitate the task of analyzing data from gravitational wave detectors.

Perturbative evolution of nonlinear initial data for binary black holes: Zerilli versus Teukolsky equation

We consider the problem of evolving nonlinear initial data in the close limit regime. For exact Misner initial data (two equal mass black holes initially at rest), metric perturbations evolved via the Zerilli equation suffer from a premature breakdown (at a proper separation of the holes $L/M\ensuremath{\approx}2.2)$ while we find that the exact Weyl scalar ${\ensuremath{\psi}}_{4}$ evolved via the Teukolsky equation keeps very good agreement with the full numerical results up to $L/M\ensuremath{\approx}3.5.$ Metric and curvature perturbations of nonrotating black holes are equivalent to first perturbative order, but the Moncrief waveform in the former case and the Weyl scalar ${\ensuremath{\psi}}_{4}$ in the latter differ when nonlinearities are present. We argue that this inequivalent behavior holds for a wider class of conformally flat initial data.