No Arabic abstract
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.
By means of three-dimensional hydrodynamic simulations with a Eulerian PPM code we investigate the formation and the properties of the accretion torus around the stellar mass black hole which originates from the merging of two neutron stars. The simulations are performed with four nested cartesian grids which allow for both a good resolution near the central black hole and a large computational volume. They include the use of a physical equation of state as well as the neutrino emission from the hot matter of the torus. The gravity of the black hole is described with a Newtonian and alternatively with a Paczynski-Wiita potential. In a post-processing step, we evaluate our models for the energy deposition by nu-nubar annihilation around the accretion torus. Our models show that nu-nubar annihilation can yield the energy to account for weak, short gamma-ray bursts, if moderate beaming is involved. In fact, the barrier of the dense baryonic gas of the torus suggests that the low-density pair-photon-plasma is beamed as axial jets into a fraction 2 delta Omega/ (4 pi) between 1/100 and 1/10 of the sky, corresponding to opening half-angles of roughly ten to several tens of degrees. Thus gamma-burst energies of 10^{50}--10^{51} erg seem within the reach of our models (if the source is interpreted as radiating isotropically), corresponding to luminosities around 10^{51} erg/s for typical burst durations of 0.1--1 s. Gravitational capture of radiation by the black hole, redshift and ray bending do not reduce the jet energy significantly. Effects associated with the Kerr character of the rapidly rotating black hole, however, could increase the gamma-burst energy considerably, and effects due to magnetic fields might even be required to get the energies for long complex gamma-ray bursts.
The central engine of Gamma Ray Bursts may live much longer than the duration of the prompt emission. Some evidence of it comes from the presence of strong precursors, post-cursors, and X-ray flares in a sizable fraction of bursts. Additional evidence comes from the fact that often the X-ray and the optical afterglow light curves do not track one another, suggesting that they are two different emission components. The typical steep-flat-steep behavior of the X-ray light curve can be explained if the same central engine responsible for the main prompt emission continues to be active for a long time, but with a decreasing power. The early X-ray afterglow emission is then the extension of the prompt emission, originating at approximately the same location, and is not due to forward shocks. If the bulk Lorentz factor Gamma is decreasing in time, the break ending the shallow phase can be explained, since at early times Gamma is large, and we see only a fraction of the emitting area. Later, when Gamma decreases, we see an increasing fraction of the emitting surface up to the time when Gamma ~ 1/theta_j. This time ends the shallow phase of the X-ray light curve. The origin of the late prompt emission can be the accretion of the fall-back material, with an accretion rate dot M proportional to t^(-5/3). The combination of this late prompt emission with the flux produced by the standard forward shock can explain the great diversity of the optical and the X-ray light curves.
Short gamma-ray bursts may originate from the merger of double neutron stars (NS) or that of a black hole (BH) and an NS. We propose that the bright X-ray flare related to the central engine reactivity may hint a BH-NS merger, since such a merger can provide more fall-back materials and therefore a more massive accretion disk than the NS-NS merger. Based on the observed 49 short bursts with Swift/X-ray Telescope follow-up observations, we find that three bursts have bright X-ray flares, among which three flares from two bursts are probably related to the central engine reactivity. We argue that these two bursts may originate from the BH-NS merger rather than the NS-NS merger. Our suggested link between the central engine-powered bright X-ray flare and the BH-NS merger event can be checked by the future gravitational wave detections from advanced LIGO and Virgo.
Observations of gravitational waves and their electromagnetic counterparts may soon uncover the existence of coalescing compact binary systems formed by a stellar-mass black hole and a neutron star. These mergers result in a remnant black hole, possibly surrounded by an accretion disk. The mass and spin of the remnant black hole depend on the properties of the coalescing binary. We construct a map from the binary components to the remnant black hole using a sample of numerical-relativity simulations of different mass ratios $q$, (anti-)aligned dimensionless spins of the black hole $a_{rm BH}$, and several neutron star equations of state. Given the binary total mass, the mass and spin of the remnant black hole can therefore be determined from the three parameters $(q,a_{rm BH},Lambda)$, where $Lambda$ is the tidal deformability of the neutron star. Our models also incorporate the binary black hole and test-mass limit cases and we discuss a simple extension for generic black hole spins. We combine the remnant characterization with recent population synthesis simulations for various metallicities of the progenitor stars that generated the binary system. We predict that black-hole-neutron-star mergers produce a population of remnant black holes with masses distributed around $7M_odot$ and $9M_odot$. For isotropic spin distributions, nonmassive accretion disks are favoured: no bright electromagnetic counterparts are expected in such mergers.
The first locations of short gamma-ray bursts (GRBs) in elliptical galaxies suggest they are produced by the mergers of double neutron star (DNS) binaries in old stellar populations. Globular clusters, where the extreme densities of very old stars in cluster cores create and exchange compact binaries efficiently, are a natural environment to produce merging NSs. They also allow some short GRBs to be offset from their host galaxies, as opposed to DNS systems formed from massive binary stars which appear to remain in galactic disks. Starting with a simple scaling from the first DNS observed in a galactic globular, which will produce a short GRB in ~300My, we present numerical simulations which show that ~10-30% of short GRBs may be produced in globular clusters vs. the much more numerous DNS mergers and short GRBs predicted for galactic disks. Reconciling the rates suggests the disk short GRBs are more beamed, perhaps by both the increased merger angular momentum from the DNS spin-orbit alignment (random for the DNS systems in globulars) and a larger magnetic field on the secondary NS.