Neutron star crust and outer core equation of state from chiral effective field theory with quantified uncertainties

TL;DR

This paper develops a Bayesian framework using Gaussian processes to quantify uncertainties in the neutron star crust and outer core equation of state derived from chiral effective field theory.

Contribution

It introduces a 2D Gaussian process trained on many-body perturbation theory results to efficiently evaluate the EOS with uncertainty quantification.

Findings

The method accurately benchmarks EFT uncertainties for neutron matter.

Constructs EOSs for neutron star crust consistent with chiral EFT results.

Analyzes the phase diagram of neutron-rich matter including surface and Coulomb effects.

Abstract

We study the order-by-order expansion of the energy per particle of asymmetric nuclear matter up to twice saturation density in chiral effective field theory (EFT) within a Bayesian framework. For this, we develop a two-dimensional Gaussian process (2D GP) that is trained using many-body perturbation theory results based on chiral two- and three-nucleon interactions from leading to next-to-next-to-next-to-leading order (N$^3$LO). This allows for an efficient evaluation of the equation of state (EOS) and thermodynamic derivatives with EFT truncation uncertainties. After benchmarking our 2D GP against Bayesian uncertainties for pure neutron matter and symmetric matter, we study the energy per particle, pressure, and chemical potentials of neutron star matter in $\beta$-equilibrium including EFT uncertainties. We investigate the phase diagram of neutron-rich matter from neutron- to proton-drip and to the uniform phase, including surface and Coulomb corrections. Based on this, we construct EOSs for the inner crust of neutron stars that are consistent with the chiral EFT results for uniform matter at N$^3$LO.