skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: New Neutron Star Equation of State with Quark–Hadron Crossover
We present a much improved equation of state for neutron star matter, QHC19, with a smooth crossover from the hadronic regime at lower densities to the quark regime at higher densities. We now use the Togashi et al.equation of state, a generalization of the Akmal–Pandharipande–Ravenhall equation of state of uniform nuclear matter, in the entire hadronic regime; the Togashi equation of state consistently describes nonuniform as well as uniform matter, and matter at beta equilibrium without the need for an interpolation between pure neutron and symmetric nuclear matter. We describe the quark matter regime at higher densities with the Nambu–Jona–Lasinio model, now identifying tight constraints on the phenomenological universal vector repulsion between quarks and the pairing interaction between quarks arising from the requirements of thermodynamic stability and causal propagation of sound. The resultant neutron star properties agree very well with the inferences of the LIGO/Virgo collaboration, from GW170817, of the pressure versus baryon density, neutron star radii, and tidal deformabilities. The maximum neutron star mass allowed by QHC19 is 2.35 solar masses, consistent with all neutron star mass determinations.  more » « less
Award ID(s):
1714042
PAR ID:
10126214
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
The Astrophysical journal
Volume:
885
Issue:
42
ISSN:
2471-4259
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The nonlocal three-flavor Nambu-Jona-Lasinio model is used to study quark deconfinement in the cores of neutron stars (NSs). The quark-hadron phase transition is modeled using both the Maxwell construction and the Gibbs construction. For the Maxwell construction, we find that all NSs with core densities beyond the phase transition density are unstable. Therefore, no quark matter cores would exist inside such NSs. The situation is drastically different if the phase transition is treated as a Gibbs transition, resulting in stable NSs whose stellar cores are a mixture of hadronic matter and deconfined quarks. The largest fractions of quarks achieved in the quark-hadron mixed phase are around 50%. No choice of parametrization or composition leads to a pure quark matter core. The inclusion of repulsive vector interactions among the quarks is crucial since the equation of state (EoS) in the quark-hadron mixed phase is significantly softer than that of the pure hadronic phase. 
    more » « less
  2. Abstract: With recent advances in astronomical observations, major progress has been made in determining the pressure of neutron star matter at high density. This pressure is constrained by the neutron star deformability, determined from gravitational waves emitted in a neutron-star merger, and the mass-radii relation of two neutron stars, determined from a new X-ray observatory on the International Space Station. Previous studies have relied on nuclear theory calculations to constrain the equation of state at low density. Here we use a combination of constraints composed of three astronomical observations and twelve nuclear experimental constraints that extend over a wide range of densities. A Bayesian inference framework is then used to obtain a comprehensive nuclear equation of state. This data-centric result provides benchmarks for theoretical calculations and modeling of nuclear matter and neutron stars. Furthermore, it provides insights into the microscopic degrees of freedom of the nuclear matter equation of state and on the composition of neutron stars and their cooling via neutrino radiation. 
    more » « less
  3. Neutron stars were first posited in the early thirties and discovered as pulsars in late sixties; however, only recently are we beginning to understand the matter they contain. This talk describes the continuing development of a consistent picture of the liquid interiors of neutron stars, driven by four advances: observations of heavy neutron stars with masses in the range of two solar masses; inferences of masses and radii simultaneously for an increasing number of neutron stars in low mass X-ray binaries, and ongoing determinations via the NICER observatory; the observation of the binary neutron star merger, GW170817, through gravitational waves as well as across the electromagnetic spectrum; and an emerging understanding in QCD of how nuclear matter can turn into deconfined quark matter in the interior. We describe the modern quark-hadron crossover equation of state, QHC18 and now QHC19, and the corresponding neutron stars, which agree well with current observations. 
    more » « less
  4. Recent experimental and ab initio theory investigations of the 208Pb neutron skin thickness have the potential to inform the neutron star equation of state. In particular, the strong correlation between the 208Pb neutron skin thickness and the pressure of neutron matter at normal nuclear densities leads to modified predictions for the radii, tidal deformabilities, and moments of inertia of typical 1.4M⊙ neutron stars. In the present work, we study the relative impact of these recent analyses of the 208Pb neutron skin thickness on bulk properties of neutron stars within a Bayesian statistical analysis. Two models for the equation of state prior are employed in order to highlight the role of the highly uncertain high-density equation of state. From our combined Bayesian analysis of nuclear theory, nuclear experiment, and observational constraints on the dense matter equation of state, we find at the 90% credibility level R1.4=12.36−0.73+0.38 km for the radius of a 1.4M⊙ neutron star, R2.0=11.96−0.71+0.94 km for the radius of a 2.0M⊙ neutron star, Λ1.4=440−144+103 for the tidal deformability of a 1.4M⊙ neutron star, and I1.338=1.425−0.146+0.074×1045gcm2 for the moment of inertia of PSR J0737-3039A whose mass is 1.338M⊙. 
    more » « less
  5. 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