Extending OpenKIM with an Uncertainty Quantification Toolkit for Molecular Modeling

Abstract

Atomistic simulations are an important tool in materials modeling. Interatomic potentials (IPs) are at the heart of such molecular models, and the accuracy of a model's predictions depends strongly on the choice of IP. Uncertainty quantification (UQ) is an emerging tool for assessing the reliability of atomistic simulations. The Open Knowledgebase of Interatomic Models (OpenKIM) is a cyberinfrastructure project whose goal is to collect and standardize the study of IPs to enable transparent, reproducible research. Part of the OpenKIM framework is the Python package, KIM-based Learning-Integrated Fitting Framework (KLIFF), that provides tools for fitting parameters in an IP to data. This paper introduces a UQ toolbox extension to KLIFF. We focus on two sources of uncertainty: variations in parameters and inadequacy of the functional form of the IP. Our implementation uses parallel-tempered Markov chain Monte Carlo (PTMCMC), adjusting the sampling temperature to estimate the uncertainty due to the functional form of the IP. We demonstrate on a Stillinger--Weber potential that makes predictions for the atomic energies and forces for silicon in a diamond configuration. Finally, we highlight some potential subtleties in applying and using these tools with recommendations for practitioners and IP developers.

References in corpus (19)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• emcee: The MCMC Hammer
• Interatomic potentials from first-principles calculations: the force-matching method
• Parallel Tempering: Theory, Applications, and New Perspectives
• Interatomic Potential for Silicon Defects and Disordered Phases
• Dynamic temperature selection for parallel-tempering in Markov chain Monte Carlo simulations
• Parameter Space Compression Underlies Emergent Theories and Predictive Models
• Potfit: effective potentials from ab-initio data
• Bayesian Error Estimation in Density Functional Theory
• Why are nonlinear fits so challenging?
• The sloppy model universality class and the Vandermonde matrix
• Force-matched embedded-atom method potential for niobium
• Improvements to the Levenberg-Marquardt algorithm for nonlinear least-squares minimization
• The limitations of model-based experimental design and parameter estimation in sloppy systems
• Bridging Mechanistic and Phenomenological Models of Complex Biological Systems
• The parameters uncertainty inflation fallacy
• KLIFF: A framework to develop physics-based and machine learning interatomic potentials
• Uncertainty quantification for classical effective potentials: an extension to potfit
• A KIM-compliant potfit for fitting sloppy interatomic potentials: Application to the EDIP model for silicon