# Simulations and some results

## Simulations and some results

NS-NS and NS-BH mergers have been suggested as sources of $`\gamma `$-ray bursts at cosmological distances (e.g., ) and of r-process nuclei . These merger events are estimated to occur at rates of the order $`10^5`$ per year per galaxy (e.g., ).

By means of three-dimensional hydrodynamic simulations with a Eulerian PPM code we have investigated the dynamical phases of the coalescence . The simulations were performed with four nested cartesian grids which allow for both a good resolution near the central black hole and a large computational volume. In a post-processing step, we evaluated our models for the energy deposition by $`\nu \overline{\nu }`$ annihilation around the merging objects.

| Table 1. Results of merger models. Energies are given in units of 1foe = $`10^{51}\mathrm{erg}`$. | | | | | | | | | | | | |
| --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- | --- |
| model | masses | spin | $`\mathrm{\Delta }M_\mathrm{d}`$ | $`t_{\mathrm{sim}}`$ | $`L_\nu `$ | $`a_{\mathrm{BH}}`$ | $`q_\nu `$ | $`q_{\nu \overline{\nu }}`$ | $`\dot{E}_{\nu \overline{\nu }}`$ | $`t_{\mathrm{acc}}`$ | $`E_\nu `$ | $`E_{\nu \overline{\nu }}`$ |
|  | ($`M_{}`$) |  | ($`\frac{M_{}}{100}`$) | $`\left(\mathrm{ms}\right)`$ | ($`100\frac{\mathrm{foe}}{\mathrm{s}}`$) |  | % | % | ($`\frac{\mathrm{foe}}{\mathrm{s}}`$) | $`\left(\mathrm{ms}\right)`$ | $`\left(\mathrm{foe}\right)`$ | $`\left(\mathrm{foe}\right)`$ |
| NS-NS | 1.2–1.2 | solid | 2.0 | 10 | 1.5 |  |  | 1.0 | 1.4 |  |  |  |
| NS-NS | 1.2–1.8 | solid | 3.8 | 13 | 2.0 |  |  | 1.3 | 2.6 |  |  |  |
| NS-NS | 1.6–1.6 | anti | 0.01 | 10 | 4.0 |  |  | 1.5 | 6.0 |  |  |  |
| NS-NS | 1.6–1.6 | none | 0.23 | 10 | 5.0 |  |  | 1.8 | 9.1 |  |  |  |
| NS-NS | 1.6–1.6 | solid | 2.4 | 10 | 3.3 |  |  | 2.1 | 7.0 |  |  |  |
| BH-AD | 3.0–0.26 | solid |  | 15 | 1.2 | 0.42 | 1.3 | 0.4 | 0.49 | 68 | 8 | 0.03 |
| BH-NS | 2.5–1.6 | solid | 0.2 | 10 | 8.0 | 0.76 | 3.0 | 2.8 | 22.4 | 40 | 32 | 0.9 |
| BH-NS | 5.0–1.6 | solid | 5.6 | 15 | 5.0 | 0.33 | 2.8 | 1.8 | 9.0 | 56 | 28 | 0.5 |
| BH-NS | 10.0–1.6 | solid | 10.0 | 10 | 2.5 | 0.14 | 2.8 | 0.9 | 2.3 | 82 | 21 | 0.2 |

Results for some simulations are listed in Table 1. The NS and accretion disk (AD) masses are baryonic masses, the spin parameter indicates whether no spin was assumed for the NS’s (“none”), or the NS spin was taken parallel (“solid”) or anti-parallel (“anti”) to the orbital angular momentum. $`\mathrm{\Delta }M_\mathrm{d}`$ is the dynamically ejected mass, $`t_{\mathrm{sim}}`$ the time at the end of the simulation, $`L_\nu `$ the typical total neutrino luminosity, $`a_{\mathrm{BH}}=Jc/(GM_{\mathrm{BH}}^2)`$ the relativistic rotation parameter of the black hole, $`q_\nu =L_\nu /(\dot{M}c^2)`$ the conversion efficiency of rest-mass energy of the gas swallowed by the black hole into $`\nu `$ energy, $`q_{\nu \overline{\nu }}=E_{\nu \overline{\nu }}/E_\nu `$ the conversion efficiency of $`\nu `$ energy into $`e^\pm `$ pairs, $`\dot{E}_{\nu \overline{\nu }}`$ the integral rate of energy deposition by $`\nu \overline{\nu }`$ annihilation, $`t_{\mathrm{acc}}`$ the estimated duration of the accretion by the BH, $`E_\nu `$ the energy released in neutrinos during $`t_{\mathrm{acc}}`$, and $`E_{\nu \overline{\nu }}`$ the corresponding energy converted into $`e^\pm `$ pairs.

We find that during the dynamical phase of the merging between about $`10^4M_{}`$ and $`10^2M_{}`$ of cool gas can be ejected (Fig. 1). The exact number depends very sensitively on the assumed neutron star spin(s) and on the masses of the merging binary components (Table 1). This is of potential interest for heavy-element nucleosynthesis .

Figure 1 sketches the result that during the dynamical interaction of the binary components (1st phase) low-entropy (dense and cold) matter is ejected from the ends of long spiral arms which are formed in the orbital plane. During the merging the matter is heated by shocks and viscous dissipation, and intense neutrino fluxes are emitted before the gas is accreted into the central black hole. Energy transfer by these neutrinos drives a high-entropy gas outflow (“neutrino wind”) off the surface of the accretion disk. At the same time, energy deposition by $`\nu \overline{\nu }`$ annihilation in the essentially baryon-free funnel around the rotation axis powers relativistically expanding $`e^\pm \gamma `$ jets which can give rise to $`\gamma `$-ray bursts.

Our simulations confirm suggestions that these events could be the sources of a subclass of short gamma-ray bursts at cosmological distances (Fig. 2). BH-NS mergers yield between 10 and 100 times more energy than NS-NS mergers (Table 1). The annihilation of neutrinos and antineutrinos (Fig. 3) which are emitted from the hot matter that is accreted into the (forming) black hole, yields total energies $`E_{\nu \overline{\nu }}`$ up to $`10^{51}\mathrm{erg}`$ in our models. It is therefore enough efficient to explain observable burst luminosities $`L_\gamma E_{\nu \overline{\nu }}/(f_\mathrm{\Omega }t_\gamma )`$ up to several $`10^{53}\mathrm{erg}\mathrm{s}^1`$ for burst durations $`t_\gamma 0.1`$–1 s, if the $`\gamma `$ emission is beamed in two moderately focussed jets into a fraction $`f_\mathrm{\Omega }=2\delta \mathrm{\Omega }/(4\pi )1/100`$$`1/10`$ of the sky. Such a modest amount of beaming (jet opening half-angles between about ten and several ten degrees) has to be expected because the dense gas in the equatorial plane of the accreting black hole prevents relativistic expansion of the pair-plasma jets outside a low-density funnel along the system axis. More energy for the burst can be available if Kerr effects of the rotating black hole are taken into account.

The evolution tracks of compact binaries shown in Fig. 2 depend on the NS and BH masses, on the total angular momentum in the system, and on the stiffness of the nuclear equation of state. Energetic relativistic pair-plasma jets can be produced when a hot accretion disk forms around the BH, because high neutrino fluxes are emitted over a time much longer than the dynamical timescale of the system. This allows for high efficiencies $`q_\nu `$ and $`q_{\nu \overline{\nu }}`$ and a large total energy conversion by $`\nu \overline{\nu }e^+e^{}`$. Gamma-ray bursts from these scenarios are accompanied by the ejection of nucleosynthesis products in considerable amounts, the emission of large numbers of neutrinos with MeV energies, and the production of strong gravitational-wave signals.

Acknowledgements.  MR is grateful for support by a PPARC Advanced Fellowship, HTJ acknowledges support by the DFG on grant “SFB 375 für Astro-Teilchenphysik”.