No Arabic abstract
(Abridged) In this paper we present a compilation of results from our most advanced neutron star merger simulations, including a description of the employed numerical procedures and a more complete overview over a large number of computed models. The three-dimensional hydrodynamic simulations were done with a code based on the Piecewise Parabolic Method with up to five levels of nested Cartesian grids. The simulations are basically Newtonian, but gravitational-wave emission and the corresponding back-reaction are taken into account. The use of a physical nuclear equation of state allows us to follow the thermodynamic history of the stellar medium and to compute the energy and lepton number loss due to the emission of neutrinos. The computed models differ concerning the neutron star masses and mass ratios, the neutron star spins, the numerical resolution expressed by the cell size of the finest grid and the number of grid levels, and the calculation of the temperature from the solution of the entropy equation instead of the energy equation. Our simulations show that the details of the gravitational-wave emission are still sensitive to the numerical resolution, even in our highest-quality calculations. The amount of mass which can be ejected from neutron star mergers depends strongly on the angular momentum of the system. Our results do not support the initial conditions of temperature and proton-to-nucleon ratio assumed in recent work for producing a solar r-process pattern for nuclei around and above the A approx 130 peak. The improved models confirm our previous conclusion that gamma-ray bursts are not powered by neutrino emission during the dynamical phase of the merging of two neutron stars.
Three-dimensional hydrodynamical, Newtonian calculations of the coalescence of equal-mass binary neutron stars are performed, including a physical high-density equation of state and a treatment of the neutrino emission of the heated matter. The total neutrino luminosity climbs to a maximum value of 1--$1.5cdot 10^{53}$~erg/s of which 90--95% originate from the toroidal gas cloud surrounding the very dense core formed after the merging. When the neutrino luminosities are highest, $ ubar u$-annihilation deposits about 0.2--0.3% of the emitted neutrino energy in the immediate neighborhood of the merger, and the maximum integral energy deposition rate is 3--$4cdot 10^{50}$~erg/s. Since the $3,M_{odot}$ core of the merged object will most likely collapse into a black hole within milliseconds, the energy that can be pumped into a pair-photon fireball is insufficient by a factor of about 1000 to explain $gamma$-ray bursts at cosmological distances with an energy of the order of $10^{51}/(4pi)$~erg/steradian. Analytical estimates show that the additional energy provided by the annihilation of $ ubar u$ pairs emitted from a possible accretion torus of $sim 0.1,M_{odot}$ around the central black hole is still more than a factor of 10 too small, unless focussing of the fireball into a jet-like expansion plays an important role. About $10^{-4}$--$10^{-3}$~$M_odot$ of material lost during the neutron star merging and swept out from the system in a neutrino-driven wind might be a site for nucleosythesis. Aspects of a possible r-processing in these ejecta are discussed.
Precision mass spectrometry of neutron-rich nuclei is of great relevance for astrophysics. Masses of exotic nuclides impose constraints on models for the nuclear interaction and thus affect the description of the equation of state of nuclear matter, which can be extended to describe neutron-star matter. With knowledge of the masses of nuclides near shell closures, one can also derive the neutron-star crustal composition. The Penning-trap mass spectrometer ISOLTRAP at CERN-ISOLDE has recently achieved a breakthrough measuring the mass of 82Zn, which allowed constraining neutron-star crust composition to deeper layers (Wolf et al., PRL 110, 2013). We perform a more detailed study on the sequence of nuclei in the outer crust of neutron stars with input from different nuclear models to illustrate the sensitivity to masses and the robustness of neutron-star models. The dominant role of the N=50 and N=82 closed neutron shells for the crustal composition is confirmed.
Recently exploratory studies were performed on the possibility of constraining the neutron star equation of state (EOS) using signals from coalescing binary neutron stars, or neutron star-black hole systems, as they will be seen in upcoming advanced gravitational wave detectors such as Advanced LIGO and Advanced Virgo. In particular, it was estimated to what extent the combined information from multiple detections would enable one to distinguish between different equations of state through hypothesis ranking or parameter estimation. Under the assumption of zero neutron star spins both in signals and in template waveforms and considering tidal effects to 1 post-Newtonian (1PN) order, it was found that O(20) sources would suffice to distinguish between a hard, moderate, and soft equation of state. Here we revisit these results, this time including neutron star tidal effects to the highest order currently known, termination of gravitational waveforms at the contact frequency, neutron star spins, and the resulting quadrupole-monopole interaction. We also take the masses of neutron stars in simulated sources to be distributed according to a relatively strongly peaked Gaussian, as hinted at by observations, but without assuming that the data analyst will necessarily have accurate knowledge of this distribution for use as a mass prior. We find that especially the effect of the latter is dramatic, necessitating many more detections to distinguish between different EOS and causing systematic biases in parameter estimation, on top of biases due to imperfect understanding of the signal model pointed out in earlier work. This would get mitigated if reliable prior information about the mass distribution could be folded into the analyses.
We discuss new limits on masses and radii of compact stars and we conclude that they can be interpreted as an indication of the existence of two classes of stars: normal compact stars and ultra-compact stars. We estimate the critical mass at which the first configuration collapses into the second.
We discuss three topics related to the neutron star (NS) mass spectrum. At first we discuss the possibility to form low-mass ($ M stackrel{<}{sim} 1 M_{odot}$) and suggest this is possible only due to fragmentation of rapidly rotating proto-NSs. Such low-mass NSs should have very high spatial velocities which could allow identification. A critical assessment of this scenario is given. Secondly, we discuss mass growth due to accretion for NSs in close binary systems. With the help of numerical population synthesis calculations we derive the mass spectrum of massive ($M > 1.8 M_{odot}$) NSs. Finally, we discuss the role of the mass spectrum in population studies of young cooling NSs. We formulate a kind of {it mass constraint} which can be helpful, in our opinion, in discussing different competive models of the thermal evolution of NSs.