We sample the posterior probability distributions of the low-energy constants (LECs) in $\mathrm{\Delta }$-full chiral effective field theory $\left(\chi \mathrm{EFT}\right)$ up to third order. We use eigenvector continuation for fast and accurate emulation of the likelihood and Hamiltonian Monte Carlo to draw effectively independent samples from the posteriors. Our Bayesian inference is conditioned on the Granada database of neutron-proton ($np$) cross sections and polarizations. We use priors grounded in $\chi \mathrm{EFT}$ assumptions and a Roy-Steiner analysis of pion-nucleon scattering data. We model correlated EFT truncation errors using a two-feature Gaussian process, and find correlation lengths for $np$ scattering energies and angles in the ranges 45–83 MeV and 24–39 degrees, respectively. These correlations yield a nondiagonal covariance matrix and reduce the number of independent scattering data with factors of 8 and 4 at the second and third chiral orders, respectively. The relatively small difference between the second- and third-order predictions in $\mathrm{\Delta }$-full $\chi \mathrm{EFT}$ suppresses the marginal variance of the truncation error and the effects of its correlation structure. Our results are particularly important for analyzing the predictive capabilities in ab initio nuclear theory.

• Received 27 April 2023
• Revised 7 March 2024
• Accepted 11 April 2024

DOI:https://doi.org/10.1103/PhysRevC.109.064003

Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI. Funded by Bibsam.

Published by the American Physical Society