An experimental study of Cobalt (II) complexation in Cl and H2S-bearing hydrothermal solutions

The speciation of cobalt (II) in Cl⁻ and H₂S-bearing solutions was investigated spectrophotometrically at temperatures of 200, 250, and 300 °C and a pressure of 100 bars, and by measuring the solubility of cobaltpentlandite at temperatures of 120-300 °C and variable pressures of H₂S. From the results of these experiments, it is evident that CoHS⁺ and predominate in the solutions except at 150 °C, for which the dominant chloride complex is CoCl₃⁻. The logarithms of the stability constant for CoHS⁺ show moderate variation with temperature, decreasing from 6.24 at 120 °C to 5.84 at 200 °C, and increasing to 6.52 at 300 °C. Formation constants for chloride species increase smoothly with temperature and at 300°C their logarithms reach 8.33 for , 6.44 for CoCl₃⁻, 4.94 to 5.36 for , and 2.42 for CoCl⁺. Calculations based on the composition of a model hydrothermal fluid (Ksp-Mu-Qz, KCl = 0.25 m, NaCl = 0.75 m, ΣS = 0.3 m) suggest that at temperatures ?200 °C, cobalt occurs dominantly as CoHS⁺, whereas at higher temperatures the dominant species is .     

Experimental determination of synthetic NdPO4 monazite end-member solubility in water from 21°C to 300°C: implications for rare earth element mobility in crustal fluids

The solubility of synthetic NdPO₄ monazite end-member has been determined experimentally from 21 to 300°C in aqueous solutions at pH = 2, and at 21°C and pH = 2 for GdPO₄. Measurements were performed in batch reactors, with regular solution sampling for pH measurement, rare earths and phosphorous analysis by inductively coupled plasma mass spectrometry (ICP-MS) coupled with a desolvation system. Scanning electron microscopy (SEM) and X-ray photoelectron spectroscopy (XPS) were employed to check that no reprecipitation of secondary phases occurred and that the mineral surfaces remained those of a monazite. Coupled with speciation calculations, measured solution compositions permitted the determination of NdPO₄ and GdPO₄ solubility products which are in general agreement with previous experimental determination on rhabdophane at 25°C, but showing that monazite is more than two orders of magnitude less soluble than inferred on the basis of previous thermodynamic estimates. The temperature evolution from 21 to 300°C of the equilibrium constant (K) of the NdPO₄ monazite end-member dissolution reaction given by:

     

A spectrophotometric study of neodymium(III) complexation in sulfate solutions at elevated temperatures

The formation constants of neodymium complexes in sulfate solutions have been determined spectrophotometrically at temperatures of 30-250 °C and a pressure of 100 bars. The dominant species in the solution are NdSO₄⁺ and Nd(SO₄)₂⁻, with the latter complex being more important at higher temperature. Equilibrium constants were calculated for the following reactions:

     

A potentiometric study of the stability of aqueous yttrium-acetate complexes from 25 to 175 °C and 1-1000 bar

The stability of yttrium-acetate (Y-Ac) complexes in aqueous solution was determined potentiometrically at temperatures 25-175 °C (at P_s) and pressures 1-1000 bar (at 25 and 75 °C). Measurements were performed using glass H⁺-selective electrodes in potentiometric cells with a liquid junction. The species YAc²⁺ and were found to dominate yttrium aqueous speciation in experimental solutions at 25-100 °C (log [Ac⁻] < −1.5, pH < 5.2), whereas at 125, 150 and 175 °C introduction of into the Y-Ac speciation model was necessary. The overall stability constants β_n were determined for the reaction

     

Zinc complexation in aqueous sulfide solutions: Determination of the stoichiometry and stability of complexes via ZnS(cr) solubility measurements at 100 °C and 150 bars

The solubility of ZnS_(cr) was measured at 100 °C, 150 bars in sulfide solutions as a function of sulfur concentration (m(S_total) = 0.02-0.15) and acidity (pH_t = 2-11). The experiments were conducted using a Ti flow-through hydrothermal reactor enabling the sampling of large volumes of solutions at experimental conditions, with the subsequent concentration and determination of trace quantities of Zn. Prior to the experiments, a long-term in situ conditioning of the solid phase was performed in order to attain the reproducible Zn concentrations (i.e. solubilities). The ZnS_(cr) solubility product was monitored in the course of the experiment. The following species were found to account for Zn speciation in solution: Zn²⁺ (pH_t < 3), (pH_t 3-4.5), (pH_t 5-8), and ZnS(HS)⁻ (pH_t > 8) (pH_t predominance regions are given for m(S_total) = 0.1). Solubility data collected in this study at pH_t > 3 were combined with the ZnS_(cr) solubility product determined at lower pH to yield the following equilibrium constants (t = 100 °C, P = 150 bars):

     

Uranophane dissolution and growth in CaCl2-SiO2(aq) test solutions

The dissolution and growth of uranophane [Ca(UO₂)₂(SiO₃OH)₂·5H₂O] have been examined in Ca- and Si-rich test solutions at low temperatures (20.5 ± 2.0 °C) and near-neutral pH (∼6.0). Uranium-bearing experimental solutions undersaturated and supersaturated with uranophane were prepared in matrices of ∼10⁻² M CaCl₂ and ∼10⁻³ M SiO₂(aq). The experimental solutions were reacted with synthetic uranophane and analyzed periodically over 10 weeks. Interpretation of the aqueous solution data permitted extraction of a solubility constant for the uranophane dissolution reaction and standard state Gibbs free energy of formation for uranophane ( kJ mol⁻¹).     

Reaction pathways for quartz dissolution determined by statistical and graphical analysis of macroscopic experimental data

In light of recent work on the reactivity of specific sites on large (hydr)oxo-molecules and the evolution of surface topography during dissolution, we examined the ability to extract molecular-scale reaction pathways from macroscopic dissolution and surface charge measurements of powdered minerals using an approach that involved regression of multiple datasets and statistical graphical analysis of model fits. The test case (far-from-equilibrium quartz dissolution from 25 to 300 °C, pH 1-12, in solutions with [Na⁺] ? 0.5 M) avoids the objections to this goal raised in these recent studies. The strategy was used to assess several mechanistic rate laws, and was more powerful in distinguishing between models than the statistical approaches employed previously. The best-fit model included three mechanisms—two involving hydrolysis of Si centers by H₂O next to neutral (>Si-OH⁰) and deprotonated (>Si-O⁻) silanol groups, and one involving hydrolysis of Si centers by OH⁻. The model rate law is

     

Stability of chloridogold(I) complexes in aqueous solutions from 300 to 600°C and from 500 to 1800 bar

The solubility of gold has been measured in the system H₂O+H₂+HCl+NaCl+NaOH at temperatures from 300 to 600°C and pressures from 500 to 1800 bar in order to determine the stability and stoichiometry of chloride complexes of gold(I) in hydrothermal solutions. The experiments were carried out in a flow-through autoclave system. This approach permitted the independent determination of the concentrations of all critical aqueous components in solution for the determination of the stability and stoichiometry of gold(I) complexes. The solubilities (i.e. total dissolved gold) were in the range 9.9 × 10⁻⁹ to 3.26 × 10⁻⁵ mol kg⁻¹ (0.002-6.42 mg kg⁻¹) in solutions of total dissolved chloride between 0.150 and 1.720 mol kg⁻¹, total dissolved sodium between 0.000 and 0.975 mol kg⁻¹ and total dissolved hydrogen between 4.34 × 10⁻⁶ and 7.87 × 10⁻⁴ mol kg⁻¹. A nonlinear least squares treatment of the data demonstrates that the solubility of gold in chloride solutions is accurately described by the reactions,

     

Experimental study of gold-hydrosulphide complexing in aqueous solutions at 350-500°C, 500 and 1000 bars using mineral buffers

The solubility of gold was measured in KCl solutions (0.001-0.1 m) at near-neutral to weakly acidic pH in the presence of the K-feldspar-muscovite-quartz, andalusite-muscovite-quartz, and pyrite-pyrrhotite-magnetite buffers at temperatures 350 to 500°C and pressures 0.5 and 1 kbar. These mineral buffers were used to simultaneously constrain pH, f(S₂), and f(H₂). The experiments were performed using a CORETEST flexible Ti-cell rocking hydrothermal reactor enabling solution sampling at experimental conditions. Measured log m(Au) (mol/kg H₂O) ranges from −7.5 at weakly acid pH to −5.9 in near-neutral solutions, and increases slightly with temperature. Gold solubility in weakly basic and near-neutral solutions decreases with decreasing pH at all temperatures, which implies that Au(HS)₂ is the dominant Au species in solution. In more acidic solutions, solubility is independent of pH. Comparison of the experimentally measured solubilities with literature values for Au hydrolysis constants demonstrates that at 350°C dominates Au aqueous speciation at the weakly acidic pH and f(S₂)/f(H₂) conditions imposed by the pyrite-pyrrhotite-magnetite buffer. In contrast, at temperatures >400°C becomes less important and predominates in weakly acid solutions. Solubility data collected in this study were used to calculate the following equilibrium reaction constants:

     

An experimental study of magnesite dissolution rates at neutral to alkaline conditions and 150 and 200 °C as a function of pH, total dissolved carbonate concentration, and chemical affinity

Steady-state magnesite dissolution rates were measured in mixed-flow reactors at 150 and 200 °C and 4.6 < pH < 8.4, as a function of ionic strength (0.001 M ? I ? 1 M), total dissolved carbonate concentration (10⁻⁴ M < ΣCO₂ < 0.1 M), and distance from equilibrium. Rates were found to increase with increasing ionic strength, but decrease with increasing temperature from 150 to 200 °C, pH, and aqueous CO₃²⁻ activity. Measured rates were interpreted using the surface complexation model developed by Pokrovsky et al. (1999a) in conjunction with transition state theory (Eyring, 1935). Within this formalism, magnesite dissolution rates are found to be consistent with

     

Complexation of Cu in Hydrothermal NaCl Brines: Ab initio molecular dynamics and energetics

Chloride complexation of Cu⁺ controls the solubility of copper(I) oxide and sulfide ore minerals in hydrothermal and diagenetic fluids. Solubility measurements and optical spectra of high temperature CuCl solutions have been interpreted as indicating the formation of CuCl, , and complexes. However, no other monovalent cation forms tri- and tetrachloro complexes. EXAFS spectra of high temperature Cu-Cl solutions, moreover, appear to show only CuCl and complexes at T > 100 °C. To reconcile these results, I investigated the nature and stability of Cu-Cl complexes using ab initio cluster calculations and ab initio (Car-Parrinello) molecular dynamics simulations for CuCl-NaCl-H₂O systems at 25 to 450 °C. Ab initio molecular dynamic simulations of 1 m CuCl in a 4 m Cl solution give a stable complex at 25 °C over 4 ps but show that the third Cl⁻ is weakly bound. When the temperature is increased along the liquid-vapour saturation curve to 125 °C, the complex dissociates into and Cl⁻; only forms at 325 °C and 1 kbar. Even in a 15.6 m Cl brine at 450 °C, only the complex forms over a 4 ps simulation run.Cluster calculations with a static dielectric continuum solvation field (COSMO) were used in an attempt directly estimate free energies of complex formation in aqueous solution. Consistent with the MD simulations, the complex is slightly stable at 25 °C but decreases in stability with decreasing dielectric constant (ε). The complex is predicted to be unstable at 25 °C and becomes increasingly unstable with decreasing dielectric constant. In hydrothermal fluids (ε < 30) both the and complexes are unstable to dissociation into and Cl⁻.The results obtained here are at odds with recent equations of state that predict and complexes are the predominant species in hydrothermal brines. In contrast, I predict that only complexes will be significant at T > 125 °C, even in NaCl-saturated brines. The high-temperature (T > 125 °C) optical spectra of CuCl solutions and solubility measurements of Cu minerals in Cl-brines need to be reinterpreted in terms of only the CuCl and complexes.     

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.     

Effect of pH and organic ligands on the kinetics of smectite dissolution at 25 °C

Forward dissolution rates of Na-Montmorillonite (Wyoming) SWy-2 smectite (Ca_0.06Na_0.56)[Al_3.08Fe(III)_0.38Mg_0.54] [Si_7.93 Al_0.07]O₂₀(OH)₄ were measured at 25 °C in a mixed-flow reactor equipped with interior dialysis compartment (6-8 kDa membrane) as a function of pH (1-12), dissolved carbonate (0.5-10 mM), phosphate (10⁻⁵ to 0.03 M), and nine organic ligands (acetate, oxalate, citrate, EDTA, alginate, glucuronic acid, 3,4-dihydroxybenzoic acid, gluconate, and glucosamine) in the concentration range from 10⁻⁵ to 0.03 M. In organic-free solutions, the Si-based rates decrease with increasing pH at 1 ? pH ? 8 with a slope close to −0.2. At 9 ? pH ? 12, the Si-based rates increase with a slope of ∼0.3. In contrast, non-stoichiometric Mg release weakly depends on pH at 1 ? pH ? 12 and decreases with increasing pH. The empirical expression describing Si-release rates [R, mol/cm²/s] obtained in the present study at 25 °C, I = 0.01 M is given by

     

Molybdic acid ionisation under hydrothermal conditions to 300 °C

This UV spectrophotometric study was aimed at providing precise, experimentally derived thermodynamic data for the ionisation of molybdic acid (H₂MoO₄) from 30 to 300 °C and at equilibrium saturated vapour pressures. The determination of the equilibrium constants and associated thermodynamic parameters were facilitated by spectrophotometric measurements using a specially designed high temperature optical Ti-Pd flow-through cell with silica glass windows.The following van’t Hoff isochore equations describe the temperature dependence of the first and second ionisation constants of molybdic acid up to 300 °C:

     

Antimony speciation in saline hydrothermal fluids: A combined X-ray absorption fine structure spectroscopy and solubility study

Solubility of senarmontite (Sb₂O₃, cubic) in pure water and NaCl-HCl aqueous solutions, and local atomic structure around antimony in these fluids were characterized using in situ X-ray absorption fine structure (XAFS) spectroscopy at temperatures to 450 °C and pressures to 600 bars. These experiments were performed using a new X-ray cell which allows simultaneous measurement of the absolute concentration of the absorbing element in the fluid, and atomic environment around the absorber. Results show that aqueous Sb(III) speciation is dominated by the complex in pure water, mixed Sb-hydroxide-chloride complexes in acidic NaCl-HCl solutions (2 m NaCl-0.1 m HCl), and by Sb-chloride species in concentrated HCl solutions (3.5 m HCl). Interatomic Sb-O and Sb-Cl distances in these complexes range from 1.96 to 1.97 Å and from 2.37 to 2.47 Å, respectively. These structural data, together with senarmontite solubility determined from XAFS spectra, were complemented by batch-reactor measurements of senarmontite and stibnite (Sb₂S₃, rhombic) solubilities over a wide range of HCl and NaCl concentrations from 300 to 400 °C. Analysis of the whole dataset shows that Sb(III) speciation in high-temperature moderately acid (pH > 2-3) Cl-rich fluids is dominated by mixed hydroxy-chloride species like Sb(OH)₂Cl° and Sb(OH)₃Cl⁻, but other species containing two or three Cl atoms appear at higher acidities and moderate temperatures (?300 °C). Calculations using stability constants retrieved in this study indicate that mixed hydroxy-chloride complexes control antimony transport in saline high-temperature ore fluids at acidic conditions. Such species allow for a more effective Sb partitioning into the vapor phase during boiling and vapor-brine separation processes occurring in magmatic-hydrothermal systems. Antimony hydroxy-chloride complexes are however minor in the neutral low- to moderate-temperature solutions (?250-300 °C) typical of Sb deposits formation; the antimony speciation in these systems is dominated by Sb(OH)₃ and potentially Sb-sulfide species.     

Synthesis, characterization and thermochemistry of a Pb-jarosite

The enthalpy of formation from the elements of a well-characterized synthetic Pb-jarosite sample corresponding to the chemical formula (H₃O)_0.74Pb_0.13Fe_2.92(SO₄)₂(OH)_5.76(H₂O)_0.24 was measured by high temperature oxide melt solution calorimetry. This value ( = −3695.9 ± 9.7 kJ/mol) is the first direct measurement of the heat of formation for a lead-containing jarosite. Comparison to the thermochemical properties of hydronium jarosite and plumbojarosite end-members strongly suggests the existence of a negative enthalpy of mixing possibly related to the nonrandom distribution of Pb²⁺ ions within the jarosite structure. Based on these considerations, the following thermodynamic data are proposed as the recommended values for the enthalpy of formation from the elements of the ideal stoichiometric plumbojarosite Pb_0.5Fe₃(SO₄)₂(OH)₆:  = −3118.1 ± 4.6 kJ/mol,  = −3603.6 ± 4.6 kJ/mol and S° = 376.6 ± 4.5 J/(mol K). These data should prove helpful for the calculation of phase diagrams of the Pb-Fe-SO₄-H₂O system and for estimating the solubility product of pure plumbojarosite. For illustration, the evolution of the estimated solubility product of ideal plumbojarosite as a function of temperature in the range 5-45 °C was computed (Log(K_sp) ranging from −24.3 to −26.2). An Eh-pH diagram is also presented.     

Dawsonite synthesis and reevaluation of its thermodynamic properties from solubility measurements: Implications for mineral trapping of CO2

Over the last decade, a significant research effort has focused on determining the feasibility of sequestering large amounts of CO₂ in deep, permeable geologic formations to reduce carbon dioxide emissions to the atmosphere. Most models indicate that injection of CO₂ into deep sedimentary formations will lead to the formation of various carbonate minerals, including the common phases calcite (CaCO₃), dolomite (CaMg(CO₃)₂), magnesite (MgCO₃), siderite (FeCO₃), as well as the far less common mineral, dawsonite (NaAlCO₃(OH)₂). Nevertheless, the equilibrium and kinetics that control the precipitation of stable carbonate minerals are poorly understood and few experiments have been performed to validate computer codes that model CO₂ sequestration.In order to reduce this uncertainty we measured the solubility of synthetic dawsonite according to the equilibrium: , from under- and oversaturated solutions at 50-200 °C in basic media at 1.0 mol · kg⁻¹ NaCl. The solubility products (Q_s) obtained were extrapolated to infinite dilution to obtain the solubility constants (. Combining the fit of these values and fixing  at 25 °C, which was derived from the calorimetric data of Ferrante et al. [Ferrante, M.J., Stuve, J.M., and Richardson, D.W., 1976. Thermodynamic data for synthetic dawsonite. U.S. Bureau of Mines Report Investigation, 8129, Washington, D.C., 13p.], the following thermodynamic parameters for the dissolution of dawsonite were calculated at 25 °C: , and . Subsequently, we were able to derive values for the Gibbs energy of formation (, enthalpy of formation ( and entropy ( of dawsonite. These results are within the combined experimental uncertainties of the values reported by Ferrante et al. (1976). Predominance diagrams are presented for the dawsonite/boehmite and dawsonite/bayerite equilibria at 100 °C in the presence of a saline solution with and without silica-containing minerals.     

Dissolution kinetics of synthetic Na-smectite. An integrated experimental approach

The effect of pH and Gibbs energy on the dissolution rate of a synthetic Na-montmorillonite was investigated by means of flow-through experiments at 25 and 80 °C at pH of 7 and 9. The dissolution reaction took place stoichiometrically at 80 °C, whereas at 25 °C preferential release of Mg over Si and Al was observed. The TEM-EDX analyses (transmission electronic microscopy with quantitative chemical analysis) of the dissolved synthetic phase at 25 °C showed the presence of newly formed Si-rich phases, which accounts for the Si deficit. At low temperature, depletion of Si concentration was attributed to incongruent clay dissolution with the formation of detached Si tetrahedral sheets (i.e., alteration product) whereas the Al behaviour remains uncertain (e.g., possible incorporation into Al-rich phases). Hence, steady-state rates were based on the release of Mg. Ex situ AFM measurements were used to investigate the variations in reactive surface area. Accordingly, steady-state rates were normalized to the initial edge surface area (11.2 m² g⁻¹) and used to propose the dissolution rate law for the dissolution reactions as a function of ΔG_r at 25 °C and pH∼9:

     

The effect of sulfur on the partitioning of Ni and other first-row transition elements between olivine and silicate melt

The effect of sulfur dissolved as sulfide (S²⁻) in silicate melts on the activity coefficients of NiO and some other oxides of divalent cations (Ca, Cr, Mn, Fe and Co) has been determined from olivine/melt partitioning experiments at 1400 °C in six melt compositions in the system CaO-MgO-Al₂O₃-SiO₂ (CMAS), and in derivatives of these compositions at 1370 °C, obtained from the six CMAS compositions by substituting Fe for Mg (FeCMAS). Amounts of S²⁻ were varied from zero to sulfide saturation, reaching 4100 μg g⁻¹ S in the most sulfur-rich silicate melt. The sulfide solubilities compare reasonably well with those predicted from the parameterization of the sulfide capacity of silicate melts at 1400 °C of O’Neill and Mavrogenes (2002), although in detail systematic deviations indicate that a more sophisticated model may improve the prediction of sulfide capacities.The results show a barely discernible effect of S²⁻ in the silicate melt on Fe, Co and Ni partition coefficients, and also surprisingly, a tiny but resolvable effect on Ca partitioning, but no detectable effect on Cr, Mn or some other lithophile incompatible elements (Sc, Ti, V, Y, Zr and Hf). Decreasing Mg# of olivine (reflecting increasing FeO in the system) has a significant influence on the partitioning of several of the divalent cations, particularly Ca and Ni. We find a remarkably systematic correlation between and the ionic radius of M²⁺, where M = Ca, Cr, Mn, Fe, Co or Ni, which is attributable to a simple relationship between size mismatch and excess free energies of mixing in Mg-rich olivine solid solutions.Neither the effect of S²⁻ nor of Mg#^ol is large enough by an order of magnitude to account for the reported variations of obtained from electron microprobe analyses of olivine/glass pairs from mid-ocean ridge basalts (MORBs). Comparing these MORB glass analyses with the Ni-MgO systematics of MORB from other studies in the literature, which were obtained using a variety of analytical techniques, shows that these electron microprobe analyses are anomalous. We suggest that the reported variation of with S content in MORB is an analytical artifact.Mass balance of melt and olivine compositions with the starting compositions shows that dissolved S²⁻ depresses the olivine liquidus of haplobasaltic silicate melts by 5.8 × 10⁻³ (±1.3 × 10⁻³) K per μg g⁻¹ of S²⁻, which is negligible in most contexts. We also present data for the partitioning of some incompatible trace elements (Sc, Ti, Y, Zr and Hf) between olivine and melt. The data for Sc and Y confirm previous results showing that and decrease with increasing SiO₂ content of the melt. Values of average 0.01 with most falling in the range 0.005-0.015. Zr and Hf are considerably more incompatible than Ti in olivine, with and about 10⁻³. The ratio / is well constrained at 0.611 ± 0.016.