A Blind Structure Prediction of Ice XIV

Experimental evidence for the existence of a second partially-ordered phase of ice VI

Ice exhibits extraordinary structural variety in its polymorphic structures. The existence of a new form of diversity in ice polymorphism has recently been debated in both experimental and theoretical studies, questioning whether hydrogen-disordered ice can transform into multiple hydrogen-ordered phases, contrary to the known one-to-one correspondence between disordered ice and its ordered phase. Here, we report a high-pressure phase, ice XIX, which is a second hydrogen-partially-ordered phase of ice VI. We demonstrate that disordered ice undergoes different manners of hydrogen ordering, which are thermodynamically controlled by pressure in the case of ice VI. Such multiplicity can appear in all disordered ice, and it widely provides a research approach to deepen our knowledge, for example of the crucial issues of ice: the centrosymmetry of hydrogen-ordered configurations and potentially induced (anti-)ferroelectricity. Ultimately, this research opens up the possibility of completing the phase diagram of ice.

Prediction of properties from first principles with quantitative accuracy: six representative ice phases

Based on a high-level MP2 theory with the fragment approach, the crystal structure, vibration spectra and phase transitions of six representative ice phases (II, VI, VII, VIII, IX, and XV) are predicted.

A new theoretical model for hexagonal ice, Ih(d), from first principles investigations

Due to their great importance in science, technology, and the life sciences, water and ice have been extensively investigated over many years. In particular, hexagonal ice Ih has been of great interest since it is the most common form of ice, and several modifications, Ih(a), Ih(b) and Ih(c) are known, whose structural details are still under discussion. In this study, we present an alternative theoretical model, called Ih(d), for the hexagonal ice modification in space group P63/mmc (no. 194), based on first-principles calculations that have been performed using DFT-LDA, GGA-PBE, and hybrid B3LYP and PBE0 functionals.

Thermodynamic and kinetic isotope effects on the order-disorder transition of ice XIV to ice XII

Isotope effects accompanying the order-disorder transition of ice XIV to ice XII are studied using calorimetry, X-ray diffraction, and dielectric spectroscopy. Particular emphasis is placed on the impact of the cooling rate applied during high-pressure production and during ambient-pressure recooling on the degree of hydrogen order in the low-temperature ice XIV phase. For specimens from D2O, ordering is harder to achieve in the sense that despite smaller cooling rates, the degree of order is less than in crystals produced from H2O. The degree of ordering can be quantified in terms of the Pauling entropy using calorimetry and manifests itself in structural and dynamical features that were examined using X-ray diffraction and dielectric spectroscopy, respectively. In hydrogen chloride doped samples, H/D substitution was found to slow down the dipolar dynamics up to about 30-fold and shifts the order-disorder transition by 4-6 K. By contrast to earlier assumptions it is possible to reach a high degree of ordering also at ambient pressure, provided the cooling rate is small enough. That is, at ambient pressure, orthorhombic stress slows down the dipolar reorientation near the ordering transition by a factor of 300-2000 for H2O and 30-100 for D2O samples. Furthermore, by long-term storage of our samples at 77 K we have reached surprisingly large increases in degree of order. For the D2O samples we observed an unprecedented high order, corresponding to more than 45% of the Pauling entropy.

DFT Simulations of the Vibrational Spectrum and Hydrogen Bonds of Ice XIV

It is always a difficult task to assign the peaks recorded from a vibrational spectrum. Herein, we explored a new pathway of density functional theory (DFT) simulation to present three kinds of spectra of ice XIV that can be referenced as inelastic neutron scattering (INS), infrared (IR), and Raman experimental spectrum. The INS spectrum is proportional to the phonon density of states (PDOS) while the photon scattering signals reflect the normal vibration frequencies near the Brillouin zone (BZ) center. Based on good agreements with the experimental data, we identified the relative frequency and made scientific assignments through normal vibration modes analysis. The two hydrogen bond (H-bond) peaks among the ice phases from INS were discussed and the dynamic process of the H-bond vibrations was found to be classified into two basic modes. We deduced that two H-bond modes are a general rule among the ice family and more studies are ongoing to investigate this subject.

Mapping uncharted territory in ice from zeolite networks to ice structures

Ice is one of the most extensively studied condensed matter systems. Yet, both experimentally and theoretically several new phases have been discovered over the last years. Here we report a large-scale density-functional-theory study of the configuration space of water ice. We geometry optimise 74,963 ice structures, which are selected and constructed from over five million tetrahedral networks listed in the databases of Treacy, Deem, and the International Zeolite Association. All prior knowledge of ice is set aside and we introduce "generalised convex hulls" to identify configurations stabilised by appropriate thermodynamic constraints. We thereby rediscover all known phases (I-XVII, i, 0 and the quartz phase) except the metastable ice IV. Crucially, we also find promising candidates for ices XVIII through LI. Using the "sketch-map" dimensionality-reduction algorithm we construct an a priori, navigable map of configuration space, which reproduces similarity relations between structures and highlights the novel candidates. By relating the known phases to the tractably small, yet structurally diverse set of synthesisable candidate structures, we provide an excellent starting point for identifying formation pathways.

Modeling Polymorphic Molecular Crystals with Electronic Structure Theory

Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.

Two Dimensional Ice from First Principles: Structures and Phase Transitions

Despite relevance to disparate areas such as cloud microphysics and tribology, major gaps in the understanding of the structures and phase transitions of low-dimensional water ice remain. Here, we report a first principles study of confined 2D ice as a function of pressure. We find that at ambient pressure hexagonal and pentagonal monolayer structures are the two lowest enthalpy phases identified. Upon mild compression, the pentagonal structure becomes the most stable and persists up to ∼2  GPa, at which point the square and rhombic phases are stable. The square phase agrees with recent experimental observations of square ice confined within graphene sheets. This work provides a fresh perspective on 2D confined ice, highlighting the sensitivity of the structures observed to both the confining pressure and the width.

Communication: On the stability of ice 0, ice i, and I(h)

Using ab initio methods, we examine the stability of ice 0, a recently proposed tetragonal form of ice implicated in the homogeneous freezing of water [J. Russo, F. Romano, and H. Tanaka, Nat. Mater. 13, 670 (2014)]. Vibrational frequencies are computed across the complete Brillouin Zone using Density Functional Theory (DFT), to confirm mechanical stability and quantify the free energy of ice 0 relative to ice I(h). The robustness of this result is tested via dispersion corrected semi-local and hybrid DFT, and Quantum Monte-Carlo calculation of lattice energies. Results indicate that popular molecular models only slightly overestimate the stability of ice zero. In addition, we study all possible realisations of proton disorder within the ice zero unit cell, and identify the ground state as ferroelectric. Comparisons are made to other low density metastable forms of ice, suggesting that the ice i structure [C. J. Fennel and J. D. Gezelter, J. Chem. Theory Comput. 1, 662 (2005)] may be equally relevant to ice formation.

Proton Ordering of Cubic Ice Ic: Spectroscopy and Computer Simulations

Several proton-disordered crystalline ice structures are known to proton order at sufficiently low temperatures, provided that the right preparation procedure is used. For cubic ice, ice Ic, however, no proton ordering has been observed so far. Here, we subject ice Ic to an experimental protocol similar to that used to proton order hexagonal ice. In situ FT-IR spectroscopy carried out during this procedure reveals that the librational band of the spectrum narrows and acquires a structure that is observed neither in proton-disordered ice Ic nor in ice XI, the proton-ordered variant of hexagonal ice. On the basis of vibrational spectra computed for ice Ic and four of its proton-ordered variants using classical molecular dynamics and ab initio simulations, we conclude that the features of our experimental spectra are due to partial proton ordering, providing the first evidence of proton ordering in cubic ice. We further find that the proton-ordered structure with the lowest energy is ferroelectric, while the structure with the second lowest energy is weakly ferroelectric. Both structures fit the experimental spectral similarly well such that no unique assignment of proton order is possible based on our results.

Hydrogen bonds and van der waals forces in ice at ambient and high pressures

The first principles methods, density-functional theory and quantum Monte Carlo, have been used to examine the balance between van der Waals (vdW) forces and hydrogen bonding in ambient and high-pressure phases of ice. At higher pressure, the contribution to the lattice energy from vdW increases and that from hydrogen bonding decreases, leading vdW to have a substantial effect on the transition pressures between the crystalline ice phases. An important consequence, likely to be of relevance to molecular crystals in general, is that transition pressures obtained from density-functional theory exchange-correlation functionals which neglect vdW forces are greatly overestimated.

Ice XV: a new thermodynamically stable phase of ice

A new phase of ice, named ice XV, has been identified and its structure determined by neutron diffraction. Ice XV is the hydrogen-ordered counterpart of ice VI and is thermodynamically stable at temperatures below approximately 130 K in the 0.8 to 1.5 GPa pressure range. The regions of stability in the medium pressure range of the phase diagram have thus been finally mapped, with only hydrogen-ordered phases stable at 0 K. The ordered ice XV structure is antiferroelectric (P1), in clear disagreement with recent theoretical calculations predicting ferroelectric ordering (Cc).

Liquid water and ices: understanding the structure and physical properties

A review of the structure and some properties of condensed phases of water is given. Since the discovery of the polymorphism of crystalline ice (beginning of the twentieth century), 15 ice modifications have been found and their structures have been determined. If we do not take into consideration proton ordering or disordering, nine distinct crystalline ice modifications in which water molecules retain their individuality are known. In the tenth, ice X, there are no H(2)O molecules. It contains ions (or atoms) of oxygen and hydrogen. The structure of all these modifications is described and information about their fields of stability and about the transition between them is given. It is emphasized that there are ice modifications which are metastable at any temperature and pressure (ices Ic, IV and XII), and many modifications can exist as metastable phases beyond their fields of stability. The ability of water to exist in metastable states is one of its remarkable properties. Several amorphous ice modifications (all of them are metastable) are known. Brief information about their properties and transitions between them is given. At the end of the 1960s the conception of the water structure as a three-dimensional hydrogen-bonded network was conclusively formed. Discovery of the polymorphism of amorphous ices awakened interest in the heterogeneity of the water network. Structural and dynamical heterogeneity of liquid water is discussed in detail. Computer simulation showed that the diffusion coefficient of water molecules in dense regions of the network is lower than in the loose regions, while an increase of density of the entire network gives rise to an increase of diffusion coefficient. This finding contradicts the conceptions associated with the primitive two-state models and can be explained from pressure dependences of melting temperature and of homogeneous nucleation temperature. A brief discussion of the picture of molecular motions in liquid water based on experiment and on computer simulation is given. This picture is still very incomplete. The most fascinating idea that was put forward during the last 20 years was the second critical point conjecture. It is still not clear whether this conjecture corresponds to reality.

The melting point of hexagonal ice (Ih) is strongly dependent on the quadrupole of the water models

It is shown that quadrupolar interactions play a determinant role in the melting temperatures of common water models and that there is a simple relationship between the strength of the quadrupolar forces and the position of the negative charge; our conclusion is that acceptable predictions for the melting temperature can only be obtained when the negatively charged site is shifted from the oxygen atom towards the hydrogens.