Toward stable 3D numerical evolutions of black-hole spacetimes

Three dimensional (3D) numerical evolutions of static black holes with excision are presented. These evolutions extend to about 8000M, where M is the mass of the black hole. This degree of stability is achieved by using growth-rate estimates to guide the fine tuning of the parameters in a multi-parameter family of symmetric hyperbolic representations of the Einstein evolution equations. These evolutions were performed using a fixed gauge in order to separate the intrinsic stability of the evolution equations from the effects of stability-enhancing gauge choices.Comment: 4 pages, 5 figures. To appear in Phys. Rev. D. Minor additions to text for clarification. Added short paragraph about inner boundary dependenc

Using Full Information When Computing Modes of Post-Newtonian Waveforms From Inspiralling Compact Binaries in Circular Orbit

The increasing sophistication and accuracy of numerical simulations of compact binaries (especially binary black holes) presents the opportunity to test the regime in which post-Newtonian (PN) predictions for the emitted gravitational waves are accurate. In order to confront numerical results with those of post-Newtonian theory, it is convenient to compare multipolar decompositions of the two waveforms. It is pointed out here that the individual modes can be computed to higher post-Newtonian order by examining the radiative multipole moments of the system, rather than by decomposing the 2.5PN polarization waveforms. In particular, the dominant (l = 2, m = 2) mode can be computed to 3PN order. Individual modes are computed to as high a post-Newtonian order as possible given previous post-Newtonian results.Comment: 15 page

Numerical simulations of neutron star-black hole binaries in the near-equal-mass regime

Simulations of neutron star-black hole (NSBH) binaries generally consider black holes with masses in the range $(5-10)M_\odot$, where we expect to find most stellar mass black holes. The existence of lower mass black holes, however, cannot be theoretically ruled out. Low-mass black holes in binary systems with a neutron star companion could mimic neutron star-neutron (NSNS) binaries, as they power similar gravitational wave (GW) and electromagnetic (EM) signals. To understand the differences and similarities between NSNS mergers and low-mass NSBH mergers, numerical simulations are required. Here, we perform a set of simulations of low-mass NSBH mergers, including systems compatible with GW170817. Our simulations use a composition and temperature dependent equation of state (DD2) and approximate neutrino transport, but no magnetic fields. We find that low-mass NSBH mergers produce remnant disks significantly less massive than previously expected, and consistent with the post-merger outflow mass inferred from GW170817 for moderately asymmetric mass ratio. The dynamical ejecta produced by systems compatible with GW170817 is negligible except if the mass ratio and black hole spin are at the edge of the allowed parameter space. That dynamical ejecta is cold, neutron-rich, and surprisingly slow for ejecta produced during the tidal disruption of a neutron star : $v\sim (0.1-0.15)c$. We also find that the final mass of the remnant black hole is consistent with existing analytical predictions, while the final spin of that black hole is noticeably larger than expected -- up to $\chi_{\rm BH}=0.84$ for our equal mass case

Gravitational radiation reaction in compact binary systems: Contribution of the quadrupole-monopole interaction

The radiation reaction in compact spinning binaries on eccentric orbits due to the quadrupole-monopole interaction is studied. This contribution is of second post-Newtonian order. As result of the precession of spins the magnitude $L$ of the orbital angular momentum is not conserved. Therefore a proper characterization of the perturbed radial motion is provided by the energy $E$ and angular average $\bar{L}$. As powerful computing tools, the generalized true and eccentric anomaly parametrizations are introduced. Then the secular losses in energy and magnitude of orbital angular momentum together with the secular evolution of the relative orientations of the orbital angular momentum and spins are found for eccentric orbits by use of the residue theorem. The circular orbit limit of the energy loss agrees with Poisson's earlier result.Comment: accepted for publication in Phys. Rev.

Spin effects in gravitational radiation backreaction II. Finite mass effects

A convenient formalism for averaging the losses produced by gravitational radiation backreaction over one orbital period was developed in an earlier paper. In the present paper we generalize this formalism to include the case of a closed system composed from two bodies of comparable masses, one of them having the spin S. We employ the equations of motion given by Barker and O'Connell, where terms up to linear order in the spin (the spin-orbit interaction terms) are kept. To obtain the radiative losses up to terms linear in the spin, the equations of motion are taken to the same order. Then the magnitude L of the angular momentum L, the angle kappa subtended by S and L and the energy E are conserved. The analysis of the radial motion leads to a new parametrization of the orbit. From the instantaneous gravitational radiation losses computed by Kidder the leading terms and the spin-orbit terms are taken. Following Apostolatos, Cutler, Sussman and Thorne, the evolution of the vectors S and L in the momentary plane spanned by these vectors is separated from the evolution of the plane in space. The radiation-induced change in the spin is smaller than the leading-order spin terms in the momentary angular momentum loss. This enables us to compute the averaged losses in the constants of motion E, L and L_S=L cos kappa. In the latter, the radiative spin loss terms average to zero. An alternative description using the orbital elements a,e and kappa is given. The finite mass effects contribute terms, comparable in magnitude, to the basic, test-particle spin terms in the averaged losses.Comment: 12 pages, 1 figure, Phys.Rev.D15, March, 199

Gravity Waves, Chaos, and Spinning Compact Binaries

Spinning compact binaries are shown to be chaotic in the Post-Newtonian expansion of the two body system. Chaos by definition is the extreme sensitivity to initial conditions and a consequent inability to predict the outcome of the evolution. As a result, the spinning pair will have unpredictable gravitational waveforms during coalescence. This poses a challenge to future gravity wave observatories which rely on a match between the data and a theoretical template.Comment: Final version published in PR

Spin effects in gravitational radiation backreaction III. Compact binaries with two spinning components

The secular evolution of a spinning, massive binary system in eccentric orbit is analyzed, expanding and generalizing our previous treatments of the Lense-Thirring motion and the one-spin limit. The spin-orbit and spin-spin effects up to the 3/2 post-Newtonian order are considered, both in the equations of motion and in the radiative losses. The description of the orbit in terms of the true anomaly parametrization provides a simple averaging technique, based on the residue theorem, over eccentric orbits. The evolution equations of the angle variables characterizing the relative orientation of the spin and orbital angular momenta reveal a speed-up effect due to the eccentricity. The dissipative evolutions of the relevant dynamical and angular variables is presented in the form of a closed system of differential equations.Comment: 10 pages, 1 figur

Evolution systems for non-linear perturbations of background geometries

The formulation of the initial value problem for the Einstein equations is at the heart of obtaining interesting new solutions using numerical relativity and still very much under theoretical and applied scrutiny. We develop a specialised background geometry approach, for systems where there is non-trivial a priori knowledge about the spacetime under study. The background three-geometry and associated connection are used to express the ADM evolution equations in terms of physical non-linear deviations from that background. Expressing the equations in first order form leads naturally to a system closely linked to the Einstein-Christoffel system, introduced by Anderson and York, and sharing its hyperbolicity properties. We illustrate the drastic alteration of the source structure of the equations, and discuss why this is likely to be numerically advantageous.Comment: 12 pages, 3 figures, Revtex v3.0. Revised version to appear in Physical Review

Spin-spin effects in radiating compact binaries

The dynamics of a binary system with two spinning components on an eccentric orbit is studied, with the inclusion of the spin-spin interaction terms appearing at the second post-Newtonian order. A generalized true anomaly parametrization properly describes the radial component of the motion. The average over one radial period of the magnitude of the orbital angular momentum $\bar{L}$ is found to have no nonradiative secular change. All spin-spin terms in the secular radiative loss of the energy and magnitude of orbital angular momentum are given in terms of $\bar{L}$ and other constants of the motion. Among them, self-interaction spin effects are found, representing the second post-Newtonian correction to the 3/2 post-Newtonian order Lense-Thirring approximation.Comment: 12 pages, to appear in Phys. Rev.

3D simulations of linearized scalar fields in Kerr spacetime

We investigate the behavior of a dynamical scalar field on a fixed Kerr background in Kerr-Schild coordinates using a 3+1 dimensional spectral evolution code, and we measure the power-law tail decay that occurs at late times. We compare evolutions of initial data proportional to f(r) Y_lm(theta,phi) where Y_lm is a spherical harmonic and (r,theta,phi) are Kerr-Schild coordinates, to that of initial data proportional to f(r_BL) Y_lm(theta_BL,phi), where (r_BL,theta_BL) are Boyer-Lindquist coordinates. We find that although these two cases are initially almost identical, the evolution can be quite different at intermediate times; however, at late times the power-law decay rates are equal.Comment: 12 pages, 9 figures, revtex4. Major revision: added figures, added subsection on convergence, clarified discussion. To appear in Phys Rev