A Bayesian Approach for Parameter Estimation and Prediction using a Computationally Intensive Model

Abstract

Bayesian methods have been very successful in quantifying uncertainty in physics-based problems in parameter estimation and prediction. In these cases, physical measurements y are modeled as the best fit of a physics-based model $\eta (\theta )$ where $\theta$ denotes the uncertain, best input setting. Hence the statistical model is of the form $y=\eta (\theta )+\epsilon$, where $\epsilon$ accounts for measurement, and possibly other error sources. When non-linearity is present in $\eta (\cdot )$, the resulting posterior distribution for the unknown parameters in the Bayesian formulation is typically complex and non-standard, requiring computationally demanding computational approaches such as Markov chain Monte Carlo (MCMC) to produce multivariate draws from the posterior. While quite generally applicable, MCMC requires thousands, or even millions of evaluations of the physics model $\eta (\cdot )$. This is problematic if the model takes hours or days to evaluate. To overcome this computational bottleneck, we present an approach adapted from Bayesian model calibration. This approach combines output from an ensemble of computational model runs with physical measurements, within a statistical formulation, to carry out inference. A key component of this approach is a statistical response surface, or emulator, estimated from the ensemble of model runs. We demonstrate this approach with a case study in estimating parameters for a density functional theory (DFT) model, using experimental mass/binding energy measurements from a collection of atomic nuclei. We also demonstrate how this approach produces uncertainties in predictions for recent mass measurements obtained at Argonne National Laboratory (ANL).

19 pages, 10 figures, 1 table; Invited paper for the Journal of Physics G: Nuclear and Particle Physics focus section entitled "Enhancing the interaction between nuclear experiment and theory through information and statistics". In press

References in corpus (5)

• Nuclear Energy Density Optimization
• Nuclear energy density optimization: Large deformations
• Cosmic Calibration
• First Results from the CARIBU Facility: Mass Measurements on the r-Process Path
• Error Analysis in Nuclear Density Functional Theory

Cited by in corpus (28)

• Applying Bayesian parameter estimation to relativistic heavy-ion collisions: simultaneous characterization of the initial state and quark-gluon plasma medium
• Machine Learning in Nuclear Physics
• Future of Nuclear Fission Theory
• A data-driven analysis for the temperature and momentum dependence of the heavy quark diffusion coefficient in relativistic heavy-ion collisions
• Uncertainty Quantification for Nuclear Density Functional Theory and Information Content of New Measurements
• Estimating initial state and quark-gluon plasma medium properties using a hybrid model with nucleon substructure calibrated to $p$-Pb and Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV
• Bayesian parameter estimation for effective field theories
• Quantifying properties of hot and dense QCD matter through systematic model-to-data comparison
• Exploring QCD matter in extreme conditions with Machine Learning
• Constraints on rapidity-dependent initial conditions from charged particle pseudorapidity densities and two-particle correlations
• Error Analysis in Nuclear Density Functional Theory
• Propagation of Statistical Uncertainties of Skyrme Mass Models to Simulations of $r$-Process Nucleosynthesis
• Statistical aspects of nuclear mass models
• Spectral likelihood expansions for Bayesian inference
• Determination of Quark-Gluon-Plasma Parameters from a Global Bayesian Analysis
• Estimating parameter uncertainty in binding-energy models by the frequency-domain bootstrap
• Wave-function-based emulation for nucleon-nucleon scattering in momentum space
• Bayesian tomography with prior-knowledge-based parametrization and surrogate modeling
• Examination of nucleon distribution with Bayesian imaging for isobar collisions
• Optimization and Supervised Machine Learning Methods for Fitting Numerical Physics Models without Derivatives
• Building Surrogate Models of Nuclear Density Functional Theory with Gaussian Processesand Autoencoders
• Uncertainty propagation within the UNEDF models
• Calibration of Energy Density Functionals with Deformed Nuclei
• On model emulation and closure tests for 3+1D relativistic heavy-ion collisions
• A Fast and Calibrated Computer Model Emulator: An Empirical Bayes Approach
• Constructing a simulation surrogate with partially observed output
• Shell quenching in nuclear charge radii based on Monte Carlo dropout Bayesian neural network
• Multipole responses in fissioning nuclei and their uncertainties