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 inspiral 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 index G=2 and G=5/3) 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 (the disc mass is approximately 0.2 solar masses) accretion torus around the spinning (Kerr) black hole. The rotation axis is clear of baryons (less than 1.e-5 solar masses within 10 degrees) 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.