We construct an efficient emulator for two-body scattering observables using the general (complex) Kohn variational principle and trial wave functions derived from eigenvector continuation. The emulator simultaneously evaluates an array of Kohn varia
tional principles associated with different boundary conditions, which allows for the detection and removal of spurious singularities known as Kohn anomalies. When applied to the $K$-matrix only, our emulator resembles the one constructed by Furnstahl et al. [Phys. Lett. B 809, 135719] although with reduced numerical noise. After a few applications to real potentials, we emulate differential cross sections for $^{40}$Ca$(n,n)$ scattering based on a realistic optical potential and quantify the model uncertainties using Bayesian methods. These calculations serve as a proof of principle for future studies aimed at improving optical models.
We combine Newtons variational method with ideas from eigenvector continuation to construct a fast & accurate emulator for two-body scattering observables. The emulator will facilitate the application of rigorous statistical methods for interactions
that depend smoothly on a set of free parameters. Our approach begins with a trial $K$ or $T$ matrix constructed from a small number of exact solutions to the Lippmann--Schwinger equation. Subsequent emulation only requires operations on small matrices. We provide several applications to short-range potentials with and without the Coulomb interaction and partial-wave coupling. It is shown that the emulator can accurately extrapolate far from the support of the training data. When used to emulate the neutron-proton cross section with a modern chiral interaction as a function of 26 free parameters, it reproduces the exact calculation with negligible error and provides an over 300x improvement in CPU time.
We discuss high-order calculations in perturbative effective field theory for fermions at low energy scales. The Fermi-momentum or $k_{rm F} a_s$ expansion for the ground-state energy of the dilute Fermi gas is calculated to fourth order, both in cut
off regularization and in dimensional regularization. For the case of spin one-half fermions we find from a Bayesian analysis that the expansion is well-converged at this order for ${| k_{rm F} a_s | lesssim 0.5}$. Further, we show that Pad{e}-Borel resummations can improve the convergence for ${| k_{rm F} a_s | lesssim 1}$. Our results provide important constraints for nonperturbative calculations of ultracold atoms and dilute neutron matter.
Born in the aftermath of core collapse supernovae, neutron stars contain matter under extraordinary conditions of density and temperature that are difficult to reproduce in the laboratory. In recent years, neutron star observations have begun to yiel
d novel insights into the nature of strongly interacting matter in the high-density regime where current theoretical models are challenged. At the same time, chiral effective field theory has developed into a powerful framework to study nuclear matter properties with quantified uncertainties in the moderate-density regime for modeling neutron stars. In this article, we review recent developments in chiral effective field theory and focus on many-body perturbation theory as a computationally efficient tool for calculating the properties of hot and dense nuclear matter. We also demonstrate how effective field theory enables statistically meaningful comparisons between nuclear theory predictions, nuclear experiments, and observational constraints on the nuclear equation of state.
The nuclear symmetry energy is a key quantity in nuclear (astro)physics. It describes the isospin dependence of the nuclear equation of state (EOS), which is commonly assumed to be almost quadratic. In this work, we confront this standard quadratic e
xpansion of the EOS with explicit asymmetric nuclear-matter calculations based on a set of commonly used Hamiltonians including two- and three-nucleon forces derived from chiral effective field theory. We study, in particular, the importance of non-quadratic contributions to the symmetry energy, including the non-analytic logarithmic term introduced by Kaiser [Phys.~Rev.~C textbf{91}, 065201 (2015)]. Our results suggest that the quartic contribution to the symmetry energy can be robustly determined from the various Hamiltonians employed, and we obtain 1.00(8) MeV (or 0.55(8) MeV for the potential part) at saturation density, while the logarithmic contribution to the symmetry energy is relatively small and model-dependent. We finally employ the meta-model approach to study the impact of the higher-order contributions on the neutron-star crust-core transition density, and find a small 5% correction.
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 per
turbation theory calculations with 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 multitask 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 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 equation of state, such as $chi$EFT truncation errors, to derived quantities. We apply this framework to state-of-the-art many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the $chi$EFT expansion. This produces the first statistically robust 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 study ground-state energies and charge radii of closed-shell medium-mass nuclei based on novel chiral nucleon-nucleon (NN) and three-nucleon (3N) interactions, with a focus on exploring the connections between finite nuclei and nuclear matter. To
this end, we perform in-medium similarity renormalization group (IM-SRG) calculations based on chiral interactions at next-to-leading order (NLO), N$^2$LO, and N$^3$LO, where the 3N interactions at N$^2$LO and N$^3$LO are fit to the empirical saturation point of nuclear matter and to the triton binding energy. Our results for energies and radii at N$^2$LO and N$^3$LO overlap within uncertainties, and the cutoff variation of the interactions is within the EFT uncertainty band. We find underbound ground-state energies, as expected from the comparison to the empirical saturation point. The radii are systematically too large, but the agreement with experiment is better. We further explore variations of the 3N couplings to test their sensitivity in nuclei. While nuclear matter at saturation density is quite sensitive to the 3N couplings, we find a considerably weaker dependence in medium-mass nuclei. In addition, we explore a consistent momentum-space SRG evolution of these NN and 3N interactions, exhibiting improved many-body convergence. For the SRG-evolved interactions, the sensitivity to the 3N couplings is found to be stronger in medium-mass nuclei.
Using effective field theory methods, we calculate for the first time the complete fourth-order term in the Fermi-momentum or $k_{rm F} a_s$ expansion for the ground-state energy of a dilute Fermi gas. The convergence behavior of the expansion is exa
mined for the case of spin one-half fermions and compared against quantum Monte-Carlo results, showing that the Fermi-momentum expansion is well-converged at this order for $| k_{rm F} a_s | lesssim 0.5$.
