Fast & rigorous predictions for A=6 nuclei with Bayesian posterior sampling
Abstract
We make ab initio predictions for the A = 6 nuclear level scheme based on two- and three-nucleon interactions up to next-to-next-to-leading order in chiral effective field theory (). We utilize eigenvector continuation and Bayesian methods to quantify uncertainties stemming from the many-body method, the truncation, and the low-energy constants of the nuclear interaction. The construction and validation of emulators is made possible via the development of JupiterNCSM -- a new M-scheme no-core shell model code that uses on-the-fly Hamiltonian matrix construction for efficient, single-node computations up to Nmax = 10 for 6Li. We find a slight underbinding of 6He and 6Li, although consistent with experimental data given our theoretical error bars. As a result of incorporating a correlated -truncation errors we find more precise predictions (smaller error bars) for separation energies: Sd(6Li) = 0.89 0.44 MeV, S2n(6He) = 0.20 0.60 MeV, and for the beta decay Q-value: Qβ-(6He) = 3.71 0.65 MeV. We conclude that our error bars can potentially be reduced further by extending the model space used by JupiterNCSM.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.