Neutron diffraction study of the structure of deuterated ice VIII to 10 GPa

Halogen molecular modifications at high pressure: the case of iodine

Metallization and dissociation are key transformations in diatomic molecules at high densities particularly significant for modeling giant planets. Using X-ray absorption spectroscopy and atomistic modeling, we demonstrate that in halogens, the formation of a connected molecular structure takes place at pressures well below metallization. Here we show that the iodine diatomic molecule first elongates by ∼0.007 Å up to a critical pressure of Pc ∼ 7 GPa, developing bonds between molecules. Then its length continuously decreases with pressure up to 15-20 GPa. Universal trends in halogens are shown and allow us to predict for chlorine a pressure of 42 ± 8 GPa for molecular bond-length reversal. Our findings contribute to tackling the molecule invariability paradigm in diatomic molecular phases at high pressures and may be generalized to other abundant diatomic molecules in the universe, including hydrogen.

Temperature-dependent kinetic pathways featuring distinctive thermal-activation mechanisms in structural evolution of ice VII

Ice amorphization, low- to high-density amorphous (LDA-HDA) transition, as well as (re)crystallization in ice, under compression have been studied extensively due to their fundamental importance in materials science and polyamorphism. However, the nature of the multiple-step "reverse" transformation from metastable high-pressure ice to the stable crystalline form under reduced pressure is not well understood. Here, we characterize the rate and temperature dependence of the structural evolution from ice VII to ice I recovered at low pressure (∼5 mTorr) using in situ time-resolved X-ray diffraction. Unlike previously reported ice VII (or ice VIII)→LDA→ice I transitions, we reveal three temperature-dependent successive transformations: conversion of ice VII into HDA, followed by HDA-to-LDA transition, and then crystallization of LDA into ice I. Significantly, the temperature-dependent characteristic times indicate distinctive thermal activation mechanisms above and below 110-115 K for both ice VIII-to-HDA and HDA-to-LDA transitions. Large-scale molecular-dynamics calculations show that the structural evolution from HDA to LDA is continuous and involves substantial movements of the water molecules at the nanoscale. The results provide a perspective on the interrelationship of polyamorphism and unravel its underpinning complexities in shaping ice-transition kinetic pathways.

The role of water in Earth's mantle

Geophysical observations suggest that the transition zone is wet locally. Continental and oceanic sediment components together with the basaltic and peridotitic components might be transported and accumulated in the transition zone. Low-velocity anomalies at the upper mantle–transition zone boundary might be caused by the existence of dense hydrous magmas. Water can be carried farther into the lower mantle by the slabs. The anomalous Q and shear wave regions locating at the uppermost part of the lower mantle could be caused by the existence of fluid or wet magmas in this region because of the water-solubility contrast between the minerals in the transition zone and those in the lower mantle. δ-H solid solution AlO₂H–MgSiO₄H₂ carries water into the lower mantle. Hydrogen-bond symmetrization exists in high-pressure hydrous phases and thus they are stable at the high pressures of the lower mantle. Thus, the δ-H solid solution in subducting slabs carries water farther into the bottom of the lower mantle. Pyrite FeO₂H_x is formed due to a reaction between the core and hydrated slabs. This phase could be a candidate for the anomalous regions at the core–mantle boundary.

Direct observation of symmetrization of hydrogen bond in δ-AlOOH under mantle conditions using neutron diffraction

At ambient pressure, the hydrogen bond in materials such as ice, hydrates, and hydrous minerals that compose the Earth and icy planets generally takes an asymmetric O-H···O configuration. Pressure significantly affects this configuration, and it is predicted to become symmetric, such that the hydrogen is centered between the two oxygen atoms at high pressure. Changes of physical properties of minerals relevant to this symmetrization have been found; however, the atomic configuration around this symmetrization has remained elusive so far. Here we observed the pressure response of the hydrogen bonds in the aluminous hydrous minerals δ-AlOOH and δ-AlOOD by means of a neutron diffraction experiment. We find that the transition from P2₁nm to Pnnm at 9.0 GPa, accompanied by a change in the axial ratios of δ-AlOOH, corresponds to the disorder of hydrogen bond between two equivalent sites across the center of the O···O line. Symmetrization of the hydrogen bond is observed at 18.1 GPa, which is considerably higher than the disorder pressure. Moreover, there is a significant isotope effect on hydrogen bond geometry and transition pressure. This study indicates that disorder of the hydrogen bond as a precursor of symmetrization may also play an important role in determining the physical properties of minerals such as bulk modulus and seismic wave velocities in the Earth's mantle.

When immiscible becomes miscible-Methane in water at high pressures

At low pressures, the solubility of gases in liquids is governed by Henry's law, which states that the saturated solubility of a gas in a liquid is proportional to the partial pressure of the gas. As the pressure increases, most gases depart from this ideal behavior in a sublinear fashion, leveling off at pressures in the 1- to 5-kbar (0.1 to 0.5 GPa) range with solubilities of less than 1 mole percent (mol %). This contrasts strikingly with the well-known marked increase in solubility of simple gases in water at high temperature associated with the critical point (647 K and 212 bar). The solubility of the smallest hydrocarbon, the simple gas methane, in water under a range of pressure and temperature is of widespread importance, because it is a paradigmatic hydrophobe and occurs widely in terrestrial and extraterrestrial geology. We report measurements up to 3.5 GPa of the pressure dependence of the solubility of methane in water at 100°C-well below the latter's critical temperature. Our results reveal a marked increase in solubility between 1 and 2 GPa, leading to a state above 2 GPa where the maximum solubility of methane in water exceeds 35 mol %.

New Raman measurements for H2O ice VII in the range of 300 cm-1 to 4000 cm-1 at pressures up to 120 GPa

Raman spectroscopic measurements for H₂O ice VII have been conducted to 120 GPa at 300 K in the spectroscopic range of 300-4000 cm^-1. Both moissanite and diamond anvils were used for the experiments. This overcomes the problems of overlapping spectra between the diamond anvil and sample, which had prevented the observation of the stretching modes at pressures higher than ∼23 GPa in all previous measurements. The new results reveal many bands which have not been reported before. The pressure dependences of the Raman modes show anomalous changes at 13-15, ∼27, ∼44, ∼60, and 90 GPa, implying possible structural changes at these pressures. The new results demonstrate that the predicted symmetric hydrogen bond phase X transition does not occur below 120 GPa.

Hydrogen bonding and delocalization in the ELF analysis approach

Delocalization in hydrogen bonded systems is revealed by ELF along the interaction line.

Brillouin scattering of H2O ice to megabar pressures

The sound velocity in polycrystalline ice was measured as a function of pressure at room temperature to 100 GPa, through the phase field of ice VII and crossing the ice X transition, by Brillouin scattering in order to examine the elasticity, compression mechanism, and structural transitions in this pressure range. In particular, we focused on previously proposed phase transitions below 60 GPa. Throughout this pressure range, we find no evidence for anomalous changes in compressibility, and the sound velocities and elastic moduli do not exhibit measurable discontinuous shifts with pressure. Subtle changes in the pressure dependence of the bulk modulus at intermediate pressures can be attributed to high shear stresses at these compressions. The C(11) and C(12) moduli are consistent with previously reported results to 40 GPa and increase monotonically at higher pressures.

High-pressure crystallography

The ability of pressure to change inter-atomic distances strongly leads to a wide range of pressure-induced phenomena at high pressures: for example metallisation, amorphisation, superconductivity and polymerisation. Key to understanding these phenomena is the determination of the crystal structure using x-ray or neutron diffraction. The tools necessary to compress matter above 1 million atmospheres (1 Megabar or 100 GPa) were established by the mid 1970s, but it is only since the early 1990s that we have been able to determine the detailed crystal structures of materials at such pressures. In this chapter I briefly review the history of high-pressure crystallography, and describe the techniques used to obtain and study materials at high pressure. Recent crystallographic studies of elements are then used to illustrate what is now possible using modern detectors and synchrotron sources. Finally, I speculate as to what crystallographic studies might become possible over the next decade.

Hydrogen bond topology and the ice VII/VIII and Ih/XI proton ordering phase transitions

Ice Ih, ordinary ice at atmospheric pressure, is a proton-disordered crystal that when cooled under special conditions is believed to transform to ferroelectric proton-ordered ice XI, but this transformation is still subject to controversy. Ice VII, also proton disordered throughout its region of stability, transforms to proton-ordered ice VIII upon cooling. In contrast to the ice Ih/XI transition, the VII/VIII transition and the crystal structure of ice VIII are well characterized. In order to shed some light on the ice Ih proton ordering transition, we present the results of periodic electronic density functional theory calculations and statistical simulations. We are able to describe the small energy differences among the innumerable H-bond configurations possible in a large simulation cell by using an analytic theory to extrapolate from electronic DFT calculations on small unit cells to cells large enough to approximate the thermodynamic limit. We first validate our methods by comparing our predictions to the well-characterized ice VII/VIII proton ordering transition, finding agreement with respect to both the transition temperature and structure of the low-temperature phase. For ice Ih, our results indicate that a proton-ordered phase is attainable at low temperatures, the structure of which is in agreement with the experimentally proposed ferroelectric structure. The predicted transition temperature of is in qualitative agreement with the observed transition at on KOH-doped ice samples.

Computer simulation of the 13 crystalline phases of ice

As a reference for follow-up studies toward more accurate model parametrizations, we performed molecular-dynamics and Monte Carlo simulations for all known crystalline phases of ice, as described by the simple point-charge/extended and TIP4P water models. We started from the measured structures, densities, and temperatures, and carried out classical canonical simulations for all these arrangements. All simulated samples were cooled down close to 0 K to facilitate the comparison with theoretical estimates. We determined configurational internal energies as well as pressures, and monitored how accurately the measured configurations were preserved during the simulations. While these two models predicted very similar thermophysical and structural properties for water at ambient conditions, the predicted features for the corresponding ice polymorphs may differ significantly.

Anharmonic Raman spectra in high-pressure ice from ab initio simulations

We calculate from ab initio molecular dynamics the Raman scattering of high-pressure ice. To this effect we apply a new method based on the Berry phase theory of polarization. Our results are in agreement with recent and difficult experiments and are compatible with a picture in which ice VII is a proton-disordered system and in ice X the hydrogen bond is symmetric.

Compression of Ice to 210 Gigapascals: Infrared Evidence for a Symmetric Hydrogen-Bonded Phase

Protonated and deuterated ices (H2O and D2O) compressed to a maximum pressure of 210 gigapascals at 85 to 300 kelvin exhibit a phase transition at 60 gigapascals in H2O ice (70 gigapascals in D2O ice) on the basis of their infrared reflectance spectra determined with synchrotron radiation. The transition is characterized by soft-mode behavior of the nu3 O-H or O-D stretch below the transition, followed by a hardening (positive pressure shift) above it. This behavior is interpreted as the transformation of ice phase VII to a structure with symmetric hydrogen bonds. The spectroscopic features of the phase persisted to the maximum pressures (210 gigapascals) of the measurements, although changes in vibrational mode coupling were observed at 150 to 160 gigapascals.

Pressure-Tuned Fermi Resonance in Ice VII

Fermi resonance was observed between the OH stretch and the overtone of the OH bending modes of HDO molecules contaminated in phase VII of D(2)O ice over the pressure range from 17 to 30 gigapascals. An anharmonic coupling constant, which is related to the potential energy surface on which hydrogen-bonded protons oscillate, was found to range around 50 wave numbers through the resonant pressure range. Its experimentally obtained magnitude and pressure-insensitive behavior will be useful for theoretical studies of the potential energy surface and hence of the nature of hydrogen bonding in ice.