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: Bayes goes fast: Uncertainty quantification for a covariant energy density functional emulated by the reduced basis method
A covariant energy density functional is calibrated using a principled Bayesian statistical framework informed by experimental binding energies and charge radii of several magic and semi-magic nuclei. The Bayesian sampling required for the calibration is enabled by the emulation of the high-fidelity model through the implementation of a reduced basis method (RBM)—a set of dimensionality reduction techniques that can speed up demanding calculations involving partial differential equations by several orders of magnitude. The RBM emulator we build—using only 100 evaluations of the high-fidelity model—is able to accurately reproduce the model calculations in tens of milliseconds on a personal computer, an increase in speed of nearly a factor of 3,300 when compared to the original solver. Besides the analysis of the posterior distribution of parameters, we present model calculations for masses and radii with properly estimated uncertainties. We also analyze the model correlation between the slope of the symmetry energy L and the neutron skin of 48 Ca and 208 Pb. The straightforward implementation and outstanding performance of the RBM makes it an ideal tool for assisting the nuclear theory community in providing reliable estimates with properly quantified uncertainties of physical observables. Such uncertainty quantification tools will become essential given the expected abundance of data from the recently inaugurated and future experimental and observational facilities.  more » « less
Award ID(s):
2004601
PAR ID:
10433202
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Frontiers in Physics
Volume:
10
ISSN:
2296-424X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We introduce a model of hadronization based on invertible neural networks that faithfully reproduces a simplified version of the Lund string model for meson hadronization. Additionally, we introduce a new training method for normalizing flows, termed MAGIC, that improves the agreement between simulated and experimental distributions of high-level (macroscopic) observables by adjusting single-emission (microscopic) dynamics. Our results constitute an important step toward realizing a machine-learning based model of hadronization that utilizes experimental data during training. Finally, we demonstrate how a Bayesian extension to this normalizing-flow architecture can be used to provide analysis of statistical and modeling uncertainties on the generated observable distributions. 
    more » « less
  2. A modern power system is characterized by an increasing penetration of wind power, which results in large uncertainties in its states. These uncertainties must be quantified properly; otherwise, the system security may be threatened. Facing this challenge, we propose a cost-effective, data-driven approach for the probabilistic load-margin assessment problem. Using actual wind data, a kernel-density-estimator is applied to infer the nonparametric wind speed distributions, which are further merged into the framework of a vine copula. The latter enables us to simulate complex multivariate and highly dependent model inputs with a variety of bivariate copulae that precisely represent the tail dependence in the correlated samples. Furthermore, to reduce the prohibitive computational time of the traditional Monte-Carlo simulations processing a large amount of samples, we propose to use a nonparametric, Gaussian-process-emulator-based reduced-order model to replace the original complicated continuation power-flow model through a Bayesian learning framework. To accelerate the convergence rate of this Bayesian algorithm, a truncated polynomial chaos surrogate is developed, which serves as a highly efficient, parametric Bayesian prior. This emulator allows us to execute the time-consuming continuation power-flow solver at the sampled values with a negligible computational cost. Simulation results reveal the impressive performances of the proposed method. 
    more » « less
  3. 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
  4. e. With recent advances in online sensing technology and high-performance computing, structural health monitoring (SHM) has begun to emerge as an automated approach to the real-time conditional monitoring of civil infrastructure. Ideal SHM strategies detect and characterize damage by leveraging measured response data to update physics-based finite element models (FEMs). When monitoring composite structures, such as reinforced concrete (RC) bridges, the reliability of FEM based SHM is adversely affected by material, boundary, geometric, and other model uncertainties. Civil engineering researchers have adapted popular artificial intelligence (AI) techniques to overcome these limitations, as AI has an innate ability to solve complex and ill-defined problems by leveraging advanced machine learning techniques to rapidly analyze experimental data. In this vein, this study employs a novel Bayesian estimation technique to update a coupled vehicle-bridge FEM for the purposes of SHM. Unlike existing AI based techniques, the proposed approach makes intelligent use of an embedded FEM model, thus reducing the parameter space while simultaneously guiding the Bayesian model via physics-based principles. To validate the method, bridge response data is generated from the vehicle-bridge FEM given a set of “true” parameters and the bias and standard deviation of the parameter estimates are analyzed. Additionally, the mean parameter estimates are used to solve the FEM model and the results are compared against the results obtained for “true” parameter values. A sensitivity study is also conducted to demonstrate methods for properly formulating model spaces to improve the Bayesian estimation routine. The study concludes with a discussion highlighting factors that need to be considered when leveraging experimental data to update FEMs of concrete structures using AI techniques. 
    more » « less
  5. Ferromagnetic resonance (FMR) is a broadly used dynamical measurement used to characterize a wide range of magnetic materials. Applied research and development on magnetic thin film materials is growing rapidly alongside a growing commercial appetite for magnetic memory and computing technologies. The ability to execute high-quality, fast FMR surveys of magnetic thin films is needed to meet the demanding throughput associated with rapid materials exploration and quality control. Here, we implement optimal Bayesian experimental design software developed by [McMichael et al. J. Res. Natl. Inst. Stand. Technol. 126, 126002 (2021)] in a vector network analyzer-FMR setup to demonstrate an unexplored opportunity to accelerate FMR measurements. A systematic comparison is made between the optimal Bayesian measurement and the conventional measurement. Reduced uncertainties in the linewidth and resonance frequency ranging from 40% to 60% are achieved with the Bayesian implementation for the same measurement duration. In practical terms, this approach reaches a target uncertainty of ±5 MHz for the linewidth and ±1 MHz for the resonance frequency in 2.5× less time than the conventional approach. As the optimal Bayesian approach only decreases random errors, we evaluate how large systematic errors may limit the full advantage of the optimal Bayesian approach. This approach can be used to deliver gains in measurement speed by a factor of 3 or more and as a software add-on has the flexibility to be added on to any FMR measurement system to accelerate materials discovery and quality control measurements, alike. 
    more » « less