Abstract:Effective potentials are an essential ingredient of classical molecular dynamics (MD) simulations. Little is understood of the consequences of representing the complex energy landscape of an atomic configuration by an effective potential or force field containing considerably fewer parameters. The probabilistic potential ensemble method has been implemented in the potfit force matching code. This introduces uncertainty quantification into the interatomic potential generation process. Uncertainties in the effective potential are propagated through MD to obtain uncertainties in quantities of interest (QoI), which are a measure of the confidence in the model predictions. We demonstrate the technique using three potentials for nickel: two simple pair potentials, Lennard-Jones and Morse, and a local density dependent embedded atom method potential. A potential ensemble fit to density functional theory (DFT) reference data is constructed for each potential to calculate the uncertainties in lattice constants, elastic constants and thermal expansion. We quantitatively illustrate the cases of poor model selection and fit, highlighted by the uncertainties in the quantities calculated. This shows that our method can capture the effects of the error incurred in QoI resulting from the potential generation process without resorting to comparison with experiment or DFT, which is an essential part to assess the predictive power of MD simulations. Our understanding of the physics underlying material properties relies on verification from computational models of materials and molecules. Such materials simulations also allow us to predict new properties and structures which can then be reproduced experimentally. In order to facilitate the modelling, interatomic potentials have long been used to circumvent the limitations in timescale and simulation size of costly first principles calculations by specifying the energy dependence on the atomic positions. This functional representation is a key constituent of molecular dynamics (MD) simulations, where the quality of the output relies predominantly on the potential employed. Currently, the systematic error incurred in using an interatomic potential is generally unknown, as is the resulting effect on quantities of interest (QoI) it is used to predict, therein forming the motivation for this work. The general idea behind interatomic potentials is that the energy of a collection of atoms can be represented by an explicit functional form or model, dependent on the atomic separation. These analytic potential functions encode the physics into the system and contain a limited number of additional parameters which are adjusted to reproduce desired quantities. The first potentials used intuitive functional forms, fitted to experimental data; however new potentials are frequently fitted to ab-initio data such as atomic forces, energies and stresses. We use the potfit open source force matching code to fit interatomic potentials to density functional theory (DFT) data, and subsequently quantify the uncertainty incurred in simulations using the potential, in lieu of a first principles approach. There has been a significant recent effort in the literature to quantify uncertainty in this multiscale modelling process, with Bayesian frameworks proposed for a variety of interatomic models and force fields [ 1 – 3 ]. There has also been work toward the quantification of uncertainty due to the potential fitting reference set [ 4 ], as well as a proposed framework to efficiently propagate parameter uncertainties to MD outputs [ 5 ]. More specifically, quantification of parameter uncertainty for single potentials has been undertaken in a handful of cases [ 6 – 9 ]. We add to the growing body of uncertainty quantification (UQ) work in potential development and application by providing an open source implementation of the framework in [ 1 ] for use in future potential development projects. We have implemented a new module in potfit which adds a potential ensemble functionality to the potential fitting workflow. For potentials fitted with potfit , the corresponding ensemble informs of the parameter uncertainties, given the associated reference set. Section 2.1 introduces the potfit code, followed by an outline of the potential ensemble method in section 2.2 . Implementation specific details are outlined in section 3 . We then demonstrate in section 4 how an ensemble can be propagated through MD simulations, illustrating the incurred uncertainties in both equilibrium and non-equilibrium QoI: the equil -Abstract Truncated-

Potfit: effective potentials from ab-initio data

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

Transferability of Zr-Zr interatomic potentials

HessFit: A Toolkit to Derive Automated Force Fields from Quantum Mechanical Information

Efficient Calculation of the Lattice Thermal Conductivity by Atomistic Simulations with Ab Initio Accuracy

Searching for Globally Optimal Functional Forms for Inter-Atomic Potentials Using Parallel Tempering and Genetic Programming

Atomic permutationally invariant polynomials for fitting molecular force fields

QDπ: A Quantum Deep Potential Interaction Model for Drug Discovery

Machine Learning Interatomic Potentials for Amorphous Zeolitic Imidazolate Frameworks

Simulation of crack propagation in alumina with ab-initio based polarizable force field

A Hessian-Based Assessment of Atomic Forces for Training Machine Learning Interatomic Potentials

ACEpotentials.jl: A Julia Implementation of the Atomic Cluster Expansion

Gaussian Approximation Potentials: theory, software implementation and application examples

Efficient Approximations of Complete Interatomic Potentials for Crystal Property Prediction

Interatomic potentials of the binary transition metal systems and some applications in materials physics

Adaptive-precision potentials for large-scale atomistic simulations

GAFit: A general-purpose, user-friendly program for fitting potential energy surfaces

Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set

Hyperactive Learning (HAL) for Data-Driven Interatomic Potentials

Machine learned Force-Fields for an ab-initio Quality Description of Metal-Organic Frameworks