skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Hydrogen bond symmetrization and high-spin to low-spin transition of ε-FeOOH at the pressure of Earth’s lower mantle
Abstract We focus on the ferric end-member of phase H: ε-FeOOH using density functional theory at the PBEsol+U level. At 300 K, we find that ε-FeOOH undergoes a hydrogen bond symmetrization at 37 GPa and a sharp high-spin to low-spin transition at 45 GPa. We find excellent agreement with experimental measurements of the equation of state, lattice parameters, atomic positions, vibrational frequencies, and optical properties as related to the band gap, which we find to be finite and small, decreasing with pressure. The hydrogen bond symmetrization transition is neither first-nor second-order, with no discontinuity in volume or any of the elastic moduli. Computed IR and Raman frequencies and intensities show that vibrational spectroscopy may provide the best opportunity for locating the hydrogen bond symmetrization transition experimentally. We find that ε-FeOOH is highly anisotropic in both longitudinal- and shear-wave velocities at all pressures, with the shear wave velocity varying with propagation and polarization direction by as much as 24% at zero pressure and 43% at 46 GPa. The shear and bulk elastic moduli increase by 18% across the high-spin to low-spin transition.  more » « less
Award ID(s):
2223935
PAR ID:
10535919
Author(s) / Creator(s):
;
Publisher / Repository:
Mineralogical Society of America
Date Published:
Journal Name:
American Mineralogist
Volume:
108
Issue:
12
ISSN:
0003-004X
Page Range / eLocation ID:
2209 to 2218
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract Constraining the accommodation, distribution, and circulation of hydrogen in the Earth's interior is vital to our broader understanding of the deep Earth due to the significant influence of hydrogen on the material and rheological properties of minerals. Recently, a great deal of attention has been paid to the high-pressure polymorphs of FeOOH (space groups P21nm and Pnnm). These structures potentially form a hydrogen-bearing solid solution with AlOOH and phase H (MgSiO4H2) that may transport water (OH–) deep into the Earth's lower mantle. Additionally, the pyrite-type polymorph (space group Pa3 of FeOOH), and its potential dehydration have been linked to phenomena as diverse as the introduction of hydrogen into the outer core (Nishi et al. 2017), the formation of ultralow-velocity zones (ULVZs) (Liu et al. 2017), and the Great Oxidation Event (Hu et al. 2016). In this study, the high-pressure evolution of FeOOH was re-evaluated up to ~75 GPa using a combination of synchrotron-based X-ray diffraction (XRD), Fourier transform infrared spectroscopy (FTIR), and optical absorption spectroscopy. Based on these measurements, we report three principal findings: (1) pressure-induced changes in hydrogen bonding (proton disordering or hydrogen bond symmetrization) occur at substantially lower pressures in ε-FeOOH than previously reported and are unlikely to be linked to the high-spin to low-spin transition; (2) ε-FeOOH undergoes a 10% volume collapse coincident with an isostructural Pnnm → Pnnm transition at approximately 45 GPa; and (3) a pressure-induced band gap reduction is observed in FeOOH at pressures consistent with the previously reported spin transition (40 to 50 GPa). 
    more » « less
  2. Abstract The 3.65 Å phase [MgSi(OH)6] is a hydrous phase that is predicted to be stable in a simplified MgO-SiO2-H2O (MSH) ternary system at pressures exceeding 9 GPa. Along cold subduction zones, it is likely to transport water, bound in its crystalline lattice, into the Earth’s interior. The 3.65 Å phase consists of Mg and Si octahedral sites attached to the hydroxyl group that forms a hydrogen bond and is predicted to undergo pressure-induced symmetrization of the hydrogen bond. Therefore, in this study, we investigate the high-pressure behavior of the 3.65 Å phase using Raman spectroscopy. We have conducted five distinct compressions up to ~60 GPa using two different pressure-transmitting media—alcohol mixture and neon. At ambient conditions, we identified vibrational modes using complementary first-principles simulations based on density functional perturbation theory. Upon compression, we note that the first derivative of the vibrational modes in the lattice region stiffens, i.e., b1lattice > 0. In contrast, the hydroxyl region softens, i.e., b1OH > 0. This is indicative of the strengthening of hydrogen bonding upon compression. We noticed a significant broadening of vibrational modes related to hydroxyl groups that are indicative of proton disorder. However, within the maximum pressures explored in this study, we did not find evidence for pressure-induced symmetrization of the hydrogen bonds. We used the pressure derivative of the vibrational modes to determine the ratio of the bulk moduli and their pressure derivative. We note that the smaller bulk moduli of hydrous phases compared to the major mantle phases are compensated by significantly larger pressure derivatives of the bulk moduli for the hydrous phases. This leads to a significant reduction in the elasticity contrast between hydrous and major mantle phases. Consequently, the detection of the degree of mantle hydration is likely to be challenging at greater depths. 
    more » « less
  3. null (Ed.)
    Abstract The high-pressure phases of oxyhydroxides (δ-AlOOH, ε-FeOOH, and their solid solution), candidate components of subducted slabs, have wide stability fields, thus potentially influencing volatile circulation and dynamics in the Earth’s lower mantle. Here, we report the elastic wave velocities of δ-(Al,Fe)OOH (Fe/(Al + Fe) = 0.13, δ-Fe13) to 79 GPa, determined by nuclear resonant inelastic X-ray scattering. At pressures below 20 GPa, a softening of the phonon spectra is observed. With increasing pressure up to the Fe 3+ spin crossover (~ 45 GPa), the Debye sound velocity ( v D ) increases. At higher pressures, the low spin δ-Fe13 is characterized by a pressure-invariant v D . Using the equation of state for the same sample, the shear-, compressional-, and bulk-velocities ( v S , v P , and v Φ ) are calculated and extrapolated to deep mantle conditions. The obtained velocity data show that δ-(Al,Fe)OOH may cause low- v Φ and low- v P anomalies in the shallow lower mantle. At deeper depths, we find that this hydrous phase reproduces the anti-correlation between v S and v Φ reported for the large low seismic velocity provinces, thus serving as a potential seismic signature of hydrous circulation in the lower mantle. 
    more » « less
  4. Abstract The post-stishovite transition is a classic pseudo-proper typed ferroelastic transition with a symmetry-breaking spontaneous strain. This transition has been studied using high-pressure spontaneous strains, optic modes, and elastic moduli (Cij) based on the Landau modeling, but its atomistic information and structural distortion remain poorly understood. Here we have conducted synchrotron single-crystal X-ray diffraction measurements on stishovite crystals up to 75.3 GPa in a diamond-anvil cell. Analysis of the data reveals atomic positions, bond lengths, bond angles, and variations of SiO6 octahedra across the transition at high pressure. Our results show that the O coordinates split at ~51.4 GPa, where the apical and equatorial Si-O bond lengths cross over, the SiO6 octahedral distortion vanishes, and the SiO6 octahedra start to rotate about the c axis. Moreover, distortion mode analysis shows that an in-plane stretching distortion (GM1+ mode) occurs in the stishovite structure at high pressure while a rotational distortion (GM2+ mode) becomes dominant in the post-stishovite structure. These results are used to correlate with elastic moduli and Landau parameters (symmetry-breaking strain e1–e2 and order parameter Q) to provide atomistic insight into the ferroelastic transition. When the bond lengths of two Si-O bonds are equal due to the contribution from the GM1+ stretching mode, C11 converges with C12, and the shear wave VS1[110] polarizing along [110] and propagating along [110] vanishes. Values of e1–e2 and Q are proportional to the SiO6 rotation angle from the occurrence of the GM1+ rotational mode in the post-stishovite structure. Our results on the pseudo-proper type transition are also compared with that for the proper type in albite and improper type in CaSiO3 perovskite. The symmetry-breaking strain, in all these types of transitions, arises as the primary effect from the structural angle (such as SiO6 rotation or lattice constant angle) and its relevant distortion mode in the low-symmetry ferroelastic phase. 
    more » « less
  5. Using first-principles calculations, this study evaluates the structure, equation of state, and elasticity of three compositions of phase D up to 75 GPa: (1) the magnesium endmember [MgSi2O4(OH)2], (2) the aluminum endmember [Al2SiO4(OH)2], and (3) phase D with 50% Al-substitution [AlMg0.5Si1.5O4(OH)2]. We find that the Mg-endmember undergoes hydrogen-bond symmetrization and that this symmetrization is linked to a 22% increase in the bulk modulus of phase D, in agreement with previous studies. Al2SiO4(OH)2 also undergoes hydrogen-bond symmetrization, but the concomitant increase in bulk modulus is only 13%—a significant departure from the 22% increase of the Mg-endmember. Additionally, Al-endmember phase D is denser (2%–6%), less compressible (6%–25%), and has faster compressional (6%–12%) and shear velocities (12%–15%) relative to its Mg-endmember counterpart. Finally, we investigated the properties of phase D with 50% Al-substitution [AlMg0.5Si1.5O4(OH)2], and found that the hydrogen-bond symmetrization, equation of state parameters, and elastic constants of this tie-line composition cannot be accurately modeled by interpolating the properties of the Mg- and Al-endmembers. 
    more » « less