Thermodynamics of plagioclase II: Temperature evolution of the spontaneous strain at the I\bar 1 \leftrightarrow P\bar 1 phase transition in anorthite

The temperature dependence of the lattice parameters of pure anorthite with high Al/Si order reveals the predicted tricritical behaviour of the \(I\bar 1 \leftrightarrow P\bar 1\) phase transition at T _c ^* =510 K. The spontaneous strain couples to the order parameter Q° as x _i∝S _xQ ⁱ Q°² with S _xQ ¹ =4.166×10^?3, S _xQ ² =0.771×10^?3, S _xQ ³ =?7.223×10^?3 for the diagonal elements. The temperature dependence of Q° is $$Q^{\text{o}} = \left( {1 - \frac{T}{{510}}} \right)^\beta ,{\text{ }}\beta = \tfrac{{\text{1}}}{{\text{4}}}$$ A strong dependence of T _c ^* , S _xQ ⁱ and β is predicted for Al/Si disordered anorthite.     

The P–V–T equation of state of CaPtO3 post-perovskite

Orthorhombic post-perovskite CaPtO₃ is isostructural with post-perovskite MgSiO₃, a deep-Earth phase stable only above 100 GPa. Energy-dispersive X-ray diffraction data (to 9.4 GPa and 1,024 K) for CaPtO₃ have been combined with published isothermal and isobaric measurements to determine its P–V–T equation of state (EoS). A third-order Birch–Murnaghan EoS was used, with the volumetric thermal expansion coefficient (at atmospheric pressure) represented by α(T) = α₀ + α₁(T). The fitted parameters had values: isothermal incompressibility, $ K_{{T_{0} }} $  = 168.4(3) GPa; $ K_{{T_{0} }}^{\prime } $  = 4.48(3) (both at 298 K); $ \partial K_{{T_{0} }} /\partial T $  = ?0.032(3) GPa K^?1; α₀ = 2.32(2) × 10^?5 K^?1; α₁ = 5.7(4) × 10^?9 K^?2. The volumetric isothermal Anderson–Grüneisen parameter, δ _T , is 7.6(7) at 298 K. $ \partial K_{{T_{0} }} /\partial T $ for CaPtO₃ is similar to that recently reported for CaIrO₃, differing significantly from values found at high pressure for MgSiO₃ post-perovskite (?0.0085(11) to ?0.024 GPa K^?1). We also report axial P–V–T EoS of similar form, the first for any post-perovskite. Fitted to the cubes of the axes, these gave $ \partial K_{{aT_{0} }} /\partial T $  = ?0.038(4) GPa K^?1; $ \partial K_{{bT_{0} }} /\partial T $  = ?0.021(2) GPa K^?1; $ \partial K_{{cT_{0} }} /\partial T $  = ?0.026(5) GPa K^?1, with δ _T  = 8.9(9), 7.4(7) and 4.6(9) for a, b and c, respectively. Although $ K_{{T_{0} }} $ is lowest for the b-axis, its incompressibility is the least temperature dependent.     

Carbon Dioxide in igneous petrogenesis: I

A number of experimental CO₂ solubility data for silicate and aluminosilicate melts at a variety of P- T conditions are consistent with solution of CO₂ in the melt by polymer condensation reactions such as SiO _4(m ^4? +CO_2(v)+Si _n O _3n+1(m) ⁽²ⁿ⁺¹⁾ ?Si _n+1O _3n+4(m) ^(2n+4)? +CO _3(m ^)2? . For various metalsilicate systems the relative solubility of CO₂ should depend markedly on the relative Gibbs free change of reaction. Experimental solubility data for the systems Li₂O-SiO₂, Na₂O-SiO₂, K₂O-SiO₂, CaO-SiO₂, MgO-SiO₂ and other aluminosilicate melts are in complete accord with predictions based on Gibbs Free energies of model polycondesation reactions. A rigorous thermodynamic treatment of published P- T-wt.% CO₂ solubility data for a number of mineral and natural melts suggests that for the reaction CO_2(m) ? CO_2(v)

1. CO₂-melt mixing may be considered ideal (i.e., { \(a_{{\text{CO}}_{\text{2}} }^m = X_{{\text{CO}}_{\text{2}} }^m \) );
2. \(\bar V_{{\text{CO}}_{\text{2}} }^m \) , the partial molal volume of CO₂ in the melt, is approximately equal to 30 cm³ mole^?1 and independent of P and T;
3. Δ C _p ⁰ is approximately equal to zero in the T range 1,400° to 1,650 °C and
4. enthalpies and entropies of the dissolution reaction depend on the ratio of network modifiers to network builders in the melt. Analytic expressions which relate the CO₂ content of a melt to P, T, and \(f_{{\text{CO}}_{\text{2}} } \) for andesite, tholeiite and olivine melilite melts of the form

$$\ln X_{{\text{CO}}_{\text{2}} }^m = \ln f_{{\text{CO}}_{\text{2}} } - \frac{A}{T} - B - \frac{C}{T}(P - 1)$$ have been determined. Regression parameters are (A, B, C): andesite (3.419, 11.164, 0.408), tholeiite (14.040, 5.440,0.393), melilite (9.226, 7.860, 0.352). The solubility equations are believed to be accurate in the range 3<P<30 kbar and 1,100°<T<1,650 °C. A series of CO₂ isopleth diagrams for a wide range of T and P are drawn for andesitic, tholeiitic and alkalic melts.     

Vibrational mode analysis and heat capacity calculation of K2SiSi3O9-wadeite

The phonon dispersions and vibrational density of state (VDoS) of the K₂SiSi₃O₉-wadeite (Wd) have been calculated by the first-principles method using density functional perturbation theory. The vibrational frequencies at the Brillouin zone center are in good correspondence with the Raman and infrared experimental data. The calculated VDoS was then used in conjunction with a quasi-harmonic approximation to compute the isobaric heat capacity (C _P ) and vibrational entropy ( $S_{298}^{0}$ ), yielding C _P (T) = 469.4(6) ? 2.90(2) × 10³  T ^?0.5 ? 9.5(2) × 10⁶  T ^?2 + 1.36(3) × 10⁹  T ^?3 for the T range of 298–1,000 K and $S_{298}^{0}$  = 250.4 J mol^?1 K^?1. In comparison, these thermodynamic properties were calculated by a second method, the classic Kieffer’s lattice vibrational model. On the basis of the vibrational mode analysis facilitated by the first-principles simulation result, we developed a new Kieffer’s model for the Wd phase. This new Kieffer’s model yielded C _P (T) = 475.9(6) ? 3.15(2) × 10³  T ^?0.5 – 8.8(2) × 10⁶  T ^?2 + 1.31(3) × 10⁹  T ^?3 for the T range of 298–1,000 K and $S_{298}^{0}$  = 249.5(40) J mol^?1 K^?1, which are in good agreement both with the results from our first method containing the component of the first-principles calculation and with some calorimetric measurements in the literature.     

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

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

Ephesite,Na(LiAl2) [Al2Si2O10] (OH)2: I. thermal stability and standard state thermodynamic properties

Ephesite, Na(LiAl₂) [Al₂Si₂O₁₀] (OH)₂, has been synthesized for the first time by hydrothermal treatment of a gel of requisite composition at 300≦T(° C)≦700 and \(P_{H_2 O}\) upto 35 kbar. At \(P_{H_2 O}\) between 7 and 35 kbar and above 500° C, only the 2M₁ polytype is obtained. At lower temperatures and pressures, the 1M polytype crystallizes first, which then inverts to the 2M₁ polytype with increasing run duration. The X-ray diffraction patterns of the 1M and 2M₁ poly types can be indexed unambiguously on the basis of the space groups C2 and Cc, respectively. At its upper thermal stability limit, 2M₁ ephesite decomposes according to the reaction (1) $$\begin{gathered} {\text{Na(LiAl}}_{\text{2}} {\text{) [Al}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{{\text{10}}} {\text{] (OH)}}_{\text{2}} \hfill \\ {\text{ephesite}} \hfill \\ {\text{ = Na[AlSiO}}_{\text{4}} {\text{] + LiAl[SiO}}_{\text{4}} {\text{] + }}\alpha {\text{ - Al}}_{\text{2}} {\text{O}}_{\text{3}} {\text{ + H}}_{\text{2}} {\text{O}} \hfill \\ {\text{nepheline }}\alpha {\text{ - eucryptite corundum}} \hfill \\ \end{gathered}$$ Five reversal brackets for (1) have been established experimentally in the temperature range 590–750° C, at \(P_{H_2 O}\) between 400 and 2500 bars. The equilibrium constant, K, for this reaction may be expressed as (2) $$log K{\text{ = }}log f_{{\text{H}}_{\text{2}} O}^* = 7.5217 - 4388/T + 0.0234 (P - 1)T$$ where \(f_{H_2 O}^* = f_{H_2 O} (P,T)/f_{H_2 O}^0\) (1,T), with T given in degrees K, and P in bars. Combining these experimental data with known thermodynamic properties of the decomposition products in (1), the following standard state (1 bar, 298.15 K) thermodynamic data for ephesite were calculated: H _f,298.15 ⁰ =-6237372 J/mol, S _298.15 ⁰ =300.455 J/K·mol, G _298.15 ⁰ =-5851994 J/mol, and V _298.15 ⁰ =13.1468 J/bar·mol.     

Bimetasomatic zoning in the CaO-MgO-SiO2-H2O-CO2 system: Experiments with the use of natural rock samples

Experiments reproducing the development of bimetasomatic zoning in the CaO-MgO-SiO₂-H₂O-CO₂ system were conducted at elevated P-T parameters with the use of samples of naturally occurring quartzdolomite and calcite-serpentinite rocks. In order to maintain mass transfer exclusively via the diffusion-controlled mechanism, we used the method of the ensured compaction of the cylindrical sample surface with a thin-walled gold tube. In the course of the experiments, a single diopside zone ～2.5 × 10^?5 m thick was obtained at the quartz-dolomite interface at T = 600°C, $P_{H_2 O + CO_2 } $ = 200 MPa, and $X_{CO_2 } $ = 0.5 for 25–40 days and a succession of metasomatic zones at T = 750°C, $P_{H_2 O + CO_2 } $ = 300 MPa, and $X_{CO_2 } $ = 0.4 for 48 days. The metasomatic zones were as follows (listed in order from quartz to dolomite): wollastonite ‖ diopside ‖ tremolite ‖ calcite + forsterite; with the average width of the diopside zone equal to ～1.3 × 10^?5 m and the analogous part of the wollastonite zone equal to ～2.6 × 10^?5 m. Two zones (listed in order from calcite to serpentine) diopside and diopside-forsterite (the average widths of these zones were ～6 × 10^?4 and ～8 × 10^?4 m, respectively) were determined to develop at contact between serpentine and calcite during experiments that lasted 124 days at T = 500°C, $P_{H_2 O + CO_2 } $ = 200 MPa, and $X_{CO_2 } $ = 0.2–0.4. In the former and latter situations, the growth rate of the zoning ranged between 3.1 × 10^?12 and 1.2 × 10^?11 m/s and between 5.6 × 10^?11 and 7.5 × 10^?11 m/s, respectively. The higher growth rate in the latter case can be explained by the higher water mole fraction in the fluid, with this water released during serpentinite decomposition in the experiments. The development of the only diopside zone in the experiments modeling the interaction of quartz and dolomite at T = 600–650°C and $P_{H_2 O + CO_2 } $ = 200 MPa is in conflict with theoretical considerations underlain by the Korzhinskii-Fisher-Joesten model. The interaction of quartz and dolomite in the CaO-MgO-SiO₂-CO₂-H₂O system at the P-T- $X_{CO_2 } $ parameters specified above should be attended by the origin of a number of reaction zones consisting of various proportions of talc, forsterite, tremolite, diopside, and calcite. The saturation of the fluid with respect to these minerals was likely not reached, and this resulted in the degeneration of the respective stability fields in the succession of zones. Conceivably, this was related to the insufficient rates of quartz and dolomite dissolution and the relatively low diffusion rates of the dissolved species in the low-permeable medium. In the experiments with interacting calcite and serpentine, the zoning calcite ‖ diopside ‖ diopside + forsterite ‖ serpentine developed in its complete form, in agreement with the theory. Equilibrium was likely achieved in these experiments due to the higher diffusion coefficients.     

Estimating the Density and Compressibility of Seawater to High Temperatures Using the Pitzer Equations

The density and compressibility of seawater solutions from 0 to 95 °C have been examined using the Pitzer equations. The apparent molal volumes (X = V) and compressibilities (X = κ) are in the form $$ X_{\phi } = \bar{X}^{0} + A_{X} I/(1.2 \, m)\ln (1 + 1.2 \, I^{0.5} ) + \, 2{\text{RT }}m \, (\beta^{(0)X} + \beta^{(1)X} g(y) + C^{X} m) $$ where $ \bar{X}^{0} $ is the partial molal volume or compressibility, I is the ionic strength, m is the molality of sea salt, A_X is the Debye–Hückel slope for volume (X = V) or adiabatic compressibility (X = κ _s), and g(y) = (2/y ²)[1 ? (1 + y) exp(?y)] where y = 2I ^0.5. The values of the partial molal volume and compressibility ( $ \bar{X}^{0} $ ) and Pitzer parameters (β ^(0)X , β ^(1)X and C ^X ) are functions of temperature in the form $$ Y^{X} = \sum_{i} a_{i} (T-T_{\text{R}} )^{i} $$ where a _i are adjustable parameters, T is the absolute temperature in Kelvin, and T _R = 298.15 K is the reference temperature. The standard errors of the seawater fits for the specific volumes and adiabatic compressibilities are 5.35E?06 cm³ g^?1 and 1.0E?09 bar^?1, respectively. These equations can be combined with similar equations for the osmotic coefficient, enthalpy and heat capacity to define the thermodynamic properties of sea salt to high temperatures at one atm. The Pitzer equations for the major components of seawater have been used to estimate the density and compressibility of seawater to 95 °C. The results are in reasonable agreement with the measured values (0.010E?03 g cm^?3 for density and 0.050E?06 bar^?1 for compressibility) from 0 to 80 °C and salinities from 0 to 45 g kg^?1. The results make it possible to estimate the density and compressibility of all natural waters of known composition over a wide range of temperature and salinity.     

On the nonstoichiometry and point defects of olivine

This paper presents the point-defect thermodynamics for fayalite and olivine solid solutions (Fe _x Mg_1?x )₂SiO₄. By means of thermogravimetry, the metal-to-oxygen ratio of these silicates has been determined as a function of oxygen potential, compositionx and temperature. Experiments were performed in the range of 1,000° C≦T≦1,280° C and 0.2≦x≦1.0. It is found that V _Me ^″ , Fe _Me ^· and the associate {Fe′ _Si ^′ Fe _Me ^· } are the majority defects. With this knowledge it is possible to calculate the nonstoichiometry at given temperature as a function of \(p_{O_2 } \) and \(a_{SiO_2 } \) . The cation vacancy concentration shows a \(p_{O_2 }^{1/5} \) -dependence (forx≧0.2) and increases at givenT and \(p_{O_2 } \) almost exponentially with compositionx. In the composition range studied here, the silicates show an oxygen excess, and FeO is more soluble in the olivine than SiO₂.     

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.     

Kinetic rate laws derived from order parameter theory IV: Kinetics of Al,Si disordering in Na-feldspars

The kinetic rate laws of Al-Si disordering under dry conditions (T = 1353K, 1253 K, 1223 K, 1183 K) and in the presence of water (p = 1 kbar, T = 1023 K, 1073 K, 1103 K) were studied both experimentally and theoretically. A gradual change of the degree of order was found under dry conditions. For intermediate degrees of order broad distributions of the order parameter Q _od occur. The variations of Q _od are correlated with structural modulations as observed in the transmission electron microscope. The time evolution of the mean value of Q _od can be well described by the rate law: $$\frac{{dQ_{od} }}{{dt}} = - \frac{\gamma }{{RT}}\exp \sum\limits_{i = 1}^n {X_i^2 } \left[ {\frac{{ - (G_a^0 + \varepsilon (\Delta Q_{od} )^2 )}}{{RT}}} \right]\frac{{dG}}{{dQ_{od} }}$$ with the excess Gibbs energy G and G _a ⁰ = 433.8 kJ/mol, ?= -27.4 kJ/mol, γ = 1.687 · 10¹⁴ h ^?1. Under wet conditions, two processes were found which occur simultaneously. Firstly, some material renucleated with the equilibrium degree of order. Secondly, the bulk of the material transformed following the same rate law as under dry conditions but with the reduced activation energy G _a ⁰ = 332.0 kJ/mol and ? = -43.0 kJ/ mol, γ = 1.047 · 10¹³ h^?1. The applicability of the kinetic theory is discussed and some ideas for the analysis of geological observations are evolved.     

Thermodynamic data from redox reactions at high temperatures. II. The MnO-0Mn3O4 oxygen buffer,and implications for the thermodynamic properties of MnO and Mn3O4

The \(\mu _{O_2 } \) defined by the reaction 6 MnO+O₂ =2 Mn₃O₄ has been determined from 917 to 1,433 K using electrochemical cells (with calcia-stabilized zirconta, CSZ) of the type: Steady emfs were achieved rapidly at all temperatures on both increasing and decreasing temperature, indicating that the MnO-Mn₃O₄ oxygen buffer equilibrates relatively easily. It therefore makes a useful alternative choice in experimental petrology to Fe₂O₃-Fe₃O₄ for buffering oxygen potentials at oxidized values. The results are (in J/mol, temperature in K, reference pressure 1 bar); \(\mu _{O_2 } \) (±200)=-563,241+1,761.758T-220.490T inT+0.101819T ² with an uncertainty of ±200 J/mol. Third law analysis of these data, including a correction for the deviations in stoichiometry of MnO, impliesS _298.15 for Mn₃O₄ of 166.6 J/K · mol, which is 2.5 J/K · mol higher than the calorimetric determination of Robie and Hemingway (1985). The low value of the calorimetric entropy may be due to incomplete ordering of the magnetic spins. The third law value of Δ _r H _298.15 ⁰ is-450.09 kJ/mol, which is significantly different from the calorimetric value of-457.5±3.4 kJ/mol, calculated from Δ _f H _298.15 ⁰ of MnO and Mn₃O₄, implying a small error in one or both of these latter.     

An experimental study of Ti and Zr partitioning among zircon, rutile, and granitic melt

In order to evaluate the effect of trace and minor elements (e.g., P, Y, and the REEs) on the high-temperature solubility of Ti in zircon (zrc), we conducted 31 experiments on a series of synthetic and natural granitic compositions [enriched in TiO₂ and ZrO₂; Al/(Na + K) molar ~1.2] at a pressure of 10 kbar and temperatures of ~1,400 to 1,200 °C. Thirty of the experiments produced zircon-saturated glasses, of which 22 are also saturated in rutile (rt). In seven experiments, quenched glasses coexist with quartz (qtz). SiO₂ contents of the quenched liquids range from 68.5 to 82.3 wt% (volatile free), and water concentrations are 0.4–7.0 wt%. TiO₂ contents of the rutile-saturated quenched melts are positively correlated with run temperature. Glass ZrO₂ concentrations (0.2–1.2 wt%; volatile free) also show a broad positive correlation with run temperature and, at a given T, are strongly correlated with the parameter (Na + K + 2Ca)/(Si·Al) (all in cation fractions). Mole fraction of ZrO₂ in rutile $ \left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) $ in the quartz-saturated runs coupled with other 10-kbar qtz-saturated experimental data from the literature (total temperature range of ~1,400 to 675 °C) yields the following temperature-dependent expression: $ {\text{ln}}\left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) + {\text{ln}}\left( {a_{{{\text{SiO}}_{2} }} } \right) = 2.638(149) - 9969(190)/T({\text{K}}) $ , where silica activity $ a_{{{\text{SiO}}_{2} }} $ in either the coexisting silica polymorph or a silica-undersaturated melt is referenced to α-quartz at the P and T of each experiment and the best-fit coefficients and their uncertainties (values in parentheses) reflect uncertainties in T and $ \mathop X\nolimits_{{{\text{ZrO}}_{2} }}^{\text{rt}} $ . NanoSIMS measurements of Ti in zircon overgrowths in the experiments yield values of ~100 to 800 ppm; Ti concentrations in zircon are positively correlated with temperature. Coupled with values for $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ for each experiment, zircon Ti concentrations (ppm) can be related to temperature over the range of ~1,400 to 1,200 °C by the expression: $ \ln \left( {\text{Ti ppm}} \right)^{\text{zrc}} + \ln \left( {a_{{{\text{SiO}}_{2} }} } \right) - \ln \left( {a_{{{\text{TiO}}_{2} }} } \right) = 13.84\left( {71} \right) - 12590\left( {1124} \right)/T\left( {\text{K}} \right) $ . After accounting for differences in $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ , Ti contents of zircon from experiments run with bulk compositions based on the natural granite overlap with the concentrations measured on zircon from experiments using the synthetic bulk compositions. Coupled with data from the literature, this suggests that at T ≥ 1,100 °C, natural levels of minor and trace elements in “granitic” melts do not appear to influence the solubility of Ti in zircon. Whether this is true at magmatic temperatures of crustal hydrous silica-rich liquids (e.g., 800–700 °C) remains to be demonstrated. Finally, measured $ D_{\text{Ti}}^{{{\text{zrc}}/{\text{melt}}}} $ values (calculated on a weight basis) from the experiments presented here are 0.007–0.01, relatively independent of temperature, and broadly consistent with values determined from natural zircon and silica-rich glass pairs.     

Using chemical analysis and modeling to enhance the understanding of soil solution of some calcareous soils

The chemistry of soil solutions can be altered by human activities, due to the intense agricultural and husbandry, leading to leaching of nutrients and subsequently elevating ground water levels. Multivariate statistical and inverse geochemical modeling techniques were used to determine the main factors controlling soil solution chemistry of calcareous soils. In this research, a total of 21 calcareous soils was characterized and assessed for soil solution using soil column. The major cations in the studied soil solutions were in the decreasing order as Ca²⁺ > Mg²⁺ > Na⁺ > K⁺. The anions were also arranged in decreasing order as HCO $ _{3}^{ - } $  > Cl $ ^{ - } $  > SO $ _{4}^{2 - } $  > NO $ _{3}^{ - } $ . Concentrations of NO $ _{3}^{ - } $ , P, and K⁺ in soil solutions were in the range of 6.8–307.5 mg l^?1 (mean 63.2 mg l^?1), 5.0–10.4 mg l^?1 (mean 5.9 mg l^?1), and 2.8–54.6 mg l^?1 (mean 11.3 mg l^?1), respectively. Results suggest that the concentration of P in the soil solutions could be primarily controlled by the solubility of dicalcium phosphate dihydrate and dicalcium phosphate. Interactions between soil properties and observed solubility of nutrients were described, and put into empirical multivariate formulations. Obtained equations contained electrical conductivity (EC) as a key factor in determining nutrients solubility. Inverse geochemical modeling of soil solution using PHREEQC indicates the dissolution of calcite, anhydrite, halite, CO₂ (g), N₂ (g), and hydroxyapatite, and precipitation of sulfur. Cation exchange between Ca²⁺, Mg²⁺, K⁺ and Na⁺ occurred with Mg²⁺ and K⁺ into the solution, and Ca²⁺ and Na⁺ out of the solution. Determination of soil solution will improve soil management in the area, and preventing groundwater deterioration.     

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

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

The high-temperature heat capacity of natural calcite (CaCO3)

The heat capacity (C _P) of a natural sample of calcite (CaCO₃) has been measured from 350 to 775 K by differential scanning calorimetry (DSC). Heat capacities determined for a powdered sample and a single-crystal disc are in close agreement and have a total uncertainty of ±1 percent. The following equation for the heat capacity of calcite from 298 to 775 K was fit by least squares to the experimental data and constrained to join smoothly with the low-temperature heat capacity data of Staveley and Linford (1969) (C _P in J mol^?1 K^?1, T in K): $$\begin{gathered} C_p = - 184.79 + 0.32322T - 3,688,200T^{ - 2} \hfill \\ {\text{ }} - (1.2974{\text{ }} \times {\text{ 10}}^{ - {\text{4}}} )T^2 + 3,883.5T^{ - 1/2} \hfill \\ \end{gathered} $$ Combining this equation with the S ₂₉₈ ⁰ value from Staveley and Linford (1969), entropies for calcite are calculated and presented to 775 K. A simple method of extrapolating the heat capacity function of calcite above 775 K is presented. This method provides accurate entropies of calcite for high-temperature thermodynamic calculations, as evidenced by calculation of the equilibrium: CaCO_{3 (s)}=CaO_(s)+CO_{2 (s)}.     

The ferric-ferrous ratio of natural silicate liquids equilibrated in air

Results of chemical analyses of glasses produced in 46 melting experiments in air at 1,350° C and 1,450° C on rocks ranging in composition from nephelinite to rhyolite have been combined with other published data to obtain an empirical equation relating in \((X_{{\text{Fe}}_{\text{2}} {\text{O}}_{\text{3}} }^{{\text{liq}}} /X_{{\text{FeO}}}^{{\text{liq}}} )\) to T, \(\ln f_{{\text{O}}_{\text{2}} } \) and bulk composition. The whole set of experimental data range over 1,200–1,450° C and oxygen fugacities of 10^?9.00 to 10^?0.69 bars, respectively. The standard errors of temperature and \(\log _{10} f_{{\text{O}}_{\text{2}} } \) predictions from this equation are 52° C and 0.5 units, respectively, for 186 experiments.     

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} $$      

On the least-squares fit of small and great circles to spherically projected orientation data

Given the direction cosines a _i ^′ = (a ₁ ⁱ , a ₂ ⁱ , a ₃ ⁱ )corresponding to a set of pspherically projected fabric poles, an initial estimate x′ = (x₁, x₂, x₃, x₄)for the angular radius x₄,and direction cosines of the center of the least-squares small circle which minimizes the sum of the squares of the angular residuals $$r = \sum\limits_p {\left[ {x_4 - \cos ^{ - 1} \left( {a_1^i x_1 + a_2^i x_2 + a_3^i x_3 } \right)} \right]} ^2 $$ can be iteratively improved by taking x_j+1 = x_j + Δxwhere x_j is the value of xat the jth iteration and $$\Delta x = - H_j^{ - 1} \left[ {q_j + x_j \left( {x'_j H_j^{ - 1} x_j } \right)\left( {q_j - x'_j H_j^{ - 1} q_j } \right)} \right],$$ where As an initial approximation for xwe have found it convenient to ignore the fact that the data are constrained to lie on the surface of the reference sphere and to use the parameters of a least-squares plane through the given poles. Generalization of this approach to fitting variously constrained great and small circles is easily made. The relative merits of differently constrained fits to the same data can be tested approximately if it is assumed that the errors in the location of the poles are isotropic and normally distributed. It is thus possible to statistically assess the relative significance of conflicting structural models which predict different geometrical patterns of fabric elements.