Towards shell model interactions with credible uncertainties

Abstract

Background: The nuclear shell model is a powerful framework for predicting nuclear structure observables, but relies on interaction matrix elements fit to experimental data as its inputs. Extending the shell model's applicability, particularly toward dripline nuclei, requires efficient fitting methods and credible uncertainty quantification. Traditional approaches face computational challenges and may underestimate uncertainties. Purpose: We develop and test a framework combining eigenvector continuation and Markov chain Monte Carlo to efficiently fit shell model interaction matrix elements and quantify their uncertainties. Methods: Eigenvector continuation is used to emulate shell model calculations, reducing computational costs. The emulator enables Markov chain Monte Carlo sampling to optimize interaction matrix elements and rigorously assess parametric uncertainties. The framework is benchmarked using the USDB interaction in the sd shell. Results: The emulator reproduces the USDB interaction with negligible error, validating its use in shell model fitting applications. However, we find that to obtain credible predictive intervals, the model defect of the shell model itself, rather than experimental or emulator error, must be taken into account in order to obtain credible uncertainties. Conclusions: The proposed framework provides an efficient and rigorous approach for fitting shell model interactions and quantifying uncertainties. Further, the normality assumption used in the past appears sufficient to describe the distribution of interaction matrix elements. However, it is crucial to account for model correlations to avoid underestimating uncertainties.

Supplemental material provided with Physical Review publication. 16 pages, 9 figures

References in corpus (18)

• emcee: The MCMC Hammer
• The Shell Model as Unified View of Nuclear Structure
• Effective interaction for pf-shell nuclei
• Non-Empirical Interactions for the Nuclear Shell Model: An Update
• The merging of the islands of inversion at N=20 and N=28
• Eigenvector continuation with subspace learning
• Eigenvector Continuation as an Efficient and Accurate Emulator for Uncertainty Quantification
• Ab initio multi-shell valence-space Hamiltonians and the island of inversion
• The Neutron-Rich Edge of the Nuclear Landscape. Experiment and Theory
• Ab initio effective interactions for sd-shell valence nucleons
• Shell-model coupled-cluster method for open-shell nuclei
• Evolution of the N=20 and 28 Shell Gaps and 2-particle-2-hole states in the FSU Interaction
• Perturbative Approach to Effective Shell-Model Hamiltonians and Operators
• Uncertainty-quantified phenomenological optical potentials for single-nucleon scattering
• Constructing approximate shell-model wavefunctions by eigenvector continuation
• Uncertainty quantification of an empirical shell-model interaction using principal component analysis
• Uncertainty quantification in nuclear shell model
• Improving the predictive power of empirical shell-model Hamiltonians