skip to main content


Title: An ab-initio study on the thermodynamics of disulfide, sulfide, and bisulfide incorporation into apatite and the development of a more comprehensive temperature, pressure, pH, and composition-dependent model for ionic substitution in minerals
The mineral apatite, Ca10(PO4)6(F,OH,Cl)2, incorporates sulfur (S) during crystallization from S-bearing hydrothermal fluids and silicate melts. Our previous studies of natural and experimental apatite demonstrate that the oxidation state of S in apatite varies systematically as a function of oxygen fugacity (fO2). The S oxidation states –1 and –2 were quantitatively identified in apatite crystallized from reduced, S-bearing hydrothermal fluids and silicate melts by using sulfur K-edge X‑ray absorption near-edge structure spectroscopy (S-XANES) where S6+/ΣS in apatite increases from ~0 at FMQ-1 to ~1 at FMQ+2, where FMQ refers to the fayalite-magnetite-quartz fO2 buffer. In this study, we employ quantum-mechanical calculations to investigate the atomistic structure and energetics of S(-I) and S(-II) incorporated into apatite and elucidate incorporation mechanisms. One S(-I) species (disulfide, S22−) and two S(-II) species (bisulfide, HS−, and sulfide, S2−) are investigated as possible forms of reduced S species in apatite. In configuration models for the simulation, these reduced S species are positioned along the c-axis channel, originally occupied by the column anions F, Cl, and OH in the end-member apatites. In the lowest-energy configurations of S-incorporated apatite, disulfide prefers to be positioned halfway between the mirror planes at z = 1/4 and 3/4. In contrast, the energy-optimized bisulfide is located slightly away from the mirror planes by ~0.04 fractional units in the c direction. The energetic stability of these reduced S species as a function of position along the c-axis can be explained by the geometric and electrostatic constraints of the Ca and O planes that constitute the c-axis channel. The thermodynamics of incorporation of disulfide and bisulfide into apatite are evaluated by using solid-state reaction equations where the apatite host and a solid S-bearing source phase (pyrite and Na2S2(s) for disulfide; troilite and Na2S(s) for sulfide) are the reactants, and the S-incorporated apatite and an anion sink phase are the products. The Gibbs free energy (ΔG) is lower for incorporation with Na-bearing phases than with Fe-bearing phases, which is attributed to the higher energetic stability of the iron sulfide minerals as a source phase for S than the sodium sulfide phases. The thermodynamics of incorporation of reduced S are also evaluated by using reaction equations involving dissolved disulfide and sulfide species [HnS2(aq)(2–n) and HnS(aq)(2–n); n = 0, 1, and 2] as a source phase. The ΔG of S-incorporation increases for fluorapatite and chlorapatite and decreases for hydroxylapatite as these species are protonated (i.e., as n changes from 0 to 2). These thermodynamic results demonstrate that the presence of reduced S in apatite is primarily controlled by the chemistry of magmatic and hydrothermal systems where apatite forms (e.g., an abundance of Fe; solution pH). Ultimately, our methodology developed for evaluating the thermodynamics of S incorporation in apatite as a function of temperature, pH, and composition is highly applicable to predicting the trace and volatile element incorporation in minerals in a variety of geological systems. In addition to solid-solid and solid-liquid equilibria treated here at different temperatures and pH, the methodology can be easily extended also to different pressure conditions by just performing the quantum-mechanical calculations at elevated pressures.  more » « less
Award ID(s):
1924142
PAR ID:
10357798
Author(s) / Creator(s):
Date Published:
Journal Name:
The American mineralogist
ISSN:
1945-3027
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The mineral apatite, Ca10(PO4)6(F,OH,Cl)2, incorporates sulfur (S) during crystallization from S-bearing hydrothermal fluids and silicate melts. Our previous studies of natural and experimental apatite demonstrate that the oxidation state of S in apatite varies systematically as a function of oxygen fugacity (fO2). The S oxidation states –1 and –2 were quantitatively identified in apatite crystallized from reduced, S-bearing hydrothermal fluids and silicate melts by using sulfur K-edge X-ray absorption near-edge structure spectroscopy (S-XANES) where S 6+/ΣS in apatite increases from ~0 at FMQ-1 to ~1 at FMQ+2, where FMQ refers to the fayalite-magnetite-quartz fO2 buffer. In this study, we employ quantum-mechanical calculations to investigate the atomistic structure and energetics of S(-I) and S(-II) incorporated into apatite and elucidate incorporation mechanisms.

    One S(-I) species (disulfide, S22−) and two S(-II) species (bisulfide, HS−, and sulfide, S2−) are investigated as possible forms of reduced S species in apatite. In configuration models for the simulation, these reduced S species are positioned along the c-axis channel, originally occupied by the column anions F, Cl, and OH in the end-member apatites. In the lowest-energy configurations of S-incorporated apatite, disulfide prefers to be positioned halfway between the mirror planes at z = 1/4 and 3/4. In contrast, the energy-optimized bisulfide is located slightly away from the mirror planes by ~0.04 fractional units in the c direction. The energetic stability of these reduced S species as a function of position along the c-axis can be explained by the geometric and electrostatic constraints of the Ca and O planes that constitute the c-axis channel.

    The thermodynamics of incorporation of disulfide and bisulfide into apatite is evaluated by using solid-state reaction equations where the apatite host and a solid S-bearing source phase (pyrite and Na2S2(s) for disulfide; troilite and Na2S(s) for sulfide) are the reactants, and the S-incorporated apatite and an anion sink phase are the products. The Gibbs free energy (ΔG) is lower for incorporation with Na-bearing phases than with Fe-bearing phases, which is attributed to the higher energetic stability of the iron sulfide minerals as a source phase for S than the sodium sulfide phases. The thermodynamics of incorporation of reduced S is also evaluated by using reaction equations involving dissolved disulfide and sulfide species [HnS(aq)(2−n) and HnS(aq)(2−n); n = 0, 1, and 2] as a source phase. The ΔG of S-incorporation increases for fluorapatite and chlorapatite, and decreases for hydroxylapatite, as these species are protonated (i.e., as n changes from 0 to 2). These thermodynamic results demonstrate that the presence of reduced S in apatite is primarily controlled by the chemistry of magmatic and hydrothermal systems where apatite forms (e.g., an abundance of Fe; solution pH). Ultimately, our methodology developed for evaluating the thermodynamics of S incorporation in apatite as a function of temperature, pH, and composition is highly applicable to predicting the trace and volatile element incorporation in minerals in a variety of geological systems. In addition to solid-solid and solid-liquid equilibria treated here at different temperatures and pH, the methodology can be easily extended to different pressure conditions by just performing the quantum-mechanical calculations at elevated pressures.

     
    more » « less
  2. Most known porphyry Cu±Au deposits are associated with moderately oxidized and sulfur-rich, calc-alkaline to mildly alkalic arc-related magmas in the Phanerozoic. In contrast, sodium-enriched tonalite–trondhjemite–granodiorite–diorite (TTG) magmas predominant in the Archean are hypothesized to be unoxidized and sulfur-poor, which together preclude porphyry Cu deposit formation. Here, we test this hypothesis by interrogating the causative magmas for the ∼2·7 Ga TTG-related Côté Gold, St-Jude, and Clifford porphyry-type Cu±Au deposit settings in the Neoarchean southern Abitibi subprovince. New and previously published geochronological results constrain the age of emplacement of the causative magmas at ∼2·74 Ga, ∼2·70 Ga, and∼2·69 Ga, respectively. The dioritic and trondhjemitic magmas associated with Côté Gold and St-Jude evolved along a plagioclase-dominated fractionation trend, in contrast to amphibole-dominated fractionation for tonalitic magma at Clifford. Analyses of zircon grains from the Côté Gold, St-Jude, and Clifford igneous rocks yielded εHf(t)±SD values of 4·5±0·3, 4·2±0·6, and 4·3±0·4, and δ18O±SD values of 5·40±0·11  , 3·91±0·13  , and 4·83±0·12  , respectively. These isotopic signatures indicate that, although these magmas are mantle-sourced with minimal crustal contamination, for the St- Jude and Clifford settings the magmas or their sources may have undergone variable alteration by heated seawater or meteoric fluids. Primary barometric minerals (i.e. zircon, amphibole, apatite, and magnetite–ilmenite) that survived variable alteration and metamorphism (up to greenschist facies) were used for estimating fO2 of the causative magmas. Estimation of magmatic fO2 values, reported relative to the fayalite–magnetite–quartz buffer as  FMQ, using zircon geochemistry indicates that the fO2 values of the St-Jude, Côté Gold, and Clifford magmas increase from  FMQ –0·3±0·6 to  FMQ +0·8±0·4 and to  FMQ +1·2±0·4, respectively. In contrast, amphibole chemistry yielded systematically higher fO2 values of  FMQ +1·6±0·3 and  FMQ +2·6±0·1 for Côté Gold and Clifford, respectively, which are consistent with previous studies that indicate that amphibole may overestimate the fO2 of intrusive rocks by up to 1 log unit. Micro X-ray absorption near edge structure (μ-XANES) spectrometric determination of sulfur (i.e. S6+/ S) in primary apatite yielded ≥ FMQ−0·3 and FMQ+1·4–1·8 for St-Jude and Clifford, respectively. The magnetite–ilmenite mineral pairs from the Clifford tonalite yielded  FMQ +3·3±1·3 at equilibrium temperatures of 634±21 ◦C, recording the redox state of the late stage of magma crystallization. Electron probe microanalyses revealed that apatite grains from Clifford are enriched in S (up to 0·1 wt%) relative to those of Côté Gold and St-Jude (below the detection limit), which is attributed to either relatively oxidized or sulfur-rich features of the Clifford tonalite. We interpret these results to indicate that the deposits at Côté Gold and Clifford formed from mildly (∼ FMQ +0·8±0·4) to moderately (∼ FMQ +1·5) oxidized magmas where voluminous early sulfide saturation was probably limited, whereas the St-Jude deposit represents a rare case whereby the ingress of externally derived hydrothermal fluids facilitated metal fertility in a relatively reduced magma chamber (∼ FMQ +0). Furthermore, we conclude that variable modes of formation for these deposits and, in addition, the apparent rarity of porphyry-type Cu–Au deposits in the Archean may be attributed to either local restriction of favorable metallogenic conditions, and/or preservation, or an exploration bias. 
    more » « less
  3. Abstract Recycling of oxidized sulfur from subducting slabs to the mantle wedge provides simultaneous explanations for the elevated oxygen fugacity (fO2) in subduction zones, their high hydrothermal and magmatic sulfur outputs, and the enriched sulfur isotopic signatures (i.e., δ34S > 0‰) of these outputs. However, a quantitative understanding of the abundance and speciation of sulfur in slab fluids consistent with high pressure experiments is lacking. Here we analyze published experimental data for anhydrite solubility in H2O‐NaCl solutions to calibrate a high‐pressure aqueous speciation model of sulfur within the framework of the deep earth water model. We characterize aqueous complexes, required to account for the high experimental anhydrite solubilities. We then use this framework to predict the speciation and solubility of sulfur in chemically complex fluids in equilibrium with model subducting mafic and ultramafic lithologies, from 2 to 3 GPa and 400 to 800°C at logfO2from FMQ‐2 to FMQ+4. We show that sulfate complexes of calcium and sodium markedly enhance the stability of sulfate in moderately oxidized fluids in equilibrium with pyrite atfO2conditions of FMQ+1 to +2, causing large sulfur isotope fractionations up to 10‰ in the fluid relative to the slab. Such fluids could impart oxidized, sulfur‐rich and high δ34S signatures to the mantle wedge that are ultimately transferred to arc magmas, without the need to invoke34S‐rich subducted lithologies.

     
    more » « less
  4. The large range in oxidation states of sulfur (-II to +VI) provides it with a large oxidation potential in rocks, even at relatively low concentrations. Most importantly, the transition from sulfide to sulfate species in rocks and silicate melts occurs in the same approximate fO2 region (for a given temperature) as the transition from ferrous to ferric iron, and reduced S species can coexist with oxidized Fe and vice versa. The result is a large potential for reactions involving sulfur to oxidize or reduce Fe in silicate minerals, since Fe only occurs in two oxidation states (+II and +III). In order for sulfur to be released during slab dehydration, sulfur in sulfide must be converted into an easily dissolved species, such as SO42− or H2S, through either oxidation or reduction. We propose that oxidation of sulfur in sulfide follows the generalized reaction: 8Fe3+SiaOb(OH)c +S2− = 8Fe2+SidOe +SO42− +(H2O)f (1) In this type of reaction, sulfur participates in the dehydration of greenschist- or blueschist-facies hydrous silicates during transition to the eclogite facies: ferric Fe in Fe-bearing silicates (chlorite, amphibole, epidote) is reduced to ferrous Fe in anhydrous ferromagnesian silicates (pyroxene, garnet). At the same time, the reaction consumes sulfide by oxidation of S2− to produce SO42−, which is readily dissolved in the fluid produced during dehydration. Additionally, a similar redox reaction could oxidize sulfur by reducing ferric Fe in oxides. It is important to note that one mole of S has the same redox potential as 8 moles of Fe. The molar ratio of 8 moles of Fe per 1 mole of S translates to a mass ratio of approximately 14; therefore, small concentrations of sulfur can have a large impact on reduction/oxidation of the silicate assemblage. Our observations show that sulfide minerals that can be identified as primary or related to the peak metamorphic stage are rare in eclogites and restricted to inclusions in garnet, consistent with reaction (1). Thermodynamic modeling is currently underway to assess the influence of sulfur on the phase equilibria of silicate phases during high pressure metamorphism. 
    more » « less
  5. Abstract

    Secondary minerals in martian nakhlites provide a powerful tool for investigating the nature, composition, and duration of aqueous activity in the martian crust. Northwest Africa (NWA) 998 crystallized early from the nakhlite magmatic source and has evidence of minimal signatures of the late hydrothermal alteration event that altered the nakhlites. Using FIB‐TEM techniques to study a cumulus apatite grain in NWA 998, we report the first evidence of a submicron‐scale vein consisting of fluorapatite and an SiO2‐rich phase. Fluorapatite grew epitaxially on the walls of an opened cleavage plane of host F‐bearing chlorapatite and the SiO2‐rich phase filled the center of the vein. The presence of nanoporosity and nanometer‐scale amorphous material and the sharp interface between the vein and the host apatite indicate the vein represents a coupled dissolution–reprecipitation process that generated apatite of a different composition that was more stable with the fluid. Using experimental data and diffusion coefficients of Cl in apatite from the literature, we conclude that the vein was caused by a low temperature (~300°C), slightly acidic, F‐, Si‐rich, aqueous fluid that acted as a closed system. Based on the characteristics of the vein (formation by rapid injection of fluid) and the fluid (composition, temperature, pH), and the lack of terrestrial weathering products in our SEM and TEM images, we infer that the vein is pre‐terrestrial in origin. Our observations support the hypothesis that the heat source triggering a hydrothermal system was a low‐shock velocity impact and rule out a magmatic origin. Finally, the vein could have formed from a late‐stage fluid different from that reported in other nakhlites, but formation during the same magmatic event by, for example, a less evolved fluid might also be plausible.

     
    more » « less