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

Title: Exploring the free energy surface using ab initio molecular dynamics

Abstract

Efficient exploration of the configuration space and identification of metastable structures are challenging from both computational as well as algorithmic perspectives. Here, we extend the recently proposed orderparameter aided temperature accelerated sampling schemes to efficiently and systematically explore free energy surfaces, and search for metastable states and reaction pathways within the framework of density functional theory based molecular dynamics. The sampling method is applied to explore the relevant parts of the configuration space in prototypical materials SiO 2 and Ti to identify the different metastable structures corresponding to different phases in these materials. In addition, we use the string method in collective variables to study the melting pathways in the high pressure cotunnite phase of SiO 2 and the hcp to fcc phase transition in Ti.

Authors:
 [1];  [1];  [1]
  1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC-22)
OSTI Identifier:
1262186
Alternate Identifier(s):
OSTI ID: 1249678
Report Number(s):
LLNL-JRNL-674389
Journal ID: ISSN 0021-9606; JCPSA6
Grant/Contract Number:
AC52-07NA27344
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Journal of Chemical Physics
Additional Journal Information:
Journal Volume: 144; Journal Issue: 16; Journal ID: ISSN 0021-9606
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English
Subject:
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SUPERFLUIDITY; Lagrangian mechanics; Free energy; Melting; Molecular dynamics

Citation Formats

Samanta, Amit, Morales, Miguel A., and Schwegler, Eric. Exploring the free energy surface using ab initio molecular dynamics. United States: N. p., 2016. Web. doi:10.1063/1.4945653.
Samanta, Amit, Morales, Miguel A., & Schwegler, Eric. Exploring the free energy surface using ab initio molecular dynamics. United States. doi:10.1063/1.4945653.
Samanta, Amit, Morales, Miguel A., and Schwegler, Eric. Fri . "Exploring the free energy surface using ab initio molecular dynamics". United States. doi:10.1063/1.4945653. https://www.osti.gov/servlets/purl/1262186.
@article{osti_1262186,
title = {Exploring the free energy surface using ab initio molecular dynamics},
author = {Samanta, Amit and Morales, Miguel A. and Schwegler, Eric},
abstractNote = {Efficient exploration of the configuration space and identification of metastable structures are challenging from both computational as well as algorithmic perspectives. Here, we extend the recently proposed orderparameter aided temperature accelerated sampling schemes to efficiently and systematically explore free energy surfaces, and search for metastable states and reaction pathways within the framework of density functional theory based molecular dynamics. The sampling method is applied to explore the relevant parts of the configuration space in prototypical materials SiO2 and Ti to identify the different metastable structures corresponding to different phases in these materials. In addition, we use the string method in collective variables to study the melting pathways in the high pressure cotunnite phase of SiO2 and the hcp to fcc phase transition in Ti.},
doi = {10.1063/1.4945653},
journal = {Journal of Chemical Physics},
number = 16,
volume = 144,
place = {United States},
year = {Fri Apr 22 00:00:00 EDT 2016},
month = {Fri Apr 22 00:00:00 EDT 2016}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

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

Save / Share:
  • Cited by 3
  • Cited by 10
  • The ab initio/classical free energy perturbation (ABC-FEP) method proposed previously by Wood et al. [J. Chem. Phys. 110, 1329 (1999)] uses classical simulations to calculate solvation free energies within an empirical potential model, then applies free energy perturbation theory to determine the effect of changing the empirical solute-solvent interactions to corresponding interactions calculated from ab initio methods. This approach allows accurate calculation of solvation free energies using an atomistic description of the solvent and solute, with interactions calculated from first principles. Results can be obtained at a feasible computational cost without making use of approximations such as a continuum solventmore » or an empirical cavity formation energy. As such, the method can be used far from ambient conditions, where the empirical parameters needed for approximate theories of solvation may not be available. The sources of error in the ABC-FEP method are the approximations in the ab initio method, the finite sample of configurations, and the classical solvent model. This article explores the accuracy of various approximations used in the ABC-FEP method by comparing to the experimentally well-known free energy of hydration of water at two state points (ambient conditions, and 973.15 K and 600 kg/m3). The TIP4P-FQ model [J. Chem. Phys. 101, 6141 (1994)] is found to be a reliable solvent model for use with this method, even at supercritical conditions. Results depend strongly on the ab initio method used: a gradient-corrected density functional theory is not adequate, but a localized MP2 method yields excellent agreement with experiment. Computational costs are reduced by using a cluster approximation, in which ab initio pair interaction energies are calculated between the solute and up to 60 solvent molecules, while multi-body interactions are calculated with only a small cluster (5 to 12 solvent molecules). Sampling errors for the ab initio contribution to solvation free energies are {+-}2 kJ/mol or less when 50-200 configurations are used. Using the largest clusters and most accurate ab initio methods, ABC-FEP predicts hydration free energies of water at both state points that agree with equations of state, within the sampling error. These results are the first calculation of a free energy of solvation at extreme conditions from a fully atomistic model with ab initio methods. (c) 2000 American Institute of Physics.« less
  • We present an ab initio approach for evaluating a free energy profile along a reaction coordinate by combining logarithmic mean force dynamics (LogMFD) and first-principles molecular dynamics. The mean force, which is the derivative of the free energy with respect to the reaction coordinate, is estimated using density functional theory (DFT) in the present approach, which is expected to provide an accurate free energy profile along the reaction coordinate. We apply this new method, first-principles LogMFD (FP-LogMFD), to a glycine dipeptide molecule and reconstruct one- and two-dimensional free energy profiles in the framework of DFT. The resultant free energy profilemore » is compared with that obtained by the thermodynamic integration method and by the previous LogMFD calculation using an empirical force-field, showing that FP-LogMFD is a promising method to calculate free energy without empirical force-fields.« less