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.

We perform three-dimensional relativistic hydrodynamical simulations of the coalescence of strange stars (SSs) and explore the possibility to decide on the strange matter hypothesis by means of gravitational-wave (GW) measurements. Selfbinding of strange quark matter (SQM) and the generally more compact stars yield features that clearly distinguish SS from neutron star (NS) mergers, e.g. hampering tidal disruption during the plunge of quark stars. Furthermore, instead of forming dilute halos around the remnant as in the case of NS mergers, the coalescence of SSs results in a differentially rotating hypermassive object with a sharp surface layer surrounded by a geometrically thin, clumpy high-density SQM disk. We also investigate the importance of including non-zero temperature equations of state (EoSs) in NS and SS merger simulations. In both cases we find a crucial sensitivity of the dynamics and outcome of the coalescence to thermal effects, which, e.g., determine the outer remnant structure and the delay time of the dense remnant core to black hole collapse. For comparing and classifying the GW signals, we use a number of characteristic quantities like the maximum frequency during inspiral or the dominant frequency of oscillations of the postmerger remnant. In general, these frequencies are higher for SS mergers. If not, additional features of the GW luminosity spectrum may help to discriminate coalescence events of the different types. Future GW measurements may thus help to decide on the existence of SQM stars. (abridged)

We determine the Galactic production rate of strangelets as a canonical input to calculations of the measurable cosmic ray flux of strangelets by performing simulations of strange star mergers and combining the results with recent estimates of stellar binary populations. We find that the flux depends sensitively on the bag constant of the MIT bag model of QCD and disappears for high values of the bag constant and thus more compact strange stars. In the latter case strange stars could coexist with ordinary neutron stars as they are not converted by the capture of cosmic ray strangelets. An unambiguous detection of an ordinary neutron star would then not rule out the strange matter hypothesis.

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.

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)

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.

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.

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.

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.

With the exception of the Big Bang, responsible for 1,2H, 3,4He, and 7Li, stars act as sources for the composition of the interstellar medium. Cosmic rays are related to the latter and very probably due to acceleration of the mixed interstellar medium by shock waves from supernova remnants. Thus, the understanding of the abundance evolution in the interstellar medium and especially the enrichment of heavy elements, as a function of space and time, is essential. It reflects the history of star formation and the lifetimes of the diverse contributing stellar objects. Therefore, the understanding of the endpoints of stellar evolution is essential as well. These are mainly planetary nebulae and type II/Ib/Ic supernovae as evolutionary endpoints of single stars, but also events in binary systems can contribute, like e.g. supernovae of type Ia, novae and possibly X-ray bursts and neutron star or neutron star - black hole mergers. Despite many efforts, a full and self-consistent understanding of supernovae (the main contributors to nucleosynthesis in galaxies) is not existing, yet. Their fingerprints, however, seen either in spectra, lightcurves, radioactivities/decay gamma-rays or in galactic evolution, can help to constrain the composition of their ejecta and related model uncertainties.

We present a new 3D SPH code which solves the general relativistic field + hydrodynamics equations in the conformally flat approximation. Several test cases are considered to test different aspects of the code. We finally apply then the code to the coalescence of a neutron star binary system. The neutron stars are modeled by a polytropic equation of state (EoS) with adiabatic indices $\Gamma=2.0$, $\Gamma=2.6$ and $\Gamma=3.0$. We calculate the gravitational wave signals, luminosities and frequency spectra by employing the quadrupole approximation for emission and back reaction in the slow motion limit. In addition, we consider the amount of ejected mass.

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.