Emulators for scarce and noisy data II: Application to auxiliary-field diffusion Monte Carlo for neutron matter

Abstract: Understanding the equation of state (EOS) of pure neutron matter is necessary for interpreting multi-messenger observations of neutron stars. Reliable data analyses of these observations require well-quantified uncertainties for the EOS input, ideally propagating uncertainties from nuclear interactions directly to the EOS. This, however, requires calculations of the EOS for a prohibitively larger number of nuclear Hamiltonians, solving the nuclear many-body problem for each one. Quantum Monte Carlo methods, such as Auxiliary field diffusion Monte Carlo (AFDMC), provide precise and accurate results for the neutron matter EOS but they are very computationally expensive, making them unsuitable for the fast evaluations necessary for uncertainty propagation. Here, we employ parametric matrix models to develop fast emulators for AFDMC calculations of neutron matter, and use them to directly propagate uncertainties of coupling constants in the Hamiltonian to the EOS. As these uncertainties include estimates of the EFT truncation uncertainty, this approach provides robust uncertainty estimates for use in astrophysical data analyses. This work will enable novel applications such as using astrophysical observations to put constraints on coupling constants for nuclear interactions.