Elastic moduli ($C_{ij}$) of single-crystal stishovite and post-stishovite are determined using Brillouin light scattering, impulsive stimulated light scattering, and x-ray diffraction up to 70 GPa. The $C_{12}$ of stishovite converges with the $C_{11}$ at $\sim 55\mathrm{GPa}$, where the transverse wave $V_{S1}$ propagating along [110] also vanishes. Landau modeling of the $C_{ij}$, $B_{1g}$ optic mode, and lattice parameters reveals a pseudoproper type ferroelastic post-stishovite transition. The transition would cause peculiar anomalies in $V_S$ and Poisson’s ratio in silica-bearing subducting slabs in the mid-lower mantle.

Elastic properties of Si _0.95 H _0.21 O ₂ with hydrogarnet substitution (4H ⁺ = Si ⁴⁺ ) across the post‐stishovite transition (28–42 GPa) are determined up to 70 GPa using ab initio calculations and a pseudo‐proper type Landau model. At 28 GPa, elastic coefficients C ₁₁ and C ₁₂ converge, and the average shear and compressional velocity ( V _S and V _P ) decrease by a maximum of 25.5% and 5.2%, respectively. Hydrogarnet substitution reduces ambient elastic moduli and sound velocities, and shifts shear softening to lower pressure. 2%–13% Si _0.95 H _0.21 O ₂ may cause a V _S anomaly of −0.5% to −2.6% at 700–820 km depth, explaining low V _S layers beneath North America and the European Alps. Additionally, 20 vol % SiO ₂ in subducted basalt, with decreasing water content from 3.2 wt% to zero, could cause a V _S anomaly of up to −7(4) % from 700 to 1,900 km depth, aligning with seismic scatterers identified in some subduction regions.

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.

The metastable olivine (Ol) wedge hypothesis assumes that Ol may exist as a metastable phase at the P conditions of the mantle transition zone (MTZ) and even deeper regions due to inhibition of the phase transitions from Ol to wadsleyite and ringwoodite caused by low T in the cold subducting slabs. It is commonly invoked to account for the stagnation of the descending slabs, deep focus earthquakes and other geophysical observations. In the last few years, several new structures with the forsterite (Fo) composition, namely Fo-II, Fo-III and Fo-IV, were either experimentally observed or theoretically predicted at very low T conditions. They may have important impacts on the metastable Ol wedge hypothesis. By performing first-principles calculations, we have systematically examined their crystallographic characteristics, elastic properties and dynamic stabilities from 0 to 100 GPa, and identified the Fo-III phase as the most likely metastable phase to occur in the cold slabs subducted to the depths equivalent to the lower part of the MTZ (below the ~600 km depth) and even the lower mantle. As disclosed by our theoretical simulations, the Fo-III phase is a post-spinel phase (space group Cmc2₁), has all cations in sixfold coordination at P < ~60 GPa, and shows dynamic stability for the entire P range from 0 to 100 GPa. Further, our static enthalpy calculations have suggested that the Fo-III phase may directly form from the Fo material at ~22 GPa (0 K), and our high-T phase relation calculations have located the Fo/Fo-III phase boundary at ~23.75 GPa (room T) with an averaged Clapeyron slope of ~−1.1 MPa/K for the T interval from 300 to 1800 K. All these calculated phase transition pressures are likely overestimated by ~3 GPa because of the GGA method used in this study. The discrepancy between our predicted phase transition P and the experimental observation (~58 GPa at 300 K) can be explained by slow reaction rate and short experimental durations. Taking into account the P-T conditions in the cold downgoing slabs, we therefore propose that the Fo-III phase, rather than the Ol, highly possibly occurs as the metastable phase in the cold slabs subducted to the P conditions of the lower part of the MTZ (below the ~600 km depth) and even the lower mantle. In addition, our calculation has showed that the Fo-III phase has higher bulk seismic velocity, and thus may make important contributions to the high seismic speeds observed in the cold slabs stagnated near the upper mantle-lower mantle boundary. Future seismic studies may discriminate the effects of the Fo-III phase and the low T. Surprisingly, the Fo-III phase will speed up, rather than slow down, the subducting process of the cold slabs, if it metastably forms from the Ol. In general, the Fo-III phase has a higher density than the warm MTZ, but has a lower density than the lower mantle, as suggested by our calculations.