Quantum path integral simulation of isotope effects in the melting temperature of ice Ih.

The isotope effect in the melting temperature of ice Ih has been studied by free energy calculations within the path integral formulation of statistical mechanics. Free energy differences between isotopes are related to the dependence of their kinetic energy on the isotope mass. The water simulations were performed by using the q-TIP4P/F model, a point charge empirical potential that includes molecular flexibility and anharmonicity in the OH stretch of the water molecule. The reported melting temperature at ambient pressure of this model (T=251 K) increases by 6.5±0.5 and 8.2±0.5 K upon isotopic substitution of hydrogen by deuterium and tritium, respectively. These temperature shifts are larger than the experimental ones (3.8 and 4.5 K, respectively). In the classical limit, the melting temperature is nearly the same as that for tritiated ice. This unexpected behavior is rationalized by the coupling between intermolecular interactions and molecular flexibility. This coupling makes the kinetic energy of the OH stretching modes larger in the liquid than in the solid phase. However, the opposite behavior is found for intramolecular modes, which display larger kinetic energy in ice than in liquid water.

[1]  P. Rossky,et al.  Model dependence of quantum isotope effects in liquid water , 1991 .

[2]  Dietmar Paschek.  Temperature dependence of the hydrophobic hydration and interaction of simple solutes: an examination of five popular water models. , 2004, The Journal of chemical physics.

[3]  C. Vega,et al.  Heat capacity of water: A signature of nuclear quantum effects. , 2010, The Journal of chemical physics.

[4]  R. A. Kuharski,et al.  A quantum mechanical study of structure in liquid H2O and D2O , 1985 .

[5]  Sidney Yip,et al.  Optimized free-energy evaluation using a single reversible-scaling simulation , 1999 .

[6]  J. Jouzel,et al.  Climate and atmospheric history of the past 420,000 years from the Vostok ice core, Antarctica , 1999, Nature.

[7]  Peter S. Pacheco.  Parallel programming with MPI , 1996 .

[8]  Mark E. Tuckerman,et al.  Molecular dynamics algorithms for path integrals at constant pressure , 1999 .

[9]  C. Chakravarty.  Path integral simulations of atomic and molecular systems , 1997 .

[10]  C. Vega,et al.  Quantum contributions in the ice phases: the path to a new empirical model for water-TIP4PQ/2005. , 2009, The Journal of chemical physics.

[11]  R. Ludwig,et al.  Quantum cluster equilibrium theory of liquids: light and heavy QCE/3-21G model water , 2000 .

[12]  Emilio Artacho,et al.  Electrons and hydrogen-bond connectivity in liquid water. , 2005, Physical review letters.

[13]  Gillan.  Quantum simulation of hydrogen in metals. , 1988, Physical review letters.

[14]  V. Buch,et al.  Simulations of H2O Solid, Liquid, and Clusters, with an Emphasis on Ferroelectric Ordering Transition in Hexagonal Ice , 1998 .

[15]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[16]  C. Herrero,et al.  Quantum path-integral study of the phase diagram and isotope effects of neon. , 2011, The Journal of chemical physics.

[17]  F. Stillinger,et al.  Molecular Dynamics Study of Liquid Water , 1971 .

[18]  Watanabe,et al.  Direct dynamical calculation of entropy and free energy by adiabatic switching. , 1990, Physical review letters.

[19]  G. Voth,et al.  Quantum effects in liquid water from an ab initio-based polarizable force field. , 2007, The Journal of chemical physics.

[20]  C. Herrero,et al.  Path integral calculation of free energies: quantum effects on the melting temperature of neon. , 2008, The Journal of chemical physics.

[21]  S. Yoo,et al.  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. , 2009, The Journal of chemical physics.

[22]  B. Berne,et al.  Quantum effects in liquid water: Path-integral simulations of a flexible and polarizable ab initio model , 2001 .

[23]  D. Ceperley.  Path integrals in the theory of condensed helium , 1995 .

[24]  John A. Zollweg,et al.  The Lennard-Jones equation of state revisited , 1993 .

[25]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[26]  M. Tuckerman,et al.  IN CLASSICAL AND QUANTUM DYNAMICS IN CONDENSED PHASE SIMULATIONS , 1998 .

[27]  Roberto Car,et al.  Nuclear quantum effects in water. , 2008, Physical review letters.

[28]  Mark E. Tuckerman,et al.  Explicit reversible integrators for extended systems dynamics , 1996 .

[29]  Julian Tirado-Rives,et al.  Potential energy functions for atomic-level simulations of water and organic and biomolecular systems. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[30]  B. Berne,et al.  On path integral Monte Carlo simulations , 1982 .

[31]  M. Parrinello,et al.  Study of an F center in molten KCl , 1984 .

[32]  M. Shiga,et al.  Calculation of heat capacities of light and heavy water by path-integral molecular dynamics. , 2005, The Journal of chemical physics.

[33]  C. Vega,et al.  A general purpose model for the condensed phases of water: TIP4P/2005. , 2005, The Journal of chemical physics.

[34]  Thomas E. Markland,et al.  Competing quantum effects in the dynamics of a flexible water model. , 2009, The Journal of chemical physics.

[35]  C. Vega,et al.  The melting temperature of the most common models of water. , 2005, The Journal of chemical physics.

[36]  Wei Zhang,et al.  An accurate and simple quantum model for liquid water. , 2006, The Journal of chemical physics.

[37]  Michael W. Mahoney,et al.  Quantum, intramolecular flexibility, and polarizability effects on the reproduction of the density anomaly of liquid water by simple potential functions , 2001 .

[38]  M. Klein,et al.  Hydrogen bonding in water. , 2003, Physical review letters.

[39]  J. A. Barker,et al.  Structure of water; A Monte Carlo calculation , 1969 .

[40]  Nikos L. Doltsinis,et al.  Quantum Simulations of complex many-body systems:from theory to algorithms , 2002 .

[41]  J. Grotendorst,et al.  Quantum Simulations of Complex Many-Body Systems: From Theory to Algorithms , 2001 .

[42]  M. Tuckerman,et al.  Path integral molecular dynamics: a computational approach to quantum statistical mechanics , 1998 .

[43]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[44]  C. Vega,et al.  Quantum effects on the maximum in density of water as described by the TIP4PQ/2005 model. , 2009, The Journal of chemical physics.

[45]  P. Kusalik,et al.  Quantum effects in light and heavy liquid water: A rigid-body centroid molecular dynamics study. , 2004, The Journal of chemical physics.

[46]  J. Morales,et al.  Path integral simulation of the free energy of (Lennard-Jones) neon , 1991 .

[47]  C. Herrero.  Isotope effects in structural and thermodynamic properties of solid neon , 2001 .