skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Model reduction methods for nuclear emulators
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
Award ID(s):
2209442 2004601 1913069
PAR ID:
10421096
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Journal of Physics G: Nuclear and Particle Physics
Volume:
49
Issue:
10
ISSN:
0954-3899
Page Range / eLocation ID:
102001
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. 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
  4. The sorption properties of [Zr 6 O 4 (OH) 4 (NH 3 + -BDC) 6 ]Cl 6 · x H 2 O ( MOR-1 ) and H 16 [Zr 6 O 16 (H 2 PATP) 4 ]Cl 8 · x H 2 O ( MOR-2 ) towards ReO 4 − and TcO 4 − were studied in detail. Both MOR-1 and MOR-2 are very effective sorbents for ReO 4 − and TcO 4 − anions, with MOR-2 showing the highest sorption capacity (up to 4.1 ± 0.4 mmol g −1 ) among the known metal organic materials. Importantly, the exceptional sorption capacity of MOR-2 is retained even under conditions simulating acidic nuclear waste. In addition, MOR-1 and MOR-2 exhibit selective luminescence ReO 4 − sensing properties, demonstrated for the first time for MOF materials. 
    more » « less
  5. Abstract Climate emulators are a powerful instrument for climate modeling, especially in terms of reducing the computational load for simulating spatiotemporal processes associated with climate systems. The most important type of emulators are statistical emulators trained on the output of an ensemble of simulations from various climate models. However, such emulators oftentimes fail to capture the “physics” of a system that can be detrimental for unveiling critical processes that lead to climate tipping points. Historically, statistical mechanics emerged as a tool to resolve the constraints on physics using statistics. We discuss how climate emulators rooted in statistical mechanics and machine learning can give rise to new climate models that are more reliable and require less observational and computational resources. Our goal is to stimulate discussion on how statistical climate emulators can further be improved with the help of statistical mechanics which, in turn, may reignite the interest of statistical community in statistical mechanics of complex systems. 
    more » « less