Newtonian hydrodynamics of the coalescence of black holes with neutron stars - IV. Irrotational binaries with a soft equation of state

Abstract

We present the results of three-dimensional hydrodynamical simulations of the final stages of in-spiral in a black hole–neutron star binary, when the separation is comparable to the stellar radius. We use a Newtonian smooth particle hydrodynamics (SPH) code to model the evolution of the system, and take the neutron star to be a polytrope with a soft (adiabatic indices and equation of state and the black hole to be a Newtonian point mass. The only non-Newtonian effect we include is a gravitational radiation back reaction force, computed in the quadrupole approximation for point masses. We use irrotational binaries as initial conditions for our dynamical simulations, which are begun when the system is on the verge of initiating mass transfer and followed for approximately 23 ms. For all the cases studied we find that the star is disrupted on a dynamical time-scale, and forms a massive accretion torus around the spinning (Kerr) black hole. The rotation axis is clear of baryons (less than 10−5 M⊙ within 10°) to an extent that would not preclude the formation of a relativistic fireball capable of powering a cosmological gamma-ray burst. Some mass (the specific amount is sensitive to the stiffness of the equation of state) may be dynamically ejected from the system during the coalescence and could undergo r-process nucleosynthesis. We calculate the waveforms, luminosities and energy spectra of the gravitational radiation signal, and show how they reflect the global outcome of the coalescence process.