Unified nonparametric equation-of-state inference from the neutron-star crust to perturbative-QCD densities

TL;DR

This paper develops a unified nonparametric framework to incorporate perturbative QCD constraints into neutron-star equation-of-state inference, improving modeling consistency across densities from the crust to core.

Contribution

It introduces a unified EOS model that seamlessly integrates pQCD constraints as priors, avoiding arbitrary density choices and enhancing the physical consistency of neutron-star matter modeling.

Findings

pQCD constraints soften the EOS at high densities

No significant impact on maximum neutron-star mass observed

Speed of sound in neutron-star cores is bounded by ~0.5 at 90% confidence

Abstract

Perturbative quantum chromodynamics (pQCD), while valid only at densities exceeding those found in the cores of neutron stars, could provide constraints on the dense-matter equation of state (EOS). In this work, we examine the impact of pQCD information on the inference of the EOS using a nonparametric framework based on Gaussian processes (GPs). We examine the application of pQCD constraints through a "pQCD likelihood," and verify the findings of previous works; namely, a softening of the EOS at the central densities of the most massive neutron stars and a reduction in the maximum neutron-star mass. Although the pQCD likelihood can be easily integrated into existing EOS inference frameworks, this approach requires an arbitrary selection of the density at which the constraints are applied. The EOS behavior is also treated differently on either side of the chosen density. To mitigate these issues, we extend the EOS model to higher densities, thereby constructing a "unified" description of the EOS from the neutron-star crust to densities relevant for pQCD. In this approach the pQCD constraints effectively become part of the prior. Since the EOS is unconstrained by any calculation or data between the densities applicable to neutron stars and pQCD, we argue for maximum modeling flexibility in that regime. We compare the unified EOS with the traditional pQCD likelihood, and although we confirm the EOS softening, we do not see a reduction in the maximum neutron-star mass or any impact on macroscopic observables. Though residual model dependence cannot be ruled out, we find that pQCD suggests the speed of sound in the densest neutron-star cores has already started decreasing toward the asymptotic limit; we find that the speed of sound squared at the center of the most massive neutron star has an upper bound of $\sim 0.5$ at the $90\%$ level.