Bayesian parameter estimation in $\chi$EFT using Hamiltonian Monte Carlo

Abstract

The number of low-energy constants (LECs) in chiral effective field theory ($\chi$EFT) grows rapidly with increasing chiral order, necessitating the use of Markov chain Monte Carlo techniques for sampling their posterior probability density function. For this we introduce a Hamiltonian Monte Carlo (HMC) algorithm and sample the LEC posterior up to next-to-next-to-leading order (NNLO) in the two-nucleon sector of $\chi$EFT. We find that the sampling efficiency of HMC is three to six times higher compared to an affine-invariant sampling algorithm. We analyze the empirical coverage probability and validate that the NNLO model yields predictions for two-nucleon scattering data with largely reliable credible intervals, provided that one ignores the leading order EFT expansion parameter when inferring the variance of the truncation error. We also find that the NNLO truncation error dominates the error budget.

22 pages, 12 figures

References in corpus (6)

• Array Programming with NumPy
• Chiral effective field theory and nuclear forces
• Global sensitivity analysis of bulk properties of an atomic nucleus
• Rigorous constraints on three-nucleon forces in chiral effective field theory from fast and accurate calculations of few-body observables
• Fast & accurate emulation of two-body scattering observables without wave functions
• Power counting in chiral effective field theory and nuclear binding

Cited by in corpus (14)

• Machine Learning in Nuclear Physics
• What is ab initio in nuclear theory?
• Model reduction methods for nuclear emulators
• Bayesian probability updates using Sampling/Importance Resampling: Applications in nuclear theory
• Inference of the low-energy constants in $\Delta$-full chiral effective field theory including a correlated truncation error
• Bayesian Analysis of $\chi$EFT at Leading Order in a Modified Weinberg Power Counting Approach
• Perturbative computations of neutron-proton scattering observables using renormalization-group invariant $\chi$EFT up to N${}^3$LO
• Bayesian approach for many-body uncertainties in nuclear structure: Many-body perturbation theory for finite nuclei
• Shell quenching in nuclear charge radii based on Monte Carlo dropout Bayesian neural network
• Posterior predictive distributions of neutron-deuteron cross sections
• Bayesian method for fitting the low-energy constants in chiral perturbation theory
• Bayesian estimation of the low-energy constants up to fourth order in the nucleon-nucleon sector of chiral effective field theory
• Meson-baryon scattering lengths without annihilation diagrams to order $p^4$ in heavy baryon chiral perturbation theory
• Chiral interactions and superfluidity in the calcium isotopic chain