Error correction in multi-fidelity molecular dynamics simulations using functional uncertainty quantification

We use functional, Fréchet, derivatives to quantify how thermodynamic outputs of a molecular dynamics (MD) simulation depend on the potential used to compute atomic interactions. Our approach quantifies the sensitivity of the quantities of interest with respect to the input functions as opposed to i...

Full description

Saved in:
Bibliographic Details
Published inJournal of computational physics Vol. 334
Main Authors Reeve, Samuel Temple, Strachan, Alejandro
Format Journal Article
LanguageEnglish
Published United States 01.04.2017
Subjects
CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS
COMPUTERIZED SIMULATION
CORRECTIONS
CRYSTALS
ERRORS
FORECASTING
FREE ENERGY
INTERACTIONS
LIQUIDS
MOLECULAR DYNAMICS METHOD
PHASE SPACE
POTENTIAL ENERGY
QUANTUM MECHANICS
SENSITIVITY
THERMODYNAMICS
TWO-BODY PROBLEM
Online AccessGet full text

Cover

More Information
Summary:We use functional, Fréchet, derivatives to quantify how thermodynamic outputs of a molecular dynamics (MD) simulation depend on the potential used to compute atomic interactions. Our approach quantifies the sensitivity of the quantities of interest with respect to the input functions as opposed to its parameters as is done in typical uncertainty quantification methods. We show that the functional sensitivity of the average potential energy and pressure in isothermal, isochoric MD simulations using Lennard–Jones two-body interactions can be used to accurately predict those properties for other interatomic potentials (with different functional forms) without re-running the simulations. This is demonstrated under three different thermodynamic conditions, namely a crystal at room temperature, a liquid at ambient pressure, and a high pressure liquid. The method provides accurate predictions as long as the change in potential can be reasonably described to first order and does not significantly affect the region in phase space explored by the simulation. The functional uncertainty quantification approach can be used to estimate the uncertainties associated with constitutive models used in the simulation and to correct predictions if a more accurate representation becomes available.
ISSN:0021-9991
1090-2716
DOI:10.1016/J.JCP.2016.12.039