Short Gamma-Ray Bursts from Binary Neutron Star Mergers

We present the results from new relativistic hydrodynamic simulations of binary neutron star mergers using realistic non-zero temperature equations of state. We vary several unknown parameters in the system such as the neutron star (NS) masses, their spins and the nuclear equation of state. The results are then investigated with special focus on the post-merger torus-remnant system. Observational implications on the Gamma-ray burst (GRB) energetics are discussed and compared with recent observations

Quark Matter in Neutron Star Mergers

Binary neutron star mergers are expected to be one of the most promising source of gravitational waves (GW) for the network of laser interferometric and bar detectors becoming operational in the next few years. The merger wave signal is expected to be sensitive to the interior structure of the neutron star (NS). The structure of high density phases of matter is under current experimental investigation in heavy-ion collisions. We investigate the dependence of the merger process and its GW signal on the presence of quarks in these phases by performing numerical simulations, where the smoothed particle hydrodynamics (SPH) method and the conformally flat approximation for the 3-geometry in general relativistic gravity are implemented.Comment: 4 Pages, 4 Figures, Proc. Nuclei in the Cosmos 7, 200

Testing Approximations of Thermal Effects in Neutron Star Merger Simulations

We perform three-dimensional relativistic hydrodynamical calculations of neutron star mergers to assess the reliability of an approximate treatment of thermal effects in such simulations by combining an ideal-gas component with zero-temperature, micro-physical equations of state. To this end we compare the results of simulations that make this approximation to the outcome of models with a consistent treatment of thermal effects in the equation of state. In particular we focus on the implications for observable consequences of merger events like the gravitational-wave signal. It is found that the characteristic gravitational-wave oscillation frequencies of the post-merger remnant differ by about 50 to 250 Hz (corresponding to frequency shifts of 2 to 8 per cent) depending on the equation of state and the choice of the characteristic index of the ideal-gas component. In addition, the delay time to black hole collapse of the merger remnant as well as the amount of matter remaining outside the black hole after its formation are sensitive to the description of thermal effects.Comment: 10 pages, 6 figures, 9 eps files; revised with minor additions due to referee comments; accepted by Phys.Rev.

Relativistic neutron star merger simulations with non-zero temperature equations of state I. Variation of binary parameters and equation of state

An extended set of binary neutron star (NS) merger simulations is performed with an approximative conformally flat treatment of general relativity to systematically investigate the influence of the nuclear equation of state (EoS), the neutron star masses, and the NS spin states prior to merging. We employ the two non-zero temperature EoSs of Shen et al. (1998a,b) and Lattimer & Swesty (1991). In addition, we use the cold EoS of Akmal et al. (1998) with a simple ideal-gas-like extension according to Shibata & Taniguchi (2006), and an ideal-gas EoS with parameters fitted to the supernuclear part of the Shen-EoS. We estimate the mass sitting in a dilute high-angular momentum ``torus'' around the future black hole (BH). The dynamics and outcome of the models is found to depend strongly on the EoS and on the binary parameters. Larger torus masses are found for asymmetric systems (up to ~0.3 M_sun for a mass ratio of 0.55), for large initial NSs, and for a NS spin state which corresponds to a larger total angular momentum. We find that the postmerger remnant collapses either immediately or after a short time when employing the soft EoS of Lattimer& Swesty, whereas no sign of post-merging collapse is found within tens of dynamical timescales for all other EoSs used. The typical temperatures in the torus are found to be about 3-10 MeV depending on the strength of the shear motion at the collision interface between the NSs and thus depending on the initial NS spins. About 10^{-3}-10^{-2} M_sun of NS matter become gravitationally unbound during or right after the merging process. This matter consists of a hot/high-entropy component from the collision interface and (only in case of asymmetric systems) of a cool/low-entropy component from the spiral arm tips. (abridged)Comment: 20 pages, 15 figures, accepted for publication in A&A, included changes based on referee comment

Dynamical non-axisymmetric instabilities in rotating relativistic stars

We present new results on dynamical instabilities in rapidly rotating neutron-stars. In particular, using numerical simulations in full General Relativity, we analyse the effects that the stellar compactness has on the threshold for the onset of the dynamical bar-mode instability, as well as on the appearance of other dynamical instabilities. By using an extrapolation technique developed and tested in our previous study [1], we explicitly determine the threshold for a wide range of compactnesses using four sequences of models of constant baryonic mass comprising a total of 59 stellar models. Our calculation of the threshold is in good agreement with the Newtonian prediction and improves the previous post-Newtonian estimates. In addition, we find that for stars with sufficiently large mass and compactness, the m=3 deformation is the fastest growing one. For all of the models considered, the non-axisymmetric instability is suppressed on a dynamical timescale with an m=1 deformation dominating the final stages of the instability. These results, together with those presented in [1], suggest that an m=1 deformation represents a general and late-time feature of non-axisymmetric dynamical instabilities both in full General Relativity and in Newtonian gravity.Comment: To appear on CQG, NFNR special issue. 16 pages, 5 color figures, movies from http://www.fis.unipr.it/numrel/BarMode/ResearchBarMode.htm

The influence of quark matter at high densities on binary neutron star mergers

We consider the influence of potential quark matter existing at high densities in neutron star (NS) interiors on gravitational waves (GWs) emitted in a binary NS merger event. Two types of equations of state (EoSs) at zero temperature are used - one describing pure nuclear matter and the other nuclear matter with a phase transition to quark matter at very high densities. Binary equilibrium sequences close to the innermost stable circular orbit (ISCO) are calculated to determine the GW frequencies just before the merger. It is found that the effects of the EoSs begin to play a role when gravitational masses are larger than M∞≃ 1.5 M⊙. The difference in the GW frequency at the ISCO increases by up to ≃10 per cent for the maximum mass permitted by the EoSs. We then perform three-dimensional hydrodynamic simulations for each EoS while varying the initial mass and determine the characteristic GW frequencies of the merger remnant. The implications of the presence of quark matter show up mainly in the collapse behaviour of the merger remnant. If the collapse does not take place immediately after the merger, we find a phase difference between the two EoSs in the post-merger GW signal. We also compare the GW frequencies emitted by the remnant of the merger to values obtained from simulations using a polytropic EoS and find an imprint of the non-constant adiabatic index of our EoSs. All calculations are based on the conformally flat approximation to general relativity and the GW signal from the merger simulation is extracted up to quadrupole orde

Torus Formation in Neutron Star Mergers and Well-Localized Short Gamma-Ray Bursts

Merging neutron stars (NSs) are hot candidates for the still enigmatic sources of short gamma-ray bursts (GRBs). If the central engines of the huge energy release are accreting relic black holes (BHs) of such mergers, it is important to understand how the properties of the BH-torus systems, in particular disc masses and mass and rotation rate of the compact remnant, are linked to the characterizing parameters of the NS binaries. For this purpose we present relativistic smoothed particle hydrodynamics simulations with conformally flat approximation of the Einstein field equations and a physical, non-zero temperature equation of state. Thick disc formation is highlighted as a dynamical process caused by angular momentum transfer through tidal torques during the merging process of asymmetric systems or in the rapidly spinning triaxial post-merger object. Our simulations support the possibility that the first well-localized short and hard GRBs 050509b, 050709, 050724, 050813 have originated from NS merger events and are powered by neutrino-antineutrino annihilation around a relic BH-torus system. Using model parameters based on this assumption, we show that the measured GRB energies and durations lead to estimates for the accreted masses and BH mass accretion rates which are compatible with theoretical expectations. In particular, the low energy output and short duration of GRB 050509b set a very strict upper limit of less than 100 ms for the time interval after the merging until the merger remnant has collapsed to a BH, leaving an accretion torus with a small mass of only about 0.01 solar masses. This favors a (nearly) symmetric NS+NS binary with a typical mass as progenitor system.Comment: 12 pages, 9 figures, high-resolution color figures available on request; accepted by MNRA

Gravitational waves from relativistic neutron star mergers with nonzero-temperature equations of state

We analyze the gravitational wave (GW) emission from our recently published set of relativistic neutron star (NS) merger simulations and determine characteristic signal features that allow one to link GW measurements to the properties of the merging binary stars. We find that the distinct peak in the GW energy spectrum that is associated with the formation of a hypermassive merger remnant has a frequency that depends strongly on the properties of the nuclear equation of state (EoS) and on the total mass of the binary system, whereas the mass ratio and the NS spins have a weak influence. If the total mass can be determined from the inspiral chirp signal, the peak frequency of the postmerger signal is a sensitive indicator of the EoS.Comment: 5 pages, 4 figures, revised version accepted for publication in PR

The Influence of Quark Matter at High Densities on Binary Neutron Star Mergers

We consider the influence of potential quark matter existing at high densities in neutron star interiors on gravitational waves (GW) emitted in a binary neutron star merger event. Two types of equations of state (EoS) at zero temperatures are used, one describing pure nuclear matter, the other nuclear matter with a phase transition to quark matter at very high densities. Binary equilibrium sequences close to the innermost stable circular orbit (ISCO) are calculated to determine the GW frequencies just before merger. It is found that EoS effects begin to play a role for gravitational masses larger than $M_\infty\simeq1.5M_\odot$. The difference in the gravitational wave frequency at the ISCO grows to up to $\simeq 10$ for the maximal allowed mass given by the EoSs used. Then, we perform 3D hydrodynamic simulations for each EoS varying the initial mass and determine the characteristic GW frequencies of the merger remnants. The implications of quark matter show up mainly in a different collapse behaviour of the merger remnant. If the collapse does not take place immediately after merger, we find a phase difference between two EoS's in the post-merger GW signal. We also compare the GW frequencies emitted by the merger remnant to values from simulations using a polytropic EoS and find an imprint of the non-constant adiabatic index of our EoSs. All calculations are based on the conformally flat (CF) approximation to general relativity and the GW signal from the merger simulation is extracted up to quadrupole order.Comment: 13 pages, 8 figure

Post Newtonian SPH

We introduce an adaptation of the well known Tree+SPH numerical scheme to Post Newtonian (PN) hydrodynamics and gravity. Our code solves the (0+1+2.5)PN equations. These equations include Newtonian hydrodynamics and gravity (0PN), the first order relativistic corrections to those (1PN) and the lowest order gravitational radiation terms (2.5PN). We test various aspects of our code using analytically solvable test problems. We then proceed to study the 1PN effects on binary neutron star coalescence by comparing calculations with and without the 1PN terms. We find that the effect of the 1PN terms is rather small. The largest effect arises with a stiff equation of state for which the maximum rest mass density increases. This could induce black hole formation. The gravitational wave luminosity is also affected.Comment: 28 pages, 13 figures, revised version published in Ap