Some constraints on the thermodynamic mixing properties of Fe-Mg orthopyroxenes and olivines

Existing data on the temperature and composition dependence of the Fe²⁺-Mg²⁺ distribution between Fe-Mg olivine and orthopyroxene, the intra-crystalline distribution of Fe²⁺ and Mg²⁺ between M1 and M2 sites in orthopyroxene, and macroscopic activity-composition relations in olivine and orthopyroxene are shown to be inconsistent with generally accepted thermodynamic formulations which assume that the non-configurational Gibbs energy of orthopyroxene is independent of the degree of long-range ordering of Fe²⁺ and Mg⁺ between M1 and M2 sites. These data are interpreted in terms of the constraints they provide on the size of Bragg-Williams type energy, entropy, and volume terms for olivine and orthopyroxene. The apparent equilibrium constant for Fe-Mg exchange between olivine and orthopyroxene is shown to be a potentially useful ‘geothermometer’ for olivine-orthopyroxene assemblages with olivines with mole fraction of Fe₂SiO₄ component less than 0.2 or greater than 0.6. A provisional calibration of this ‘geothermometer’ is presented.     

Variation in ferrous iron-magnesium distribution coefficients of metamorphic pyroxenes from Quairading,Western Australia

Distribution coefficients (K _D·Fe ⁺⁺ _–Mg) calculated for orthopyroxene-clinopyroxene pairs from 12 basic granulites of the Quairading district, Western Australia, range from 1.87 for magnesian pyroxenes (Opx Mg value=78.1) to 1.70 for iron-rich varieties (Opx Mg value 37.7). Field and petrographic evidence indicates that these pyroxenes have probably reached equilibrium within a narrow temperature range. In order to account for the observed variation of K _D values it is suggested that one (or both) of the pyroxene structures is not the ideal Fe⁺⁺-Mg solution proposed in the thermodynamic model of the pyroxene equilibrium exchange. After consideration of the geometry of the pyroxene cation sites, the relative bond energies of each site (especially crystal field stabilization energy) and structural ordering a model is proposed to explain the non-ideal behaviour of Fe⁺⁺-Mg in the pyroxene system. The distribution pattern in low-iron pyroxenes will probably show Fe⁺⁺ favouring the M _2·Opx site; competing unfavourably with Mg⁺⁺ for the M _1·Cpxsite; and probably excluded by Mg⁺⁺ from the M _1·Opxsite. As the iron content of the system increases the M _2·Opxsite will begin to become saturated with Fe⁺⁺ and this ion will enter the M _1·Opx site. Further increases in the iron content of the system will cause the Fe⁺⁺-Mg distribution to depend increasingly on the relative attraction of the M ₁ sites of both pyroxenes. Of these sites Fe⁺⁺ should show preference to the more distorted M _1·Cpxsite. The distribution coefficient reflects this swing towards a relative enrichment of Fe⁺⁺ in the clinopyroxene by decreasing regularly with increasing iron content. It is likely that this downward trend will not become evident until the iron content of the M _2·Opx site reaches saturation. This would explain why the K _D values for the magnesian pyroxenes remain practically unchanged until the orthopyroxene Mg value is approximately 60; from here on the iron-rich pyroxenes show a rapid decrease in K _D value with increasing iron content.The Ca content of the pyroxenes is also significant since the Quairading pyroxenes show a marked increase in mutual solubility with increasing iron content. Calcium taken into the orthopyroxene structure will enter the M ₂ site ahead of Fe⁺⁺ so that this site will reach Fe⁺⁺ saturation at a lower iron content than when the orthopyroxene is Ca-free.The application of K _D values to the regional study of basic granulites, particularly when establishing relative temperature zones on the basis of K _D variation, should only be attempted when pyroxenes which extend over a wide range of Fe-Mg content are available.     

The composition,structure, and stability of guinier-preston zones in lunar and terrestrial orthopyroxene

Lunar and terrestrial orthopyroxenes (Mg,Fe,Ca)₂Si₂O₆ contain varying abundances of coherent, Ca-enriched Guinier-Preston (G.P.) zones. G.P. zones 5–6 unit cells thick have been found in one lunar sample whereas all other examples (lunar and terrestrial) are only one unit-cell-thick. Electron diffraction maxima from the larger lunar G.P. zones indicate that d ₁₀₀=18.52 Å whereas, d ₁₀₀=18.2 Å for the host. This increase in the a direction corresponds to an increase in calcium content in the G.P. zones over that of the host of ～25 mol% Ca₂Si₂O₆. Diffraction patterns of the hk0 net from an area containing G.P. zones show extra spots (h=2n+1) not observed in the host orthopyroxene (Pbca), that violate the a-glide of the host. The G.P. zones, therefore, have space group Pbc2₁ if it is assumed that the c-glide of pyroxene is retained and the space group of the G.P. zone is a subgroup of Pbca. The loss of the a-glide in the G.P. zones results in 4 distinct silica chains and 4 distinct cation sites M1A, M1B, M2A, M2B; by symmetry, equivalent M2A or M2B sites are clustered together in only one-half of the unit cell. As one-fourth of the divalent cations in the G.P. zones are calcium, ordering of Ca on M2A or M2B would produce a zone 9 Å thick extended parallel to (100) with the composition of Ca(Mg,Fe)Si₂O₆, but constrained by the host to the structure of orthopyroxene. This zone and the Ca-poor half-unit-cell then constitute an 18 Å thick G.P. zone. Heating experiments of varying duration indicate that the zones become unstable with respect to the host orthopyroxene at ～950°C for Wo_0.6 and ～1,050°C for Wo_2.5. The zones are interpreted in terms of the pyroxene subsolidus as a metastable phase having either a solvus relationship with orthopyroxene or originating as a distinct phase. The size, distribution, composition and structure of G.P. zones may be an important indicator of the low-temperature thermal history of orthopyroxene.     

Intracrystalline relationships in olivine,orthopyroxene, clinopyroxene and spinel from a suite of spinel lherzolite xenoliths from Mt. Noorat,Victoria, Australia

A detailed crystal chemical study of coexisting olivine, orthopyroxene, clinopyroxene and spinel from selected Victorian (Australia) lherzolite suites was carried out by means of single crystal x-ray diffraction and electron probe microanalysis to obtain actual site occupancies. The aim of this study was primarily to characterise the intracrystalline configurations and related cation ordering on sites in major mantle constituents. The results demonstrate that cation ordering on sites is subject to distinctive crystallographic controls which depend on the petrological evolution of the suite. Mg-Fe²⁺ ordering in M1–M2 pyroxene sites depends on variations of the smaller cations, mainly Al^vi, Ti⁴⁺, Fe³⁺, and related configurations of M 1. Pressuresensitive Al^vi is crucial to Fe²⁺, the more ordered clinopyroxene showing high Al^vi configurations which tend to exclude the larger bivalent cations and yield small polyhedral volumes for M 1, M 2, T sites and the unit cell. Conversely, the coexisting orthopyroxene, characterised by lower Al^vi configuration and higher M 1 and unit cell volumes, is relatively more disordered. Olivine is consistent with the coexisting clinopyroxene, the more disordered crystals coexisting with more disordered clinopyroxene, while Al-Mg order in the coexisting spinel shows the reverse relationship. Estimated temperatures of apparent equilibration based on current geothermometers are not considered realistic. Assumptions of ideal cation mixing on sites in pyroxene and spinel are not supported.     

The site of Fe in Fe-bearing MgSiO3 enstatite and perovskite. a theoretical-, x-ray multiple-scattering study at Fe K-edge

Two polycrystalline-, Fe-bearing MgSiO₃ enstatite and perovskite have been probed by x-ray absorption near edge structure (XANES) spectroscopy at the Fe K-edge under ambient conditions. The perovskite sample was synthesized at 260 kbar and 1973 K in a multianvil apparatus. The experimental XANES spectrum has been compared to ab-initio-, x-ray multiple-scattering calculations (Feff 6 code). Calculations confirm that the Fe K-edge arises mainly from multiple scattering involving the first shell of oxygen neighbors around Fe. In Fe-enstatite, these calculations are consistent with Fe²⁺ as substituted in the M2 site of this orthopyroxene, in good agreement with crystal structure refinements and previous XANES studies. In perovskite, Feff 6 suggests that Fe is likely to be substituted to Mg within the (8+4)-coordinated sites of that perovskite. No evidences for 6-coordinated Fe were found. These results are consistent with a previous anharmonic analysis of the extended x-ray absorption fine structure (EXAFS) study that evidenced the presence of 8-coordinated Fe in the same perovskite sample.     

Calorimetric determination of the enthalpy of Mg-Fe ordering in orthopyroxene

The enthalpy of Mg-Fe ordering in En₅₀Fs₅₀ orthopyroxene was measured using the transposed temperature drop calorimetric method. Heat effects associated with two consecutive drops were recorded. In the first drop, synthetic orthopyroxene samples equilibrated at 823?K, 0.1?MPa and a f?O₂ of the WI buffer were dropped from 823?K into the calorimeter, which was held at 1173?K. The measured heat effect corresponds to the enthalpy change due to the heat capacity of the sample from 823 to 1173?K and to the enthalpy associated with the (dis)ordering of Mg and Fe²⁺. In the second drop, the samples, with an Fe-Mg order corresponding to 1173?K, were dropped again from 823 to 1173?K. From the difference of the heat effects measured in the two experiments, the enthalpy of disordering associated with the temperature change from 823 to 1173?K was calculated to be ?1.73±0.04 J mol^?1. The observed enthalpy corresponds to a change in the mole fraction of iron on the M2 site, ΔX _Fe,M2=?0.096 ± 0.001, which leads to of ΔH ⁰ _exch of 18.0 ± 0.4 kJ mol^?1 for the exchange reaction: The degree of Fe-Mg order was characterized by ⁵⁷Fe Mössbauer resonance spectroscopy. In order to minimize the error due to the thickness of the absorber, the iron concentration of the absorber was reduced step by step from 5 to 1 mg?Fe?cm^?2. The iron distribution extrapolated to zero thickness was used for the calculations of the enthalpy of exchange reaction.     

Ion migration in MgSiO3-perovskite and olivine by molecular dynamics calculations

We applied molecular dynamics (MD) simulations to finding likely paths of atomic migration of the Mg ion in forsterite (Mg₂SiO₄) and the oxygen ion in MgSiO₃-perovskite to better understand atomic diffusion in minerals. Our simulations show that there exist two routes of Mg migration in the forsterite structure, that is, paths between the M1 and M2 sites and between the M1 and M1 sites. In the MgSiO₃-perovskite structure, some oxygen ions migrate to the next sites all together through the O1 vacant site showing co-operative movements. The O ions are relatively mobile mainly along the b axis in the perovskite structure. Meta-stable sites are often present between a stable site and another stable site on atomic migration. In spite of many assumptions, our MD simulations may show likely paths of atomic migration in crystal.     

Thermodynamics of element partitioning: (1) Systematics of transition metals in crystalline and molten silicates and (2) Defect chemistry and “the Henry's law problem”

The distribution of Mg, Mn, Fe, Co and Ni among olivine, orthopyroxene, calcic clinopyroxene and liquid can be described by exchange reactions of the form: M (phase A) + Mg (phase B) = M (phase B) + Mg (phase A). The thermochemical data predict the observed partitioning within the limits of error of both sets of data when assumptions of ideal solid solution are used, except for Mn-Mg exchange. Because ΔS and ΔV for these exchange reactions are generally small, K_D varies relatively little with temperature and pressure, although individual ion distribution coefficients (D values) are more sensitive to these variables and to changing liquid composition.Trace elements present at the ppm level can enter crystals in normal lattice sites, can participate in point defect equilibria and can enter various defect sites present metastably in the crystal as a result of its P,T history. Metastable defect equilibria, in which the number of sites (grain boundaries, surfaces, intergrowths, etc.) is fixed, can lead to enhanced trace element solubility in the crystal at very low concentrations and thus to apparent deviations from Henry's law.     

Importance of considerations of mixing properties in establishing an internally consistent thermodynamic database: thermochemistry of minerals in the system Mg2SiO4-Fe2SiO4-SiO2

A thermodynamic solution model is developed for minerals whose compositions lie in the two binary systems Mg₂SiO₄-Fe₂SiO₄ and Mg₂Si₂O₆-Fe₂Si₂O₆. The formulation makes explicit provision for nonconvergent ordering of Fe²⁺ and Mg²⁺ between M1 and M2 sites in orthopyroxenes and non-zero Gibbs energies of reciprocal ordering reactions in both olivine and orthopyroxene. The calibration is consistent with (1) constraints provided by available experimental and natural data on the Fe-Mg exchange reaction between olivine and orthopyroxene ± quartz, (2) site occupancy data on orthopyroxenes including both crystallographic refinements and Mössbauer spectroscopy, (3) enthalpy of solution data on olivines and orthopyroxenes and enthalpy of disordering data on orthopyroxene, (4) available data on the temperature and ordering dependence of the excess volume of orthopyroxene solid solutions, and (5) direct activity-composition determinations of orthopyroxene and olivine solid solutions at elevated temperatures. Our analysis suggests that the entropies of the exchange [Mg(M2)Fe(M1)Fe(M2)Mg(M1)] and reciprocal ordering reactions [Mg(M2)Mg(M1)+ Fe(M2)Fe(M1)Fe(M2)Mg(M1)+Mg(M2)Fe(M1)] cannot differ significantly (± 1 cal/K) from zero over the temperature range of calibration (400°–1300° C). Consideration of the mixing properties of olivine-orthopyroxene solid solutions places tight constraints on the standard state thermodynamic quantities describing Fe-Mg exchange reactions involving olivine, orthopyroxene, pyralspite garnets, aluminate spinels, ferrite spinels and biotite. These constraints are entirely consistent with the standard state properties for the phases-quartz,-quartz, orthoenstatite, clinoenstatite, protoenstatite, fayalite, ferrosilite and forsterite which were deduced by Berman (1988) from an independent analysis of phase equilibria and calorimetric data. In conjunction with these standard state properties, the solution model presented in this paper provides a means of evaluating an internally consistent set of Gibbs energies of mineral solid solutions in the system Mg₂SiO₄-Fe₂SiO₄-SiO₂ over the temperature range 0–1300° C and pressure interval 0.001–50 kbars. As a consequence of our analysis, we find that the excess Gibbs energies associated with mixing of Fe and Mg in (Fe, Mg)₂SiO₄ olivines, (Fe, Mg)₃Al₂Si₃O₁₂ garnets, (Fe, Mg)Al₂O₄ and (Fe, Mg)Fe₂O₄ spinels, and K(Mg, Fe)₃AlSi₃O₁₀(OH)₂ biotites may be satisfactory described, on a macroscopic basis, with symmetric regular solution type parameters having values of 4.86±0.12 (olivine), 3.85±0.09 (garnet), 1.96±0.13 (spinel), and 3.21±0.29 kcals/gfw (biotite). Applications of the proposed solution model demonstrate the sensitivity of petrologic modeling to activity-composition relations of olivine-orthopyroxene solutions. We explore the consequences of estimating the activity of silica in melts forming in the mantle and we develop a graphical geothermometer/geobarometer for metamorphic assemblages of olivine+orthopyroxene+quartz. Quantitative evaluation of these results suggests that accurate and realistic estimates of silica activity in melts derived from mantle source regions,P-T paths of metamorphism and other intensive variables of petrologic interest await further refinements involving the addition of trace elements (Al³⁺ and Fe³⁺) to the thermodynamic formulation for orthopyroxenes.     

Phosphate bonding configuration on ferrihydrite based on molecular orbital calculations and XANES fingerprinting

Sorption of phosphate by Fe(III)- and Al(III)-(hydr)oxide minerals regulates the mobility of this potential water pollutant in the environment. The objective of this research was to determine the molecular configuration of phosphate bound on ferrihydrite at pH 6 by interpreting P K-edge XANES spectra in terms of bonding mode. XANES and UV-visible absorption spectra for aqueous Fe(III)-PO₄ solutions (Fe/P molar ratio = 0-2.0) provided experimental trends for energies of P(3p)-O(2p) and Fe(3d)-O(2p) antibonding molecular orbitals. Molecular orbitals for Fe(III)-PO₄ or Al(III)-PO₄ complexes in idealized monodentate or bidentate bonding mode were generated by conceptual bonding arguments, and Extended-Hückel molecular orbital computations were used to understand and assign XANES spectral features to bound electronic states. The strong white line at the absorption edge in P K-edge XANES spectra for Fe-PO₄ or Al-PO₄ systems is attributable to an electronic transition from a P 1s atomic orbital into P(3p)-O(2p) or P(3p)-O(2p)-Al(3p) antibonding molecular orbitals, respectively. For Fe-PO₄ systems, a XANES peak at 2-5 eV below the edge was assigned to a P 1s electron transition into Fe(4p)-O(2p) antibonding molecular orbitals. Similarly, a shoulder on the low-energy side of the white line for variscite corresponds to a transition into Al(3p)-O(2p) orbitals. In monodentate-bonded phosphate, Fe-O bonding is optimized and P-O bonding is weakened, and the converse is true of bidentate-bonded phosphate. These differences explained an inverse correlation between energies of P(3p)-O(2p) and Fe(3d)-O(2p) antibonding molecular orbitals consistent with a monodentate-to-bidentate transition in aqueous Fe(III)-PO₄ solutions. The intensity of the XANES pre-edge feature in Fe(III)-bonded systems increased with increasing number of Fe(III)-O-P bonds. Based on the similarity of intensity and splitting of the pre-edge feature for phosphate sorbed on ferrihydrite at 750 mmol/kg at pH 6 and aqueous Fe-PO₄ solutions containing predominantly bidentate complexes, XANES results indicated that phosphate adsorbed on ferrihydrite was predominantly a bidentate-binuclear surface complex.     

Fe2+-Mg order-disorder in orthopyroxene: equilibrium fractionation between the octahedral sites and thermodynamic analysis

The equilibrium intracrystalline distribution coefficient, k _D ^*, of Fe^* (i.e. Fe²⁺ + Mn) and Mg between the M1 and M2 sites of three natural nearly binary Fe²⁺-Mg orthopyroxene crystals (Fs₁₄, Fs₁₅ and Fs₄₉) were determined by annealing experiments at several temperatures between 550 and 1000 °C and single crystal X-ray structure refinements. In addition, the X-ray data of an orthopyroxene crystal (Fs₂₃), which were collected earlier by Molin et al. (1991) between 700 and 1000 °C, were re-refined. The data were processed through two different refinement programs (SHELXL-93 and RFINE90) using both unit and individual weights and also both ionic and atomic scattering factors. The calculated site occupancies were found to agree within their estimated standard errors. However, the use of ionic scattering factors led to significantly better goodness of fit and agreement index, and smaller standard deviations of the site occupancies than those obtained from the use of atomic scattering factors. Furthermore, the weighted refinements yielded significantly smaller standard deviations of the site occupancies than the unweighted refinements even when the same set of reflections was used in the two procedures. The site occupancy data from this study were combined with selected published data to develop expressions of k _D ^* as a function of temperature and composition. Calculation of the excess configurational entropy, ΔS ^XS, suggests that orthopyroxene should be treated as a two parameter symmetric solution instead of as a “simple mixture”. The calculated ΔS ^XS values and the excess Gibbs free energy of mixing suggested by available cation exchange data lead to a slightly negative enthalpy of mixing in the orthopyroxene solid solution. Received: 25 August 1998 / Accepted: 10 March 1999     

Computer simulations of the interactions of the (0 1 2) and (0 0 1) surfaces of jarosite with Al, Cd, Cu and Zn

Jarosite is an important mineral on Earth, and possibly on Mars, where it controls the mobility of iron, sulfate and potentially toxic metals. Atomistic simulations have been used to study the incorporation of Al³⁺, and the M²⁺ impurities Cd, Cu and Zn, in the (0 1 2) and (0 0 1) surfaces of jarosite. The calculations show that the incorporation of Al on an Fe site is favorable on all surfaces in which terminal Fe ions are exposed, and especially on the (0 0 1) [Fe₃(OH)₃]⁶⁺ surface. Incorporation of Cd, Cu or Zn on a K site balanced by a K vacancy is predicted to stabilize the surfaces, but calculated endothermic solution energies and the high degree of distortion of the surfaces following incorporation suggest that these substitutions will be limited. The calculations also suggest that incorporation of Cd, Cu and Zn on an Fe site balanced by an OH vacancy, or by coupled substitution on both K and Fe sites, is unfavorable, although this might be compensated for by growth of a new layer of jarosite or goethite, as predicted for bulk jarosite. The results of the simulations show that surface structure will exert an influence on uptake of impurities in the order Cu > Cd > Zn, with the most favorable surfaces for incorporation being (0 1 2) [KFe(OH)₄]⁰ and (0 0 1) [Fe₃(OH)₃]⁶⁺.     

Symplectitic cordierite-orthopyroxene-garnet assemblages as products of contact metamorphism of pre-existing basement granulites in the Vredefort Structure,South Africa,and their relations to pseudotachylite

Former idioblastic garnet crystals of an Archaean granulite with Mg/(Mg+Fe+Mn)-ratios (= M) near 0.40 were converted to symplectitic pseudomorphs consisting mainly of cordierite (M=0.61), orthopyroxene (M=0.40), and a relic garnet with a new composition (M=0.18), during a static metamorphism which is related in time to the Vredefort event. On the basis of experimental data for the continuous reaction garnet+quartz=orthopyroxene+cordierite the conditions of metamorphism were near 5 kb, 700 °C. Orthopyroxenes crystallized initially as minute, myrmekitic grains with metastable excess Al contents up to 13 mol% Al₂O₃ and, through grain growth and Al exsolution, transformed into dense aggregates of coarser crystals with equilibrated compositions near 3 mol% Al₂O₃. In the absence of free silica hercynitic spinel appears as an additional phase that coexists with cordierite, orthopyroxene, and garnet of more magnesian compositions than in the case of silica saturation.Pseudotachylite veins crosscutting the hornfelsed granulite consist of pyroxene of variable Al contents, two feldspars, opaques, and quartz, and are finely recrystallized. There is textural and mineral chemical evidence from both the Al-contents of the pyroxenes and the Mg/Fe distribution among the phases adjacent to the pseudotachylite that these veins were emplaced late during the static metamorphism, which was probably caused by a local, deep-seated magma diapir. On this basis, the meteorite impact hypothesis for the origin of the Vredefort Structure is considered less likely than a mechanism of internal origin which can be made responsible for the close link in time and space between static metamorphism and the high strain rate deformation that produced the pseudotachylite.     

Oxygen isotope evidence for sorption of molecular oxygen to pyrite surface sites and incorporation into sulfate in oxidation experiments

Experiments were conducted to investigate (i) the rate of O-isotope exchange between SO₄ and water molecules at low pH and surface temperatures typical for conditions of acid mine drainage (AMD) and (ii) the O- and S-isotope composition of sulfates produced by pyrite oxidation under closed and open conditions (limited and free access of atmospheric O₂) to identify the O source/s in sulfide oxidation (water or atmospheric molecular O₂) and to better understand the pyrite oxidation pathway. An O-isotope exchange between SO₄ and water was observed over a pH range of 0–2 only at 50 °C, whereas no exchange occurred at lower temperatures over a period of 8 a. The calculated half-time of the exchange rate for 50 °C (pH = 0 and 1) is in good agreement with former experimental data for higher and lower temperatures and excludes the possibility of isotope exchange for typical AMD conditions (T  25 °C, pH  3) for decades.Pyrite oxidation experiments revealed two dependencies of the O-isotope composition of dissolved sulfates: O-isotope values decreased with longer duration of experiments and increasing grain size of pyrite. Both changes are interpreted as evidence for chemisorption of molecular O₂ to pyrite surface sites. The sorption of molecular O₂ is important at initial oxidation stages and more abundant in finer grained pyrite fractions and leads to its incorporation in the produced SO₄. The calculated bulk contribution of atmospheric O₂ in the dissolved SO₄ reached up to 50% during initial oxidation stages (first 5 days, pH 2, fine-grained pyrite fraction) and decreased to less than 20% after about 100 days. Based on the direct incorporation of molecular O₂ in the early-formed sulfates, chemisorption and electron transfer of molecular O₂ on S sites of the pyrite surface are proposed, in addition to chemisorption on Fe sites. After about 10 days, the O of all newly-formed sulfates originates only from water, indicating direct interaction of hydroxyls from water with S at the anodic S pyrite surface site. Then, the role of molecular O₂ is as proposed in previous studies: acting as electron acceptor only at the cathodic Fe pyrite surface site for oxidation of Fe(II) to Fe(III).     

Temperature-jump induced cation exchange kinetics in (Co<Subscript>0.1</Subscript>Mg<Subscript>0.9</Subscript>)<Subscript>2</Subscript>SiO<Subscript>4</Subscript> olivine: an in situ optical spectroscopic study

In the olivine crystal structure, cations are distributed over two inequivalent octahedral sites, M1 and M2. Kinetics of cation exchange between the two octahedral sites in (Co_0.1Mg_0.9)₂SiO₄ single crystal have been studied in the temperature range from 600 to 800°C by monitoring the time evolution of the absorbance of Co²⁺ ions in M1 or M2 sites using optical spectroscopy after rapid temperature jumps. It was found from such temperature-jump induced relaxation experiments that with increasing temperature the absorbance of Co²⁺ ions in the M1 site decreases while that in the M2 site increases. This indicates a tendency of Co²⁺ cations to populate the M2 site with increasing temperatures and vice versa. The experimental relaxation data can be modeled using a triple exponential equation based on theoretical analysis. Activation energies of 221 ± 4 and 213 ± 10 kJ/mol were derived from relaxation experiments on the M2 site and M1 site, respectively, for the cation exchange processes in (Co_0.1Mg_0.9)₂SiO₄ olivine. Implications for cation diffusion at low temperatures are discussed.     

The influence of crystal field stabilization energy on Fe2+ partitioning in paragenetic minerals

In order to explore possible quantitative relations between crystal field stabilization energy, CFSE, and partitioning behaviour of the 3d ⁶-configured Fe²⁺ ion, a suite of 29 paragenetic rock-forming minerals from 12 high-grade metamorphic rock samples of the Ukrainian shield, including the parageneses garnet/orthopyroxene/clinopyroxene (2x), orthopyroxene/clinopyroxene, garnet/clinopyroxene, garnet/orthopyroxene/biotite, garnet/biotite, garnet/cordierite, garnet/cordierite/biotite, garnet/orthopyroxene/clinopyroxene/Ca-amphibole, Ca-amphibole/biotite (retrograde), was studied by electron microprobe analysis to obtain the respective K _D Fe^{2+ (Ph1/Ph2)} values and by polarized single crystal electronic absorption spectroscopy to evaluate the respective CFSE_Fe2+ values. Other than in the case of Cr³⁺, a clear quantitative relation between K _D ^(Ph1/Ph2) and the ΔCFSE^(Ph1/Ph2) was only observed when geometrical factors, mainly the volume of crystallographic sites and ionic radii of ions competing in the partitioning process, are similar in the respective two paragenetic phases to within 15–20%. In such cases, the ΔCFSE_Fe2+ contribution to K _D ^(Ph1/Ph2) amounts to 0.1 to 0.2 log K _D per 100 cm⁻¹ΔCFSE. The conclusion is that ΔCFSE_Fe2+ plays only a secondary role after geometrical factors, in the partitioning behaviour of Fe²⁺. The reason for this is seen in the facts that, compared to the 3d ³-configured Cr³⁺ ion, CFSE of the 3d ⁶-configured Fe²⁺ amounts only to 20–25%, and that the former ion enters only octahedral sites with similar geometrical properties in the paragenetic mineral phases. Received: 17 November 1998 / Accepted: 28 June 1999     

The bonding of ferrous iron to sulfur and oxygen in tetrahedral coordination: A comparative study using SCF Xα scattered wave molecular orbital calculations

Molecular quantum mechanical calculations have been performed on high-spin ferrous iron tetrahedrally coordinated to sulfur and oxygen, respectively. The molecular orbital energies obtained from the calculations are compared with experimental optical and X-ray emission spectra. Good agreement was found between calculated and experimental spectral transition energies for the optical absorption spectra of Fe²⁺ in sphalerite, of Fe²⁺ in FeAl₂O₄, staurolite and (Zn, Fe)O, and for the FeKβ X-ray emission spectra of FeCr₂O₄. This both clarified interpretation of the spectra and established the validity of the calculations. Distinct differences occur in the molecular orbital structures of the sulfide and oxide clusters. In the sulfide, the crystal field type (mainly Fe 3d) molecular orbitals lie within the nonbonding (mainly S 3p) orbitais in energy, whereas in the oxide, they lie well above the 02p nonbonding orbitals. This also results in a wider valence band in the oxide than in the sulfide. The crystal field type (Fe 3d) molecular orbitais have more ligand character in the sulfide than the oxide and the chalcophilic properties of iron are partly attributed to this observation.     

Oxygen-exchange pathways in aluminum polyoxocations

Using molecular dynamics simulations and electronic structure methods, we postulate a mechanism to explain the complicated reactivity trends that are observed for oxygen isotope exchange reactions between sites in aluminum polyoxocations of the ε-Keggin type and bulk solution. Experimentally, the molecules have four nonequivalent oxygens that differ considerably in reactivity both within a molecule, and between molecules in the series: Al₁₃, GaAl₁₂, and GeAl₁₂ [MO₄Al₁₂(OH)₂₄(H₂O)₁₂ⁿ⁺(aq); with M = Al(III) for Al₁₃, n = 7; M = Ga(III) for GaAl₁₂, n = 7; M = Ge(IV) for GeAl₁₂, n = 8]. We find that a partly dissociated, metastable intermediate molecule of expanded volume is necessary for exchange of both sets of μ₂-OH and that the steady-state concentration of this intermediate reflects the bond strengths between the central metal and the μ₄-O. Thus the central metal exerts extraordinary control over reactions at hydroxyl bridges, although these are three bonds away.This mechanism not only explains the reactivity trends for oxygen isotope exchange in μ₂-OH and η-OH₂ sites in the ε-Keggin aluminum molecules, but also explains the observation that the reactivities of minerals tend to reflect the presence of highly coordinated oxygens, such as the μ₄-O in boehmite, α-, and γ-Al₂O₃ and their Fe(III) analogs. The partial dissociation of these highly coordinated oxygens, coupled with simultaneous activation and displacement of neighboring metal centers, may be a fundamental process by which metals atoms undergo ligand exchanges at mineral surfaces.     

Crystal orbital contributions to the pyrrhotite valence band with XPS evidence for weak Fe–Fe π bond formation

 Synchrotron excited X-ray photoelectron spectra (SXPS) of hexagonal pyrrhotite reveal three distinct Fe 3d-derived photopeaks within its outer valence band. The t _2gα band (majority spin) is centered at about 2.5 eV, the e _g α band at about 1.0 eV and the t _2gβ (minority spin) contribution at about 0.25 eV. From these data the ligand field splitting energy is 1.5 (±0.2) eV and the majority spin pairing energy is 2.25 (±0.2) eV. These are the first such XPS measurements for this mineral. S 3p-derived bonding and non-bonding bands are identified, with the former centred at about 6.5 eV and the latter near 4.5 eV. The XPS results are remarkably consistent with SCF-Xα scattered wave molecular orbital calculations. Although the calculations and the collected spectra are consistent, they differ from a recent interpretation of the pyrrhotite valence band. An explanation for the discrepant results is provided. Auger resonant enhancement of Fe 3d photopeaks at 60 eV photon energy results in the t _2gα emission (at 2.5 eV) being strongly enhanced and broader than the t _2gβ emission (0.25 eV). The explanation of these observations requires the presence of weak Fe–Fe π and π^* crystal (molecular) orbitals located near 2.5 eV, and separated by no more than about 0.5 eV. The π-bonded crystal orbitals are derived from weak mixing of adjacent Fe t _2g atomic orbitals along the c crystallographic axis. Received: 15 June 2000 / Accepted: 11 June 2001