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.     

The experimental investigation of the reactions MnCO3+SiO2=MnSiO37#x002B;CO2 and MnSiO3+MnCO3=Mn2SiO4+CO2 in CO2/H2O gasmixtures at a total pressure of 500 bars

The equilibrium constants for the reaction (2) Rhodochrosite + Quartz=Pyroxmangite+CO₂ obtained are:logK(₂)(bars)= $$\begin{gathered}{\text{log}}f_{co_2 } = - \frac{{(9862 \pm 102)}}{T} \hfill \\+ (15.887 \pm 0.220) + (0.1037 \pm 0.0020)\frac{{P - 1}}{T} \hfill \\\end{gathered} $$ and for the reaction (3) Rhodochrosite+Pyroxmangite=Tephroite+CO₂: logK(₃)(bars)= $$\begin{gathered}{\text{log}}f_{co_2 } = - \frac{{(6782 \pm 205)}}{T} \hfill \\+ (11.296 \pm 0.304) + (0.0835 \pm 0.0030)\frac{{P - 1}}{T} \hfill \\\end{gathered} $$ The present data lie within reasonable limits of error of the values calculated from previous experimental results at P _tot = 2000 bars.     

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.     

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”.     

An experimental study of the partitioning of Fe and Mg between garnet and orthopyroxene

The partitioning of Fe and Mg between garnet and aluminous orthopyroxene has been experimentally investigated in the pressure-temperature range 5–30 kbar and 800–1,200° C in the FeO-MgO-Al₂O₃-SiO₂ (FMAS) and CaO-FeO-MgO-Al₂O₃-SiO₂ (CFMAS) systems. Within the errors of the experimental data, orthopyroxene can be regarded as macroscopically ideal. The effects of Calcium on Fe-Mg partitioning between garnet and orthopyroxene can be attributed to non-ideal Ca-Mg interactions in the garnet, described by the interaction term:W _CaMg ^ga -W _CaFe ^ga =1,400±500 cal/mol site. Reduction of the experimental data, combined with molar volume data for the end-member phases, permits the calibration of a geothermometer which is applicable to garnet peridotites and granulites: $$T(^\circ C) = \left\{ {\frac{{3,740 + 1,400X_{gr}^{ga} + 22.86P(kb)}}{{R\ln K_D + 1.96}}} \right\} - 273$$ with $$K_D = {{\left\{ {\frac{{Fe}}{{Mg}}} \right\}^{ga} } \mathord{\left/ {\vphantom {{\left\{ {\frac{{Fe}}{{Mg}}} \right\}^{ga} } {\left\{ {\frac{{Fe}}{{Mg}}} \right\}}}} \right. \kern-\nulldelimiterspace} {\left\{ {\frac{{Fe}}{{Mg}}} \right\}}}$$ and $$X_{gr}^{ga} = (Ca/Ca + Mg + Fe)^{ga} .$$ The accuracy and precision of this geothermometer are limited by largerelative errors in the experimental and natural-rock data and by the modest absolute variation inK _D with temperature. Nevertheless, the geothermometer is shown to yield reasonable temperature estimates for a variety of natural samples.     

Oxygen isotope fractionation between rutile and water

Synthetic rutile-water fractionations (1000 ln α) at 775, 675, and 575° C were found to be ?2.8, ?3.5, and ?4.8, respectively. Partial exchange experiments with natural rutile at 575° C and with synthetic rutile at 475° C failed to yield reliable fractionations. Isotopic fractionation within the range 575–775° C may be expressed as follows: 1 $$1000\ln \alpha ({\rm T}i{\rm O}_{2 } - H_2 O) = - 4.1 \frac{{10^6 }}{{T_{k^2 } }} + 0.96$$ . Combined with previously determined quartz-water fractionations, the above data permit calibration of the quartz-rutile geothermometer: 1 $$1000\ln \alpha ({\text{S}}i{\rm O}_{2 } - Ti{\rm O}_{2 } ) = 6.6 \frac{{10^6 }}{{T_{k^2 } }} - 2.9$$ . When applied to B-type eclogites from Europe, as an example, the latter equation yields a mean equilibration temperature of 565° C.     

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 as derived from order parameter theory II: Interpretation of experimental data by laplace-transformation,the relaxation spectrum,and kinetic gradient coupling between two order parameters

A unifying theory of kinetic rate laws, based on order parameter theory, is presented. The time evolution of the average order parameter is described by $$\langle Q\rangle \propto \smallint P(x)e^{^{^{^{^{^{^{ - xt} } } } } } } dx = L(P)$$ where t is the time, x is the effective inverse susceptibility, and L indicates the Laplace transformation. The probability function P(x) can be determined from experimental data by inverse Laplace transformation. Five models are presented:

1. Polynomial distributions of P(x) lead to Taylor expansions of 〈Q〉 as $$\langle Q\rangle = \frac{{\rho _1 }}{t} + \frac{{\rho _2 }}{{t^2 }} + ...$$
2. Gaussian distributions (e.g. due to defects) lead to a rate law $$\langle Q\rangle = e^{ - x_0 t} e^{^{^{^{^{\frac{1}{2}\Gamma t^2 } } } } } erfc\left( {\sqrt {\frac{\Gamma }{2}} t} \right)$$ where x ₀ is the most probable inverse time constant, Γ is the Gaussian line width and erfc is the complement error integral.
3. Maxwell distributions of P are equivalent to the rate law 〈Q〉∝e^?k√t .
4. Pseudo spin glasses possess a logarithmic rate law 〈Q〉∝lnt.
5. Power laws with P(x)=x ^a lead to a rate law: ln〈Q〉=-(α + 1) ln t.

The power spectra of Q are shown for Gaussian distributions and pseudo spin glasses. The mechanism of kinetic gradient coupling between two order parameters is evaluated.     

A possible relation between the masses of black holes in galactic nuclei and the parameters of the host galaxies

Data on about forty virialized galaxy clusters with bright central galaxies, for which both the galactic velocity dispersion (?? _gal) and the stellar velocity dispersion in the brightest galaxies (??*) are measured, have been used to obtain several approximate relations between ?? _gal, ??*, the absolute B magnitude of the brightest central galaxyM _B ^BCG , and the mass of the central massive black holeM _BH: $\begin{gathered} \log \sigma _* = (0.12 \pm 0.14)\log \sigma _{gal} + (2.1 \pm 0.4), \hfill \\ \log \sigma _* = - (0.15 \pm 0.02)M_B^{BCG} + (0.85 \pm 0.5), \hfill \\ \log M_{BH} = 0.51\log \sigma _{gal} + 7.28. \hfill \\ \end{gathered} $ . These relations can be used to derive crude estimates ofMBH in the nuclei of the brightest galaxies using the parameters of the both host galaxies and the host galaxy clusters. The last relation above confirms earlier suggestions of a quadratic relation between the masses of the coronas of the host systems and the masses their central objects: M _hg ^halo ?? M _cent ² . The relations obtained are consistent with the common evolution of subsystems with different scales and masses formed in the process of hierarchical clustering.     

The influence of fluid unmixing on cation exchange between plagioclases and aqueous chloride solutions at 700° C, 1 kbar

Experimental exchanges between plagioclases (synthesized from gels) and aqueous solutions (0.5N–8N) were carried out according to the reaction $$\begin{gathered} 2NaA1Si_3 O_8 + CaC1_2 \hfill \\ \leftrightarrow CaA1_2 Si_2 O_8 + 4SiO_2 + 2NaC1. \hfill \\ \end{gathered}$$ Distribution coefficients defined by $$K_D = \frac{{X_{An} }}{{(X_{Ab} )^2 }}\frac{{(X_{NaC1} )^2 }}{{X_{CaC1_2 } }}$$ were determined at 700° C and 1 kbar. From previous studies it is known that variations in the concentration of the aqueous solutions have no influence upon K _D if the fluid is a single phase. In this study, variation of K _D with the concentration of the solutions is interpreted as the result of fluid unmixing to vapour and brine phases. This implies boiling of CaCl₂-NaCl-H₂O fluids analogous to that known for the system NaCl-H₂O. Experimental data permit calculation of the compositions of vapours and estimation of those of the brines for fluids in which Ca/Na<0.5. Boiling has an effect upon the exchange between feldspars and solutions (metasomatism) and must be considered when determining the activity coefficients.     

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.     

The role of the minor substitutions in the crystal structure of natural challacolloite, KPb2Cl5, and hephaistosite, TlPb2Cl5, from Vulcano (Aeolian Archipelago, Italy)

Crystals of challacolloite, KPb₂Cl₅, and hephaistosite, TlPb₂Cl₅, from volcanic sublimates formed on the crater rim of the “La Fossa Crater” at Vulcano, Aeolian Archipelago, Italy, were investigated. Chemical compositions were ${\left( {{\text{K}}_{{0.93}} {\text{Tl}}_{{0.02}} } \right)}_{{\Sigma = 0.95}} {\text{Pb}}_{{2.04}} {\left( {{\text{Cl}}_{{4.90}} {\text{Br}}_{{0.11}} } \right)}_{{\Sigma = 5.01}} $ and ${\text{Tl}}_{{0.94}} {\text{Pb}}_{{2.01}} {\left( {{\text{Cl}}_{{4.91}} {\text{Br}}_{{0.14}} } \right)}_{{\Sigma = 5.05}} $ , respectively. Single crystal X-ray measurements showed monoclinic symmetry for both phases, space group P2₁/c, with the following unit-cell parameters: a = 8.8989(4), b = 7.9717(5), c = 12.5624(8) Å, β = 90.022(4)°, V = 891.2(1) ų, Z = 4 (challacolloite) and a = 9.0026(6), b = 7.9723(6), c = 12.5693(9) Å, β = 90.046(4)°, V = 902.1(1) ų, Z = 4 (hephaistosite). The structure refinements converge to R = 3.99% and R = 3.86%, respectively. The effects of Br?Cl and K?Tl substitutions on the structure of these natural compounds have been discussed.     

The thermodynamic properties of reciprocal solid solutions

A multisite solid solution of the type (A, B) (X, Y) has the four possible components AX, AY, BX, BY. Taking the standard state to be the pure phase at the pressure and temperature of interest, the mixing of these components is shown not to be ideal unless the condition: $$\Delta G^0 = (\mu _{AX}^0 + \mu _{BY}^0 - \mu _{AY}^0 - \mu _{BX}^0 = 0$$ applies. Even for the case in which mixing on each of the individual sublattices is ideal, ΔG ⁰ contributes terms of the following form to the activity coefficients of the constituent components: $$RT\ln \gamma _{AX} = - X_{B_1 } X_{Y_2 } \Delta G^0$$ (X _Ji refers to the atomic fraction of J on sublattice i). The above equation, which assumes complete disorder on (A, B) sites and on (X, Y) sites is extended to the general n-component case. Methods of combining the “cross-site” or reciprocal terms with non-ideal terms for each of the individual sites are also described. The reciprocal terms appear to be significant in many geologically important solid solutions, and clinopyroxene, garnet and spinel solid solutions are all used as examples. Finally, it is shown that the assumption of complete disorder only applies under the condition: $$\Delta G^0 \ll zn_1 RT$$ where z is the number of nearest-neighbour (X, Y) sites around A and n ₁ is the number of (A, B) sites in the formula unit. If ΔG ⁰ is relatively large, then substantial short range oder must occur and the activity coefficient is given by (ignoring individual site terms): $$\gamma _{AX} = \left( {\frac{{1 - X'_{Y_2 } }}{{1 - X_{Y_2 } }}} \right)^{zn_1 }$$ where X′_Y2 is the equilibrium atomic fraction of Y atoms surrounding A atoms in the structure. The ordered model may be developed for multicomponent solutions and individual site interactions added, but numerical methods are needed to solve the simultaneous equations involved.     

Peculiarities of radio pulsars with emission outside the radio range

Parameters of 100 radio pulsars detected outside the radio range (he pulsars) are compared with those of pulsars radiating only in the radio (n pulsars). The periods of he pulsars are, on average, appreciably shorter than those of n pulsars: 〈P〉 = 0.10 and 0.56 s, respectively. The distribution of the magnetic field at the light cylinder is shifted toward higher magnetic fields for the pulsars with high-energy radiation, compared to the distribution for pulsars radiating only in the radio. The magnetic fields at the light cylinder are 〈B _lc〉 = 9×10³ G for he radio pulsars, and 〈B_lc〉 = 56 G formost purely radio pulsars. This suggests the generation of high-energy nonthermal radiation in radio pulsars at the peripheries of their magnetospheres. The distribution of the spin-energy loss rate dE/dt is uniform for he pulsars, and is characterized by a higher average value \(\left( {\left\langle {\log \frac{{dE}} {{dt}}} \right\rangle = 35.53} \right) \) , compared to n pulsars, \(\left( {\left\langle {\log \frac{{dE}} {{dt}}} \right\rangle = 32.60} \right) \) . The spatial distribution of he pulsars is nonuniform: they form two well separated clouds.     

Experimental investigation of the metamorphism of siliceous dolomites

For the reaction: 1 diopside+3 dolomite ?2 forsterite+4 calcite+2 CO₂ (14) the following P _total?T-brackets have been determined experimentally in the presence of a gasphase consisting of 90 mole%CO₂ and 10 mole%H₂O∶1 kb, 544°±20° C; 3kb, 638°±15° C; 5kb, 708°±10° C; 10kb, 861°±10° C. The determination was carried out with well defined synthetic minerals in the starting mixture. The MgCO₃-contents of the magnesian calcites formed by the reaction in equilibrium with dolomite agree very well with the calcite-dolomite miscibility gap, which can be recalculated from the activities and the activity coefficients of MgCO₃ as given by Gordon and Greenwood (1970). The equilibrium constant K _14b was calculated with respect to the reference pressure P ₀=1 bar using the experimentally determined \(P_{total} TX_{CO_2 }\) brackets, the activities of MgCO₃ and CaCO₃ (Gordon and Greenwood 1970; Skippen 1974) and the fugacities of CO₂ Holloway (1977) considering the correction of Flowers (1979). Results are plotted as function of the absolute reciprocal temperature in Fig. 1. For the temperature range of 530° to 750° C the following linear expression can be given for the natural logarithm of K_14b: (g) $$[ln K_{14b} ]_T^P = - \frac{{18064.43}}{{T\left( {^\circ K} \right)}} + 38.58 + \frac{{0.308(P - 1 bar)}}{{T\left( {^\circ K} \right)}}$$ where P is the total pressure in bars and T the temperature in degrees Kelvin. Combining Equation (g) with the activities of MgCO₃ and CaCO₃ gives the equilibrium fugacity \(f_{CO_2 }\) : (i) $$[ln f_{CO_2 } ]_T^P = - \frac{{11635.44}}{{T\left( {^\circ K} \right)}} + 21.09 + \frac{{0.154(P - 1 bar)}}{{T\left( {^\circ K} \right)}}$$ Equation (i) and the fugacities of CO₂ permit to calculate the equilibrium data in terms of \(P_{CO_2 }\) and T (see Fig. 3) or P _total, T and \(X_{CO_2 }\) (see Fig. 5). Combining the \(P_{total} TX_{CO_2 }\) equilibrium data of the above reaction with those of the previously investigated reaction (Metz 1976): 1 tremolite+11 dolomite ?8 forsterite+13 calcite+9 CO₂+1 H₂O yields the stability conditions of the four-mineral assemblage: diopside+calcian dolomite+forsterite +magnesian calcite and the stability conditions of the five-mineral assemblage: tremolite+calcian dolomite+forsterite +magnesian calcite+diopside both shown in Fig. 6. Since these assemblages are by no means rare in metamorphic siliceous dolomites (Trommsdorff 1972; Suzuki 1977; Puhan 1979) the data of Fig. 6 can be used to determine the pressure of metamorphism and to estimate the composition of the CO₂-H₂O fluid provided the temperature of the metamorphic event was determined using the calcite-dolomite geothermometer.     

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.     

Apsidal motion in the eclipsing binary GG Ori

High-precision WBVR photoelectric observations of the eclipsing binary GG Ori (B9.5V+B9.5V), which has an eccentric orbit (e=0.22), were carried out in 1988–2001 at the Moscow and high-altitude Tian-Shan Observatories of the Sternberg Astronomical Institute. The aim of these observations was investigation of the apsidal motion of the system. Analysis of the resulting 12-year series of observations enabled us for the first time to accurately (to within 11%) measure the rate of rotation of the orbit $\dot \omega _{obs} = 0.046 \pm 0.005^\circ /yr$ and to appreciably improve estimates of the photometric and absolute parameters. The observed value of $\dot \omega _{obs}$ is 28% higher than the theoretical prediction of $\dot \omega _{th} = \dot \omega _{cl} + \dot \omega _{rel} = 0.036 \pm 0.001^\circ /yr$ . The relativistic part of the apsidal motion in this system $\dot \omega _{rel}$ is a factor of 2.5 greater than the classical term $\dot \omega _{cl}$ due to the tidal and rotational deformations of the components. The interstellar extinction in the direction of the star (at a distance of r=425 pc) is very large (A _v =1.75 ^m ). A number of recently published results (in particular, the conclusion that the components of this eclipsing binary are young) are confirmed.     

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.     

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.     

Calcic amphibole equilibria and a new amphibole-plagioclase geothermometer

There is currently a dearth of reliable thermobarometers for many hornblende and plagioclase-bearing rocks such as granitoids and amphibolites. A semi-empirical thermodynamic evaluation of the available experimental data on amphibole+plagioclase assemblages leads to a new thermometer based on the Al^iv content of amphibole coexisting with plagioclase in silica saturated rocks. The principal exchange vector in amphiboles as a function of temperature in both the natural and experimental studies is \(\left( {Na\square _{ - 1} } \right)^A \left( {AlSi_{ - 1} } \right)^{T1}\) . We have analysed the data using 3 different amphibole activity models to calibrate the thermometer reactions 1. $$1. Edenite + 4 Quartz = Tremolite + Albite$$ 2. $$2. Pargasite + 4 Quartz = Hornblende + Albite.$$ The equilibrium relation for both (1) and (2) leads to the proposed new thermometer $$T = \frac{{0.677P - 48.98 + Y}}{{ - 0.0429 - 0.008314 ln K}} and K = \left( {\frac{{Si - 4}}{{8 - Si}}} \right)X_{Ab}^{Plag} ,$$ where Si is the number of atoms per formula unit in amphiboles, with P in kbar and T in K; the term Y represents plagioclase non-ideality, RTlnγ_ab, from Darken's Quadratic formalism (DQF) with Y=0 for X _ab>0.5 and Y=-8.06+25.5(1-X _ab)² for X _ab<0.5. The best fits to the data were obtained by assuming complete coupling between Al on the T1 site and Na in the A site of amphibole, and the standard deviation of residuals in the fit is ±38°C. The thermometer is robust to ferric iron recalculation procedures from electron probe data and should yield temperatures of equilibration for hornblende-plagioclase assemblages with uncertainties of around ±75° C for rocks equilibrated at temperatures in the range 500°–1100° C. The thermometer should only be used in this temperature range and for assemblages with plagioclase less calcic than An₉₂ and with amphiboles containing less than 7.8 Si atoms pfu. Good results have been attained on natural examples from greenschist to granulite facies metamorphic rocks as well as from a variety of mafic to acid intrusive and extrusive igneous rocks. Our analysis shows that the pressure dependence is poorly constrained and the equilibria are not suitable for barometry.