TY - JOUR
T1 - Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional
AU - Piaggi, Pablo M.
AU - Panagiotopoulos, Athanassios Z.
AU - Debenedetti, Pablo G.
AU - Car, Roberto
N1 - Publisher Copyright:
©
PY - 2021/5/11
Y1 - 2021/5/11
N2 - Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.
AB - Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.
UR - http://www.scopus.com/inward/record.url?scp=85101440719&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85101440719&partnerID=8YFLogxK
U2 - 10.1021/acs.jctc.1c00041
DO - 10.1021/acs.jctc.1c00041
M3 - Article
C2 - 33835819
AN - SCOPUS:85101440719
SN - 1549-9618
VL - 17
SP - 3065
EP - 3077
JO - Journal of Chemical Theory and Computation
JF - Journal of Chemical Theory and Computation
IS - 5
ER -