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

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.