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.

arXiv comment22 pages, 12 figures

References in corpus (26)

• Array Programming with NumPy
• emcee: The MCMC Hammer
• MCMC using Hamiltonian dynamics
• MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
• Modern Theory of Nuclear Forces
• Rank-normalization, folding, and localization: An improved $\hat{R}$ for assessing convergence of MCMC
• Chiral effective field theory and nuclear forces
• Effective Field Theory for Few-Nucleon Systems
• Nuclear effective field theory: status and perspectives
• Semilocal momentum-space regularized chiral two-nucleon potentials up to fifth order
• Uncertainties of Predictions from Parton Distribution Functions I: the Lagrange Multiplier Method
• Quantifying truncation errors in effective field theory
• Roy-Steiner-equation analysis of pion-nucleon scattering
• Quantifying Correlated Truncation Errors in Effective Field Theory
• Bayesian truncation errors in chiral effective field theory: nucleon-nucleon observables
• Uncertainty analysis and order-by-order optimization of chiral nuclear interactions
• Coarse grained Potential analysis of neutron-proton and proton-proton scattering below pion production threshold
• Global sensitivity analysis of bulk properties of an atomic nucleus
• Bayesian parameter estimation for effective field theories
• Reconciling threshold and subthreshold expansions for pion-nucleon scattering
• Eigenvector Continuation as an Efficient and Accurate Emulator for Uncertainty Quantification
• Exploring Bayesian parameter estimation for chiral effective field theory using nucleon-nucleon phase shifts
• Rigorous constraints on three-nucleon forces in chiral effective field theory from fast and accurate calculations of few-body observables
• Partial Wave Analysis of Nucleon-Nucleon Scattering below pion production threshold
• 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