A Linear-Nonlinear Formulation of Einstein Equations for the Two-Body Problem in General Relativity

A formulation of Einstein equations is presented that could yield advantages in the study of collisions of binary compact objects during regimes between linear-nonlinear transitions. The key idea behind this formulation is a separation of the dynamical variables into i) a fixed conformal 3-geometry, ii) a conformal factor possessing nonlinear dynamics and iii) transverse-traceless perturbations of the conformal 3-geometry.Comment: 7 pages, no figure

Irrotational binary neutron stars in quasiequilibrium

We report on numerical results from an independent formalism to describe the quasi-equilibrium structure of nonsynchronous binary neutron stars in general relativity. This is an important independent test of controversial numerical hydrodynamic simulations which suggested that nonsynchronous neutron stars in a close binary can experience compression prior to the last stable circular orbit. We show that, for compact enough stars the interior density increases slightly as irrotational binary neutron stars approach their last orbits. The magnitude of the effect, however, is much smaller than that reported in previous hydrodynamic simulations.Comment: 4 pages, 2 figures, revtex, accepted for publication in Phys. Rev.

Relativistic Models for Binary Neutron Stars with Arbitrary Spins

We introduce a new numerical scheme for solving the initial value problem for quasiequilibrium binary neutron stars allowing for arbitrary spins. The coupled Einstein field equations and equations of relativistic hydrodynamics are solved in the Wilson-Mathews conformal thin sandwich formalism. We construct sequences of circular-orbit binaries of varying separation, keeping the rest mass and circulation constant along each sequence. Solutions are presented for configurations obeying an n=1 polytropic equation of state and spinning parallel and antiparallel to the orbital angular momentum. We treat stars with moderate compaction ((m/R) = 0.14) and high compaction ((m/R) = 0.19). For all but the highest circulation sequences, the spins of the neutron stars increase as the binary separation decreases. Our zero-circulation cases approximate irrotational sequences, for which the spin angular frequencies of the stars increases by 13% (11%) of the orbital frequency for (m/R) = 0.14 ((m/R) = 0.19) by the time the innermost circular orbit is reached. In addition to leaving an imprint on the inspiral gravitational waveform, this spin effect is measurable in the electromagnetic signal if one of the stars is a pulsar visible from Earth.Comment: 21 pages, 14 figures. A few explanatory sentences added and some typos corrected. Accepted for publication in Phys. Rev.

Black Hole Boundary Conditions and Coordinate Conditions

This paper treats boundary conditions on black hole horizons for the full 3+1D Einstein equations. Following a number of authors, the apparent horizon is employed as the inner boundary on a space slice. It is emphasized that a further condition is necessary for the system to be well posed; the ``prescribed curvature conditions" are therefore proposed to complete the coordinate conditions at the black hole. These conditions lead to a system of two 2D elliptic differential equations on the inner boundary surface, which coexist nicely to the 3D equation for maximal slicing (or related slicing conditions). The overall 2D/3D system is argued to be well posed and globally well behaved. The importance of ``boundary conditions without boundary values" is emphasized. This paper is the first of a series. This revised version makes minor additions and corrections to the previous version.Comment: 13 pages LaTeX, revtex. No figure

Corotating and irrotational binary black holes in quasi-circular orbits

A complete formalism for constructing initial data representing black-hole binaries in quasi-equilibrium is developed. Radiation reaction prohibits, in general, true equilibrium binary configurations. However, when the timescale for orbital decay is much longer than the orbital period, a binary can be considered to be in quasi-equilibrium. If each black hole is assumed to be in quasi-equilibrium, then a complete set of boundary conditions for all initial data variables can be developed. These boundary conditions are applied on the apparent horizon of each black hole, and in fact force a specified surface to be an apparent horizon. A global assumption of quasi-equilibrium is also used to fix some of the freely specifiable pieces of the initial data and to uniquely fix the asymptotic boundary conditions. This formalism should allow for the construction of completely general quasi-equilibrium black hole binary initial data.Comment: 13 pages, no figures, revtex4; Content changed slightly to reflect fact that regularized shift solutions do satisfy the isometry boundary condition

Potential for ill-posedness in several 2nd-order formulations of the Einstein equations

Second-order formulations of the 3+1 Einstein equations obtained by eliminating the extrinsic curvature in terms of the time derivative of the metric are examined with the aim of establishing whether they are well posed, in cases of somewhat wide interest, such as ADM, BSSN and generalized Einstein-Christoffel. The criterion for well-posedness of second-order systems employed is due to Kreiss and Ortiz. By this criterion, none of the three cases are strongly hyperbolic, but some of them are weakly hyperbolic, which means that they may yet be well posed but only under very restrictive conditions for the terms of order lower than second in the equations (which are not studied here). As a result, intuitive transferences of the property of well-posedness from first-order reductions of the Einstein equations to their originating second-order versions are unwarranted if not false.Comment: v1:6 pages; v2:7 pages, discussion extended, to appear in Phys. Rev. D; v3: typos corrected, published versio

Relativistic dust disks and the Wilson-Mathews approach

Treating problems in full general relativity is highly complex and frequently approximate methods are employed to simplify the solution. We present comparative solutions of a infinitesimally thin relativistic, stationary, rigidly rotating disk obtained using the full equations and the approximate approach suggested by Wilson & Mathews. We find that the Wilson-Mathews method has about the same accuracy as the first post-Newtonian approximation.Comment: 4 Pages, 5 eps-figures, uses revtex.sty. Submitted to PR

On the Circular Orbit Approximation for Binary Compact Objects In General Relativity

One often-used approximation in the study of binary compact objects (i.e., black holes and neutron stars) in general relativity is the instantaneously circular orbit assumption. This approximation has been used extensively, from the calculation of innermost circular orbits to the construction of initial data for numerical relativity calculations. While this assumption is inconsistent with generic general relativistic astrophysical inspiral phenomena where the dissipative effects of gravitational radiation cause the separation of the compact objects to decrease in time, it is usually argued that the timescale of this dissipation is much longer than the orbital timescale so that the approximation of circular orbits is valid. Here, we quantitatively analyze this approximation using a post-Newtonian approach that includes terms up to order ({Gm/(rc^2)})^{9/2} for non-spinning particles. By calculating the evolution of equal mass black hole / black hole binary systems starting with circular orbit configurations and comparing them to the more astrophysically relevant quasicircular solutions, we show that a minimum initial separation corresponding to at least 6 (3.5) orbits before plunge is required in order to bound the detection event loss rate in gravitational wave detectors to < 5% (20%). In addition, we show that the detection event loss rate is > 95% for a range of initial separations that include all modern calculations of the innermost circular orbit (ICO).Comment: 10 pages, 12 figures, revtex

Improved numerical stability of stationary black hole evolution calculations

We experiment with modifications of the BSSN form of the Einstein field equations (a reformulation of the ADM equations) and demonstrate how these modifications affect the stability of numerical black hole evolution calculations. We use excision to evolve both non-rotating and rotating Kerr-Schild black holes in octant and equatorial symmetry, and without any symmetry assumptions, and obtain accurate and stable simulations for specific angular momenta J/M of up to about 0.9M.Comment: 13 pages, 11 figures, 1 typo in Eq. (20) correcte

Illustrating Stability Properties of Numerical Relativity in Electrodynamics

We show that a reformulation of the ADM equations in general relativity, which has dramatically improved the stability properties of numerical implementations, has a direct analogue in classical electrodynamics. We numerically integrate both the original and the revised versions of Maxwell's equations, and show that their distinct numerical behavior reflects the properties found in linearized general relativity. Our results shed further light on the stability properties of general relativity, illustrate them in a very transparent context, and may provide a useful framework for further improvement of numerical schemes.Comment: 5 pages, 2 figures, to be published as Brief Report in Physical Review