Structural characterization of ice XIX as the second polymorph related to ice VI

Direct Visualization of Molecular Stacking in Quasi-2D Hexagonal Ice

Understanding ice nucleation and growth is of great interest to researchers due to its importance in the biological, cryopreservation, and environmental fields. However, microstructural investigations of ice on the molecular scale are still lacking. In this paper, a simple method is proposed to prepare quasi-2-dimensional ice Ih films, which have been characterized via cryogenic transmission electron microscope. The intersecting stacking faults of basal (BSF) and prismatic (PSF) types have been directly visualized and resolved with a notable first-time report of PSF in ice Ih. Moreover, the possible growth pathways of BSF, namely, the Ic phase, were elucidated by the theoretical calculations and the chair conformation of H2O molecules. This study offers valuable insights that can enhance researchers' understanding of the growth kinetics of crystalline ice.

Kinetic isotope effects on hydrogen/deuterium disordering and ordering in ice crystals: A Raman and dielectric study of ice VI, XV, and XIX

Ice XIX and ice XV are both partly hydrogen-ordered counterparts to disordered ice VI. The ice XIX → XV transition represents the only order-to-order transition in ice physics. Using Raman and dielectric spectroscopies, we investigate the ambient-pressure kinetics of the two individual steps in this transition in real time (of hours), that is, ice XIX → transient ice VI (the latter called VI‡) and ice VI‡ → ice XV. Hydrogen-disordered ice VI‡ appears intermittent between 101 and 120 K, as inferred from the appearance and subsequent disappearance of the ice VI Raman marker bands. A comparison of the rate constants for the H2O ices reported here with those from D2O samples [Thoeny et al., J. Chem. Phys. 156, 154507 (2022)] reveals a large kinetic isotope effect for the ice XIX decay, but a much smaller one for the ice XV buildup. An enhancement of the classical overbarrier rate through quantum tunneling for the former can provide a possible explanation for this finding. The activation barriers for both transitions are in the 18-24 kJ/mol range, which corresponds to the energy required to break a single hydrogen bond. These barriers do not show an H/D isotope effect and are the same, no matter whether they are derived from Raman scattering or from dielectric spectroscopy. These findings favor the notion that a dipolar reorientation, involving the breakage of a hydrogen bond, is the rate determining step at the order-to-order transition.

Configurational entropy of ice XIX and its isotope effect

Ice XIX is a partly hydrogen-ordered polymorph related to disordered ice VI, similar to ice XV. We here investigate the order-order-disorder sequence ice XIX→ice XV→ice VI based on calorimetry at ambient pressure both for D2O and H2O-ice XIX. From these data we extract configurational entropy differences between ice XIX, ice XV and ice VI. This task is complex because, unlike for all other ices, the order-disorder transition from ice XIX to ice VI takes place in two steps via ice XV. Even more challenging, these two steps take place in an overlapping manner, so that careful separation of slow kinetics is necessary. This is evidenced best by changing the heating rate in calorimetry experiments: For fast heating experiments the second step, disordering of ice XV, is suppressed because the first step, formation of ice XV from ice XIX, is too slow. The transient state ice VI‡ that is initially produced upon ice XIX decay then does not have enough time to convert to ice XV, but remains disordered all along. In order to tackle the challenge to determine the entropy difference between ice XIX and VI as well as the entropy difference between ice XV and VI we employ two different approaches that allow assessing the impact of kinetics on the entropy change. "Single peak integration" defines a kinetically limited result, but "combined peak integration" allows estimation of the true thermodynamic values. Our best estimate for the true value shows ice XIX to be much more ordered than ice XV (25 ± 3% vs 9 ± 4% of the Pauling entropy). For D2Oice XIX samples we obtain 28% of order, but only when a small number of fast H-isotope defects are used. In the second part we use these results to estimate the location of the ice XIX phase boundary both for protiated and deuterated ice XIX. The initial Clapeyron slope at ambient pressure is determined from the combination of neutron powder diffraction volume differences and calorimetry entropy differences data to be 21 K GPa-1 with an order-disorder transition temperature To-d(0.0 GPa) = 103 ± 1 K. An in situ bracketing experiment at 1.8 GPa yields To-d(1.8 GPa) = 116 ± 3 K, i.e., the phase boundary slope flattens at higher pressures. These data allow us to determine the region of thermodynamic stability of ice XIX in the phase diagram and to explain the surprising isotope shift reversal at 1.6 GPa compared to 0.0 GPa, i.e., why D2O-ice XIX disorders at lower temperatures than H2O-ice XIX at 1.6 GPa, but at higher temperatures at ambient pressures.

Formation of a two-dimensional helical square tube ice in hydrophobic nanoslit using the TIP5P water model

In extreme and nanoconfinement conditions, the tetrahedral arrangement of water molecules is challenged, resulting in a rich and new phase behavior unseen in bulk phases. The unique phase behavior of water confined in hydrophobic nanoslits has been previously observed, such as the formation of a variety of two-dimensional (2D) ices below the freezing temperature. The primary identified 2D ice phase, termed square tube ice (STI), represents a unique arrangement of water molecules in 2D ice, which can be viewed as an array of 1D ice nanotubes stacked in the direction parallel to the confinement plane. In this study, we report the molecular dynamics (MD) simulations evidence of a novel 2D ice phase, namely, helical square tube ice (H-STI). H-STI is characterized by the stacking of helical ice nanotubes in the direction parallel to the confinement plane. Its structural specificity is evident in the presence of helical square ice nanotubes, a configuration unseen in both STI and single-walled ice nanotubes. A detailed analysis of the hydrogen bonding strength showed that H-STI is a 2D ice phase diverging from the Bernal-Fowler-Pauling ice rules by forming only two strong hydrogen bonds between adjacent molecules along its helical ice chain. This arrangement of strong hydrogen bonds along ice nanotube and weak bonds between the ice nanotube shows a similarity to quasi-one-dimensional van der Waals materials. Ab initio molecular dynamics simulations (over a 30 ps) were employed to further verify H-STI's stability at 1 GPa and temperature up to 200 K.

GenIce-core: Efficient algorithm for generation of hydrogen-disordered ice structures

Ice is different from ordinary crystals because it contains randomness, which means that statistical treatment based on ensemble averaging is essential. Ice structures are constrained by topological rules known as the ice rules, which give them unique anomalous properties. These properties become more apparent when the system size is large. For this reason, there is a need to produce a large number of sufficiently large crystals that are homogeneously random and satisfy the ice rules. We have developed an algorithm to quickly generate ice structures containing ions and defects. This algorithm is provided as an independent software module that can be incorporated into crystal structure generation software. By doing so, it becomes possible to simulate ice crystals on a previously impossible scale.

Thermodynamically Stable Intermediate in the Course of Hydrogen Ordering from Ice V to Ice XIII

Even though many partially ordered ices are known, it remains elusive to understand and categorize them. In this study, we study the ordering from ice V to XIII using calorimetry at ambient pressure and discover that the transition takes place via an intermediate that is thermodynamically stable at 113-120 K. Our isothermal ordering approach allows us to highlight the distinction of this intermediate from ice V and XIII, where there are clear differences both in terms of enthalpy and ordering kinetics. We suggest that the approach developed in the present work can also reveal the nature of partially ordered forms in the hydrogen order-disorder series of other ice phases.

Close-Packed Ices in Nanopores

Water molecules in any of the ice polymorphs organize themselves into a perfect four-coordinated hydrogen-bond network at the expense of dense packing. Even at high pressures, there seems to be no way to reconcile the ice rules with the close packing. Here, we report several close-packed ice phases in carbon nanotubes obtained from molecular dynamics simulations of two different water models. Typically they are in plastic states at high temperatures and are transformed into the hydrogen-ordered ice, keeping their close-packed structures at lower temperatures. The close-packed structures of water molecules in carbon nanotubes are identified with those of spheres in a cylinder. We present design principles of hydrogen-ordered, close-packed structures of ice in nanotubes, which suggest many possible dense ice forms with or without nonzero polarization. In fact, some of the simulated ices are found to exhibit ferroelectric ordering upon cooling.

Slightly Hydrogen-Ordered State of Ice IV Evidenced by In Situ Neutron Diffraction

Ice IV is a metastable high-pressure phase of ice in which the water molecules exhibit orientational disorder. Although orientational ordering is commonly observed for other ice phases, it has not been reported for ice IV. We conducted in situ powder neutron diffraction experiments for DCl-doped D2O ice IV to investigate its hydrogen ordering. We found abrupt changes in the temperature derivative of unit-cell volume, dV/dT, at ∼120 K, and revealed a slightly ordered structure at low temperatures based on the Rietveld method. The occupancy of the D1 site deviates from 0.5 in particular; it increased when samples were cooled at higher pressures and reached 0.174(14) at 2.38 GPa, 58 K. Our results evidence the presence of a low-symmetry hydrogen-ordered state corresponding to ice IV. It seems, however, difficult to experimentally access the completely ordered phase corresponding to ice IV by slow cooling at high pressure.

Simultaneous measurements of volume, pressure, optical images, and crystal structure with a dynamic diamond anvil cell: A real-time event monitoring system

The dynamic diamond anvil cell (dDAC) technique has attracted great interest because it possibly provides a bridge between static and dynamic compression studies with fast, repeatable, and controllable compression rates. The dDAC can be a particularly useful tool to study the pathways and kinetics of phase transitions under dynamic pressurization if simultaneous measurements of physical quantities are possible as a function of time. We here report the development of a real-time event monitoring (RTEM) system with dDAC, which can simultaneously record the volume, pressure, optical image, and structure of materials during dynamic compression runs. In particular, the volume measurement using both Fabry-Pérot interferogram and optical images facilitates the construction of an equation of state (EoS) using the dDAC in a home-laboratory. We also developed an in-line ruby pressure measurement (IRPM) system to be deployed at a synchrotron x-ray facility. This system provides simultaneous measurements of pressure and x-ray diffraction in low and narrow pressure ranges. The EoSs of ice VI obtained from the RTEM and the x-ray diffraction data with the IRPM are consistent with each other. The complementarity of both RTEM and IRPM systems will provide a great opportunity to scrutinize the detailed kinetic pathways of phase transitions using dDAC.

The hydrogen-bond network in sodium chloride tridecahydrate: analogy with ice VI

The structure of a recently found hyperhydrated form of sodium chloride (NaCl·13H2O and NaCl·13D2O) has been determined by in situ single-crystal neutron diffraction at 1.7 GPa and 298 K. It has large hydrogen-bond networks and some water molecules have distorted bonding features such as bifurcated hydrogen bonds and five-coordinated water molecules. The hydrogen-bond network has similarities to ice VI in terms of network topology and disordered hydrogen bonds. Assuming the equivalence of network components connected by pseudo-symmetries, the overall network structure of this hydrate can be expressed by breaking it down into smaller structural units which correspond to the ice VI network structure. This hydrogen-bond network contains orientational disorder of water molecules in contrast to the known salt hydrates. An example is presented here for further insights into a hydrogen-bond network containing ionic species.

Theoretical Prediction of the Anti-Icing Activity of Two-Dimensional Ice I

Two-dimensional (2D) ice I is atomic-level ice that is composed of two interlocked atomic layers saturated with hydrogen bonds. It has recently been experimentally observed, but its properties have yet to be clarified. Accordingly, we theoretically studied the hydrophobic properties of 2D ice I. On the contrary, a simulation of a hydrogen fluoride molecule on a 2D ice surface manifested that it destroyed the 2D ice structure and connected new hydrogen bonds with water molecules. Investigations of the interfacial effect between 2D and three-dimensional (3D) ice films indicated that the network structure of 2D ice was not destroyed by a 3D ice surface, as the former was saturated with hydrogen bonds. However, the surface of 3D ice reorganized to form as many hydrogen bonds as possible. Thus, the 2D ice film was hydrophobic and inhibited the growth of 3D ice. This shows that if 2D ice can be produced on an industrial scale, it can be used as an anti-3D-icing agent under low temperatures.

Melting curves of ice polymorphs in the vicinity of the liquid-liquid critical point

The possible existence of a liquid-liquid critical point in deeply supercooled water has been a subject of debate due to the challenges associated with providing definitive experimental evidence. The pioneering work by Mishima and Stanley [Nature 392, 164-168 (1998)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected liquid-liquid transition and its associated critical point. Based on the continuous or discontinuous changes in the slope of the melting curves, Mishima [Phys. Rev. Lett. 85, 334 (2000)] suggested that the liquid-liquid critical point lies between the melting curves of ice III and ice V. We explore this conjecture using molecular dynamics simulations with a machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII and find that all of them are supercritical and do not intersect the liquid-liquid transition locus. We also find a pronounced, yet continuous, change in the slope of the melting lines upon crossing of the liquid locus of maximum compressibility. Finally, we analyze the literature in light of our findings and conclude that the scenario in which the melting curves are supercritical is favored by the most recent computational and experimental evidence. Although the preponderance of evidence is consistent with the existence of a second critical point in water, the behavior of ice polymorph melting lines does not provide strong evidence in support of this viewpoint, according to our calculations.

Medium-density amorphous ice

Amorphous ices govern a range of cosmological processes and are potentially key materials for explaining the anomalies of liquid water. A substantial density gap between low-density and high-density amorphous ice with liquid water in the middle is a cornerstone of our current understanding of water. However, we show that ball milling "ordinary" ice Ih at low temperature gives a structurally distinct medium-density amorphous ice (MDA) within this density gap. These results raise the possibility that MDA is the true glassy state of liquid water or alternatively a heavily sheared crystalline state. Notably, the compression of MDA at low temperature leads to a sharp increase of its recrystallization enthalpy, highlighting that H₂O can be a high-energy geophysical material.

Mateverse, the Future Materials Science Computation Platform Based on Metaverse

Currently, computational materials science involves human-computer interaction through coding in software or neural networks. There is still no direct way for human intelligence endorsement. The digitalization of human intelligence should be the ultimate goal for many disciplines. In materials science, human intelligence is still irreplaceable from machine learning techniques, where humans can deal with complex correlations in the real world. We design the framework of Mateverse, a materials science computation platform based on Metaverse, which unifies human intelligence, experiment data, and theoretical simulations. In Mateverse, we intensively study the properties of H₂O, including the liquid and solid phases. We show that we can optimize a new water force field (which we name TIP4P-Meta) directly from the interactions between human and visible properties of H₂O. This force field is validated to be better than the conventional water model, and new ice polymorphs can be generated. We believe our platform can provide valuable hints in the paradigm upgrade in future computational materials science development.

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.

Infrared Analysis and Effect of Nitrogen and Nitrous Oxide on the Glass Transition of Methanol Cryofilms

Methanol plays an important role in studying the structure and dynamics of hydrogen bonds in alcohols and other physiologically important compounds in the condensed state. The physical vapor deposition method and two-beam interferometry make it possible to study the structure of polyatomic molecules, the nature and character of intermolecular interactions, and the internal structure of various compounds. Thus, it becomes possible to analyze changes in the internal structure of substances near the points of their phase transitions and glass transitions at ultralow temperatures. In this work, we have studied the effect of nitrogen and nitrous oxide on the dependence of refractive indices of methanol on temperature and on structural-phase transformations upon heating from 16 to 160 K.

Water adsorption and dynamics on graphene and other 2D materials: Computational and experimental advances

The interaction of water and surfaces, at molecular level, is of critical importance for understanding processes such as corrosion, friction, catalysis and mass transport. The significant literature on interactions with single crystal metal surfaces should not obscure unknowns in the unique behaviour of ice and the complex relationships between adsorption, diffusion and long-range inter-molecular interactions. Even less is known about the atomic-scale behaviour of water on novel, non-metallic interfaces, in particular on graphene and other 2D materials. In this manuscript, we review recent progress in the characterisation of water adsorption on 2D materials, with a focus on the nano-material graphene and graphitic nanostructures; materials which are of paramount importance for separation technologies, electrochemistry and catalysis, to name a few. The adsorption of water on graphene has also become one of the benchmark systems for modern computational methods, in particular dispersion-corrected density functional theory (DFT). We then review recent experimental and theoretical advances in studying the single-molecular motion of water at surfaces, with a special emphasis on scattering approaches as they allow an unparalleled window of observation to water surface motion, including diffusion, vibration and self-assembly.

Optical Studies of Thin Films of Cryocondensed Mixtures of Water and Admixture of Nitrogen and Argon

The interaction of host molecules with water molecules is of primary importance in astrophysical and atmospheric studies. Water-binding interactions continue to attract a broad interest in various fields, especially those related to the formation of assembly structures. Using the physical vapor deposition (PVD) method and a two-beam interferometer with a wavelength of 406 nm, the refractive indices of thin films of a water and nitrogen (argon) mixture were calculated in the range from 15 to 35 K. The results of temperature transformations of the obtained films from a two-beam interferometer, and thermal desorption characteristics from the temperature of condensation to the temperature of evaporation of water (15-180 K), are presented. The relationship between the signal of the interferometer, the refractive index, and the film thickness during glass transition is demonstrated.

Accurate crystal structure of ice VI from X-ray diffraction with Hirshfeld atom refinement

Water is an essential chemical compound for living organisms, and twenty of its different crystal solid forms (ices) are known. Still, there are many fundamental problems with these structures such as establishing the correct positions and thermal motions of hydrogen atoms. The list of ice structures is not yet complete as DFT calculations have suggested the existence of additional and - to date - unknown phases. In many ice structures, neither neutron diffraction nor DFT calculations nor X-ray diffraction methods can easily solve the problem of hydrogen atom disorder or accurately determine their anisotropic displacement parameters (ADPs). Here, accurate crystal structures of H2O, D2O and mixed (50%H2O/50%D2O) ice VI obtained by Hirshfeld atom refinement (HAR) of high-pressure single-crystal synchrotron and laboratory X-ray diffraction data are presented. It was possible to obtain O-H/D bond lengths and ADPs for disordered hydrogen atoms which are in good agreement with the corresponding single-crystal neutron diffraction data. These results show that HAR combined with X-ray diffraction can compete with neutron diffraction in detailed studies of polymorphic forms of ice and crystals of other hydrogen-rich compounds. As neutron diffraction is relatively expensive, requires larger crystals which can be difficult to obtain and access to neutron facilities is restricted, cheaper and more accessible X-ray measurements combined with HAR can facilitate the verification of the existing ice polymorphs and the quest for new ones.

Machine Learning Methods for Multiscale Physics and Urban Engineering Problems

We present an overview of four challenging research areas in multiscale physics and engineering as well as four data science topics that may be developed for addressing these challenges. We focus on multiscale spatiotemporal problems in light of the importance of understanding the accompanying scientific processes and engineering ideas, where "multiscale" refers to concurrent, non-trivial and coupled models over scales separated by orders of magnitude in either space, time, energy, momenta, or any other relevant parameter. Specifically, we consider problems where the data may be obtained at various resolutions; analyzing such data and constructing coupled models led to open research questions in various applications of data science. Numeric studies are reported for one of the data science techniques discussed here for illustration, namely, on approximate Bayesian computations.

The Emergence of 2D MXenes Based Zn-Ion Batteries: Recent Development and Prospects

Rechargeable zinc-ion batteries (ZIBs) with exceptional theoretical capacity have garnered significant interest in large-scale electrochemical energy storage devices due to their low cost, abundant material, inherent safety, high specific energy, and ecofriendly nature. Metal carbides/nitrides, known as MXenes, have emerged as a large family of 2D transition metal carbides or carbonitrides with excellent properties, e.g., high electrical conductivity, large surface functional groups (e.g., F, O, and OH), low energy barriers for the diffusion of electrolyte ions with wide interlayer spaces. After a decade of effort, significant development has been achieved in the synthesis, properties, and applications of MXenes. Thus, it has opened up various exciting opportunities to construct advanced MXene-based nanostructures for ZIBs with excellent specific energy and power. Herein, this review summarizes the advances across multiple synthesis routes, related properties, morphological and structural characteristics, and chemistries of MXenes for ZIBs. The recent development of MXene-based electrodes is introduced, and electrolytes for ZIBs are elucidated in detail. MXene-based rocking chair ZIBs, strategies to enhance the performance of MXene-based cathodes, suppress the dendrites in MXene-based anodes, and MXene-based flexible ZIBs are pointed out. A rational design and modification of the MXenes as well as the production of composites with metal oxides exhibits promise in solving issues and enhancing the electrochemical performance of ZIBs. Finally, the present challenges and future prospects for MXene-based ZIBs are discussed.

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.

Calorimetric Investigation of Hydrogen-Atom Sublattice Transitions in the Ice VI/XV/XIX Trio

Ice XIX represents the latest discovery of ice polymorphs and exists in the medium pressure range near 1–2 GPa. Ice XIX is a partially hydrogen-ordered phase, by contrast to its disordered mother phase ice VI, which shares the same oxygen-atom network with ice XIX. Ice XIX differs in terms of the ordering of the hydrogen-atom sublattice, and hence the space group, from its hydrogen-ordered sibling ice XV, which also features the same type of oxygen network. Together, ice VI, XV, and XIX form the only known trio of ice polymorphs, where polymorphic transformations from order to order, order to disorder, and disorder to order are possible, which also compete with each other depending on the thermodynamic path taken and the cooling/heating rates employed. These transitions in the H-sublattice have barely been investigated, so we study here the unique triangular relation in the ice VI/XV/XIX trio based on calorimetry experiments. We reveal the following key features for H-sublattice transitions: (i) upon cooling ice VI, domains of ice XV and XIX develop simultaneously, where pure ice XV forms at ≤0.85 GPa and pure ice XIX forms at ≥1.60 GPa, (ii) ice XIX transforms into ice XV via a transient disordered state, (iii) ice XV recooled at ambient pressure features a complex domain structure, possibly containing an unknown H-ordered polymorph, (iv) recooled ice XV partly transforms back into ice XIX at 1.80 GPa, and (v) partial deuteration slows down domain reordering strongly. These findings not only are of interest in understanding possible hydrogen-ordering and -disordering processes in the interior of icy moons and planets but, more importantly, also provide a challenging benchmark for our understanding and parameterizing many-body interactions in H-bonded networks.

Metastability of Liquid Water Freezing into Ice VII under Dynamic Compression

The ubiquitous nature and unusual properties of water have motivated many studies on its metastability under temperature- or pressure-induced phase transformations. Here, nanosecond compression by a high-power laser is used to create the nonequilibrium conditions where liquid water persists well into the stable region of ice VII. Through our experiments, as well as a complementary theoretical-computational analysis based on classical nucleation theory, we report that the metastability limit of liquid water under nearly isentropic compression from ambient conditions is at least 8 GPa, higher than the 7 GPa previously reported for lower loading rates.

Advances in Atomic Force Microscopy: Imaging of Two- and Three-Dimensional Interfacial Water

Interfacial water is closely related to many core scientific and technological issues, covering a broad range of fields, such as material science, geochemistry, electrochemistry and biology. The understanding of the structure and dynamics of interfacial water is the basis of dealing with a series of issues in science and technology. In recent years, atomic force microscopy (AFM) with ultrahigh resolution has become a very powerful option for the understanding of the complex structural and dynamic properties of interfacial water on solid surfaces. In this perspective, we provide an overview of the application of AFM in the study of two dimensional (2D) or three dimensional (3D) interfacial water, and present the prospect and challenges of the AFM-related techniques in experiments and simulations, in order to gain a better understanding of the physicochemical properties of interfacial water.

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.

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.

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.

Structure and nature of ice XIX

Ice is a material of fundamental importance for a wide range of scientific disciplines including physics, chemistry, and biology, as well as space and materials science. A well-known feature of its phase diagram is that high-temperature phases of ice with orientational disorder of the hydrogen-bonded water molecules undergo phase transitions to their ordered counterparts upon cooling. Here, we present an example where this trend is broken. Instead, hydrochloric-acid-doped ice VI undergoes an alternative type of phase transition upon cooling at high pressure as the orientationally disordered ice remains disordered but undergoes structural distortions. As seen with in-situ neutron diffraction, the resulting phase of ice, ice XIX, forms through a Pbcn-type distortion which includes the tilting and squishing of hexameric clusters. This type of phase transition may provide an explanation for previously observed ferroelectric signatures in dielectric spectroscopy of ice VI and could be relevant for other icy materials.

Dynamics & Spectroscopy with Neutrons-Recent Developments & Emerging Opportunities

This work provides an up-to-date overview of recent developments in neutron spectroscopic techniques and associated computational tools to interrogate the structural properties and dynamical behavior of complex and disordered materials, with a focus on those of a soft and polymeric nature. These have and continue to pave the way for new scientific opportunities simply thought unthinkable not so long ago, and have particularly benefited from advances in high-resolution, broadband techniques spanning energy transfers from the meV to the eV. Topical areas include the identification and robust assignment of low-energy modes underpinning functionality in soft solids and supramolecular frameworks, or the quantification in the laboratory of hitherto unexplored nuclear quantum effects dictating thermodynamic properties. In addition to novel classes of materials, we also discuss recent discoveries around water and its phase diagram, which continue to surprise us. All throughout, emphasis is placed on linking these ongoing and exciting experimental and computational developments to specific scientific questions in the context of the discovery of new materials for sustainable technologies.

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.