An official website of the United States government
Abstract The field of model order reduction (MOR) is growing in importance due to its ability to extract the key insights from complex simulations while discarding computationally burdensome and superfluous information. We provide an overview of MOR methods for the creation of fast & accurate emulators of memory- and compute-intensive nuclear systems, focusing on eigen-emulators and variational emulators. As an example, we describe how ‘eigenvector continuation’ is a special case of a much more general and well-studied MOR formalism for parameterized systems. We continue with an introduction to the Ritz and Galerkin projection methods that underpin many such emulators, while pointing to the relevant MOR theory and its successful applications along the way. We believe that this guide will open the door to broader applications in nuclear physics and facilitate communication with practitioners in other fields.
more »
« less
BUQEYE guide to projection-based emulators in nuclear physics
The BUQEYE collaboration (Bayesian Uncertainty Quantification: Errors in Your effective field theory) presents a pedagogical introduction to projection-based, reduced-order emulators for applications in low-energy nuclear physics. The term emulator refers here to a fast surrogate model capable of reliably approximating high-fidelity models. As the general tools employed by these emulators are not yet well-known in the nuclear physics community, we discuss variational and Galerkin projection methods, emphasize the benefits of offline-online decompositions, and explore how these concepts lead to emulators for bound and scattering systems that enable fast and accurate calculations using many different model parameter sets. We also point to future extensions and applications of these emulators for nuclear physics, guided by the mature field of model (order) reduction. All examples discussed here and more are available as interactive, open-source Python code so that practitioners can readily adapt projection-based emulators for their own work.
more »
« less
- PAR ID:
- 10420991
- Date Published:
- Journal Name:
- Frontiers in Physics
- Volume:
- 10
- ISSN:
- 2296-424X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Gravitational-wave observations of binary neutron-star (BNS) mergers have the potential to revolutionize our understanding of the nuclear equation of state (EOS) and the fundamental interactions that determine its properties. However, Bayesian parameter estimation frameworks do not typically sample over microscopic nuclear-physics parameters that determine the EOS. One of the major hurdles in doing so is the computational cost involved in solving the neutron-star structure equations, known as the Tolman–Oppenheimer–Volkoff (TOV) equations. In this paper, we explore approaches to emulating solutions for the TOV equations: multilayer perceptrons (MLPs), Gaussian processes, and a data-driven variant of the reduced basis method (RBM). We implement these emulators for three different parameterizations of the nuclear EOS, each with a different degree of complexity represented by the number of model parameters. We find that our MLP-based emulators are generally more accurate than the other two algorithms, whereas the RBM results in the largest speedup with respect to the full high-fidelity TOV solver. We employ these emulators for a simple parameter inference using a potentially loud BNS observation and show that the posteriors predicted by our emulators are in excellent agreement with those obtained from the full TOV solver.more » « less
-
Penning-trap mass spectrometry in atomic and nuclear physics has become a well-established and reliable tool for the determination of atomic masses. In combination with short-lived radioactive nuclides it was first introduced at ISOLTRAP at the Isotope Mass Separator On-Line facility (ISOLDE) at CERN. Penning traps have found new applications in coupling to other production mechanisms, such as in-flight production and separation systems. The applications in atomic and nuclear physics range from nuclear structure studies and related precision tests of theoretical approaches to description of the strong interaction to tests of the electroweak Standard Model, quantum electrodynamics and neutrino physics, and applications in nuclear astrophysics. The success of Penning-trap mass spectrometry is due to its precision and accuracy, even for low ion intensities (i.e., low production yields), as well as its very fast measurement cycle, enabling access to short-lived isotopes. The current reach in relative mass precision goes beyond δ m/ m=10 −8 , the half-life limit is as low as a few milliseconds, and the sensitivity is on the order of one ion per minute in the trap. We provide a comprehensive overview of the techniques and applications of Penning-trap mass spectrometry in nuclear and atomic physics.more » « less
-
We discuss emulators from the ab initio symmetry-adapted no-core shell-model framework for studying the formation of alpha clustering and collective properties without effective charges. We present a new type of an emulator, one that utilizes the eigenvector continuation technique but is based on the use of symplectic symmetry considerations. This is achieved by using physically relevant degrees of freedom, namely, the symmetry-adapted basis, which exploits the almost perfect symplectic symmetry in nuclei. Specifically, we study excitation energies, point-proton root-mean-square radii, along with electric quadrupole moments and transitions for 6 Li and 12 C. We show that the set of parameterizations of the chiral potential used to train the emulators has no significant effect on predictions of dominant nuclear features, such as shape and the associated symplectic symmetry, along with cluster formation, but slightly varies details that affect collective quadrupole moments, asymptotic normalization coefficients, and alpha partial widths up to a factor of two. This makes these types of emulators important for further constraining the nuclear force for high-precision nuclear structure and reaction observables.more » « less
-
Elastic continuum mechanical models are widely used to compute deformations due to pressure changes in buried cavities, such as magma reservoirs. In general, analytical models are fast but can be inaccurate as they do not correctly satisfy boundary conditions for many geometries, while numerical models are slow and may require specialized expertise and software. To overcome these limitations, we trained supervised machine learning emulators (model surrogates) based on parallel partial Gaussian processes which predict the output of a finite element numerical model with high fidelity but >1,000× greater computational efficiency. The emulators are based on generalized nondimensional forms of governing equations for finite non‐dipping spheroidal cavities in elastic halfspaces. Either cavity volume change or uniform pressure change boundary conditions can be specified, and the models predict both surface displacements and cavity (pore) compressibility. Because of their computational efficiency, using the emulators as numerical model surrogates can greatly accelerate data inversion algorithms such as those employing Bayesian Markov chain Monte Carlo sampling. The emulators also permit a comprehensive evaluation of how displacements and cavity compressibility vary with geometry and material properties, revealing the limitations of analytical models. Our open‐source emulator code can be utilized without finite element software, is suitable for a wide range of cavity geometries and depths, includes an estimate of uncertainties associated with emulation, and can be used to train new emulators for different source geometries.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3389/fphy.2022.1092931
