Enhancing the formation of ionic defects to study the ice Ih/XI transition with molecular dynamics simulations

Abstract

Ice Ih, the common form of ice in the biosphere, contains proton disorder. Its proton-ordered counterpart, ice XI, is thermodynamically stable below 72 K. However, even below this temperature the formation of ice XI is kinetically hindered and experimentally it is obtained by doping ice with KOH. Doping creates ionic defects that promote the migration of protons and the associated change in proton configuration. In this article, we mimic the effect of doping in molecular dynamics simulations using a bias potential that enhances the formation of ionic defects. The recombination of the ions thus formed proceeds through fast migration of the hydroxide and results in the jump of protons along a hydrogen bond loop. This provides a physical and expedite way to change the proton configuration, and to accelerate diffusion in proton configuration space. A key ingredient of this approach is a machine learning potential trained with density functional theory data and capable of modeling molecular dissociation. We exemplify the usefulness of this idea by studying the order-disorder transition using an appropriate order parameter to distinguish the proton environments in ice Ih and XI. We calculate the changes in free energy, enthalpy, and entropy associated with the transition. Our estimated entropy agrees with experiment within the error bars of our calculation.

17 pages, 10 figures

References in corpus (21)

• Escaping free-energy minima
• Maximally-localized generalized Wannier functions for composite energy bands
• PLUMED 2: New feathers for an old bird
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• On representing chemical environments
• Deep Potential Molecular Dynamics: a scalable model with the accuracy of quantum mechanics
• DeePMD-kit: A deep learning package for many-body potential energy representation and molecular dynamics
• Ab initio theory and modeling of water
• A Variational Approach to Enhanced Sampling and Free Energy Calculations
• Enhanced Sampling in the Well-Tempered Ensemble
• Why does hydronium diffuse faster than hydroxide in liquid water?
• Deep neural network for the dielectric response of insulators
• Raman Spectrum and Polarizability of Liquid Water from Deep Neural Networks
• Monte Carlo simulation of ice models
• Calculation of phase diagrams in the multithermal-multibaric ensemble
• A microscopic description of acid-base equilibrium
• Residual Entropy of Ordinary Ice from Multicanonical Simulations
• Configurational entropy of ice from thermodynamic integration
• Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations
• Electronic and Nuclear Quantum Effects on the Ice XI/Ice Ih Phase Transition
• Proton strings and rings in atypical nucleation of ferroelectricity in ice

Cited by in corpus (2)

• Deep Potentials for Materials Science
• Homogeneous ice nucleation in an ab initio machine learning model of water