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

Title: Communication: Improved ab initio molecular dynamics by minimally biasing with experimental data

Authors:
 [1]; ORCiD logo [2]; ORCiD logo [3]; ORCiD logo [3]
  1. Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 S Ellis Ave., Chicago, Illinois 60637, USA; Department of Chemical Engineering, University of Rochester, Rochester, New York 14627, USA
  2. Leadership Computing Facility, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, Illinois 60439, USA
  3. Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 S Ellis Ave., Chicago, Illinois 60637, USA
Publication Date:
Research Org.:
Argonne National Lab. (ANL), Argonne, IL (United States)
Sponsoring Org.:
US Department of the Navy, Office of Naval Research (ONR); Argonne National Laboratory - Argonne Leadership Computing Facility; University of Chicago - Materials Research Science & Engineering Center (MRSEC)
OSTI Identifier:
1429896
DOE Contract Number:
AC02-06CH11357
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Chemical Physics; Journal Volume: 146; Journal Issue: 4
Country of Publication:
United States
Language:
English

Citation Formats

White, Andrew D., Knight, Chris, Hocky, Glen M., and Voth, Gregory A. Communication: Improved ab initio molecular dynamics by minimally biasing with experimental data. United States: N. p., 2017. Web. doi:10.1063/1.4974837.
White, Andrew D., Knight, Chris, Hocky, Glen M., & Voth, Gregory A. Communication: Improved ab initio molecular dynamics by minimally biasing with experimental data. United States. doi:10.1063/1.4974837.
White, Andrew D., Knight, Chris, Hocky, Glen M., and Voth, Gregory A. Sat . "Communication: Improved ab initio molecular dynamics by minimally biasing with experimental data". United States. doi:10.1063/1.4974837.
@article{osti_1429896,
title = {Communication: Improved ab initio molecular dynamics by minimally biasing with experimental data},
author = {White, Andrew D. and Knight, Chris and Hocky, Glen M. and Voth, Gregory A.},
abstractNote = {},
doi = {10.1063/1.4974837},
journal = {Journal of Chemical Physics},
number = 4,
volume = 146,
place = {United States},
year = {Sat Jan 28 00:00:00 EST 2017},
month = {Sat Jan 28 00:00:00 EST 2017}
}
  • Cited by 1
  • Here, the interface between a solid and a complex multi-component liquid forms a unique reaction environment whose structure and composition can significantly deviate from either bulk or liquid phase and is poorly understood due the innate difficulty to obtain molecular level information. Feldspar minerals, as typified by the Ca-end member Anorthite, serve as prototypical model systems to assess the reactivity and ion mobility at solid/water-bearing supercritical fluid (WBSF) interfaces due to recent X-ray based measurements that provide information on water-film formation, and cation vacancies at these surfaces. Using density functional theory based molecular dynamics, which allows the evaluation of reactivitymore » and condensed phase dynamics on equal footing, we report on the structure and dynamics of water nucleation and surface aggregation, carbonation and Ca mobilization under geologic carbon sequestration scenarios (T = 323 K and P = 90 bar). We find that water has a strong enthalpic preference for aggregation on a Ca-rich, O-terminated anorthite (001) surface, but entropy strongly hinders the film formation at very low water concentrations. Carbonation reactions readily occur at electron-rich terminal Oxygen sites adjacent to cation vacancies, when in contact with supercritical CO 2. Cation vacancies of this type can form readily in the presence of a water layer that allows for facile and enthalpicly favorable Ca 2+ extraction and solvation. Apart from providing unprecedented molecular level detail of a complex three component (mineral, water and scCO 2) system), this work highlights the ability of modern capabilities of AIMD methods to begin to qualitatively and quantitatively address structure and reactivity at solid-liquid interfaces of high chemical complexity. This work was supported by the US Department of Energy, Office of Fossil Energy (M.-S. L., B. P. M. and V.-A. G.) and the Office of Basic Energy Science, Division of Chemical Sciences, Geosciences and Biosciences (R.R.), and performed at the Pacific Northwest National Laboratory (PNNL). PNNL is a multi-program national laboratory operated for DOE by Battelle. Computational resources were provided by PNNL’s Platform for Institutional Computing (PIC), the W. R. Wiley Environmental Molecular Science Laboratory (EMSL), a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research located at PNNL and the National Energy Research Scientific Computing Center (NERSC) at Lawrence Berkeley National Laboratory.« less
  • Ab initio molecular dynamics (AIMD) simulations of the hydration shells surrounding the Zn2+ ion are reported for temperatures near 300oC. Simulations using a combined ab initio and classical molecular dynamics (AIMD/MM) approach are also carried out. Both simulations are done with 64 solvating water molecules (~15 ps). The hydration structure predicted from both simulations is found to agree very well with known results from X-ray data. The 1st hydration shell contains six water molecules in an octahedral structure with the hydrogen atoms oriented away from the Zn2+ ion. The six waters in the 1st shell are located at an averagemore » distance of 2.44Å. A 2nd hydration shell is observed at 4.59Å. Beyond these shells, the bonding pattern substantially returns to the tetrahedral structure of bulk water. No exchanges are seen between the 1st and 2nd hydrations shells, however many water transfers between the 2nd and outer hydrations shells are observed to occur on a picosecond (ps) time scale via dissociative and associative mechanisms. In general, it is found that the AIMD and AIMD/MM simulations give nearly identical results for structural parameters, EXAFS spectra, and exchange dynamics. These results suggest that AIMD/MM can be used to extend the particle scale and time scale of AIMD simulations of highly charged ions in solution.« less