Mineralogic and sulfur isotopic effects accompanying oxidation of pyrite in millimolar solutions of hydrogen peroxide at temperatures from 4 to 150 °C

Oxidation of pyrite by hydrogen peroxide (H₂O₂) at millimolar levels has been studied from 4 to 150 °C in order to evaluate isotopic effects potentially associated with radiolytic oxidation of pyrite. Gaseous, aqueous, and solid phases were collected and measured following sealed-tube experiments that lasted from 1 to 14 days. The dominant gaseous product was molecular oxygen. No volatile sulfur species were recovered from any experiment. Sulfate was the only aqueous sulfur species detected in solution, with sulfite and thiosulfate below the detection limits. X-ray diffraction patterns and images from scanning electron microscopy reveal solid residues composed primarily of hydrated ferric iron sulfates and sporadic ferric-ferrous iron sulfates. Hematite was detected only in solid residue produced during high temperature experiments. Elemental sulfur and/or polysulfides are inferred to be form on reacting pyrite surface based on extraction with organic solvents. Pyrite oxidation by H₂O₂ increases in rate with increasing H₂O₂concentration, pyrite surface area, and temperature. Rates measured in sealed-tube experiments at 25°C, for H₂O₂ concentration of 2 × 10⁻³ M are 8.8 × 10⁻⁹ M/m²/sec, which are higher than previous estimates. A combination of reactive oxygen species from H₂O₂ decomposition products and reactive iron species from pyrite dissolution is inferred to aggressively oxidize the receding pyrite surface. Competing oxidants with temperature-dependent oxidation efficiencies results in multiple reaction mechanisms for different temperatures and surface conditions. Sulfur isotope values of remaining pyrite were unchanged during the experiments, but showed distinct enrichment of ³⁴S in produced sulfate and depletion in elemental sulfur. The Δ_{sulfate-pyrite} and Δ_{elemental sulfur-pyrite} was +0.5 to +1.5‰ and was −0.2 to −1‰, respectively. Isotope data from high-temperature experiments indicate an additional ³⁴S-depleted sulfur fraction, with up to 4‰ depletion of ³⁴S, in the hematite. Sulfur isotope trends were not influenced by H₂O₂ concentration, temperature, or reaction time. Results of this study indicate that radiolytically produced oxidants, such as hydrogen peroxide and hydroxyl radicals, could efficiently oxidize pyrite in an otherwise oxygen-limited environment. Although H₂O₂ is generally regarded as being of minor geochemical significance on Earth, the H₂O₂ molecule plays a pivotal role in Martian atmospheric and soil chemistry. Additional experimental and field studies are needed to characterize sulfur and oxygen isotope systematics during radiolytical oxidation of metallic sulfides and elemental sulfur.     

Experimental measurement of the solubility of bismuth phases in water vapor from 220 °C to 300 °C: Implications for ore formation

Preliminary measurements were carried out of the solubility of the O_2-buffering assemblage bismuth + bismite (Bi₂O₃) in aqueous liquid–vapor and vapor-only systems at temperatures of 220, 250 and 300 °C. All experiments were carried out in Ti reaction vessels and were designed such that the Bi solids were contained in a silica tube that prevented contact with liquid water at any time during the experiment. Two blank (no Bi solids present) liquid–vapor experiments at 220 °C yielded Bi concentrations (±1σ) in the condensed liquid of 0.22 ± 0.02 mg/L, whereas the solubility measurements at this temperature yielded an average value of approximately 6 ± 9 mg/L, with replicate experiments ranging from 0.3 to 26 mg/L. Although the 6 mg/L value is associated with a considerable degree of uncertainty, the experiments do indicate transport of Bi through the vapor phase. Measured Bi concentrations in the condensed liquid at 250 °C were in the same range as those at 220 °C, whereas those at 300 °C were significantly lower (i.e., all below the blank value). Vapor-only experiments necessarily contained much smaller initial volumes of water, thereby making the results more susceptible to contamination. Single blank runs at 220 and 300 °C yielded Bi concentrations of 82 and 16 mg/L, respectively. Measured concentrations (±1σ) of Bi in the vapor-only solubility experiments at 220 °C were 235 ± 78 mg/L for an initial water volume of 0.5 mL, and at 300 °C were 56 ± 30 mg/L and 33 ± 21 for initial water volumes of 1 and 2 mL, respectively, suggesting strong preferential partitioning of Bi into the vapor. The results indicate a negative dependence of Bi solubility on temperature, but are inconclusive with respect to the dependence of Bi solubility on water density or fugacity.     

Arsenous acid ionisation in aqueous solutions from 25 to 300 °C

The ultraviolet spectra of dilute, aqueous arsenic (III)-containing solutions have been measured from 25 to 300 °C at the saturated vapour pressure. From these measurements, the equilibrium constant was obtained for the reaction

     

Dissolution of illite in saline-acidic solutions at 25 °C

The dissolution rate of illite, a common clay mineral in Australian soils, was studied in saline-acidic solutions under far from equilibrium conditions. The clay fraction of Na-saturated Silver Hill illite (K_1.38Na_0.05)(Al_2.87Mg_0.46Fe³⁺_0.39Fe²⁺_0.28Ti_0.07)[Si_7.02Al_0.98]O₂₀(OH)₄ was used for this study. The dissolution rates were measured using flow-through reactors at 25 ± 1 °C, solution pH range of 1.0-4.25 (H₂SO₄) and at two ionic strengths (0.01 and 0.25 M) maintained using NaCl solution. Illite dissolution rates were calculated from the steady state release rates of Al and Si. The dissolution stoichiometry was determined from Al/Si, K/Si, Mg/Si and Fe/Si ratios. The release rates of cations were highly incongruent during the initial stage of experiments, with a preferential release of Al and K over Si in majority of the experiments. An Al/Si ratio >1 was observed at pH 2 and 3 while a ratio close to the stoichiometric composition was observed at pH 1 and 4 at the higher ionic strength. A relatively higher K⁺ release rate was observed at I = 0.25 in 2-4 pH range than at I = 0.01, possibly due to ion exchange reaction between Na⁺ from the solution and K⁺ from interlayer sites of illite. The steady state release rates of K, Fe and Mg were higher than Si over the entire pH range investigated in the study. From the point of view of the dominant structural cations (Si and Al), stoichiometric dissolution of illite occurred at pH 1-4 in the higher ionic strength experiments and at pH ?3 for the lower ionic strength experiments. The experiment at pH 4.25 and at the lower ionic strength exhibited lower R_Al (dissolution rate calculated from steady state Al release) than R_Si (dissolution rate calculated from steady state Si release), possibly due to the adsorption of dissolved Al as the output solutions were undersaturated with respect to gibbsite. The dissolution of illite appears to proceed with the removal of interlayer K followed by the dissolution of octahedral cations (Fe, Mg and Al), the dissolution of Si is the limiting step in the illite dissolution process. A dissolution rate law showing the dependence of illite dissolution rate on proton concentration in the acid-sulfate solutions was derived from the steady state dissolution rates and can be used in predicting the impact of illite dissolution in saline acid-sulfate environments. The fractional reaction orders of 0.32 (I = 0.25) and 0.36 (I = 0.01) obtained in the study for illite dissolution are similar to the values reported for smectite. The dissolution rate of illite is mainly controlled by solution pH and no effect of ionic strength was observed on the dissolution rates.     

Antimonous acid protonation/deprotonation equilibria in hydrothermal solutions to 300 °C

The ultraviolet spectra of dilute aqueous solutions of antimony (III) have been measured from 25 to 300 °C at the saturated vapour pressure. From these measurements, equilibrium constants were obtained for the following reactions:

H₃SbO₃⁰ ? H⁺ + H₂SbO₃⁻     

An experimental study of the mechanism of the replacement of magnetite by pyrite up to 300 °C

We present the results of an experimental study into the sulfidation of magnetite to form pyrite/marcasite under hydrothermal conditions (90-300 °C, vapor saturated pressures), a process associated with gold deposition in a number of ore deposits. The formation of pyrite/marcasite was studied as a function of reaction time, temperature, pH, sulfide concentration, solid-weight-to-fluid-volume ratio, and geometric surface area of magnetite in polytetrafluoroethylene-lined autoclaves (PTFE) and a titanium and stainless steel flow-through cell. Marcasite was formed only at pH_21°C <4 and was the dominant Fe disulfide at pH_21°C 1.11, while pyrite predominated at pH_21°C >2 and formed even under basic conditions (up to pH_21°C 12-13). Marcasite formation was favored at higher temperatures. Fine-grained pyrrhotite formed in the initial stage of the reaction together with pyrite in some experiments with large surface area of magnetite (grain size <125 μm). This pyrrhotite eventually gave way to pyrite. The transformation rate of magnetite to Fe disulfide increased with decreasing pH (at 120 °C; pH_120°C 0.96-4.42), and that rate of the transformation increased from 120 to 190 °C.Scanning electron microscope (SEM) imaging revealed that micro-pores (0.1-5 μm scale) existed at the reaction front between the parent magnetite and the product pyrite, and that the pyrite and/or marcasite were euhedral at pH_21°C <4 and anhedral at higher pH. The newly formed pyrite was micro-porous (0.1-5 μm); this micro-porosity facilitates fluid transport to the reaction interface between magnetite and pyrite, thus promoting the replacement reaction. The pyrite precipitated onto the parent magnetite was polycrystalline and did not preserve the crystallographic orientation of the magnetite. The pyrite precipitation was also observed on the PTFE liner, which is consistent with pyrite crystallizing from solution. The mechanism of the reaction is that of a dissolution-reprecipitation reaction with the precipitation of pyrite being the rate-limiting step relative to magnetite dissolution under mildly acidic conditions (e.g., pH_155°C 4.42).The experimental results are in good agreement with sulfide phase assemblage and textures reported from sulfidized Banded Iron Formations: pyrite, marcasite and pyrrhotite have been found to exist or co-exist in different sulfidized Banded Iron Formations, and the microtextures show no evidence of sub-μm-scale pseudomorphism of magnetite by pyrite.     

Oxygen and hydrogen isotope fractionation in serpentine-water and talc-water systems from 250 to 450 °C, 50 MPa

Oxygen and hydrogen isotope fractionation factors in the talc-water and serpentine-water systems have been determined by laboratory experiment from 250 to 450 °C at 50 MPa using the partial exchange technique. Talc was synthesized from brucite + quartz, resulting in nearly 100% exchange during reaction at 350 and 450 °C. For serpentine, D-H exchange was much more rapid than ¹⁸O-¹⁶O exchange when natural chrysotile fibers were employed in the initial charge. In experiments with lizardite as the starting charge, recrystallization to chrysotile enhanced the rate of ¹⁸O-¹⁶O exchange with the coexisting aqueous phase. Oxygen isotope fractionation factors in both the talc-water and serpentine-water systems decrease with increasing temperature and can be described from 250 to 450 °C by the relationships: 1000 ln  = 11.70 × 10⁶/T² − 25.49 × 10³/T + 12.48 and 1000 ln  = 3.49 × 10⁶/T² − 9.48 where T is temperature in Kelvin. Over the same temperature interval at 50 MPa, talc-water D-H fractionation is only weakly dependent on temperature, similar to brucite and chlorite, and can be described by the equation: 1000 ln = 10.88 × 10⁶/T² − 41.52 × 10³/T + 5.61 where T is temperature in Kelvin. Our D-H serpentine-water fractionation factors calibrated by experiment decrease with temperature and form a consistent trend with fractionation factors derived from lower temperature field calibrations. By regression of these data, we have refined and extended the D-H fractionation curve from 25 to 450 °C, 50 MPa as follows: 1000 ln  = 3.436 × 10⁶/T² − 34.736 × 10³/T + 21.67 where T is temperature in Kelvin. These new data should improve the application of D-H and ¹⁸O-¹⁶O isotopes to constrain the temperature and origin of hydrothermal fluids responsible for serpentine formation in a variety of geologic settings.     

Kinetics of Fe(III) precipitation in aqueous solutions at pH 6.0-9.5 and 25 °C

The kinetics of Fe(III) precipitation in synthetic buffered waters have been investigated over the pH range 6.0-9.5 using a combination of visible spectrophotometry, ⁵⁵Fe radiometry combined with ion-pair solvent extraction of chelated iron and numerical modeling. The rate of precipitation, which is first order with respect to both dissolved and total inorganic ferric species, varies by nearly two orders of magnitude with a maximum rate constant of 16 ± 1.5 × 10⁶ M⁻¹ s⁻¹ at a pH of around 8.0. Our results support the existence of the dissolved neutral species, Fe(OH)₃⁰, and suggest that it is the dominant precursor in Fe(III) polymerization and subsequent precipitation at circumneutral pH. The intrinsic rate constant of precipitation of Fe(OH)₃⁰ was calculated to be allowing us to predict rates of Fe(III) precipitation in the pH range 6.0-9.5. The value of this rate constant, and the variation in the precipitation rate constant over the pH range considered, are consistent with a mechanism in which the kinetics of iron precipitation are controlled by rates of water exchange in dissolved iron hydrolysis species.     

Kinetics of calcite precipitation induced by ureolytic bacteria at 10 to 20°C in artificial groundwater

The kinetics of calcite precipitation induced in response to the hydrolysis of urea by Bacillus pasteurii at different temperatures in artificial groundwater (AGW) was investigated. The hydrolysis of urea by B. pasteurii exhibited a temperature dependence with first order rate constants of 0.91 d⁻¹ at 20°C, 0.18 d⁻¹ at 15°C, and 0.09 d⁻¹ at 10°C. At all temperatures, the pH of the AGW increased from 6.5 to 9.3 in less than 1 d. Dissolved Ca²⁺ concentrations decreased in an asymptotic fashion after 1 d at 20°C and 15°C, and 2 d at 10°C. The loss of Ca²⁺ from solution was accompanied by the development of solid phase precipitates that were identified as calcite by X-ray diffraction. The onset of calcite precipitation at each temperature occurred after similar amounts of urea were hydrolyzed, corresponding to 8.0 mM NH₄⁺. Specific rate constants for calcite precipitation and critical saturation state were derived from time course data following a second-order chemical affinity-based rate law. The calcite precipitation rate constants and critical saturation states varied by less than 10% between the temperatures with mean values of 0.16 ± 0.01 μmoles L⁻¹ d⁻¹ and 73 ±3, respectively. The highest calcite precipitation rates (ca. 0.8 mmol L⁻¹ d⁻¹) occurred near the point of critical saturation. While unique time course trajectories of dissolved Ca²⁺ concentrations and saturation state values were observed at different temperatures, calcite precipitation rates all followed the same asymptotic profile decreasing with saturation state regardless of temperature. This emphasizes the fundamental kinetic dependence of calcite precipitation on saturation state, which connects the otherwise dissimilar temporal patterns of calcite precipitation that evolved under the different temperature and biogeochemical regimes of the experiments.     

Heat capacities of aqueous solutions of sodium hydroxide and water ionization up to 300 °C at 10 MPa

A commercial (Setaram C80) calorimeter has been modified to measure the heat capacities of highly caustic solutions at temperatures up to 300 °C and pressures up to 20 MPa. The improvements have allowed more accurate determination of the isobaric volumetric heat capacities of chemically aggressive liquids at high temperatures. Test measurements with aqueous solutions of sodium chloride showed a reproducibility of about ±0.1%, with an accuracy of ∼0.3% or better, over the whole temperature range. Heat capacities of aqueous solutions of sodium hydroxide at concentrations from 0.5 to 8 mol/kg were measured at temperatures from 50 to 300 °C and a pressure of 10 MPa. Apparent molar isobaric heat capacities of NaOH(aq) were calculated using densities determined previously for the same solutions by vibrating-tube densimetry. Standard state (infinite dilution) partial molar isobaric heat capacities of NaOH(aq) were obtained by extrapolation using an extended Redlich-Meyer equation. Values of the standard heat capacity change for the ionization of water up to 300 °C were derived by combining the present results with the literature data for HCl(aq) and NaCl(aq).     

An experimental study of the solubility of molybdenum in H2O and KCl-H2O solutions from 500 °C to 800 °C, and 150 to 300 MPa

The solubility of molybdenum (Mo) was determined at temperatures from 500 °C to 800 °C and 150 to 300 MPa in KCl-H₂O and pure H₂O solutions in cold-seal experiments. The solutions were trapped as synthetic fluid inclusions in quartz at experimental conditions, and analyzed by laser ablation inductively coupled plasma mass spectrometry (LA ICPMS).Mo solubilities of 1.6 wt% in the case of KCl-bearing aqueous solutions and up to 0.8 wt% in pure H₂O were found. Mo solubility is temperature dependent, but not pressure dependent over the investigated range, and correlates positively with salinity (KCl concentration). Molar ratios of ∼1 for Mo/Cl and Mo/K are derived based on our data. In combination with results of synchrotron X-ray absorption spectroscopy of individual fluid inclusions, it is suggested that Mo-oxo-chloride complexes are present at high salinity (>20 wt% KCl) and ion pairs at moderate to low salinity (<11 wt% KCl) in KCl-H₂O aqueous solutions. Similarly, in the pure H₂O experiments molybdic acid is the dominant species in aqueous solution. The results of these hydrothermal Mo experiments fit with earlier studies conducted at lower temperatures and indicate that high Mo concentrations can be transported in aqueous solutions. Therefore, the Mo concentration in aqueous fluids seems not to be the limiting factor for ore formation, whereas precipitation processes and the availability of sulfur appear to be the main controlling factors in the formation of molybdenite (MoS₂).     

Dissolution kinetics of calcite at 50-70 °C: An atomic force microscopic study under near-equilibrium conditions

Direct measurements of calcite faces were performed using in situ atomic force microscopy (AFM) to reveal the dissolution processes as a function of solution saturation state and temperature. Time-sequential AFM images demonstrated that step velocities at constant temperature increased with increasing undersaturation. The anisotropy of obtuse and acute step velocities appeared to become more significant as solutions approached equilibrium and temperature increased. At saturation state Ω > 0.02, a curvilinear boundary was formed at the intersection of two acute steps and the initially rhombohedral etch pit exhibited a nearly triangular shape. This suggests that the and steps may not belong to the calcite-aqueous solution equilibrium system. Further increase in the saturation state (Ω ? 0.3) led to a lack of etch pit formation and dissolution primarily occurred at existing steps, in accordance with Teng (2004). Analysis of step kinetics at different temperatures yielded activation energies of 25 ± 6 kJ/mol and 14 ± 13 kJ/mol for obtuse and acute steps, respectively. The inconsistencies in etch pit morphology, step anisotropy, and step activation energies from the present study with those of studies far-from-equilibrium can be explained by increased influence of the backward reaction, or growth, near-equilibrium. We propose that the backward reaction occurs preferentially at the acute-acute kink sites. The kinetics and effective activation energies of near-equilibrium calcite dissolution presented in this work provide accurate experimental data under likely CO₂ sequestration conditions, and thus are crucial to the development of robust geochemical models that predict the long-term performance of mineral-trapped CO₂.     

Kinetics and mechanism of natural fluorapatite dissolution at 25 °C and pH from 3 to 12

The dissolution rates of natural fluorapatite (FAP), Ca₁₀(PO₄)₆F₂, were measured at 25 °C in mixed-flow reactors as a function of pH from 3.0 to 11.7, and aqueous calcium, phosphorus, and fluoride concentration. After an initial preferential Ca and/or F release, stoichiometric Ca, P, and F release was observed. Measured FAP dissolution rates decrease with increasing pH at 3 ? pH ? 7, FAP dissolution rates are pH independent at 7 ? pH ? 10, and FAP dissolution rates again decrease with increasing pH at pH ? 10. Measured FAP dissolution rates are independent of aqueous Ca, P, and F concentration at pH ≈ 3 and pH ≈ 10.Apatite dissolution appears to be initiated by the relatively rapid removal from the near surface of F and the Ca located in the M1 sites, via proton for Ca exchange reactions. Dissolution rates are controlled by the destruction of this F and Ca depleted surface layer. The destruction of this layer is facilitated by the adsorption/penetration of protons into the surface at acidic conditions, and by surface hydration at neutral and basic conditions. Taking into account these two parallel mechanisms, measured fluorapatite forward dissolution rates can be accurately described using

     

Experimental study and modeling of fluid reaction paths in the quartz-kyanite ± muscovite-water system at 0.7 GPa in the 350-550 °C range: Implications for Al selective transfer during metamorphism

In order to quantify Al transfer in response to fluid-mineral equilibration under evolving metamorphic conditions, isobaric (0.7 GPa) experiments were conducted in the 350-550 °C range. Disequilibrium was induced (1) by holding initially pure water and natural minerals (kyanite + quartz ± muscovite enclosed in a perforated inner capsule) under isothermal conditions and (2) by stepwise temperature variations. In all experiments, secondary Al-bearing phases crystallized in the external tube of a “tube-in-tube” setup (SEM characterization); they are interpreted as witnesses of the evolution of the fluid composition (fluid reaction path). These reaction paths and the subsequent amount of secondary crystallizations were modeled using thermodynamic data from SUPCRT92 and estimates of both starting-mineral dissolution rates and elemental diffusion coefficients from the literature. A major result is that the amount of aluminum transferred to secondary phases is a thousand times larger than the calculated Al concentration in the fluid. Although the crystallization of Al-bearing phases was expected as a response to a temperature decrease, the stepwise temperature increase (20 °C/day) also led to aluminum transfer towards secondary phases. In the course of re-equilibration, the fluid first becomes saturated with respect to aluminosilicates and then reaches silica saturation, due to the low solubility of Al-minerals. Consequently, aluminosilicates partly recrystallize in response to a temperature increase. Crystallization of secondary Al-phases in the external tube implies that aqueous aluminum was efficiently transported from the inner capsule, even in the pure Al₂O₃-SiO₂-H₂O system. Therefore, mass balance calculations considering a constant Al reference frame, i.e., postulating Al immobility, should be regarded with caution.     

Controls on gold solubility in arc magmas: An experimental study at 1000 °C and 4 kbar

In order to (1) explain the worldwide association between epithermal gold-copper-molybdenum deposits and arc magmas and (2) test the hypothesis that adakitic magmas would be Au-specialized, we have determined the solubility of Au at 4 kbar and 1000 °C for three intermediate magmas (two adakites and one calc-alkaline composition) from the Philippines. The experiments were performed over a fO₂ range corresponding to reducing (∼NNO−1), moderately oxidizing (∼NNO+1.5) and strongly oxidizing (∼NNO+3) conditions as measured by solid Ni-Pd-O sensors. They were carried out in gold containers, the latter serving also as the source of gold, in presence of variable amounts of H₂O and, in a few additional experiments, of S. Concentrations of Au in glasses were determined by LA-ICPMS. Gold solubility in melt is very low (30-240 ppb) but increases with fO₂ in a way consistent with the dissolution of gold as both Au¹⁺ and Au³⁺ species. In the S-bearing experiments performed at ∼NNO−1, gold solubility reaches much higher values, from ∼1200 to 4300 ppb, and seems to correlate with melt S content. No systematic difference in gold solubility is observed between the adakitic and the non-adakitic compositions investigated. Oxygen fugacity and the sulfur concentration in melt are the main parameters controlling the incorporation and concentration of gold in magmas. Certain adakitic and non-adakitic magmas have high fO₂ and magmatic S concentrations favorable to the incorporation and transport of gold. Therefore, the cause of a particular association between some arc magmas and Au-Cu-Mo deposits needs to be searched in the origin of those specialized magmas by involvement of Au- and S-rich protoliths. The subducted slab, which contains metal-rich massive sulfides, may constitute a potentially favorable protolith for the genesis of magmas specialized with respect to gold.     

A synthetic fluid inclusion study of copper solubility in hydrothermal brines from 525 to 725 °C and 0.3 to 1.7 GPa

Fluid inclusions were synthesized in a piston-cylinder apparatus under mineral-buffered conditions over a range of Cl concentration (0.29 to 11.3 mol kg⁻¹), temperature (525 to 725 °C), and pressure (0.3 to 1.7 GPa). All fluids were buffered by the mineral assemblage native copper + cuprite + talc + quartz. In situ fluid composition was determined by analysing individual fluid inclusions by LA-ICPMS and independently analysing the quench solution. The solubility data provide basic information necessary to model the high temperature behaviour of Cu in magmatic-hydrothermal systems. Copper concentrations up to ∼15 wt% were measured at 630 °C and 0.34 GPa. These results give an upper limit for Cu in natural fluids and support field-based observations of similar high Cu concentrations in fluids at near-magmatic conditions. Experimental evidence indicates that Cu⁺ may form neutral chloride complexes with the general stoichiometry with n up to 4, though n ? 2 is typical for the majority of the experimental conditions. At high pressure (>∼0.5 GPa) there is evidence that hydroxide species, e.g., CuOH⁰, become increasingly important and may predominate over copper(I)-chloride complexes. The roles of fluid mixing, cooling and decompression in ore-forming environments are also discussed.     

O and H diffusion in uraninite: Implications for fluid-uraninite interactions, nuclear waste disposal, and nuclear forensics

Diffusion coefficients for oxygen and hydrogen were determined from a series of natural uraninite-H₂O experiments between 50 and 700 °C. Under hydrous conditions there are two diffusion mechanisms: (1) an initial extremely fast-path diffusion mechanism that overprinted the oxygen isotopic composition of the entire crystals regardless of temperature and (2) a slower volume-diffusive mechanism dominated by defect clusters that displace or eject nearest neighbor oxygen atoms to form two interstitial sites and two partial vacancies, and by vacancy migration. Using the volume diffusion coefficients in the temperature range of 400-600 °C, diffusion coefficients for oxygen can be represented by D = 1.90e⁻⁵ exp (−123,382 J/RT) cm²/s and for temperatures between 100 and 300 °C the diffusion coefficients can be represented by D = 1.95e⁻¹⁰ exp (−62484 J/RT) cm²/s, where the activation energies for uraninite are 123.4 and 62.5 kJ/mol, respectively. Hydrogen diffusion in uraninite appears to be controlled by similar mechanisms as oxygen. Using the volume diffusion coefficients for temperatures between 50 and 700 °C, diffusion coefficients for hydrogen can be represented by D = 9.28e⁻⁶ exp (−156,528 J/RT) cm²/s for temperatures between 450 and 700 °C and D = 1.39e⁻¹⁴ exp (−34518 J/RT) cm²/s for temperatures between 50 and 400 °C, where the activation energies for uraninite are 156.5 and 34.5 kJ/mol, respectively.Results from these new experiments have implications for isotopic exchange during natural UO₂-water interactions. The exceptionally low δ¹⁸O values of natural uraninites (i.e. −32‰ to −19.5‰) from unconformity-type uranium deposits in Saskatchewan, in conjunction with theoretical and experimental uraninite-water and UO₃-water fractionation factors, suggest that primary uranium mineralization is not in oxygen isotopic equilibrium with coeval clay and silicate minerals. The low δ¹⁸O values have been interpreted as resulting from the low temperature overprinting of primary uranium mineralization in the presence of relatively modern meteoric fluids having δ¹⁸O values of ca. −18‰, despite petrographic and U-Pb isotope data that indicate limited alteration. Our data show that the anomalously low oxygen isotopic composition of the uraninite from the Athabasca Basin can be due to meteoric water overprinting under reducing conditions, and meteoric water or groundwater can significantly affect the oxygen isotopic composition of spent nuclear fuel in a geologic repository, with minimal change to the chemical composition or texture. Moreover, the rather fast oxygen and hydrogen diffusion coefficients for uraninite, especially at low temperatures, suggest that oxygen and hydrogen diffusion may impart characteristic isotopic signals that can be used to track the route of fissile material.     

Surface chemistry of kaolinite and Na-montmorillonite in aqueous electrolyte solutions at 25 and 60 °C: Experimental and modeling study

The aqueous interfacial chemistry of kaolinite and Na-montmorillonite samples was investigated by potentiometric measurements using acid/base continuous titrations and batch experiments at 25 and 60 °C. Using the batch experimental method, a continuous drift of pH was observed reflecting the mineral dissolution. Consequently, the continuous titration method appears to be the best way of studying solid surface reactions. For each clay mineral, the net proton surface excess/consumption was calculated as a function of pH and ionic strength (0.025, 0.1 and 0.5 M). At 25 °C, and according to the literature data, the pH corresponding to zero net proton consumption for montmorillonite appears to depend on ionic strength, whereas the value for kaolinite is constant and close to 5. Similar results are obtained at 60 °C, which suggests that the point of zero net proton consumption for clay minerals does not depend on temperature, at least up to 60 °C. On the other hand, the temperature rise induces a slight increase of the net proton surface excess. Finally, the diffuse double layer formalism (DDLM) is used to model the experimental data. The model involves two processes: the protonation/deprotonation of two types of edge sites (aluminol and silanol) and H⁺/Na⁺ exchange reactions on basal surfaces, while a tiny proportion of the negative structural charge remains uncompensated. This last process maintains a negative surface potential whatever the pH of the solution, which is in agreement with electrokinetic data.     

An in situ XAS study of copper(I) transport as hydrosulfide complexes in hydrothermal solutions (25-592 °C, 180-600 bar): Speciation and solubility in vapor and liquid phases

Chloride and hydrosulfide are the principal ligands assumed to govern transport of copper in hydrothermal fluids. Existing solubility experiments suggest that Cu(I)-hydrosulfide complexes are dominant compared to chloride complexes at low salinities in alkaline solutions (H₂S_(aq)/HS⁻ pH buffer), and may be important in transporting Cu in low density magmatic vapors, potentially controlling the liquid-vapor partitioning of Cu. This study provides the first in situ evidence of the solubility of copper sulfides and the nature and structure of the predominant Cu species in sulfur-containing fluids at temperatures up to 592 °C and pressures of 180-600 bar. XANES and EXAFS data show that at elevated T (?200 °C), Cu solubility occurs via a linear Cu complex. At 428 °C in alkaline solutions, Cu is coordinated by two sulfur atoms in a distorted linear coordination (angle ∼150-160°). This geometry is consistent with the species predicted by earlier solubility studies. In addition, in situ measurements of the solubility of chalcocite in 2 m NaHS solutions performed in this study are in remarkably good agreement with the solubilities calculated using available thermodynamic data for Cu(I)-hydrosulfide complexes, also supporting the interpretation of speciation in these studies and validating the extrapolation of low-T thermodynamic properties for to high P-T. Data on phase separation for the 2 m NaHS solution show that while significant amounts of copper can be partitioned into the vapor phase, there is no indication for preferential partitioning of Cu into the vapor. This is consistent with recent partitioning experiments conducted in autoclaves by Pokrovski et al. (2008a) and Simon et al. (2006). XANES data suggest that the species present in the low density phase is very similar to that present in the high density liquid, i.e., , although Cu(HS)(H₂S)⁰ cannot be excluded on the basis of XAS data.