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.     

Heat capacity and thermodynamic functions of xenotime YPO<Subscript>4</Subscript>(c) at 0–1600 K

The heat capacity of xenotime YPO₄(c) was measured by adiabatic calorimetry at 4.78–348.07 K. Our experimental and literature data on H ⁰(T)-H ⁰(298.15 K) of Y orthophosphate were utilized to derive the C _p ⁰(T) function of xenotime at 0–1600 K, which was then used to calculate the values of thermodynamic functions: entropy, enthalpy change, and reduced Gibbs energy. These functions assume the following values at 298.15 K: C _p ⁰ (298.15 K) = 99.27 ± 0.02 J K⁻¹ mol⁻¹, S ⁰(298.15 K) = 93.86 ± 0.08 J K⁻¹ mol⁻¹, H ⁰(298.15 K) − H ⁰(0) = 15.944 ± 0.005 kJ mol⁻¹, Φ⁰(298.15 K) = 40.38 ± 0.08 J K⁻¹ mol⁻¹. The value of the free energy of formation Δ _f G ⁰(YPO₄, 298.15 K) is −1867.9 ± 1.7 kJ mol⁻¹.     

Low-temperature thermodynamic properties of disordered zeolites of the natrolite group

 The heat capacity of paranatrolite and tetranatrolite with a disordered distribution of Al and Si atoms has been measured in the temperature range of 6–309 K using the adiabatic calorimetry technique. The composition of the samples is represented with the formula (Na_1.90K_0.22Ca_0.06)[Al_2.24Si_2.76O₁₀]·nH₂O, where n=3.10 for paranatrolite and n=2.31 for tetranatrolite. For both zeolites, thermodynamic functions (vibrational entropy, enthalpy, and free energy function) have been calculated. At T=298.15 K, the values of the heat capacity and entropy are 425.1 ± 0.8 and 419.1 ±0.8 J K⁻¹ mol⁻¹ for paranatrolite and 381.0 ± 0.7 and 383.2 ± 0.7 J K⁻¹ mol⁻¹ for tetranatrolite. Thermodynamic functions for tetranatrolite and paranatrolite with compositions corrected for the amount of extraframework cations and water molecules have also been calculated. The calculation for tetranatrolite with two water molecules and two extraframework cations per formula yields: C _p (298.15)=359.1 J K⁻¹ mol⁻¹, S(298.15) −S(0)=362.8 J K⁻¹ mol⁻¹. Comparing these values with the literature data for the (Al,Si)-ordered natrolite, we can conclude that the order in tetrahedral atoms does not affect the heat capacity. The analysis of derivatives dC/dT for natrolite, paranatrolite, and tetranatrolite has indicated that the water- cations subsystem within the highly hydrated zeolite may become unstable at temperatures above 200 K. Received: 30 July 2001 / Accepted: 15 November 2001     

The reaction MgCr2O4 + SiO2 = Cr2O3 + MgSiO3 and the free energy of formation of magnesiochromite (MgCr2O4)

Low-temperature heat capacity measurements for MgCr₂O₄ have only been performed down to 52 K, and the commonly quoted third-law entropy at 298 K (106 J K⁻¹ mol⁻¹) was obtained by empirical extrapolation of these measurements to 0 K without considering the magnetic or electronic ordering contributions to the entropy. Subsequent magnetic measurements at low temperature reveal that the Néel temperature, at which magnetic ordering of the Cr³⁺ ions in MgCr₂O₄ occurs, is at ∼15 K. Hence a substantial contribution to the entropy of MgCr₂O₄ has been missed. We have determined the position of the near-univariant reaction MgCr₂O₄+SiO₂=MgSiO₃+Cr₂O₃. The reaction, which has a small positive slope in P-T space, has been bracketed at 100 K intervals between 1273 and 1773 K by reversal experiments. The reaction is extremely sluggish, and lengthy run times with a flux (H₂O, BaO-B₂O₃ or K₂O-B₂O₃) are needed to produce tight reversal brackets. The results, combined with assessed thermodynamic data for Cr₂O₃, MgSiO₃ and SiO₂, give the entropy and enthalpy of formation of MgCr₂O₄ spinel. As expected, our experimental results are not in good agreement with the presently available thermodynamic data. We obtain Δ _f H ^∘ ₂₉₈=−1759.2±1.5 kJ mol⁻¹ and S ^∘ ₂₉₈=122.1±1.0 J K⁻¹ mol⁻¹ for MgCr₂O₄. This entropy is some 16 J K⁻¹ mol⁻¹ more than the calorimetrically determined value, and implies a value for the magnetic entropy of MgCr₂O₄ consistent with an effective spin quantum number (S') for Cr³⁺ of 1/2 rather than the theoretical 3/2, indicating, as in other spinels, spin quenching. Received: 9 May 1997 / Accepted: 28 July 1997     

Heat capacity and thermodynamic functions of pretulite ScPO<Subscript>4</Subscript>(c) at 0–1600 K

The heat capacity of synthetic pretulite ScPO₄(c) was measured by adiabatic calorimetry within a temperature range of 12.13–345.31 K, and the temperature dependence of the pretulite heat capacity at 0–1600 K was derived from experimental and literature data on H ⁰(T)-H ⁰(298.15 K) for Sc orthophosphate. This dependence was used to calculate the values of the following thermodynamic functions: entropy, enthalpy change, and reduced Gibbs energy. They have the following values at 298.15 K: C _p ⁰ (298.15 K) = 97.45 ± 0.06 J K⁻¹ mol⁻¹, S ⁰(298.15 K) = 84.82 ± 0.18 J K⁻¹ mol⁻¹, H ⁰(298.15 K)-H ⁰(0) = 14.934 ± 0.016 kJ mol⁻¹, and Φ ⁰(298.15 K) = 34.73 ± 0.19 J K⁻¹mol⁻¹. The enthalpy of formation Δ _f H ⁰(ScPO₄, 298.15 K) = − 1893.6 ± 8.4 kJ mol⁻¹.     

Thermodynamic functions of eskolaite Cr<Subscript>2</Subscript>O<Subscript>3</Subscript>(c) at 0–1800 K

The heat capacity of eskolaite Cr₂O₃(c) was determined by adiabatic vacuum calorimetry at 11.99–355.83 K and by differential calorimetry at 320–480 K. Experimental data of the authors and data compiled from the literature were applied to calculate the heat capacity, entropy, and the enthalpy change of Cr₂O₃ within the temperature range of 0–1800 K. These functions have the following values at 298.15 K: C _p ⁰ (298.15) = 121.5 ± 0.2 J K⁻¹mol⁻¹, S ⁰(298.15) = 80.95 ± 0.14 J K⁻¹mol⁻¹, and H ⁰(298.15)-H ⁰(0) = 15.30±0.02 kJ mol⁻¹. Data were obtained on the transitions from the antiferromagnetic to paramagnetic states at 228–457 K; it was determined that this transition has the following parameters: Neel temperature T _N = 307 K, Δ _tr S = 6.11 ± 0.12 J K⁻¹mol⁻¹ and δ _tr H = 1.87 ± 0.04 kJ mol⁻¹.     

Heat capacity,entropy and phase equilibria of stishovite

The low-temperature heat capacity (C _P) of stishovite (SiO₂) synthesized with a multi-anvil device was measured over the range of 5–303 K using the heat capacity option of a physical properties measurement system (PPMS) and around ambient temperature using a differential scanning calorimeter (DSC). The entropy of stishovite at standard temperature and pressure calculated from DSC-corrected PPMS data is 24.94 J mol⁻¹ K⁻¹, which is considerably smaller (by 2.86 J mol⁻¹ K⁻¹) than that determined from adiabatic calorimetry (Holm et al. in Geochimica et Cosmochimica Acta 31:2289–2307, 1967) and about 4% larger than the recently reported value (Akaogi et al. in Am Mineral 96:1325–1330, 2011). The coesite–stishovite phase transition boundary calculated using the newly determined entropy value of stishovite agrees reasonably well with the previous experimental results by Zhang et al. (Phys Chem Miner 23:1–10, 1996). The calculated phase boundary of kyanite decomposition reaction is most comparable with the experimental study by Irifune et al. (Earth Planet Sci Lett 77:245–256, 1995) at low temperatures around 1,400 K, and the calculated slope in this temperature range is mostly consistent with that determined by in situ X-ray diffraction experiments (Ono et al. in Am Mineral 92:1624–1629, 2007).     

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.     

Heat capacity,entropy, and phase equilibria of dmitryivanovite

The heat capacity (C _p ) of dmitryivanovite synthesized with a cubic press was measured in the temperature range of 5–664 K using the heat capacity option of a physical properties measurement system and a differential scanning calorimeter. The entropy of dmitryivanovite at standard temperature and pressure (STP) was calculated to be 110.1 ± 1.6 J mol⁻¹ K⁻¹ from the measured C _p data. With the help of new phase equilibrium experiments done at 1.5 GPa, the phase transition boundary between krotite and dmitryivanovite was best represented by the equation: P (GPa) = −2.1825 + 0.0025 T (K). From the temperature intercept of this phase boundary and other available thermodynamic data for krotite and dmitryivanovite, the enthalpy of formation and Gibbs free energy of formation of dmitryivanovite at STP were calculated to be −2326.7 ± 2.1 and −2,208.1 ± 2.1 kJ mol⁻¹, respectively. It is also inferred that dmitryivanovite is the stable CaAl₂O₄ phase at STP and has a wide stability field at high pressures whereas the stability field of krotite is located at high temperatures and relatively low pressures. This conclusion is consistent with natural occurrences (in Ca–Al-rich inclusions) of dmitryivanovite and krotite, where the former is interpreted as the shock metamorphic product of originally present krotite.     

Thermodynamic data of the high-pressure phase Mg5Al5Si6O21(OH)7 (Mg-sursassite)

 Calorimetric and P–V–T data for the high-pressure phase Mg₅Al₅Si₆O₂₁(OH)₇ (Mg-sursassite) have been obtained. The enthalpy of drop solution of three different samples was measured by high-temperature oxide melt calorimetry in two laboratories (UC Davis, California, and Ruhr University Bochum, Germany) using lead borate (2PbO·B₂O₃) at T=700 ^∘C as solvent. The resulting values were used to calculate the enthalpy of formation from different thermodynamic datasets; they range from −221.1 to −259.4 kJ mol⁻¹ (formation from the oxides) respectively −13892.2 to −13927.9 kJ mol⁻¹ (formation from the elements). The heat capacity of Mg₅Al₅Si₆O₂₁(OH)₇ has been measured from T=50 ^∘C to T=500 ^∘C by differential scanning calorimetry in step-scanning mode. A Berman and Brown (1985)-type four-term equation represents the heat capacity over the entire temperature range to within the experimental uncertainty: C _P (Mg-sursassite) =(1571.104 −10560.89×T ^−0.5−26217890.0 ×T ⁻²+1798861000.0×T ⁻³) J K⁻¹ mol⁻¹ (T in K). The P V T behaviour of Mg-sursassite has been determined under high pressures and high temperatures up to 8 GPa and 800 ^∘C using a MAX 80 cubic anvil high-pressure apparatus. The samples were mixed with Vaseline to ensure hydrostatic pressure-transmitting conditions, NaCl served as an internal standard for pressure calibration. By fitting a Birch-Murnaghan EOS to the data, the bulk modulus was determined as 116.0±1.3 GPa, (K ^′=4), V _T,0 =446.49 ³ exp[∫(0.33±0.05) × 10⁻⁴ + (0.65±0.85)×10⁻⁸ T dT], (K _T/T) _P  = −0.011± 0.004 GPa K⁻¹. The thermodynamic data obtained for Mg-sursassite are consistent with phase equilibrium data reported recently (Fockenberg 1998); the best agreement was obtained with Δ_f H ⁰ ₂₉₈ (Mg-sursassite) = −13901.33 kJ mol⁻¹, and S ⁰ ₂₉₈ (Mg-sursassite) = 614.61 J K⁻¹ mol⁻¹. Received: 21 September 2000 / Accepted: 26 February 2001     

Structure of magnetite (Fe<Subscript>3</Subscript>O<Subscript>4</Subscript>) above the Curie temperature: a cation ordering study

A pristine magnetite (Fe₃O₄) specimen was studied by means of Neutron Powder Diffraction in the 273–1,073 K temperature range, in order to characterize its structural and magnetic behavior at high temperatures. An accurate analysis of the collected data allowed the understanding of the behavior of the main structural and magnetic features of magnetite as a function of temperature. The magnetic moments of both tetrahedral and octahedral sites were extracted by means of magnetic diffraction up to the Curie temperature (between 773 and 873 K). A change in the thermal expansion coefficient around the Curie temperature together with an increase in the oxygen coordinate value above 700 K can be observed, both features being the result of a change in the thermal expansion of the tetrahedral site. This anomaly is not related to the magnetic transition but can be explained with an intervened cation reordering, as magnetite gradually transforms from a disordered configuration into a partially ordered one. Based on a simple model which takes into account the cation-oxygen bond length, the degree of order as a function of temperature and consequently the enthalpy and entropy of the reordering process were determined. The refined values are ΔH⁰ = −23.2(1.7) kJ mol⁻¹ and ΔS⁰ = −16(2) J K⁻¹ mol⁻¹. These results are in perfect agreement with values reported in literature (Mack et al. in Solid State Ion 135(1–4):625–630, 2000; Wu and Mason in J Am Ceramic Soc 64(9):520–522, 1981).     

Heat capacity and thermodynamic properties of GdPO<Subscript>4</Subscript> in the temperature range 0–1600 K

The heat capacity of gadolinium orthophosphate (GdPO₄) measured in the temperature range 11.15–344.11 K by adiabatic calorimetry and available literature data were used to calculate its thermodynamic functions at 0–1600 K. At 298.15 K, these functions are as follows: C _p ⁰(298.15 K) = 101.85 ± 0.05 J K⁻¹ mol⁻¹, S ⁰(298.15 K) = 123.82 ± 0.18 J K⁻¹ mol⁻¹, H ⁰(298.15 K)–H ⁰(0) = 17.250 ± 0.012 kJ mol⁻¹, and Φ ⁰(298.15 K) = 65.97 ± 0.18 J K⁻¹ mol⁻¹ The calculated Gibbs free energy of formation from the elements of GdPO₄ is Δ _f G ⁰ (298.15 K) = −1844.3 ± 4.7 kJ mol⁻¹.     

Low-temperature heat capacity of magnesioferrite (MgFe2O4)

The low-temperature heat capacity of magnesioferrite (MgFe₂O₄) was measured between 1.5 K and 300 K, and thermochemical functions were derived from the results. No heat capacity anomaly was observed. From our data, we suggest a standard entropy (298.15 K) for magnesioferrite of 120.8±0.6 J mol⁻¹ K⁻¹, which is about 2.4 J mol⁻¹ K⁻¹ higher than previously reported calorimetric studies; but is in rough agreement with predictions from sets of internally consistent thermodynamic data.     

Thermodynamic properties of nickel chromite (NiCr2O4) based on adiabatic calorimetry at low temperatures

 The hitherto unknown low-temperature heat capacity of nickel chromite (NiCr₂O₄) was measured between 8 and 381 K using adiabatic calorimetry, and some thermochemical functions [C_P(T), S(T), S°₂₉₈, H_(T)−H₍₀₎] were derived from the results. The standard entropy (S°₂₉₈=140.0 ± 0.3 J mol⁻¹ K⁻¹) for nickel chromite was calculated from the results. Our calorimetric measurements indicate three major anomalies in the heat-capacity curve at temperatures between 8 and 381 K. A short literature review indicates that two of these anomalies can be accounted for, whereas an anomaly peaking at 29 K has not been reported previously.     

Thermodynamic mixing behavior of synthetic Ca-Tschermak–diopside pyroxene solid solutions: II. Heat of mixing and activity–composition relationships

The heat of mixing for the binary solid solution diopside–Ca-Tschermak was investigated at T = 980 K by lead borate solution calorimetry. A new statistical technique was applied to overcome the problem of using experimental data of various precisions. A two-parameter Margules model was fitted to the calorimetric data leading to W _CaTs–Di^H = 31.3 ± 3.4 kJ mol⁻¹ and W ^H _Di–CaTs = 2.4 ± 4.3 kJ mol⁻¹. The results are in good agreement with calorimetric data given in the literature. They agree also with enthalpy data that were extracted from phase equilibrium experiments. With configurational entropy values taken from the literature, the volume and the vibrational entropy, presented in Part I of this work, and the enthalpy data of this study, the activity–composition relationships of the diopside–Ca-Tschermak binary were calculated.     

Phase equilibrium and calorimetric study of the transition of MnTiO3 from the ilmenite to the lithium niobate structure and implications for the stability field of perovskite

The phase boundary between MnTiO₃ I (ilmenite structure) and MnTiO₃ II (lithium niobate structure) has been determined by analysis of quench products from reversal experiments in a cubic anvil apparatus at 1073–1673 K and 43–75 kbar using mixtures of MnTiO₃ I and II as starting materials. Tight brackets of the boundary give P(kbar)=121.2−0.045 T(K). Thermodynamic analysis of this boundary gives ΔH^o=5300±1000 J·mol⁻¹, ΔS^o = 1.98 ±1J·K⁻¹· mol⁻¹. The enthalpy of transformation obtained directly by transposed-temperature-drop calorimetry is 8359 ±2575 J·mol⁻¹. Possible topologies of the phase relations among the ilmenite, lithium niobate, and perovskite polymorphs are constrained using the above data and the observed (reversible with hysteresis) transformation of II to III at 298 K and 20–30 kbar (Ross et al. 1989). The observed II–III transition is likely to lie on a metastable extension of the II–III boundary into the ilmenite field. However the reversed I–II boundary, with its negative dP/ dT does represent stable equilibrium between ilmenite and lithium niobate, as opposed to the lithium niobate being a quench product of perovskite. We suggest a topology in which the perovskite occurs stably at low T and high P with a triple point (I, II, III) at or below 1073 K near 70 kbar. The I–II boundary would have a negative P-T slope while the II–III and I–III boundaries would be positive, implying that entropy decreases in the order lithium niobate, ilmenite, perovskite. The inferred positive slope of the ilmenite-perovskite transition in MnTiO₃ is different from the negative slopes in silicates and germanates. These thermochemical parameters are discussed in terms of crystal structure and lattice vibrations.     

In situ high-temperature powder X-ray diffraction study on the spinel solid solutions (Mg1?xMnx)Cr2O4

Using a conventional high-T furnace, the solid solutions between magnesiochromite and manganochromite, (Mg_1−x Mn _x )Cr₂O₄ with x = 0.00, 0.19, 0.44, 0.61, 0.77 and 1.00, were synthesized at 1,473 K for 48 h in open air. The ambient powder X-ray diffraction data suggest that the V–x relationship of the spinels does not show significant deviation from the Vegard’s law. In situ high-T powder X-ray diffraction measurements were taken up to 1,273 K at ambient pressure. For the investigated temperature range, the unit-cell parameters of the spinels increase smoothly with temperature increment, indicating no sign of cation redistribution between the tetrahedral and octahedral sites. The V–T data were fitted with a polynomial expression for the volumetric thermal expansion coefficient (a_T = a₀ + a₁ T + a₂ T ^{- 2} \alpha_{T} = a_{0} + a_{1} T + a_{2} T^{ - 2} ), which yielded insignificant a ₂ values. The effect of the composition on a ₀ is adequately described by the equation a ₀ = [17.7(8) − 2.4(1) × x] 10⁻⁶ K⁻¹, whereas that on a ₁ by the equation a ₁ = [8.6(9) + 2.1(11) × x] 10⁻⁹ K⁻².     

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.     

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.     

P–V–T equation of state of Ca3Al2Si3O12 grossular garnet

The thermoelastic parameters of synthetic Ca₃Al₂Si₃O₁₂ grossular garnet were examined in situ at high-pressure and high-temperature by energy dispersive X-ray diffraction, using a Kawai-type multi-anvil press apparatus coupled with synchrotron radiation. Measurements have been conducted at pressures up to 20 GPa and temperatures up to 1,650 K: this P, T range covered the entire high-P, T stability field of grossular garnet. The analysis of room temperature data yielded V _0,300 = 1,664 ± 2 ?³ and K ₀ = 166 ± 3 GPa for K^￠₀ K^{\prime}_{0} fixed to 4.0. Fitting of our P–V–T data by means of the high-temperature third order Birch–Murnaghan or the Mie–Grüneisen–Debye thermal equations of state, gives the thermoelastic parameters: (∂K _0,T /∂T) _P  = −0.019 ± 0.001 GPa K⁻¹ and α _0,T  = 2.62 ± 0.23 × 10⁻⁵ K⁻¹, or γ ₀ = 1.21 for fixed values q ₀ = 1.0 and θ ₀ = 823 (Isaak et al. Phys Chem Min19:106–120, 1992). From the comparison of fits from two different approaches, we propose to constrain the bulk modulus of grossular garnet and its pressure derivative to K _T0 = 166 GPa and K^￠_T0 K^{\prime}_{T0}  = 4.03–4.35. Present results are compared with previously determined thermoelastic properties of grossular-rich garnets.