We perform binary neutron star merger simulations using a newly derived set of finite-temperature equations of state in the Brueckner-Hartree-Fock approach. We point out the important and opposite roles of finite temperature and rotation for stellar stability and systematically investigate the gravitational-wave properties, matter distribution, and ejecta properties in the postmerger phase for the different cases. The validity of several universal relations is also examined and the most suitable EOSs are identified.
We investigate the use of hybrid equations of state in binary neutron-star simulations in full general relativity, where thermal effects are included in an approximate way through the adiabatic index $Gamma_{rm{th}}$. We employ a newly developed finite-temperature equation of state derived in the Brueckner-Hartree-Fock approach and carry out comparisons with the corresponding hybr
We present results from three-dimensional general relativistic simulations of binary neutron star coalescences and mergers using public codes. We considered equal mass models where the baryon mass of the two Neutron Stars (NS) is $1.4M_{odot}$, described by four different equations of state (EOS) for the cold nuclear matter (APR4, SLy, H4, and MS1; all parametrized as piecewise polytropes). We started the simulations from four different initial interbinary distances ($40, 44.3, 50$, and $60$ km), including up to the last 16 orbits before merger. That allows to show the effects on the gravitational wave phase evolution, radiated energy and angular momentum due to: the use of different EOSs, the orbital eccentricity present in the initial data and the initial separation (in the simulation) between the two stars. Our results show that eccentricity has a major role in the discrepancy between numerical and analytical waveforms until the very last few orbits, where tidal effects and missing high-order post-Newtonian coefficients also play a significant role. We test different methods for extrapolating the gravitational wave signal extracted at finite radii to null infinity. We show that an effective procedure for integrating the Newman-Penrose $psi_4$ signal to obtain the gravitational wave strain $h$ is to apply a simple high-pass digital filter to $h$ after a time domain integration, where only the two physical motivated integration constants are introduced. That should be preferred to the more common procedures of introducing additional integration constants, integrating in the frequency domain or filtering $psi_4$ before integration.
We study the properties of hot beta-stable nuclear matter using equations of state derived within the Brueckner-Hartree-Fock approach at finite temperature including consistent three-body forces. Simple and accurate parametrizations of the finite-temperature equations of state are provided. The properties of hot neutron stars are then investigated within this framework, in particular the temperature dependence of the maximum mass. We find very small temperature effects and analyze the interplay of the different contributions.
High-accuracy numerical simulations of merging neutron stars play an important role in testing and calibrating the waveform models used by gravitational wave observatories. Obtaining high-accuracy waveforms at a reasonable computational cost, however, remains a significant challenge. One issue is that high-order convergence of the solution requires the use of smooth evolution variables, while many of the equations of state used to model the neutron star matter have discontinuities, typically in the first derivative of the pressure. Spectral formulations of the equation of state have been proposed as a potential solution to this problem. Here, we report on the numerical implementation of spectral equations of state in the Spectral Einstein Code. We show that, in our code, spectral equations of state allow for high-accuracy simulations at a lower computational cost than commonly used `piecewise polytrope equations state. We also demonstrate that not all spectral equations of state are equally useful: different choices for the low-density part of the equation of state can significantly impact the cost and accuracy of simulations. As a result, simulations of neutron star mergers present us with a trade-off between the cost of simulations and the physical realism of the chosen equation of state.