Experimental determination of the stability and stoichiometry of sulphide complexes of silver(I) in hydrothermal solutions to 400°C

The solubility of silver sulphide (acanthite/argentite) has been measured in aqueous sulphide solutions between 25 and 400°C at saturated water vapour pressure and 500 bar to determine the stability and stoichiometry of sulphide complexes of silver(I) in hydrothermal solutions. The experiments were carried out in a flow-through autoclave, connected to a high-performance liquid chromatographic pump, titanium sampling loop, and a back-pressure regulator on line. Samples for silver determination were collected via the titanium sampling loop at experimental temperatures and pressures. The solubilities, measured as total dissolved silver, were in the range 1.0 × 10⁻⁷ to 1.30 × 10⁻⁴ mol kg⁻¹ (0.01 to 14.0 ppm), in solutions of total reduced sulphur between 0.007 and 0.176 mol kg⁻¹ and pH_T,p of 3.7 to 12.7. A nonlinear least squares treatment of the data demonstrates that the solubility of silver sulphide in aqueous sulphide solutions of acidic to alkaline pH is accurately described by the reactions0.5Ag₂S(s) + 0.5H₂S(aq) = AgHS(aq) K_s,1110.5Ag₂S(s) + 0.5H₂S(aq) + HS⁻ = Ag(HS)2− K_s,122Ag₂S(s) + 2HS⁻ = Ag₂S(HS)22− K_s,232where AgHS(aq) is the dominant species in acidic solutions, Ag(HS)2− under neutral pH conditions and Ag₂S(HS)22− in alkaline solutions. With increasing temperature the stability field of Ag(HS)2− increases and shifts to more alkaline pH in accordance with the change in the first ionisation constant of H₂S(aq). Consequently, Ag₂S(HS)22− is not an important species above 200°C. The solubility constant for the first reaction is independent of temperature to 300°C, with values in the range logK_s,111 = −5.79 (±0.07) to −5.59 (±0.09), and decreases to −5.92 (±0.16) at 400°C. The solubility constant for the second reaction increases almost linearly with inverse temperature from logK_s,122 = −3.97 (±0.04) at 25°C to −1.89 (±0.03) at 400°C. The solubility constant for the third reaction increases with temperature from logK_s,232 = −4.78 (±0.04) at 25°C to −4.57 (±0.18) at 200°C. All solubility constants were found to be independent of pressure within experimental uncertainties. The interaction between Ag⁺ and HS⁻ at 25°C and 1 bar to form AgHS(aq) has appreciable covalent character, as reflected in the exothermic enthalpy and small entropy of formation. With increasing temperature, the stepwise formation reactions become progressively more endothermic and are accompanied by large positive entropies, indicating greater electrostatic interaction. The aqueous speciation of silver is very sensitive to fluid composition and temperature. Below 100°C silver(I) sulphide complexes predominate in reduced sulphide solutions, whereas Ag⁺ and AgClOH⁻ are the dominant species in oxidised waters. In high-temperature hydrothermal solutions of seawater salinity, chloride complexes of silver(I) are most important, whereas in dilute hydrothermal fluids of meteoric origin typically found in active geothermal systems, sulphide complexes predominate. Adiabatic boiling of dilute and saline geothermal waters leads to precipitation of silver sulphide and removal of silver from solution. Conductive cooling has insignificant effects on silver mobility in dilute fluids, whereas it leads to quantitative loss of silver for geothermal fluids of seawater salinity.     

A spectrophotometric study of aqueous copper(I)-chloride complexes in LiCl solutions between 100 °C and 250 °C

Copper transport and deposition in highly saline hydrothermal fluids are controlled by the stability of copper(I) complexes with ligands such as chloride and hydrosulphide. However, our understanding of the behavior of copper(I)-chloride complexes at elevated temperatures and in highly saline brines is limited by the conditions of existing experimental studies where the maximum chloride concentration is 2 m. This paper presents the results of a study of copper(I)-chloride complexes at much higher chloride concentrations, 1.5 m to 9.1 m, using a UV spectrophotometric method. The UV spectra of copper(I)-bearing LiCl solutions were measured at temperatures between 100 °C and 250 °C at vapor-saturated pressures and quantitative interpretation of the spectra shows that CuCl₂⁻, CuCl₃²⁻, and CuCl₄³⁻ were present in the experimental solutions. The fitted logarithms of formation constants (log K) for CuCl₂⁻ are in good agreement with the previous results of solubility experiments reported by Xiao et al. (1998) and Liu et al. (2001). The log K values for CuCl₃²⁻ also agree with those of Liu et al. (2001) and theoretical estimates of Sverjensky et al. (1997). This study presents the first experimentally determined formation constants for CuCl₄³⁻, at temperatures greater than 25 °C, and indicates that this complex predominates at chloride concentrations greater than 5 m. Based on the new log K values generated from this study, the calculated chalcopyrite solubility in NaCl solutions indicates that in addition to cooling, fluid mixing (dilution of saline fluids) may be an important factor controlling the deposition of copper minerals from hydrothermal solutions.     

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.     

Carbonate complexation of yttrium and the rare earth elements in natural waters

Potentiometric measurements of Yttrium and Rare Earth Element (YREE) complexation by carbonate and bicarbonate indicate that the quality of carbonate complexation constants previously obtained via solvent exchange analyses are superior to characterizations obtained using solubility and adsorptive exchange analyses. The results of our analyses at 25°C are combined with the results of previous solvent exchange analyses to obtain YREE carbonate complexation constants over a wide range of ionic strength (0 ≤ I ≤3 molal). YREE carbonate complexation constants are reported for the following equilibria, M³⁺+nHCO₃⁻?M(CO₃)_n³⁻²ⁿ+nH⁺, where n = 1 or 2. Formation constants written in terms of HCO₃⁻ concentrations require only minor corrections for ion pairing relative to the corrections required for constants expressed in terms of CO₃²⁻ concentrations. Formation constants for the above complexation equilibria, _CO3^Hβ₁=[MCO₃⁺][H⁺][M³⁺]⁻¹[HCO₃⁻]⁻¹ and _CO3^Hβ₂=[M(CO₃)₂⁻][H⁺]²[M³⁺]⁻¹[HCO₃⁻]⁻², have very similar dependencies on ionic strength because the reaction MCO₃⁺+HCO₃⁻?M(CO₃)₂⁻+H⁺ is isocoulombic. Potentiometric analyses indicate that the dependence of log_CO3^Hβ₁ and log_CO3^Hβ₂ on ionic strength at 25°C is given as

(A)     

An experimental study of the solubility of baddeleyite (ZrO2) in fluoride-bearing solutions at elevated temperature

The solubility of baddeleyite (ZrO₂) and the speciation of zirconium have been investigated in HF-bearing aqueous solutions at temperatures up to 400 °C and pressures up to 700 bar. The data obtained suggest that in HF-bearing solutions zirconium is transported mainly in the form of the hydroxyfluoride species ZrF(OH)₃° and ZrF₂(OH)₂°. Formation constants determined for these species (Zr⁴⁺ + nF⁻ + mOH⁻ = ZrF_n(OH)_m°) range from 43.7 at 100 °C to 46.41 at 400 °C for ZrF(OH)₃°, and from 37.25 at 100 °C to 43.88 at 400 °C for ZrF₂(OH)₂°.Although the solubility of ZrO₂ is retrograde with respect to temperature, the measured concentrations of Zr are orders of magnitude higher than those predicted from theoretical extrapolations based on simple fluoride species (ZrF³⁺-ZrF₆²⁻). Model calculations performed for zircon show that zirconium can be transported by aqueous fluids in concentrations sufficient to account for the concentration of this metal at conditions commonly encountered in fluoride-rich natural hydrothermal systems.     

Estimates of the second dissociation constant of H2S from the surface sulfidation of crystalline sulfur

The adsorption of hydrogen sulfide (Γ_H2S) and protons (Γ_H⁺) on the surface of crystalline sulfur was investigated experimentally in H₂S-bearing solutions at temperatures of 25, 50, and 70°C, NaCl concentrations of 0.1 and 0.5 mol/dm⁻³ and log C_H⁺ values in the range −2.3 to −5. At all temperatures, the dominant process on the surface of the sulfur was deprotonation, and the average values of Γ_H2S were very close to the highest values determined for Γ_H⁺. This finding, combined with the lack of detectable proton adsorption in H₂S-free solutions, suggests that proton adsorption/desorption on the surface of sulfur occurs through formation of ≡ S − H₂S complexes in the presence of H₂S.We propose that this complexation represents sulfidation of the sulfur surface, a process analogous to hydroxylation of oxide surfaces, and that the sulfidation can be described by the reaction: ≡ S + H₂S = ≡SSH₂⁰ β° The deprotonation of the ≡ SH° complex occurs via the reaction: ≡ SSH₂⁰ = ≡SSH⁻ + H⁺ β⁻ Values of 2.9, 2.8, and 2.9 (± 0.23) were obtained for −log β⁻ at 25, 50, and 70°C, respectively. These data were employed to estimate the second dissociation constant for hydrogen sulfide in aqueous solutions using the extrapolation method proposed by Schoonen and Barnes (1988) and yielded corresponding values for the constant of 17.4 ± 0.3, 15.7, and 14.5, respectively. The value for 25°C is in very good agreement with the experimentally determined values of Giggenbach (1971) at 17 ± 0.1; Meyer et al. (1983) at 17 ± 1; Licht and Manassen (1987) at 17.6 ± 0.3; and Licht et al. (1990) at 17.1 ± 0.3.     

The hydrolysis of gold(I) in aqueous solutions to 600°c and 1500 bar

The solubility of gold has been measured in aqueous solutions at temperatures between 300 and 600°C and pressures from 500 to 1500 bar to determine the stability and stoichiometry of the hydroxy complexes of gold(I) in hydrothermal solutions. The experiments were carried out using a flow-through autoclave system. The solubilities, measured as total dissolved gold, were in the range 1.2 × 10⁻⁸ to 2.0 × 10⁻⁶ mol kg⁻¹ (0.002 to 0.40 mg kg⁻¹), in solutions of total dissolved sodium between 0.0 and 0.5 mol kg⁻¹, and total dissolved hydrogen between 4.0 × 10⁻⁶ and 4.0 × 10⁻⁴ mol kg⁻¹. At constant hydrogen molality, the solubility of gold increases with increasing temperature and decreases with increasing pressure. The solubilities were found to be independent of pH but increased with decreasing hydrogen molality at constant temperature and pressure. Consequently, gold dissolves in aqueous solutions of acidic to alkaline pH according to the reactionAu(s)+H₂O(l)=AuOH(aq)+0.5H₂(g) K_s,1The solubility constant, logK_s,1, increases with increasing temperature from a minimum of −8.76 (±0.18) at 300°C and 500 bar to a maximum of −7.50 (±0.11) at 500°C and 1500 bar and decreases to −7.61 (±0.08) at 600°C and 1500 bar. From the equilibrium solubility constant and the redox potential of gold, the formation constant to form AuOH(aq) was calculated. At 25°C the complex formation is characterised by an exothermic enthalpy and a positive entropy. With increasing temperature and decreasing pressure, the formation reaction becomes endothermic and is accompanied by a large positive entropy, indicating a greater electrostatic interaction between Au⁺ and OH⁻.     

UV-Vis spectrophotometric and XAFS studies of ferric chloride complexes in hyper-saline LiCl solutions at 25-90 °C

The speciation and thermodynamic properties of ferric chloride complexes in hydrothermal solutions and hypersaline brines are still poorly understood, despite the importance of this element as a micronutrient and ore-component. Available experimental data are limited to room temperature and relatively low chloride concentrations. This paper reports results of UV-Vis spectrophotometric and synchrotron XAFS experiments of ferric chloride complexes in chloride concentrations up to 15 m and at temperatures of 25-90 °C. Qualitative interpretation of the UV-Vis spectra shows that FeCl²⁺, FeCl₂⁺, FeCl_3(aq) and FeCl₄⁻ were present in the experimental solutions. As chloride concentrations increase, higher ligand number complexes become important with FeCl₄⁻ predominating in solutions containing more than 10 m at 25 °C. The predominance fields of FeCl_3(aq) and FeCl₄⁻ expand to lower Cl concentrations with increasing T. Both XANES and UV-Vis spectra reveal a major change in the geometry of the complex between FeCl₂⁺ and FeCl_3(aq). EXAFS data confirm that the number of chloride ligands increases with increasing chloride concentration and show that Fe³⁺, FeCl²⁺ and FeCl₂⁺ share an octahedral geometry. FeCl_3(aq) could be either tetrahedral or trigonal dipyramidal, while FeCl₄⁻ is expected to be tetrahedral. EXAFS data support a tetrahedral geometry for FeCl₄⁻, especially at 90 °C, but do not allow to distinguish between a tetrahedral or trigonal dipyramidal geometry for FeCl_3(aq) because of similar Fe-Cl distances. At room temperature, EXAFS data suggest that FeCl_3(aq) may be a mixture of octahedral and tetrahedral or trigonal dipyramidal forms.The room temperature formation constants for three ferric chloride complexes (FeCl₂⁺, FeCl_3(aq) and FeCl₄⁻) determined from the UV data are generally in good agreement with previous studies. Calculations based on the properties extrapolated to 300 °C show that hematite solubility is much higher than previously estimated, and that the high orders complexes FeCl_3(aq) and FeCl₄⁻ are important at high temperatures even in solutions with low chloride concentrations. The accuracy of these properties is limited by a poor understanding of activity-composition relationships in concentrated electrolytes, and by limitations in the available experimental techniques and extrapolation algorithms; however, the inclusion of higher order complexes in numerical models of ore transport and deposition allows for a more accurate qualitative prediction of Fe behaviour in hydrothermal and hypersaline systems.     

A spectrophotometric study of neodymium(III) complexation in chloride solutions

The formation constants of neodymium complexes in chloride solutions have been determined spectrophotometrically at temperatures of 25 to 250°C and a pressure of 50 bars. The simple ion, Nd³⁺, is dominant at 25°C, whereas NdCl²⁺ and NdCl₂⁺ are the dominant species at elevated temperatures. Equilibrium constants were calculated for the following reactions:Nd³⁺ + Cl⁻ = NdCl²⁺ β₁,Nd³⁺ + 2 · Cl⁻ = NdCl⁺₂ β₂.The values of β₁ were found to be identical within experimental error to the values reported by Gammons et al. (1996) but substantially different from those proposed by Stepanchikova and Kolonin (1999). The values of β₂ obtained in this study agree relatively well with those of Gammons et al. (1996); differences are greatest at intermediate temperature and reach a maximum of one half an order of magnitude at 200°C.Theoretical estimates of β₁ and β₂ by Haas et al. (1995) using the revised Helgeson-Kirkham-Flowers (HKF) equation of state predict lower stability of NdCl²⁺ and NdCl₂⁺ at temperatures above 150°C than determined in this study. A new fit to the HKF equation of state is therefore proposed, which yields values for β₁ and β₂ similar to those obtained experimentally.Using the formation constants reported in this study, we predict that typical seafloor hydrothermal vent fluids will contain a maximum concentration of Nd of ∼2 ppb. This value is several orders of magnitude lower than would be required to explain the levels of Nd mobility commonly reported for seafloor hydrothermal systems and suggests that other ligands may be more important than Cl in transporting rare earth elements in the Earth’s crust.     

An experimental study of the solubility and speciation of neodymium (III) fluoride in F-bearing aqueous solutions

The solubility of neodymium (III) fluoride was investigated at temperatures of 150, 200 and 250 °C, saturated water vapor pressure, and a total fluoride concentration (HF°_aq + F⁻) ranging from 2.0 × 10⁻³ to 0.23 mol/l. The results of the experiments show that Nd³⁺ and NdF²⁺ are the dominant species in solution at the temperatures investigated and were used to derive formation constants for NdF²⁺ and a solubility product for NdF₃. The solubility product of NdF₃(logK_sp=loga_Nd³⁺+3loga_F^-) is −24.4 ± 0.2, −22.8 ± 0.1, and −21.5 ± 0.2 at 250, 200 and 150 °C, respectively, and the formation constant of NdF²⁺(logβ=loga_NdF²⁺-loga_Nd³⁺-loga_F^-) is 6.8 ± 0.1, 6.2 ± 0.1, and 5.5 ± 0.2 at 250, 200 and 150 °C, respectively. The results of this study show that published theoretical predictions significantly overestimate the stability of NdF²⁺ and the solubility of NdF₃.The potential impact of the results on natural systems was evaluated for a hypothetical fluid with a composition similar to that responsible for REE mineralization in the Capitan pluton, New Mexico. In contrast to results obtained using the theoretical predictions of Haas [Haas J. R., Shock E. L., and Sassani D. C. (1995) Rare earth elements in hydrothermal systems: estimates of standard partial molal thermodynamic properties of aqueous complexes of the rare earth elements at high pressures and temperatures. Geochim. Cosmochim. Acta59, 4329-4350.], which indicate that NdF²⁺ is the dominant species in solution, calculations employing the data presented in this paper and previously published experimental data for chloride and sulfate species [Migdisov A. A., and Williams-Jones A. E. (2002) A spectrophotometric study of neodymium(III) complexation in chloride solutions. Geochim. Cosmochim. Acta66, 4311-4323; Migdisov A. A., Reukov V. V., and Williams-Jones A. E. (2006) A spectrophotometric study of neodymium(III) complexation in sulfate solutions at elevated temperatures. Geochim. Cosmochim. Acta70, 983-992.] show that neodymium chloride species predominate and that neodymium fluoride species are relatively unimportant. This suggests that accepted models for REE deposits that invoke fluoride complexation as the method of hydrothermal REE transport may need to be re-evaluated.     

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:

     

Experimental determination and analysis of the solubility of corundum in 0.1-molal CaCl2 solutions between 400 and 600°C at 0.6 to 2.0 kbar

Corundum (α-Al₂O₃) solubility was measured in 0.1-molal CaCl₂ solutions from 400 to 600°C between 0.6 and 2.0 kbar. The Al molality at 2 kbar increases from 3.1 × 10⁻⁴ at 400°C to 12.7 × 10⁻⁴ at 600°C. At 1 kbar, the solubility increases from 1.5 × 10⁻⁴m at 400°C to 3.4 × 10⁻⁴m at 600°C. These molalities are somewhat less than corundum solubility in pure H₂O (Walther, 1997) at 400°C but somewhat greater at 600°C. The distribution of species was computed considering the Al species Al(OH)₃⁰ and Al(OH)₄⁻, consistent with the solubility of corundum in pure H₂O of Walther (1997) and association constants reported in the literature. The calculated solubility was greater than that measured except at 600°C and 2.0 kbar, indicating that neutral-charged species interactions are probably important.A Setchénow model for neutral species resulted in poor fitting of the measured values at 1.0 kbar. This suggests that Al(OH)₃⁰ has a greater stability relative to Al(OH)₄⁻ than given by the models of Pokrovskii and Helgeson (1995) or Diakonov et al. (1996). The significantly lower Al molalities in CaCl₂ relative to those in NaCl solutions at the same concentration confirm the suggestions of Walther (2001) and others that NaAl(OH)₄⁰ rather than an Al-Cl complex must be significant in supercritical NaCl solutions to give the observed increase in corundum solubility with increasing NaCl concentrations.     

The solubility of nantokite (CuCl(s)) and Cu speciation in low-density fluids near the critical isochore: An in-situ XAS study

The solubility of nantokite (CuCl_(s)) and the structure of the predominant copper species in supercritical water (290-400 bar at 420 °C; 350-450 °C at 290 bar; 500 °C at 350 bar; density = 0.14-0.65 g/cm³) were investigated concurrently using synchrotron X-ray absorption spectroscopy (XAS) techniques. These conditions were chosen as they represent single phase solutions near the critical isochore, where the fluid density is intermediate of typical values for vapour and brine and is highly sensitive to even small changes in pressure. X-ray absorption near edge spectroscopy (XANES) and extended X-ray absorption spectroscopy (EXAFS) analyses show that aqueous copper occurs in a slightly distorted linear coordination in the solutions studied, with an average of 1.35(±0.3) Cl and 0.65(±0.3) O neighbours. The solubility of CuCl_(s) decreases exponentially with decreasing water density (i.e., decreasing pressure at constant temperature), in a manner similar to the solubility behaviour of salts such as NaCl in water vapour. Based on this similarity, an apparent equilibrium constant for the dissolution reaction of 0.5 ± 0.4 was calculated from a regression of the data at 420 °C, and it was determined that each Cu atom is solvated by approximately three water molecules. This indicates that under these conditions, copper solubility is controlled mainly by the structure of the second-shell hydration, which is essentially invisible to the XAS techniques used in this study.These results demonstrate that for a supercritical fluid near the critical isochore, decreasing pressure may initiate precipitation of copper even before boiling or phase separation. Such a process could be responsible for near-surface ore deposition in seafloor hydrothermal systems, where supercritical fluids experience rapid pressure changes during the transition between lithostatic and hydrostatic domains.     

Solubility of Pt and Pd sulfides and Au metal in aqueous bisulfide solutions

An experimental study of the solubility of Pt and Pd sulfides and Au metal in aqueous bisulfide solutions was conducted at temperatures from 200° to 350 °C and at saturated vapor pressure. A 500-mL Bridgemantype pressure vessel constructed of titanium, and equipped with a motor-driven magnetic stirrer was employed. The pH and the oxidation state were buffered by the coexistence of H₂S/HS^–/SO _inf4 ^sup2– . The pH at temperature was calculated to be in the range 5.91–9.43, and S was 0.3–2.2 m. Under the experimental conditions, the measured solubility of gold is about two to three orders of magnitude greater than that of either platinum and palladium, and the measured solubility of platinum is, in general, approximately equal to that of palladium, in molal units. The solubilities are found to be in the range: platinum 4–800 ppb, palladium 1–400 ppb, and gold 2–300 ppm. The solubility data can be modeled adequately using the following reactions: Au+H₂S+HH^–=Au(HS) ₂ ^– +1/2H₂ (K₁₄); PtS+HS^–+H⁺=Pt (HS) ₂ ⁰ (K₁₅); PdS+HS^–+H⁺=Pd (HS) ₂ ⁰ (K₁₆); PtS₂+H₂=Pt (HS) ₂ ⁰ (K₂₁).With equilibrium constants determined as follows (errors represent two standard deviations): Preliminary measurements of the solubilities of metallic Pt, Pd and Au as hydroxide complexes were also conducted using a second titanium pressure vessel, at temperatures of 200° to 350 °C and vapor saturation pressure, with pH and the oxidation state controlled or buffered by adding known amounts of NaOH and H₂ gas. The concentration of NaOH was in the range 0.01–1.3 m, and the partial pressure of H₂ at 200 °C was 62–275 bars, initially. Under the temperature and pressure conditions of these experiments, the solubility of platinum in 1 m NaOH solution is less than 100 ppb, that of palladium is less than 10 ppb and that of gold is less than 0.2 ppm; and in 0.01 m NaOH solutions, both Pt and Pd solubilities are less than 1 ppb. These data indicate that the contributions of hydroxide complexes to the total solubilities in the bisulfide runs, where the pH was in the range of 5.9–9.4, are negligible. The concentrations of both Pt and Pd as bisulfide complexes in the Salton Sea geothermal system predicted using the stability constants determined in this work agree very well with those values measured by McKibben et al. (1990). This calculation strongly suggests that the PGE are transported in moderately reducing, near neutral hydrothermal fluids as bisulfide complexes, as is gold. However, the much lower maximum solubility of the PGE relative to gold severely constrains models of re genesis, and may explain the relative rarity of hydrothermal PGE deposits compared to the relative abundance of hydrothermal Au deposits.     

A spectrophotometric study of aqueous Au(III) halide-hydroxide complexes at 25-80 °C

The mobility and transport of gold in low-temperature waters and brines is affected by the aqueous speciation of gold, which is sensitive in particular to pH, oxidation and halide concentrations. In this study, we use UV-Vis spectrophotometry to identify and measure the thermodynamic properties of Au(III) aqueous complexes with chloride, bromide and hydroxide. Au(III) forms stable square planar complexes with hydroxide and halide ligands. Based on systematic changes in the absorption spectra of solutions in three binary systems NaCl-NaBr, NaCl-NaOH and NaBr-NaOH at 25 °C, we derived log dissociation constants for the following mixed and end-member halide and hydroxide complexes: [AuCl₃Br]⁻, [AuCl₂Br₂]⁻, [AuBr₃Cl]⁻ and [AuBr₄]⁻; [AuCl₃(OH)]⁻, [AuCl₂(OH)₂]⁻, [AuCl(OH)₃]⁻ and [Au(OH)₄]⁻; and [AuBr₃(OH)]⁻, [AuBr₂(OH)₂]⁻ and [AuBr(OH)₃]⁻. These are the first reported results for the mixed chloride-bromide complexes. Increasing temperature to 80 °C resulted in an increase in the stability of the mixed chloride-bromide complexes, relative to the end-member chloride and bromide complexes. For the [AuCl_(4−n)(OH)_n]⁻ series of complexes (n = 0-4), there is an excellent agreement between our spectrophotometric results and previous electrochemical results of Chateau et al. [Chateau et al. (1966)]. In other experiments, the iodide ion (I⁻) was found to be unstable in the presence of Au(III), oxidizing rapidly to I₂(g) and causing Au to precipitate. Predicted Au(III) speciation indicates that Au(III) chloride-bromide complexes can be important in transporting gold in brines with high bromide-chloride ratios (e.g., >0.05), under oxidizing (atmospheric), acidic (pH < 5) conditions. Native gold solubility under atmospheric oxygen conditions is predicted to increase with decreasing pH in acidic conditions, increasing pH in alkaline conditions, increasing chloride, especially at acid pH, and increasing bromide for bromide/chloride ratios greater than 0.05. The results of our study increase the understanding of gold aqueous geochemistry, with the potential to lead to new methods for mineral exploration, hydrometallurgy and medicine.     

Experimental study of the stability of aluminate-borate complexes in hydrothermal solutions

Formation of aqueous aluminate-borate complexes was characterized at 25°C using ²⁷Al NMR spectroscopy, and at 50-200°C via measurements of gibbsite and boehmite solubility in the presence of boric acid. ²⁷Al spectra performed at pH = 9 in Al-B solution with m(B) = 0.02 show the presence of two peaks at 80.5 and 74.5 ppm which correspond to Al(OH)₄⁻ and a single Al-substituted Q¹_Al dimer, Al(OH)₃OB(OH)₂⁻, respectively. In 0.08 m and 0.2 m borate solution, a third peak appears at 68.5 ppm which can be assigned to the Q²_Al trimer Al(OH)₂O₂(B(OH)₂)₂⁻. These chemical shifts are close to those measured for Al(OH)₃OSi(OH)₃⁻ and Al(OH)₂O₂(Si(OH)₃)₂⁻ (74 and 69.5 ppm, respectively; Pokrovski et al., Min. Mag.62a (1998), 1194) which demonstrates the similar structure of Al-B and Al-Si complexes formed in alkaline solutions. Gibbsite and boehmite solubility were measured in weakly basic solutions as a function of boric acid concentration at 50°C and 78 to 200°C, respectively. Equilibrium was reached within several days at m(B) = 0.01-0.1, but more slowly at higher boron concentrations, and at 50°C and m(B) = 0.2, Al concentration increased continuously during at least 3 months as a result of the sluggish formation of Al-polyborates. The equilibrium constant of the reaction Al(OH)₄⁻ + B(OH)₃⁰_(aq) = Al(OH)₃OB(OH)₂⁻ + H₂O decreases very slowly with increasing temperature to 200°C. The log K values are 1.58 ± 0.10, 1.46 ± 0.10, 1.52 ± 0.15, and 1.25 ± 0.15 at 50, 78, 150 and 200°C, respectively, which result in the following values of the standard thermodynamic properties for this reaction: Δ_rG⁰ = −9.22 ± 3.25 kJ/mol, Δ_rH⁰ = −4.6 ± 2.5 kJ/mol, Δ_rS⁰ = 15.5 ± 6.9 J/mol K. The thermodynamic data generated in this study indicate that Al-B complexes can dominate aqueous aluminum speciation in solutions containing ≥0.7 g/L of boron at temperature to at least 400°C.     

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 .     

Kinetic and equilibrium Fe isotope fractionation between aqueous Fe(II) and Fe(III)

Equilibrium and kinetic Fe isotope fractionation between aqueous ferrous and ferric species measured over a range of chloride concentrations (0, 11, 110 mM Cl⁻) and at two temperatures (0 and 22°C) indicate that Fe isotope fractionation is a function of temperature, but independent of chloride contents over the range studied. Using ⁵⁷Fe-enriched tracer experiments the kinetics of isotopic exchange can be fit by a second-order rate equation, or a first-order equation with respect to both ferrous and ferric iron. The exchange is rapid at 22°C, ∼60-80% complete within 5 seconds, whereas at 0°C, exchange rates are about an order of magnitude slower. Isotopic exchange rates vary with chloride contents, where ferrous-ferric isotope exchange rates were ∼25 to 40% slower in the 11 mM HCl solution compared to the 0 mM Cl⁻ (∼10 mM HNO₃) solutions; isotope exchange rates are comparable in the 0 and 110 mM Cl⁻ solutions.The average measured equilibrium isotope fractionations, Δ_{Fe(III)-Fe(II)}, in 0, 11, and 111 mM Cl⁻ solutions at 22°C are identical within experimental error at +2.76±0.09, +2.87±0.22, and +2.76±0.06 ‰, respectively. This is very similar to the value measured by Johnson et al. (2002a) in dilute HCl solutions. At 0°C, the average measured Δ_{Fe(III)-Fe(II)} fractionations are +3.25±0.38, +3.51±0.14 and +3.56±0.16 ‰ for 0, 11, and 111 mM Cl⁻ solutions. Assessment of the effects of partial re-equilibration on isotope fractionation during species separation suggests that the measured isotope fractionations are on average too low by ∼0.20 ‰ and ∼0.13 ‰ for the 22°C and 0°C experiments, respectively. Using corrected fractionation factors, we can define the temperature dependence of the isotope fractionation from 0°C to 22°C as: where the isotopic fractionation is independent of Cl⁻ contents over the range used in these experiments. These results confirm that the Fe(III)-Fe(II) fractionation is approximately half that predicted from spectroscopic data, and suggests that, at least in moderate Cl⁻ contents, the isotopic fractionation is relatively insensitive to Fe-Cl speciation.     

Hydrosulfide/sulfide complexes of zinc to 250 °C and the thermodynamic properties of sphalerite

The solubility of synthetic ZnS_(cr) was measured at 25–250 °C and P = 150 bars as a function of pH in aqueous sulfide solutions (~ 0.015–0.15 m of total reduced sulfur). The solubility determinations were performed using a Ti flow-through hydrothermal reactor. The solubility of ZnS_(cr) was found to increase slowly with temperature over the whole pH range from 2 to ~ 10. The values of the Zn–S–HS complex stability constant, β, were determined for Zn(HS)₂⁰_(aq), Zn(HS)₃^?, Zn(HS)₄^2?, and ZnS(HS)^?. Based on the experimental values the Ryzhenko–Bryzgalin electrostatic model parameters for these stability constants were calculated, and the ZnS_(cr) solubility and the speciation of Zn in sulfide-containing hydrothermal solutions were evaluated. The most pronounced solubility increase, about 3 log units at m(S_total) = 0.1 for the temperatures from 25 to 250 °C, was found in acidic solutions (pH ~ 3 to 4) in the Zn(HS)₂⁰_(aq) predominance field. In weakly alkaline solutions, where Zn(HS)₃^? and Zn(HS)₄^2? are the dominant Zn–S–HS complexes, the ZnS_(cr) solubility increases by 1 log unit at the same conditions. It was found that ZnS(HS)^? and especially Zn(HS)₄^2? become less important in high temperature solutions. At 25 °C and m(S_total) = 0.1, these species dominate Zn speciation at pH > 7. At 100 °C and m(S_total) = 0.1, the maximum fraction of Zn(HS)₄^2? is only 20% of the total Zn concentration (i.e. at pH_t ~ 7.5), whereas at 350 °C and 3 <pH_t <10, the fraction of Zn(HS)₄^2? and ZnS(HS)^? is less than 0.05% and 2.5% respectively, of the total Zn concentration and Zn(HS)₂⁰ and Zn(HS)₃^? predominate. The measured equilibrium formation constants were combined with the literature data on the stability of Zn–Cl complexes in order to evaluate the concentration and speciation of Zn in chloride solutions. It was found that at acidic pH, and in more saline fluids having total chloride > 0.05 m, Zn–Cl complexes are responsible for hydrothermal Zn transport with no significant contribution of Zn–S–HS complexes. The hydrosulfide/sulfide complexes will play a more important role in lower salinity (< 0.05 m chloride) hydrothermal solutions which are characteristic of many epithermal ore depositing environments. The value of Δ_fG° (β-ZnS_(cr)) = ? 198.6 ± 0.2 kJ/mol at 25 °C was determined via solubility measurements of natural low-iron Santander (Spain) sphalerite.