High-pressure crystal chemistry of chrysoberyl,Al2BeO4: Insights on the origin of olivine elastic anisotropy

High-pressure crystal structure refinements and axial compressibilities have been determined by x-ray methods for the olivine isomorph chrysoberyl, Al₂BeO₄. Unlike silicate olivines, which are more than twice as compressible along b than along a, chrysoberyl (space group Pbnm) has nearly isotropic compressibility with β _a =1.12±0.04, β _b =1.46±0.05, and β _c =1.31±0.03 (all×10^?4 kbar^?1). The resultant bulk modulus is 2.42±0.05 Mbar, with K′ assumed to be 4. The axial compression ratios of chrysoberyl are 1.00:1.30:1.17, compared to axial compression ratios 1.00:2.02:1.60 for forsterite. These differences in compression anisotropy arise from differences in relative bond compressibilities. In chrysoberyl the average aluminum-oxygen and beryllium-oxygen bond compressibilities are similar, yielding nearly isotropic compression, but in silicate olivines octahedral cation-oxygen bonds are significantly more compressible than Si-O bonds, so that compression parallel to a is much more restricted than that parallel to b. The inherent anisotropy of the olivine structure is not, by itself, sufficient to cause anisotropic compression. It appears that in the case of olivine the distribution of cations of different valences, in conjunction with the structure type, leads to anisotropies in physical properties.     

Hydrostatic compression of ilmenite phase of ZnSiO3 and MgGeO3

Accurate measurements of cell parameters were performed on the ilmenite phases of ZnSiO₃ and MgGeO₃ using an X-ray diffraction method under hydrostatic conditions. The linear changes in cell parameter are represented by 1?a/a ₀=(1.06±0.04)×10^?4 P(kbar) and 1?c/c ₀=(2.11±0.04)×10^?4 P for ZnSiO₃, and 1?a/a ₀=(1.37±0.03)×10^?4 P and 1?c/c ₀=(2.05±0.04)×10^?4 P for MgGeO₃. A least-squares calculation using the first-order Birch-Murnaghan equation gives K _T =2.16±0.02 Mbar and K _T =1.87±0.02 Mbar for ZnSiO₃ and MgGeO₃, respectively. Elastic systematics assuming K _T V _m =constant give a predicted value K _T =2.14 Mbar for the ilmenite phase of MgSiO₃.     

High-Pressure crystal chemistry of spinel (MgAl2O4) and magnetite (Fe3O4): Comparisons with silicate spinels

High-pressure crystal structures and compressibilities have been determined by x-ray methods for MgAl₂O₄ spinel and its isomorph magnetite, Fe₃O₄. The measured bulk moduli, K, of spinel and magnetite (assuming K′=4) are 1.94±0.06 and 1.86±0.05 Mbar, respectively, in accord with previous ultrasonic determinations. The oxygen u parameter, the only variable atomic position coordinate in the spinel structure (Fd3m, Z=8), decreases with pressure in MgAl₂O₄, thus indicating that the magnesium tetrahedron is more compressible than the aluminum octahedron. In magnetite the u parameter is unchanged, and both tetrahedron and octahedron display the 1.9 Mbar bulk modulus characteristic of the entire crystal. This behavior contrasts with that of nickel silicate spinel (γ-Ni₂SiO₄), in which the u parameter increases with pressure because the silicon tetrahedron is relatively incompressible compared to the nickel octahedron.     

High-pressure polymorphism and structural transitions of norsethite,BaMg(CO3)2

In situ high-pressure investigations on norsethite, BaMg(CO₃)₂, have been performed in sequence of diamond-anvil cell experiments by means of single-crystal X-ray and synchrotron diffraction and Raman spectroscopy. Isothermal hydrostatic compression at room temperature yields a high-pressure phase transition at P _c ≈ 2.32 ± 0.04 GPa, which is weakly first order in character and reveals significant elastic softening of the high-pressure form of norsethite. X-ray structure determination reveals C2/c symmetry (Z = 4; a = 8.6522(14) Å, b = 4.9774(13) Å, c = 11.1542(9) Å, β = 104.928(8)°, V = 464.20(12) Å³ at 3.00 GPa), and the structure refinement (R ₁ = 0.0763) confirms a distorted, but topologically similar crystal structure of the so-called γ-norsethite, with Ba in 12-fold and Mg in octahedral coordination. The CO₃ groups were found to get tilted off the ab-plane direction by ~16.5°. Positional shifts, in particular of the Ba atoms and the three crystallographically independent oxygen sites, give a higher flexibility for atomic displacements, from which both the relatively higher compressibility and the remarkable softening originate. The corresponding bulk moduli are K ₀ = 66.2 ± 2.3 GPa and dK/dP = 2.0 ± 1.8 for α-norsethite and K ₀ = 41.9 ± 0.4 GPa and dK/dP = 6.1 ± 0.3 for γ-norsethite, displaying a pronounced directional anisotropy (α: β _a ^?1  = 444(53) GPa, β _c ^?1  = 76(2) GPa; γ: β _a ^?1  = 5.1(1.3) × 10³ GPa, β _b ^?1  = 193(6) GPa β _c ^?1  = 53.4(0.4) GPa). High-pressure Raman spectra show a significant splitting of several modes, which were used to identify the transformation in high-pressure high-temperature experiments in the range up to 4 GPa and 542 K. Based on the experimental series of data points determined by XRD and Raman measurements, the phase boundary of the α-to-γ-transition was determined with a Clausius–Clapeyron slope of 9.8(7) × 10^?3 GPa K^?1. An in situ measurement of the X-ray intensities was taken at 1.5 GPa and 411 K in order to identify the nature of the structural variation on increased temperatures corresponding to the previously reported transformation from α- to β-norsethite at 343 K and 1 bar. The investigations revealed, in contrast to all X-ray diffraction data recorded at 298 K, the disappearance of the superstructure reflections and the observed reflection conditions confirm the anticipated \(R\bar{3}m\) space-group symmetry. The same superstructure reflections, which disappear as temperature increases, were found to gain in intensity due to the positional shift of the Ba atoms in the γ-phase.     

High-pressure behaviour and equation of state of calcite, CaCO3

The high-pressure response of the cell parameters of calcite, CaCO₃, has been investigated by single crystal X-ray diffraction. The unit cell parameters have been refined from 0 to 1.435?GPa, and the linear and volume compressibilities have been measured as β _a =2.62(2)?×?10^?3?GPa^?1,β _c =7.94(7)?×?10^?3?GPa^?1, β _v =13.12?×?10^?3?GPa^?1. The bulk modulus has been obtained from a fit to the Birch-Murnaghan equation of state, giving K ₀=73.46?±?0.27?GPa and V ₀=367.789 ±?0.004?ų with K′=4. Combined with earlier data for magnesite, ankerite and dolomite, these data suggest that K ₀ V ₀ is a constant for the Ca-Mg rhombohedral carbonates.     

High-temperature and high-pressure equation of state for the hexagonal phase in the system NaAlSiO4 – MgAl2O4

Thermal equation of state of an Al-rich phase with Na_1.13Mg_1.51Al_4.47Si_1.62O₁₂ composition has been derived from in situ X-ray diffraction experiments using synchrotron radiation and a multianvil apparatus at pressures up to 24 GPa and temperatures up to 1,900 K. The Al-rich phase exhibited a hexagonal symmetry throughout the present pressure–temperature conditions and the refined unit-cell parameters at ambient condition were: a=8.729(1) Å, c=2.7695(5) Å, V ₀=182.77(6) Å³ (Z=1; formula weight=420.78 g/mol), yielding the zero-pressure density ρ₀=3.823(1) g/cm³ . A least-square fitting of the pressure-volume-temperature data based on Anderson’s pressure scale of gold (Anderson et al. in J Appl Phys 65:1534–543, 1989) to high-temperature Birch-Murnaghan equation of state yielded the isothermal bulk modulus K ₀=176(2) GPa, its pressure derivative K ₀ ^′ =4.9(3), temperature derivative (?K _T /?T) _P =?0.030(3) GPa K^?1 and thermal expansivity α(T)=3.36(6)×10^?5+7.2(1.9)×10^?9 T, while those values of K ₀=181.7(4) GPa, (?K _T /?T) _P =?0.020(2) GPa K^?1 and α(T)=3.28(7)×10^?5+3.0(9)×10^?9 T were obtained when K ₀ ^′ was assumed to be 4.0. The estimated bulk density of subducting MORB becomes denser with increasing depth as compared with earlier estimates (Ono et al. in Phys Chem Miner 29:527–531 2002; Vanpeteghem et al. in Phys Earth Planet Inter 138:223–230 2003; Guignot and Andrault in Phys Earth Planet Inter 143–44:107–128 2004), although the difference is insignificant (<0.6%) when the proportions of the hexagonal phase in the MORB compositions (～20%) are taken into account.     

Static elasticity of cordierite I: Effect of heavy ion irradiation on the compressibility of hydrous cordierite

The effect of ion beam irradiations on the elastic properties of hydrous cordierite was investigated by means of Raman and X-ray diffraction experiments. Oriented single crystals were exposed to swift heavy ions (Au, Bi) of various specific energies (10.0–11.1 MeV/u and 80 MeV/u), applying fluences up to 5 × 10¹³ ions/cm². The determination of unit-cell constants yields a volume strain of 3.4 × 10^?3 up to the maximum fluence, which corresponds to a compression of non-irradiated cordierite at ~480 ± 10 MPa. The unit-cell contraction is anisotropic (e ₁ = 1.4 ± 0.1 × 10^?3, e ₂ = 1.5 ± 0.1 × 10^?3, and e ₃ = 7 ± 1 × 10^?4) with the c-axis to shrink only half as much as the axes within the ab-plane. The lattice elasticity for irradiated cordierite (? = 1 × 10¹² ions/cm²) was determined from single-crystal XRD measurements in the diamond anvil cell. The fitted third-order Birch–Murnaghan equation-of-state parameters of irradiated cordierite (V ₀ = 1548.41 ± 0.16 Å³, K ₀ = 117.1 ± 1.1 GPa, ?K/?P = ?0.6 ± 0.3) reveal a 10–11 % higher compressibility compared to non-irradiated cordierite. While the higher compressibility is attributed to the previously reported irradiation-induced loss of extra-framework H₂O, the anomalous elasticity as expressed by elastic softening (β _a ^?1 , β _b ^?1 , β _c ^?1  = 397 ± 9, 395 ± 28, 308 ± 11 GPa, ?(β ^?1)/?P = ?4.5 ± 2.7, ?6.6 ± 8.4, ?5.4 ± 3.0) appears to be related to the framework stability and to be independent of the water content in the channels and thus of the ion beam exposure.     

Thermal equation of state of natural tourmaline at high pressure and temperature

Synchrotron-based in situ angle-dispersive X-ray diffraction experiments were conducted on a natural uvite-dominated tourmaline sample by using an external-heating diamond anvil cell at simultaneously high pressures and temperatures up to 18 GPa and 723 K, respectively. The angle-dispersive X-ray diffraction data reveal no indication of a structural phase transition over the P–T range of the current experiment in this study. The pressure–volume–temperature data were fitted by the high-temperature Birch–Murnaghan equation of state. Isothermal bulk modulus of K ₀ = 96.6 (9) GPa, pressure derivative of the bulk modulus of \(K_{0}^{\prime } = 12.5 \;(4)\), thermal expansion coefficient of α ₀ = 4.39 (27) × 10^?5 K^?1 and temperature derivative of the bulk modulus (?K/?T) _P  = ?0.009 (6) GPa K^?1 were obtained. The axial thermoelastic properties were also obtained with K _a0 = 139 (2) GPa, \(K_{a0}^{\prime }\) = 11.5 (7) and α _a0 = 1.00 (11) × 10^?5 K^?1 for the a-axis, and K _c0 = 59 (1) GPa, \(K_{c0}^{\prime }\) = 11.4 (5) and α _c0 = 2.41 (24) × 10^?5 K^?1 for the c-axis. Both of axial compression and thermal expansion exhibit large anisotropic behavior. Thermoelastic parameters of tourmaline in this study were also compared with that of the other two ring silicates of beryl and cordierite.     

Thermal equation of state of CaIrO<Subscript>3</Subscript> post-perovskite

The pressure–volume–temperature (P–V–T) relation of CaIrO₃ post-perovskite (ppv) was measured at pressures and temperatures up to 8.6 GPa and 1,273 K, respectively, with energy-dispersive synchrotron X-ray diffraction using a DIA-type, cubic-anvil apparatus (SAM85). Unit-cell dimensions were derived from the Le Bail full profile refinement technique, and the results were fitted using the third-order Birth-Murnaghan equation of state. The derived bulk modulus \( K_{T0} \) at ambient pressure and temperature is 168.3 ± 7.1 GPa with a pressure derivative \( K_{T0}^{\prime } \) = 5.4 ± 0.7. All of the high temperature data, combined with previous experimental data, are fitted using the high-temperature Birch-Murnaghan equation of state, the thermal pressure approach, and the Mie-Grüneisen-Debye formalism. The refined thermoelastic parameters for CaIrO₃ ppv are: temperature derivative of bulk modulus \( (\partial K_{T} /\partial T)_{P} \) = ?0.038 ± 0.011 GPa K^?1, \( \alpha K_{T} \) = 0.0039 ± 0.0001 GPa K^?1, \( \left( {\partial K_{T} /\partial T} \right)_{V} \) = ?0.012 ± 0.002 GPa K^?1, and \( \left( {\partial^{2} P/\partial T^{2} } \right)_{V} \) = 1.9 ± 0.3 × 10^?6 GPa² K^?2. Using the Mie-Grüneisen-Debye formalism, we obtain Grüneisen parameter \( \gamma_{0} \) = 0.92 ± 0.01 and its volume dependence q = 3.4 ± 0.6. The systematic variation of bulk moduli for several oxide post-perovskites can be described approximately by the relationship K _T0  = 5406.0/V(molar) + 5.9 GPa.     

Thermoelastic properties of chromium oxide Cr<Subscript>2</Subscript>O<Subscript>3</Subscript> (eskolaite) at high pressures and temperatures

A new synchrotron X-ray diffraction study of chromium oxide Cr₂O₃ (eskolaite) with the corundum-type structure has been carried out in a Kawai-type multi-anvil apparatus to pressure of 15 GPa and temperatures of 1873 K. Fitting the Birch–Murnaghan equation of state (EoS) with the present data up to 15 GPa yielded: bulk modulus (K _0,T0), 206 ± 4 GPa; its pressure derivative K′_0,T , 4.4 ± 0.8; (?K _0,T /?T) = ?0.037 ± 0.006 GPa K^?1; a = 2.98 ± 0.14 × 10^?5 K^?1 and b = 0.47 ± 0.28 × 10^?8 K^?2, where α _0,T  = a + bT is the volumetric thermal expansion coefficient. The thermal expansion of Cr₂O₃ was additionally measured at the high-temperature powder diffraction experiment at ambient pressure and α _0,T0 was determined to be 2.95 × 10^?5 K^?1. The results indicate that coefficient of the thermal expansion calculated from the EoS appeared to be high-precision because it is consistent with the data obtained at 1 atm. However, our results contradict α ₀ value suggested by Rigby et al. (Brit Ceram Trans J 45:137–148, 1946) widely used in many physical and geological databases. Fitting the Mie–Grüneisen–Debye EoS with the present ambient and high-pressure data yielded the following parameters: K _0,T0 = 205 ± 3 GPa, K′_0,T  = 4.0, Grüneisen parameter (γ ₀) = 1.42 ± 0.80, q = 1.82 ± 0.56. The thermoelastic parameters indicate that Cr₂O₃ undergoes near isotropic compression at room and high temperatures up to 15 GPa. Cr₂O₃ is shown to be stable in this pressure range and adopts the corundum-type structure. Using obtained thermoelastic parameters, we calculated the reaction boundary of knorringite formation from enstatite and eskolaite. The Clapeyron slope (with \({\text{d}}P/{\text{d}}T = - 0.014\) GPa/K) was found to be consistent with experimental data.     

P–V–T equation of state of Mn3Al2Si3O12 spessartine garnet

The thermoelastic parameters of synthetic Mn₃Al₂Si₃O₁₂ spessartine garnet were examined in situ at high pressure up to 13 GPa and high temperature up to 1,100 K, by synchrotron radiation energy dispersive X-ray diffraction within a DIA-type multi-anvil press apparatus. The analysis of room temperature data yielded K ₀ = 172 ± 4 GPa and K ₀ ^′  = 5.0 ± 0.9 when V _0,300 is fixed to 1,564.96 Å³. Fitting of P–V–T data by means of the high-temperature third-order Birch–Murnaghan EoS gives the thermoelastic parameters: K ₀ = 171 ± 4 GPa, K ₀ ^′  = 5.3 ± 0.8, (?K _0,T /?T) _P  = ?0.049 ± 0.007 GPa K^?1, a ₀ = 1.59 ± 0.33 × 10^?5 K^?1 and b ₀ = 2.91 ± 0.69 × 10^?8 K^?2 (e.g., α _0,300 = 2.46 ± 0.54 × 10^?5 K^?1). Comparison with thermoelastic properties of other garnet end-members indicated that the compression mechanism of spessartine might be the same as almandine and pyrope but differs from that of grossular. On the other hand, at high temperature, spessartine softens substantially faster than pyrope and grossular. Such softening, which is also reported for almandine, emphasize the importance of the cation in the dodecahedral site on the thermoelastic properties of aluminosilicate garnet.     

Pressure dependence of structural parameters of paragonite

The compressibility and structure of a 2M₁ paragonite with composition [Na_0.88K_0.10Ca_0.01Ba_0.01] [Al_1.97Ti_0.007Fe_0.01Mn_0.002Mg_0.006]Si_3.01Al_0.99O₁₀OH₂ were determined at pressures between 1 bar and 41 kbar, by single crystal X-ray diffraction using a Merrill-Bassett diamond anvil cell. Compressibility turned out to be largely anisotropic, linear compressibility coefficients parallel to the unit cell edges being β_a=3.5(1)·10^?4, β_b=3.6(1)·10^?4, β_c=8.3(3)·10^?4 kbar^?1 (β_a:β_b·β_c=1:1028:2.371). The isothermal bulk modulus, calculated as the reciprocal of the mean compressibility of the cell volume, was 650(20) kbar. The main features of the deformation mechanism resulting from structural refinements at pressures of 0.5, 25.4, 40.5 kbar were: –?variation in sheet thickness, showing that compression of the c parameter was mainly due to the interlayer thickness reduction from 3.07 Å at 0.5 kbar to 2.81 Å at 40.5 kbar; –?the compressibility of octahedra was greater than that of tetrahedra, the dimensional misfit between tetrahedral and octahedral sheets increased with P, so that tetrahedral rotation angel α increased from 15° at 0.5 kbar to 21.6° at 40.5 kbar; –?the basal surface corrugation (Δz) of the tetrahedral layer, due to the different dimensions of M1 and M2 octahedra and to the octahedral distortion, decreased with P (Δz=0.19 and 0.12 Å at 0.5 and 40.5 kbar respectively). Comparison of the new data on paragonite with those of a K-muscovite and a Na-rich muscovite (Comodi and Zanazzi 1995) revealed a clear trend toward decreasing of compressibility when Na substitutes for K atoms in the interlayer sites.     

Hydrostatic compression of γ-Mg2SiO4 to mantle pressures and 700 K: Thermal equation of state and related thermoelastic properties

P-V-T equations of state for the γ phase of Mg₂SiO₄ have been fitted to unit cell volumes measured under simultaneous high pressure (up 30 GPa) and high temperature (up to 700 K) conditions. The measurements were conducted in an externally heated diamond anvil cell using synchrotron x-ray diffraction. Neon was used as a pressure medium to provide a more hydrostatic pressure environment. The P-V-T data include 300 K-isothermal compression to 30 GPa, 700 K-compression to 25 GPa and some additional data in P-T space in the region 15 to 30 GPa and 300 to 700 K. The isothermal bulk modulus and its pressure derivative, determined from the isothermal compression data, are 182(3) GPa and 4.2(0.3) at T=300 K, and 171(4) GPa and 4.4(0.5) at T=700 K. Fitting all the P-V-T data to a high-temperature Murnaghan equation of state yields: K _TO=182(3.0) GPa, K _TO=4.0(0.3), ?K _T /?T)₀=?2.7(0.5)×10^?2 GPa/K and (?² K _T /?P?T)₀=5.5(5.2)×10^?4/K at the ambient condition.     

A high P-T single-crystal X-ray diffraction study of thermoelasticity of MgSiO3 orthoenstatite

P-V-T data of MgSiO₃ orthoenstatite have been measured by single-crystal X-ray diffraction at simultaneous high pressures (in excess of 4.5 GPa) and temperatures (up to 1000 K). The new P-V-T data of the orthoenstatite, together with previous compression data and thermal expansion data, are described by a modified Birch-Murnaghan equation of state for diverse temperatures. The fitted thermoelastic parameters for MgSiO₃ orthoenstatite are: thermal expansion ?α/?P with values of a=2.86(29)×10^-5 K^-1 and b=0.72(16)×10^-8 K^-2; isothermal bulk modulus K _{T o} =102.8(2) GPa; pressure derivative of bulk modulus K′=?K/?P=10.2(1.2); and temperature derivative of bulk modulus K=?K/?T=-0.037(5) GPa/K. The derived thermal Grüneisen parameter is γ _th=1.05 for ambient conditions; Anderson-Grüneisen parameter is δ _{T o} =11.6, and the pressure derivative of thermal expansion is ?α/?P=-3.5×10^-6K^-1 GPa^-1. From the P-V-T data and the thermoelastic equation of state, thermal expansions at two constant pressures of 1.5 GPa and 4.0 GPa are calculated. The resulting pressure dependence of thermal expansion is Δα/ΔP=-3.2(1)× 10^-6 K^-1 GPa^-1. The significantly large values of K′, K, δ _T and ?α/?P indicate that compression/expansion of MgSiO₃ orthoenstatite is very sensitive to changes of pressure and temperature.     

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.     

Helium measurements of pore fluids obtained from the San Andreas Fault Observatory at Depth (SAFOD, USA) drill cores

⁴He accumulated in fluids is a well established geochemical tracer used to study crustal fluid dynamics. Direct fluid samples are not always collectable; therefore, a method to extract rare gases from matrix fluids of whole rocks by diffusion has been adapted. Helium was measured on matrix fluids extracted from sandstones and mudstones recovered during the San Andreas Fault Observatory at Depth (SAFOD) drilling in California, USA. Samples were typically collected as subcores or from drillcore fragments. Helium concentration and isotope ratios were measured 4?C6 times on each sample, and indicate a bulk ⁴He diffusion coefficient of 3.5?±?1.3?×?10^?C8 cm²?s^?C1 at 21°C, compared to previously published diffusion coefficients of 1.2?×?10^?C18 cm²?s^?C1 (21°C) to 3.0?×?10^?C15 cm²?s^?C1 (150°C) in the sands and clays. Correcting the diffusion coefficient of ⁴He_water for matrix porosity (??3%) and tortuosity (??6?C13) produces effective diffusion coefficients of 1?×?10^?C8 cm^2?s^?C1 (21°C) and 1?×?10^?C7 (120°C), effectively isolating pore fluid ⁴He from the ⁴He contained in the rock matrix. Model calculations indicate that <6% of helium initially dissolved in pore fluids was lost during the sampling process. Complete and quantitative extraction of the pore fluids provide minimum in situ porosity values for sandstones 2.8?±?0.4% (SD, n?=?4) and mudstones 3.1?±?0.8% (SD, n?=?4).     

Single-crystal elastic properties of alunite, KAl3(SO4)2(OH)6

The single-crystal elastic constants of natural alunite (ideally KAl₃(SO₄)₂(OH)₆) were determined by Brillouin spectroscopy. Chemical analysis by electron microprobe gave a formula KAl₃(SO₄)₂(OH)₆. Single crystal X-ray diffraction refinement with R ₁ = 0.0299 for the unique observed reflections (|F _o| > 4σ _F) and wR ₂ = 0.0698 for all data gave a = 6.9741(3) Å, c = 17.190(2) Å, fractional positions and thermal factors for all atoms. The elastic constants (in GPa), obtained by fitting the spectroscopic data, are C ₁₁ = 181.9 ± 0.3, C ₃₃ = 66.8 ± 0.8, C ₄₄ = 42.8 ± 0.2, C ₁₂ = 48.2 ± 0.5, C ₁₃ = 27.1 ± 1.0, C ₁₄ = 5.4 ± 0.5, and C ₆₆ = ½(C ₁₁–C ₁₂) = 66.9 ± 0.3 GPa. The VRH averages of bulk and shear modulus are 63 and 49 GPa, respectively. The aggregate Poisson ratio is 0.19. The high value of the ratio C ₁₁/C ₃₃ = 2.7 and of the ratio C ₆₆/C ₄₄ = 1.6 are characteristic of an anisotropic structure with very weak interlayer interactions along the c-axis. The basal plane (001) is characterized by 0.1% longitudinal acoustic anisotropy and 0.9–1.1% shear acoustic anisotropy, which gives alunite a characteristic pseudo-hexagonal elastic behavior, and is related to the pseudo-hexagonal arrangement of the Al(O,OH)₆ octahedra in the basal layer. The elastic Debye temperature of alunite is 654 K. The large discrepancy between the elastic and heat capacity Debye temperature is also a consequence of the layered structure.