The phase diagram of ice: a quasi-harmonic study based on a flexible water model

The phase diagram of ice is studied by a quasi-harmonic approximation. The free energy of all experimentally known ice phases has been calculated with the flexible q-TIP4P/F model of water. The only exception is the high pressure ice X, in which the presence of symmetric O-H-O bonds prevents its modeling with this empirical interatomic potential. The simplicity of our approach allows us to study ice phases at state points of the T-P plane that have been omitted in previous simulations using free energy methods based on thermodynamic integration. The effect in the phase diagram of averaging the proton disorder that appears in several ice phases has been studied. It is found particularly relevant for ice III, at least for cell sizes typically used in phase coexistence simulations. New insight into the capability of the employed water model to describe the coexistence of ice phases is presented. We find that the H-ordered ices IX and XIV, as well as the H-disordered ice XII, are particularly stable for this water model. This fact disagrees with experimental data. The unexpected large stability of ice IX is a property related to the TIP4P-character of the water model. Only after omission of these three stable ice phases, the calculated phase diagram becomes in reasonable qualitative agreement to the experimental one in the T-P region corresponding to ice Ih, II, III, V, and VI. The calculation of the phase diagram in the quantum and classical limits shows that the most important quantum effect is the stabilization of ice II due to its lower zero-point energy when compared to that one of ices Ih, III, and V.Comment: 13 pages, 8 figures, 5 table

Molar volume of solid isotopic helium mixtures

Solid isotopic helium mixtures have been studied by path-integral Monte Carlo simulations in the isothermal-isobaric ensemble. This method allowed us to study the molar volume as a function of temperature, pressure, and isotopic composition. At 25 K and 0.2 GPa, the relative difference between molar volumes of isotopically-pure crystals of 3He and 4He is found to be about 3%. This difference decreases under pressure, and for 12 GPa it is smaller than 1%. For isotopically-mixed crystals, a linear relation between lattice parameters and concentrations of helium isotopes is found, in agreement with Vegard's law. The virtual crystal approximation, valid for isotopic mixtures of heavier atoms, does not give reliable results for solid solutions of helium isotopes.Comment: 7 pages, 5 figure