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

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

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.

Gauge conditions for binary black hole puncture data based on an approximate helical Killing vector

We show that puncture data for quasicircular binary black hole orbits allow a special gauge choice that realizes some of the necessary conditions for the existence of an approximate helical Killing vector field. Introducing free parameters for the lapse at the punctures we can satisfy the condition that the Komar and ADM mass agree at spatial infinity. Several other conditions for an approximate Killing vector are then automatically satisfied, and the 3-metric evolves on a timescale smaller than the orbital timescale. The time derivative of the extrinsic curvature however remains significant. Nevertheless, quasicircular puncture data are not as far from possessing a helical Killing vector as one might have expected.Comment: 11 pages, 6 figures, 2 table

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

Comparing initial-data sets for binary black holes

We compare the results of constructing binary black hole initial data with three different decompositions of the constraint equations of general relativity. For each decomposition we compute the initial data using a superposition of two Kerr-Schild black holes to fix the freely specifiable data. We find that these initial-data sets differ significantly, with the ADM energy varying by as much as 5% of the total mass. We find that all initial-data sets currently used for evolutions might contain unphysical gravitational radiation of the order of several percent of the total mass. This is comparable to the amount of gravitational-wave energy observed during the evolved collision. More astrophysically realistic initial data will require more careful choices of the freely specifiable data and boundary conditions for both the metric and extrinsic curvature. However, we find that the choice of extrinsic curvature affects the resulting data sets more strongly than the choice of conformal metric.Comment: 18 pages, 12 figures, accepted for publication in Phys. Rev.

Comparing Criteria for Circular Orbits in General Relativity

We study a simple analytic solution to Einstein's field equations describing a thin spherical shell consisting of collisionless particles in circular orbit. We then apply two independent criteria for the identification of circular orbits, which have recently been used in the numerical construction of binary black hole solutions, and find that both yield equivalent results. Our calculation illustrates these two criteria in a particularly transparent framework and provides further evidence that the deviations found in those numerical binary black hole solutions are not caused by the different criteria for circular orbits.Comment: 4 pages; to appear in PRD as a Brief Report; added and corrected reference

Conformal-thin-sandwich initial data for a single boosted or spinning black hole puncture

Sequences of initial-data sets representing binary black holes in quasi-circular orbits have been used to calculate what may be interpreted as the innermost stable circular orbit. These sequences have been computed with two approaches. One method is based on the traditional conformal-transverse-traceless decomposition and locates quasi-circular orbits from the turning points in an effective potential. The second method uses a conformal-thin-sandwich decomposition and determines quasi-circular orbits by requiring the existence of an approximate helical Killing vector. Although the parameters defining the innermost stable circular orbit obtained from these two methods differ significantly, both approaches yield approximately the same initial data, as the separation of the binary system increases. To help understanding this agreement between data sets, we consider the case of initial data representing a single boosted or spinning black hole puncture of the Bowen-York type and show that the conformal-transverse-traceless and conformal-thin-sandwich methods yield identical data, both satisfying the conditions for the existence of an approximate Killing vector.Comment: 13 pages, 2 figure

Self-organization of (001) cubic crystal surfaces

Self-organization on crystal surface is studied as a two dimensional spinodal decomposition in presence of a surface stress. The elastic Green function is calculated for a $(001)$ cubic crystal surface taking into account the crystal anisotropy. Numerical calculations show that the phase separation is driven by the interplay between domain boundary energy and long range elastic interactions. At late stage of the phase separation process, a steady state appears with different nanometric patterns according to the surface coverage and the crystal elastic constants

Dynamical Gauge Conditions for the Einstein Evolution Equations

The Einstein evolution equations have been written in a number of symmetric hyperbolic forms when the gauge fields--the densitized lapse and the shift--are taken to be fixed functions of the coordinates. Extended systems of evolution equations are constructed here by adding the gauge degrees of freedom to the set of dynamical fields, thus forming symmetric hyperbolic systems for the combined evolution of the gravitational and the gauge fields. The associated characteristic speeds can be made causal (i.e. less than or equal to the speed of light) by adjusting 14 free parameters in these new systems. And 21 additional free parameters are available, for example to optimize the stability of numerical evolutions. The gauge evolution equations in these systems are generalizations of the ``K-driver'' and ``Gamma-driver'' conditions that have been used with some success in numerical black hole evolutions.Comment: New appendix on constraint evolution adde