Newtonian hydrodynamics of the coalescence of black holes with neutron stars — III. Irrotational binaries with a stiff equation of state

Abstract

We present a numerical study of the hydrodynamics in the final stages of inspiral of a black hole-neutron star binary, when the binary separation becomes comparable to the stellar radius. We use a Newtonian three-dimensional Smooth Particle Hydrodynamics (SPH) code, and model the neutron star with a soft (adiabatic index Gamma=5/3) polytropic equation of state and the black hole as a Newtonian point mass which accretes matter via an absorbing boundary at the Schwarzschild radius. Our initial conditions correspond to tidally locked binaries in equilibrium, and we have explored configurations with different values of the mass ratio q=Mns/Mbh, ranging from q=1 to q=0.1. The dynamical evolution is followed for approximately 23 ms, and in every case studied here we find that the neutron star is tidally disrupted on a dynamical timescale, forming a dense torus around the black hole that contains a few tenths of a solar mass. A nearly baryon-free axis is present in the system throughout the coalescence, and only modest beaming of a fireball that could give rise to a gamma-ray burst would be sufficient to avoid excessive baryon contamination. We find that some mass (on the order of 0.01 to 0.001 solar masses) may be dynamically ejected from the system, and could thus contribute substantially to the amount of observed r-process material in the galaxy. We calculate the gravitational radiation waveforms and luminosity in the quadrupole approximation.