Inferring the neutron star equation of state with nuclear-physics informed semiparametric models

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 (NSs), substantial uncertainty about the dense matter equation of state (EoS) remains. In this work, we present a semiparametric EoS 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 (GP) 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 GP 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 (PSRs). We find that maximum Tolman–Oppenheimer–Volkoff masses above 3.2 M ⊙ can be supported by causal EoS compatible with nuclear constraints at low densities. We then combine information from existing observations of heavy PSR masses, gravitational waves emitted from binary NS mergers, and x-ray pulse profile modeling of millisecond PSRs within a Bayesian inference scheme using our semiparametric EoS prior. With information from all public NS Interior Composition ExploRer PSRs (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 ρ nuc ) = 1.98 − 1.08 + 2.13 × 10 34 dyn cm −2 , a radius of a 1.4 M ⊙ NS value of R 1.4 = 11.4 − 0.60 + 0.98 km, and M max = 2.31 − 0.23 + 0.35 M ⊙ at the 90% credible level (C.L).