We perform statistically rigorous uncertainty quantification (UQ) for chiral effective field theory ($\chi$EFT) applied to infinite nuclear matter up to twice nuclear saturation density. The equation of state (EOS) is based on high-order many-body perturbation theory calculations with consistent nucleon-nucleon and three-nucleon interactions up to fourth order in the $\chi$EFT expansion. From these calculations our newly developed Bayesian machine-learning approach extracts the size and smoothness properties of the correlated EFT truncation error. We then propose a novel extension that uses multi-task machine learning to reveal correlations between the EOS at different proton fractions. The inferred in-medium $\chi$EFT breakdown scale in pure neutron matter and symmetric nuclear matter is consistent with that from free-space nucleon-nucleon scattering. These significant advances allow us to provide posterior distributions for the nuclear saturation point and propagate theoretical uncertainties to derived quantities: the pressure and incompressibility of symmetric nuclear matter, the nuclear symmetry energy, and its derivative. Our results, which are validated by statistical diagnostics, demonstrate that an understanding of truncation-error correlations between different densities and different observables is crucial for reliable UQ. The methods developed here are publicly available as annotated Jupyter notebooks.

We introduce a new framework for quantifying correlated uncertainties of the infinite-matter equation of state (EOS) derived from chiral effective field theory ($\chi$EFT). Bayesian machine learning, via Gaussian Processes with physics-based hyperparameters, allows us to efficiently quantify and propagate theoretical uncertainties of the EOS, such as $\chi$EFT truncation errors, to derived quantities. We apply this framework to state-of-the-art many-body perturbation theory calculations with consistent nucleon-nucleon and three-nucleon interactions up to fourth order in the $\chi$EFT expansion. This produces the first statistically meaningful uncertainty estimates for key quantities of neutron stars. We give results up to twice nuclear saturation density for the energy per particle, pressure, and speed of sound of neutron matter, as well as for the nuclear symmetry energy and its derivative. At nuclear saturation density the predicted symmetry energy and its slope are consistent with experimental constraints.

We present an efficient Monte Carlo framework for perturbative calculations of infinite nuclear matter based on chiral two-, three-, and four-nucleon interactions. The method enables the incorporation of all many-body contributions in a straightforward and transparent way, and makes it possible to extract systematic uncertainty estimates by performing order-by-order calculations in the chiral expansion as well as the many-body expansion. The versatility of this new framework is demonstrated by applying it to chiral low-momentum interactions, exhibiting a very good many-body convergence up to fourth order. Following these benchmarks, we explore new chiral interactions up to next-to-next-to-next-to-leading order (N$^3$LO). Remarkably, simultaneous fits to the triton and to saturation properties can be achieved, while all three-nucleon low-energy couplings remain natural. The theoretical uncertainties of nuclear matter are significantly reduced when going from next-to-next-to-leading order to N$^3$LO.

We calculate the properties of isospin-asymmetric nuclear matter based on chiral nucleon-nucleon (NN) and three-nucleon (3N) interactions. To this end, we develop an improved normal-ordering framework that allows to include general 3N interactions starting from a plane-wave partial-wave-decomposed form. We present results for the energy per particle for general isospin asymmetries based on a set of different Hamiltonians, study their saturation properties, the incompressibility, symmetry energy, and also provide an analytic parametrization for the energy per particle as a function of density and isospin asymmetry.

Neutron matter is an ideal laboratory for nuclear interactions derived from chiral effective field theory since all contributions are predicted up to next-to-next-to-next-to-leading order (N$^3$LO) in the chiral expansion. By making use of recent advances in the partial-wave decomposition of three- nucleon (3N) forces, we include for the first time N$^3$LO 3N interactions in many-body perturbation theory (MBPT) up to third order and in self-consistent Green's function theory (SCGF). Using these two complementary many-body frameworks we provide improved predictions for the equation of state of neutron matter at zero temperature and also analyze systematically the many-body convergence for different chiral EFT interactions. Furthermore, we present an extension of the normal-ordering framework to finite temperatures. These developments open the way to improved calculations of neutron-rich matter including estimates of theoretical uncertainties for astrophysical applications.

We present solutions of the BCS gap equation in the channels ${}^1S_0$ and ${}^3P_2-{}^3F_2$ in neutron matter based on nuclear interactions derived within chiral effective field theory (EFT). Our studies are based on a representative set of nonlocal nucleon-nucleon (NN) plus three-nucleon (3N) interactions up to next-to-next-to-next-to-leading order (N$^3$LO) as well as local and semilocal chiral NN interactions up to N$^2$LO and N$^4$LO, respectively. In particular, we investigate for the first time the impact of subleading 3N forces at N$^3$LO on pairing gaps and also derive uncertainty estimates by taking into account results for pairing gaps at different orders in the chiral expansion. Finally, we discuss different methods for obtaining self-consistent solutions of the gap equation. Besides the widely-used quasi-linear method by Khodel et al. we demonstrate that the modified Broyden method is well applicable and exhibits a robust convergence behavior. In contrast to Khodel's method it is based on a direct iteration of the gap equation without imposing an auxiliary potential and is straightforward to implement.

We investigate asymmetric nuclear matter with two- and three-nucleon interactions based on chiral effective field theory, where three-body forces are fit only to light nuclei. Focusing on neutron-rich matter, we calculate the energy for different proton fractions and include estimates of the theoretical uncertainty. We use our ab-initio results to test the quadratic expansion around symmetric matter with the symmetry energy term, and confirm its validity for highly asymmetric systems. Our calculations are in remarkable agreement with an empirical parametrization for the energy density. These findings are very useful for astrophysical applications and for developing new equations of state.