skip to main content


Title: Realistic Finite-Temperature Effects in Neutron Star Merger Simulations
Binary neutron star mergers provide a unique probe of the dense-matter equation of state (EoS) across a wide range of parameter space, from the zero-temperature EoS during the inspiral to the high-temperature EoS following the merger. In this paper, we implement a new model for calculating parametrized finite-temperature EoS effects into numerical relativity simulations. This "M* model" is based on a two-parameter approximation of the particle effective mass and includes the leading-order effects of degeneracy in the thermal pressure and energy. We test our numerical implementation by performing evolutions of rotating single stars with zero- and non-zero temperature gradients, as well as evolutions of binary neutron star mergers. We find that our new finite-temperature EoS implementation can support stable stars over many dynamical timescales. We also perform a first parameter study to explore the role of the M* parameters in binary neutron star merger simulations. All simulations start from identical initial data with identical cold EoSs, and differ only in the thermal part of the EoS. We find that both the thermal profile of the remnant and the post-merger gravitational wave signal depend on the choice of M* parameters, but that the total merger ejecta depends only weakly on the finite-temperature part of the EoS across a wide range of parameters. Our simulations provide a first step toward understanding how the finite-temperature properties of dense matter may affect future observations of binary neutron star mergers.  more » « less
Award ID(s):
1912619
NSF-PAR ID:
10276495
Author(s) / Creator(s):
Date Published:
Journal Name:
Physical review
ISSN:
2470-0010
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    Parametric equations of state (EoSs) provide an important tool for systematically studying EoS effects in neutron star merger simulations. In this work, we perform a numerical validation of the M*-framework for parametrically calculating finite-temperature EoS tables. The framework, introduced by Raithel et al., provides a model for generically extending any cold, β-equilibrium EoS to finite temperatures and arbitrary electron fractions. In this work, we perform numerical evolutions of a binary neutron star merger with the SFHo finite-temperature EoS, as well as with the M*-approximation of this same EoS, where the approximation uses the zero-temperature, β-equilibrium slice of SFHo and replaces the finite-temperature and composition-dependent parts with the M*-model. We find that the approximate version of the EoS is able to accurately recreate the temperature and thermal pressure profiles of the binary neutron star remnant, when compared to the results found using the full version of SFHo. We additionally find that the merger dynamics and gravitational wave signals agree well between both cases, with differences of $\lesssim 1\!-\!2\,{\textrm{per cent}}$ introduced into the post-merger gravitational wave peak frequencies by the approximations of the EoS. We conclude that the M*-framework can be reliably used to probe neutron star merger properties in numerical simulations.

     
    more » « less
  2. Abstract

    The equation of state (EOS) of dense strongly interacting matter can be probed by astrophysical observations of neutron stars (NS), such as X-ray detections of pulsars or the measurement of the tidal deformability of NSs during the inspiral stage of NS mergers. These observations constrain the EOS at most up to the density of the maximum-mass configuration,nTOV, which is the highest density that can be explored by stable NSs for a given EOS. However, under the right circumstances, binary neutron star (BNS) mergers can create a postmerger remnant that explores densities abovenTOV. In this work, we explore whether the EOS abovenTOVcan be measured from gravitational-wave or electromagnetic observations of the postmerger remnant. We perform a total of 25 numerical-relativity simulations of BNS mergers for a range of EOSs and find no case in which different descriptions of the matter abovenTOVhave a detectable impact on postmerger observables. Hence, we conclude that the EOS abovenTOVcan likely not be probed through BNS merger observations for the current and next generation of detectors.

     
    more » « less
  3. ABSTRACT

    GW190425 was the second gravitational wave (GW) signal compatible with a binary neutron star (BNS) merger detected by the Advanced LIGO and Advanced Virgo detectors. Since no electromagnetic counterpart was identified, whether the associated kilonova was too dim or the localization area too broad is still an open question. We simulate 28 BNS mergers with the chirp mass of GW190425 and mass ratio 1 ≤ q ≤ 1.67, using numerical-relativity simulations with finite-temperature, composition dependent equations of state (EOS) and neutrino radiation. The energy emitted in GWs is $\lesssim 0.083\mathrm{\, M_\odot }c^2$ with peak luminosity of 1.1–$2.4\times ~10^{58}/(1+q)^2\, {\rm {erg \, s^{-1}}}$. Dynamical ejecta and disc mass range between 5 × 10−6–10−3 and 10−5–$0.1 \mathrm{\, M_\odot }$, respectively. Asymmetric mergers, especially with stiff EOSs, unbind more matter and form heavier discs compared to equal mass binaries. The angular momentum of the disc is 8–$10\mathrm{\, M_\odot }~GM_{\rm {disc}}/c$ over three orders of magnitude in Mdisc. While the nucleosynthesis shows no peculiarity, the simulated kilonovae are relatively dim compared with GW170817. For distances compatible with GW190425, AB magnitudes are always dimmer than ∼20 mag for the B, r, and K bands, with brighter kilonovae associated to more asymmetric binaries and stiffer EOSs. We suggest that, even assuming a good coverage of GW190425’s sky location, the kilonova could hardly have been detected by present wide-field surveys and no firm constraints on the binary parameters or EOS can be argued from the lack of the detection.

     
    more » « less
  4. null (Ed.)
    In this work, we discuss the dense matter equation of state (EOS) for the extreme range of conditions encountered in neutron stars and their mergers. The calculation of the properties of such an EOS involves modeling different degrees of freedom (such as nuclei, nucleons, hyperons, and quarks), taking into account different symmetries, and including finite density and temperature effects in a thermodynamically consistent manner. We begin by addressing subnuclear matter consisting of nucleons and a small admixture of light nuclei in the context of the excluded volume approach. We then turn our attention to supranuclear homogeneous matter as described by the Chiral Mean Field (CMF) formalism. Finally, we present results from realistic neutron-star-merger simulations performed using the CMF model that predict signatures for deconfinement to quark matter in gravitational wave signals. 
    more » « less
  5. Abstract

    Neutrinos are copiously emitted by neutron star mergers, due to the high temperatures reached by dense matter during the merger and its aftermath. Neutrinos influence the merger dynamics and shape the properties of the ejecta, including the resultingr-process nucleosynthesis and kilonova emission. In this work, we analyse neutrino emission from a large sample of binary neutron star merger simulations in Numerical Relativity, covering a broad range of initial masses, nuclear equation of state and viscosity treatments. We extract neutrino luminosities and mean energies, and compute quantities of interest such as the peak values, peak broadnesses, time averages and decrease time scales. We provide a systematic description of such quantities, including their dependence on the initial parameters of the system. We find that for equal-mass systems the total neutrino luminosity (several$$10^{53}{\hbox {erg}~{\hbox {s}}^{-1}}$$1053ergs-1) decreases as the reduced tidal deformability increases, as a consequence of the less violent merger dynamics. Similarly, tidal disruption in asymmetric mergers leads to systematically smaller luminosities. Peak luminosities can be twice as large as the average ones. Electron antineutrino luminosities dominate (initially by a factor of 2-3) over electron neutrino ones, while electron neutrinos and heavy flavour neutrinos have similar luminosities. Mean energies are nearly constant in time and independent on the binary parameters. Their values reflect the different decoupling temperature inside the merger remnant. Despite present uncertainties in neutrino modelling, our results provide a broad and physically grounded characterisation of neutrino emission, and they can serve as a reference point to develop more sophisticated neutrino transport schemes.

     
    more » « less