Phase equilibrium of water with hexagonal and cubic ice using the SCAN functional

Abstract

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 towards 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 ice 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.

20 pages, 9 figures

References in corpus (14)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Canonical sampling through velocity-rescaling
• Advanced capabilities for materials modelling with Quantum ESPRESSO
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• Libxc: a library of exchange and correlation functionals for density functional theory
• Competing quantum effects in the dynamics of a flexible water model
• Ab initio theory and modeling of water
• A Variational Approach to Enhanced Sampling and Free Energy Calculations
• On the accuracy of the MB-pol many-body potential for water: Interaction energies, vibrational frequencies, and classical thermodynamic and dynamical properties from clusters to liquid water and ice
• 86 PFLOPS Deep Potential Molecular Dynamics simulation of 100 million atoms with ab initio accuracy
• Homogeneous ice nucleation evaluated for several water models
• A Unified Approach to Enhanced Sampling
• Homogeneous nucleation of ice
• Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations

Cited by in corpus (18)

• Deep Potentials for Materials Science
• Homogeneous ice nucleation in an ab initio machine learning model of water
• Liquid-liquid transition in water from first principles
• Heat transport in liquid water from first-principles and deep-neural-network simulations
• Viscosity in water from first-principles and deep-neural-network simulations
• A first-principles machine-learning force field for heterogeneous ice nucleation on microcline feldspar
• ArcaNN: automated enhanced sampling generation of training sets for chemically reactive machine learning interatomic potentials
• DMC-ICE13: ambient and high pressure polymorphs of ice from Diffusion Monte Carlo and Density Functional Theory
• A Deep Potential model for liquid-vapor equilibrium and cavitation rates of water
• Recent advances in describing and driving crystal nucleation using machine learning and artificial intelligence
• NEP-MB-pol: A unified machine-learned framework for fast and accurate prediction of water's thermodynamic and transport properties
• Size and Quality of Quantum Mechanical Data Sets for Training Neural Network Force Fields for Liquid Water
• Melting curves of ice polymorphs in the vicinity of the liquid-liquid critical point
• Pressure-Induced Structural and Dielectric Changes in Liquid Water at Room Temperature
• Dynamical Heterogeneity in Supercooled Water and its Spectroscopic Fingerprints
• Quantum Anharmonic Effects in Hydrogen-Bond Symmetrization of High-Pressure Ice
• Kinetics of Stacking Order Evolution During Heterogeneous Ice Formation
• Determination of melting temperature of hexagonal ice using Lee-Yang phase transition theory