An official website of the United States government
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
Inference of neutron-star properties with unified crust-core equations of state for parameter estimation
Context.Relating different global neutron-star (NS) properties, such as tidal deformability and radius, or mass and radius, requires an equation of state (EoS). Determining the NS EoS is therefore not only the science goal of a variety of observational projects, but it also enters in the analysis process; for example, to predict a NS radius from a measured tidal deformability via gravitational waves (GW) during the inspiral of a binary NS merger. To this aim, it is important to estimate the theoretical uncertainties on the EoS, one of which is the possible bias coming from an inconsistent treatment of the low-density region; that is, the use of a so called non-unified NS crust. Aims.We propose a numerical tool allowing the user to consistently match a nuclear-physics informed crust to an arbitrary high-density EoS describing the core of the star. Methods.We introduce an inversion procedure of the EoS close to saturation density that allows users to extract nuclear-matter parameters and extend the EoS to lower densities in a consistent way. For the treatment of inhomogeneous matter in the crust, a standard approach based on the compressible liquid-drop (CLD) model approach was used in our work. A Bayesian analysis using a parametric agnostic EoS representation in the high-density region is also presented in order to quantify the uncertainties induced by an inconsistent treatment of the crust. Results.We show that the use of a fixed, realistic-but-inconsistent model for the crust causes small but avoidable errors in the estimation of global NS properties and leads to an underestimation of the uncertainties in the inference of NS properties. Conclusions.Our results highlight the importance of employing a consistent EoS in inference schemes. The numerical tool that we developed to reconstruct such a thermodynamically consistent EoS, CUTER, has been tested and validated for use by the astrophysical community.
more »
« less
- Award ID(s):
- 2116686
- PAR ID:
- 10527743
- Publisher / Repository:
- EDP Sciences
- Date Published:
- Journal Name:
- Astronomy & Astrophysics
- Volume:
- 687
- ISSN:
- 0004-6361
- Page Range / eLocation ID:
- A44
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Over the past decade, an abundance of information from neutron-star observations, nuclear experiments and theory has transformed our efforts to elucidate the properties of dense matter. However, at high densities relevant to the cores of neutron stars, substantial uncertainty about the dense matter equation of state (EoS) remains. In this work, we present a semiparametric equation of state framework aimed at better integrating knowledge across these domains in astrophysical inference. We use a Meta-model and realistic crust at low densities, and Gaussian Process extensions at high densities. Comparisons between our semiparametric framework to fully nonparametric EoS representations show that imposing nuclear theoretical and experimental constraints through the Meta-model up to nuclear saturation density results in constraints on the pressure up to twice nuclear saturation density. We also show that our Gaussian Process trained on EoS models with nucleonic, hyperonic, and quark compositions extends the range of EoS explored at high density compared to a piecewise polytropic extension schema, under the requirements of causality of matter and of supporting the existence of heavy pulsars. We find that maximum TOV masses above $$3.2 M_{\odot}$$ can be supported by causal EoS compatible with nuclear constraints at low densities. We then combine information from existing observations of heavy pulsar masses, gravitational waves emitted from binary neutron star mergers, and X-ray pulse profile modeling of millisecond pulsars within a Bayesian inference scheme using our semiparametric EoS prior. With information from all public NICER pulsars (including PSR J0030$$+$$0451, PSR J0740$$+$$6620, PSR J0437-4715, and PSR J0614-3329), we find an astrophysically favored pressure at two times nuclear saturation density of $$P(2\rho_{\rm nuc}) = 1.98^{+2.13}_{-1.08}\times10^{34}$$ dyn/cm$$^{2}$$, a radius of a $$1.4 M_{\odot}$$ neutron star value of $$R_{1.4} = 11.4^{+0.98}_{-0.60}$$\;km, and $$M_{\rm max} = 2.31_{-0.23}^{+0.35} M_{\odot}$$ at the 90\% credible level.more » « less
-
Abstract Neutron stars provide a unique laboratory for studying matter at extreme pressures and densities. While there is no direct way to explore their interior structure, X-rays emitted from these stars can indirectly provide clues to the equation of state (EOS) of the superdense nuclear matter through the inference of the star's mass and radius. However, inference of EOS directly from a star's X-ray spectra is extremely challenging and is complicated by systematic uncertainties. The current state of the art is to use simulation-based likelihoods in a piece-wise method which relies on certain theoretical assumptions and simplifications about the uncertainties. It first infers the star's mass and radius to reduce the dimensionality of the problem, and from those quantities infer the EOS. We demonstrate a series of enhancements to the state of the art, in terms of realistic uncertainty quantification and a path towards circumventing the need for theoretical assumptions to infer physical properties with machine learning. We also demonstrate novel inference of the EOS directly from the high-dimensional spectra of observed stars, avoiding the intermediate mass-radius step. Our network is conditioned on the sources of uncertainty of each star, allowing for natural and complete propagation of uncertainties to the EOS.more » « less
-
Due to the high-density nuclear matter equation of state (EOS) being as yet unknown, neutron stars (NSs) do not have a confirmed limiting “Chandrasekhar” type maximum mass. However, observations of NSs (PSR J1614-2230, PSR J0348+0432, PSR J0740+6620, PSR J0952–0607) indicate that the NS’s limiting mass, if there is any, could be well over 2M⊙. On the other hand, there seems to be an observational mass gap (of around 2.5−5M⊙ ) between the lightest black hole and the heaviest NS. Therefore, the “massive NSs” are prime candidates to fill that gap. Several NS EOSs have been proposed using both microscopic and phenomenological approaches. In this project, we look at a class of phenomenological nuclear matter EOSs—relativistic mean field models—and see what kind of NS is formed from them. We compute the max- imum mass supported by each model EOS to observe if the mass of the NS is indeed in the “massive” NS (>2M⊙) regime. We also observe the effects of including exotic particles (hyperons, Δs) in the NS EOS and how that affects the NS mass. However, only by introducing the magnetic field, i.e. for magnetized anisotro- pic NSs, we find the mass exceeding 2.5M⊙. Using tidal deformability constraints from gravitational wave observations, we place a further check on how physical the EOS and NSs are.more » « less
-
Abstract The detection of GW170817, the first neutron star-neutron star merger observed by Advanced LIGO and Virgo, and its following analyses represent the first contributions of gravitational wave data to understanding dense matter. Parameterizing the high density section of the equation of state of both neutron stars through spectral decomposition, and imposing a lower limit on the maximum mass value, led to an estimate of the stars’ radii ofkm andkm (Abbottet al2018Phys. Rev. Lett.121161101). These values do not, however, take into account any uncertainty owed to the choice of the crust low-density equation of state, which was fixed to reproduce the SLy equation of state model (Douchin and Haensel 2001Astron. Astrophys.380151). We here re-analyze GW170817 data and establish that different crust models do not strongly impact the mass or tidal deformability of a neutron star—it is impossible to distinguish between low-density models with gravitational wave analysis. However, the crust does have an effect on inferred radius. We predict the systematic error due to this effect using neutron star structure equations, and compare the prediction to results from full parameter estimation runs. For GW170817, this systematic error affects the radius estimate by 0.3 km, approximatelyof the neutron stars’ radii.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1051/0004-6361/202348402
