TY - JOUR
T1 - Why Is MP2-Water ''Cooler'' and ''Denser'' than DFT-Water?
AU - Willow, Soohaeng Yoo
AU - Zeng, Xiao Cheng
AU - Xantheas, Sotiris S.
AU - Kim, Kwang S.
AU - Hirata, So
N1 - Funding Information:
S.Y.W. and S.H. were supported by the National Science Foundation under award CHE-1361586 and by CREST, Japan Science and Technology Agency. X.C.Z. was supported by grants from the NSF (CBET-1066947 and CHE-1306326) and the Nebraska Research Initiative. S.S.X. was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences at Pacific Northwest National Laboratory. Battelle operates the Pacific Northwest National Laboratory for the U.S. Department of Energy. K.S.K was supported by Korean National Research Foundation (National Honor Scientist 2010-0020414) and by Korea Institute of Science and Technology Information (KSC-2013-C3-019).
Publisher Copyright:
© 2016 American Chemical Society.
PY - 2016/2/18
Y1 - 2016/2/18
N2 - Density functional theory (DFT) with a dispersionless generalized gradient approximation (GGA) needs much higher temperature and pressure than the ambient conditions to maintain water in the liquid phase at the correct (1 g/cm3) density during first-principles simulations. Conversely, ab initio second-order many-body perturbation (MP2) calculations of liquid water require lower temperature and pressure than DFT/GGA to keep water liquid. Here we present a unifying explanation of these trends derived from classical water simulations using a polarizable force field with different sets of parameters. We show that the different temperatures and pressures between DFT/GGA and MP2 at which the simulated water displays the experimentally observed liquid structure under the ambient conditions can be largely explained by their differences in polarizability and dispersion interaction, respectively. In DFT/GGA, the polarizability and thus the induced dipole moments and the hydrogen-bond strength are all overestimated. This hinders the rotational motion of molecules and requires a higher temperature for DFT-water to be liquid. MP2 gives a stronger dispersion interaction and thus shorter intermolecular distances than dispersionless DFT/GGA, which is why MP2-water is denser than DFT-water under the same external pressure.
AB - Density functional theory (DFT) with a dispersionless generalized gradient approximation (GGA) needs much higher temperature and pressure than the ambient conditions to maintain water in the liquid phase at the correct (1 g/cm3) density during first-principles simulations. Conversely, ab initio second-order many-body perturbation (MP2) calculations of liquid water require lower temperature and pressure than DFT/GGA to keep water liquid. Here we present a unifying explanation of these trends derived from classical water simulations using a polarizable force field with different sets of parameters. We show that the different temperatures and pressures between DFT/GGA and MP2 at which the simulated water displays the experimentally observed liquid structure under the ambient conditions can be largely explained by their differences in polarizability and dispersion interaction, respectively. In DFT/GGA, the polarizability and thus the induced dipole moments and the hydrogen-bond strength are all overestimated. This hinders the rotational motion of molecules and requires a higher temperature for DFT-water to be liquid. MP2 gives a stronger dispersion interaction and thus shorter intermolecular distances than dispersionless DFT/GGA, which is why MP2-water is denser than DFT-water under the same external pressure.
UR - http://www.scopus.com/inward/record.url?scp=84958956546&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84958956546&partnerID=8YFLogxK
U2 - 10.1021/acs.jpclett.5b02430
DO - 10.1021/acs.jpclett.5b02430
M3 - Article
C2 - 26821830
AN - SCOPUS:84958956546
SN - 1948-7185
VL - 7
SP - 680
EP - 684
JO - Journal of Physical Chemistry Letters
JF - Journal of Physical Chemistry Letters
IS - 4
ER -