(Ca,Mg)2Si2O6 clinopyroxenes: A solution model based on nonconvergent site-disorder

Experiments at high pressure and temperature indicate that excess Ca may be dissolved in diopside. If the (Ca, Mg)₂Si₂O₆ clinopyroxene solution extends to more Ca-rich compositions than CaMgSi₂O₆, macroscopic regular solution models cannot strictly be applied to this system. A nonconvergent site-disorder model, such as that proposed by Thompson (1969, 1970), may be more appropriate. We have modified Thompson's model to include asymmetric excess parameters and have used a linear least-squares technique to fit the available experimental data for Ca-Mg orthopyroxene-clinopyroxene equilibria and Fe-free pigeonite stability to this model. The model expressions for equilibrium conditions \(\mu _{{\text{Mg}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{opx}}} = \mu _{{\text{Mg}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{cpx}}} \) (reaction A) and \(\mu _{{\text{Ca}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{opx}}} = \mu _{{\text{Ca}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{cpx}}} \) (reaction B) are given by: 1 $$\begin{gathered} \Delta \mu _{\text{A}}^{\text{O}} = {\text{RT 1n}}\left[ {\frac{{(X_{{\text{Mg}}}^{{\text{opx}}} )^2 }}{{X_{{\text{Mg}}}^{{\text{M1}}} \cdot X_{{\text{Mg}}}^{{\text{M2}}} }}} \right] - \frac{1}{2}\{ W_{21} [2(X_{{\text{Ca}}}^{{\text{M2}}} )^3 - (X_{{\text{Ca}}}^{{\text{M2}}} ] \hfill \\ {\text{ + 2W}}_{{\text{22}}} [X_{{\text{Ca}}}^{{\text{M2}}} )^2 - (X_{{\text{Ca}}}^{{\text{M2}}} )^3 + \Delta {\text{G}}_{\text{*}}^{\text{0}} (X_{{\text{Ca}}}^{{\text{M1}}} \cdot X_{{\text{Ca}}}^{{\text{M2}}} )\} \hfill \\ {\text{ + W}}^{{\text{opx}}} (X_{{\text{Wo}}}^{{\text{opx}}} )^2 \hfill \\ \Delta \mu _{\text{B}}^{\text{O}} = {\text{RT 1n}}\left[ {\frac{{(X_{{\text{Ca}}}^{{\text{opx}}} )^2 }}{{X_{{\text{Ca}}}^{{\text{M1}}} \cdot X_{{\text{Ca}}}^{{\text{M2}}} }}} \right] - \frac{1}{2}\{ 2W_{21} [2(X_{{\text{Mg}}}^{{\text{M2}}} )^2 - (X_{{\text{Mg}}}^{{\text{M2}}} )^3 ] \hfill \\ {\text{ + W}}_{{\text{22}}} [2(X_{{\text{Mg}}}^{{\text{M2}}} )^3 - (X_{{\text{Mg}}}^{{\text{M2}}} )^2 + \Delta {\text{G}}_{\text{*}}^{\text{0}} (X_{{\text{Mg}}}^{{\text{M1}}} \cdot X_{{\text{Mg}}}^{{\text{M2}}} )\} \hfill \\ {\text{ + W}}^{{\text{opx}}} (X_{{\text{En}}}^{{\text{opx}}} )^2 \hfill \\ \hfill \\ \end{gathered} $$ where 1 $$\begin{gathered} \Delta \mu _{\text{A}}^{\text{O}} = 2.953 + 0.0602{\text{P}} - 0.00179{\text{T}} \hfill \\ \Delta \mu _{\text{B}}^{\text{O}} = 24.64 + 0.958{\text{P}} - (0.0286){\text{T}} \hfill \\ {\text{W}}_{{\text{21}}} = 47.12 + 0.273{\text{P}} \hfill \\ {\text{W}}_{{\text{22}}} = 66.11 + ( - 0.249){\text{P}} \hfill \\ {\text{W}}^{{\text{opx}}} = 40 \hfill \\ \Delta {\text{G}}_*^0 = 155{\text{ (all values are in kJ/gfw)}}{\text{.}} \hfill \\ \end{gathered} $$ . Site occupancies in clinopyroxene were determined from the internal equilibrium condition 1 $$\begin{gathered} \Delta G_{\text{E}}^{\text{O}} = - {\text{RT 1n}}\left[ {\frac{{X_{{\text{Ca}}}^{{\text{M1}}} \cdot X_{{\text{Mg}}}^{{\text{M2}}} }}{{X_{{\text{Ca}}}^{{\text{M2}}} \cdot X_{{\text{Mg}}}^{{\text{M1}}} }}} \right] + \tfrac{1}{2}[(2{\text{W}}_{{\text{21}}} - {\text{W}}_{{\text{22}}} )(2{\text{X}}_{{\text{Ca}}}^{{\text{M2}}} - 1) \hfill \\ {\text{ + }}\Delta G_*^0 (X_{{\text{Ca}}}^{{\text{M1}}} - X_{{\text{Ca}}}^{{\text{M2}}} ) + \tfrac{3}{2}(2{\text{W}}_{{\text{21}}} - {\text{W}}_{{\text{22}}} ) \hfill \\ {\text{ (1}} - 2X_{{\text{Ca}}}^{{\text{M1}}} )(X_{{\text{Ca}}}^{{\text{M1}}} + \tfrac{1}{2})] \hfill \\ \end{gathered} $$ where δG _E ⁰ =153+0.023T+1.2P. The predicted concentrations of Ca on the clinopyroxene Ml site are low enough to be compatible with crystallographic studies. Temperatures calculated from the model for coexisting ortho- and clinopyroxene pairs fit the experimental data to within 10° in most cases; the worst discrepancy is 30°. Phase relations for clinopyroxene, orthopyroxene and pigeonite are successfully described by this model at temperatures up to 1,600° C and pressures from 0.001 to 40 kbar. Predicted enthalpies of solution agree well with the calorimetric measurements of Newton et al. (1979). The nonconvergent site disorder model affords good approximations to both the free energy and enthalpy of clinopyroxenes, and, therefore, the configurational entropy as well. This approach may provide an example for Febearing pyroxenes in which cation site exchange has an even more profound effect on the thermodynamic properties.     

The reversed alumina contents of orthopyroxene in equilibrium with spinel and forsterite in the system MgO-Al2O3-SiO2

Equilibrium alumina contents of orthopyroxene coexisting with spinel and forsterite in the system MgO-Al₂O₃-SiO₂ have been reversed at 15 different P-T conditions, in the range 1,030–1,600° C and 10–28 kbar. The present data and three reversals of Danckwerth and Newton (1978) have been modeled assuming an ideal pyroxene solid solution with components Mg₂Si₂O₆ (En) and MgAl₂SiO₆ (MgTs), to yield the following equilibrium condition (J, bar, K): $$\begin{gathered} RT{\text{ln(}}X_{{\text{MgTs}}} {\text{/}}X_{{\text{En}}} {\text{) + 29,190}} - {\text{13}}{\text{.42 }}T + 0.18{\text{ }}T + 0.18{\text{ }}T^{1.5} \hfill \\ + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP = 0,} \hfill \\ \end{gathered} $$ where $$\begin{gathered} + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP} \hfill \\ = [0.013 + 3.34 \times 10^{ - 5} (T - 298) - 6.6 \times 10^{ - 7} P]P. \hfill \\ \end{gathered} $$ The data of Perkins et al. (1981) for the equilibrium of orthopyroxene with pyrope have been similarly fitted with the result: $$\begin{gathered} - RT{\text{ln(}}X_{{\text{MgTs}}} \cdot X_{{\text{En}}} {\text{) + 5,510}} - 88.91{\text{ }}T + 19{\text{ }}T^{1.2} \hfill \\ + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP = 0,} \hfill \\ \end{gathered} $$ where $$\begin{gathered} + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP} \hfill \\ = [ - 0.832 - 8.78{\text{ }} \times {\text{ 10}}^{ - {\text{5}}} (T - 298) + 16.6{\text{ }} \times {\text{ 10}}^{ - 7} P]{\text{ }}P. \hfill \\ \end{gathered} $$ The new parameters are in excellent agreement with measured thermochemical data and give the following properties of the Mg-Tschermak endmember: $$H_{f,970}^0 = - 4.77{\text{ kJ/mol, }}S_{298}^0 = 129.44{\text{ J/mol}} \cdot {\text{K,}}$$ and $$V_{298,1}^0 = 58.88{\text{ cm}}^{\text{3}} .$$ The assemblage orthopyroxene+spinel+olivine can be used as a geothermometer for spinel lherzolites, subject to a choice of thermodynamic mixing models for multicomponent orthopyroxene and spinel. An ideal two-site mixing model for pyroxene and Sack's (1982) expressions for spinel activities provide, with the present experimental calibration, a geothermometer which yields temperatures of 800° C to 1,350° C for various alpine peridotites and 850° C to 1,130° C for various volcanic inclusions of upper mantle origin.     

The system Fe-Si-O: Oxygen buffer calibrations to 1,500K

The five solid-phase oxygen buffers of the system Fe-Si-O, iron-wuestite (IW), wuestite-magnetite (WM), magnetite-hematite (MH), quartz-iron-fayalite (QIF) and fayalite-magnetite-quartz (FMQ) have been recalibrated at 1 atm pressure and temperatures from 800°–1,300° C, using a thermogravimetric gas mixing furnace. The oxygen fugacity, \(f_{{\text{O}}_{\text{2}} }\) was measured with a CaO-doped ZrO₂ electrode. Measurements were made also for wuestite solid solutions in order to determine the redox behavior of wuestites with O/Fe ratios varying from 1.05 to 1.17. For FMQ, additional determinations were carried out at 1 kb over a temperature range of 600° to 800° C, using a modified Shaw membrane. Results agree reasonably well with published data and extrapolations. The reaction parameters K, ΔG _r ^o , ΔH _r ^o , and ΔS _r ^o were calculated from the following log \(f_{{\text{O}}_{\text{2}} }\) /T relations (T in K): $$\begin{gathered} {\text{IW }}\log f_{{\text{O}}_{\text{2}} } = - 26,834.7/T + 6.471\left( { \pm 0.058} \right) \hfill \\ {\text{ }}\left( {{\text{800}} - 1,260{\text{ C}}} \right), \hfill \\ {\text{WM }}\log f_{{\text{O}}_{\text{2}} } = - 36,951.3/T + 16.092\left( { \pm 0.045} \right) \hfill \\ {\text{ }}\left( {{\text{1,000}} - 1,300{\text{ C}}} \right), \hfill \\ {\text{MH }}\log f_{{\text{O}}_{\text{2}} } = - 23,847.6/T + 13.480\left( { \pm 0.055} \right) \hfill \\ {\text{ }}\left( {{\text{1,040}} - 1,270{\text{ C}}} \right), \hfill \\ {\text{QIF }}\log f_{{\text{O}}_{\text{2}} } = - 27,517.5/T + 6.396\left( { \pm 0.049} \right) \hfill \\ {\text{ }}\left( {{\text{960}} - 1,140{\text{ C}}} \right), \hfill \\ {\text{FMQ }}\log f_{{\text{O}}_{\text{2}} } = - 24,441.9/T + 8.290\left( { \pm 0.167} \right) \hfill \\ {\text{ }}\left( {{\text{600}} - 1,140{\text{ C}}} \right). \hfill \\ \end{gathered}$$ These experimentally determined reaction parameters were combined with published 298 K data to determine the parameters G_f, H_f, and S_f for the phases wuestite, magnetite, hematite, and fayalite from 298 K to the temperatures of the experiments. The T? \(f_{{\text{O}}_{\text{2}} }\) data for wuestite solid solutions were used to obtain activities, excess free energies and Margules mixing parameters. The new data provide a more reliable, consistent and complete reference set for the interpretation of redox reactions at elevated temperatures in experiments and field settings encompassing the crust, mantle and core as well as extraterrestrial environments.     

Thermodynamic properties of andradite and application to skarn with coexisting andradite and hedenbergite

The enthalpy of formation of andradite (Ca₃Fe₂Si₃O₁₂) has been estimated as-5,769.700 (±5) kJ/mol from a consideration of the calorimetric data on entropy (316.4 J/mol K) and of the experimental phaseequilibrium data on the reactions: 1 $$\begin{gathered} 9/2 CaFeSi_2 O_6 + O_2 = 3/2 Ca_3 Fe_2 Si_3 O_{12} + 1/2 Fe_3 O_4 + 9/2 SiO_2 (a) \hfill \\ Hedenbergite andradite magnetite quartz \hfill \\ \end{gathered} $$ 1 $$\begin{gathered} 4 CaFeSi_2 O_6 + 2 CaSiO_3 + O_2 = 2 Ca_3 Fe_2 Si_3 O_{12} + 4 SiO_2 (b) \hfill \\ Hedenbergite wollastonite andradite quartz \hfill \\ \end{gathered} $$ 1 $$\begin{gathered} 18 CaSiO_3 + 4 Fe_3 O_4 + O_2 = 6Ca_3 Fe_2 Si_3 O_{12} (c) \hfill \\ Wollastonite magnetite andradite \hfill \\ \end{gathered} $$ 1 $$\begin{gathered} Ca_3 Fe_2 Si_3 O_{12} = 3 CaSiO_3 + Fe_2 O_3 . (d) \hfill \\ Andradite pseudowollastonite hematite \hfill \\ \end{gathered} $$ and $$log f_{O_2 } = E + A + B/T + D(P - 1)/T + C log f_{O_2 } .$$ Oxygen-barometric scales are presented as follows: $$\begin{gathered} E = 12.51; D = 0.078; \hfill \\ A = 3 log X_{Ad} - 4.5 log X_{Hd} ; C = 0; \hfill \\ B = - 27,576 - 1,007(1 - X_{Ad} )^2 - 1,476(1 - X_{Hd} )^2 . \hfill \\ \end{gathered} $$ For the assemblage andradite (Ad)-hedenbergite (Hd)-magnetite-quartz: $$\begin{gathered} E = 13.98; D = 0.0081; \hfill \\ A = 4 log(X_{Ad} / X_{Hd} ); C = 0; \hfill \\ B = - 29,161 - 1,342.8(1 - X_{Ad} )^2 - 1,312(1 - X_{Hd} )^2 . \hfill \\ \end{gathered} $$ For the assemblage andradite-hedenbergite-wollastonite-quartz: 1 $$\begin{gathered} E = 13.98;{\text{ }}D = 0.0081; \hfill \\ A = 4\log (X_{Ad} /X_{Hd} );{\text{ C = 0;}} \hfill \\ B = - 29,161 - 1,342.8(1 - X_{Ad} )^2 - 1,312(1 - X_{Hd} )^2 . \hfill \\ \end{gathered} $$ For the assemblage andradite-hedenbergite-calcitequartz: 1 $$\begin{gathered} E = - 1.69;{\text{ }}D = - 0.199; \hfill \\ A = 4\log (X_{Ad} /X_{Hd} );{\text{ C = 2;}} \hfill \\ B = - 20,441 - 1,342.8(1 - X_{Ad} )^2 - 1,312(1 - X_{Hd} )^2 . \hfill \\ \end{gathered} $$ For the assemblage andradite-hedenbergite-wollastonite-calcite: 1 $$\begin{gathered} E = - 17.36;{\text{ }}D = - 0.403; \hfill \\ A = 4\log (X_{Ad} /X_{Hd} );{\text{ C = 4;}} \hfill \\ B = - 11,720 - 1,342.8(1 - X_{Ad} )^2 - 1,312(1 - X_{Hd} )^2 \hfill \\ \end{gathered} $$ The oxygen fugacity of formation of those skarns where andradite and hedenbergite assemblage is typical can be calculated by using the above equations. The oxygen fugacity of formation of this kind of skarn ranges between carbon dioxide/graphite and hematite/magnetite buffers. It increases from the inside zones to the outside zones, and appears to decrease with the ore-types in the order Cu, Pb?Zn, Fe, Mo, W(Sn) ore deposits.     

Thermal expansion and spontaneous strain associated with the normal-incommensurate phase transition in melilites

The linear thermal expansions of åkermanite (Ca₂MgSi₂O₇) and hardystonite (Ca₂ZnSi₂O₇) have been measured across the normal-incommensurate phase transition for both materials. Least-squares fitting of the high temperature (normal phase) data yields expressions linear in T for the coefficients of instantaneous linear thermal expansion, $$\alpha _1 = \frac{1}{l}\frac{{dl}}{{dT}}$$ for åkermanite: $$\begin{gathered} \alpha _{[100]} = 6.901(2) \times 10^{ - 6} + 1.834(2) \times 10^{ - 8} T \hfill \\ \alpha _{[100]} = - 2.856(1) \times 10^{ - 6} + 11.280(1) \times 10^{ - 8} T \hfill \\ \end{gathered} $$ for hardystonite: $$\begin{gathered} \alpha _{[100]} = 15.562(5) \times 10^{ - 6} - 1.478(3) \times 10^{ - 8} T \hfill \\ \alpha _{[100]} = - 11.115(5) \times 10^{ - 6} + 11.326(3) \times 10^{ - 8} T \hfill \\ \end{gathered} $$ Although there is considerable strain for temperatures within 10° C of the phase transition, suggestive of a high-order phase transition, there appears to be a finite ΔV of transition, and the phase transition is classed as “weakly first order”.     

The concept of geochemical potential field exemplified by ore formation field

Geochemical potential field is defined as the scope within the earth’s space where a given component in a certain phase of a certain material system is acted upon by a diffusion force, depending on its spatial coordinatesX, Y andZ. The three coordinates follow the relations: $$NF_{ix} = - \frac{{\partial \mu }}{{\partial x}}, NF_{iy} = - \frac{{\partial \mu }}{{\partial y}}, NF_{iz} = - \frac{{\partial \mu }}{{\partial z}}$$ The characteristics of such a field can be summarized as: (1) The summation of geochemical potentials related to the coordinatesX, Y, Z, or pseudo-velocity head, pseudo-pressure head and pseudo-potential head of a certain component in the earth is a constant as given by $$\mu _x + \mu _y + \mu _z = c$$ or $$\mu _{x2} + \mu _{y2} + \mu _{z2} = \mu _{x1} + \mu _{y1} + \mu _{z1} $$ Derived from these relations is the principle of geochemical potential conservation. The following relations have the same physical significance: $$\mu _k + \mu _u + \mu _p = c$$ or $$\mu _{k2} + \mu _{u2} + \mu _{p2} = \mu _{k1} + \mu _{u1} + \mu _{p1} $$ (2) Geochemical potential field is a vector field quantified by geochemical field intensity which is defined as the diffusion force applied to one molecular volume (or one atomic volume) of a certain component moving from its higher concentration phase to lower concentration phase. The geochemical potential field intensity is given by $$\begin{gathered} E = - grad\mu \hfill \\ E = \frac{{RT}}{x}i + \frac{{RT}}{y}j + \frac{{RT}}{z}K \hfill \\ \end{gathered} $$ The present theory has been inferred to interpret the mechanism of formation of some tungsten ore deposits in China.     

Assemblages among tephroite,pyroxmangite, rhodochrosite,quartz: Experimental data and occurrences in the Rhetic Alps

Reactions involving the phases quartz-rhodochrosite-tephroite-pyroxmangite-fluid have been studied experimentally in the system MnO-SiO₂-CO₂-H₂O at a pressure of 2 000 bars and resulted in the following expressions 1 $$\begin{gathered} {\text{Rhodochrosite + Quartz = Pyroxmangite + CO}}_2 \hfill \\ {\text{ log}}_{{\text{10}}} K^{{\text{2000 bars}}} = - \frac{{11.765}}{T} + 18.618. \hfill \\ {\text{Rhodochrosite + Pyroxmangite = Tephroite + CO}}_2 \hfill \\ {\text{ log}}_{{\text{10}}} K^{{\text{2000 bars}}} = - \frac{{7.083}}{T} + 11.870. \hfill \\ \end{gathered}$$ which can be used to derive data for the remaining two reactions among the phases under consideration. Field data from the Alps are in agreement with the metamorphic sequence resulting from the experiments.     

Intrinsic oxygen Fugacity measurements on chromites from the Bushveld Complex and their petrogenetic significance

Oxygen Fugacity measurements were carried out on chromites from the Eastern Bushveld Complex (Maandagshoek) and are compared with former measurements on chromites from the western Bushveld Complex (Zwartkop Chrome Mine). These results together with those of Hill and Roeder (1974) yield the following conditions of formation for the massive chromitite layers: Western Bushveld Complex (Zwartkop Chrome Mine) $$\begin{gathered} Layer{\text{ }}T(^\circ C) p_{O_2 } (atm) \hfill \\ LG3{\text{ 1160}} - {\text{1234 10}}^{ - {\text{5}}} - 10^{ - 7.6} \hfill \\ LG4{\text{ 1175}} - {\text{1200 10}}^{ - 6.35} - 10^{ - 7.20} \hfill \\ LG6{\text{ 1162}} - {\text{1207 10}}^{ - 6.20} - 10^{ - 7.50} \hfill \\ \hfill \\ \end{gathered} $$ Eastern Bushveld Complex (Farm Maandagshoek) $$\begin{gathered} {\text{LXI 1115}} - {\text{1150 10}}^{ - 7.80} - 10^{ - 8.80} \hfill \\ ( = {\text{Steelpoort Seam)}} \hfill \\ {\text{LX 1125 10}}^{ - 8.25} \hfill \\ {\text{V 1120 10}}^{ - 8.55} \hfill \\ {\text{LII 1120 10}}^{ - 8.0} - 10^{ - 8.60} \hfill \\ \end{gathered} $$ The comparison of the data shows, that the chronitite layers within each particular sequence were formed under approximately identicalp _{o 2}- andT-conditions. The chromites from the western Bushveld Complex, however, were formed at higher temperatures and higher oxygen fugacities than the chromites from the eastern Bushveld Complex. Fromp _{o 2}-T-curves of disseminated chromites and the temperatures derived above, the following conditions of formation for the host rocks were obtained: Western Bushveld Complex $$T = 1200^\circ {\text{C; }}p_{{\text{o}}_{\text{2}} } = 10^{ - 7.25} - 10^{ - 7.50} $$ Eastern Bushveld Complex $$T = 1125^\circ {\text{C; }}p_{{\text{o}}_{\text{2}} } = 10^{ - 8.50} - 10^{ - 9.25} $$ Consequently, the host rocks in the Zwartkop-Chrome-Mine, were formed under higher temperatures and higher oxygen fugacities than the host rocks at Maandagshoek. The rock sequence in the Zwartkop-Chrome-Mine therefore originated in an earlier stage of the differentiation of the Bushveld magma. Comparison of the chromites from the host rocks with the chromites from massive layers supports Ulmer's (1969) thesis that an increase of the oxygen fugacity is responsible for the formation of massive chromitite layers. The values in this investigation show that increases of only about 0.5–1.0 log units are necessary to enhance chromitite layer formation.     

Non-ideal Mg-Fe binary mixing in cordierite: Constraints from experimental data on Mg-Fe partitioning in garnet and cordierite and a reformulation of garnet-cordierite geothermometer

The non-ideal regular Mg-Fe binary in cordierite has been derived through multivariate linear regression of the expressionRT InK_D +(P- 1)ΔVK ₁ ⁰ , 298 along with updated subfegular mixing parameter of almandine-pyrope solution (Hackler and Wood 1989; Berman 1990). The data base used for multivariate analyses consists of published experimental data (n = 177) on Mg-Fe partitioning between garnet and cordierite in theP-T range 650–1050°C and 4–12 K bar. The non-ideality can be approximated by temperature-dependent Margules parameters. The retrieved values of ΔH_<T> ^o and ΔH_<T> ^o of exchange reaction between garnet and cordierite and enthalpy and entropy of mixing of Mg-Fe cordierite were combined with recent quaternary (Fe-Mg-Ca-Mn) mixing data in garnet to obtain the geothermometric expressions to determine temperature (T Kelvin): $$\begin{gathered} T(WH) = 6832 + 0.031(P - 1) - \{ 166(X_{Mg}^{Gt} )^2 - 506(X_{Fe}^{Gt} )^2 + 680X_{Fe}^{Gt} X_{Mg}^{Gt} + 336(X_{Ca} + X_{Mn} ) \hfill \\ (X_{Mg} - X_{Fe} )^{Gt} - 3300X_{Ca}^{Gt} - 358X_{Mn}^{Gt} \} + 954(X_{Fe} - X_{Mg} )^{Crd} /1.987\ln K_D + 3.41 + 1.5X_{Ca}^{Gt} \hfill \\ + 1.23(X_{Fe} - X_{Mg} )^{Crd} \hfill \\ \end{gathered} $$ $$\begin{gathered} T(Br) = 6920 + 0.031(p - 1) - \{ 18(X_{Mg}^{Gt} )^2 - 296(X_{Fe}^{Gt} )^2 + 556X_{Fe}^{Gt} X_{Mg}^{Gt} - 6339X_{Ca}^{Gt} X_{Mg}^{Gt} \hfill \\ - 99(X_{Ca}^{Gt} )^2 + 4687X_{Ca}^{Gt} (X_{Mg} - X_{Fe}^{Gt} ) - 4269X_{Ca}^{Gt} X_{Fe}^{Gt} - 358X_{Mn}^{Gt} \} + 640(X_{Fe} - X_{Mg} )^{Crd} \hfill \\ + 1.90X_{Ca}^{Gt} (X_{Mg} - X_{Ca} )^{Gt} . \hfill \\ \end{gathered} $$      

Energetics of hydration of cordierite and water barometry in cordierite-granulites

The Gibbs free energy and volume changes attendant upon hydration of cordierites in the system magnesian cordierite-water have been extracted from the published high pressure experimental data at \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) =P _total, assuming an ideal one site model for H₂O in cordierite. Incorporating the dependence of ΔG and ΔV on temperature, which was found to be linear within the experimental conditions of 500°–1,000°C and 1–10,000 bars, the relation between the water content of cordierite and P, T and \(f_{{\text{H}}_{\text{2}} {\text{O}}} \) has been formulated as $$\begin{gathered} X_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{crd}}} = \hfill \\ \frac{{f_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{P, T}}} }}{{\left[ {{\text{exp}}\frac{1}{{RT}}\left\{ {64,775 - 32.26T + G_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{1, }}T} - P\left( {9 \times 10^{ - 4} T - 0.5142} \right)} \right\}} \right] + f_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{P, T}}} }} \hfill \\ \end{gathered} $$ The equation can be used to compute H₂O in cordierites at \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) <1. Our results at different P, T and partial pressure of water, assuming ideal mixing of H₂O and CO₂ in the vapour phase, are in very good agreement with the experimental data of Johannes and Schreyer (1977, 1981). Applying the formulation to determine \(X_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{crd}}} \) in the garnet-cordierite-sillimanite-plagioclase-quartz granulites of Finnish Lapland as a test case, good agreement with the gravimetrically determined water contents of cordierite was obtained. Pressure estimates, from a thermodynamic modelling of the Fe-cordierite — almandine — sillimanite — quartz equilibrium at \(P_{{\text{H}}_{\text{2}} {\text{O}}} = 0\) and \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) =P_total, for assemblages from South India, Scottish Caledonides, Daly Bay and Hara Lake areas are compatible with those derived from the garnetplagioclase-sillimanite-quartz geobarometer.     

Experimental study of Fe-Mg distribution between biotite and orthopyroxene at P=490 MPa

Partitioning of Mg and Fe between coexisting biotite and orthopyroxene has been experimentally determined at temperatures 700, 750 and 800° C and 490 MPa total pressure in the system KAlO₂-MgO-FeO-SiO₂-H₂O. Oxygen fugacity was controlled by the QFM buffer. Starting materials were synthetic minerals of differing Fe/(Fe+Mg) values. Run products were analyzed for partitioning of components by a microprobe. Orthopyroxene was established to be notably inhomogeneous, whereas biotite was essentially homogeneous. To establish equilibrium relations, statistical treatment of the results of each experiment in addition to the whole complex of experimental data was applied. The regression equations for isotherms of the Fe-Mg partitioning between the minerals studied have been obtained. As a result, the equation for a two-dimensional regression may be written as: $$\begin{gathered} Y = (A + A_1 t + A_2 t^2 )(X - X^4 ) + (B + B_1 t + B_1 t^2 )(X^2 - X^4 ) + \hfill \\ (C + C_1 t + C_1 t^2 )(X^3 - X^4 ) + X^4 {\text{ where }}Y = X_{{\text{Opx}}}^{{\text{Fe}}} ;{\text{ X}} = {\text{X}}_{{\text{Bi}}}^{{\text{Fe}}} ; \hfill \\ t = 1000/T,K, \hfill \\ \begin{array}{*{20}c} {A = {\text{ }}4.59398,} & {A_1 = - {\text{ }}8.29838,} & {A_2 = {\text{ }}4.97316,} \\ {B = - 11.13731,} & {B_1 = {\text{ }}28.19304,} & {B_2 = - 20.98240,} \\ {A = {\text{ }}8.25072,} & {C_1 = - 20.80485,} & {C_2 = {\text{ }}15.35967} \\ \end{array} \hfill \\ {\text{ }}\sigma = 0.0143{\text{ }} \hfill \\ \end{gathered}$$ . This equation enables extrapolation of partitioning isotherms over a wide range of temperatures.     

Non-ideal mixing in the phlogopite-annite binary: constraints from experimental data on Mg−Fe partitioning and a reformulation of the biotite-garnet geothermometer

The existing experimental data [Ferry and Spear 1978; Perchuk and Lavrent'eva 1983] on Mg?Fe partitioning between garnet and biotite are disparate. The underlying assumption of ideal Mg?Fe exchange between the minerals has been examined on the basis of recently available thermochemical data. Using the updated mixing parameters for the pyrope-almandine asymmetric regular solution as inputs [Ganguly and Saxena 1984; Hackler and Wood 1984], thermodynamic analysis points to non-ideal mixing in the phlogopite-annite binary in the temperature range of 550°C–950°C. The non-ideality can be approximated by a temperature-independent, one constant Margules parameter. The retrieved values for enthalpy of mixing for Mg?Fe biotites and the standard state enthalpy and entropy changes of the exchange reaction were combined with existing thermochemical data on grossular-pyrope and grossular-almandine binaries to obtain geothermometric expressions for Mg?Fe fractionation between biotite and garnet. [T in K] $$\begin{gathered} {\text{T(HW) = [20286 + 0}}{\text{.0193P - \{ 2080(X}}_{{\text{Mg}}}^{{\text{Gt}}} {\text{)}}^{\text{2}} {\text{ - 6350(X}}_{{\text{Fe}}}^{{\text{Gt}}} {\text{)}}^{\text{2}} \hfill \\ {\text{ - 13807(X}}_{{\text{Ca}}}^{{\text{Gt}}} {\text{)(1 - X}}_{{\text{Mn}}}^{{\text{Gt}}} {\text{) + 8540(X}}_{{\text{Fe}}}^{{\text{Gt}}} {\text{)(X}}_{{\text{Mg}}}^{{\text{Gt}}} {\text{)(1 - X}}_{{\text{Mn}}}^{{\text{Gt}}} {\text{)}} \hfill \\ {\text{ + 4215(X}}_{{\text{Ca}}}^{{\text{Gt}}} {\text{)(X}}_{{\text{Mg}}}^{{\text{Gt}}} {\text{ - X}}_{{\text{Fe}}}^{{\text{Gt}}} {\text{)\} + 4441}}{{{\text{(2X}}_{{\text{Mg}}}^{{\text{Bt}}} {\text{ - 1)]}}} \mathord{\left/ {\vphantom {{{\text{(2X}}_{{\text{Mg}}}^{{\text{Bt}}} {\text{ - 1)]}}} {{\text{[13}}{\text{.138}}}}} \right. \kern-\nulldelimiterspace} {{\text{[13}}{\text{.138}}}} \hfill \\ {\text{ + 8}}{\text{.3143 InK}}_{\text{D}} {\text{ + 6}}{\text{.276(X}}_{{\text{Ca}}}^{{\text{Gt}}} ){\text{(1 - X}}_{{\text{Mn}}}^{{\text{Gt}}} )] \hfill \\ {\text{T(GS) = [13538 + 0}}{\text{.0193P - \{ 837(X}}_{{\text{Mg}}}^{{\text{Gt}}} )^{\text{2}} {\text{ - 10460(X}}_{{\text{Fe}}}^{{\text{Gt}}} )^2 \hfill \\ {\text{ - 13807(X}}_{{\text{Ca}}}^{{\text{Gt}}} )(1{\text{ - X}}_{{\text{Mn}}}^{{\text{Gt}}} {\text{) + 19246(X}}_{{\text{Fe}}}^{{\text{Gt}}} ){\text{(X}}_{{\text{Mg}}}^{{\text{Gt}}} ){\text{(1 - X}}_{{\text{Mn}}}^{{\text{Gt}}} ) \hfill \\ {\text{ }}{{{\text{ + 5649(X}}_{{\text{Ca}}}^{{\text{Gt}}} ){\text{(X}}_{{\text{Mg}}}^{{\text{Gt}}} {\text{ - X}}_{{\text{Fe}}}^{{\text{Gt}}} ){\text{\} + 7972(2X}}_{{\text{Mg}}}^{{\text{Bt}}} {\text{ - 1)]}}} \mathord{\left/ {\vphantom {{{\text{ + 5649(X}}_{{\text{Ca}}}^{{\text{Gt}}} ){\text{(X}}_{{\text{Mg}}}^{{\text{Gt}}} {\text{ - X}}_{{\text{Fe}}}^{{\text{Gt}}} ){\text{\} + 7972(2X}}_{{\text{Mg}}}^{{\text{Bt}}} {\text{ - 1)]}}} {{\text{[6}}{\text{.778}}}}} \right. \kern-\nulldelimiterspace} {{\text{[6}}{\text{.778}}}} \hfill \\ {\text{ + 8}}{\text{.3143InK}}_{\text{D}} {\text{ + 6}}{\text{.276(X}}_{{\text{Ca}}}^{{\text{Gt}}} )(1{\text{ - X}}_{{\text{Mn}}}^{{\text{Gt}}} )] \hfill \\ \end{gathered} $$ The reformulated geothermometer is an improvement over existing biotite-garnet geothermometers because it reconciles the experimental data sets on Fe?Mg partitioning between the two phases and is based on updated activity-composition relationship in Fe?Mg?Ca garnet solid solutions.     

The reaction garnet+clinopyroxene+quartz =2 orthopyroxene+anorthite: A potential geobarometer for granulites

A mineralogic geobarometer based on the reaction garnet+clinopyroxene+quartz=2 orthopyroxene+anorthite is proposed. The geobarometric formulations for the Fe- and Mg- end member equilibria are $$\begin{gathered} P_{({\text{Fe}})} {\text{ }}({\text{bars}}){\text{ = 32}}{\text{.097 }}T{\text{ }} - {\text{ 26385 }} - {\text{ 22}}{\text{.79 (}}T - 848 - T1{\text{n(}}T/848{\text{))}} \hfill \\ {\text{ }} - (3.655 + 0.0138T){\text{ }}\left( {\frac{{{\text{(}}T - 848{\text{)}}^{\text{2}} }}{T}} \right) \hfill \\ {\text{ }} - {\text{(3}}{\text{.123) }}T1{\text{n }}\frac{{(a_{a{\text{n}}}^{{\text{Plag}}} )(a_{{\text{fs}}}^{{\text{P}}\ddot u{\text{x}}} )^2 }}{{(a_{{\text{alm}}}^{{\text{Gt}}} )(a_{{\text{hed}}}^{{\text{Opx}}} )}} \hfill \\ P_{({\text{Mg}})} {\text{ (bars) = 9}}{\text{.270 }}T + 4006 - 0.9305{\text{ }}(T - 848 - T1{\text{n (}}T/848{\text{)}}) \hfill \\ {\text{ }} - (1.1963{\text{ }} - {\text{ }}6.0128{\text{ x 10}}^{ - {\text{3}}} T)\left( {\frac{{(T - 848)^2 }}{T}} \right) \hfill \\ {\text{ }} - 3.489{\text{ }}T1{\text{n }}\frac{{(a_{an}^{{\text{Plag}}} ){\text{ }}(a_{{\text{ens}}}^{{\text{Opx}}} )}}{{{\text{(}}a_{{\text{pyr}}}^{{\text{Gt}}} {\text{) (}}a_{{\text{diop}}}^{{\text{Cpx}}} {\text{)}}}}. \hfill \\ \end{gathered}$$ The end member thermodynamic data have been taken from the data base of Helgeson et al. (1978) and Saxena and Erikson (1983). The activities of pyroxene components and anorthite in plagioclase have been modelled after Wood and Banno (1973) and Newton (1983) respectively. The activities of pyrope and almandine are calculated from the binary interaction parameters for garnet solid solutions proposed by Saxena and Erikson (1983). Pressures computed from these equations for fifty sets of published mineral data from several granulite areas are comparable with those obtained from dependable geobarometers. The pressure values determined from the Fe-end member equilibrium appear to be more reasonable than those from the Mg-end member reaction. It is likely that the difference in pressures computed from the Fe- and Mg-end members, ΔP ^*, have been caused by non-ideal mixing in the phases, especially in garnets.     

Stochastic analysis of bed-thickness distribution of sediments

On formation of a bed and distribution of bed thickness, A. N. Kolmogorov presented a mathematical explanation that if repetitive alternations of material accumulation and erosion form a sequence of beds, the resultant bed-thickness distribution curve takes a shape truncated by the ordinate at zero thickness. In this truncated distribution curve, its continuation and extension from positive to negative thickness represents the distribution of beds with negative thickness, that is, the depth of erosion. When a distribution curve, including both positive and negative parts, is expressed by a function f(x),the ratio \(\int_0^\infty {f(x)dx to} \int_{ - \infty }^\infty {f(x)dx} \) ,called Kolmogorov's coefficient and designated as p,is a parameter representing the degree of accumulation in the depositional environment. On the assumption that f(x)is described by the Gaussian distribution function, the coefficient pfor Permian and Pliocene sequences in central Japan was calculated. The coefficients also were obtained from published data for different types of sediments from other areas. It was determined that they are more or less different depending on their depositional environments. The calculated results are summarized as follows: $$\begin{gathered} p = 0.80 - 1.0for{\text{ }}alluvial{\text{ }}or{\text{ }}fluvial{\text{ }}deposits \hfill \\ p = 0.65 - 0.95for{\text{ }}nearshore{\text{ }}sediments \hfill \\ p = 0.55 - 0.95for{\text{ }}geosynclinal{\text{ }}sediments \hfill \\ p = 0.90 - 1.0for{\text{ }}varves \hfill \\ \end{gathered} $$ In addition, a ratio \(q = \int_0^\infty {xf(x)dx/} \int_{ - \infty }^\infty {|x|f(x)dx} \) ,called Kolmogorov's ratio in this paper, is introduced for estimating a degree of total thickness actually observed in the field relative to total thickness once present in a basin. The calculated results of Kolmogorov's ratio are as follows: $$\begin{gathered} q = 0.88 - 1.0for{\text{ }}alluvial{\text{ }}or{\text{ }}fluvial{\text{ }}deposits \hfill \\ q = 0.68 - 0.98for{\text{ }}nearshore{\text{ }}sediments \hfill \\ q = 0.55 - 0.96for{\text{ }}geosynclinal{\text{ }}sediments \hfill \\ q = 0.92 - 1.0for{\text{ }}varves \hfill \\ \end{gathered} $$ The sedimentological significance of these values is discussed.     

High temperature calorimetry of sulfide systems

The standard enthalpies of formation of FeS (troilite), FeS₂ (pyrite), Co_0.9342S, Co₃S₄ (linnaeite), Co₉S₈ (cobalt pentlandite), CoS₂ (cattierite), CuS (covellite), and Cu₂S (chalcocite) have been determined by high temperature direct reaction calorimetry at temperatures between 700 K and 1021 K. The following results are reported: $$\Delta {\rm H}_{f,FeS}^{tr} = - 102.59 \pm 0.20kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,FeS}^{py} = - 171.64 \pm 0.93kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,Co_{0.934} S} = - 99.42 \pm 1.52kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,Co_9 S_8 }^{ptl} = - 885.66 \pm 16.83kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,Co_3 S_4 }^{In} = - 347.47 \pm 7.27kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,CoS_2 }^{ct} = - 150.94 \pm 4.85kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,Cu_2 S}^{cc} = - 80.21 \pm 1.51kJ mol^{ - 1} ,$$ and $$\Delta {\rm H}_{f,CuS}^{cv} = - 53.14 \pm 2.28kJ mol^{ - 1} ,$$ The enthalpy of formation of CuFeS₂ (chalcopyrite) from (CuS+FeS) and from (Cu+FeS₂) was determined by solution calorimetry in a liquid Ni_0.60S_0.40 melt at 1100 K. The results of these measurements were combined with the standard enthalpies of formation of CuS, FeS, and FeS₂, to calculate the standard enthalpy of formation of CuFeS₂. We found \(\Delta {\rm H}_{f,CuFeS_2 }^{ccp} = - 194.93 \pm 4.84kJ mol^{ - 1}\) . Our results are compared with earlier data given in the literature; generally the agreement is good and our values agree with previous estimates within the uncertainties present in both.     

Experimental study of subsolidus phase relations and mixing properties of pyroxene and plagioclase in the system Na2O-CaO-Al2O3-SiO2

In the system Na₂O-CaO-Al₂O₃-SiO₂ (NCAS), the equilibrium compositions of pyroxene coexisting with grossular and corundum were experimentally determined at 40 different P-T conditions (1,100–1,400° C and 20.5–38 kbar). Mixing properties of the Ca-Tschermak — Jadeite pyroxene inferred from the data are (J, K): $$\begin{gathered} G_{Px}^{xs} = X_{{\text{CaTs}}} X_{{\text{Jd}}} [14,810 - 7.15T - 5,070(X_{{\text{CaTs}}} - X_{{\text{Jd}}} ) \hfill \\ {\text{ }} - 3,350(X_{{\text{CaTs}}} - X_{{\text{Jd}}} )^2 ] \hfill \\ \end{gathered} $$ The excess entropy is consistent with a complete disorder of cations in the M2 and the T site. Compositions of coexisting pyroxene and plagioclase were obtained in 11 experiments at 1,190–1,300° C/25 kbar. The data were used to infer an entropy difference between low and high anorthite at 1,200° C, corresponding to the enthalpy difference of 9.6 kJ/mol associated with the C \(\bar 1\) =I \(\bar 1\) transition in anorthite as given by Carpenter and McConnell (1984). The resulting entropy difference of 5.0 J/ mol · K places the transition at 1,647° C. Plagioclase is modeled as ideal solutions, C \(\bar 1\) and I \(\bar 1\) , with a non-first order transition between them approximated by an empirical expression (J, bar, K): $$\Delta G_T = \Delta G_{1,473} \left[ {1 - 3X_{Ab} \tfrac{{T^4 - 1,473^4 }}{{\left( {1,920 - 0.004P} \right)^4 - 1,473^4 }}} \right],$$ where $$\Delta G_{1,473} = 9,600 - 5.0T - 0.02P$$ The derived mixing properties of the pyroxene and plagioclase solutions, combined with the thermodynamic properties of other phases, were used to calculate phase relations in the NCAS system. Equilibria involving pyroxene+plagioclase +grossular+corundum and pyroxene+plagioclase +grossular+kyani te are suitable for thermobarometry. Albite is the most stable plagioclase.     

Thermodynamic properties of glaucophane new data from calorimetric and spectroscopic measurements

New data concerning glaucophane are presented. New high temperature drop calorimetry data from 400 to 800 K are used to constrain the heat capacity at high temperature. Unpublished low temperature calorimetric data are used to estimate entropy up to 900 K. These data, corrected for composition, are fitted for C _p and S to the polynomial expressions (J · mol^?1 · K^?2) for T> 298.15 K: $$\begin{gathered} C_p = 11.4209 * 10^2 - 40.3212 * 10^2 /T^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} - 41.00068 * 10^6 /T^2 \hfill \\ + 52.1113 * 10^8 /T^3 \hfill \\ \end{gathered} $$ $$\begin{gathered} S = 539 + 11.4209 * 10^2 * \left( {\ln T - \ln 298.15} \right) - 80.6424 * 10^2 \hfill \\ * \left( {T^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} - 1/\left( {298.15} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} } \right) + 20.50034 * 10^6 \hfill \\ * \left( {T^{ - 2} - 1/\left( {298.15} \right)^2 } \right) - 17.3704 * 10^8 * \left( {T^{ - 3} - \left( {1/298.15} \right)^3 } \right) \hfill \\ \end{gathered} $$ IR and Raman spectra from 50 to 3600 cm^?1 obtained on glaucophane crystals close to the end member composition are also presented. These spectroscopic data are used with other data (thermal expansion, acoustic velocities etc.) in vibrational modelling. This last method provides an independent way for the determination of the thermodynamic properties (C_p and entropy). The agreement between measured and calculated properties is excellent (less than 2% difference between 100 and 1000 K). It is therefore expected that vibrational modelling could be applied to other amphiboles for which spectroscopic data are available. Finally, the enthalpy of formation of glaucophane is calculated.     

Gold Solubility in SiO2—HCl—H2O System at 200℃：a Preliminary Assessment of the Implications of Silicification with Regard to Gold Mineralization

The complexation between gold and silica was experimentally, confirmed and calibrated at 200 °C: $$\begin{gathered} Au^ + + H_3 SiO_4^ - \rightleftharpoons AuH_3 SiO_4^0 \hfill \\ \log K_{(200^\circ C)} = 19.26 \pm 0.4 \hfill \\ \end{gathered} $$ Thermodynamic calculations show that AuH₃SiO ₄ ⁰ would be far more abundant than AuCl ₂ ^? under physicochemical conditions of geological interest, suggesting that silica is much more important than chloride as ligands for gold transport. In systems containing both sulfur and silica, AuH₃SiO ₄ ⁰ would be increasingly more important than Au (HS) ₂ ^? as the proportion of SiO₂ in the system increases. The dissolution of gold in aqueous SiO₂ solutions can be described by the reaction: $$\begin{gathered} Au + 1/4O_2 + H_4 SiO_4^0 \rightleftharpoons AuH_3 SiO_4^0 + 1/2H_2 O \hfill \\ log K_{(200^\circ C)} = 6.23 \hfill \\ \end{gathered} $$ which indicates that SiO₂ precipitation is an effective mechanism governing gold deposition, and thus explains the close association of silicification and gold mineralization.     

Dislocation-assisted diffusion of oxygen in albite

Oxygen diffusion in albite has been determined by the integrating (bulk ¹⁸O) method between 750° and 450° C, for a P _H2O of 2 kb. The original material has a low dislocation density (<10⁶ cm^?2), and its lattice diffusion coefficient (D ₁), given below, agrees well with previous determinations. A sample was deformed at high temperature and pressure to produce a uniform dislocation density of 5 × 10⁹ cm^?2. The diffusion coefficient (D _a) for this deformed material, given below, is about 0.5 and 0.7 orders of magnitude larger than D ₁ at 700° and 450° C, respectively. This enhancement is believed due to faster diffusion along the cores of dislocations. Assuming a dislocation core radius of 4 Å, the calculated pipe diffusion coefficient (D _p), given below, is about 5 orders of magnitude larger than D ₁. These results suggest that volume diffusion at metamorphic conditions may be only slightly enhanced by the presence of dislocations. $$\begin{gathered} D_1 = 9.8 \pm 6.9 \times 10^{ - 6} (cm^2 /\sec ) \hfill \\ {\text{ }} \cdot \exp [ - 33.4 \pm 0.6(kcal/mole)/RT] \hfill \\ \end{gathered} $$ $$\begin{gathered} D_a = 7.6 \pm 4.0 \times 10^{ - 6} (cm^2 /\sec ) \hfill \\ {\text{ }} \cdot \exp [ - 30.9 \pm 1.1(kcal/mole)/RT] \hfill \\ \end{gathered} $$ $$\begin{gathered} D_p \approx 1.2 \times 10^{ - 1} (cm^2 /\sec ) \hfill \\ {\text{ }} \cdot \exp [ - 29.8(kcal/mole)/RT]. \hfill \\ \end{gathered} $$      

Magnesiochloritoid: Compressibility and high pressure structure refinement

The response of magnesiochloritoid to pressure has been studied by single crystal X-ray diffraction in a diamond anvil cell, using crystals with composition Mg_1.3Fe_0.7Al₄Si₂O₁₀(OH)₄. The unit cell parameters decrease from a = 9.434 (3), b = 5.452 (2), c = 18.136 (5) Å, β = 101.42° (2) (1 bar pressure) to a = 9.370 (7), b = 5.419 (5), c = 17.88 (1) Å, β = 101.5° (1) (42 kbar pressure), following a slightly anisotropic compression pattern (linear compressibilities parallel to unit cell edges: β _a = 1.85, β _b = 1.74, β_c = 3.05 × 10^?4 kbar^?1) with a bulk modulus of 1480 kbar. Perpendicular to c, the most compressible direction, the crystal structure (space group C2/c) consists of two kinds of alternating octahedral layers connected via isolated SiO₄ tetrahedra. With increasing pressure the slightly wavy layer [Mg_1.3Fe_0.7AlO₂(OH)₄] tends to flatten. Furthermore, the octahedra in this layer, with all cations underbonded, are more compressible than the octahedra in the (A1₃O₈) layer with slightly overbonded aluminum. Comparison between high-pressure and high-temperature data yields the following equations: $$\begin{gathered} a_{P,T} = 9.434{\text{ }}{\AA} - 174 \cdot 10^{ - 5} {\text{ }}{\AA}{\text{kb}}^{{\text{ - 1}}} \cdot P \hfill \\ {\text{ }} + 9 \cdot 10^{ - 5} {\text{ }}{\AA}^\circ C^{ - 1} \cdot (T - 25^\circ C) \hfill \\ b_{P,T} = 5.452{\text{ }}{\AA} - 95 \cdot 10^{ - 5} {\text{ }}{\AA}{\text{kb}}^{{\text{ - 1}}} \cdot P \hfill \\ {\text{ }} + 5 \cdot 65 \cdot 10^{ - 5} {\text{ }}{\AA}^\circ C^{ - 1} \cdot (T - 25^\circ C) \hfill \\ c_{P,T} = 18.136{\text{ }}{\AA} - 549 \cdot 10^{ - 5} {\text{ }}{\AA}{\text{kb}}^{{\text{ - 1}}} \cdot P \hfill \\ {\text{ }} + 16 \cdot 2^{ - 5} {\text{ }}{\AA}^\circ C^{ - 1} \cdot (T - 25^\circ C) \hfill \\ \end{gathered} $$ with P in kbar and T in °C. These equations indicate that the unit cell and bond geometry of magnesiochloritoid at formation conditions do not differ greatly from those at the outcrop conditions, e.g. the calculated unitcell volume is 917.3 ų at P = 16 kbar and T=500 °C, whereas the observed volume at room conditions is 914.4 ų. In addition, they show that the specific gravity increases from formation at depth to outcrop at surface conditions.