The elastic constants of monoclinic single-crystal chrome-diopside to 1,300 K

Values of the complete adiabatic elastic tensor for single-crystal chrome-diopside (a monoclinic pyroxene mineral) are presented from 298 to 1,300 K. The data were obtained using resonant ultrasound spectroscopy (RUS). They are the first published results for the temperature T dependences of the 13 individual elastic constants C _ij of any clinopyroxene mineral. Each C _ij is appropriately described by a linear function in T throughout the range of T. Values for each (∂C _ij /∂T) _P in GPa K⁻¹ are as follows: C ₁₁, −0.0291; C ₂₂, −0.0248; C ₃₃, −0.0179; C ₄₄, −0.0103; C ₅₅, −0.0077; C ₆₆, −0.0152; C ₁₂, −0.0119; C ₁₃, −0.0064; C ₂₃, 0.0000; C ₁₅, 0.0025; C ₂₅, 0.0022; C ₃₅, −0.0046; and C ₄₆, 0.0026. Values of (∂M/∂T) _P in GPa K⁻¹, where M represents an isotropic bulk property calculated from the C _ij data, are as follows: adiabatic bulk modulus K _S , −0.0123; isothermal bulk modulus K _T , −0.0178; and shear modulus G, −0.00998. Some diopside derivatives, notably (∂K _S /∂T) _P , (∂K _T /∂T) _P , and (∂V _P /∂T) _P , where V _P is the compressional wave velocity, have smaller magnitudes than all other minerals of importance in Earth’s mantle, thus, confirming predictions from systematics studies. We find several dimensionless quantities for this monoclinic mineral have normal values compared to other mantle minerals. Further, αK _T (α is the volume coefficient of thermal expansion) for diopside is approximately independent of both T and volume V at elevated temperature, so its equation of state is accurately expressed in simplified form.     

Triclinic elastic constants for low albite

The full set of elastic constants for plagioclase end-member phase albite (NaAlSi₃O₈) is reported for the first time. Velocities of surface acoustic waves (both Rayleigh and pseudo-surface waves) were measured using impulsively stimulated light scattering on polished surfaces having six different orientations (three normal to the Cartesian axes and three lying on diagonals). Data were inverted and results tested using several non-linear optimization techniques. Compliance moduli determined under hydrostatic compression provided additional constraints and reduced covariance in the reported constants. The Cartesian coordinate system associated with the constants (using the unit cell) has the y-axis parallel to the crystal b axis, the x-axis parallel to a* (perpendicular to b and c) and the z-axis consistent with a right-handed coordinate system. The values of the moduli C₁₁, C₁₂, C₁₃, C₁₄, C₁₅, C₁₆, C₂₂, C₂₃, C₂₄, C₂₅, C₂₆, C₃₃, C₃₄, C₃₅, C₃₆, C₄₄, C₄₅, C₄₆, C₅₅, C₅₆, C₆₆ and their 2σ uncertainties (in parentheses) are, respectively, 69.1(0.6), 34.0(0.7), 30.8(0.5), 5.1(0.1), −2.4(0.1), −0.9(0.1), 183.5(2.7), 5.5(2.2), −3.9(0.5), −7.7(0.7), −5.8(0.7), 179.5(2.3), −8.7(0.4), 7.1(0.6), −9.8(0.6), 24.9(0.1), −2.4(0.1), −7.2(0.1), 26.8 (0.2), 0.5(0.1), 33.5(0.2). These constants differ significantly from the previously reported pseudo-monoclinic constants that were based on velocity measurements on polysynthetic twinned crystal aggregates. Differences are consistent with systematic errors in the earlier study associated with sparse data and the presence of cracks and other imperfections.     

The pressure and temperature dependence of the elastic properties of single-crystal fayalite Fe2SiO4

The nine adiabatic elastic stiffness constants of synthetic single-crystal fayalite, Fe₂SiO₄, were measured as functions of pressure (range, 0 to 1.0 GPa) and temperature (range, 0 to 40° C) using the pulse superposition ultrasonic method. Summary calculated results for a dense fayalite polycrystalline aggregate, based on the HS average of our single-crystal data, are as follows: V_p = 6.67 km/s; V_s = 3.39km/s; K= 127.9 GPa; μ = 50.3 GPa; (?K/?P)_T = 5.2; (?μ/?P)_T=1.5;(?K/?T)_P= ?0.030 GPa/K;and,(?/?T)_P =-0.013 GPa/K (the pressure and temperature data are referred to 25° C and 1 atm, respectively). Accuracy of the single-crystal results was maintained by numerous cross and redundancy checks. Compared to the single-crystal elastic properties of forsterite, Mg₂SiO₄, the fayalite stiffness constants, as well as their pressure derivatives, are lower for each of the on-diagonal (C _ij for which i=j) values, and generally higher for the off-diagonal (C _ij for which i≠j) data. As a result, the bulk moduli (K) and dK/dP for forsterite and fayalite are very similar, but the rigidity modulus (μ) and dμ/dP for polycrystalline fayalite are much lower than their forsterite counterparts. The bulk compression properties derived from this study are very consistent with the static-compression x-ray results of Yagi et al. (1975). The temperature dependence of the bulk modulus of fayalite is somewhat greater (in a negative sense) than that of forsterite. The rigidity dependencies are almost equivalent. Over the temperature range relevant to this study, the elastic property results are generally consistent with the data of Sumino (1978), which were obtained using the RPR technique. However, some of the compressional modes are clearly discrepant. The elastic constants of fayalite appear to be less consistent with a theoretical HCP model (Leibfried 1955) than forsterite, reflecting the more covalent character of the Fe-O bonding in the former.     

Elasticity of single-crystal aragonite by Brillouin spectroscopy

The elastic constants of natural single-crystal aragonite (CaCO₃) have been measured by Brillouin spectroscopy at ambient conditions. The elastic constants C₁₁, C₂₂, C₃₃, C₄₄, C₅₅, C₆₆, C₁₂, C₁₃ and C₂₃ are 171.1±1.0, 110.1±0.9, 98.4±1.2, 39.3±0.6, 24.2±0.4, 40.2±0.6, 60.3±1.0, 27.8±1.6 and 41.9±2.0 GPa, respectively, for aragonite. The linear compressibilities of the a-, b- and c-axis for aragonite at ambient conditions were derived from our measured data to be 3.0±0.2, 4.2±0.2 and 7.3±0.6×10^–3 GPa^–1, respectively. The aggregate bulk and shear moduli for aragonite using the Voigt-Reuss-Hill (VRH) scheme are thus calculated to be 68.9±1.4 and 35.8±0.2 GPa, respectively. The value of bulk modulus is in remarkable contrast to the literature value of 46.9 GPa measured almost a century ago. Our new datum, however, is closer to that derived from recent atomistic simulation and static compression studies.     

Elasticity of monticellite

The adiabatic single-crystal elastic moduli of a natural sample of monticellite (CaMgSiO₄), with the olivine crystal structure, have been measured under ambient conditions using Brillouin spectroscopy. From the single-crystal moduli the aggregate bulk and shear moduli are calculated to be K _s=106±1 and = 55.2±0.4 GPa, respectively. These results are consistent with a systematic decrease in bulk modulus with increasing molar volume among the olivine-structured silicates. The longitudinal moduli decrease in the order c ₁₁>c ₃₃>c₂₂, indicating that the structure is stiffest along the a axis and most compliant along the b axis. This relationship among the longitudinal moduli holds for all silicate and germanate olivines, and is thus inferred to result from the topology of the olivine crystal structure. However, the moduli obtained in this study are at variance with previous conclusions concerning deviations from the Cauchy relations (e.g. c ₁₂=c ₄₄). For monticellite, off-diagonal shear moduli of the c ₁₂-type are uniformly greater than pure shear moduli such as c ₄₄. Similar behavior is found in pyroxenes such as diopside. The relative magnitudes of shear and off-diagonal moduli are not, therefore, a diagnostic chemical signature in minerals with complex crystal structures.Department of Applied Physics, Stanford University, Stanford CA, USA     

Elastic properties of the incommensurate phase of akermanite,Ca2MgSi2O7

The crystal structure of akermanite, Ca₂Mg-Si₂O₇, consists of mixed tetrahedral sheets formed by [MgO₄] tetrahedra and [Si₂O₇] groups interleaved along the c axis with Ca²⁺ ions in eight-fold coordination. Above 358 K, the structure is tetragonal , and below it is incommensurate with modulations parallel to [110] and . The elastic stiffness moduli, C ^ij of the incommensurate phase at room temperature were measured from wave velocities in the 20–75 MHz carrier frequency range by the ultrasonic phase comparison method using optically clear synthetic single crystal plates (3×3×2 mm) oriented parallel to (100), (001), (110) and (101) planes. The C _ij values (GPa) are: C ₁₁ 159.40, C ₃₃ 149.43, C ₄₄ 30.26, C ₆₆ 58.10, C ₁₂ 76.58 and C ₁₃ 57.80. In (010) and (001) planes, the compressional modulus, V ²(L) from the longitudinal wave, L is considerably larger than the shear moduli, V ²(T₁, T ₂) both from the in-plane and perpendicular-to-plane shear waves, T ₁ and T ₂. The relatively small values of the shear moduli indicate the ease of tetrahedral rotations in response to in-plane and perpendicular-to-plane shears and may provide preconditions for structural changes involving shear-type atomic movements.     

Elasticity of the pyrope (Mg3Al2Si3O12)-majorite (MgSiO3) garnets solid solution

Dense isotropic polycrystalline specimens of majorite-rich garnets (Py₁₀₀, Py₆₂Mj₃₈, Py₅₀Mj₅₀, Py₂₁Mj₇₉ and Mj₁₀₀) along the pyrope (Mg₃Al₂Si₃O₁₂ = Py₁₀₀)-majorite (MgSiO₃ = Mj₁₀₀) join were fabricated in a 2000-ton uniaxial split-sphere anvil apparatus (USSA-2000) at pressures from 10 to 18.5 GPa and temperatures from 1200 to 1850 °C, within their stability fields in runs of 2–4-h duration, using hot-pressing techniques developed by Gwanmesia et al. (1993). These specimens are single-phased, fine-grained (≤5 mm), free of microcracks, and have bulk densities greater than 99% of the corresponding single-crystal X-ray density. Elastic compressional (P) and shear (S) wave velocities were determined at room pressure and temperature for these polycrystalline garnet specimens by phase comparison ultrasonic interferometry. For Mj₁₀₀, the P and S wave velocities are within 1% of the Hashin-Shtrikman averages calculated from the single crystal elastic moduli measured by Brillouin spectroscopy. Both the elastic bulk modulus (K) and the shear modulus (G) decrease continuously with increasing majorite content from pyrope garnet (Py₁₀₀) to pure majorite garnet (Mj₁₀₀). The compositional dependence of K and G are given by K = 172.3 (40) − 0.085X, and G = 91.6 (10) − 0.038X, where X = mol% majorite), respectively, indicating that substitution of Si for Mg and Al decreases both K and G by about 5% along the solid solution series. Received: 25 March 1999 / Accepted: 12 July 1999     

Elasticity of superhydrous B

The adiabatic single-crystal elastic moduli of superhydrous B, Mg₁₀Si₃O₁₄(OH)₄, have been measured at ambient conditions using Brillouin spectroscopy. This material is the first hydrous phase found to be stable at the extreme conditions of 20 GPa and 1400 °C. The single-crystal moduli, in GPa, are: C ₁₁=280.0±1.5, C ₂₂=307.4±1.6, C ₃₃=293.4±1.4, C ₄₄=90.0±1.1, C ₅₅=99.2±0.8, C ₆₆=89.6±0.6, C ₁₂=66.1±2.2, C ₁₃=105.6±2.6, C ₂₃=81.8±2.6. With aggregate elastic properties of K _VRH =154.0±4.2 and μ _VRH =97.0±0.7 GPa, superhydrous B is approximately 16% suffer than forsterite and 20% softer than magnesium silicate spinel; it is also considerably more elastically isotropic than forsterite. The single-crystal moduli are compared to those of forsterite, magnesium silicate spinel and periclase, materials that are both structurally and compositionally similar to superhydrous B. The longitudinal moduli of superhydrous B and forsterite follow similar trends and appear to be dominated by the incompressibility and rotation of silicon tetrahedra. The shear and off-diagonal moduli more closely resemble those of periclase and spinel and may reflect the properties inherent to layers of magnesium octahedra.     

Elasticity of an upper mantle clinopyroxene

The ambient pressure elastic properties of a natural clinopyroxene (C2/c symmetry) from Kilbourne Hole, NM have been determined. In terms of end-members, diopside (CaMgSi₂O₆), hedenbergite (CaFeSi₂O₆), jadeite (NaAlSi₂O₆), cosmochlor (NaCrSi₂O₆), and Mg-Tschermak (MgAl(AlSi)O₆), its composition is Di₇₂He₉Jd₃Cr₃Ts₁₂. The analytic density, based on chemistry and cell parameters is 3.327 (0.003) g/cm³. The elastic constants [c₁₁, c₁₂, c₁₃, c₁₅, c₂₂, c₂₃, c₂₅, c₃₃, c₃₅, c₄₄, c₄₆, c₅₅, c₆₆] are [273.8 (0.9), 83.5 (1.3), 80.0 (1.1), 9.0 (0.6), 183.6 (0.9), 59.9 (1.6), 9.5 (1.0), 229.5 (0.9), 48.1 (0.6), 76.5 (0.9), 8.4 (0.8), 73.0 (0.4), 81.6 (1.0)] GPa where uncertainties are reported at the 95% confidence level. The aggregate (mean of Hashin-Strikman bounds) adiabatic bulk modulus is 117.2 (0.7) GPa, and the shear modulus is 72.2 (0.2) GPa. Although measured moduli are broadly consistent with trends in elasticity versus atomic volume, deviations from the systematics would produce significant (percent level) changes in calculated velocities for candidate mantle mineral assemblages. The compositional dependence of elasticity for several clinopyroxenes is explored on the basis of just the Di+He and Jd+Ts mole fractions. The bulk modulus lies within experimental uncertainties of the linear mixture of end-member properties while the shear modulus deviates by 3%. Received: 29 September 1997 / Revised, accepted: 4 March 1998     

Elastic properties of a single-crystal forsterite Mg2SiO4, up to 1,200 K

Elastic moduli of forsterite were measured between 300 and 1,200 K (? 1.6 times the Debye temperature) by the Rectangular Parallelepiped Resonance method. All the moduli decrease regularly with temperature. A summary of the results is as follows:

Pressure and temperature dependence of the elastic properties of synthetic MnO

The adiabatic elastic stiffness constants of synthetic single-crystal MnO were measured in this study using pulse superposition interferometry. Data were obtained up to 1.0 GPa in pressure and over the temperature range 273 to 473 K. As a result, we were able to determine the complete set of second-order stiffness moduli (C _ij ^s ) and their pressure and temperature derivatives, as well as higher-order properties for selected modes. Relevant results for the adiabatic bulk modulus are: K ^s=155.1±0.8 GPa; (K^s/P)_T=4.70±0.13; and, (K ^s/T)_P= -0.0203±0.0009 GPa/K. Our results for the second-order moduli are generally consistent with the data from previous studies. However, relative to the estimated uncertainties, small and systematic discrepancies appear to characterize the data set. The available evidence indicates that the differences result from microstructural variations (in particular, microcracks and Mn₃O₄ inclusions) between the synthetic MnO specimens used in different investigations. The pure shear mode C ₄₄ exhibits anomalous soft-mode behavior with both temperature (the ambient derivative is positive) and pressure (the ambient derivative is negative). In both cases the C ₄₄ data trends appear to primarily reflect the influence of Mn-Mn magnetoelastic interactions associated with the onset of a paramagnetic-antiferromagnetic (PM AFM) phase transition.     

High-pressure single-crystal X-ray diffraction and synchrotron Mössbauer study of monoclinic ferrosilite

The phase and spin transitions in single-crystal monoclinic ferrosilite, FeSiO₃, were investigated using X-ray diffraction and Mössbauer spectroscopy up to lower-mantle pressures and room temperature in a helium pressure medium. Using single-crystal X-ray diffraction, we measured the equation of state of ferrosilite up to about 43 GPa. We observed a P2₁/c-to-C2/c phase transition between 1.5 and 1.7 GPa and a phase transition from C2/c to a distinct P2₁/c structure between 30 and 34 GPa. With time-domain Mössbauer spectroscopy, we determined the hyperfine parameters of ferrous iron up to 95 GPa. The phase transitions were correlated with discontinuities in Mössbauer spectral features. We observed the onset of high-spin-to-low-spin transitions in the M1 and M2 sites at ～37 GPa and ～74 GPa, respectively. Understanding the electronic structure of iron in a well-characterized single crystal of ferrosilite may help interpret the behavior of iron in complex dense silicate phases.     

Elastic properties of pyrope

Brillouin spectroscopy was used to measure the single crystal elastic properties of a pure synthetic pyrope and a natural garnet containing 89.9 mol% of the pyrope end member (Mg₃Al₂Si₃O₁₂). The elastic moduli, c _ij , of the two samples are entirely consistent and agree with previous estimates of the elastic properties of pyrope based upon the moduli of solid solutions. Our results indicate that the elastic moduli of pyrope end-member are c ₁₁=296.2±0.5, c ₁₂=111.1±0.6, c ₄₄=91.6±0.3, K_s=172.8±0.3, =92.0±0.2, all in units of GPa. These results differ by several percent from those reported previously for synthetic pyrope, but are based upon a much larger data set. Although the hydrous components of the two samples from the present study are substantially different, representing both dry and saturated samples, we find no discernable effect of structurally bound water on the elastic properties. This is due to the small absolute solubility of water in pyrope, as compared with other garnets such as grossular.     

Elasticity of CaTiO3–CaSiO3 perovskites

Polycrystalline specimens in the CaTiO₃–CaSiO₃ perovskite system have been hot-pressed in a 2000-ton uniaxial split-sphere apparatus (USSA-2000) at pressures up to 15 GPa and temperature of 1550°C, for the compositions CaTiO₃, Ca(Ti_0.75Si_0.25)O₃, Ca(Ti_0.5Si_0.5)O₃. For the specimens with the bulk densities within 1% of the X-ray density, compressional and shear wave velocity measurements have been conducted using ultrasonic interferometry. The measured adiabatic bulk moduli (K _s ) for the CaTiO₃ and Ca(Ti_0.5Si_0.5)O₃ perovskites are 175(1) and 188(1) GPa and shear moduli (G) of 106(1) and 109(1) GPa. In situ X-ray diffraction studies at high pressure and temperature resulted in isothermal values for K ₀ of 170(5) and 185(5) GPa, respectively. For the unquenchable CaSiO₃ perovskite, elasticity theory and systematics were used to predict K ₀=212(7) GPa and G ₀=112(5) GPa; this shear modulus is 37% less than that for (Mg,Fe)SiO₃ perovskite, suggesting that CaSiO₃ perovskite cannot be ignored in modeling the composition of the Earth’s lower mantle. Received: 27 June 1997 / Revised, accepted: 25 November 1997     

Elasticity of TiO2 rutile to 1800 K

The ambient pressure elastic properties of single-crystal TiO₂ rutile are reported from room temperature (RT) to 1800 K, extending by more than 1200 oK the maximum temperature for which rutile elasticity data are available. The magnitudes of the temperature derivatives decrease with increasing temperature for five of the six adiabatic elastic moduli (C _ij ). At RT, we find (units, GPa): C ₁₁=268(1); C ₃₃=484(2); C ₄₄=123.8(2); C ₆₆=190.2(5); C ₂₃=147(1); and C ₁₂=175(1). The temperature derivatives (units, GPa K⁻¹) at RT are: (∂C ₁₁/∂T) _P =−0.042(5); (∂C ₃₃/∂T) _P =−0.087(6); (∂C ₄₄/∂T) _P =−0.0187(2); (∂C ₆₆/∂T) _P =−0.067(2); (∂C ₂₃/∂T) _P =−0.025; and (∂C ₁₂/∂T) _P −0.048(5). The values for K _S (adiabatic bulk modulus) and μ (isotropic shear modulus) and their temperature derivatives are K _S =212(1) GPa; μ=113(1) GPa; (∂K _S /∂T) _P =−0.040(4) GPa K⁻¹; and (∂μ/∂T) _P =−0.018(1) GPa K⁻¹. We calculate several dimensionless parameters over a large temperature range using our new data. The unusually high values for the Anderson-Gròneisen parameters at room temperature decrease with increasing temperature. At high T, however, these parameters are still well above those for most other oxides. We also find that for TiO₂, anharmonicity, as evidenced by a non-zero value of [∂ln (K _T )/∂lnV] _T , is insignificant at high T, implying that for the TiO₂ analogue of stishovite, thermal pressure is independent of volume (or pressure). Systematic relations indicate that ∂² K _S /∂T∂P is as high as 7×10⁻⁴ K⁻¹ for rutile, whereas ∂²μ/∂T∂P is an order of magnitude less. Received: 19 September 1997 / Revised, accepted: 27 February 1998     

The nonlinear anomalous lattice elasticity associated with the high-pressure phase transition in spodumene: a high-precision static compression study

The high-pressure behavior of the lattice elasticity of spodumene, LiAlSi₂O₆, was studied by static compression in a diamond-anvil cell up to 9.3 GPa. Investigations by means of single-crystal XRD and Raman spectroscopy within the hydrostatic limits of the pressure medium focus on the pressure ranges around ~3.2 and ~7.7 GPa, which have been reported previously to comprise two independent structural phase transitions. While our measurements confirm the well-established first-order C2/c–P2₁/c transformation at 3.19 GPa (with 1.2% volume discontinuity and a hysteresis between 0.02 and 0.06 GPa), both unit-cell dimensions and the spectral changes observed in high-pressure Raman spectra give no evidence for structural changes related to a second phase transition. Monoclinic lattice parameters and unit-cell volumes at in total 59 different pressure points have been used to re-calculate the lattice-related properties of spontaneous strain, volume strain, and the bulk moduli as a function of pressure across the transition. A modified Landau free energy expansion in terms of a one component order parameter has been developed and tested against these experimentally determined data. The Landau solution provides a much better reproduction of the observed anomalies than any equation-of-state fit to data sets truncated below and above P _tr, thus giving Landau parameters of K ₀ = 138.3(2) GPa, K′ = 7.46(5), λ _V  = 33.6(2) GPa, a = 0.486(3), b = −29.4(6) GPa and c = 551(11) GPa.     

Equation of state, elasticity, and shear strength of pyrite under high pressure

 Physical properties including the equation of state, elasticity, and shear strength of pyrite have been measured by a series of X-ray diffraction in diamond-anvil cells at pressures up to 50 GPa. A Birch–Murnaghan equation of state fit to the quasihydrostatic pressure–volume data obtained from laboratory X-ray source/film techniques yields a quasihydrostatic bulk modulus K _0T =133.5 (±5.2) GPa and bulk modulus first pressure derivative K ^′ _0T =5.73 (±0.58). The apparent equation of state is found to be strongly dependent on the stress conditions in the sample. The stress dependency of the high-pressure properties is examined with anisotropic elasticity theory from subsequent measurements of energy-dispersive radial diffraction experiments in the diamond-anvil cell. The calculated values of K _0T depend largely upon the angle ψ between the diffracting plane normal and the maximum stress axis. The uniaxial stress component in the sample, t=σ₃−σ₁, varies with pressure as t=−3.11+0.43P between 10 and 30 GPa. The pressure derivatives of the elastic moduli dC ₁₁/dP=5.76 (±0.15), dC ₁₂/dP=1.41 (±0.11) and dC ₄₄/dP=1.92 (±0.06) are obtained from the diffraction data assuming previously reported zero-pressure ultrasonic data (C ₁₁=382 GPa, C ₁₂=31 GPa, and C ₄₄=109 GPa). Received: 21 December 2000 / Accepted: 11 July 2001     

A theoretical investigation of the relative stabilities of Fe-free clinozoisite and orthozoisite

 The relative stabilities of orthozoisite, Ca₂Al₃[O|OH|Si₂O₇|SiO₄], space group Pnma, and the monoclinic polymorph, clinozoisite, space group P2₁/m, have been investigated using calculations based on density functional theory. It is found that orthozoisite is more stable than clinozoisite by about 1 kJ mol⁻¹ at zero pressure in the athermal limit. The bulk moduli of the two polymorphs have been calculated to be B_ortho=117.5(1.7) GPa and B_clino=136(4) GPa. Received: 20 March 2000 / Accepted: 26 February 2001     

Elasticity of a beryllium silicate (phenacite: Be2SiO4)

The adiabatic single-crystal elastic moduli of a beryllium silicate (phenacite: Be₂SiO₄, trigonal, have been determined at atmospheric pressure and 22° C by Brillouin spectroscopy. The elastic stiffness moduli in gigapascals are: C ₁₁=341.9 C ₃₃=391.0 C ₄₄=91.4 C ₆₆= 96.9 C ₁₂=148.0 C ₁₃=136.0 C ₁₄= 0.1 C ₁₅= 3.5Overall, the elastic stiffness moduli for phenacite parallel and perpendicular to the c axis are comparable (i.e., it is almost cubic in its elastic signature). The elastic moduli can be rationalized in terms of division of the structure into two types of coordination polyhedra (1Si+2Be) with slightly different stiffnesses, which are linked to form a three dimensional framework. Values of the isothermal bulk modulus and the linear compressibilities, as determined from hydrostatic compression experiments of Hazen and Au (1986), are in good agreement with those obtained here. Combining the two studies indicates a low pressure derivative of the bulk modulus for phenacite.     

Single-crystal elastic properties of chondrodite: implications for water in the upper mantle

Chondrodite, a member of the humite group of minerals, forms by hydration of olivine and is stable over a range of temperatures and pressures that includes a portion of the uppermost mantle. We have measured the single crystal elastic properties of a natural chondrodite specimen at ambient conditions using Brillouin spectroscopy. The isotropic aggregate bulk (K) and shear (μ) moduli calculated from the single-crystal elastic moduli, C_ij, are: K_S=118.4(16) GPa and μ=75.6(7) GPa. A comparison of the structures and elasticity of olivine and chondrodite indicate that the replacement of O with (OH,F) in M²⁺O₆ octahedra has a small effect on the elasticity of humite-group minerals. The slightly diminished elastic moduli of humite-group minerals (as compared to olivine) are likely caused by a smaller ratio of strong structural elements (SiO₄ tetrahedra) to weaker octahedra, and perhaps a more flexible geometry of edge-sharing MO₄(O,OH,F)₂ octahedra. In contrast to the humite-olivine group minerals, the incorporation of water into garnets and spineloids leads to a more substantial decrease in the elastic properties of these minerals. This contrasting behavior is due to formation of O₄H₄ tetrahedra and vacant hydroxyl-bearing octahedra in the garnets and spineloids, respectively. Therefore, the mechanism of incorporation of H/OH into mineral phases, not only degree of hydration, should be taken into account when estimating the effect of water on the elastic properties of minerals. The bulk elastic wave velocities of chondrodite and olivine are very similar. If humite-like incorporation of OH is predominant in the upper mantle, then the reaction of OH with olivine will have a minor or possibly no detectable effect on seismic velocities. Thus, it may be difficult to distinguish chondrodite-bearing rocks from “anhydrous” mantle on the basis of seismically determined velocities for the Earth. Received: 25 February 1998 / Revised, accepted: 18 August 1998