Emulating \emph{ab initio} computations of infinite nucleonic matter

W. G. Jiang,C. Forssén,T. Djärv,G. Hagen

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

2024-06-16

Abstract:We construct efficient emulators for the \emph{ab initio} computation of the infinite nuclear matter equation of state. These emulators are based on the subspace-projected coupled-cluster method for which we here develop a new algorithm called small-batch voting to eliminate spurious states that might appear when emulating quantum many-body methods based on a non-Hermitian Hamiltonian. The efficiency and accuracy of these emulators facilitate a rigorous statistical analysis within which we explore nuclear matter predictions for $> 10^6$ different parametrizations of a chiral interaction model with explicit $\Delta$-isobars at next-to-next-to leading order. Constrained by nucleon-nucleon scattering phase shifts and bound-state observables of light nuclei up to \nuc{4}{He}, we use history matching to identify non-implausible domains for the low-energy coupling constants of the chiral interaction. Within these domains we perform a Bayesian analysis using sampling/importance resampling with different likelihood calibrations and study correlations between interaction parameters, calibration observables in light nuclei, and nuclear matter saturation properties.

Nuclear Theory

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to simulate the equation of state (EOS) of infinite nuclear matter by constructing efficient emulators. Specifically, the researchers developed a new algorithm - small - batch voting - to eliminate spurious states that may occur when using non - Hermitian Hamiltonians for quantum many - body method simulations. These emulators can efficiently and accurately predict the properties of nuclear matter under more than 10^6 different chiral interaction model parameterizations. In addition, through the history matching technique, the researchers were able to identify the non - unreasonable domains of low - energy coupling constants (LECs), and on this basis, conduct Bayesian analysis to study the correlations between different interaction parameters, as well as the relationships between these parameters and the calibration observables of light nuclei and the saturation properties of nuclear matter.

The key contributions of the paper are as follows:
1. **Developed the small - batch voting algorithm**: It is used to eliminate spurious states that occur during the simulation process, improving the accuracy and reliability of the emulators.
2. **Constructed efficient emulators**: These models are based on the subspace - projected coupled - cluster method (SPCC) and can maintain high - precision predictions of the nuclear matter equation of state while significantly reducing the computational cost.
3. **Conducted large - scale statistical analysis**: By analyzing more than 10^6 different chiral interaction model parameterizations, the multi - dimensional parameter space of nuclear matter properties was explored, providing a new perspective for understanding the basic properties of nuclear matter.
4. **Applied history matching and Bayesian analysis**: Combined with experimental data, the reasonable parameter range was narrowed down through the history matching technique, and the correlations and uncertainties between parameters were further studied through Bayesian analysis.

In summary, this paper has solved multiple challenges in the calculation of the nuclear matter equation of state by developing new algorithms and techniques, providing important tools and methods for research in the field of nuclear physics.