Greedy Emulators for Nuclear Two-Body Scattering

Abstract

Applications of reduced basis method emulators are increasing in low-energy nuclear physics because they enable fast and accurate sampling of high-fidelity calculations, enabling robust uncertainty quantification. In this paper, we develop, implement, and test two model-driven emulators based on (Petrov-)Galerkin projection using the prototypical test case of two-body scattering with the Minnesota potential and a more realistic local chiral potential. The high-fidelity scattering equations are solved with the matrix Numerov method, a reformulation of the popular Numerov recurrence relation for solving special second-order differential equations as a linear system of coupled equations. A novel error estimator based on reduced-space residuals is applied to an active learning approach (a greedy algorithm) to choosing training samples ("snapshots") for the emulator and contrasted with a proper orthogonal decomposition (POD) approach. Both approaches allow for computationally efficient offline-online decompositions, but the greedy approach requires much fewer snapshot calculations. These developments set the groundwork for emulating scattering observables based on chiral nucleon-nucleon and three-nucleon interactions and optical models, where computational speed-ups are necessary for Bayesian uncertainty quantification. Our emulators and error estimators are widely applicable to linear systems.

minor changes, close to published version; 24 pages, 11 figures, and 2 tables

References in corpus (17)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Local chiral effective field theory interactions and quantum Monte Carlo applications
• Global sensitivity analysis of bulk properties of an atomic nucleus
• pyMOR - Generic Algorithms and Interfaces for Model Order Reduction
• Rigorous constraints on three-nucleon forces in chiral effective field theory from fast and accurate calculations of few-body observables
• Efficient emulators for scattering using eigenvector continuation
• Model reduction methods for nuclear emulators
• Fast & accurate emulation of two-body scattering observables without wave functions
• Eigenvector Continuation and Projection-Based Emulators
• Power counting in chiral effective field theory and nuclear binding
• Fast emulation of quantum three-body scattering
• Self-learning Emulators and Eigenvector Continuation
• A simple and efficient numerical scheme to integrate non-local potentials
• Emulating \emph{ab initio} computations of infinite nucleonic matter
• Surrogate models for quantum spin systems based on reduced order modeling
• Wave-function-based emulation for nucleon-nucleon scattering in momentum space
• Reduced basis surrogates for quantum spin systems based on tensor networks

Cited by in corpus (6)

• Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description
• Non-Hermitian Quantum Mechanics Approach for Extracting and Emulating Continuum Physics Based on Bound-State-Like Calculations
• Deformation, halo, and bubble structure: A paradigm shift of exotic phenomena in light to medium mass nuclei
• Star Log-extended eMulation: a method for efficient computation of the Tolman-Oppenheimer-Volkoff equations
• Active learning emulators for nuclear two-body scattering in momentum space
• Toward a Unified Understanding of the Dense Matter Equation of State