Communication: Hypothetical ultralow-density ice polymorphs

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.

Ice structures assembled from cubic water clusters of D2d and S4 symmetry

The study of self-assembly processes is of key importance for fundamental science and modern technologies. Cubic water clusters of D2d and S4 symmetry show great potential as building blocks for self-assembly. The objective of this paper is to construct possible ice structures formed by hydrogen bonding of these very stable water clusters. A number of such structures are herein presented, including quasi-2D and 3D ices as well as spatial layered and tubular ices. The energetics and structure of many configurations differing in the arrangement of hydrogen atoms in hydrogen bonds have been studied. It was established that the proton disorder of all such ices is of island type. The residual entropy of these ices is equal to ln(3)/4 in dimensionless form. For layered structures formed by the stacking of multiple bilayers, the determining role of the van der Waals interactions is shown. Note that, for all considered ices, the lowest-energy configurations are formed only by clusters of D2d symmetry.

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.

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.

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.

Novel Algorithm to Generate Hydrogen-Disordered Ice Structures

We propose an efficient algorithm for generating hydrogen-disordered ice networks utilizing graph theory and the topological characteristic of the network. The computational efficiency with the new algorithm is much higher than the conventional ones developed by Rahman and Stillinger and Buch et al. The difference in the computational time between our algorithm and either of the two conventional ones increases with increasing the system size.

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.

Hydrogen Inter-Cage Hopping and Cage Occupancies inside Hydrogen Hydrate: Molecular-Dynamics Analysis

The inter-cage hopping in a type II clathrate hydrate with different numbers of H2 and D2 molecules, from 1 to 4 molecules per large cage, was studied using a classical molecular dynamics simulation at temperatures of 80 to 240 K. We present the results for the diffusion of these guest molecules (H2 or D2) at all of the different occupations and temperatures, and we also calculated the activation energy as the energy barrier for the diffusion using the Arrhenius equation. The average occupancy number over the simulation time showed that the structures with double and triple large-cage H2 occupancy appeared to be the most stable, while the small cages remained with only one guest molecule. A Markov model was also calculated based on the number of transitions between the different cage types.

Computational Prediction of Novel Ice Phases: A Perspective

Although computational prediction of new ice phases is a niche field in water science, the scientific subject itself is representative of two important areas in physical chemistry, namely, statistical thermodynamics and molecular simulations. The prediction of a variety of novel ice phases has also attracted general public interest since the 1980s. In particular, the prediction of low-dimensional ice phases has gained momentum since the confirmation of a number of low-dimensional "computer ice" phases in the laboratory over the past decade. In this Perspective, the research advancements in computational prediction of novel ice phases over the past few years are reviewed. Particular attention is placed on new ice phases whose physical properties or dimensional structures are distinctly different from conventional bulk ices. Specific topics include the (i) formation of superionic ices, (ii) electrofreezing of water under high pressure and in a high external electric field, (iii) prediction of low-density porous ice at strongly negative pressure, (iv) ab initio computational study of two-dimensional (2D) ice under nanoscale confinement, and (v) 2D ices formed on a solid surface near ambient temperature without nanoscale confinement. Clearly, the formation of most of these novel ice phases demands certain extreme conditions. Ongoing challenges and new opportunities for predicting new ice phases from either classical molecular dynamics simulation or high-level ab initio computation are discussed.

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.

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.

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.

Advances in the experimental exploration of water's phase diagram

Water's phase diagram displays enormous complexity with currently 17 experimentally confirmed polymorphs of ice and several more predicted computationally. For almost 120 years, it has been a stomping ground for scientific discovery, and ice research has often been a trailblazer for investigations into a wide range of materials-related phenomena. Here, the experimental progress of the last couple of years is reviewed, and open questions as well as future challenges are discussed. The specific topics include (i) the polytypism and stacking disorder of ice I, (ii) the mechanism of the pressure amorphization of ice I, (iii) the emptying of gas-filled clathrate hydrates to give new low-density ice polymorphs, (iv) the effects of acid/base doping on hydrogen-ordering phase transitions as well as (v) the formation of solid solutions between salts and the ice polymorphs, and the effect this has on the appearance of the phase diagram. In addition to continuing efforts to push the boundaries in terms of the extremes of pressure and temperature, the exploration of the "chemical" dimensions of ice research appears to now be a newly emerging trend. It is without question that ice research has entered a very exciting era.

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.

Adsorption of Kinetic Hydrate Inhibitors on Growing Surfaces: A Molecular Dynamics Study

We investigate the mechanism of a typical kinetic hydrate inhibitor (KHI), polyvinylcaprolactam (PVCap), which has been applied to prevent hydrate plugs from forming in gas pipe lines, using molecular dynamics simulations of crystal growth of ethylene oxide hydrate. Water-soluble ethylene oxide is chosen as a guest species to avoid problems associated with the presence of the gas phase in the simulation cell such as slow crystal growth. A PVCap dodecamer adsorbs irreversibly on the hydrate surface which grows at supercooling of 3 K when the hydrophobic part of two pendent groups are trapped in open cages at the surface. The amide hydrogen bonds make no contribution to the adsorption. PVCap can adsorb on various crystallographic planes of sI hydrate. This is in contrast to antifreeze proteins, each of which prefers a specific plane of ice. The trapped PVCap gives rise to necessarily the concave surface of the hydrate. The crystal growth rate decreases with increasing surface curvature, indicating that the inhibition by PVCap is explained by the Gibbs-Thomson effect.

Assessing Intermolecular Interactions in Guest-Free Clathrate Hydrate Systems

Recently, empty hydrate structures sI, sII, sH, and others have been proposed as low-density ice structures by both experimental observations and computer simulations. Some of them have been synthesized in the laboratory, which motivates further investigations on the stability of such guest-free clathrate structures. Using semiempirical and ab initio-based water models, as well as dispersion-corrected density functional theory approaches, we predict their stability, including cooperative many-body effects, in comparison with reference data from converged wave function-based DF-MP2 electronic structure calculations. We show that large basis sets and counterpoise corrections are required to improve convergence in the interaction/binding energies for such systems. Therefore, extrapolation schemes based on triple/quadruple and quadruple/quintuple ζ quality basis sets are used to reach high accuracy. Eleven different water structures corresponding to dodecahedron, edge sharing, face sharing, and fused cubes, as a part of the WATER27 database, as well as cavities from the sI, sII, and sH clathrate hydrates formed by 20, 24, 28, and 36 water molecules, are employed, and new benchmark energies are reported. Using these benchmark sets of interaction energies, we assess the performance of both analytical models and direct DFT calculations for such clathrate-like systems. In particular, seven popular water models (TIP4P/ice, TIP4P/2005, q-TIP4P/F, TTM2-F, TTM3-F, TTM4-F, and MB-pol) available in the literature, and nine density functional approximations (3 meta-GGAs, 3 hybrids, and 3 range separated functionals) are used to investigate their accuracy. By including dispersion corrections, our results show that errors in the interaction energies are reduced for most of the DFT functionals. Despite the difficulties faced by current water models and DFT functionals to accurately describe the interactions in such water systems, we found some general trends that could serve to extend their applicability to larger systems.