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

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

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.

Solving the Darwin problem in the first post-Newtonian approximation of general relativity

We analytically calculate the equilibrium sequence of the corotating binary stars of incompressible fluid in the first post-Newtonian(PN) approximation of general relativity. By calculating the total energy and total angular momentum of the system as a function of the orbital separation, we investigate the innermost stable circular orbit for corotating binary(we call it ISCCO). It is found that by the first PN effect, the orbital separation of the binary at the ISCCO becomes small with increase of the compactness of each star, and as a result, the orbital angular velocity at the ISCCO increases. These behaviors agree with previous numerical works.Comment: 33 pages, revtex, 4 figures(eps), accepted for publication in Phys. Rev.

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

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.

Ruling Out Chaos in Compact Binary Systems

We investigate the orbits of compact binary systems during the final inspiral period before coalescence by integrating numerically the second-order post-Newtonian equations of motion. We include spin-orbit and spin-spin coupling terms, which, according to a recent study by Levin [J. Levin, Phys. Rev. Lett. 84, 3515 (2000)], may cause the orbits to become chaotic. To examine this claim, we study the divergence of initially nearby phase-space trajectories and attempt to measure the Lyapunov exponent gamma. Even for systems with maximally spinning objects and large spin-orbit misalignment angles, we find no chaotic behavior. For all the systems we consider, we can place a strict lower limit on the divergence time t_L=1/gamma that is many times greater than the typical inspiral time, suggesting that chaos should not adversely affect the detection of inspiral events by upcoming gravitational-wave detectors.Comment: 8 pages, 4 figures, submitted to Phys. Rev. Let

Up-down instability of binary black holes in numerical relativity

Binary black holes with spins that are aligned with the orbital angular momentum do not precess. However, post-Newtonian calculations predict that "up-down" binaries, in which the spin of the heavier (lighter) black hole is aligned (antialigned) with the orbital angular momentum, are unstable when the spins are slightly perturbed from perfect alignment. This instability provides a possible mechanism for the formation of precessing binaries in environments where sources are preferentially formed with (anti) aligned spins. In this paper, we present the first full numerical relativity simulations capturing this instability. These simulations span $\sim 100$ orbits and $\sim 3$-$5$ precession cycles before merger, making them some of the longest numerical relativity simulations to date. Initialized with a small perturbation of $1^{\circ}$-$10^{\circ}$, the instability causes a dramatic growth of the spin misalignments, which can reach $\sim 90^{\circ}$ near merger. We show that this leaves a strong imprint on the subdominant modes of the gravitational wave signal, which can potentially be used to distinguish up-down binaries from other sources. Finally, we show that post-Newtonian and effective-one-body approximants are able to reproduce the unstable dynamics of up-down binaries extracted from numerical relativity

Solving Einstein's Equations With Dual Coordinate Frames

A method is introduced for solving Einstein's equations using two distinct coordinate systems. The coordinate basis vectors associated with one system are used to project out components of the metric and other fields, in analogy with the way fields are projected onto an orthonormal tetrad basis. These field components are then determined as functions of a second independent coordinate system. The transformation to the second coordinate system can be thought of as a mapping from the original ``inertial'' coordinate system to the computational domain. This dual-coordinate method is used to perform stable numerical evolutions of a black-hole spacetime using the generalized harmonic form of Einstein's equations in coordinates that rotate with respect to the inertial frame at infinity; such evolutions are found to be generically unstable using a single rotating coordinate frame. The dual-coordinate method is also used here to evolve binary black-hole spacetimes for several orbits. The great flexibility of this method allows comoving coordinates to be adjusted with a feedback control system that keeps the excision boundaries of the holes within their respective apparent horizons.Comment: Updated to agree with published versio

2.5PN kick from black-hole binaries in circular orbit: Nonspinning case

Using the Multipolar post-Minskowskian formalism, we compute the linear momentum flux from black-hole binaries in circular orbits and having no spins. The total linear momentum flux contains various types of instantaneous (which are functions of the retarded time) and hereditary (which depends on the dynamics of the binary in the past) terms both of which are analytically computed. In addition to the inspiral contribution, we use a simple model of plunge to compute the kick or recoil accumulated during this phase.Comment: To appear in Proceedings of "Relativity and Gravitation - 100 Years after Einstein in Prague" Ed. J. Bicak (2013