The extrapolation of elastic moduli to high pressure and temperature

The elastic constants of a crystal under stress, defined as the second derivative of the crystal free energy with respect to strain, require a correction related to the static pressure at non-zero pressures. The corrections required for the elastic constants calculated by the free energy minimisation code PARAPOCS are described and tested by comparison with the elastic constants calculated numerically by applying small stresses in the appropriate orientations to simulated crystals of fluorite, forsterite, α-quartz and albite. The corrected elastic constants are then used to investigate the extrapolation of the bulk and shear moduli (and hence also the seismic wave velocities V _p and V _s) of β-spinel and forsterite to upper mantle pressures. A Murnaghan equation, thirdorder Eulerian finite strain equation, second order polynomial equation and a logistic equation were all fitted to the simulated bulk and shear moduli between 0 and 3 GPa pressure. The parameters derived for these equations are used to extrapolate the bulk and shear moduli to 14 GPa and the results are compared to the simulated high pressure moduli. Over this pressure range, the second order polynomial provides the best extrapolation of the bulk modulus, but the use of the logistic equation results in the best extrapolation of the shear modulus.     

Anomalous elasticity of talc at high pressures: Implications for subduction systems

Talc is a layered hydrous silicate mineral that plays a vital role in transporting water into Earth’s interior and is crucial for explaining geophysical observations in subduction zone settings. In this study, we explored the structure, equation of state, and elasticity of both triclinic and monoclinic talc under high pressures up to 18 GPa using first principles simulations based on density functional theory corrected for dispersive forces. Our results indicate that principal components of the full elastic constant tensor C₁₁ and C₂₂, shear components C₆₆, and several off-diagonal components show anomalous pressure dependence. This non-monotonic pressure dependence of elastic constant components is likely related to the structural changes and is often manifested in a polytypic transition from a low-pressure polytype talc-I to a high-pressure polytype talc-II. The polytypic transition of talc occurs at pressures within its thermodynamic stability. However, the bulk and shear elastic moduli show no anomalous softening. Our study also shows that talc has low velocity, extremely high anisotropy, and anomalously high V_P/V_S ratio, thus making it a potential candidate mineral phase that could readily explain unusually high V_P/V_S ratio and large shear wave splitting delays as observed from seismological studies in many subduction systems.     

Structural change of GeO2 under pressure

Pressure-induced amorphization of α-quartz type GeO₂ was studied with a newly developed X-ray diffraction system which consists of a 4-circle goniometer and a curved position sensitive detector. Single-crystal diffraction was measured under pressurs up to 7.3 GPa at room temperature in order to investigate pretransitional phenomena. Diffraction intensity and line width of the diffraction profiles showed no remarkable change up to 5.9 GPa. However, no sharp diffraction line was observed at pressures over 6.5 GPa. The bulk modulus at 0.1 MPa and its pressure derivative of α-quartz type GeO₂ were determined to be K _T =32.8(3.3) GPa and K′ _T =6.0(2.0), respectively. In situ microscopic observations of the amorphization transformation was also performed. The large volume change due to amorphization was observed and estimated to be about 10%.     

Shock-induced formation mechanism of seifertite in shergottites

The Martian meteorites Shergotty, Zagami and Dhofar 378 have been re-investigated in order to elucidate the shock-induced formation of seifertite. The occurrence of orthorhombic seifertite (α-PbO₂ structured SiO₂) has been confirmed for the mesostasis of Shergotty and Zagami by transmission electron microscopy with lattice parameters of a = 4.05(1) Å, b = 5.05(1) Å and c = 4.45(1) Å. Seifertite crystals are interpreted as shock-induced transformation products occurring together with maskelynite of both plagioclase and alkali-feldspar composition in a largely preserved eutectic crystallisation texture. Shock-induced microstructures in accessory minerals demonstrate that these regions cannot have been completely re-molten. No further features indicating shock-pressures above ~30 GPa are detected. Hence, seifertite must have been formed below its stability field by a fast solid-state process. Significantly higher shock-pressures of Dhofar 378 indicate an inhibition of a potential seifertite crystallisation by resulting high post-shock temperatures. Crystallographic considerations reveal that a direct formation of seifertite from a high-pressure derivate of cristobalite is possible without breaking any silicon-oxygen bonds. Important implications arise from the existence of such a non-equilibrium pathway. Inferring shock-pressures from metastably formed phases appears implausible, and the transition pressure could be even below 30 GPa. Furthermore, the transformation product is determined by the precursor phase. Epitaxial intergrowth with other silica high-pressure polymorphs should be induced by certain features of the precursor, for example, planar defects, or heterogeneous strain conditions. Due to symmetrical considerations, seifertite will get amorphous during a potential back-transformation, which provides an explanation for the formation of numerous amorphous lamellae.     

Test of the vibrational modelling for the λ-type transitions: Application to the α-β quartz transition

Vibrational modelling is at the present time the only known way to predict the heat capacities of the Earth's mantle minerals at high-pressure and high-temperature. To test the validity of this method for λ-type transitions, we have applied it to the α-β quartz transition (T ₀=846±1 K). Raman spectra of quartz were recorded up to 900 K. Measured frequency shifts of the α-quartz Raman modes were then used in conjunction with available high-pressure Raman data to calculate intrinsic mode anharmonicity, through the parameter a _i=(?Lnv_i/?T)_v. Vibrational modelling of the heat capacity at constant volume, Cv, and at constant pressure, Cp, including anharmonic corrections deduced from the a _i parameters, are compared to experimental data. Taking into account the soft-mode associated to the α-β quartz transition, the model reproduces the excess of Cp related to the transition. Then, this study confirms that detecting a soft-mode from vibrational data allows one to predict λ-type transitions. However, when modelling the thermodynamic properties, the contribution of a soft-mode cannot be established from spectroscopic data. Therefore, one needs first to determine this contribution in order to predict the heat capacities of Earth's mantle minerals displaying λ-type transitions. In α-quartz, this contribution has been determined as 0.007% of the total number of the optic modes in the model of the density of states.     

The high-pressure phase diagram of synthetic epsomite (MgSO4·7H2O and MgSO4·7D2O) from ultrasonic and neutron powder diffraction measurements

We present an ultrasonic and neutron powder diffraction study of crystalline MgSO₄·7H₂O (synthetic epsomite) and MgSO₄·7D₂O under pressure up to ~3 GPa near room temperature and up to ~2 GPa at lower temperatures. Both methods provide complementary data on the phase transitions and elasticity of magnesium sulphate heptahydrate, where protonated and deuterated counterparts exhibit very similar behaviour and properties. Under compression in the declared pressure intervals, we observed three different sequences of phase transitions: between 280 and 295 K, phase transitions occurred at approximately 1.4, 1.6, and 2.5 GPa; between 240 and 280 K, only a single phase transition occurred; below 240 K, there were no phase transformations. Overall, we have identified four new phase fields at high pressure, in addition to that of the room-pressure orthorhombic structure. Of these, we present neutron powder diffraction data obtained in situ in the three phase fields observed near room temperature. We present evidence that these high-pressure phase fields correspond to regions where MgSO₄·7H₂O decomposes to a lower hydrate by exsolving water. Upon cooling to liquid nitrogen temperatures, the ratio of shear modulus G to bulk modulus B increases and we observe elastic softening of both moduli with pressure, which may be a precursor to pressure-induced amorphization. These observations may have important consequences for modelling the interiors of icy planetary bodies in which hydrated sulphates are important rock-forming minerals, such as the large icy moons of Jupiter, influencing their internal structure, dynamics, and potential for supporting life.     

两种不同燧石中α-石英 β-石英相变的研究

 利用差热分析和高温 X-射线衍射分析的实验手段,对两种不同燧石的α-石英⇔β-石英的相变反应进行了研究。结果表明,它们在相变热效应和晶格参数随温度的变化规律等方面存在着明显的差异。这种差异性的存在是由于两者内部α-石英颗粒的粒度、双晶发育程度以及燧石本身的致密度不同所造成。     

Pressure induced phase transition in Pb<Subscript>6</Subscript>Bi<Subscript>2</Subscript>S<Subscript>9</Subscript>

The crystal structure of Pb₆Bi₂S₉ is investigated at pressures between 0 and 5.6 GPa with X-ray diffraction on single-crystals. The pressure is applied using diamond anvil cells. Heyrovskyite (Bbmm, a = 13.719(4) Å, b = 31.393(9) Å, c = 4.1319(10) Å, Z = 4) is the stable phase of Pb₆Bi₂S₉ at ambient conditions and is built from distorted moduli of PbS-archetype structure with a low stereochemical activity of the Pb²⁺ and Bi³⁺ lone electron pairs. Heyrovskyite is stable until at least 3.9 GPa and a first-order phase transition occurs between 3.9 and 4.8 GPa. A single-crystal is retained after the reversible phase transition despite an anisotropic contraction of the unit cell and a volume decrease of 4.2%. The crystal structure of the high pressure phase, β-Pb₆Bi₂S₉, is solved in Pna2 ₁ (a = 25.302(7) Å, b = 30.819(9) Å, c = 4.0640(13) Å, Z = 8) from synchrotron data at 5.06 GPa. This structure consists of two types of moduli with SnS/TlI-archetype structure in which the Pb and Bi lone pairs are strongly expressed. The mechanism of the phase transition is described in detail and the results are compared to the closely related phase transition in Pb₃Bi₂S₆ (lillianite).     

Elasticity of CaSnO3 perovskite

The elastic properties of CaSnO₃ perovskite have been measured by both ultrasonic interferometry and single-crystal X-ray diffraction at high pressures. The single-crystal diffraction data collected using a diamond-anvil cell show that CaSnO₃ perovskite does not undergo any phase transitions at pressures below 8.5?GPa at room temperature. Ultrasonic measurements in the multianvil press to a maximum pressure of ～8?GPa at room temperature yielded S- and P-wave velocity data as a function of pressure. For a third-order Birch-Murnaghan EoS the adiabatic elastic moduli and their pressure derivatives determined from these velocity data are K _S0=167.2±3.1?GPa, K ^′ _S0=4.89±0.17, G ₀=89.3±1.0?GPa, G ^′ ₀=0.90±0.02. The quoted uncertainties include contributions from uncertainties in both the room pressure length and density of the specimen, as well as uncertainties in the pressure calibration of the multianvil press. Because the sample is a polycrystalline specimen, this value of K _S0 represents an upper limit to the Reuss bound (conditions of uniform stress) on the elastic modulus of CaSnO₃ perovskite. If the value of αγT is assumed to be 0.01, the value of K _S0 corresponds to K _T0=165.5±3.1?GPa. The 10 P-V data obtained by single-crystal diffraction were fit with a third-order Birch–Murnaghan equation-of-state to obtain the parameters V ₀=246.059±0.013 ų, K _T0=162.6±1.0?GPa, K ^′ _T0=5.6±0.3. Because single-crystal measurements under hydrostatic conditions are made under conditions of uniform stress, they yield bulk moduli equivalent to the Reuss bound on a polycrystalline specimen. The results from the X-ray and ultrasonic experiments are therefore consistent. The bulk modulus of CaSnO₃ perovskite lies above the linear trend of K ₀ with inverse molar volume, previously determined for Ca perovskites. This prevents an estimation of the bulk modulus of CaSiO₃ perovskite by extrapolation. However, our value of G ₀ for CaSnO₃ perovskite combined with values for CaTiO₃ and CaGeO₃ forms a linear trend of G ₀ with octahedral tilt angle. This allows a lower bound of 150?GPa to be placed on the shear modulus of CaSiO₃ by extrapolation.     

A high-pressure neutron diffraction study of the ferroelastic phase transition in RbCaF3

The fluoroperovskite phase RbCaF₃ has been investigated using high-pressure neutron powder diffraction in the pressure range ~0–7.9 GPa at room temperature. It has been found to undergo a first-order high-pressure structural phase transition at ~2.8 GPa from the cubic aristotype phase to a hettotype phase in the tetragonal space group I4/mcm. This transition, which also occurs at ~200 K at ambient pressure, is characterised by a linear phase boundary and a Clapeyron slope of 2.96 × 10^?5 GPa K^?1, which is in excellent agreement with earlier, low-pressure EPR investigations. The bulk modulus of the high-pressure phase (49.1 GPa) is very close to that determined for the low-pressure phase (50.0 GPa), and both are comparable with those determined for the aristotype phases of CsCdF₃, TlCdF₃, RbCdF₃, and KCaF₃. The evolution of the order parameter with pressure is consistent with recent modifications to Landau theory and, in conjunction with polynomial approximations to the pressure dependence of the lattice parameters, permits the pressure variation of the bond lengths and angles to be predicted. On entering the high-pressure phase, the Rb–F bond lengths decrease from their extrapolated values based on a third-order Birch–Murnaghan fit to the aristotype equation of state. By contrast, the Ca–F bond lengths behave atypically by exhibiting an increase from their extrapolated magnitudes, resulting in the volume and the effective bulk modulus of the CaF₆ octahedron being larger than the cubic phase. The bulk moduli for the two component polyhedra in the tetragonal phase are comparable with those determined for the constituent binary fluorides, RbF and CaF₂.     

High-temperature thermal expansion and elasticity of calcium-rich garnets

We present new high temperature elasticity data on two grossular garnet specimens. One specimen is single-crystal, of nearly endmember grossular, the other is polycrystalline with about 22% molar andradite. Our data extend the high temperature regime for which any garnet elasticity data are available from 1000 to 1350 K and the compositional range of temperature data to near endmember grossular. We also present new data on the thermal expansivity of calcium-rich garnet. We find virtually no discernable differences in the temperatureT derivatives at ambient conditions of the isotropic bulkK _S and shearμ moduli when comparing our results between these two specimens. These calcium-rich garnets have the lowest values of ¦(?K _S /?T) _P ¦ = (1.47,1.49) x 10^-2GPa/K, and among the highest values of ¦(?μ/?T) _P ¦ = 1.25 x 10^-2GPa/K, when compared with other garnets. Small, but measurable, nonlinear temperature dependences of most of the elastic moduli are observed. Several dimensionless parameters are computed with the new data and used to illustrate the effects of different assumptions on elastic equations of state extra-polated to high temperatures. We discuss how dimensionless parameters and other systematic considerations can be useful in estimating the temperature dependence of some properties of garnet phases for which temperature data are not yet available. While we believe it is premature to quantitatively predict the temperature variation ofK _S andμ for majorite garnets, our results have bearing on the amount of diopside required to explain the shear velocity gradients in Earth's transition zone.     

Compressibility and crystal structure of sillimanite, Al2SiO5, at high pressure

The unit-cell dimensions and crystal structure of sillimanite at various pressures up to 5.29 GPa have been refined from single-crystal X-ray diffraction data. As pressure increases, a and b decrease linearly, whereas c decreases nonlinearly with a slightly positive curvature. The axial compression ratios at room pressure are β_a:β_b:β_c=1.22:1.63:1.00. Sillimanite exhibits the least compressibility along c, but the least thermal expansivity along a (Skinner et al. 1961; Winter and Ghose 1979). The bulk modulus of sillimanite is 171(1) GPa with K′=4 (3), larger than that of andalusite (151 GPa), but smaller than that of kyanite (193 GPa). The bulk moduli of the [Al1O₆], [Al2O₄], and [SiO₄] polyhedra are 162(8), 269(33), and 367(89) GPa, respectively. Comparison of high-pressure data for Al₂SiO₅ polymorphs reveals that the [SiO₄] tetrahedra are the most rigid units in all these polymorphic structures, whereas the [AlO₆] octahedra are most compressible. Furthermore, [AlO₆] octahedral compressibilities decrease from kyanite to sillimanite, to andalusite, the same order as their bulk moduli, suggesting that [AlO₆] octahedra control the compression of the Al₂SiO₅ polymorphs. The compression of the [Al1O₆] octahedron in sillimanite is anisotropic with the longest Al1-OD bond shortening by ～1.9% between room pressure and 5.29 GPa and the shortest Al1-OB bond by only 0.3%. The compression anisotropy of sillimanite is primarily a consequence of its topological anisotropy, coupled with the compression anisotropy of the Al-O bonds within the [Al1O₆] octahedron.     

Solid carbon at high pressure: Electrical resistivity and phase transition

The electrical resistivity of polycrystalline graphite and amorphous carbon are measured at high pressures and room temperature. The results show that the resistivity of these carbon phases decreases with increasing pressure below 19 GPa. The pressure dependence of the resistivity (dln?/dP) is determined to be-0.037 GPa^?1 for the polycrystalline graphite and-0.039 GPa^?1 for the amorphous carbon. A phase transition, interpreted as the transformation to hexagonal diamond phase, is observed in the polycrystalline graphite at about 20 GPa but not in the amorphous carbon. The present experimental results support the assumption that this phase transition is martensitic in nature.     

Fe2N-type SiO2 from shocked quartz

The modified niccolite structure (Fe₂N-type) of SiO₂, previously found in diamond anvil experiments at 35 to 40 GPa, was formed in a porous mixture of crystalline α-quartz and copper powder at shock pressures estimated at 12 to 27 GPa. It is suggested that quartz can invert during shock compression not only to coesite, stishovite and an amorphous or glass phase of silica, but also to Fe₂N-type SiO₂, depending upon the shock history.     

Equation of state of adamite up to 11 GPa: a synchrotron X-ray diffraction study

The compression behavior of natural adamite [Zn₂AsO₄OH] has been investigated up to 11.07 GPa at room temperature utilizing in situ angle-dispersive X-ray diffraction and a diamond anvil cell. No phase transition has been observed within the pressure range investigated. A third-order Birch–Murnaghan equation of state fitted to all of the data points yielded V ₀ = 430.1(4) Å³, K ₀ = 80(3) GPa, K′ ₀ = 1.9(5). The K ₀ was obtained as 69(1) GPa when K′ ₀ was fixed at 4. Analysis of axial compressible moduli shows the intense compression anisotropy of adamite: K _a0 = 37(3) GPa, K _b0 = 153(6) GPa, K _c0 = 168(8) GPa; hence, a axis is the most compressible and the compressibility of b and c axis is comparable. Furthermore, the comparisons among the compressional properties of adamite, libethenite (Cu₂PO₄OH, also belongs to olivenite group), and andalusite (Al₂SiO₄O has the similar structure with adamite) at high pressure were made.     

High-pressure Raman spectroscopic study of Mn2GeO4, Ca2GeO4, Ca2SiO4, and CaMgGeO4 olivines

We present a Raman spectroscopic study of the structural modifications of several olivines at high pressures and ambient temperature. At high pressures, the following modifications in the Raman spectra are observed: 1)?in Mn₂GeO₄, between 6.7 and 8.6?GPa the appearance of weak bands at 560 and 860?cm^?1; between 10.6 and 23?GPa, the progressive replacement of the olivine spectrum by the spectrum of a crystalline high pressure phase; upon decompression, the inverse sequence of transformations is observed with some hysteresis in the transformation pressures; this sequence may be interpreted as the progressive transformation of the olivine to a spinelloid where Ge tetrahedra are polymerized, and then to a partially inverse spinel; 2)?in Ca₂SiO₄, the olivine transforms to larnite between 1.9 and 2.1?GPa; larnite is observed up to the maximum pressure of 24?GPa and it can partially back-transform to olivine during decompression; 3)?in Ca₂GeO₄, the olivine transforms to a new structure between 6.8 and 8?GPa; the vibrational frequencies of the new phase suggest that the phase transition involves an increase of the Ca coordination number and that Ge tetrahedra are isolated; this high pressure phase is observed up to the maximum pressure of 11?GPa; during decompression, it transforms to a disordered phase below 5?GPa; 4)?in CaMgGeO₄, no significant modification of the olivine spectrum is observed up to 15?GPa; between 16 and 26?GPa, broadening of some peaks and the appearance of a weak broad feature at 700–900?cm^?1 suggests a progressive amorphization of the structure; near 27?GPa, amorphization is complete and an amorphous phase is quenched down to ambient pressure; this unique behaviour is interpreted as the result of the incompatibilities in the high pressure behaviour of the Ca and Mg sublattices in the olivine structure.     

Structural and vibrational behaviour of fluorapatite with pressure. Part I: in situ single-crystal X-ray diffraction investigation

Using single-crystal X-ray diffraction from a diamond anvil cell, the compressibility of a synthetic fluorapatite was determined up to about 7?GPa. The compression pattern was anisotropic, with greater change along a than c. Unit cell parameters varied linearly with β _a =3.32(8)?10^?3 and β _c =2.40(5)?10^?3 GPa^?1, giving a ratio β _a :β _c =1.38:1. Data fitted with a third-order Birch-Murnaghan EOS yielded a bulk modulus of K ₀=93(4)?GPa with K′=5.8(1.8). The evolution of the crystal structure of fluorapatite was analysed using data collected at room pressure, at 3.04 and 4.72?GPa. The bulk modulus of phosphate tetrahedron is about three times greater than the bulk modulus of calcium polyhedra. The values were 270(10), 100(4) and 86(3) GPa for P, Ca1 (nine-coordinated) and Ca2 (seven-coordinated) respectively. While the calcium polyhedra became more regular with pressure, the distortion of the phosphate tetrahedron remained unchanged. The size of the channel extending along the [001] direction represented the most compressible direction. The Ca2–Ca2 distance decreased from 3.982 to 3.897?Å on compression from 0.0001 to 4.72?GPa. The anisotropic compressional pattern may be understood in terms of the greater compressibility of the channel size over the polyhedral units. The reduction of the channel volume was measured by the evolution of the trigonal prism, having the Ca2–Ca2–Ca2 triangle as its base and the c lattice parameter as its height. This prism volume changed from 47.3?ų at room pressure to 44.78?ų at 4.72?GPa. Its relatively high bulk moduli, 86(3) GPa, indicated that the channel did not collapse with pressure and the apatite structure could remain stable at very high pressure.     

Elasticity of CaTiO3, SrTiO3 and BaTiO3 perovskites up to 3.0 Gpa: the effect of crystallographic structure

Elasticity of CaTiO₃, SrTiO₃ and BaTiO₃ perovskites has been experimentally investigated as a function of pressure up to 3.0 GPa in a liquid-medium piston cylinder apparatus using a high precision ultrasonic interferometric technique. Specimens used are hot-pressed fine-grained (3–10 μm) polycrystalline aggregates with low porosity (<1.5%). Compressional and shear wave velocities and their pressure derivatives have been measured. The results are compared with previous studies on other perovskites and the role of structural transitions is examined. We find that the role of Ti-O₆ polyhedral tilting (such as observed in CaTiO₃) is small in the sense that a single well-defined general trend exists in perovskites with a wide range of tilting angles, although there is suggestion that cubic perovskites have slightly higher bulk modulus than orthorhombic perovskites. In contrast, cation-anion displacement that changes crystal symmetry from cubic to tetragonal in BaTiO₃ has very large effects on elasticity. This distortion significantly reduces the bulk modulus (but not much the shear modulus) and results in an unusually large pressure derivative of bulk modulus (dK/dP～10). A large change in elasticity in BaTiO₃ associated with the structural transition (without a significant volume change) is a clear example of the breakdown of the Birch's law between densities and elastic wave velocities.     

Equation of state of Fe3+-bearing phase-X

We investigated the high-pressure behaviour of Fe³⁺-bearing hydrous phase-X, (K_1.307Na_0.015)(Mg_1.504Fe _0.373 ³⁺ Al_0.053Ti _0.004 ⁴⁺ )Si₂O₇H_0.36, up to 34?GPa at room temperature by synchrotron X-ray powder diffraction. The lattice parameters behave anisotropically, with the [001] direction stiffer than [100]. In the 10^?4 to 22?GPa pressure range, the axial bulk moduli are K _0a ?=?112(3) GPa and K′?=?4, and K _0c ?=?158(2) GPa and K′?=?4, and the anisotropy of the lattice parameters is β_0c :β_0a ?=?0.71:1. The cell volumes are fitted by a second-order Birch–Murnaghan equation of state giving a bulk modulus of K ₀?=?127(1) GPa and K′?=?4 in the same pressure range. After 22?GPa, a discontinuity in volume and lattice parameters can be recognized. Sample did not become amorphous up to 34?GPa. The coupled substitution K?+?Mg?=?[]?+?Fe³⁺ has only a limited influence on the bulk modulus and structural stability of phase-X.