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: Water in the crystal structure of CaSiO3 perovskite
Abstract While the water storage capacities of the upper 700 km depths of the mantle have been constrained by high-pressure experiments and diamond inclusion studies, the storage capacity of the lower mantle remains controversial. A recent high-pressure experimental study on CaSiO3 perovskite, which is the third most abundant mineral in the lower mantle, reported possible storage of H2O up to a few weight percent. However, the substitution mechanism for H in this phase remains unknown. We have conducted a series of density functional theory calculations under static-lattice conditions and high pressures to elucidate hydration mechanisms at the atomic scale. All of the possible dodecahedral (Ca2+ → 2H+) and octahedral (Si4+ → 4H+) substitution configurations for a tetragonal perovskite lattice have very small energy differences, suggesting the coexistence of multiples of H configurations in CaSiO3 perovskite at mantle pressures and temperatures. The dodecahedral substitutions decrease the bulk modulus, resulting in a smaller unit-cell volume of hydrous CaSiO3 perovskite under pressure, consistent with the experimental observations. Although the octahedral substitutions also decrease the bulk modulus, they increase the unit-cell volume at 1 bar. The H atoms substituted in the dodecahedral sites develop much less hydrogen bonding with O atoms, leading to a large distortion in the neighboring SiO6 octahedra. Such distortion may be responsible for the non-cubic peak splittings observed in experiments on hydrous CaSiO3 perovskite. Our calculated infrared spectra suggest that the observed broad OH modes in CaSiO3 perovskite can result from the existence of multiples of H configurations in the phase. Combined with the recent experimental results, our study suggests that CaSiO3 can be an important mineral phase to consider for the H2O storage in the lower mantle.  more » « less
Award ID(s):
2019565
PAR ID:
10335204
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
American Mineralogist
Volume:
107
Issue:
4
ISSN:
0003-004X
Page Range / eLocation ID:
631 to 641
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Phase egg, [AlSiO3(OH)], is an aluminosilicate hydrous mineral that is thermodynamically stable in lithological compositions represented by Al2O3-SiO2-H2O (ASH) ternary, i.e., a simplified ternary for the mineralogy of subducted sediments and continental crustal rocks. High-pressure and high-temperature experiments on lithological compositions resembling hydrated sedimentary layers in subducting slabs show that phase egg is stable up to pressures of 20–30 GPa, which translates to the transition zone to lower mantle depths. Thus, phase egg is a potential candidate for transporting water into the Earth’s mantle transition zone. In this study, we use first-principles simulations based on density functional theory to explore the pressure dependence of crystal structure and how it influences energetics and elasticity. Our results indicate that phase egg exhibits anomalous behavior of the pressure dependence of the elasticity at mantle transition zone depths (~15 GPa). Such anomalous behavior in the elasticity is related to changes in the hydrogen bonding O-H···O configurations, which we delineate as a transition from a low-pressure to a high-pressure structure of phase egg. Full elastic constant tensors indicate that phase egg is very anisotropic resulting in a maximum anisotropy of compressional wave velocity, AvP ≈ 30% and of shear wave velocity, AvS ≈ 17% at zero pressures. Our results also indicate that the phase egg has one of the fastest bulk sound velocities (vP and vS) compared to other hydrous aluminous phases in the ASH ternary, which include topaz-OH, phase Pi, and d-AlOOH. However, the bulk sound velocity of phase egg is slower than that of stishovite. At depths corresponding to the base of mantle transition zone, phase egg decomposes to a mixture of d-AlOOH and stishovite. The changes in compressional DvP and shear DvS velocity associated with the decomposition is ~0.42% and –1.23%, respectively. Although phase egg may be limited to subducted sediments, it could hold several weight percentages of water along a normal mantle geotherm. 
    more » « less
  2. Abstract Davemaoite (CaSiO3 perovskite) is considered the third most abundant phase in the pyrolytic lower mantle and the second most abundant phase in the subducted mid-ocean ridge basalt (MORB). During the partial melting of the pyrolytic upper mantle, incompatible titanium (Ti) becomes enriched in the basaltic magma, forming Ti-rich MORB. Davemaoite is considered an important Ti-bearing mineral in subducted slabs by forming a Ca(Si,Ti)O3 solid solution. However, the crystal structure and compressibility of Ca(Si,Ti)O3 perovskite solid solution at relevant pressure and temperature conditions had not been systematically investigated. In this study, we investigated the structure and equations of state of Ca(Si0.83Ti0.17)O3 and Ca(Si0.75Ti0.25)O3 perovskites at room temperature up to 82 and 64 GPa, respectively, by synchrotron X-ray diffraction (XRD). We found that both Ca(Si0.83Ti0.17)O3 and Ca(Si0.75Ti0.25)O3 perovskites have a tetragonal structure up to the maximum pressures investigated. Based on the observed data and compared to pure CaSiO3 davemaoite, both Ca(Si0.83Ti0.17)O3 and Ca(Si0.75Ti0.25)O3 perovskites are expected to be less dense up to the core-mantle boundary (CMB), and specifically ~1–2% less dense than CaSiO3 davemaoite in the pressure range of the transition zone (15–25 GPa). Our results suggest that the presence of Ti-bearing davemaoite phases may result in a reduction in the average density of the subducting slabs, which in turn promotes their stagnation in the lower mantle. The presence of low-density Ti-bearing davemaoite phases and subduction of MORB in the lower mantle may also explain the seismic heterogeneity in the lower mantle, such as large low shear velocity provinces (LLSVPs). 
    more » « less
  3. Abstract The influence of Al substitution on the elastic properties of stishovite and its transition to post-stishovite is of great importance for interpreting the seismic wave velocities of subducted mid-ocean ridge basalt (MORB) within the mantle transition zone and the lower mantle. However, atomistic mechanisms of Al substitution effects on the transition and its associated elasticity remain debated. Here synchrotron single-crystal X-ray diffraction measurements have been performed at room temperature on Al1.3-SiO2 (1.3 mol% Al in the chemical formula of Si0.965(3)Al0.041(1)O2H0.017(4)) and Al2.1-SiO2 (2.1 mol% Al in Si0.948(2)Al0.064(1)O2H 0.018(3)) crystals in diamond anvil cells with Boehler-Almax designed anvils up to 38.0 GPa and 28.5 GPa, respectively. Refinements of the diffraction patterns show that a transformation from stishovite (space group P42/mnm; No. 136) to CaCl2-typed post-stishovite (space group Pnnm; No. 58) is accompanied by splitting of O coordinates. The Al substitution in stishovite results in a faster decrease in the O coordinate, softer apical (Si,Al)-O bonds, and a softer and less distorted (Si,Al)O6 octahedron under compression. This leads to reduced adiabatic bulk modulus (KS), shear modulus (G), shear wave velocity (VS), and compressional wave velocity (VP) in the stishovite phase, explaining seismic wave perturbations in the mantle transition zone. Together with Raman data, Landau theory modeling shows that Al substitution increases the order parameter and excess free energy, stabilizing the post-stishovite phase at lower pressures. Correlation between elasticity and octahedral distortion index (D) reveals that at certain D, the Al substitution reduces KS, G, VS, and VP of the stishovite phase while increasing G, VS, and VP of the post-stishovite phase. Importantly, the maximum shear reduction is slightly enhanced at D = 0.00620(9) at the transition point. Our results help explain the seismically observed small-scale VS anomalies beneath subduction regions in the shallow lower mantle where Al,H-bearing stishovite undergoes the post-stishovite transition. 
    more » « less
  4. 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
  5. 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