Quantum path integral simulation of isotope effects in the melting temperature of ice Ih
arXiv:1108.3718 · doi:10.1063/1.3503764
Abstract
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 K 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 K 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.
14 pages, 13 figures, 5 tables
References in corpus (5)
- Competing quantum effects in the dynamics of a flexible water model
- Nuclear quantum effects in water
- Quantum contributions in the ice phases: the path to a new empirical model for water -- TIP4PQ/2005
- Path-integral molecular dynamics simulation of 3C-SiC
- Quantum path-integral study of the phase diagram and isotope effects of neon
Cited by in corpus (21)
- Ab initio thermodynamics of liquid and solid water
- Simulation and understanding of quantum crystals
- Anomalous Nuclear Quantum Effects in Ice
- Direct Measurement of Competing Quantum Effects on the Kinetic Energy of Heavy Water upon Melting
- Quantum-mechanical exploration of the phase diagram of water
- Quasi-harmonic approximation of thermodynamic properties of ice Ih, II, and III
- Isotope effects in ice Ih: A path-integral simulation
- Anharmonic nuclear motion and the relative stability of hexagonal and cubic ice
- The phase diagram of water from quantum simulations
- Kinetic energy of protons in ice Ih and water: a path integral study
- Simulation of quantum zero-point effects in water using a frequency-dependent thermostat
- A study of the influence of isotopic substitution on the melting point and temperature of maximum density of water by means of path integral simulations of rigid models
- Quantum Nature of the Hydrogen Bond from Ambient Conditions down to Ultra-low Temperatures
- The phase diagram of ice Ih, II, and III: a quasi-harmonic study
- High-density amorphous ice: A path-integral simulation
- Path-integral simulation of ice Ih: The effect of pressure
- Path-integral simulation of ice VII: Pressure and temperature effects
- Topological characterization of crystalline ice structures from coordination sequences
- Isotopic effects in structural properties of graphene
- Water at negative pressure: Nuclear quantum effects
- h-CMD: An efficient hybrid fast centroid and quasi-centroid molecular dynamics method for the simulation of vibrational spectra