Black Hole - Neutron Star Mergers as Central Engines of Gamma-Ray Bursts

Hydrodynamic simulations of the merger of stellar mass black hole - neutron star binaries (BH/NS) are compared with mergers of binary neutron stars (NS/NS). The simulations are Newtonian, but take into account the emission and backreaction of gravitational waves. The use of a physical nuclear equation of state allows us to include the effects of neutrino emission. For low neutron star to black hole mass ratios the neutron star transfers mass to the black hole during a few cycles of orbital decay and subsequent widening before finally being disrupted, whereas for ratios near unity the neutron star is already distroyed during its first approach. A gas mass between about 0.3 and about 0.7 solar masses is left in an accretion torus around the black hole and radiates neutrinos at a luminosity of several 10^{53} erg/s during an estimated accretion time scale of about 0.1 s. The emitted neutrinos and antineutrinos annihilate into electron-positron pairs with efficiencies of 1-3% percent and rates of up to 2*10^{52} erg/s, thus depositing an energy of up to 10^{51} erg above the poles of the black hole in a region which contains less than 10^{-5} solar masses of baryonic matter. This could allow for relativistic expansion with Lorentz factors around 100 and is sufficient to explain apparent burst luminosities of up to several 10^{53} erg/s for burst durations of approximately 0.1-1 s, if the gamma emission is collimated in two moderately focussed jets in a fraction of about 1/100-1/10 of the sky.Comment: 8 pages, LaTex, 4 postscript figures, 2 tables. ApJ Letters, accepted; revised and shortened version, Fig. 2 change

A Solution to the Protostellar Accretion Problem

Accretion rates of order 10^-8 M_\odot/yr are observed in young protostars of approximately a solar mass with evidence of circumstellar disks. The accretion rate is significantly lower for protostars of smaller mass, approximately proportional to the second power of the stellar mass, \dot{M}_accr\propto M^2. The traditional view is that the observed accretion is the consequence of the angular momentum transport in isolated protostellar disks, controlled by disk turbulence or self--gravity. However, these processes are not well understood and the observed protostellar accretion, a fundamental aspect of star formation, remains an unsolved problem. In this letter we propose the protostellar accretion rate is controlled by accretion from the large scale gas distribution in the parent cloud, not by the isolated disk evolution. Describing this process as Bondi--Hoyle accretion, we obtain accretion rates comparable to the observed ones. We also reproduce the observed dependence of the accretion rate on the protostellar mass. These results are based on realistic values of the ambient gas density and velocity, as inferred from numerical simulations of star formation in self--gravitating turbulent clouds.Comment: 4 pages, 2 figures, ApJ Letters, in pres

Numerical Models of Binary Neutron Star System Mergers. I.: Numerical Methods and Equilibrium Data for Newtonian Models

The numerical modeling of binary neutron star mergers has become a subject of much interest in recent years. While a full and accurate model of this phenomenon would require the evolution of the equations of relativistic hydrodynamics along with the Einstein field equations, a qualitative study of the early stages on inspiral can be accomplished by either Newtonian or post-Newtonian models, which are more tractable. In this paper we offer a comparison of results from both rotating and non-rotating (inertial) frame Newtonian calculations. We find that the rotating frame calculations offer significantly improved accuracy as compared with the inertial frame models. Furthermore, we show that inertial frame models exhibit significant and erroneous angular momentum loss during the simulations that leads to an unphysical inspiral of the two neutron stars. We also examine the dependence of the models on initial conditions by considering initial configurations that consist of spherical neutron stars as well as stars that are in equilibrium and which are tidally distorted. We compare our models those of Rasio & Shapiro (1992,1994a) and New & Tohline (1997). Finally, we investigate the use of the isolated star approximation for the construction of initial data.Comment: 32 pages, 19 gif figures, manuscript with postscript figures available at http://www.astro.sunysb.edu/dswesty/docs/nspap1.p

Neutrino and anti-neutrino transport in accretion disks

We numerically solve the one dimensional Boltzmann equation of the neutrino and anti-neutrino transport in accretion disks and obtain the fully energy dependent and direction dependent neutrino and anti-neutrino emitting spectra, under condition that the distribution of the mass density,temperature and chemical components are given. Then, we apply the resulting neutrino and anti-neutrino emitting spectra to calculate the corresponding annihilation rate of neutrino pairs above the neutrino dominated accretion disk and find that the released energy resulting from the annihilation of neutrino pairs can not provide sufficient energy for the most energetic short gamma ray bursts whose isotropic luminosity can be as high as $10^{52}$ ergs/s unless the high temperature zone where the temperature is beyond 10 MeV can stretch over 200 km in the disk. We also compare the resulting luminosity of neutrinos and anti-neutrinos with the results from the two commonly used approximate treatment of the neutrino and anti-neutrino luminosity: the Fermi-Dirac black body limit and a simplified model of neutrino transport, i.e., the gray body model, and find that both of them overestimate the neutrino/anti-neutrino luminosity and their annihilation rate greatly. Additionally, as did in Sawyer (2003), we also check the validity of the two stream approximation, and find that it is a good approximation to high accuracy.Comment: Phys. Rev. D in press, 15 preprint papers, 5 figure