skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

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

Authors:
ORCiD logo;
Publication Date:
Sponsoring Org.:
USDOE
OSTI Identifier:
1396717
Grant/Contract Number:
B615461
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 334; Journal Issue: C; Related Information: CHORUS Timestamp: 2017-10-04 15:15:56; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English

Citation Formats

Reeve, Samuel Temple, and Strachan, Alejandro. Error correction in multi-fidelity molecular dynamics simulations using functional uncertainty quantification. United States: N. p., 2017. Web. doi:10.1016/j.jcp.2016.12.039.
Reeve, Samuel Temple, & Strachan, Alejandro. Error correction in multi-fidelity molecular dynamics simulations using functional uncertainty quantification. United States. doi:10.1016/j.jcp.2016.12.039.
Reeve, Samuel Temple, and Strachan, Alejandro. Sat . "Error correction in multi-fidelity molecular dynamics simulations using functional uncertainty quantification". United States. doi:10.1016/j.jcp.2016.12.039.
@article{osti_1396717,
title = {Error correction in multi-fidelity molecular dynamics simulations using functional uncertainty quantification},
author = {Reeve, Samuel Temple and Strachan, Alejandro},
abstractNote = {},
doi = {10.1016/j.jcp.2016.12.039},
journal = {Journal of Computational Physics},
number = C,
volume = 334,
place = {United States},
year = {Sat Apr 01 00:00:00 EDT 2017},
month = {Sat Apr 01 00:00:00 EDT 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/j.jcp.2016.12.039

Citation Metrics:
Cited by: 2works
Citation information provided by
Web of Science

Save / Share:
  • 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 undermore » 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.« less
  • In this study we investigate the free energy barrier associated with the dissociation of strong acids, XH (HBr, HCl and HNO3) deprotonation, and subsequent formation of ionpairs, X–___H3O+ in the vicinity of the air-water interface. We will show that the free energy for acid dissociation for HCl and HNO3 show significant differences at the air-water than under bulk solvation conditions producing a picture where at the interface associated molecular species can be stable. For the strongest acid we consider, HBr the more traditional picture of acids is preserved in the vicinity of the air-water interface. Our results have implications formore » our understanding of acids, and their surface tensions at the air-water interface.« less
  • Vibrational infrared (IR) spectra of gas-phase O–H⋅⋅⋅O methanol clusters up to pentamer are simulated using self-consistent-charge density functional tight-binding method using two distinct methodologies: standard normal mode analysis and Fourier transform of the dipole time-correlation function. The twofold simulations aim at the direct critical assignment of the C–H stretching region of the recently recorded experimental spectra [H.-L. Han, C. Camacho, H. A. Witek, and Y.-P. Lee, J. Chem. Phys. 134, 144309 (2011)]. Both approaches confirm the previous assignment (ibid.) of the C–H stretching bands based on the B3LYP/ANO1 harmonic frequencies, showing that ν{sub 3}, ν{sub 9}, and ν{sub 2} C–Hmore » stretching modes of the proton-accepting (PA) and proton-donating (PD) methanol monomers experience only small splittings upon the cluster formation. This finding is in sharp discord with the assignment based on anharmonic B3LYP/VPT2/ANO1 vibrational frequencies (ibid.), suggesting that some procedural faults, likely related to the breakdown of the perturbational vibrational treatment, led the anharmonic calculations astray. The IR spectra based on the Fourier transform of the dipole time-correlation function include new, previously unaccounted for physical factors such as non-zero temperature of the system and large amplitude motions of the clusters. The elevation of temperature results in a considerable non-homogeneous broadening of the observed IR signals, while the presence of large-amplitude motions (methyl group rotations and PA-PD flipping), somewhat surprisingly, does not introduce any new features in the spectrum.« less
  • Abstract not provided.