The effects of organic acids on the dissolution of silicate minerals: A case study of oxalate catalysis of kaolinite dissolution

Most studies agree that the dissolution rate of aluminosilicates in the presence of oxalic and other simple carboxylic acids is faster than the rate with non-organic acid under the same pH. However, the mechanisms by which organic ligands enhance the dissolution of minerals are in debate. The main goal of this paper was to study the mechanism that controls the dissolution rate of kaolinite in the presence of oxalate under far from equilibrium conditions (−29 < ΔG_r < −18 kcal mol⁻¹). Two types of experiments were performed: non-stirred flow-through dissolution experiments and batch type adsorption isotherms. All the experiments were conducted at pH 2.5-3.5 in a thermostatic water-bath held at a constant temperature of 25.0, 50.0 or 70.0 ± 0.1 °C. Kaolinite dissolution rates were obtained based on the release of silicon and aluminum at steady state. The results show good agreement between these two estimates of kaolinite dissolution rate. At constant temperature, there is a general trend of increase in the overall dissolution rate as a function of the total concentration of oxalate in solution. The overall kaolinite dissolution rates in the presence of oxalate was up to 30 times faster than the dissolution rate of kaolinite at the same temperature and pH without oxalate as was observed in our previous study. Therefore, these rate differences are related to differences in oxalate and aluminum concentrations. Within the experimental variability, the oxalate adsorption at 25, 50, and 70 °C showed the same dependence on the sum of the activities of oxalate and bioxalate in solution. The change of oxalate concentration on the kaolinite surface (C_s,ox) as a function of the sum of the activities of the oxalate and bioxalate in solution may be described by the general adsorption isotherm:

     

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⁻¹).     

Kinetics and mechanism of polythionate oxidation to sulfate at low pH by O2 and Fe

Polythionates (S_xO₆²⁻) are important in redox transformations involving many sulfur compounds. Here we investigate the oxidation kinetics and mechanisms of trithionate and tetrathionate oxidation between pH 0.4 and pH 2 in the presence of Fe³⁺ and/or oxygen. In these solutions, Fe³⁺ plus oxygen oxidizes tetrathionate and trithionate at least an order of magnitude faster than oxygen alone. Kinetic measurements, coupled with density functional calculations, suggest that the rate-limiting step for tetrathionate oxidation involves Fe³⁺ attachment, followed by electron density shifts that result in formation of a sulfite radical and S₃O₃⁰ derivatives. The overall reaction orders for trithionate and tetrathionate are fractional due to rearrangement reactions and side reactions between reactants and intermediate products. The pseudo-first order rate coefficients for tetrathionate range from 10⁻¹¹ s⁻¹ at 25°C to 10⁻⁸ s⁻¹ at 70°C, compared to 2 × 10⁻⁷ s⁻¹ at 35 °C for trithionate. The apparent activation energy (E_A) for tetrathionate oxidation at pH 1.5 is 104.5 ± 4.13 kJ/mol. A rate law at pH 1.5 and 70°C between 0.5 and 5 millimolar [Fe³⁺] is of the form:

     

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:

     

Kaolinite dissolution and precipitation kinetics at 22 °C and pH 4

Dissolution and precipitation rates of low defect Georgia kaolinite (KGa-1b) as a function of Gibbs free energy of reaction (or reaction affinity) were measured at 22 °C and pH 4 in continuously stirred flowthrough reactors. Steady state dissolution experiments showed slightly incongruent dissolution, with a Si/Al ratio of about 1.12 that is attributed to the re-adsorption of Al on to the kaolinite surface. No inhibition of the kaolinite dissolution rate was apparent when dissolved aluminum was varied from 0 and 60 μM. The relationship between dissolution rates and the reaction affinity can be described well by a Transition State Theory (TST) rate formulation with a Temkin coefficient of 2

     

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 .     

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.     

Oxygen isotope fractionation between apatite-bound carbonate and water determined from controlled experiments with synthetic apatites precipitated at 10-37 °C

The oxygen isotope fractionation between the structural carbonate of inorganically precipitated hydroxyapatite (HAP) and water was determined in the range 10-37 °C. Values of 1000 ln α() are linearly correlated with inverse temperature (K) according to the following equation: 1000 ln α() = 25.19 (±0.53)·T⁻¹ − 56.47 (±1.81) (R² = 0.998). This fractionation equation has a slightly steeper slope than those already established between calcite and water ( [O’Neil et al., 1969] and [Kim and O’Neil, 1997]) even though measured fractionations are of comparable amplitude in the temperature range of these experimental studies. It is consequently observed that the oxygen isotope fractionation between apatite carbonate and phosphate increases from about 7.5‰ up to 9.1‰ with decreasing temperature from 37 °C to 10 °C. A compilation of δ¹⁸O values of both phosphate and carbonate from modern mammal teeth and bones confirms that both variables are linearly correlated, despite a significant scattering up to 3.5‰, with a slope close to 1 and an intercept corresponding to a 1000 ln α() value of 8.1‰. This apparent fractionation factor is slightly higher or close to the fractionation factor expected to be in the range 7-8‰ at the body temperature of mammals.     

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

     

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.     

Carbon solubility in core melts in a shallow magma ocean environment and distribution of carbon between the Earth’s core and the mantle

The solubility of carbon in Fe and Fe-5.2 wt.% Ni melts, saturated with graphite, determined by electron microprobe analysis of quenched metal melts was 5.8 ± 0.1 wt.% at 2000 °C, 6.7 ± 0.2 wt.% at 2200 °C, and 7.4 ± 0.2 wt.% at 2410 °C at 2 GPa, conditions relevant for core/mantle differentiation in a shallow magma ocean. These solubilities are slightly lower than low-pressure literature values and significantly beneath calculated values for even higher pressures [e.g., Wood B. J. (1993) Carbon in the core. Earth Planet. Sci. Lett.117, 593-607]. The trend of C solubility versus temperature for Fe-5.2 wt.% Ni melt, within analytical uncertainties, is similar to or slightly lower (∼0.2-0.4 wt.%) than that of pure Fe. Carbon content of core melts and residual mantle silicates derived from equilibrium batch or fractional segregation of core liquids and their comparison with our solubility data and carbon content estimate of the present day mantle, respectively, constrain the partition coefficient of carbon between silicate and metallic melts, in a magma ocean. For the entire range of possible bulk Earth carbon content from chondritic to subchondritic values, of 10⁻⁴ to 1 is derived. But for ∼1000 ppm bulk Earth carbon, is between 10⁻² and 1. Using the complete range of possible for a magma ocean at ∼2200 °C, we predict maximum carbon content of the Earth’s core to be ∼6-7 wt.% and a preferred value of 0.25 ± 0.15 wt.% for a bulk Earth carbon concentration of ∼1000 ppm.     

Experimental study of brucite dissolution and precipitation in aqueous solutions: surface speciation and chemical affinity control

Dissolution and precipitation rates of brucite (Mg(OH)₂) were measured at 25°C in a mixed-flow reactor as a function of pH (2.5 to 12), ionic strength (10⁻⁴ to 3 M), saturation index (−12 < log Ω < 0.4) and aqueous magnesium concentrations (10⁻⁶ to 5·10⁻⁴ M). Brucite surface charge and isoelectric point (pH_IEP) were determined by surface titrations in a limited residence time reactor and electrophoretic measurements, respectively. The pH of zero charge and pH_IEP were close to 11. A two-pK, one site surface speciation model which assumes a constant capacitance of the electric double layer (5 F/m²) and lack of dependence on ionic strength predicts the dominance of >MgOH₂⁺ species at pH < 8 and their progressive replacement by >MgOH° and >MgO⁻ as pH increases to 10-12. Rates are proportional to the square of >MgOH₂⁺ surface concentration at pH from 2.5 to 12. In accord with surface speciation predictions, dissolution rates do not depend on ionic strength at pH 6.5 to 11. Brucite dissolution and precipitation rates at close to equilibrium conditions obeyed TST-derived rate laws. At constant saturation indices, brucite precipitation rates were proportional to the square of >MgOH₂⁺ concentration. The following rate equation, consistent with transition state theory, describes brucite dissolution and precipitation kinetics over a wide range of solution composition and chemical affinity:

     

Surface properties, solubility and dissolution kinetics of bamboo phytoliths

Although phytoliths, constituted mainly by micrometric opal, exhibit an important control on silicon cycle in superficial continental environments, their thermodynamic properties and reactivity in aqueous solution are still poorly known. In this work, we determined the solubility and dissolution rates of bamboo phytoliths collected in the Réunion Island and characterized their surface properties via electrophoretic measurements and potentiometric titrations in a wide range of pH. The solubility product of “soil” phytoliths ( at 25 °C) is equal to that of vitreous silica and is 17 times higher than that of quartz. Similarly, the enthalpy of phytoliths dissolution reaction is close to that of amorphous silica but is significantly lower than the enthalpy of quartz dissolution. Electrophoretic measurements yield isoelectric point pH_IEP = 1.2 ± 0.1 and 2.5 ± 0.2 for “soil” (native) and “heated” (450 °C heating to remove organic matter) phytoliths, respectively. Surface acid-base titrations allowed generation of a 2-pK surface complexation model. Phytoliths dissolution rates, measured in mixed-flow reactors at far from equilibrium conditions at 2 ? pH ? 12, were found to be intermediate between those of quartz and vitreous silica. The dissolution rate dependence on pH was modeled within the concept of surface coordination theory using the equation:

     

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:

     

The kinetics of iodide oxidation by the manganese oxide mineral birnessite

The kinetics of iodide (I⁻) and molecular iodine (I₂) oxidation by the manganese oxide mineral birnessite (δ-MnO₂) was investigated over the pH range 4.5-6.25. I⁻ oxidation to iodate proceeded as a two-step reaction through an I₂ intermediate. The rate of the reaction varied with both pH and birnessite concentration, with faster oxidation occurring at lower pH and higher birnessite concentration. The disappearance of I⁻ from solution was first order with respect to I⁻ concentration, pH, and birnessite concentration, such that −d[I⁻]/dt = k[I⁻][H⁺][MnO₂], where k, the third order rate constant, is equal to 1.08 ± 0.06 × 10⁷ M⁻² h⁻¹. The data are consistent with the formation of an inner sphere I⁻ surface complex as the first step of the reaction, and the adsorption of I⁻ exhibited significant pH dependence. Both I₂, and to a lesser extent, sorbed to birnessite. The results indicate that iodine transport in mildly acidic groundwater systems may not be conservative. Because of the higher adsorption of the oxidized I species I₂ and , as well as the biophilic nature of I₂, redox transformations of iodine must be taken into account when predicting I transport in aquifers and watersheds.     

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

     

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.     

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: