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
Calculated Elasticity of Al-Bearing Phase D
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
- Award ID(s):
- 1725673
- PAR ID:
- 10413113
- Date Published:
- Journal Name:
- Minerals
- Volume:
- 12
- Issue:
- 8
- ISSN:
- 2075-163X
- Page Range / eLocation ID:
- 922
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The transport of hydrogen into Earth's deep interior may have an impact on lower mantle dynamics as well as on the seismic signature of subducted material. Due to the stability of the hydrous phasesδ‐AlOOH (delta phase), MgSiO2(OH)2(phase H), andε‐FeOOH at high temperatures and pressures, their solid solutions may transport significant amounts of hydrogen as deep as the core‐mantle boundary. We have constrained the equation of state, including the effects of a spin crossover in the Fe3+atoms, of (Al, Fe)‐phase H: Al0.84Fe3+0.07Mg0.02Si0.06OOH, using powder X‐ray diffraction measurements to 125 GPa, supported by synchrotron Mössbauer spectroscopy measurements on (Al, Fe)‐phase H andδ‐(Al, Fe)OOH. The changes in spin state of Fe3+in (Al, Fe)‐phase H results in a significant decrease in bulk sound velocity and occurs over a different pressure range (48–62 GPa) compared withδ‐(Al, Fe)OOH (32–40 GPa). Changes in axial compressibilities indicate a decrease in the compressibility of hydrogen bonds in (Al, Fe)‐phase H near 30 GPa, which may be associated with hydrogen bond symmetrization. The formation of (Al, Fe)‐phase H in subducted oceanic crust may contribute to scattering of seismic waves in the mid‐lower mantle (∼1,100–1,550 km). Accumulation of 1–4 wt.% (Al, Fe)‐phase H could reproduce some of the seismic signatures of large, low seismic‐velocity provinces. Our results suggest that changes in the electronic structure of phases in the (δ‐AlOOH)‐(MgSiO2(OH)2)‐(ε‐FeOOH) solid solution are sensitive to composition and that the presence of these phases in subducted oceanic crust could be seismically detectable throughout the lower mantle.more » « less
-
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
-
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
-
null (Ed.)The degree to which the Earth’s mantle stores and cycles water in excess of the storage capacity of nominally anhydrous minerals is dependent upon the stability of hydrous phases under mantle-relevant pressures, temperatures, and compositions. Two hydrous phases, phase D and phase H, are stable to the pressures and temperatures of the Earth’s lower mantle, suggesting that the Earth’s lower mantle may participate in the cycling of water. We build on our prior work of density functional theory calculations on phase H with the stability, structure, and bonding of hydrous phases D, and we predict the aluminum partitioning with H in the Al 2 O 3 -SiO 2 -MgO-H 2 O system. We address the solid solutions through a statistical sampling of site occupancy and calculation of the partition function from the grand canonical ensemble. We show that each phase has a wide solid solution series between MgSi 2 O 6 H 2 -Al 2 SiO 6 H 2 and MgSiO 4 H 2 -2 δ AlOOH + SiO 2 , in which phase H is more aluminum rich than phase D at a given bulk composition. We predict that the addition of Al to both phases D and H stabilizes each phase to higher temperatures through additional configurational entropy. While we have shown that phase H does not exhibit symmetric hydrogen bonding at high pressure, we report here that phase D undergoes a gradual increase in the number of symmetric H-bonds beginning at ∼30 GPa, and it is only ∼50% complete at 60 GPa.more » « less