Molecular Dynamics Simulations of Liquid Water Using High Quality Meta-GGA Functionals

Overview

Journal Chem Sci

Publisher Royal Society of Chemistry

Specialty Chemistry

Date 2018 Aug 30

PMID 30155200

Citations 22

Authors

Luis Ruiz Pestana

Narbe Mardirossian

Martin Head-Gordon

Teresa Head-Gordon

Affiliations

Soon will be listed here.

Abstract

We have used molecular dynamics (AIMD) to characterize water properties using two meta-generalized gradient approximation (meta-GGA) functionals, M06-L-D3 and B97M-rV, and compared their performance against a standard GGA corrected for dispersion, revPBE-D3, at ambient conditions (298 K, and 1 g cm or 1 atm). Simulations of the equilibrium density, radial distribution functions, self-diffusivity, the infrared spectrum, liquid dipole moments, and characterizations of the hydrogen bond network show that all three functionals have overcome the problem of the early AIMD simulations that erroneously found ambient water to be highly structured, but they differ substantially among themselves in agreement with experiment on this range of water properties. We show directly using water cluster data up through the pentamer that revPBE-D3 benefits from a cancellation of its intrinsic functional error by running classical trajectories, whereas the meta-GGA functionals are demonstrably more accurate and would require the simulation of nuclear quantum effects to realize better agreement with all cluster and condensed phase properties.

Citing Articles

Oxygen Vacancies on Hydrated Anatase (101) Surfaces: Insights from Classical and Ab Initio Molecular Dynamics Simulations.

Grote F, Lyubartsev A Nanomaterials (Basel). 2025; 15(5).

PMID: 40072167 PMC: 11901595. DOI: 10.3390/nano15050364.

A Robust and Versatile QM/MM Interface for Molecular Dynamics in GROMOS.

Poliak P, Bleiziffer P, Pultar F, Riniker S, Oostenbrink C J Comput Chem. 2025; 46(5):e70053.

PMID: 39918182 PMC: 11804165. DOI: 10.1002/jcc.70053.

The Effect of Surface Oxygen Coverage on the Oxygen Evolution Reaction over a CoFeNiCr High-Entropy Alloy.

Yuan G, Pestana L Nanomaterials (Basel). 2024; 14(12).

PMID: 38921934 PMC: 11206725. DOI: 10.3390/nano14121058.

Hydration of -Hydroxyurea from Ab Initio Molecular Dynamics Simulations.

Balicki M, Smiechowski M Molecules. 2024; 29(11).

PMID: 38893311 PMC: 11173572. DOI: 10.3390/molecules29112435.

To Pair or not to Pair? Machine-Learned Explicitly-Correlated Electronic Structure for NaCl in Water.

ONeill N, Shi B, Fong K, Michaelides A, Schran C J Phys Chem Lett. 2024; 15(23):6081-6091.

PMID: 38820256 PMC: 11181334. DOI: 10.1021/acs.jpclett.4c01030.

References

Lee H, Tuckerman M . Dynamical properties of liquid water from ab initio molecular dynamics performed in the complete basis set limit. J Chem Phys. 2007; 126(16):164501. DOI: 10.1063/1.2718521. View

Kuo I, Mundy C . An ab initio molecular dynamics study of the aqueous liquid-vapor interface. Science. 2004; 303(5658):658-60. DOI: 10.1126/science.1092787. View

Head-Gordon T, Johnson M . Tetrahedral structure or chains for liquid water. Proc Natl Acad Sci U S A. 2006; 103(21):7973-7. PMC: 1472414. DOI: 10.1073/pnas.0510593103. View

Perdew , Burke , Ernzerhof . Generalized Gradient Approximation Made Simple. Phys Rev Lett. 1996; 77(18):3865-3868. DOI: 10.1103/PhysRevLett.77.3865. View

Forster-Tonigold K, Gross A . Dispersion corrected RPBE studies of liquid water. J Chem Phys. 2014; 141(6):064501. DOI: 10.1063/1.4892400. View

Skinner L, Huang C, Schlesinger D, Pettersson L, Nilsson A, Benmore C . Benchmark oxygen-oxygen pair-distribution function of ambient water from x-ray diffraction measurements with a wide Q-range. J Chem Phys. 2013; 138(7):074506. DOI: 10.1063/1.4790861. View

Liu H, Wang Y, Bowman J . Quantum calculations of the IR spectrum of liquid water using ab initio and model potential and dipole moment surfaces and comparison with experiment. J Chem Phys. 2015; 142(19):194502. DOI: 10.1063/1.4921045. View

Dahlke E, Truhlar D . Assessment of the pairwise additive approximation and evaluation of many-body terms for water clusters. J Phys Chem B. 2006; 110(22):10595-601. DOI: 10.1021/jp061039e. View

Albaugh A, Demerdash O, Head-Gordon T . An efficient and stable hybrid extended Lagrangian/self-consistent field scheme for solving classical mutual induction. J Chem Phys. 2015; 143(17):174104. DOI: 10.1063/1.4933375. View

10.

Klimes J, Michaelides A . Perspective: Advances and challenges in treating van der Waals dispersion forces in density functional theory. J Chem Phys. 2012; 137(12):120901. DOI: 10.1063/1.4754130. View

11.

Zhang C, Wu J, Galli G, Gygi F . Structural and Vibrational Properties of Liquid Water from van der Waals Density Functionals. J Chem Theory Comput. 2015; 7(10):3054-61. DOI: 10.1021/ct200329e. View

12.

Marshall M, Burns L, Sherrill C . Basis set convergence of the coupled-cluster correction, δ(MP2)(CCSD(T)): best practices for benchmarking non-covalent interactions and the attendant revision of the S22, NBC10, HBC6, and HSG databases. J Chem Phys. 2011; 135(19):194102. DOI: 10.1063/1.3659142. View

13.

Miller T, Manolopoulos D . Quantum diffusion in liquid water from ring polymer molecular dynamics. J Chem Phys. 2005; 123(15):154504. DOI: 10.1063/1.2074967. View

14.

Dion M, Rydberg H, Schroder E, Langreth D, Lundqvist B . van der Waals density functional for general geometries. Phys Rev Lett. 2004; 92(24):246401. DOI: 10.1103/PhysRevLett.92.246401. View

15.

Wang L, Martinez T, Pande V . Building Force Fields: An Automatic, Systematic, and Reproducible Approach. J Phys Chem Lett. 2015; 5(11):1885-91. PMC: 9649520. DOI: 10.1021/jz500737m. View

16.

Goedecker , Teter , HUTTER . Separable dual-space Gaussian pseudopotentials. Phys Rev B Condens Matter. 1996; 54(3):1703-1710. DOI: 10.1103/physrevb.54.1703. View

17.

Yoo S, Zeng X, Xantheas S . On the phase diagram of water with density functional theory potentials: The melting temperature of ice I(h) with the Perdew-Burke-Ernzerhof and Becke-Lee-Yang-Parr functionals. J Chem Phys. 2009; 130(22):221102. DOI: 10.1063/1.3153871. View

18.

Cohen A, Mori-Sanchez P, Yang W . Insights into current limitations of density functional theory. Science. 2008; 321(5890):792-4. DOI: 10.1126/science.1158722. View

19.

Fanourgakis G, Apra E, Xantheas S . High-level ab initio calculations for the four low-lying families of minima of (H2O)20. I. Estimates of MP2/CBS binding energies and comparison with empirical potentials. J Chem Phys. 2004; 121(6):2655-63. DOI: 10.1063/1.1767519. View

20.

Matsumoto M, Saito S, Ohmine I . Molecular dynamics simulation of the ice nucleation and growth process leading to water freezing. Nature. 2002; 416(6879):409-13. DOI: 10.1038/416409a. View