Thermal expansion of gehlenite, Ca2Al[AlSiO7], and the related aluminates LnCaAl[Al2O7] with Ln = Tb, Sm

The thermal expansion of gehlenite, Ca₂Al[AlSiO₇], (up to T=830 K), TbCaAl[Al₂O₇] (up to T=1100 K) and SmCaAl[Al₂O₇] (up to T=1024 K) has been determined. All compounds are of the melilite structure type with space group Thermal expansion data were obtained from in situ X-ray powder diffraction experiments in-house and at HASYLAB at the Deutsches Elektronen Synchrotron (DESY) in Hamburg (Germany). The thermal expansion coefficients for gehlenite were found to be: α₁=7.2(4)×10⁻⁶×K⁻¹+3.6(7)×10⁻⁹ΔT×K⁻² and α₃=15.0(1)×10⁻⁶×K⁻¹. For TbCaAl[Al₂O₇] the respective values are: α₁=7.0(2)×10⁻⁶×K⁻¹+2.0(2)×10⁻⁹ΔT×K⁻² and α₃=8.5(2)×10⁻⁶×K⁻¹+2.0(3)×10⁻⁹ΔT×K⁻², and the thermal expansion coefficients for SmCaAl[Al₂O₇] are: α₁=6.9(2)×10⁻⁶×K⁻¹+1.7(2)×10⁻⁹ΔT×K⁻² and α₃=9.344(5)×10⁻⁶×K⁻¹. The expansion mechanisms of the three compounds are explained in terms of structural trends obtained from Rietveld refinements of the crystal structures of the compounds against the powder diffraction patterns. No structural phase transitions have been observed. While gehlenite behaves like a ‘proper’ layer structure, the aluminates show increased framework structure behavior. This is most probably explained by stronger coulombic interactions between the tetrahedral conformation and the layer-bridging cations due to the coupled substitution (Ca²⁺+Si⁴⁺)–(Ln ³⁺+Al³⁺) in the melilite-type structure. This article has been mistakenly published twice. The first and original version of it is available at .     

The thermal behaviour of γ-CaSO<Subscript>4</Subscript>

Thermal behaviour of γ-anhydrite (γ-CaSO₄, soluble anhydrite) has been investigated in situ real-time using laboratory parallel-beam X-ray powder diffraction data. Thermal expansion has been analysed from 303 to 569 K with temperature steps of 4 K. Lattice parameters and volume were fitted with a second-order polynomial to calculate thermal expansion coefficients. Thermal expansion of γ-anhydrite is anisotropic being larger along the c axis. Within the 343–383 K thermal range, γ-anhydrite has been found to partially re-hydrate to bassanite CaSO₄·0.5H₂O. At 455 K the transformation γ-CaSO₄ → β-CaSO₄, insoluble anhydrite, starts reaching completion at 653 K.     

Thermal expansion of gehlenite, Ca2Al[AlSiO7], and the related aluminates LnCaAl[Al2O7] with Ln=Tb, Sm

The thermal expansion of gehlenite, Ca₂Al[AlSiO₇], (up to T=830 K), TbCaAl[Al₂O₇] (up to T=1,100 K) and SmCaAl[Al₂O₇] (up to T=1,024 K) has been determined. All compounds are of the melilite structure type with space group Thermal expansion data was obtained from in situ X-ray powder diffraction experiments in-house and at HASYLAB at the Deutsches Elektronen Synchrotron (DESY) in Hamburg (Germany). The thermal expansion coefficients for gehlenite were found to be: α₁=7.2(4)×10⁻⁶ K⁻¹+3.6(7)×10⁻⁹ΔT K⁻² and α₃=15.0(1)×10⁻⁶ K⁻¹. For TbCaAl[Al₂O₇] the respective values are: α₁=7.0(2)×10⁻⁶ K⁻¹+2.0(2)×10⁻⁹ΔT K⁻² and α₃=8.5(2)×10⁻⁶ K⁻¹+2.0(3)×10⁻⁹ΔT K⁻², and the thermal expansion coefficients for SmCaAl[Al₂O₇] are: α₁=6.9(2)× 10⁻⁶ K⁻¹+1.7(2)×10⁻⁹ΔT K⁻² and α₃=9.344(5)×10⁻⁶ K⁻¹. The expansion-mechanisms of the three compounds are explained in terms of structural trends obtained from Rietveld refinements of the crystal structures of the compounds against the powder diffraction patterns. No structural phase transitions have been observed. While gehlenite behaves like a ’proper’ layer structure, the aluminates show increased framework structure behaviour. This is most probably explained by stronger coulombic interactions between the tetrahedral conformation and the layer-bridging cations due to the coupled substitution (Ca²⁺+Si⁴⁺)-(Ln ³⁺+Al³⁺) in the melilite-type structure. Electronic Supplementary Material Supplementary material is available for this article at     

Low T heat capacity measurements and new entropy data for titanite (sphene): implications for thermobarometry of high-pressure rocks

The accepted standard state entropy of titanite (sphene) has been questioned in several recent studies, which suggested a revision from the literature value 129.3 ± 0.8 J/mol K to values in the range of 110–120 J/mol K. The heat capacity of titanite was therefore re-measured with a PPMS in the range 5 to 300 K and the standard entropy of titanite was calculated as 127.2 ± 0.2 J/mol K, much closer to the original data than the suggested revisions. Volume parameters for a modified Murgnahan equation of state: V _P,T  = V ₂₉₈° × [1 + a°(T − 298) − 20a°(T − 298)] × [1 – 4P/(K ₂₉₈ × (1 – 1.5 × 10⁻⁴ [T − 298]) + 4P)]^1/4 were fit to recent unit cell determinations at elevated pressures and temperatures, yielding the constants V ₂₉₈° = 5.568 J/bar, a° = 3.1 × 10⁻⁵ K⁻¹, and K = 1,100 kbar. The standard Gibbs free energy of formation of titanite, −2456.2 kJ/mol (∆H°_f = −2598.4 kJ/mol) was calculated from the new entropy and volume data combined with data from experimental reversals on the reaction, titanite + kyanite = anorthite + rutile. This value is 4–11 kJ/mol less negative than that obtained from experimental determinations of the enthalpy of formation, and it is slightly more negative than values given in internally consistent databases. The displacement of most calculated phase equilibria involving titanite is not large except for reactions with small ∆S. Re-calculated baric estimates for several metamorphic suites yield pressure differences on the order of 2 kbar in eclogites and 10 kbar for ultra-high pressure titanite-bearing assemblages.     

Excess heat capacity and entropy of mixing along the chlorapatite–fluorapatite binary join

The heat capacity at constant pressure, C _p, of chlorapatite [Ca₅(PO₄)₃Cl – ClAp], and fluorapatite [Ca₅(PO₄)₃F – FAp], as well as of 12 compositions along the chlorapatite–fluorapatite join have been measured using relaxation calorimetry [heat capacity option of the physical properties measurement system (PPMS)] and differential scanning calorimetry (DSC) in the temperature range 5–764 K. The chlor-fluorapatites were synthesized at 1,375–1,220°C from Ca₃(PO₄)₂ using the CaF₂–CaCl₂ flux method. Most of the chlor-fluorapatite compositions could be measured directly as single crystals using the PPMS such that they were attached to the sample platform of the calorimeter by a crystal face. However, the crystals were too small for the crystal face to be polished. In such cases, where the sample coupling was not optimal, an empirical procedure was developed to smoothly connect the PPMS to the DSC heat capacities around ambient T. The heat capacity of the end-members above 298 K can be represented by the polynomials: C _p^ClAp = 613.21 − 2,313.90T ^−0.5 − 1.87964 × 10⁷ T ⁻² + 2.79925 × 10⁹ T ⁻³ and C _p^FAp = 681.24 − 4,621.73 × T ^−0.5 − 6.38134 × 10⁶ T ⁻² + 7.38088 × 10⁸ T ⁻³ (units, J mol⁻¹ K⁻¹). Their standard third-law entropy, derived from the low-temperature heat capacity measurements, is S° = 400.6 ± 1.6 J mol⁻¹ K⁻¹ for chlorapatite and S° = 383.2 ± 1.5 J mol⁻¹ K⁻¹ for fluorapatite. Positive excess heat capacities of mixing, ΔC _p^ex, occur in the chlorapatite–fluorapatite solid solution around 80 K (and to a lesser degree at 200 K) and are asymmetrically distributed over the join reaching a maximum of 1.3 ± 0.3 J mol⁻¹ K⁻¹ for F-rich compositions. They are significant at these conditions exceeding the 2σ-uncertainty of the data. The excess entropy of mixing, ΔS ^ex, at 298 K reaches positive values of 3–4 J mol⁻¹ K⁻¹ in the F-rich portion of the binary, is, however, not significantly different from zero across the join within its 2σ-uncertainty.     

The dehydration process of gypsum under high pressure

The effects of pressure on the dehydration of gypsum materials were investigated up to 633 K and 25 GPa by using Raman spectroscopy and synchrotron X-ray diffraction with an externally heated diamond anvil cell. At 2.5 GPa, gypsum starts to dehydrate around 428 K, by forming bassanite, CaSO₄ hemihydrate, which completely dehydrates to γ-anhydrite at 488 K. All the sulphate modes decrease linearly between 293 and 427 K with temperature coefficients ranging from −0.119 to −0.021 cm⁻¹ K⁻¹, where an abrupt change in the ν₃ mode and in the OH-stretching region indicates the beginning of dehydration. Increasing the temperature to 488 K, the OH-stretching modes completely disappear, marking the complete dehydration and formation of γ-anhydrite. Moreover, the sample changes from transparent to opaque to transparent again during the dehydration sequence gypsum-bassanite-γ-anhydrite, which irreversibly transforms to β-anhydrite form at 593 K. These data compared with the dehydration temperature at room pressure indicate that the dehydration temperature increases with pressure with a ΔP/ΔT slope equal to 230 bar/K. Synchrotron X-ray diffraction experiments show similar values of temperature and pressure for the first appearance of bassanite. Evidence of phase transition from β-anhydrite structure to the monazite type was observed at about 2 GPa under cold compression. On the other hand at the same pressure (2 GPa and 633 K), β-anhydrite was found, indicating a positive Clausis-Clayperon slope of the transition. This transformation is completely reversible as showed by the Raman spectra on the sample recovered after phase transition.     

Interaction of Freshly Precipitated Silica Gel with Aqueous Silicic Acid Solutions under Ambient and Near Neutral pH-conditions: A Detailed Analysis of Linear Rate Law

Interaction of freshly precipitated silica gel with aqueous solutions was studied at laboratory batch experiments under ambient and near neutral pH-conditions. The overall process showed excellent reversibility: gel growth could be considered as an opposite process to dissolution and a linear rate law could be applied to experimental data. Depending on the used rate law form, the resulting rate constants were sensitive to errors in parameters/variables such as gel surface area, equilibrium constants, Si-fluxes, and reaction quotients. The application of an Integrated Exponential Model appeared to be the best approach for dissolution data evaluation. It yielded the rate constants k _dissol ∼ (4.50 ± 0.68) × 10⁻¹² and k _growth ∼ (2.58 ± 0.39) × 10⁻⁹ mol m⁻² s⁻¹ for zero ionic strength. In contrast, a Differential Model gave best results for growth data modeling. It yielded the rate constants k _dissol ∼ (1.14 ± 0.44) × 10⁻¹¹ and k _growth ∼ (6.08 ± 2.37) × 10⁻⁹ mol m⁻² s⁻¹ for higher ionic strength (I ∼ 0.04 to 0.11 mol L⁻¹). The found silica gel solubility at zero ionic strength was somewhat lower than the generally accepted value. Based on the and standard Gibbs free energy of silica gel formation was calculated as and −850,318 ± 20 J mol⁻¹, respectively. Activation energies for silica gel dissolution and growth were determined as and respectively. An universal value for growth of any silica polymorph, is not consistent with the value for silica gel growth, which questions the hypothesis about one unique activated complex controlling the silica polymorph growth.     

Thermoelasticity and high-<Emphasis Type="Italic">T</Emphasis> behaviour of anthophyllite

The thermoelastic behaviour of anthophyllite has been determined for a natural crystal with crystal-chemical formula ^ANa_0.01 ^B(Mg_1.30Mn_0.57Ca_0.09Na_0.04) ^C(Mg_4.95Fe_0.02Al_0.03) ^T(Si_8.00)O₂₂ ^W(OH)₂ using single-crystal X-ray diffraction to 973 K. The best model for fitting the thermal expansion data is that of Berman (J Petrol 29:445–522, 1988) in which the coefficient of volume thermal expansion varies linearly with T as α _V,T  = a ₁ + 2a ₂ (T − T ₀): α₂₉₈ = a ₁ = 3.40(6) × 10⁻⁵ K⁻¹, a ₂ = 5.1(1.0) × 10⁻⁹ K⁻². The corresponding axial thermal expansion coefficients for this linear model are: α _a ,₂₉₈ = 1.21(2) × 10⁻⁵ K⁻¹, a _2,a  = 5.2(4) × 10⁻⁹ K⁻²; α _b ,₂₉₈ = 9.2(1) × 10⁻⁶ K⁻¹, a _2,b  = 7(2) × 10⁻¹⁰ K⁻². α _c ,₂₉₈ = 1.26(3) × 10⁻⁵ K⁻¹, a _2,c  = 1.3(6) × 10⁻⁹ K⁻². The thermoelastic behaviour of anthophyllite differs from that of most monoclinic (C2/m) amphiboles: (a) the ε ₁ − ε ₂ plane of the unit-strain ellipsoid, which is normal to b in anthophyllite but usually at a high angle to c in monoclinic amphiboles; (b) the strain components are ε ₁ ≫ ε ₂ > ε ₃ in anthophyllite, but ε ₁ ~ ε ₂ ≫ ε ₃ in monoclinic amphiboles. The strain behaviour of anthophyllite is similar to that of synthetic C2/m ^ANa ^B(LiMg) ^CMg₅ ^TSi₈ O₂₂ ^W(OH)₂, suggesting that high contents of small cations at the B-site may be primarily responsible for the much higher thermal expansion ⊥(100). Refined values for site-scattering at M4 decrease from 31.64 epfu at 298 K to 30.81 epfu at 973 K, which couples with similar increases of those of M1 and M2 sites. These changes in site scattering are interpreted in terms of Mn ↔ Mg exchange involving M1,2 ↔ M4, which was first detected at 673 K.     

Temperature evolution between 50 K and 320 K of the thermal expansion tensor of gypsum derived from neutron powder diffraction data

The unit cell parameters, extracted from Rietveld analysis of neutron powder diffraction data collected between 4.2 K and 320 K, have been used to calculate the temperature evolution of the thermal expansion tensor for gypsum for 50 ≤ T ≤ 320 K. At 300 K the magnitudes of the principal axes are α ₁₁  = 1.2(6) × 10⁻⁶ K⁻¹, α ₂₂  = 36.82(1) × 10⁻⁶ K⁻¹ and α ₃₃  = 25.1(5) × 10⁻⁶ K⁻¹. The maximum axis, α ₂₂ , is parallel to b, and using Institution of Radio Engineers (IRE) convention for the tensor orthonormal basis, the axes α ₁₁ and α ₃₃ have directions equal to (−0.979, 0, 0.201) and (0.201, 0, 0.979) respectively. The orientation and temperature dependent behaviour of the thermal expansion tensor is related to the crystal structure in the I2/a setting. Received 12 February 1998 / Revised, accepted 19 October 1998     

Trace element partitioning between orthopyroxene and anhydrous silicate melt on the lherzolite solidus from 1.1 to 3.2 GPa and 1,230 to 1,535°C in the model system Na<Subscript>2</Subscript>O–CaO–MgO–Al<Subscript>2</Subscript>O<Subscript>3</Subscript>–SiO<Subscript>2</Subscript>

We determined experimentally the Nernst distribution coefficient between orthopyroxene and anhydrous silicate melt for trace elements i in the system Na₂O–CaO–MgO–Al₂O₃–SiO₂ (NCMAS) along the dry model lherzolite solidus from 1.1 GPa/1,230°C up to 3.2 GPa/1,535°C in a piston cylinder apparatus. Major and trace element composition of melt and orthopyroxene were determined with a combination of electron microprobe and ion probe analyses. We provide partitioning data for trace elements Li, Be, B, K, Sc, Ti, V, Cr, Co, Ni, Rb, Sr, Y, Zr, Nb, Cs, Ba, La, Ce, Sm, Nd, Yb, Lu, Hf, Ta, Pb, U, and Th. The melts were chosen to be boninitic at 1.1 and 2.0 GPa, picritic at 2.3 GPa and komatiitic at 2.7 and 3.2 GPa. Orthopyroxene is Tschermakitic with 8 mol% Mg-Tschermaks MgAl[AlSiO₆] at 1.1 GPa while at higher pressure it has 18–20 mol%. The rare earth elements show a continuous, significant increase in compatibility with decreasing ionic radius from D _La^opx−melt ∼ 0.0008 to D _Lu^opx−melt ∼ 0.15. For the high-field-strength elements compatibility increases from D _Th^opx−melt ∼ 0.001 through D _Nb^opx−melt ∼ 0.0015, D _U^opx−melt ∼ 0.002, D _Ta^opx−melt ∼ 0.005, D _Zr^opx−melt ∼ 0.02 and D _Hf^opx−melt ∼ 0.04 to D _Ti^opx−melt ∼ 0.14. From mathematical and graphical fits we determined best-fit values for D ₀^M1, D ₀^M2, r ₀^M1, r ₀^M2, E ₀^M1, and E ₀^M2 for the two different M sites in orthopyroxene according to the lattice strain model and calculated the intracrystalline distribution between M1 and M2. Our data indicate extreme intracrystalline fractionation for most elements in orthopyroxene; for the divalent cations D _i ^M2−M1 varies by three orders of magnitude between D _Co^M2−M1 = 0.00098–0.00919 and D _Ba^M2−M1 = 2.3–28. Trivalent cations Al and Cr almost exclusively substitute on M1 while the other trivalent cations substitute on M2; D _La^M2−M1 reaches extreme values between 6.5 × 10⁷ and 1.4 × 10¹⁶. Tetravalent cations Ti, Hf, and Zr almost exclusively substitute on M1 while U and Th exclusively substitute on M2. Our new comprehensive data set can be used for polybaric-polythermal melting models along the Earth’s mantle solidus. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Thermal expansion behavior of Ce2Zr2O7 up to 898 K in conjunction with structural analyses by neutron diffraction

The thermal expansion of cubic pyrochlore Ce₂Zr₂O₇ has been measured from room temperature to 898 K on polycrystalline material in conjunction with structural analyses using neutron diffraction. This compound has a thermal expansion coefficient in line with the other comparable lanthanoide pyrochlore oxides. The coefficient can be expressed as α(T) = 8.418 × 10⁻⁶ + 0.9861 × 10⁻⁹ × T. The structural refinements performed for each measured temperature showed a comparable linear evolution of the Ce–O/Zr–O distances (within 0.57%).     

Kinetics of heterogeneous bubble nucleation in rhyolitic melts: implications for the number density of bubbles in volcanic conduits and for pumice textures

We performed decompression experiments to simulate the ascent of a phenocryst-bearing rhyolitic magma in a volcanic conduit. The starting materials were bubble-free rhyolites water-saturated at 200 MPa–800°C under oxidizing conditions: they contained 6.0 wt% dissolved H₂O and a dense population of hematite crystals (8.7 ± 2 × 10⁵ mm⁻³). Pressure was decreased from the saturation value to a final value ranging from 99 to 20 MPa, at constant temperature (800°C); the rate of decompression was either 1,000 or 27.8 kPa/s. In all experiments, we observed a single event of heterogeneous bubble nucleation beginning at a pressure P _N equal to 63 ± 3 MPa in the 1,000 kPa/s series, and to 69 ± 1 MPa in the 27.8 kPa/s series. Below P _N, the degree of water supersaturation in the liquid rapidly decreased to a few 0.1 wt%, the nucleation rate dropped, and the bubble number density (BND) stabilized to a value strongly sensitive to decompression rate: 80 mm⁻³ at 27.8 kPa/s vs. 5,900 mm⁻³ at 1,000 kPa/s. This behaviour is like the behavior formerly described in the case of homogeneous bubble nucleation in the rhyolite-H₂O system and in numerical simulations of vesiculation in ascending magmas. Similar degrees of water supersaturation were measured at 27.8 and 1,000 kPa/s, implying that a faster decompression rate does not result in a larger departure from equilibrium. Our experimental results imply that BNDs in acid to intermediate magmas ascending in volcanic conduits will depend on both the decompression rate and the number density of phenocrysts, especially the number density of magnetite microphenocrysts (1–100 mm⁻³), which is the only mineral species able to reduce significantly the degree of water supersaturation required for bubble nucleation. Very low BNDs (≈1 mm⁻³) are predicted in the case of effusive eruptions ( ≈ 0.1 kPa/s). High BNDs (up to 10⁷ mm⁻³) and bimodal bubble size distributions are expected in the case of explosive eruptions: (1) a relatively small number density of bubbles (1–100 mm⁻³) will first nucleate in the lower part of the conduit ( ≈ 10 kPa/s), either at high pressure on magnetite or at lower pressure on quartz and feldspar (or by homogeneous nucleation in the liquid) and (2) then, extreme decompression rates near the fragmentation level ( ≈ 10³ kPa/s) will trigger a major nucleation event leading to the multitude of small bubbles, typically a few micrometers to a few tens of micrometers in diameter, which characterizes most silicic pumices.     

Behavior of epidote at high pressure and high temperature: a powder diffraction study up to 10 GPa and 1,200 K

The thermo-elastic behavior of a natural epidote [Ca_1.925 Fe_0.745Al_2.265Ti_0.004Si_3.037O₁₂(OH)] has been investigated up to 1,200 K (at 0.0001 GPa) and 10 GPa (at 298 K) by means of in situ synchrotron powder diffraction. No phase transition has been observed within the temperature and pressure range investigated. P–V data fitted with a third-order Birch–Murnaghan equation of state (BM-EoS) give V ₀ = 458.8(1)ų, K _T0 = 111(3) GPa, and K′ = 7.6(7). The confidence ellipse from the variance–covariance matrix of K _T0 and K′ from the least-square procedure is strongly elongated with negative slope. The evolution of the “Eulerian finite strain” vs “normalized stress” yields Fe(0) = 114(1) GPa as intercept values, and the slope of the regression line gives K′ = 7.0(4). The evolution of the lattice parameters with pressure is slightly anisotropic. The elastic parameters calculated with a linearized BM-EoS are: a ₀ = 8.8877(7) Å, K _T0(a) = 117(2) GPa, and K′(a) = 3.7(4) for the a-axis; b ₀ = 5.6271(7) Å, K _T0(b) = 126(3) GPa, and K′(b) = 12(1) for the b-axis; and c ₀ = 10.1527(7) Å, K _T0(c) = 90(1) GPa, and K’(c) = 8.1(4) for the c-axis [K _T0(a):K _T0(b):K _T0(c) = 1.30:1.40:1]. The β angle decreases with pressure, β_P(°) = β_P0 −0.0286(9)P +0.00134(9)P ² (P in GPa). The evolution of axial and volume thermal expansion coefficient, α, with T was described by the polynomial function: α(T) = α₀ + α₁ T ^−1/2. The refined parameters for epidote are: α₀ = 5.1(2) × 10⁻⁵ K⁻¹ and α₁ = −5.1(6) × 10⁻⁴ K^1/2 for the unit-cell volume, α₀(a) = 1.21(7) × 10⁻⁵ K⁻¹ and α₁(a) = −1.2(2) × 10⁻⁴ K^1/2 for the a-axis, α₀(b) = 1.88(7) × 10⁻⁵ K⁻¹ and α₁(b) = −1.7(2) × 10⁻⁴ K^1/2 for the b-axis, and α₀(c) = 2.14(9) × 10⁻⁵ K⁻¹ and α₁(c) = −2.0(2) × 10⁻⁴ K^1/2 for the c-axis. The thermo-elastic anisotropy can be described, at a first approximation, by α₀(a): α₀(b): α₀(c) = 1 : 1.55 : 1.77. The β angle increases continuously with T, with β_T(°) = β_T0 + 2.5(1) × 10⁻⁴ T + 1.3(7) × 10⁻⁸ T ². A comparison between the thermo-elastic parameters of epidote and clinozoisite is carried out.     

Thermal expansion and dehydroxylation of phengite micas

Phengite samples (2M ₁ and 3T politypes) and a synthetic end-member muscovite specimen were studied by in situ high-temperature synchrotron radiation X-ray diffraction. The measured volume thermal expansion of 2M ₁ phengite (<α _V> ≈ 36.6 × 10⁻⁶ K⁻¹) was systematically greater than <α _V> of the 3T polytype (≈33.3 × 10⁻⁶ K⁻¹). A positive linear correlation between the average thermal expansion on (001) plane and the mean tetrahedral rotation angle at ambient condition is proposed on the ground of new measurements and literature data. Dehydroxylation processes were observed in 2M ₁, starting at 1,000 K in 3T at 800 and 945 K in synthetic muscovite. Rietveld refinements allowed a determination of structural variations upon heating of phengite samples and their dehydroxylate phases. The phengite structure expands by regularizing the tetrahedral sheet and by reducing the bond length differences between the outer and inner coordination shell of the interlayer site. The dehydroxylate phase derived from 2M ₁ is characterized by fivefold polyhedra in the low temperature form as a consequence of two OH groups reacting to form H₂O + O (residual). The dehydroxylate exhibits an increase of the cation–cation distances along the M–Or–M bonds with respect to low-temperature phengite structures. For the 3T phase, we were unable to achieve completion of dehydroxylation. The refined structural model of the dehydroxylate phase shows two hydroxyl sites, but at a short distance from one another. This result suggests that the dehydroxylation reaction did not proceed to completion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Optical spectroscopic study of synthetic NaScSi<Subscript>2</Subscript>O<Subscript>6</Subscript>–CaNiSi<Subscript>2</Subscript>O<Subscript>6</Subscript> pyroxenes at normal and high pressures

Six synthetic NaScSi₂O₆–CaNiSi₂O₆ pyroxenes were studied by optical absorption spectroscopy. Five of them of intermediate (Na_1−x , Ca _x )(Sc_1−x , Ni _x )Si₂O₆ compositions show spectra typical of Ni²⁺ in octahedral coordination, more precise Ni²⁺ at the M1 site of the pyroxene structure. The common feature of all spectra is three broad absorption bands with maxima around 8,000, 13,000 and 24,000 cm⁻¹ assigned to ³ A _2g → ³ T _2g, ³ A _2g → ³ T _1g and →³ T _1g (³ P) electronic spin-allowed transitions of ^VINi²⁺. A weak narrow peak at ∼14,400 cm⁻¹ is assigned to the spin-forbidden ³ A _2g → ¹ T _2g (¹ D) transition of Ni²⁺. Under pressure the spin-allowed bands shift to higher energies and change in intensity. The octahedral compression modulus, calculated from the shift of the ³ A _2g → ³ T _2g band in the (Na_0.7Ca_0.3)(Sc_0.7Ni_0.3)Si₂O₆ pyroxene is evaluated as 85±20 GPa. The Racah parameter B of Ni²⁺(M1) is found gradually changing from ∼919 cm⁻¹ at ambient pressure to ∼890 cm⁻¹ at 6.18 GPa. The Ni end-member pyroxene [(Ca_0.93 Ni_0.07)NiSi₂O₆] has a spectrum different from all others. In addition to the above mentioned bands of Ni²⁺(M1) it displays several new relatively intense and broad extra bands, which were attributed to electronic transitions of Ni²⁺ at the M2 site. In difference to CaO₈ polyhedron geometry of an eightfold coordination, Ni²⁺(M2)O₈ polyhedra are assumed to be relatively large distorted octahedra. Due to different distortions and different compressibilities of the M1 and M2 sites the Ni²⁺(M1)- and Ni²⁺(M2)-bands display rather different pressure-induced behaviors, becoming more resolved in the high-pressure spectra than in that measured at atmospheric pressure. The octahedral compression modulus of Ni²⁺(M1) in this end-member pyroxene is evaluated as 150 ± 25 GPa, which is noticeably larger than in Ni_0.3 pyroxene. This is due to a smaller size and, thus, a stiffer character of Ni²⁺(M1)O₆ octahedron in the (Ca_0.93Ni_0.07)NiSi₂O₆ pyroxene compared to (Na_0.7Ca_0.3)(Sc_0.7Ni_0.3)Si₂O₆.

Single-crystal EPR and DFT studies of a [BO<Subscript>4</Subscript>]<Superscript>0</Superscript> center in datolite: electronic structure,formation mechanism and implications

A natural datolite CaBSiO₄(OH) (Bergen Hill, NJ, USA), before and after gamma-ray irradiation (up to ~70 kGy), has been investigated by single-crystal and powder electron paramagnetic resonance (EPR) spectroscopy from 10 to 295 K. EPR spectra of gamma-ray-irradiated datolite show the presence of a boron-associated oxygen hole center (BOHC) and an atomic hydrogen center (H⁰), both of which grow with increasing radiation dose. The principal g and A(¹¹B) values of the BOHC at 10 K are: g ₁ = 2.04817(3), g ₂ = 2.01179(2), g ₃ = 2.00310(2), A ₁ = −0.401(7) mT, A ₂ = −0.906(2) mT, A ₃ = −0.985(2) mT, with the orientations of the g ₁ and A ₁ axes approximately along the B–OH bond direction. These experimental results suggest that the BOHC represents hole trapping on the hydroxyl oxygen atom after the removal of the proton (i.e. a [BO₄]⁰ center): via a reaction O₃BOH → O₃BO^· + H⁰, where · denotes the unpaired electron. Density functional theory (DFT) calculations (CRYSTAL06, B3PW, all-electron basis sets, and 1 × 2 × 2 supercell) support the proposed structural model and yield the following ¹¹B hyperfine coupling constants: A ₁ = −0.429 mT, A ₂ = −0.901 mT, A ₃ = −0.954 mT, in excellent agreement with the experimental results. The [BO₄]⁰ center undergoes the onset of thermal decay at ~200°C and is completely annealed out at 375°C but can be restored readily by gamma-ray irradiation. Isothermal annealing experiments show that the [BO₄]⁰ center exhibits a second-order thermal decay with an activation energy of 0.96 eV. The confirmation of the [BO₄]⁰ center (and its formation from the O₃BOH precursor) in datolite has implications for not only understanding of BOHCs in alkali borosilicate glasses but also their applications to nuclear waste disposal.     

High-pressure X-ray diffraction study of ε-FeOOH

An in situ synchrotron X-ray diffraction study was carried out on ε-FeOOH at room temperature up to a pressure of 8.6 GPa using the energy-dispersive method. The linear compressibility was determined to be β _a  = 1.69(3) × 10⁻³ GPa⁻¹, β _b  = 2.86(6) × 10⁻³ GPa⁻¹, and β _c  = 1.73(5) × 10⁻³ GPa⁻¹. The b-axis of the unit cell is more compressible than the a and c axes. The pressure–volume data were fitted to a third-order Birch–Murnaghan equation of state. The best fit was found using a room temperature isothermal bulk modulus of K ₀ = 126(3) GPa and its pressure derivative K′ = 10(1).     

EPR and magnetism of the nanostructured natural carbonaceous material shungite

The X-band EPR and magnetic susceptibility in the temperature range 4.2–300 K study of the shungite-I, natural nanostructured material from the deposit of Shunga are reported. Obtained results allow us to assign the EPR signal to conduction electrons, estimate their number, N _P, and evaluate the Pauli paramagnetism contribution to shungite susceptibility. A small occupation (~5%) of the localized nonbonding π states in the zigzag edges of the open-ended graphene-like layers and/or on σ (sp ^2+x ) orbitals in the curved parts of the shungite globules has been also revealed. The observed temperature dependence of the EPR linewidth can be explained by the earlier considered interaction of conduction π electrons with local phonon modes associated with the vibration of peripheral carbon atoms of the open zigzag-type edges and with peripheral carbon atoms cross-linking different nanostructures. The relaxation time T ₂ and diffusion time T _D are found to have comparable values (2.84 × 10⁻⁸ and 1.73 × 10⁻⁸ s at 5.2 K, respectively), and similar dependence on temperature. The magnetic measurements have revealed the suppression of orbital diamagnetism due to small amount of large enough fragments of the graphene layers.     

Thermal equations of state of dioctahedral micas on the join muscovite–paragonite

 Powder diffraction measurements at simultaneous high pressure and temperature on samples of 2M1 polytype of muscovite (Ms) and paragonite (Pg) were performed at the beamline ID30 of ESRF (Grenoble), using the Paris-Edinburgh cell. The bulk moduli of Ms, calculated from the least-squares fitting of V–P data on each isotherm using a second-order Birch–Murnaghan EoS, were: 57.0(6), 55.1(7), 51.1(7) and 48.9(5) GPa on the isotherms at 298, 573, 723 and 873 K, respectively. The value of (∂K _T /∂T) was −0.0146(2) GPa K⁻¹. The thermal expansion coefficient α varied from 35.7(3) × 10⁻⁶ K⁻¹ at P ambient to 20.1(3) × 10⁻⁶ K⁻¹ at P = 4 GPa [(∂α/∂P) _T = −3.9(1) × 10⁻⁶ GPa⁻¹ K⁻¹]. The corresponding values for Pg on the isotherms at 298, 723 and 823 K were: bulk moduli 59.9(5), 55.7(6) and 53.8(7) GPa, (∂K _T /∂T) −0.0109(1) GPa K⁻¹. The thermal expansion coefficient α varied from 44.1(2) × 10⁻⁶ K⁻¹ at P ambient to 32.5(2) × 10⁻⁶ K⁻¹ at P = 4 GPa [(∂α/∂P) _T = −2.9(1) × 10⁻⁶ GPa⁻¹ K⁻¹]. Thermoelastic coefficients showed that Pg is stiffer than Ms; Ms softens more rapidly than Pg upon heating; thermal expansion is greater and its variation with pressure is smaller in Pg than in Ms. Received: 28 January 2002 / Accepted: 5 April 2002     

X-ray diffraction study of arsenopyrite at high pressure

The high-pressure X-ray diffraction study of a natural arsenopyrite was investigated up to 28.2 GPa using in situ angle-dispersive X-ray diffraction and a diamond anvil cell at National Synchrotron Light Source, Brookhaven National Laboratory. The 16:3:1 methanol–ethanol–water mixture was used as a pressure-transmitting medium. Pressures were measured using the ruby-fluorescence method. No phase change has been observed up to 28.2 GPa. The isothermal equation of state (EOS) was determined. The values of K ₀, and K′ ₀ refined with a third-order Birch–Murnaghan EOS are K ₀ = 123(9) GPa, and K′ ₀ = 5.2(8). Furthermore, we confirm that the linear compressibilities (β) along a, b and c directions of arsenopyrite is elastically isotropic (β _a  = 6.82 × 10⁻⁴, β _b  = 6.17 × 10⁻⁴ and β _c  = 6.57 × 10⁻⁴ GPa⁻¹).