An official website of the United States government
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
Toward Accelerated Nuclear-physics Parameter Estimation from Binary Neutron Star Mergers: Emulators for the Tolman–Oppenheimer–Volkoff Equations
Abstract Gravitational-wave observations of binary neutron-star (BNS) mergers have the potential to revolutionize our understanding of the nuclear equation of state (EOS) and the fundamental interactions that determine its properties. However, Bayesian parameter estimation frameworks do not typically sample over microscopic nuclear-physics parameters that determine the EOS. One of the major hurdles in doing so is the computational cost involved in solving the neutron-star structure equations, known as the Tolman–Oppenheimer–Volkoff (TOV) equations. In this paper, we explore approaches to emulating solutions for the TOV equations: multilayer perceptrons (MLPs), Gaussian processes, and a data-driven variant of the reduced basis method (RBM). We implement these emulators for three different parameterizations of the nuclear EOS, each with a different degree of complexity represented by the number of model parameters. We find that our MLP-based emulators are generally more accurate than the other two algorithms, whereas the RBM results in the largest speedup with respect to the full high-fidelity TOV solver. We employ these emulators for a simple parameter inference using a potentially loud BNS observation and show that the posteriors predicted by our emulators are in excellent agreement with those obtained from the full TOV solver.
more »
« less
- Award ID(s):
- 2412341
- PAR ID:
- 10611555
- Publisher / Repository:
- American Astronomical Society
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 974
- Issue:
- 2
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 285
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Binary neutron star mergers (NSMs) have been confirmed as one source of the heaviest observable elements made by the rapid neutron-capture ( r -) process. However, modeling NSM outflows—from the total ejecta masses to their elemental yields—depends on the unknown nuclear equation of state (EOS) that governs neutron star structure. In this work, we derive a phenomenological EOS by assuming that NSMs are the dominant sources of the heavy element material in metal-poor stars with r -process abundance patterns. We start with a population synthesis model to obtain a population of merging neutron star binaries and calculate their EOS-dependent elemental yields. Under the assumption that these mergers were responsible for the majority of r -process elements in the metal-poor stars, we find parameters representing the EOS for which the theoretical NSM yields reproduce the derived abundances from observations of metal-poor stars. For our proof-of-concept assumptions, we find an EOS that is slightly softer than, but still in agreement with, current constraints, e.g., by the Neutron Star Interior Composition Explorer, with R 1.4 = 12.25 ± 0.03 km and M TOV = 2.17 ± 0.03 M ⊙ (statistical uncertainties, neglecting modeling systematics).more » « less
-
Kilonova Emissions from Neutron Star Merger Remnants: Implications for the Nuclear Equation of StateAbstract Multimessenger observations of binary neutron star mergers can provide valuable information on the nuclear equation of state (EOS). Here, we investigate the extent to which electromagnetic observations of the associated kilonovae allow us to place constraints on the EOS. For this, we use state-of-the-art three-dimensional general-relativistic magnetohydrodynamics simulations and detailed nucleosynthesis modeling to connect properties of observed light curves to properties of the accretion disk, and hence, the EOS. Using our general approach, we use multimessenger observations of GW170817/AT2017gfo to study the impact of various sources of uncertainty on inferences of the EOS. We constrain the radius of a 1.4M⊙neutron star to lie within 10.30 ≤R1.4≤ 13.0 km and the maximum mass to beMTOV≤ 3.06M⊙.more » « less
-
Abstract This study investigates the properties of symmetric and asymmetric nuclear matter using the relativistic Brueckner–Hartree–Fock formalism, examining both zero and finite temperatures up to 70 MeV. Employing the full Dirac space, we incorporate three Bonn potentials (A, B, and C), which account for meson masses, coupling strengths, cutoff parameters, and form factors. The calculated properties of asymmetric nuclear matter form the basis for constructing equation-of-state (EOS) models tailored for neutron stars. These models, in turn, enable the computation of bulk properties for nonrotating, uniformly rotating, and differentially rotating neutron stars. Notably, the EOS models studied in this paper are sufficiently versatile to accommodate the mass of the most massive neutron star ever detected, PSR J0952–0607, estimated to be 2.35 ± 0.17M⊙. Furthermore, they yield masses and radii for PSR J0030+451 that align with the confidence intervals established for this pulsar.more » « less
-
null (Ed.)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
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.3847/1538-4357/ad737c
