Equations of state of ice VI and ice VII at high pressure and high temperature

Metastable water at several compression rates and its freezing kinetics into ice VII

Water can be dynamically over-compressed well into the stability field of ice VII. Whether water then transforms into ice VII, vitreous ice or a metastable novel crystalline phase remained uncertain. We report here the freezing of over-compressed water to ice VII by time-resolved X-ray diffraction. Quasi-isothermal dynamic compression paths are achieved using a dynamic-piezo-Diamond-Anvil-Cell, with programmable pressure rise time from 0.1 ms to 100 ms. By combining the present data set with those obtained on various ns-dynamical platforms, a complete evolution of the solidification pressure of metastable water versus the compression rate is rationalized within the classical nucleation theory framework. Also, when crystallization into ice VII occurs in between 1.6 GPa and 2.0 GPa, that is in the stability field of ice VI, a structural evolution over few ms is then observed into a mixture of ice VI and ice VII that seems to resolve apparent contradictions between previous results.

Oxidation of iron by giant impact and its implication on the formation of reduced atmosphere in the early Earth

Giant impact-driven redox processes in the atmosphere and magma ocean played crucial roles in the evolution of Earth. However, because of the absence of rock records from that time, understanding these processes has proven challenging. Here, we present experimental results that simulate the giant impact-driven reactions between iron and volatiles (H2O and CO2) using x-ray free electron laser (XFEL) as fast heat pump and structural probe. Under XFEL pump, iron is oxidized to wüstite (FeO), while volatiles are reduced to H2 and CO. Furthermore, iron oxidation proceeds into formation of hydrides (γ-FeHx) and siderite (FeCO3), implying redox boundary near 300-km depth. Through quantitative analysis on reaction products, we estimate the volatile and FeO budgets in bulk silicate Earth, supporting the Theia hypothesis. Our findings shed light on the fast and short-lived process that led to reduced atmosphere, required for the emergence of prebiotic organic molecules in the early Earth.

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.

High-temperature and high-pressure plastic phase of ice at the boundary of liquid water and ice VII

Simultaneous high-temperature and high-pressure studies reveal phase transformation of bulk liquid water to an ice-VII-like structure having an eight coordination. It was demonstrated through this numerical study that the observed high-temperature and high-pressure phase of water obtained upon shock compression and equilibration has high rotational diffusion and thereby the hydrogen dynamics of these crystal structures are significantly complex compared with ice VII. The current work provides new characterization methods for the numerically observed plastic crystal phase of ice at the boundary of the liquid water and ice VII phases in which the molecules have a defined lattice position but rotate freely. It is anticipated that the present work will provide important data and guide new theoretical and experimental investigations in the search for plastic crystal phases of water. The power spectra plots of bulk liquid water subjected to different temperature and pressure conditions have also been presented in this numerical study, demonstrating significant differences between these high-temperature and high-pressure shock-equilibrated phases and those of pure ice VII at 10 GPa and liquid water at ambient temperature and pressure, as well as at elevated pressures and temperatures.

First-Order Phase Transformation at Constant Volume: A Continuous Transition?

We describe a first-order phase transition of a simple system in a process where the volume is kept constant. We show that, unlike what happens when the pressure is constant, (i) the transformation extends over a finite temperature (and pressure) range, (ii) each and every extensive potential (internal energy U, enthalpy H, Helmholtz energy F, and Gibbs energy G), and the entropy S is continuous across the transition, and (iii) the constant-volume heat capacity does not diverge during the transition and only exhibits discrete jumps. These non-intuitive results highlight the importance of controlling the correct variables in order to distinguish between continuous and discontinuous transitions. We apply our results to describe the transition between ice VI and liquid water using thermodynamic information available in the literature and also to show that a first-order phase transition driven in isochoric condition can be used as the operating principle of a mechanical actuator.

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.

Immiscibility in N2-H2O solids up to 140 GPa

Nitrogen and water are very abundant in nature; however, the way they chemically react at extreme pressure-temperature conditions is unknown. Below 6 GPa, they have been reported to form clathrate compounds. Here, we present Raman spectroscopy and x-ray diffraction studies in the H₂O-N₂ system at high pressures up to 140 GPa. We find that clathrates, which form locally in our diamond cell experiments above 0.3 GPa, transform into a fine grained state above 6 GPa, while there is no sign of formation of mixed compounds. We point out size effects in fine grained crystallites, which result in peculiar Raman spectra in the molecular regime, but x-ray diffraction shows no additional phase or deviation from the bulk behavior of familiar solid phases. Moreover, we find no sign of ice doping by nitrogen, even in the regimes of stability of nonmolecular nitrogen.

Interferometric measurements of refractive index and dispersion at high pressure

We describe a high precision interferometer system to measure the pressure dependence of the refractive index and its dispersion in the diamond anvil cell (DAC). The reflective Fabry-Perot fringe patterns created by both a white light and a monochromatic beam are recorded to determine both the sample thickness and its index at the laser wavelength and to characterize the dispersion in the visible range. Advances in sample preparation, optical setup, and data analysis enable us to achieve [Formula: see text] random uncertainty, demonstrated with an air sample, a factor of five improvement over the best previous DAC measurement. New data on [Formula: see text] liquid water and ice VI up to 2.21 GPa at room temperature illustrate how higher precision measurements of the index and its optical dispersion open up new opportunities to reveal subtle changes in the electronic structure of water at high pressure.

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.

High-pressure preference for reduced water content in porous zinc aspartate hydrates

The zinc aspartate (ZnAsp2) complex, a common dietary supplement, preferentially crystallizes as the dihydrate (ZnAsp2·2H2O) from aqueous solution. Under normal conditions the dihydrate easily transforms into the sesquihydrate (ZnAsp2·1.5H2O). The dihydrate crystal structure is triclinic, space group P1, and the sesquihydrate is monoclinic, space group C2/c. However, their structures are closely related and similarly consist of zinc aspartate ribbons parallel to pores accommodating water molecules. These porous structures can breathe water molecules in and out depending on the temperature and air humidity. High pressure above 50 MPa favours the sesquihydrate, as shown by recrystallizations under pressure and compressibility measured by single-crystal X-ray diffraction up to 4 GPa. This preference is explained by the reduced volume of the sesquihydrate and water compressed separately, compared with the dihydrate. The sesquihydrate undergoes an isostructural phase transition when the voids collapse at 0.8 GPa, whereas no phase transitions occur in the dihydrate, because its pores are supported by increased water content.

Local-Basis-Function Equation of State for Ice VII–X to 450 GPa at 300 K

Helmholtz energy of ice VII–X is determined in a pressure regime extending to 450 GPa at 300 K using local-basis-functions in the form of b-splines. The new representation for the equation of state is embedded in a physics-based inverse theory framework of parameter estimation. Selected pressures as a function of volume from 14 prior experimental studies and two theoretical studies constrain the behavior of Helmholtz energy. Separately measured bulk moduli, not used to construct the representation, are accurately replicated below about 20 GPa and above 60 GPa. In the intermediate range of pressure, the experimentally determined moduli are larger and have greater scatter than values predicted using the Helmholtz representation. Although systematic error in the determination of elastic moduli is possible and likely, the alternative hypothesis is a slow relaxation time associated with changes in proton mobility or the ice VII to X transition. A correlation is observed between anomalies in the pressure derivative of the predicted bulk modulus and previously suggested higher-order phase transitions. Improved determinations of elastic properties at high pressure would allow refinement of the current equation of state. More generally, the current method of data assimilation is broadly applicable to other materials in high-pressure studies and for investigations of planetary interiors.

Shock growth of ice crystal near equilibrium melting pressure under dynamic compression

Crystal growth and morphological transitions are crucial for fundamental science and wide applications. Nevertheless, their mechanisms under local nonequilibrium growth condition are unclear due to severe interference of thermal and mass transports on the interplay between thermodynamic driving force and interface kinetics. Here, we reveal the origin of the pressure-induced 2D shock growth of ice VI crystal by using dynamic compression, in which a dimensional transition from 3D to 2D is observed. Unlike generally expected, the 2D shock growth occurs from 3D crystal edges rather than from its corners upon fast compression, even near equilibrium growth condition. This is due to similar interface structure to the crystal edge plane facilitating the fast interface kinetics under local nonequilibrium growth.

Crystal growth is governed by an interplay between macroscopic driving force and microscopic interface kinetics at the crystal–liquid interface. Unlike the local equilibrium growth condition, the interplay becomes blurred under local nonequilibrium, which raises many questions about the nature of diverse crystal growth and morphological transitions. Here, we systematically control the growth condition from local equilibrium to local nonequilibrium by using an advanced dynamic diamond anvil cell (dDAC) and generate anomalously fast growth of ice VI phase with a morphological transition from three- to two-dimension (3D to 2D), which is called a shock crystal growth. Unlike expected, the shock growth occurs from the edges of 3D crystal along the (112) crystal plane rather than its corners, which implies that the fast compression yields effectively large overpressure at the crystal–liquid interface, manifesting the local nonequilibrium condition. Molecular dynamics (MD) simulation reproduces the faster growth of the (112) plane than other planes upon applying large overpressure. Moreover, the MD study reveals that the 2D shock crystal growth originates from the similarity of the interface structure between water and the (112) crystal plane under the large overpressure. This study provides insight into crystal growth under dynamic compressions, which makes a bridge for the unknown behaviors of crystal growth between under static and dynamic pressure conditions.

Rapid freezing of water under dynamic compression

Understanding the behavior of materials at extreme pressures is a central issue in fields like aerodynamics, astronomy, and geology, as well as for advancing technological grand challenges such as inertial confinement fusion. Dynamic compression experiments to probe high-pressure states often encounter rapid phase transitions that may cause the materials to behave in unexpected ways, and understanding the kinetics of these phase transitions remains an area of great interest. In this review, we examine experimental and theoretical/computational efforts to study the freezing kinetics of water to a high-pressure solid phase known as ice VII. We first present a detailed analysis of dynamic compression experiments in which water has been observed to freeze on sub-microsecond time scales to ice VII. This is followed by a discussion of the limitations of currently available molecular and continuum simulation methods in modeling these experiments. We then describe how our phase transition kinetics models, which are based on classical nucleation theory, provide a more physics-based framework that overcomes some of these limitations. Finally, we give suggestions on future experimental and modeling work on the liquid-ice VII transition, including an outline of the development of a predictive multiscale model in which molecular and continuum simulations are intimately coupled.

Compression Freezing Kinetics of Water to Ice VII

Time-resolved x-ray diffraction (XRD) of compressed liquid water shows transformation to ice VII in 6 nsec, revealing crystallization rather than amorphous solidification during compression freezing. Application of classical nucleation theory indicates heterogeneous nucleation and one-dimensional (e.g., needlelike) growth. These first XRD data demonstrate rapid growth kinetics of ice VII with implications for fundamental physics of diffusion-mediated crystallization and thermodynamic modeling of collision or impact events on ice-rich planetary bodies.

Melting dynamics of ice in the mesoscopic regime

How does a crystal melt? How long does it take for melt nuclei to grow? The melting mechanisms have been addressed by several theoretical and experimental works, covering a subnanosecond time window with sample sizes of tens of nanometers and thus suitable to determine the onset of the process but unable to unveil the following dynamics. On the other hand, macroscopic observations of phase transitions, with millisecond or longer time resolution, account for processes occurring at surfaces and time limited by thermal contact with the environment. Here, we fill the gap between these two extremes, investigating the melting of ice in the entire mesoscopic regime. A bulk ice I _h or ice VI sample is homogeneously heated by a picosecond infrared pulse, which delivers all of the energy necessary for complete melting. The evolution of melt/ice interfaces thereafter is monitored by Mie scattering with nanosecond resolution, for all of the time needed for the sample to reequilibrate. The growth of the liquid domains, over distances of micrometers, takes hundreds of nanoseconds, a time orders of magnitude larger than expected from simple H-bond dynamics.

Isostructural crystal hydrates of rare-earth metal oxalates at high pressure: from strain anisotropy to dehydration

The crystal structures of a series of isostructural rare-earth metal oxalates, (REE)2(C2O4)3·10H2O (REE=Sm, Y) and a 1:1 YSm(C2O4)3·10H2O solid solution, have been studied in situ by single-crystal X-ray diffraction and optical microscopy. The structures were followed from ambient pressure to 6 GPa in a DAC with paraffin as the hydrostatic fluid. Bulk compressibilities, anisotropic lattice strain on hydrostatic compression and the corresponding changes in the atomic coordinates were followed. Discontinuities/sharp changes in the slopes of the pressure dependences of volume and selected cell parameters have been observed for yttrium-containing salts at ~3.5 GPa. This may be related to the re-distribution of water molecules within the crystal structure. Y2(C2O4)3·10H2O undergoes a partial dehydration at 1 GPa, forming monoclinic Y2(C2O4)3·6H2O as single-crystalline inclusions in the original phase.

Probing ice VII crystallization from amorphous NaCl–D2O solutions at gigapascal pressures

The high density amorphous solution NaCl·10.2D₂O crystallises at 260 K as almost pure ice VII during annealing at gigapascal pressures.

Structural characterization of eutectic aqueous NaCl solutions under variable temperature and pressure conditions

The structure of amorphous NaCl solutions produced by fast quenching is studied as a function of pressure, up to 4 GPa, by combined neutron diffraction experiments and classical molecular dynamics simulations. Similarly to LiCl solutions the system amorphizes at ambient pressure in a dense phase structurally similar to the e-HDA phase in pure water. The measurement of the static structure factor as a function of pressure allowed us to validate a new polarizable force field developed by Tazi et al., 2012, never tested under non-ambient conditions. We infer from simulations that the hydration shells of Na(+) cations form well defined octahedra composed of both H2O molecules and Cl(-) anions at low pressure. These octahedra are gradually broken by the seventh neighbour moving into the shell of first neighbours yielding an irregular geometry. In contrast to LiCl solutions and pure water, the system does not show a polyamorphic transition under pressure. This confirms that the existence of polyamorphism relies on the tetrahedral structure of water molecules, which is broken here.