New porous water ice metastable at atmospheric pressure obtained by emptying a hydrogen-filled ice

Harnessing Nuclear Magnetic Resonance Spectroscopy to Decipher Structure and Dynamics of Clathrate Hydrates in Confinement: A Perspective

This perspective outlines recent developments in the field of NMR spectroscopy, enabling new opportunities for in situ studies on bulk and confined clathrate hydrates. These hydrates are crystalline ice-like materials, built up from hydrogen-bonded water molecules, forming cages occluding non-polar gaseous guest molecules, including CH4, CO2 and even H2 and He gas. In nature, they are found in low-temperature and high-pressure conditions. Synthetic confined versions hold immense potential for energy storage and transportation, as well as for carbon capture and storage. Using previous studies, this report highlights static and magic angle spinning NMR hardware and strategies enabling the study of clathrate hydrate formation in situ, in bulk and in nano-confinement. The information obtained from such studies includes phase identification, dynamics, gas exchange processes, mechanistic studies and the molecular-level elucidation of the interactions between water, guest molecules and confining interfaces.

Large-cage occupation and quantum dynamics of hydrogen molecules in sII clathrate hydrates

Hydrogen clathrate hydrates are ice-like crystalline substances in which hydrogen molecules are trapped inside polyhedral cages formed by the water molecules. Small cages can host only a single H2 molecule, while each large cage can be occupied by up to four H2 molecules. Here, we present a neutron scattering study on the structure of the sII hydrogen clathrate hydrate and on the low-temperature dynamics of the hydrogen molecules trapped in its large cages, as a function of the gas content in the samples. We observe spectral features at low energy transfer (between 1 and 3 meV), and we show that they can be successfully assigned to the rattling motion of a single hydrogen molecule occupying a large water cage. These inelastic bands remarkably lose their intensity with increasing the hydrogen filling, consistently with the fact that the probability of single occupation (as opposed to multiple occupation) increases as the hydrogen content in the sample gets lower. The spectral intensity of the H2 rattling bands is studied as a function of the momentum transfer for partially emptied samples and compared with three distinct quantum models for a single H2 molecule in a large cage: (i) the exact solution of the Schrödinger equation for a well-assessed semiempirical force field, (ii) a particle trapped in a rigid sphere, and (iii) an isotropic three-dimensional harmonic oscillator. The first model provides good agreement between calculations and experimental data, while the last two only reproduce their qualitative trend. Finally, the radial wavefunctions of the three aforementioned models, as well as their potential surfaces, are presented and discussed.

Assisted Self-Assembly of Nanoporous Ices via Carbon Nanomaterial Templates

Self-assembly is a widely used synthetic method in nanoscience to assemble well-organized structures. Self-assembly processes usually occur in a water solvent environment. However, the self-assembly of water molecules is rarely studied. Herein, we show a strategy to fabricate porous ice via carbon nanomaterial-assisted self-assembly. Diverse frameworks of nanoporous ice are formed by using orthorhombic and tetragonal arrays of carbon nanotubes or carbon-atom chains as templates. In contrast to many bulk ices discovered in nature, nanoporous ices are shown to be stable only under negative pressure. Hence, nanoporous ices cannot be produced through the direct nucleation of water at negative pressure. The template-assisted self-assembly method is shown to be the most effective method to fabricate nanoporous ice in quantity. Several key factors for the self-assembly of nanoporous ices are identified, including proper gap spacings in the carbon nanomaterial template and suitable interactions between water and the carbon nanomaterials.

Analysing the stability of He-filled hydrates: how many He atoms fit in the sII crystal?

Clathrate hydrates have the ability to encapsulate atoms and molecules within their cavities, and thus they could be potentially large storage capacity materials. The present work studies the multiple cage occupancy effects in the recently discovered He@sII crystal. On the basis of previous theoretical and experimental findings, the stability of He(1/1)@sII, He(1/4)@sII and He(2/4)@sII crystals was analysed in terms of structural, mechanical and energetic properties. For this purpose, first-principles DFT/DFT-D computations were performed by using both semi-local and non-local functionals, not only to elucidate which configuration is the most energetically favoured, but also to scrutinize the relevance of the long-range dispersion interactions in these kinds of compounds. We have encountered that dispersion interactions play a fundamental role in describing the underlying interactions, and different tendencies were observed depending on the choice of the functional. We found that PW86PBE-XDM shows the best performance, while the non-local functionals tested here were not able to correctly account for them. The present results reveal that the most stable configuration is the one presenting singly occupied D cages and tetrahedrally occupied H cages (He(1/4)@sII) in line with the experimental observation.

Observation of the most H2-dense filled ice under high pressure

Hydrogen hydrates are among the basic constituents of our solar system's outer planets, some of their moons, as well Neptune-like exo-planets. The details of their high-pressure phases and their thermodynamic conditions of formation and stability are fundamental information for establishing the presence of hydrogen hydrates in the interior of those celestial bodies, for example, against the presence of the pure components (water ice and molecular hydrogen). Here, we report a synthesis path and experimental observation, by X-ray diffraction and Raman spectroscopy measurements, of the most H[Formula: see text]-dense phase of hydrogen hydrate so far reported, namely the compound 3 (or C[Formula: see text]). The detailed characterisation of this hydrogen-filled ice, based on the crystal structure of cubic ice I (ice I[Formula: see text]), is performed by comparing the experimental observations with first-principles calculations based on density functional theory and the stochastic self-consistent harmonic approximation. We observe that the extreme (up to 90 GPa and likely beyond) pressure stability of this hydrate phase is due to the close-packed geometry of the hydrogen molecules caged in the ice I[Formula: see text] skeleton.

Computational investigations of stable multiple-cage-occupancy He clathrate-like hydrostructures

One of the several possibilities offered by the interesting clathrate hydrates is the opportunity to encapsulate several atoms or molecules, in such a way that more efficient storage materials could be explored or new molecules that otherwise do not exist could be created. These types of applications are receiving growing attention from technologists and chemists, given the future positive implications that they entail. In this context, we investigated the multiple cage occupancy of helium clathrate hydrates, to establish stable novel hydrate structures or ones similar to those predicted previously by experimental and theoretical studies. To this purpose, we analyzed the feasibility of including an increased number of He atoms inside the small (D) and large (H) cages of the sII structure through first-principles properly assessed density functional approaches. On the one hand, we have computed energetic and structural properties, in which we examined the guest-host and guest-guest interactions in both individual and two-adjacent clathrate-like sII cages by means of binding and evaporation energies. On the other hand, we have carried out a thermodynamical analysis on the stability of such He-containing hydrostructures in terms of changes in enthalpy, ΔH, Gibbs free energy, ΔG, and entropy, ΔS, during their formation process at various temperature and pressure values. In this way, we have been able to make a comparison with experiments, reaffirming the ability of computational DFT approaches to describe such weak guest-host interactions. In principle, the most stable structure involves the encapsulation of one and four He atoms inside the D and H sII cages, respectively; however, more He atoms could be entrapped under lower temperature and/or higher pressure thermodynamic conditions. We foresee such accurate computational quantum chemistry approaches contributing to the current emerging machine-learning model development.

Computational Analysis of Vibrational Spectra of Hydrogen Bonds in sII and sH Gas Hydrates

The amount of energy in natural gas hydrates is thought to be equivalent to twice that of all other fossil fuels combined. However, economic and safe energy recovery has remained a challenge till now. To develop a novel method of breaking the hydrogen bonds (HBs) surrounding the trapped gas molecules, we investigated the vibrational spectra of the HBs of gas hydrates with structure types II and H. Two models of 576-atom propane-methane sII hydrate and 294-atom neohexane-methane sH hydrate were built. A first-principles density functional theory (DFT) method was employed using the CASTEP package. The simulated spectra were in good agreement with the experimental data. Compared with the partial phonon density of states of guest molecules, we confirmed that the experimental infrared absorption peak in the terahertz region mainly arose from HB vibrations. By removing the components of guest molecules, we found that the theory of two kinds of hydrogen bond vibrational modes applies. The use of a terahertz laser to enable resonance absorption of HBs (at about 6 THz, to be tested) may therefore lead to the rapid melting of clathrate ice and release of guest molecules.

Nanoporous ices: an emerging class in the water/ice family

The history of scientific research on diverse ice structures dates back to more than a century. To date, 20 three-dimensional crystalline ice phases (ice I-ice XX) have been identified in the laboratory, among which ice XVI and ice XVII belong to a class of low-density nanoporous ices. Nanoporous ices can also be viewed as a special class of porous materials or water ice, as they possess a relatively high fraction of nano-cavities and/or nano-channels built into the hydrogen-bonded water framework. As such, like the prototypical class of porous materials (e.g., MOFs and COFs), nanoporous ices can be named as water oxygen-vertex frameworks (WOFs). Because of their large surface-to-volume ratio, WOFs may be potential media for gas storage, gas purification and separation. They may be applied to the biomedical field owing to their excellent biocompatibility. The field of porous ices is still emerging, as many porous ice structures that are predicted to be stable by computer simulations require future experimental confirmation. For future theoretical/computational studies, as the machine-learning method becomes an increasingly popular research tool in the material science and chemical science fields, more reliable porous ice structures and phase diagrams will be predicted with the development of more accurate machine-learning force fields.

The pressure induced phase diagram of double-layer ice under confinement: a first-principles study

Here, we present double-layer ice confined within various carbon nanotubes (CNTs) using state-of-the-art pressure induced (-5 GPa to 5 GPa) dispersion corrected density functional theory (DFT) calculations. We find that the double-layer ice exhibits remarkably rich and diverse phase behaviors as a function of pressure with varying CNT diameters. The lattice cohesive energies for various pure double layer ice phases follow the order of hexagonal > pentagonal > square tube > hexagonal-close-packed (HCP) > square > buckled-rhombic (b-RH). The confinement width was found to play a crucial role in the square and square tube phases in the intermediate pressure range of about 0-1 GPa. Unlike the phase transition in pure bilayer ice structures, the relative enthalpies demonstrate that the pentagonal phase, rather than the hexagonal structure, is the most stable ice polymorph at ambient pressure as well as in a deep negative pressure region, whereas the b-RH phase dominates under high pressure. The relatively short O⋯O distance of b-RH demonstrates the presence of a strong hydrogen bonding network, which is responsible for stabilizing the system.

Delving into guest-free and He-filled sI and sII clathrate hydrates: a first-principles computational study

The dynamics of the formation of a specific clathrate hydrate as well as its thermodynamic transitions depend on the interactions between the trapped molecules and the host water lattice. The molecular-level understanding of the different underlying processes benefits not only the description of the properties of the system, but also allows the development of multiple technological applications such as gas storage, gas separation, energy transport, etc. In this work we investigate the stability of periodic crystalline structures, such as He@sI and He@sII clathrate hydrates by first-principles computations. We consider such host water networks interacting with a guest He atom using selected density functional theory approaches, in order to explore the effects on the encapsulation of a light atom in the sI/sII crystals, by deriving all energy components (guest-water, water-water, guest-guest). Structural properties and energies were first computed by structural relaxations of the He-filled and empty sI/sII unit cells, yielding lattice and compressibility parameters comparable to experimental and theoretical values available for those hydrates. According to the results obtained, the He enclathration in the sI/sII unit cells is a stabilizing process, and both He@sI and He@sII clathrates, considering single cage occupancy, are predicted to be stable whatever the XDM or D4 dispersion correction applied. Our results further reveal that despite the weak underlying interactions the He encapsulation has a rather notable effect on both lattice parameters and energetics, with the He@sII being the most energetically favorable in accord with recent experimental observations.

Comparative Analysis of the Hydrogen Bond Vibrations of Ice XII

It is difficult to theoretically study the vibrational spectrum of hydrogen-disordered ice XII compared with its hydrogen-ordered counterpart, ice XIV. We constructed a 24-molecule supercell of ice XII to mimic its real structure. We focused on hydrogen bond (HB) vibrational modes in the translation band using first-principles density functional theory (DFT). Our simulated results were in good agreement with inelastic neutron scattering experiments. We found that the optical vibrational modes of HBs are composed of three main components. These are cluster vibrations in the lowest-frequency region, four-bond HB vibrations in the highest-frequency region, and two-bond modes in between. Although the experimentally recorded curve of ice XII is smooth in the translation region, our results support the proposal that two types of intrinsic HB vibrational modes are common in the ice family.

Comparison of the phase transitions of high-pressure phases of ammonium fluoride and ice at ambient pressure

The phase diagrams of water and ammonium fluoride (NH₄F) display some interesting parallels. Several crystalline NH₄F phases have isostructural ice counterparts and one of the famous anomalies of water, the fact that the liquid is denser than ice Ih, is also found for NH₄F. Here, we investigate the phase transitions of the pressure-quenched high-pressure phases of NH₄F upon heating at ambient pressure with x-ray diffraction and calorimetry, and we compare the results with the corresponding ices. NH₄F II transforms to NH₄F Isd, which is a stacking-disordered variant of the stable hexagonal NH₄F Ih polymorph. Heating NH₄F III gives a complex mixture of NH₄F II and NH₄F Isd, while some NH₄F III remains initially. Complete conversion to NH₄F Isd is achieved above ∼220 K. The NH₄F II obtained from NH₄F III persists to much higher temperatures compared to the corresponding pressure-quenched NH₄F II. Quantification of the stacking disorder in NH₄F Isd reveals a more sluggish conversion to NH₄F Ih for NH₄F Isd from NH₄F III. In general, the presence of stress and strain in the samples appears to have pronounced effects on the phase transition temperatures. NH₄F shows a complete lack of amorphous forms at low temperatures either upon low-temperature compression of NH₄F Ih or heating NH₄F III at ambient pressure. The amorphous forms of ice are often used to explain the anomalies of water. It will, therefore, be interesting to explore if liquid NH₄F displays more water-like anomalies despite the apparent lack of amorphous forms at low temperatures.

Computer simulation of hypothetical hydrogen ordered structure of ice XIX

A new ice phase, ice XIX, was discovered in 2018, and is the second hydrogen-ordered polymorph of hydrogen-disordered ice VI. The first hydrogen-ordered polymorph of ice VI is ice XV,...

Assessment of DFT approaches in noble gas clathrate-like clusters: stability and thermodynamics

We have assessed the performance and accuracy of different wavefunction-based electronic structure methods, such as DFMP2 and domain-based local pair-natural orbital (DLPNO-CCSD(T)), as well as a variety of density functional theory (DFT) approaches on He@(H₂O)_N cage systems. We have selected representative clathrate-like structures corresponding to the building blocks present in each of the sI, sII and sH natural gas clathrate hydrates, and we have carefully studied the interaction between a He atom with each of their individual cages. We reported well-converged DFMP2 and DLPNO-CCSD(T) reference data, together with interaction and cohesive energies of four different density functionals (two GGA, revPBE and PW86PBE, and two hybrids, B3LYP and PBE0), including diverse dispersion correction schemes (D3(0), D3(BJ), D4 and XDM) for both He-filled and empty clathrate-like cages. After the analysis of the results, we came to the conclusion that the PW86PBE functional, with both XDM and D4 corrections, and the PBE0-D4 functional present reasonably adequate approaches to describe the guest-host noncovalent interactions that take place in such He clathrate hydrates. Taking into account that the He@sII is the only helium clathrate that scientists have been able to synthesize recently, we have performed a thermodynamic study on the individual 5¹² and 5¹²6⁴ cages present in the sII crystal. We determined the change in enthalpy, ΔH, and in Gibbs free energy, ΔG, at various temperatures and pressures, and we found out that in the range of experimental conditions the reactions associated with the encapsulation of the He atom inside the cages are exothermic and spontaneous. Finally, we highlighted the importance of an accurate description of the interaction in He@water mixtures, as a crucial component in construction of reliable data-driven models.

Capturing Individual Hydrogen Bond Strengths in Ices via Periodic Local Vibrational Mode Theory: Beyond the Lattice Energy Picture

Local stretching force constants derived from periodic local vibrational modes at the vdW-DF2 density functional level have been employed to quantify the intrinsic hydrogen bond strength of 16 ice polymorphs, ices Ih, II, III, IV, V, VI, VII, VIII, IX, XI, XII, XIII, XIV, XV, XVII, and XIX, that are stable under ambient to elevated pressures. Based on this characterization on 1820 hydrogen bonds, relationships between local stretching force constants and structural parameters such as hydrogen bond length and angle were identified. Moreover, different bond strength distributions, from uniform to inhomogeneous, were observed for the 16 ices and could be explained in relation to different local structural elements within ices, that is, rings, that consist of different hydrogen bond types. In addition, criteria for the classification of hydrogen bonds as strong, intermediate, and weak were introduced. The latter was used to explore a different dimension of the water-ice phase diagram. These findings will provide important guidelines for assessing the credibility of new ice structures.

On the role of intermolecular vibrational motions for ice polymorphs. III. Mode characteristics associated with negative thermal expansion

Low-pressure ice forms, such as hexagonal and cubic ice, expand on cooling below temperature 60 K. This negative thermal expansivity has been explored in terms of phonon frequency modulation with varying volume and attributed to the negative Grüneisen parameters unique mostly to tetrahedrally coordinated substances. However, an underlying mechanism for the negative Grüneisen parameters has not been known except some schematic analyses. We investigate in this study the characteristics of the intermolecular vibrational modes whose Grüneisen parameters are negative by examining the individual vibrational modes rigorously. It is found that the low frequency modes below 100 cm^-1, which we explicitly show are mostly bending motions of three hydrogen-bonded molecules, necessarily accompany elongation of the hydrogen bond length at peak amplitudes compared with that at the equilibrium position in executing the vibrational motions. The elongation gives rise to a decrease in the repulsive interaction while an increase in the Coulombic one. The decrease in the repulsive interaction is relaxed substantially by expansion due to its steep slope against molecular separation compared with the sluggish increase in the Coulombic one, and therefore, the negative Grüneisen parameters are obtainable. This scenario is tested against some variants of cubic ice with various water potential models. It is demonstrated that four interaction-site models are suitable to describe the intermolecular vibrations and the thermal expansivity because of the moderate tendency to favor the tetrahedral coordination.

Mining of Effective Local Order Parameters to Classify Ice Polymorphs

Order parameters make it possible to quantify the degree of structural ordering in a material and thus to apply as the reaction coordinates during the free-energy analysis of phase or structure transitions. Furthermore, order parameters are useful in determining the local structures of molecular groups during transition stages. However, identifying or developing local order parameters (LOPs) that are sensitive for specific materials and phases is a non-trivial task. In this study, the ability of LOPs to classify the solid and liquid structures of water at coexistence or triple points is investigated with the aid of supervised machine learning. The classification accuracy of a total of 179,738,433 combinations of 493 LOPs is automatically and systematically compared for water structures at the ice Ih-Ic-liquid coexistence point and the ice III-V-liquid and ice V-VI-liquid triple points. The optimal sets of two LOPs are found for each point, and sets of three LOPs are suggested for better accuracy.

Formation of porous ice frameworks at room temperature

Bulk crystalline ices with ultralow densities have been demonstrated to be thermodynamically metastable at negative pressures. However, the direct formation of these bulk porous ices from liquid water at negative pressures is extremely challenging. Inspired by approaches toward porous media based on host-guest chemistry, such as metal-organic frameworks and covalent organic frameworks, we herein demonstrate via molecular dynamics simulations that a class of ultralow-density porous ices with upright channels can be formed spontaneously from liquid water at 300 K with the assistance of carbon nanotube arrays. We refer to these porous ice polymorphs as water oxygen-vertex frameworks (WOFs). Notably, our simulations revealed that the liquid-WOF phase transition is first-order and occurs at room temperature. All the WOFs exhibited the unique structural feature that they can be regarded as assemblies of nanoribbons of hexagonal bilayer ice (2D ice I) at their armchair or zigzag edges. Based on density functional theory calculations, a comprehensive phase diagram of the WOFs was constructed considering both the thermodynamic and thermal stabilities of the porous ices at negative pressures. Like other types of porous media, these WOFs may be applicable to gas storage, purification, and separation. Moreover, these biocompatible porous ice networks may be exploited as medical-related carriers.

Insight into Liquid Polymorphism from the Complex Phase Behavior of a Simple Model

We systematically explored the phase behavior of the hard-core two-scale ramp model suggested by Jagla [Phys. Rev. E 63, 061501 (2001)PRESCM1539-375510.1103/PhysRevE.63.061501] using a combination of the nested sampling and free energy methods. The sampling revealed that the phase diagram of the Jagla potential is significantly richer than previously anticipated, and we identified a family of new crystalline structures, which is stable over vast regions in the phase diagram. We showed that the new melting line is located at considerably higher temperature than the boundary between the low- and high-density liquid phases, which was previously suggested to lie in a thermodynamically stable region. The newly identified crystalline phases show unexpectedly complex structural features, some of which are shared with the high-pressure ice VI phase.

Searching for local order parameters to classify water structures at triple points

The diversity of ice polymorphs is of interest in condensed-matter physics, engineering, astronomy, and biosphere and climate studies. In particular, their triple points are critical to elucidate the formation of each phase and transitions among phases. However, an approach to distinguish their molecular structures is lacking. When precise molecular geometries are given, order parameters are often computed to quantify the degree of structural ordering and to classify the structures. Many order parameters have been developed for specific or multiple purposes, but their capabilities have not been exhaustively investigated for distinguishing ice polymorphs. Here, 493 order parameters and their combinations are considered for two triple points involving the ice polymorphs ice III-V-liquid and ice V-VI-liquid. Supervised machine learning helps automatic and systematic searching of the parameters. For each triple point, the best set of two order parameters was found that distinguishes three structures with high accuracy. A set of three order parameters is also suggested for better accuracy.

Comparative Analysis of Hydrogen-Bonding Vibrations of Ice VI

It is difficult to investigate the hydrogen-bonding dynamics of hydrogen-disordered ice VI. Here, we present a comparative method based on our previous study of its counterpart hydrogen-ordered phase, ice XV. The primitive cell of ice XV is a 10 molecule unit, and the vibrational normal modes were analyzed individually. We constructed an 80 molecule supercell of ice VI to mimic the periodic unit and performed first-principles density functional theory calculations. As the two vibrational spectra show almost identical features, we compared the molecular translation vibrations. Inspired by the phonon analysis of ice XV, we found that the vibrational modes in the translation band of ice VI are classifiable into three groups. The lowest-strength vibration modes represent vibrations between two sublattices that lack hydrogen bonding. The highest-strength vibration modes represent the vibration of four hydrogen bonds of one molecule. The middle-strength vibration modes mainly represent the molecular vibrations of only two hydrogen bonds. Although there are many overlapping stronger and middle modes, there are only two main peaks in the inelastic neutron scattering (INS) spectra. This work explains the origin of the two main peaks in the far-infrared region of ice VI and illustrates how to analyze a hydrogen-disordered ice structure.

Ammonium fluoride's analogy to ice: Possibilities and limitations

Ammonium fluoride, NH₄F, is often seen as an analog to ice, with several of its solid phases closely resembling known ice phases. While its ionic and hydrogen-ordered nature puts topological constraints on the ice-like network structures it can form, it is not clear what consequences these constraints have for NH₄F compound formation and evolution. Here, we explore computationally the reach and eventual limits of the ice analogy for ammonium fluoride. By combining data mining of known and hypothetical ice networks with crystal structure prediction and density functional calculations, we explore the high-pressure phase diagram of NH₄F and host-guest compounds of its hydrides. Pure NH₄F departs from ice-like behavior above 80 GPa with the emergence of close-packed ionic structures. The predicted stability of NH₄F hydrides shows that NH₄F can act as a host to small guest species, albeit in a topologically severely constraint configuration space. Finally, we explore the binary NH₃-HF chemical space, where we find candidate structures for several unsolved polyfluoride phases; among them is the chemical analog to H₂O₂ dihydrate.

The everlasting hunt for new ice phases

Water ice exists in hugely different environments, artificially or naturally occurring ones across the universe. The phase diagram of crystalline phases of ice is still under construction: a high-pressure phase, ice XIX, has just been reported but its structure remains ambiguous.

Formation of hot ice caused by carbon nanobrushes. II. Dependency on the radius of nanotubes

Stable crystalline structures of confined water can be different from bulk ice. In Paper I [T. Yagasaki et al., J. Chem. Phys. 151, 064702 (2019)] of this study, it was shown, using molecular dynamics (MD) simulations, that a zeolite-like ice structure forms in nanobrushes consisting of (6,6) carbon nanotubes (CNTs) when the CNTs are located in a triangle arrangement. The melting temperature of the zeolite-like ice structure is much higher than the melting temperature of ice I_h when the distance between the surfaces of CNTs is ∼0.94 nm, which is the best spacing for the bilayer structure of water. In this paper, we perform MD simulations of nanobrushes of CNTs that are different from (6,6) CNTs in radius. Several new porous ice structures form spontaneously in the MD simulations. A stable porous ice forms when the radius of its cavities matches the radius of the CNTs well. All cylindrical porous ice structures found in this study can be decomposed into a small number of structural blocks. We provide a new protocol to classify cylindrical porous ice crystals on the basis of this decomposition.

Identification of synthesisable crystalline phases of water – a prototype for the challenges of computational materials design

We discuss the identification of experimentally realisable crystalline phases of water to outline and contextualise some of the diverse building blocks of a computational materials design process.

Novel Hydrogen Clathrate Hydrate

We report a new hydrogen clathrate hydrate synthesized at 1.2 GPa and 298 K documented by single-crystal x-ray diffraction, Raman spectroscopy, and first-principles calculations. The oxygen sublattice of the new clathrate hydrate matches that of ice II, while hydrogen molecules are in the ring cavities, which results in the trigonal R3c or R3[over ¯]c space group (proton ordered or disordered, respectively) and the composition of (H_{2}O)_{6}H_{2}. Raman spectroscopy and theoretical calculations reveal a hydrogen disordered nature of the new phase C_{1}^{'}, distinct from the well-known ordered C_{1} clathrate, to which this new structure transforms upon compression and/or cooling. This new clathrate phase can be viewed as a realization of a disordered ice II, unobserved before, in contrast to all other ordered ice structures.

He Inclusion in Ice-like and Clathrate-like Frameworks: A Benchmark Quantum Chemistry Study of Guest-Host Interactions

Energetics and structural properties of selected type and size He@hydrate frameworks, e.g., from regular structured ice channels to clathrate-like cages, are presented from first-principles quantum chemistry methods. The scarcity of information on He@hydrates makes such complexes challenging targets, while their computational study entails an interesting and arduous task. Some of them have been synthesized in the laboratory, which motivates further investigations on their stability. Hence, the main focus is to examine the performance and accuracy of different wave function-based electronic structure methods, such as MP2, CCSD(T), their explicitly correlated (F12) and domain-based local pair-natural orbital (DLPNO) analogs, as well as modern and conventional density functional theory (DFT) approaches, and analytical model potentials available. Different structures are considered, starting from the "simplest system" formed by a noble gas atom (such as He) and one water molecule, followed by the study of the "fundamental units" present in all ice-like and clathrate-like frameworks (such as pentamers and hexamers) and finally the description of interactions in the "building blocks" of three-dimensional (3D) ice channels (e.g., horizontal and perpendicular ice II and I_h) and clathrate-like cages, such as the 5¹² present in the most common sI, sII, and sH clathrate-hydrate structures. The idea is to provide well-converged DLPNO-CCSD(T) and DFMP2/CBS reference datasets that in turn are used to validate how DFT functionals (in total, 29 approaches from generalized-gradient approximation (GGA), meta-GGA, to hybrid and range-separated functionals, including dispersion correction treatments, were checked) and analytical semiempirical/ab initio-based potentials perform compared with high-level alternatives. Within all tested approaches, those best-performing were identified and classified. Most of the DFT/DFT-D functionals, as well as available analytical pairwise model potentials, face difficulties in describing both hydrogen-bonded water frameworks and dispersion bound He-water interactions. Including dispersion corrections yields an overall well-balanced performance for LCωPBE-D3BJ and PBE0-D4 functionals. Such benchmark datasets can benefit research into the development of new cheminformatics models, as can serve to guide and cross-check methodologies, lending increased predicted power to future molecular simulations for investigating the role of structures and phase transitions from nanoscale clusters to macroscopic crystalline structures.

Cubic ice Ic without stacking defects obtained from ice XVII

Amongst the more than 18 different forms of water ice, only the common hexagonal phase and the cubic phase are present in nature on Earth. Nonetheless, it is now widely recognized that all samples of 'cubic ice' discovered so far do not have a fully cubic crystal structure but instead are stacking-disordered forms of ice I (namely, ice Isd), which contain both hexagonal and cubic stacking sequences of hydrogen-bonded water molecules. Here, we describe a method to obtain large quantities of cubic ice Ic with high structural purity. Cubic ice Ic is formed by heating a powder of D₂O ice XVII obtained from annealing of pristine C₀ hydrate samples under dynamic vacuum. Neutron diffraction experiments performed on two different instruments and Raman spectroscopy measurements confirm the structural purity of the cubic ice, Ic. These findings contribute to a better understanding of ice I polymorphism and the existence of the two natural ice forms.

On the role of intermolecular vibrational motions for ice polymorphs. II. Atomic vibrational amplitudes and localization of phonons in ordered and disordered ices

We investigate the vibrational amplitudes and the degree of the phonon localization in 19 ice forms, both crystalline and amorphous, by a quasi-harmonic approximation with a reliable classical intermolecular interaction model for water. The amplitude in the low pressure ices increases with compression, while the opposite trend is observed in the medium and high pressure ices. The amplitude of the oxygen atom does not differ from that of hydrogen in low pressure ices apart from the contribution from the zero-point vibrations. This is accounted for by the coherent but opposite phase motions in the mixed translational and rotational vibrations. A decoupling of translation-dominant and rotation-dominant motions significantly reduces the vibrational amplitudes in any ice form. The amplitudes in ice III are found to be much larger than any other crystalline ice form. In order to investigate the vibrational mode characteristics, the moment ratio of the atomic displacements for individual phonon modes, called the inverse participation ratio, is calculated and the degree of the phonon localization in crystalline and amorphous ices is discussed. It is found that the phonon modes in the hydrogen-ordered ice forms are remarkably spread over the entire crystal having propagative or diffusive characteristic, while many localized modes appear at the edges of the vibrational bands, called dissipative modes, in the hydrogen-disordered counterparts. The degree of localization is little pronounced in low density amorphous and high density amorphous due to disordering of oxygen atoms.

Ice Ic without stacking disorder by evacuating hydrogen from hydrogen hydrate

Water freezes below 0 °C at ambient pressure ordinarily to ice I_h, with hexagonal stacking sequence. Under certain conditions, ice with a cubic stacking sequence can also be formed, but ideal ice I_c without stacking-disorder has never been formed until recently. Here we demonstrate a route to obtain ice I_c without stacking-disorder by degassing hydrogen from the high-pressure form of hydrogen hydrate, C₂, which has a host framework isostructural with ice I_c. The stacking-disorder free ice I_c is formed from C₂ via an intermediate amorphous or nano-crystalline form under decompression, unlike the direct transformations occurring in ice XVI from neon hydrate, or ice XVII from hydrogen hydrate. The obtained ice I_c shows remarkable thermal stability, until the phase transition to ice I_h at 250 K, originating from the lack of dislocations. This discovery of ideal ice I_c will promote understanding of the role of stacking-disorder on the physical properties of ice as a counter end-member of ice I_h.

Metastable cubic ice has been identified in several conditions relevant to geo and astrochemistry, but was always characterized by stacking disorder. Here the authors synthesize a hydrogen hydrate and degas hydrogen, obtaining pure non-defected cubic ice, observed by X-ray and neutron diffraction.

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.

How flexible is the water molecule structure? Analysis of crystal structures and the potential energy surface

High level ab initio calculations predicted a possibility for energetically low-cost (±1 kcal mol⁻¹) change of the bond angle and bond lengths in wide range,from 96.4° to 112.8° and from 0.930 Å to 0.989 Å, respectively.

Translational and Rotational Diffusion in Liquid Water at Very High Pressure: A Simulation Study

Translational and rotational diffusion coefficients of liquid water have been computed from molecular dynamics simulation with a recent polarizable potential at 298, 400, and 550 K at very high pressure. At 298 K, the model reproduces the initial increase and the occurrence of a maximum for the translational and rotational diffusion coefficients when the pressure increases. At 400 and 550 K, translational and rotational diffusion coefficients are found to monotonically decrease when pressure increases in the gigapascal range, with the translational coefficient decreasing faster than the rotational one. At 400 K, such an evolution of the rotational diffusion coefficient contrasts with quasielastic neutron scattering results predicting a near independence of the rotational diffusion with a pressure increase above ≃0.5 GPa. An interpretation is proposed to explain this discrepancy. The pressure dependence of the translation-rotation coupling is analyzed. The anisotropy of rotational diffusion is investigated by computing the rotational diffusion tensor in a molecular system of axes and the reorientational correlation times of rank 1 and rank 2 of the inertia axes and of the OH bond vector. Deviation of the simulation data with respect to the predictions of the isotropic Debye model of rotational diffusion are quantified and can be used to estimate experimental rotational diffusion coefficients from experimental reorientational correlation times.

Nucleation Curve of Carbon Dioxide Hydrate from a Linear Cooling Ramp Method

Formation of gas hydrate is a first-order phase transition that starts with nucleation. Understanding of nucleation is of interest to many in chemical and petroleum industries, as nucleation, while beneficial in many chemical processes, is detrimental in flow assurance of oil and natural gas pipelines. A primary difficulty in the investigation of gas hydrate nucleation has been the inability of researchers to compare nucleation rates of gas hydrates across various systems of different scales and complexities, which in turn has been limiting the ability of researchers to study the nucleation process itself. In this study, a first-generation high-pressure automated lag time apparatus (HP-ALTA MkI) was used to determine the nucleation curve of structure I (sI) - forming carbon dioxide hydrate. The instrument subjected a quiescent water sample of well-defined dimensions to a large number of linear cooling ramps under isobaric conditions, and detected and recorded carbon dioxide hydrate formation temperature distributions. A survival curve was constructed from the measured ensemble, and a nucleation curve was derived from the survival curve using the empirical model-independent method we had previously reported. The nucleation rate of carbon dioxide hydrate was found to be significantly greater than that of pure methane hydrate or that of natural gas hydrate over the entire range of subcooling investigated. We provide a new physical interpretation of an experimentally determined nucleation curve and, by doing so, solve one of the outstanding puzzles of the HP-ALTA technology.

Phase stability of the ice XVII-based CO2 chiral hydrate from molecular dynamics simulations

We computed the phase diagram of CO₂ hydrates at high pressure (HP), from 0.3 to 20 kbar, by means of molecular dynamics simulations. The two CO₂ hydrates known to occur in this pressure range are the cubic structure I (sI) clathrate and the HP hydrate, whose water framework is the recently discovered ice XVII. We investigated the stability of both hydrates upon heating (melting) as well as the phase changes upon compression. The CO₂-filled ice XVII is found to be more stable than the sI clathrate and than the mixture of ice VI and dry ice at pressure values ranging from 6 to 18 kbar and in a wide temperature range, although a phenomenological correction suggests that the stability should more realistically range from 6.5 to 13.5 kbar. Our simulation results support the current hypothesis that the HP hydrate is stable at temperatures above the melting curve of ice VI.

Observation of methane filled hexagonal ice stable up to 150 GPa

Gas clathrates in which water cages enclose the guest gas molecules naturally form on our planet, for instance on ocean floors. Related methane hydrates are also expected to be present under very high pressures (10 to 200 GPa) in giant planetary interiors such as Uranus or Neptune. However, the stability of such structures at these pressures is currently debated. Joint Raman spectroscopy and ab initio simulations show the stability of a high-pressure methane hydrate phase up to at least 150 GPa, the highest pressure explored to date for such compounds. The structure of this phase and the complex transition mechanism from the known MH-III phase are detailed and shown to be in accordance with all known experimental observations.

Gas hydrates consist of hydrogen-bonded water frameworks enclosing guest gas molecules and have been the focus of intense research for almost 40 y, both for their fundamental role in the understanding of hydrophobic interactions and for gas storage and energy-related applications. The stable structure of methane hydrate above 2 GPa, where CH₄ molecules are located within H₂O or D₂O channels, is referred to as methane hydrate III (MH-III). The stability limit of MH-III and the existence of a new high-pressure phase above 40 to 50 GPa, although recently conjectured, remain unsolved to date. We report evidence for a further high-pressure, room-temperature phase of the CH₄–D₂O hydrate, based on Raman spectroscopy in diamond anvil cell and ab initio molecular dynamics simulations including nuclear quantum effects. Our results reveal that a methane hydrate IV (MH-IV) structure, where the D₂O network is isomorphic with ice Ih, forms at ∼40 GPa and remains stable up to 150 GPa at least. Our proposed MH-IV structure is fully consistent with previous unresolved X-ray diffraction patterns at 55 GPa [T. Tanaka et al., J. Chem. Phys. 139, 104701 (2013)]. The MH-III → MH-IV transition mechanism, as suggested by the simulations, is complex. The MH-IV structure, where methane molecules intercalate the tetrahedral network of hexagonal ice, represents the highest-pressure gas hydrate known up to now. Repulsive interactions between methane and water dominate at the very high pressure probed here and the tetrahedral topology outperforms other possible arrangements in terms of space filling.

An ultralow-density porous ice with the largest internal cavity identified in the water phase diagram

The recent back-to-back findings of low-density porous ice XVI and XVII have rekindled the century-old field of the solid-state physics and chemistry of water. Experimentally, both ice XVI and XVII crystals can be produced by extracting guest atoms or molecules enclosed in the cavities of preformed ice clathrate hydrates. Herein, we examine more than 200 hypothetical low-density porous ices whose structures were generated according to a database of zeolite structures. Hitherto unreported porous EMT ice, named according to zeolite nomenclature, is identified to have an extremely low density of 0.5 g/cm3 and the largest internal cavity (7.88 Å in average radius). The EMT ice can be viewed as dumbbell-shaped motifs in a hexagonal close-packed structure. Our first-principles computations and molecular dynamics simulations confirm that the EMT ice is stable under negative pressures and exhibits higher thermal stability than other ultralow-density ices. If all cavities are fully occupied by hydrogen molecules, the EMT ice hydrate can easily outperform the record hydrogen storage capacity of 5.3 wt % achieved with sII hydrogen hydrate. Most importantly, in the reconstructed temperature-pressure (T-P) phase diagram of water, the EMT ice is located at deeply negative pressure regions below ice XVI and at higher temperature regions next to FAU. Last, the phonon spectra of empty-sII, FAU, EMT, and other zeolite-like ice structures are computed by using the dispersion corrected vdW-DF2 functional. Compared with those of ice XI (0.93 g/cm3), both the bending and stretching vibrational modes of the EMT ice are blue-shifted due to their weaker hydrogen bonds.

Salt- and gas-filled ices under planetary conditions

In recent years, evidence has emerged that solid water can contain substantial amounts of guest species, such as small gas molecules-in gas hydrate structures-or ions-in salty ice structures-and that these 'filled' ice structures can be stable under pressures of tens of Gigapascals and temperatures of hundreds of Kelvins. The inclusion of guest species can strongly modify the density, vibrational, diffusive and conductivity properties of ice under high pressure, and promote novel exotic properties. In this review, we discuss our experimental findings and molecular dynamics simulation results on the structures formed by salt- and gas-filled ices, their unusual properties, and the unexpected dynamical phenomena observed under pressure and temperature conditions relevant for planetary interiors modelling. This article is part of the theme issue 'The physics and chemistry of ice: scaffolding across scales, from the viability of life to the formation of planets'.

Nanosecond X-ray diffraction of shock-compressed superionic water ice

Since Bridgman's discovery of five solid water (H₂O) ice phases¹ in 1912, studies on the extraordinary polymorphism of H₂O have documented more than seventeen crystalline and several amorphous ice structures^2,3, as well as rich metastability and kinetic effects^4,5. This unique behaviour is due in part to the geometrical frustration of the weak intermolecular hydrogen bonds and the sizeable quantum motion of the light hydrogen ions (protons). Particularly intriguing is the prediction that H₂O becomes superionic^6-12-with liquid-like protons diffusing through the solid lattice of oxygen-when subjected to extreme pressures exceeding 100 gigapascals and high temperatures above 2,000 kelvin. Numerical simulations suggest that the characteristic diffusion of the protons through the empty sites of the oxygen solid lattice (1) gives rise to a surprisingly high ionic conductivity above 100 Siemens per centimetre, that is, almost as high as typical metallic (electronic) conductivity, (2) greatly increases the ice melting temperature^7-13 to several thousand kelvin, and (3) favours new ice structures with a close-packed oxygen lattice^13-15. Because confining such hot and dense H₂O in the laboratory is extremely challenging, experimental data are scarce. Recent optical measurements along the Hugoniot curve (locus of shock states) of water ice VII showed evidence of superionic conduction and thermodynamic signatures for melting¹⁶, but did not confirm the microscopic structure of superionic ice. Here we use laser-driven shockwaves to simultaneously compress and heat liquid water samples to 100-400 gigapascals and 2,000-3,000 kelvin. In situ X-ray diffraction measurements show that under these conditions, water solidifies within a few nanoseconds into nanometre-sized ice grains that exhibit unambiguous evidence for the crystalline oxygen lattice of superionic water ice. The X-ray diffraction data also allow us to document the compressibility of ice at these extreme conditions and a temperature- and pressure-induced phase transformation from a body-centred-cubic ice phase (probably ice X) to a novel face-centred-cubic, superionic ice phase, which we name ice XVIII^2,17.

Room temperature electrofreezing of water yields a missing dense ice phase in the phase diagram

Water can freeze into diverse ice polymorphs depending on the external conditions such as temperature (T) and pressure (P). Herein, molecular dynamics simulations show evidence of a high-density orthorhombic phase, termed ice χ, forming spontaneously from liquid water at room temperature under high-pressure and high external electric field. Using free-energy computations based on the Einstein molecule approach, we show that ice χ is an additional phase introduced to the state-of-the-art T–P phase diagram. The χ phase is the most stable structure in the high-pressure/low-temperature region, located between ice II and ice VI, and next to ice V exhibiting two triple points at 6.06 kbar/131.23 K and 9.45 kbar/144.24 K, respectively. A possible explanation for the missing ice phase in the T–P phase diagram is that ice χ is a rare polarized ferroelectric phase, whose nucleation/growth occurs only under very high electric fields.

Water can crystallize in different ice polymorphs according to temperature and pressure conditions. Here the authors predict by molecular dynamics simulations a new ice phase spontaneously forming at room temperature under high pressure and high electric field.

Energetics of ice nucleation in mesoporous titania using positron annihilation spectroscopy

The low temperature behavior of water and kinetics of ice nucleation in titania mesopores have been probed by positron annihilation lifetime spectroscopy as a function of pore filling. It is revealed that water undergoes complete freezing at around 220 K when more than 50% of the pore volume is filled and such freezing is hindered at lower hydration levels. A model describing progressive trapping of positronium by ice nuclei in liquid water during the phase transition is employed to estimate the energy associated with the nucleation under confinement. It is observed that the energy for ice nucleation in confinement is less than the activation energy for nucleation in bulk water because of the surface assisted nucleation inside the pore. Interestingly, energy for nucleation is seen to decrease with the lowering of hydration level and ascribed to the curtailed hydrogen bonding network of water at lower pore filling.