The structural behavior of ferric and ferrous iron in aluminosilicate glass near meta-aluminosilicate joins

Iron-57 resonant absorption Mössbauer spectroscopy was used to describe the redox relations and structural roles of Fe³⁺ and Fe²⁺ in meta-aluminosilicate glasses. Melts were formed at 1500 °C in equilibrium with air and quenched to glass in liquid H₂O with quenching rates exceeding 200 °C/s. The aluminosilicate compositions were NaAlSi₂O₆, Ca_0.5AlSi₂O₆, and Mg_0.5AlSi₂O₆. Iron oxide was added in the form of Fe₂O₃, NaFeO₂, CaFe₂O₄, and MgFe₂O₄ with total iron oxide content in the range ∼0.9 to ∼5.6 mol% as Fe₂O₃. The Mössbauer spectra, which were deconvoluted by assuming Gaussian distributions of the hyperfine field, are consistent with one absorption doublet of Fe²⁺ and one of Fe³⁺. From the area ratios of the Fe²⁺ and Fe³⁺ absorption doublets, with corrections for differences in recoil-fractions of Fe³⁺ and Fe²⁺, the Fe³⁺/ΣFe is positively correlated with increasing total iron content and with decreasing ionization potential of the alkali and alkaline earth cation. There is a distribution of hyperfine parameters from the Mössbauer spectra of these glasses. The maximum in the isomer shift distribution function of Fe³⁺, δ_Fe³⁺, ranges from about 0.25 to 0.49 mm/s (at 298 K relative to Fe metal) with the quadrupole splitting maximum, Δ_Fe³⁺, ranging from ∼1.2 to ∼1.6 mm/s. Both δ_Fe³⁺ and δ_Fe²⁺ are negatively correlated with total iron oxide content and Fe³⁺/ΣFe. The dominant oxygen coordination number Fe³⁺ changes from 4 to 6 with decreasing Fe³⁺/ΣFe. The distortion of the Fe³⁺-O polyhedra of the quenched melts (glasses) decreases as the Fe³⁺/ΣFe increases. These polyhedra do, however, coexist with lesser proportions of polyhedra with different oxygen coordination numbers. The δ_Fe²⁺ and Δ_Fe²⁺ distribution maxima at 298 K range from ∼0.95 to 1.15 mm/s and 1.9 to 2.0 mm/s, respectively, and decrease with increasing Fe³⁺/ΣFe. We suggest that these hyperfine parameter values for the most part are more consistent with Fe²⁺ in a range of coordination states from 4- to 6-fold. The lower δ_Fe²⁺-values for the most oxidized melts are consistent with a larger proportion of Fe²⁺ in 4-fold coordination compared with more reduced glasses and melts.     

Aluminum speciation, vibrational entropy and short-range order in calcium aluminosilicate glasses

The heat capacity and vibrational entropy of a calcium aluminate and three peraluminous calcium aluminosilicate glasses have been determined from 2 to 300 K by heat-pulse relaxation calorimetry. Together with previous adiabatic data for six other glasses in the system CaO-Al₂O₃-SiO₂, these results have been used to determine partial molar heat capacities and entropies for five species namely, SiO₂, CaO and three different sorts of Al₂O₃ in which Al is 4-, 5- and 6-fold coordinated by oxygen. Given the determining role of oxygen coordination on low-temperature heat capacity, the composition independent entropies found for SiO₂ and CaO indicate that short-range order around Si and Ca is not sensitive to aluminum speciation up to the highest fraction of 25% observed for ^VAl by NMR spectroscopy. Because of the higher room-temperature vibrational entropy of ^IVAl₂O₃ (72.8 J/mol K) compared to ^VAl₂O₃ (48.5 J/mol K), temperature-induced changes from ^IVAl to ^VAl give rise to a small negative contribution of the order of 1 J/mol K to the partial molar configurational heat capacity of Al₂O₃ in melts. Near 0 K, pure SiO₂ glass distinguishes itself by the importance of the calorimetric boson peak. On a g atom basis, the maximum of this peak varies with the composition of calcium aluminosilicate glasses by a factor of about 2. It does not show smooth variations, however, either as a function of SiO₂ content, at constant CaO/Al₂O₃ ratio, or as a function of Al₂O₃ content, at constant SiO₂ content.     

The oxidation state of sulfur in synthetic and natural glasses determined by X-ray absorption spectroscopy

Sulfur K-edge X-ray absorption near edge structure (XANES) spectra were recorded for experimental glasses of various compositions prepared at different oxygen fugacities (fO₂) in one-atmosphere gas-mixing experiments at 1400 °C. This sample preparation method only results in measurable S concentrations under either relatively reduced (log fO₂ < −9) or oxidised (log fO₂ > −2) conditions. The XANES spectra of the reduced samples are characterised by an absorption edge crest at 2476.4 eV, typical of S²⁻. In addition, spectra of Fe-bearing compositions exhibit a pronounced absorption edge shoulder. Spectra for all the Fe-free samples are essentially identical, as are the spectra for the Fe-bearing compositions, despite significant compositional variability within each group. The presence of a sulfide phase, such as might exsolve on cooling, can be inferred from a pre-edge feature at 2470.5 eV.The XANES spectra of the oxidised samples are characterised by an intense transition at 2482.1 eV, typical of the sulfate anion SO₄²⁻. Sulfite (SO₃²⁻) has negligible solubility in silicate melts at low pressures. The previous identification of sulfite species in natural glass samples is attributed to an artefact of the analysis (photoreduction of S⁶⁺). S⁴⁺ does, however, occur unambiguously with S⁶⁺ in Fe-free and Fe-poor compositions prepared in equilibrium with CaSO₄ at 4-16 kbar, and when buffered with Re/ReO₂ at 10 kbar. Solubility of S⁴⁺ thus requires partial pressures of SO₂ considerably in excess of 1 bar. A number of experiments were undertaken in an attempt to access intermediate fO₂s more applicable to terrestrial volcanism. Although these were largely unsuccessful, S²⁻ and S⁶⁺ were found to coexist in some samples that were not in equilibrium with the imposed fO₂.The XANES spectra of natural olivine-hosted melt inclusions and submarine glasses representative of basalts at, or close to, sulfide saturation show mainly dissolved S²⁻, but with minor sulfate, and additionally a peak at 2469.5 eV, which, although presumably due to immiscible sulfide, is 1 eV lower than that typical of FeS. These sulfate and sulfide-related peaks disappear with homogenisation of the inclusions by heating to 1200 °C followed by rapid quenching, suggesting that both these features are a result of cooling under natural conditions. The presence of small amounts of sulfate in otherwise reduced basaltic magmas may be explained by the electron exchange reaction: S²⁻ + 8Fe³⁺ = S⁶⁺ + 8Fe²⁺, which is expected to proceed strongly to the right with decreasing temperature. This reaction would explain why S²⁻ and S⁶⁺ are frequently found together despite the very limited fO₂ range over which they are thermodynamically predicted to coexist. The S XANES spectra of water-rich, highly oxidised, basaltic inclusions hosted in olivine from Etna and Stromboli confirm that nearly all S is dissolved as sulfate, explaining their relatively high S contents.     

Sulfur K-edge XANES analysis of natural and synthetic basaltic glasses: Implications for S speciation and S content as function of oxygen fugacity

XANES analyses at the sulfur K-edge were used to determine the oxidation state of S species in natural and synthetic basaltic glasses and to constrain the fO₂ conditions for the transition from sulfide (S²⁻) to sulfate (S⁶⁺) in silicate melts. XANES spectra of basaltic samples from the Galapagos spreading center, the Juan de Fuca ridge and the Lau Basin showed a dominant broad peak at 2476.8 eV, similar to the spectra obtained from synthetic sulfide-saturated basalts and pyrrhotite. An additional sharp peak at 2469.8 eV, similar to that of crystalline sulfides, was present in synthetic glasses quenched from hydrous melts but absent in anhydrous glasses and may indicate differences in sulfide species with hydration or presence of minute sulfide inclusions exsolved during quenching. The XANES spectra of a basalt from the 1991 eruption of Mount Pinatubo, Philippines, and absarokitic basalts from the Cascades Range, Oregon, USA, showed a sharp peak at 2482.8 eV, characteristic of synthetic sulfate-saturated basaltic glasses and crystalline sulfate-bearing minerals such as hauyne. Basaltic samples from the Lamont Seamount, the early submarine phase of Kilauea volcano and the Loihi Seamount showed unequivocal evidence of the coexistence of S²⁻ and S⁶⁺ species, emphasizing the relevance of S⁶⁺ to these systems. XANES spectra of basaltic glasses synthesized in internally-heated pressure vessels and equilibrated at fO₂ ranging from FMQ − 1.4 to FMQ + 2.7 showed systematic changes in the features related to S²⁻ and S⁶⁺ with changes in fO₂. No significant features related to sulfite (S⁴⁺) species were observed. These results were used to construct a function that allows estimates of S⁶⁺/ΣS from XANES data. Comparison of S⁶⁺/ΣS data obtained by S Kα shifts measured with electron probe microanalysis (EPMA), S⁶⁺/ΣS obtained from XANES spectra, and theoretical considerations show that data obtained from EPMA measurements underestimate S⁶⁺/ΣS in samples that are sulfate-dominated (most likely because of photo-reduction effects during analysis) whereas S⁶⁺/ΣS from XANES provide a close match to the expected theoretical values. The XANES-derived relationship for S⁶⁺/ΣS as a function of fO₂ indicates that the transition from S²⁻ to S⁶⁻ with increasing fO₂ occurs over a narrower interval than what is predicted by the EPMA-derived relationship. The implications for natural systems is that small variation of fO₂ above FMQ + 1 will have a large effect on S behavior in basaltic systems, in particular regarding the amount of S that can be transported by basaltic melts before sulfide saturation can occur.     

Structure of H2O-saturated peralkaline aluminosilicate melt and coexisting aluminosilicate-saturated aqueous fluid determined in-situ to 800 °C and ∼800 MPa

The structure of H₂O-saturated silicate melts and of silicate-saturated aqueous solutions, as well as that of supercritical silicate-rich aqueous liquids, has been characterized in-situ while the sample was at high temperature (to 800 °C) and pressure (up to 796 MPa). Structural information was obtained with confocal microRaman and with FTIR spectroscopy. Two Al-bearing glasses compositionally along the join Na₂O•4SiO₂-Na₂O•4(NaAl)O₂-H₂O (5 and 10 mol% Al₂O₃, denoted NA5 and NA10) were used as starting materials. Fluids and melts were examined along pressure-temperature trajectories of isochores of H₂O at nominal densities (from PVT properties of pure H₂O) of 0.85 g/cm³ (NA10 experiments) and 0.86 g/cm³ (NA5 experiments) with the aluminosilicate + H₂O sample contained in an externally-heated, Ir-gasketed hydrothermal diamond anvil cell.Molecular H₂O (H₂O°) and OH groups that form bonds with cations exist in all three phases. The OH/H₂O° ratio is positively correlated with temperature and pressure (and, therefore, fugacity of H₂O, f_H2O) with (OH/H₂O°)^melt > (OH/H₂O°)^fluid at all pressures and temperatures. Structural units of Q³, Q², Q¹, and Q⁰ type occur together in fluids, in melts, and, when outside the two-phase melt + fluid boundary, in single-phase liquids. The abundance of Q⁰ and Q¹ increases and Q² and Q³ decrease with f_H2O. Therefore, the NBO/T (nonbridging oxygen per tetrahedrally coordination cations), of melt is a positive function of f_H2O. The NBO/T of silicate in coexisting aqueous fluid, although greater than in melt, is less sensitive to f_H2O.The melt structural data are used to describe relationships between activity of H₂O and melting phase relations of silicate systems at high pressure and temperature. The data were also combined with available partial molar configurational heat capacity of Qⁿ-species in melts to illustrate how these quantities can be employed to estimate relationships between heat capacity of melts and their H₂O content.     

Spectroscopic analysis (FTIR, Raman) of water in mafic and intermediate glasses and glass inclusions

Micro-Raman spectroscopy, even though a very promising technique, is not still routinely applied to analyse H₂O in silicate glasses. The accuracy of Raman water determinations critically depends on the capability to predict and take into account both the matrix effects (bulk glass composition) and the analytical conditions on band intensities. On the other hand, micro-Fourier transform infrared spectroscopy is commonly used to measure the hydrous absorbing species (e.g., hydroxyl OH⁻ and molecular H₂O) in natural glasses, but requires critical assumptions for the study of crystal-hosted glasses. Here, we quantify for the first time the matrix effect of Raman external calibration procedures for the quantification of the total H₂O content (H₂O_T = OH⁻ + H₂O_m) in natural silicate glasses. The procedures are based on the calibration of either the absolute (external calibration) or scaled (parameterisation) intensity of the 3550 cm⁻¹ band. A total of 67 mafic (basanite, basalt) and intermediate (andesite) glasses hosted in olivines, having between 0.2 and 4.8 wt% of H₂O, was analysed. Our new dataset demonstrates, for given water content, the height (intensity) of Raman H₂O_T band depends on glass density, reflectance and water environment. Hence this matrix effect must be considered in the quantification of H₂O by Raman spectroscopy irrespective of the procedure, whereas the parameterisation mainly helps to predict and verify the self-consistency of the Raman results. In addition, to validate the capability of the micro-Raman to accurately determine the H₂O content of multicomponent aluminosilicate glasses, a subset of 23 glasses was analysed by both micro-Raman and micro-FTIR spectroscopy using the band at 3550 cm⁻¹. We provide new FTIR absorptivity coefficients (ε₃₅₅₀) for basalt (62.80 ± 0.8 L mol⁻¹ cm⁻¹) and basanite (43.96 ± 0.6 L mol⁻¹ cm⁻¹). These values, together with an exhaustive review of literature data, confirm the non-linear decline of the FTIR absorptivity coefficient (ε₃₅₅₀) as the glass depolymerisation increases. We demonstrate the good agreement between micro-FTIR and micro-Raman determination of H₂O in silicate glasses when the matrix effects are properly considered.     

Sulfur speciation and network structural changes in sodium silicate glasses: Constraints from NMR and Raman spectroscopy

Information about the state of sulfur in silicate melts and glasses is important in both earth sciences and materials sciences. Because of its variety of valence states from S²⁻ (sulfide) to S⁶⁺ (sulfate), the speciation of sulfur dissolved in silicate melts and glasses is expected to be highly dependent on the oxygen fugacity. To place new constraint on this issue, we have synthesized sulfur-bearing sodium silicate glasses (quenched melts) from starting materials containing sulfur of different valence states (Na₂SO₄, Na₂SO₃, Na₂S₂O₃ and native S) using an internally heated gas pressure vessel, and have applied electron-induced SKα X-ray fluorescence, micro-Raman and NMR spectroscopic techniques to probe their structure. The wavelength shift of SKα X-rays revealed that the differences in the valence state of sulfur in the starting compounds are largely retained in the synthesized sulfur-bearing glasses, with a small reduction for more oxidized samples. The ²⁹Si MAS NMR spectra of all the glasses contain no peaks attributable to the SiO_4-nS_n (with n > 0) linkages. The Raman spectra are consistent with the coexistence of sodium sulfate (Na₂SO₄) species and one or more types of more reduced sulfur species containing S-S linkages in all the sulfur-bearing silicate glasses, with the former dominant in glasses produced from Na₂SO₄-doped starting materials, and the latter more abundant in more reduced glasses. The ²⁹Si MAS NMR and Raman spectra also revealed changes in the silicate network structure of the sulfur-bearing glasses, which can be interpreted in terms of changes in the chemical composition and sulfur speciation.     

Iron-bearing silicate melts: Relations between pressure and redox equilibria

Mossbauer spectroscopy has been used to determine the redox equilibria of iron and structure of quenched melts on the composition join Na₂Si₂O₅-Fe₂O₃ to 40 kbar pressure at 1400° C. The Fe³⁺/ΣFe decreases with increasing pressure. The ferric iron appears to undergo a gradual coordination transformation from a network-former at 1 bar to a network-modifier at higher (≧10 kbar) pressure. Ferrous iron is a network-modifier in all quenched melts. Reduction of Fe³⁺ to Fe²⁺ and coordination transformation of remaining Fe³⁺ result in depolymerization of the silicate melts (the ratio of nonbridging oxygens per tetrahedral cations, NBO/T, increases). It is suggested that this pressure-induced depolymerization of iron-bearing silicate liquids results in increasing NBO/T of the liquidus minerals. Furthermore, this depolymerization results in a more rapid pressure-induced decrease in viscosity and activation energy of viscous flow of iron-bearing silicate melts than would be expected for iron-free silicate melts with similar NBO/T.     

Influence of glass polymerisation and oxidation on micro-Raman water analysis in alumino-silicate glasses

The development of an accurate analytical procedure for determination of dissolved water in complex alumino-silicate glasses via micro-Raman analysis requires the assessment of the spectra topology dependence on glass composition. We report here a detailed study of the respective influence of bulk composition, iron oxidation state and total water content on the absolute and relative intensities of the main Raman bands related to glass network vibrations (LF: ∼490 cm⁻¹; HF: ∼960 cm⁻¹) and total water stretching (H₂O_T: ∼3550 cm⁻¹) in natural glasses. The evolution of spectra topology was examined in (i) 33 anhydrous glasses produced by the re-melting of natural rock samples, which span a very large range of polymerisation degree (NBO/T from 0.00 to 1.16), (ii) 2 sets of synthetic anhydrous basaltic glasses with variable iron oxidation state (Fe³⁺/Fe_T from 0.05 to 0.87), and (iii) 6 sets of natural hydrous glasses (C_H2OT from 0.4 to 7.0 wt%) with NBO/T varying from 0.01 to 0.76.In the explored domain of water concentration, external calibration procedure based on the H₂O_T band height is matrix-independent but its accuracy relies on precise control of the focusing depth and beam energy on the sample. Matrix-dependence strongly affects the internal calibrations based on H₂O_T height scaled to that of LF or HF bands but its effect decreases from acid (low NBO/T, SM) to basic (high NBO/T, SM) glasses. Structural parameters such as NBO/T (non-bridging oxygen per tetrahedron) and SM (sum of structural modifiers) describe the matrix-dependence better than simple compositional parameters (e.g. SiO₂, Na₂O + K₂O). Iron oxidation state has only a minor influence on band topology in basalts and is thus not expected to significantly affect the Raman determinations of water in mafic (e.g. low SiO₂, iron-rich) glasses. Modelling the evolution of the relative band height with polymerisation degree allows us to propose a general equation to predict the dissolved water content in natural glasses:

     

Effects of the degree of polymerization on the structure of sodium silicate and aluminosilicate glasses and melts: An O NMR study

Revealing the atomic structure and disorder in oxide glasses, including sodium silicates and aluminosilicates, with varying degrees of polymerization, is a challenging problem in high-temperature geochemistry as well as glass science. Here, we report ¹⁷O MAS and 3QMAS NMR spectra for binary sodium silicate and ternary sodium aluminosilicate glasses with varying degrees of polymerization (Na₂O/SiO₂ ratio and Na₂O/Al₂O₃ ratio), revealing in detail the extent of disorder (network connectivity and topological disorder) and variations of NMR parameters with the glass composition. In binary sodium silicate glasses [Na₂O-k(SiO₂)], the fraction of non-bridging oxygens (NBOs, Na-O-Si) increases with the Na₂O/SiO₂ ratio (k), as predicted from the composition. The ¹⁷O isotropic chemical shifts (¹⁷O δ_iso) for both bridging oxygen (BO) and NBO increase by about 10-15 ppm with the SiO₂ content (for k = 1-3). The quadrupolar coupling products of BOs and NBOs also increase with the SiO₂ content. These trends suggest that both NBOs and BOs strongly interact with Na; therefore, the Na distributions around BOs and NBOs are likely to be relatively homogenous for the glass compositions studied here, placing some qualitative limits on the extent of segregation of alkali channels from silica-enriched regions as suggested by modified random-network models. The peak width (in the isotropic dimension) and thus bond angle and length distributions of Si-O-Si and Na-O-Si increase with the SiO₂ content, indicating an increase in the topological disorder with the degree of polymerization. In the ternary aluminosilicate glasses [Na₂O]_x[Al₂O₃]_1−xSiO₂, the NBO fraction decreases while the Al-O-Si and Al-O-Al fractions apparently increase with increasing Al₂O₃ content. The variation of oxygen cluster populations suggests that deviation from “Al avoidance” is more apparent near the charge-balanced join (Na/Al = 1). The Si-O-Si fraction, which is closely related to the activity coefficient of silica, would decrease with increasing Al₂O₃ content at a constant mole fraction of SiO₂. Therefore, the activity of silica may decrease from depolymerized binary silicates to fully polymerized sodium aluminosilicate glasses at a constant mole fraction of SiO₂.     

The effect of water activity on the oxidation and structural state of Fe in a ferro-basaltic melt

Experimental investigations have been performed at T = 1200°C, P = 200 MPa and fH₂ corresponding to H₂O-MnO-Mn₃O₄ and H₂O-QFM redox buffers to study the effect of H₂O activity on the oxidation and structural state of Fe in an iron-rich basaltic melt. The analysis of Mössbauer and Fe K-edge X-ray absorption nearedge structure (XANES) spectra of the quenched hydrous ferrobasaltic glasses shows that the Fe³⁺/ΣFe ratio of the glass is directly related to aH₂O in a H₂-buffered system and, consequently, to the prevailing oxygen fugacity (through the reaction of water dissociation H₂O ↔ H₂ + 1/2 O₂). However, water as a chemical component of the silicate melt has an indistinguishable effect on the redox state of iron at studied conditions. The experimentally obtained relationship between fO₂ and Fe³⁺/Fe²⁺ in the hydrous ferrobasaltic melt can be adequately predicted in the investigated range by the existing empiric and thermodynamic models. The ratio of ferric and ferrous Fe is proportional to the oxygen fugacity to the power of ∼0.25 which agrees with the theoretical value from the stoichiometry of the Fe redox reaction (FeO + ¼ O₂ = FeO_1.5). The mean centre shifts for Fe²⁺ and Fe³⁺ absorption doublets in Mössbauer spectra show little change with increasing Fe³⁺/ΣFe, suggesting no significant change in the type of iron coordination. Similarly, XANES preedge spectra indicate a mixed (C3h, Td, and Oh, i.e., 5-, 4-, and sixfold) coordination of Fe in hydrous basaltic glasses.     

Effect of water on the heat capacity of polymerized aluminosilicate glasses and melts

The effect of water on heat capacity has been determined for four series of hydrated synthetic aluminosilicate glasses and supercooled liquids close to albite, phonolite, trachyte, and leucogranite compositions. Heat capacities were measured at atmospheric pressure by differential scanning calorimetry for water contents between 0 and 4.9 wt % from 300 K to about 100 K above the glass transition temperature (T_g). The partial molar heat capacity of water in polymerized aluminosilicate glasses, which can be considered as independent of composition, is (J/mol K). In liquids containing at least 1 wt % H₂O, the partial molar heat capacity of water is about 85 J/mol K. From speciation data, the effects of water as hydroxyl groups and as molecular water have tentatively been estimated, with partial molar heat capacities of 153 ± 18 and 41 ± 14 J/mol K, respectively. In all cases, water strongly increases the configurational heat capacity at T_g and exerts a marked depressing effect on T_g, in close agreement with the results of viscosity experiments on the same series of glasses. Consistent with the Adam and Gibbs theory of relaxation processes, the departure of the viscosity of hydrous melts from Arrhenian variations correlates with the magnitude of configurational heat capacities.     

The experimental calibration of the iron isotope fractionation factor between pyrrhotite and peralkaline rhyolitic melt

A first experimental study was conducted to determine the equilibrium iron isotope fractionation between pyrrhotite and silicate melt at magmatic conditions. Experiments were performed in an internally heated gas pressure vessel at 500 MPa and temperatures between 840 and 1000 °C for 120-168 h. Three different types of experiments were conducted and after phase separation the iron isotope composition of the run products was measured by MC-ICP-MS. (i) Kinetic experiments using ⁵⁷Fe-enriched glass and natural pyrrhotite revealed that a close approach to equilibrium is attained already after 48 h. (ii) Isotope exchange experiments—using mixtures of hydrous peralkaline rhyolitic glass powder (∼4 wt% H₂O) and natural pyrrhotites (Fe_1 − xS) as starting materials— and (iii) crystallisation experiments, in which pyrrhotite was formed by reaction between elemental sulphur and rhyolitic melt, consistently showed that pyrrhotite preferentially incorporates light iron. No temperature dependence of the fractionation factor was found between 840 and 1000 °C, within experimental and analytical precision. An average fractionation factor of Δ ⁵⁶Fe/⁵⁴Fe_{pyrrhotite-melt} = −0. 35 ± 0.04‰ (2SE, n = 13) was determined for this temperature range. Predictions of Fe isotope fractionation between FeS and ferric iron-dominated silicate minerals are consistent with our experimental results, indicating that the marked contrast in both ligand and redox state of iron control the isotope fractionation between pyrrhotite and silicate melt. Consequently, the fractionation factor determined in this study is representative for the specific Fe²⁺/ΣFe ratio of our peralkaline rhyolitic melt of 0.38 ± 0.02. At higher Fe²⁺/ΣFe ratios a smaller fractionation factor is expected. Further investigation on Fe isotope fractionation between other mineral phases and silicate melts is needed, but the presented experimental results already suggest that even at high temperatures resolvable variations in the Fe isotope composition can be generated by equilibrium isotope fractionation in natural magmatic systems.     

CO2 solubility in Martian basalts and Martian atmospheric evolution

To understand possible volcanogenic fluxes of CO₂ to the Martian atmosphere, we investigated experimentally carbonate solubility in a synthetic melt based on the Adirondack-class Humphrey basalt at 1-2.5 GPa and 1400-1625 °C. Starting materials included both oxidized and reduced compositions, allowing a test of the effect of iron oxidation state on CO₂ solubility. CO₂ contents in experimental glasses were determined using Fourier transform infrared spectroscopy (FTIR) and Fe³⁺/Fe^T was measured by Mössbauer spectroscopy. The CO₂ contents of glasses show no dependence on Fe³⁺/Fe^T and range from 0.34 to 2.12 wt.%. For Humphrey basalt, analysis of glasses with gravimetrically-determined CO₂ contents allowed calibration of an integrated molar absorptivity of 81,500 ± 1500 L mol⁻¹ cm⁻² for the integrated area under the carbonate doublet at 1430 and 1520 cm⁻¹. The experimentally determined CO₂ solubilities allow calibration of the thermodynamic parameters governing dissolution of CO₂ vapor as carbonate in silicate melt, K_II, (Stolper and Holloway, 1988) as follows: , ΔV⁰ = 20.85 ± 0.91 cm³ mol⁻¹, and ΔH⁰ = −17.96 ± 10.2 kJ mol⁻¹. This relation, combined with the known thermodynamics of graphite oxidation, facilitates calculation of the CO₂ dissolved in magmas derived from graphite-saturated Martian basalt source regions as a function of P, T, and f_O2. For the source region for Humphrey, constrained by phase equilibria to be near 1350 °C and 1.2 GPa, the resulting CO₂ contents are 51 ppm at the iron-wüstite buffer (IW), and 510 ppm at one order of magnitude above IW (IW + 1). However, solubilities are expected to be greater for depolymerized partial melts similar to primitive shergottite Yamato 980459 (Y 980459). This, combined with hotter source temperatures (1540 °C and 1.2 GPa) could allow hot plume-like magmas similar to Y 980459 to dissolve 240 ppm CO₂ at IW and 0.24 wt.% of CO₂ at IW + 1. For expected magmatic fluxes over the last 4.5 Ga of Martian history, magmas similar to Humphrey would only produce 0.03 and 0.26 bars from sources at IW and IW + 1, respectively. On the other hand, more primitive magmas like Y 980459 could plausibly produce 0.12 and 1.2 bars at IW and IW + 1, respectively. Thus, if typical Martian volcanic activity was reduced and the melting conditions cool, then degassing of CO₂ to the atmosphere may not be sufficient to create greenhouse conditions required by observations of liquid surface water. However, if a significant fraction of Martian magmas derive from hot and primitive sources, as may have been true during the formation of Tharsis in the late Noachian, that are also slightly oxidized (IW + 1.2), then significant contribution of volcanogenic CO₂ to an early Martian greenhouse is plausible.     

Temperature independent thermal expansivities of calcium aluminosilicate melts between 1150 and 1973 K in the system anorthite-wollastonite-gehlenite (An-Wo-Geh): A density model

The thermal expansivities of 10 compositions from within the anorthite-wollastonite-gehlenite (An-Wo-Geh) compatibility triangle have been investigated using a combination of calorimetry and dilatometry on the glassy and liquid samples. The volumes at room temperature were derived from densities measured using the Archimedean buoyancy method. For each sample, density was measured at 298 K using glass that had a cooling-heating history of 10-10 K min⁻¹. The thermal expansion coefficient of the glass from 298 K to the glass transition interval was measured by a dilatometer and the heat capacity was measured using a differential scanning calorimeter from 298 to 1135 K. The thermal expansion coefficient and the heat flow were determined at a heating rate of 10 K min⁻¹ on glasses which were previously cooled at 10 K min⁻¹. Supercooled liquid density, molar volume and molar thermal expansivities were indirectly determined by combining differential scanning calorimetric and dilatometric measurements assuming that the kinetics of enthalpy and shear relaxation are equivalent. The data obtained on supercooled liquids were compared to high-temperature predictions from the models of (Lange, R.A., Carmichael, I.S.E., 1987. Densities of Na₂O-K₂O-CaO-MgO-FeO-Fe₂O₃-Al₂O₃-TiO₂-SiO₂ liquids: New measurements and derived partial molar properties. Geochim. Cosmochim. Acta51, 2931-2946; Courtial, P., Dingwell, D.B., 1995. Nonlinear composition dependence of molar volume of melts in the CaO-Al₂O₃-SiO₂ system. Geochim. Cosmochim. Acta59 (18), 3685-3695; Lange, R.A., 1997. A revised model for the density and thermal expansivity of K₂O-Na₂O-CaO-MgO-Al₂O₃-SiO₂ liquids from 700 to 1900 K: extension to crustal magmatic temperatures. Contrib. Mineral. Petrol.130, 1-11). The best linear fit combines the supercooled liquid data presented in this study and the high temperature data calculated using the Courtial and Dingwell (1995) model. This dilatometric/calorimetric method of determining supercooled liquid molar thermal expansivity greatly increases the temperature range accessible for thermal expansion. It represents a substantial increase in precision and understanding of the thermodynamics of calcium aluminosilicate melts. This enhanced precision demonstrates clearly the temperature independence of the melt expansions in the An-Wo-Geh system. This contrasts strongly with observations for neighboring system such as anorthite-diopside and raises the question of the compositional/structural origins of temperature dependence of thermal expansivity in multicomponent silicate melts.     

The effect of titanium and phosphorus on ferric/ferrous ratio in silicate melts: an experimental study

The effect of TiO₂ and P₂O₅ on the ferric/ferrous ratio in silicate melts was investigated in model silicate melts at air conditions in the temperature range 1,400–1,550 °C at 1-atm total pressure. The base composition of the anorthite–diopside eutectic composition was modified with 10 wt % Fe₂O₃ and variable amounts of TiO₂ (up to 30 wt %) or P₂O₅ (up to 20 wt %). Some compositions also contained higher SiO₂ concentrations to compare the role of SiO₂, TiO₂, and P₂O₅ on the Fe³⁺/Fe²⁺ ratio. The ferric/ferrous ratio in experimental glasses was analyzed using a wet chemical technique with colorimetric detection of ferrous iron. It is shown that at constant temperature, an increase in SiO₂, TiO₂, and P₂O₅ content results in a decrease in the ferric/ferrous ratio. The effects of TiO₂ and SiO₂ on the Fe³⁺/Fe²⁺ ratio was found to be almost identical. In contrast, adding P₂O₅ was found to decrease ferric/ferrous ratio much more effectively than adding silica. The results were compared with the predictions from the published empirical equations forecasting Fe³⁺/Fe²⁺ ratio. It was demonstrated that the effects of TiO₂ are minor but that the effects of P₂O₅ should be included in models to better describe ferric/ferrous ratio in phosphorus-bearing silicate melts. Based on our observations, the determination of the prevailing fO₂ in magmas from the Fe³⁺/Fe²⁺ ratio in natural glasses using empirical equations published so far is discussed critically.     

Molecular H2O in armenite, BaCa2Al6Si9O30·2H2O, and epididymite, Na2Be2Si6O15·H2O: Heat capacity, entropy and local-bonding behavior of confined H2O in microporous silicates

Armenite, ideal formula BaCa₂Al₆Si₉O₃₀·2H₂O, and its dehydrated analog BaCa₂Al₆Si₉O₃₀ and epididymite, ideal formula Na₂Be₂Si₆O₁₅·H₂O, and its dehydrated analog Na₂Be₂Si₆O₁₅ were studied by low-temperature relaxation calorimetry between 5 and 300 K to determine the heat capacity, C_p, behavior of their confined H₂O. Differential thermal analysis and thermogravimetry measurements, FTIR spectroscopy, electron microprobe analysis and powder Rietveld refinements were undertaken to characterize the phases and the local environment around the H₂O molecule.The determined structural formula for armenite is Ba_0.88(0.01)Ca_1.99(0.02)Na_0.04(0.01)Al_5.89(0.03)Si_9.12(0.02)O₃₀·2H₂O and for epididymite Na_1.88(0.03)K_0.05(0.004)Na_0.01(0.004)Be_2.02(0.008)Si_6.00(0.01)O₁₅·H₂O. The infrared (IR) spectra give information on the nature of the H₂O molecules in the natural phases via their H₂O stretching and bending vibrations, which in the case of epididymite only could be assigned. The powder X-ray diffraction data show that armenite and its dehydrated analog have similar structures, whereas in the case of epididymite there are structural differences between the natural and dehydrated phases. This is also reflected in the lattice IR mode behavior, as observed for the natural phases and the H₂O-free phases. The standard entropy at 298 K for armenite is S° = 795.7 ± 6.2 J/mol K and its dehydrated analog is S° = 737.0 ± 6.2 J/mol K. For epididymite S° = 425.7 ± 4.1 J/mol K was obtained and its dehydrated analog has S° = 372.5 ± 5.0 J/mol K. The heat capacity and entropy of dehydration at 298 K are Δ = 3.4 J/mol K and ΔS^rxn = 319.1 J/mol K and Δ = −14.3 J/mol K and ΔS^rxn = 135.7 J/mol K for armenite and epididymite, respectively. The H₂O molecules in both phases appear to be ordered. They are held in place via an ion-dipole interaction between the H₂O molecule and a Ca cation in the case of armenite and a Na cation in epididymite and through hydrogen-bonding between the H₂O molecule and oxygen atoms of the respective silicate frameworks. Of the three different H₂O phases ice, liquid water and steam, the C_p behavior of confined H₂O in both armenite and epididymite is most similar to that of ice, but there are differences between the two silicates and from the C_p behavior of ice. Hydrogen-bonding behavior and its relation to the entropy of confined H₂O at 298 K is analyzed for various microporous silicates.The entropy of confined H₂O at 298 K in various silicates increases approximately linearly with increasing average wavenumber of the OH-stretching vibrations. The interpretation is that decreased hydrogen-bonding strength between a H₂O molecule and the silicate framework, as well as weak ion-dipole interactions, results in increased entropy of H₂O. This results in increased amplitudes of external H₂O vibrations, especially translations of the molecule, and they contribute strongly to the entropy of confined H₂O at T < 298 K.     

Structural control on bulk melt properties: Single and double quantum Si NMR spectroscopy on alkali-silicate glasses

The structure of 21 binary potassium, rubidium and cesium silicate glasses (in the range 15-50 mol% alkali oxide) was analyzed by ²⁹Si single quantum and double quantum MAS NMR spectroscopy. Their glass transition temperatures (T_g) were measured by calorimetry. The chemical shifts and the relative abundance of Qⁿ species correlate with the cationic field strength (Z/r) of the network modifier. A correlation is observed between T_g and the inverse of the entropy of mixing of the different Qⁿ species, which is explained in the framework of the Adam-Gibbs relaxation theory. At high alkali content, up to 44% of the SiO₄ tetrahedra are part of three-membered rings. At a given alkali content, the abundance of these rings increases with increasing cation size. The abundance of three-membered rings in K-silicate melts correlates with a temperature and a non-linear composition dependence of the heat capacity. It is also a possible cause for the anomalous volumetric behavior of potassium silicate glasses.     

Structure of hydrous aluminosilicate glasses along the diopside-anorthite join: A comprehensive one- and two-dimensional H and Al NMR study

We have taken a systematic approach utilizing advanced solid-state NMR techniques to gain new insights into the controversial issue concerning the dissolution mechanisms of water in aluminosilicate melts (glasses). A series of quenched anhydrous and hydrous (∼2 wt% H₂O) glass samples along the diopside (Di, CaMgSi₂O₆)—anorthite (An, CaAl₂Si₂O₈) join with varying An components (0, 20, 38, 60, 80, and 100 mol %) have been studied. A variety of NMR techniques, including one-dimensional (1D) ¹H and ²⁷Al MAS NMR, and ²⁷Al → ¹H cross-polarization (CP) MAS NMR, as well as two-dimensional (2D) ¹H double-quantum (DQ) MAS NMR, ²⁷Al triple-quantum (3Q) MAS NMR, and ²⁷Al → ¹H heteronuclear correlation NMR (HETCOR) and 3QMAS/HETCOR NMR, have been applied. These data revealed the presence of SiOH, free OH ((Ca,Mg)OH) and AlOH species in the hydrous glasses, with the last mostly interconnected with Si and residing in the more polymerized parts of the structure. Thus, there are no fundamental differences in water dissolution mechanisms for Al-free and Al-bearing silicate melts (glasses), both involving two competing processes: the formation of SiOH/AlOH that is accompanied by the depolymerization of the network structure, and the formation of free OH that has an opposite effect. The latter is more important for depolymerized compositions corresponding to mafic and ultramafic magmas.Aluminum is dominantly present in four coordination (Al^IV), but a small amount of five-coordinate Al (Al^V) is also observed in all the anhydrous and hydrous glasses. Furthermore, six-coordinate Al (Al^VI) is also present in most of the hydrous glasses. As Al of higher coordinations are favored by high pressure, Al^VIOH and Al^VOH may become major water species at higher pressures corresponding to those of the Earth’s mantle.     

Kinetics of biotite dissolution and Fe behavior under low O2 conditions and their implications for Precambrian weathering

Biotite dissolution experiments were carried out to better understand the dissolution kinetics and Fe behavior under low O₂ conditions, and to give an insight into the Precambrian weathering. Mineral dissolution with a continuous flow-through reactor was employed at 25 °C for up to 65 days varying partial pressure of atmospheric oxygen (PO₂), pH (6.86 and 3.01) and Fe content in mineral (1.06 and 0.11 mol of Fe per O₁₀(OH,F)₂ for biotite and phlogopite, respectively) independently for the examination of their effects on biotite dissolution. Low PO₂ conditions were achieved in a newly developed glove box (PO₂ ? 6 × 10⁻⁴ atm; referred to as anoxic conditions), which was compared to the present, ambient air conditions (0.2 atm of PO₂; oxic conditions). The biotite dissolution rate was slightly faster under anoxic conditions at pH 6.86 while it was not affected by PO₂ at pH 3.01. There was no direct effect of Fe content on dissolution rate at pH 6.86 while there was a small difference in dissolution rate between biotite and phlogopite at pH 3.01. The 1.5 order-of-magnitude faster release rate of Fe under anoxic conditions for biotite dissolution at pH 6.86 resulted from the difference in ratio of Fe³⁺ precipitates remaining in the reactor to Fe dissolved (about 60% and 100% under anoxic and oxic conditions, respectively), which is caused mainly by the difference in PO₂. The results infer that the Fe²⁺ and Fe³⁺ contents in the Paleoproterozoic paleosols, fossil weathering profiles, are reflected by atmospheric oxygen levels at the time of weathering.