Bayesian, frequentist, and information geometric approaches to parametric uncertainty quantification of classical empirical interatomic potentials

Abstract

In this paper, we consider the problem of quantifying parametric uncertainty in classical empirical interatomic potentials (IPs) using both Bayesian (Markov Chain Monte Carlo) and frequentist (profile likelihood) methods. We interface these tools with the Open Knowledgebase of Interatomic Models and study three models based on the Lennard-Jones, Morse, and Stillinger--Weber potentials. We confirm that IPs are typically sloppy, i.e., insensitive to coordinated changes in some parameter combinations. Because the inverse problem in such models is ill-conditioned, parameters are unidentifiable. This presents challenges for traditional statistical methods, as we demonstrate and interpret within both Bayesian and frequentist frameworks. We use information geometry to illuminate the underlying cause of this phenomenon and show that IPs have global properties similar to those of sloppy models from fields such as systems biology, power systems, and critical phenomena. IPs correspond to bounded manifolds with a hierarchy of widths, leading to low effective dimensionality in the model. We show how information geometry can motivate new, natural parameterizations that improve the stability and interpretation of uncertainty quantification analysis and further suggest simplified, less-sloppy models.

References in corpus (4)

• The sloppy model universality class and the Vandermonde matrix
• Force-matched embedded-atom method potential for niobium
• Error Analysis in Nuclear Density Functional Theory
• KLIFF: A framework to develop physics-based and machine learning interatomic potentials

Cited by in corpus (6)

• KLIFF: A framework to develop physics-based and machine learning interatomic potentials
• Data efficiency and extrapolation trends in neural network interatomic potentials
• Uncertainty Quantification in Atomistic Simulations of Silicon using Interatomic Potentials
• Annealing for prediction of grand canonical crystal structures: Efficient implementation of n-body atomic interactions
• An information-matching approach to optimal experimental design and active learning
• Inverse design of bespoke interatomic potentials via active learning by information-matching