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

     

Coupled phase and aqueous species equilibrium of the H2O-CO2-NaCl-CaCO3 system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to saturation of halite

A model is developed for the calculation of coupled phase and aqueous species equilibrium in the H₂O-CO₂-NaCl-CaCO₃ system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to saturation of halite. The vapor-liquid-solid (calcite, halite) equilibrium together with the chemical equilibrium of H⁺, Na⁺, Ca²⁺, , Ca(OH)⁺, OH⁻, Cl⁻, , , CO_2(aq) and CaCO_3(aq) in the aqueous liquid phase as a function of temperature, pressure, NaCl concentrations, CO_2(aq) concentrations can be calculated, with accuracy close to those of experiments in the stated T-P-m range, hence calcite solubility, CO₂ gas solubility, alkalinity and pH values can be accurately calculated. The merit and advantage of this model is its predictability, the model was generally not constructed by fitting experimental data.One of the focuses of this study is to predict calcite solubility, with accuracy consistent with the works in previous experimental studies. The resulted model reproduces the following: (1) as temperature increases, the calcite solubility decreases. For example, when temperature increases from 273 to 373 K, calcite solubility decreases by about 50%; (2) with the increase of pressure, calcite solubility increases. For example, at 373 K changing pressure from 10 to 500 bar may increase calcite solubility by as much as 30%; (3) dissolved CO₂ can increase calcite solubility substantially; (4) increasing concentration of NaCl up to 2 m will increase calcite solubility, but further increasing NaCl solubility beyond 2 m will decrease its solubility.The functionality of pH value, alkalinity, CO₂ gas solubility, and the concentrations of many aqueous species with temperature, pressure and NaCl_(aq) concentrations can be found from the application of this model. Online calculation is made available on www.geochem-model.org/models/h2o_co2_nacl_caco3/calc.php.     

The system H2O-NaCl. Part II: Correlations for molar volume, enthalpy, and isobaric heat capacity from 0 to 1000 °C, 1 to 5000 bar, and 0 to 1 XNaCl

A set of correlations for the volumetric properties and enthalpies of phases in the system H₂O-NaCl as a function of temperature, pressure, and composition has been developed that yields accurate values from 0 to 1000 °C, 1 to 5000 bar, and 0 to 1 X_NaCl. The volumetric properties of all fluid phases from low-density vapor to hydrous salt melts and single-phase binary fluids at high pressures and temperatures, can be described by a simple equation

     

Dissociation constants and speciation in aqueous Li2SO4 and K2SO4 from measurements of electrical conductance to 673 K and 29 MPa

The electrical conductivities of aqueous solutions of Li₂SO₄ and K₂SO₄ have been measured at 523-673 K at 20-29 MPa in dilute solutions for molalities up to 2 × 10⁻² mol kg⁻¹. These conductivities have been fitted to the conductance equation of Turq, Blum, Bernard, and Kunz with a consensus mixing rule and mean spherical approximation activity coefficients. In the temperature interval 523-653 K, where the dielectric constant, ε, is greater than 14, the electrical conductance data can be fitted by a solution model which includes ion association to form , , and , where M is Li or K. The adjustable parameters of this model are the first and second dissociation constants of the M₂SO₄. For the 673 K and 300 kg m⁻³ state point where the Coulomb interactions are the strongest (dielectric constant, ε = 5), models with more extensive association give good fits to the data. In the case of the Li₂SO₄ model, including the multi-ion associate, , gave an extremely good fit to the conductance data.     

A thermodynamic model for the prediction of phase equilibria and speciation in the H2O-CO2-NaCl-CaCO3-CaSO4 system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to halite saturation

A thermodynamic model is developed for the calculation of both phase and speciation equilibrium in the H₂O-CO₂-NaCl-CaCO₃-CaSO₄ system from 0 to 250 °C, and from 1 to 1000 bar with NaCl concentrations up to the saturation of halite. The vapor-liquid-solid (calcite, gypsum, anhydrite and halite) equilibrium together with the chemical equilibrium of H⁺,Na⁺,Ca²⁺, , , and CaSO_4(aq) in the aqueous liquid phase as a function of temperature, pressure and salt concentrations can be calculated with accuracy close to the experimental results.Based on this model validated from experimental data, it can be seen that temperature, pressure and salinity all have significant effects on pH, alkalinity and speciations of aqueous solutions and on the solubility of calcite, halite, anhydrite and gypsum. The solubility of anhydrite and gypsum will decrease as temperature increases (e.g. the solubility will decrease by 90% from 360 K to 460 K). The increase of pressure may increase the solubility of sulphate minerals (e.g. gypsum solubility increases by about 20-40% from vapor pressure to 600 bar). Addition of NaCl to the solution may increase mineral solubility up to about 3 molality of NaCl, adding more NaCl beyond that may slightly decrease its solubility. Dissolved CO₂ in solution may decrease the solubility of minerals. The influence of dissolved calcite on the solubility of gypsum and anhydrite can be ignored, but dissolved gypsum or anhydrite has a big influence on the calcite solubility. Online calculation is made available on www.geochem-model.org/model.     

Gold and copper partitioning in magmatic-hydrothermal systems at 800 °C and 100 MPa

Porphyry-type ore deposits sometimes contain fluid inclusion compositions consistent with the partitioning of copper and gold into vapor relative to coexisting brine at the depositional stage. However, this has not been reproduced experimentally at magmatic conditions. In an attempt to determine the conditions under which copper and gold may partition preferentially into vapor relative to brine at temperatures above the solidus of granitic magmas, we performed experiments at 800 °C, 100 MPa, oxygen fugacity () buffered by Ni-NiO, and fixed at either 3.5 × 10⁻² by using intermediate solid solution-pyrrhotite, or 1.2 × 10⁻⁴ by using intermediate solid solution-pyrrhotite-bornite. The coexisting vapor (∼3 wt.% NaCl eq.) and brine (∼68 wt.% NaCl eq.) were composed initially of NaCl + KCl + HCl + H₂O, with starting HCl set to <1000 μg/g in the aqueous mixture. Synthetic vapor and brine fluid inclusions were trapped at run conditions and subsequently analyzed by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Our experiments demonstrate that copper and gold partitioned strongly into the magmatic volatile phase(s) (MVP) (i.e., vapor or brine) relative to a silicate melt over the entire imposed range of . Nernst style partition coefficients between coexisting brine (b) and melt (m), D^b/m (±1σ), range from 3.6(±2.2) × 10¹ to 4(±2) × 10² for copper and from 1.2(±0.6) × 10² to 2.4(±2.4) × 10³ for gold. Partition coefficients between coexisting vapor (v) and melt, D^v/m range from 2.1 ± 0.7 to 18 ± 5 and 7(±3) × 10¹ to 1.6(±1.6) × 10² for copper and gold, respectively. Partition coefficients for all experiments between coexisting brine and vapor, D^b/v (±1σ), range from 7(±2) to 1.0(±0.4) × 10² and 1.7(±0.2) to 15(±2) for copper and gold, respectively. Observed average D^b/v at an of 1.2 × 10⁻⁴ were elevated, 95(±5) and 15 ± 1 for copper and gold, respectively, relative to those at the higher of 3.5 × 10⁻² where D^b/v were 10(±5) for copper and 7(±6) for gold. Thus, there is an inverse relationship between the and the D^b/v for both copper and gold with increasing resulting in a decrease in the D^b/v signifying increased importance of the vapor phase for copper and gold transport. This suggests that copper and gold may complex with volatile S-species as well as Cl-species at magmatic conditions, however, none of the experiments of our study at 800 °C and 100 MPa had a D^b/v ? 1. We did not directly determine speciation, but infer the existence of some metal-sulfur complexes based on the reported data. We suggest that copper and gold partition preferentially into the brine in most instances at or above the wet solidus. However, in most systems, the mass of vapor is greater than the mass of brine, and vapor transport of copper and gold may become more important in the magmatic environment at higher , lower , or near the critical point in a salt-water system. A D^b/v ? 1 at subsolidus hydrothermal conditions may also occur in response to changes in temperature, , , and/or acidity.Additionally, both copper and gold were observed to partition into intermediate solid solution and bornite much more strongly than into vapor, brine or silicate melt. This suggests that, although vapor and brine are both efficient at removing copper and gold from a silicate melt, the presence of Cu-Fe sulfides can sequester a substantial portion of the copper and gold contained within a silicate melt if the Cu-Fe sulfides are abundant.     

The partitioning of molybdenum(VI) between aqueous liquid and vapour at temperatures up to 370 °C

We have conducted experiments to evaluate the vapour-liquid fractionation of Mo(VI) in the system MoO₃-NH₃-H₂O at 300-370 °C and saturated vapour pressure, using a two-chamber autoclave that allows separate trapping of the vapour and liquid. The measured total Mo concentrations in each phase were used to calculate a distribution coefficient, , which increases as the density of the vapour approaches that of the liquid, and is greater than one for pH ? 4. Molybdenum speciation in the vapour is described by a single complex, MoO₃H₂O. By contrast, thermodynamic modeling of the distribution of Mo species in the liquid indicates that bimolybdate (HMoO₄⁻) is the dominant aqueous species at the conditions of our experiments, and that molybdate (MoO₄²⁻) and molybdic acid (H₂MoO₄⁰) are present in smaller quantities. As vapour-liquid fractionation occurs between neutral species, it is governed by the reaction H₂MoO₄⁰(aq) = MoO₃ · H₂O(g). Fractionation is therefore controlled by the concentration of H₂MoO₄⁰ in the liquid, which increases with increasing temperature and decreasing pH. Owing to the pH dependence of , it cannot be used to describe Mo fractionation in aqueous vapour-liquid systems with compositions different than those of this study. We have therefore calculated a composition-independent (Henry’s Law) constant, , for each experimental point, using the measured total Mo concentration in the vapour and the modeled concentration of H₂MoO₄⁰ in the liquid. This constant may be applied to aqueous vapour-liquid systems of known liquid composition to estimate the concentration of Mo in a vapour for which little chemical information is available, and thereby supplement the available fractionation data for natural porphyry-forming systems. The results of this study demonstrate that at conditions typical of natural porphyry ore-forming systems, a significant amount of molybdenum fractionates into the vapour over the liquid, and the vapour may transport quantities of Mo in excess of that in the liquid at pH conditions below those of the muscovite-microcline reaction boundary.     

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):

     

An XAS study of the structure and thermodynamics of Cu(I) chloride complexes in brines up to high temperature (400 °C, 600 bar)

The transport and deposition of copper in saline hydrothermal fluids are controlled by the stability of copper(I) complexes with ligands such as chloride. Despite their role in the formation of most hydrothermal copper deposits, the nature and stability of Cu(I) chloride complexes in highly saline brines remains controversial. We present new X-ray absorption data (P = 600 bar, T = 25-400 °C, salinity up to 17.2 m Cl), which indicate that the linear (x = 1, 2) complexes are stable up to supercritical conditions. Distorted trigonal planar complexes predominate at room temperature and at high salinity (>3 m LiCl): subtle changes in the XANES spectrum with increasing salinity may reflect geometric distortions of this complex. Similar changes were observed in UV-Vis data [Liu, W., Brugger, J., McPhail, D.C., Spiccia, L., 2002. A spectrophotometric study of aqueous copper(I) chloride complexes in LiCl solutions between 100 °C and 250 °C. Geochim. Cosmochim. Acta66, 3615-3633], and were erroneously interpreted as a new species, . Our XAS data and ab-initio XANES calculations show that this tetrahedral species is not present to any significant degree in our solutions. The stability of the complexe decreases with increasing temperature; under supercritical conditions and in brines under magmatic-hydrothermal conditions (e.g., 15.58 m Cl, 400 °C, 600 bar), only the linear Cu(I) chloride complexes were observed. This result and the instability of the complex are also consistent with the recent ab-initio molecular dynamic calculations of Sherman [Sherman D. M.(2007) Complexation of Cu⁺ in hydrothermal NaCl brines: ab-initio molecular dynamics and energetics. Geochim. Cosmochim. Acta71, 714-722]. This study illustrates the power of the quantitative nature of XANES and EXAFS measurements for deciphering the speciation of weak transition metal complexes up to magmatic-hydrothermal conditions.The systematic XANES data are used to retrieve the formation constant for at 150 °C, which is in good agreement with the reinterpretation of the UV-Vis data of Liu et al. (Liu et al., 2002). At high temperatures (?400 °C), the solubility of chalcopyrite in equilibrium with hematite-magnetite-pyrite and K-feldspar-muscovite-quartz calculated with the new properties is lower than that calculated using the previous model, and the calculated solubilities are at the lower end of the range of values measured in brine inclusions from porphyry copper systems.     

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

     

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:

     

Solubility measurements of synthetic neodymium monazite as a function of temperature at 2 kbars, and aqueous neodymium speciation in equilibrium with monazite

The solubility of synthetic NdPO₄ monazite end-member was experimentally determined from 300 up to 800 °C, at 2000 bars in pure water, and in aqueous chloride or phosphate solutions. Both the classical weight-loss method and a new method based on isotope dilution coupled with thermal ionization mass spectrometer were used. In the range of temperature studied monazite showed a prograde solubility from 10^−5.4 m at 300 °C up to 10^−2.57 m at 800 °C. Experiments in H₂O-H₃PO₄-NaCl-HCl solutions suggested Nd(OH)₃⁰ was the major species that was formed at high temperature and pressure. The equilibrium constants (log K) for the reaction:

     

Synthesis and properties of ammoniojarosites prepared with iron-oxidizing acidophilic microorganisms at 22-65 °C

Ammoniojarosite [(NH₄,H₃O)Fe₃(OH)₆(SO₄)₂], a poorly soluble basic ferric sulfate, was produced by microbiological oxidation of ferrous sulfate at pH 2.0-3.0 over a range of concentrations (5.4-805 mM) and temperatures (22-65 °C). Ammoniojarosites were also produced by chemical (abiotic) procedures in parallel thermal (36-95 °C) experiments. At 36 °C, schwertmannite [ideally Fe₈O₈(OH)₆(SO₄)] was the only solid product formed at <10 mM concentrations. Between 11.5 and 85.4 mM , a mixed product of ammoniojarosite and schwertmannite precipitated, as identified by X-ray diffraction. In excess of 165 mM , ammoniojarosite was the only solid phase produced. An increase in the incubation temperature using thermoacidophiles at 45 and 65 °C accelerated the formation of ammoniojarosite in culture solutions containing 165 mM . Both the biogenic and chemical ammoniojarosites were yellow (2Y-4Y in Munsell hue), low surface area (<1 m²/g), well crystalline materials with average c_o and a_o unit cell parameters of 17.467 ± 0.048 Å and 7.330 ± 0.006 Å, respectively. Strong positive correlations were observed between unit cell axial ratios (c_o/a_o) and increasing synthesis temperature in both biotic and abiotic systems. All samples were N deficient compared to stoichiometric ammoniojarosite, and both chemical and X-ray data indicated partial replacement of by H₃O⁺ to form solid solutions with 0.14-0.24 mole H₃O⁺ per formula unit. The morphology of the biogenic jarosites included aggregated discs, pseudo-cubic crystals and botryoidal particles, whereas the chemical specimens prepared at 36-95 °C were composed of irregular crystals with angular edges. Morphological information may thus be useful to evaluate environmental parameters and mode of formation. The data may also have application in predicting phase boundary conditions for Fe(III) precipitation in biogeochemical processes and treatment systems involving acid sulfate waters.     

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:

     

Prediction of the vapor-liquid distribution constants for volatile nonelectrolytes in water up to its critical temperature

The distribution of solutes between coexisting liquid and vapor phases of water can be expressed by the distribution constant, K_D, defined as

     

Dissolution kinetics of fosteritic olivine at 90-150 °C including effects of the presence of CO2

The steady state dissolution rate of San Carlos olivine [Mg_1.82Fe_0.18 SiO₄] in dilute aqueous solutions was measured at 90, 120, and 150 °C and pH ranging from 2 to 12.5. Dissolution experiments were performed in a stirred flow-through reactor, under either a nitrogen or carbon dioxide atmosphere at pressures between 15 and 180 bar. Low pH values were achieved either by adding HCl to the solution or by pressurising the reactor with CO₂, whereas high pH values were achieved by adding LiOH. Dissolution was stoichiometric for almost all experiments except for a brief start-up period. At all three temperatures, the dissolution rate decreases with increasing pH at acidic to neutral conditions with a slope of close to 0.5; by regressing all data for 2 ? pH ? 8.5 and 90 °C ? T ? 150 °C together, the following correlation for the dissolution rate in CO₂-free solutions is obtained:

     

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

     

Boric acid adsorption on humic acids: Ab initio calculation of structures, stabilities, B NMR and B,B isotopic fractionations of surface complexes

Boric acid, B(OH)₃, forms complexes in aqueous solution with a number of bidentate O-containing ligands, HL⁻, where H₂L is C₂O₄H₂ (oxalic acid), C₃O₄H₄ (malonic acid), C₂H₆O₂ (ethylene glycol), C₆H₆O₂ (catechol), C₁₀H₈O₂ (dioxynaphthalene) and C₂O₃H₄ (glycolic acid). McElligott and Byrne [McElligott, S., Byrne, R.H., 1998. Interaction of and in seawater: Formation of . Aquat. Geochem.3, 345-356.] have also found B(OH)₃ to form an aqueous complex with . Recently Lemarchand et al. [Lemarchand, E., Schott, J., Gaillardeet, J., 2005. Boron isotopic fractionation related to boron sorption on humic acid and the structure of surface complexes formed. Geochim. Cosmochim. Acta69, 3519-3533] have studied the formation of surface complexes of B(OH)₃ on humic acid, determining ¹¹B NMR shifts and fitted values of formation constants, and ¹¹B, ¹⁰B isotope fractionations for a number of surface complexation models. Their work helps to clarify both the nature of the interaction of boric acid with the functional groups in humic acid and the nature of some of these coordinating sites on the humic acid. The determination of isotope fractionations may be seen as a form of vibrational spectroscopy, using the fractionating element as a local probe of the vibrational spectrum. We have calculated quantum mechanically the structures, stabilities, vibrational spectra, ¹¹B NMR spectra and ¹¹B,¹⁰B isotope fractionations of a number of complexes B(OH)₂L⁻ formed by reactions of the type: