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.     

High-pressure neutron diffraction study on H–D isotope effects in brucite

A neutron powder diffraction study of hydrogenated and deuterated brucite was conducted at ambient temperature and at pressures up to 9 GPa, using a Paris–Edinburgh high-pressure cell at the WAND instrument of the ORNL High Flux Isotope Reactor. The two materials were synthesized by the same method and companion measurements of neutron diffraction were conducted under the same conditions. Our refinement results show that the lattice-parameters of the a axis, parallel to the sheets of Mg–O octahedra, decrease only slightly with pressure with no effect of H–D substitution. However, the c axis of Mg(OD)₂ is shorter and may exhibit greater compressibility with pressure than that of Mg(OH)₂. Consequently, the unit-cell volume of deuterated brucite is slightly, but systematically smaller than that of hydrogenated brucite. When fitted to a third-order Birch–Murnaghan equation in terms of the normalized unit-cell volume, values of the bulk modulus for hydrogenated and deuterated brucite (K ₀ = 39.0 ± 2.8 and 40.4 ± 1.3 GPa, respectively) are, however, indistinguishable from each other within the experimental errors. The measured effect of H–D substitution on the unit-cell volume also demonstrates that brucite (and other hydrous minerals) preferentially incorporate deuterium over hydrogen under pressure, suggesting that the distribution of hydrogen isotopes in deep-earth conditions may differ significantly from that in near-surface environments.     

High-pressure properties of diaspore,AlO(OH)

The structural compression mechanism and compressibility of diaspore, AlO(OH), were investigated by in situ single-crystal synchrotron X-ray diffraction at pressures up to 7 GPa using the diamond-anvil cell technique. Complementary density functional theory based model calculations at pressures up to 40 GPa revealed additional information on the pressure-dependence of the hydrogen-bond geometry and the vibrational properties of diaspore. A fit of a second-order Birch–Murnaghan equation of state to the p–V data resulted in the bulk modulus B ₀ = 150(3) GPa and B ₀ = 150.9(4) GPa for the experimental and theoretical data, respectively, while a fit of a third-order Birch–Murnaghan equation of state resulted in B ₀ = 143.7(9) GPa with its pressure derivative B′ = 4.4(6) for the theoretical data. The compression is anisotropic, with the a-axis being most compressible. The compression of the crystal structure proceeds mainly by bond shortening, and particularly by compression of the hydrogen bond, which crosses the channels of the crystal structure in the (001) plane, in a direction nearly parallel to the a-axis, and hence is responsible for the pronounced compression of this axis. While the hydrogen bond strength increases with pressure, a symmetrisation is not reached in the investigated pressure range up to 40 GPa and does not seem likely to occur in diaspore even at higher pressures. The stretching frequencies of the O–H bond decrease approximately linearly with increasing pressure, and therefore also with increasing O–H bond length and decreasing hydrogen bond length. Electronic Supplementary Material The online version of this article () contains supplementary material, which is available to authorized users.     

Elasticity of CaIrO<Subscript>3</Subscript> with perovskite and post-perovskite structure

Compression behaviors of CaIrO₃ with perovskite (Pv) and post-perovskite (pPv) structures have been investigated up to 31.0(1.0) and 35.3(1) GPa at room temperature, respectively, in a diamond-anvil cell with hydrostatic pressure media. CaIrO₃ Pv and pPv phases were compressed with the axial compressibility of β _a > β _c > β _b and β _b > β _a > β _c, respectively and no phase transition was observed in both phases up to the highest pressure in the present study. The order of axial compressibility for pPv phase is consistent with the crystallographic consideration for layer structured materials and previous experimental results. On the other hand, Pv phase shows anomalous compression behavior in b axis, which exhibit constant or slightly expanded above 13 GPa, although the applied pressure remained hydrostatic. Volume difference between Pv and pPv phases was gradually decreased with increasing pressure and this is consistent with the results of theoretical study based on the ab initio calculation. Present results, combined with theoretical study, suggest that these complicate compression behaviors in CaIrO₃ under high pressure might be caused by the partially filled electron of Ir⁴⁺. Special attention must be paid in case of using CaIrO₃ as analog materials to MgSiO₃, although CaIrO₃ exhibits interesting physical properties under high pressure.     

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     

Effects of Fe/Mg on the compressibility of synthetic wadsleyite: β-(Mg1-xFex)2SiO4 (x≤0.25)

Four crystals of synthetic wadsleyite, -(Mg,Fe)₂SiO₄, were mounted together in one diamond-anvil cell for the determination of unit-cell parameters as a function of pressure. The Fe/(Fe+Mg) are 0.00, 0.08, 0.16, and 0.25 (the most iron-rich stable composition). Unit-cell refinements were made at 12 pressures up to 4.5 GPa. No phase transitions were observed and all crystals remained dimensionally orthorhombic. Of the three axes, c is the most compressible (0.000239(3) GPa^-1), whereas compressibilities of a and b are both about 30% less. The Fe content has no systematic effect on volume or linear compressibilities. Bulk moduli, based on a Birch-Murnaghan equation of state (K assumed to be 4.00) are 160(3), 169(3), 164(2), and 165(3) GPa for the four crystals in order of increasing Fe. Substitution of Fe for Mg, therefore, does not appear to have a systematic effect on bulk modulus. Other factors, especially Fe³⁺/Fe²⁺ and other deviations from the strict Mg₂SiO₄-Fe₂SiO₄ binary, may have a greater influence on compressibility.     

Equation of state of hexagonal aluminous phase in basaltic composition to 63 GPa at 300 K

A laser-heated diamond-anvil cell that is capable of operating up to a pressure of 63 GPa, with X-ray diffraction facilities using a synchrotron radiation source at the SPring-8, has been developed to observe the compressibility of a hexagonal aluminous phase, [K_0.15Na_1.66Ca_0.11Mg_1.29Fe²⁺ _0.86Al_3.13Ti_0.09Si_1.98] _Σ9.27O₁₂. The hexagonal aluminous phase is a potassium host mineral from the subducted oceanic crust in the Earth's lower mantle. A sample was heated using a YAG laser at each pressure increment to relax the deviatoric stress in the sample. X-ray diffraction measurements were carried out at 300 K using an angle-dispersive technique. Pressure was measured using an internal platinum pressure calibrant. The observed unit-cell volumes were used to obtain a third-order Birch–Murnaghan equation of state: unit-cell volume V _o=185.94(±16) ų, density ρ _o=4.145 g/cm³, and bulk modulus K _o=198(±3) GPa when the first pressure is derivative of the bulk modulus K ^′ _o is fixed to 4. The density of hexagonal aluminous phase is lower than that of coexisting Mg-perovskite in the subducted oceanic crust.     

Unquenchable high-pressure perovskite polymorphs of MnSnO3 and FeTiO3

New high-pressure orthorhombic (GdFeO₃-type) perovskite polymorphs of MnSnO₃ and FeTiO₃ have been observed using in situ powder X-ray diffraction in a diamond-anvil cell with synchrotron radiation. The materials are produced by the compression of the lithium niobate polymorphs of MnSnO₃ and FeTiO₃ at room temperature. The lithium niobate to perovskite transition occurs reversibly at 7 GPa in MnSnO₃, with a volume change of -1.5%, and at 16 GPa in FeTiO₃, with a volume change of -2.8%. Both transitions show hysteresis at room temperature. For MnSnO₃ perovskite at 7.35 (8) GPa, the orthorhombic cell parameters are a=5.301 (2) A, b=5.445 (2) Å, c=7.690 (8) Å and V= 221.99 (15) ų. Volume compression data were collected between 7 and 20 GPa. The bulk modulus calculated from the compression data is 257 (18) GPa in this pressure region. For FeTiO₃ perovskite at 18.0 (5) GPa, cell parameters are a=5.022 (6) Å, b=5.169 (5) Å, c=7.239 (9) Å and V= 187.94 (36) ų. Based on published data on the quench phases, the FeTiO₃ perovskite breaks down to a rocksalt + baddelyite mixture of FeO and TiO₂ at 23 GPa. This is the first experimental verification of the pressure-induced breakdown of a perovskite to simple oxides.     

In situ determination of crystal structure for high pressure phase of ZrO2 using a diamond anvil and single crystal X-ray diffraction method

Combining a miniature diamond-anvil pressure cell with a single crystal four-circle diffractometer, the crystal structure of a synthetic ZrO₂ has been studied in situ up to 51 kbar at room temperature. The space group of the unquenchable orthorhombic high pressure phase is Pbcm. The directions of the b and c axes are preserved through the transition and the transformation is displacive. The coordination configurations of the Zr atoms and oxygen atoms are the same in the high pressure and low pressure phases. The orthorhombic high pressure phase has a higher entropy than that of low pressure monoclinic phase.     

Single-crystal hydrostatic compression of synthetic pyrope, almandine, spessartine, grossular and andradite garnets at high pressures

The compression of synthetic pyrope Mg₃Al₂ (SiO₄)₃, almandine Fe₃Al₂(SiO₄)₃, spessartine Mn₃Al₂ (SiO₄)₃ grossular Ca₃Al₂(SiO₄)₃ and andradite Ca₃Fe₂ (SiO₄)₃ was studied by loading the crystals together in a diamond anvil cell. The unit-cell parameters were determined as a function of pressure by X-ray diffraction up to 15 GPa using neon as a pressure transmitting medium. The unit-cell parameters of pyrope and almandine were measured up to 33 and 21 GPa, respectively, using helium as a pressure medium. The bulk moduli, K _{T 0}, and their first pressure derivatives, K _{T 0} ^′, were simultaneously determined for all five garnets by fitting the volume data to a third order Birch-Murnaghan equation of state. Both parameters can be further constrained through a comparison of volume compressions between pairs of garnets, giving for K _{T 0} and K _{T 0} ^′ 171(2) GPa and 4.4(2) for pyrope, 185(3) GPa and 4.2(3) for almandine, 189(1) GPa and 4.2 for spessartine, 175(1) GPa and 4.4 for grossular and 157(1) GPa and 5.1 for andradite, where the K _{T 0} ^′ are fixed in the case of spessartine, grossular and andradite. Direct comparisons of the unit-cell volumes determined at high pressures between pairs of garnets reveal anomalous compression behavior for Mg²⁺ in the 8-fold coordinated triangular dodecahedron in pyrope. This agrees with previous studies concerning the compression behaviors of Mg²⁺ in 6-fold coordinated polyhedra at high pressures. The results show that simple bulk modulus–volume systematics are not obeyed by garnets. Received: 29 July 1998 / Revised, accepted: 7 April 1999     

The high-pressure C2/c–P21/c phase transition along the LiAlSi2O6–LiGaSi2O6 solid solution

Two synthetic single-crystals with composition Li(Al_0.53Ga_0.47)Si₂O₆ and LiGaSi₂O₆ and space group C2/c at room conditions have been studied under pressure by means of X-ray diffraction using a diamond anvil cell. The unit-cell parameters were determined at 12 and 10 different pressures up to P = 8.849 and P = 7.320 GPa for Li(Al_0.53Ga_0.47)Si₂O₆ and LiGaSi₂O₆, respectively. The sample with mixed composition shows a C2/c to P2₁/c phase transformation between 1.814 and 2.156 GPa, first-order in character. The transition is characterised by a large and discontinuous decrease in the unit-cell volume and by the appearance of the b-type reflections (h + k = odd) typical of the primitive symmetry. The Ga end-member shows the same C2/c to P2₁/c transformation at a pressure between 0.0001 and 0.39 GPa. The low-pressure value at which the transition occurred did not allow collecting any data in the C2/c pressure stability field except that on room pressure. Our results compared with those relative to spodumene (LiAlSi₂O₆, Arlt and Angel 2000a) indicate that the substitution of Al for Ga at the M1 site of Li-clinopyroxenes strongly affects the transition pressure causing a decrease from 3.17 GPa (spodumene) to less than 0.39 GPa (LiGaSi₂O₆) and decreases the volume discontinuity at the transition. As already found for other compounds, the C2/c low-pressure phases are more rigid than the P2₁ /c high-pressure ones. Moreover, the increase of the M1 cation radius causes a decrease in the bulk modulus K _T0. The axial compressibility among the Li-bearing clinopyroxenes indicates that the c axis is the most rigid for the C2/c phases while it becomes the most compressible for the P2₁ /c phases.     

Room-temperature equation of state of Li<Subscript>2</Subscript>VOSiO<Subscript>4</Subscript> up to 8.5 GPa

High-pressure single-crystal X-ray diffraction measurements of lattice parameters of the compound Li₂VOSiO₄, which crystallises with a natisite-type structure, has been carried out to a pressure of 8.54(5) GPa at room temperature. Unit-cell volume data were fitted with a second-order Birch-Murnaghan EoS (BM-EoS), simultaneously refining V ₀ and K ₀ using the data weighted by the uncertainties in V. The bulk modulus is K ₀ = 99(1) GPa, with K′ fixed to 4. Refinements of third order equations-of-state yielded values of K′ that did not differ significantly from 4. The compressibility of the unit-cell is strongly anisotropic with the c axis (K ₀(c) = 49.7 ± 0.5 GPa) approximately four times more compressible than the a axis (K ₀(a) = 195 ± 3 GPa).     

High-pressure behaviour of Ca-rich C 2 /c clinopyroxenes along the join diopside-enstatite (CaMgSi2O6-Mg2Si2O6)

 An in situ high-pressure (HP) X-ray diffraction investigation of synthetic diopside and of the Ca_0.8Mg_1.2Si₂O₆ clinopyroxene (Di₈₀En₂₀) was performed up to respectively P=40.8 and 15.1 GPa, using high brilliance synchrotron radiation. The compression of the cell parameters is markedly anisotropic, with β_b ⋙ β_c > β_a > β_asinβ for any pressure range and for both diopside and Di₈₀En₂₀. The compressibility along the crystallographic axes decreases significantly with pressure and is higher in Di₈₀En₂₀ than in diopside. The β cell parameter decreases as well with pressure, at a higher rate in Di₈₀En₂₀. The cell volume decreases at almost the same rate for the two compositions, since in diopside a higher compression along a* occurs. A change in the mechanism of deformation at P higher than about 5–10 GPa is suggested for both compositions from the analysis of the strain induced by compression. In diopside at lower pressures, the deformation mainly occurs, at a similar rate, along the b axis and at a direction 145° from the c axis on the (0 1 0) plane. At higher pressures, instead, the deformation occurs mostly along the b axis. In Di₈₀En₂₀ the orientation of the strain axes is the same as in diopside. The substitution of Ca with Mg in the M2 site induces at a given pressure a higher deformation on (0 1 0) with respect to diopside, but a similar change in the compressional behaviour is found. Changes in the M2 polyhedron with pressure can explain the above compressional behaviour. A third-order Birch-Murnaghan equation of state was fit to the retrieved volumes, with K=105.1(9) GPa, K′=6.8(1) for diopside and K=107.3(1.4) GPa, K′=5.7(3) for Di₈₀En₂₀; the same equation can be applied for any pressure range. The elasticity of diopside is therefore not significantly affected by Mg substitution into the M2 site, in contrast to the significant stiffening occurring for Ca substitution into Mg-rich orthopyroxenes. Received: 3 January 2000 / Accepted: 21 May 2000     

The compressibility of a natural apatite

In-situ synchrotron X-ray diffraction (XRD) experiments of a natural apatite with the formula of Ca₅(PO₄)₃F_0.94Cl_0.06 were carried out using a diamond anvil cell and angle-dispersive technique at Photon Factory (PF), Japan. Pressure–volume data were collected up to 7.12 GPa at 300 K. The pressures were determined from the ruby fluorescence spectra shift. The unit-cell parameters and volume decreased systematically with increasing pressure, and a reliable isothermal bulk modulus and its pressure derivative were obtained in this study. The third-order Birch–Murnaghan equation of state yielded the isothermal bulk modulus of K_T=91.5(38) GPa, its pressure derivative K_T= 4.0(11), and the zero-pressure volume V₀=524.2(3) ų.     

Effect of incorporation of iron and aluminum on the thermoelastic properties of magnesium silicate perovskite

In situ X-ray diffraction measurements of Fe- and Al-bearing MgSiO₃-rich perovskite (FeAl-Pv), which was synthesized from a natural orthopyroxene, were performed at pressures of 19–32 GPa and temperatures of 300–1,500 K using a combination of a Kawai-type apparatus with eight sintered-diamond anvils and synchrotron radiation. Two runs were performed using a high-pressure cell with two sample chambers, and both MgSiO₃ perovskite (Mg-Pv) and FeAl-Pv were synthesized simultaneously in the same cell. Thus we were able to measure specific volumes (V/V ₀) of Mg-Pv and FeAl-Pv at the same P−T conditions. At all the measurement conditions, values of the specific volume of FeAl-Pv are consistent with those of Mg-Pv within 2 Standard Deviation, strongly suggesting that effect of incorporation of iron and aluminum on the thermoelastic properties of magnesium silicate perovskite is undetectable in this composition, pressure, and temperature range. Two additional runs were performed using a high-pressure cell that has one sample chamber and unit-cell volumes of FeAl-Pv were measured at pressures and temperatures up to 32 GPa and 1,500 K, respectively. All the unit-cell volume data of FeAl-Pv perovskite were fitted to the high temperature Birch–Murnaghan equation of state and a complete set of thermoelastic parameters of this perovskite was determined with an assumption of K′ _300,0 = 4. The determined parameters are K _300,0 = 243(3) GPa, (∂K _T,0/∂T) _P = −0.030(8) GPa/K, a ₀ = 2.78(18) × 10⁻⁵ K⁻¹, and b ₀ = 0.88(28) × 10⁻⁸ K⁻², where a ₀ and b ₀ are the coefficients of the following expression describing the zero-pressure thermal expansion: α _T,0 = a ₀ + b ₀ T. The equation-of-state parameters of FeAl-Pv are in good agreement with those of MgSiO₃ perovskite at the conditions corresponding to the uppermost part of the lower mantle.     

High-pressure behavior of MnGeO3 and CdGeO3 perovskites and the post-perovskite phase transition

The stability and high-pressure behavior of perovskite structure in MnGeO₃ and CdGeO₃ were examined on the basis of in situ synchrotron X-ray diffraction measurements at high pressure and temperature in a laser-heated diamond-anvil cell. Results demonstrate that the structural distortion of orthorhombic MnGeO₃ perovskite is enhanced with increasing pressure and it undergoes phase transition to a CaIrO₃-type post-perovskite structure above 60 GPa at 1,800 K. A molar volume of the post-perovskite phase is smaller by 1.6% than that of perovskite at equivalent pressure. In contrast, the structure of CdGeO₃ perovskite becomes less distorted from the ideal cubic perovskite structure with increasing pressure, and it is stable even at 110 GPa and 2,000 K. These results suggest that the phase transition to post-perovskite is induced by a large distortion of perovskite structure with increasing pressure.     

High-pressure metallization and electronic-magnetic propertiesof hexagonal cubanite (CuFe2S3)

^?57Fe Mössbauer studies at room temperature and temperature-dependent resistance studies have been performed on a natural specimen of cubanite (CuFe₂S₃) in a diamond-anvil cell at pressures up to ～10 GPa. An insulator-metal phase transition occurs in the range 3.4–5.8 GPa coinciding with a previously observed structural transition from an orthorhombic to a hexagonal NiAs (B8) structure. The room temperature data shows that the metallization process concurs with a gradual transition from a magnetically ordered phase at low pressure to a nonmagnetic or paramagnetic phase at high-pressure. The change in magnetic behaviour at the structural transition may be attributed to a reduction of the Fe-S-Fe superexchange angle formed by edge-sharing octahedra occurring in the high-pressure phase. The non-magnetic or paramagnetic metallic phase at high pressure is retained upon decompression to ambient pressure-temperature conditions, indicative of substantial hysteresis associated with the pressure driven orthorhombic→hexagonal structural transition. The pressure evolution of both the ⁵⁷Fe Mössbauer hyperfine interaction parameters and resistance behaviour is consistent with the transition from mixed-valence character in the low pressure orthorhombic structure to that of extended-electron delocalization in the hexagonal phase at high-pressure.     

On the elastic behaviour of zeolite mordenite: a synchrotron powder diffraction study

The high-pressure elastic behaviour of a synthetic zeolite mordenite, Na₆Al_6.02Si_42.02O₉₆·19H₂O [a=18.131(2), b=20.507(2), c=7.5221(5) Å, space group Cmc2₁], has been investigated by means of in situ synchrotron X-ray powder diffraction up to 5.68 GPa. No phase transition has been observed within the pressure range investigated. Axial and volume bulk moduli have been calculated using a truncated second-order Birch–Murnaghan equation-of-state (II-BM-EoS). The refined elastic parameters are: V ₀=2801(11) ų, K _T0= 41(2) GPa for the unit-cell volume; a ₀=18.138(32) Å, K _T0(a)=70(8) GPa for the a-axis; b ₀=20.517(35) Å, K _T0(b)=29(2) GPa for the b-axis and c ₀=7.531(5) Å, K _T0(c)=38(1) GPa for the c-axis [K _T0(a): K _T0(b): K _T0(c)=2.41:1.00:1.31]. Axial and volume Eulerian finite strain versus “normalized stress” plots (fe–Fe plot) show an almost linear trend and the weighted linear regression through the data points yields the following intercept values: Fe(0)=39(4) GPa for V; Fe _a (0)=65(18) GPa for a; Fe _b (0)=28(3) GPa for b; Fe _c (0)=38(2) GPa for c. The magnitudes of the principal Lagrangian unit-strain coefficients, between 0.47 GPa (the lowest HP-data point) and each measured P>0.47 GPa, were calculated. The unit-strain ellipsoid is oriented with ε₁ || b, ε₂ || c, ε₃ || a and |ε₁|> |ε₂|> |ε₃|. Between 0.47 and 5.68 GPa the relationship between the unit-strain coefficient is ε₁: ε₂: ε₃=2.16:1.81:1.00. The reasons of the elastic anisotropy are discussed.An erratum to this article can be found at     

High-pressure behaviour along the jadeite NaAlSi2O6–aegirine NaFeSi2O6 solid solution up to 10 GPa

The volume variation as a function of pressure along the jadeite–aegirine solid solution was determined at room temperature up to pressures between 6.5 and 9.7 GPa by single-crystal X-ray diffraction. The unit-cell volumes collected at room pressure for the different compositions indicate a slight deviation from linearity along the join. The pressure–volume data have been fitted using a third-order Birch-Murnaghan equation of state (BM3-EoS). The bulk modulus, K _T0, varies from 134.0(7) GPa for pure jadeite to 116.1(5) GPa for pure aegirine. Its evolution with composition along the join is not linear and can be described by the following second order polynomial:

Equation of state and pressure-induced structural changes in mirabilite (Na2SO4·10H2O) determined from ab initio density functional theory calculations

We have carried out ab initio calculations using density functional theory to determine the bulk elastic properties of mirabilite, Na₂SO₄·10H₂O, and to obtain information on structural trends caused by the application of high pressure up to ~60 GPa. We have found that there are substantial isosymmetric discontinuous structural re-organisations at ~7.7 and ~20 GPa caused by changes in the manner in which the sodium cations are coordinated by water molecules. The low-pressure and intermediate-pressure phases both have sodium in sixfold coordination but in the high-pressure phase the coordination changes from sixfold to sevenfold. These coordination changes force a re-arrangement of the hydrogen-bond network in the crystal. The trend is towards a reduction in the number of hydrogen bonds donated to the sulphate group (from twelve down to six over the range 0–60 GPa) and an increase in hydrogen bonding amongst the Na-coordinated water molecules and the two interstitial water molecules. Ultimately, we observe proton transfers from the interstitial waters (forming OH^? ions) to two of the Na-coordinated waters (forming a pair of H₃O⁺ ions). The equation of state in the athermal limit of the low-pressure phase of mirabilite, parameterised by fitting an integrated form of the third-order Birch-Murnaghan expression to the calculated energy as a function of unit-cell volume, yields the zero-pressure unit-cell volume, V ₀ = 1468.6(9) ų, the incompressibility, K ₀ = 22.21(9) GPa, and the first pressure derivative K ₀′ = (?K/?P)₀ = 5.6(1).