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.     

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.     

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

Oxidation of Sb(III) to Sb(V) by O2 and H2O2 in aqueous solutions

The rates of Sb(III) oxidation by O₂ and H₂O₂ were determined in homogeneous aqueous solutions. Above pH 10, the oxidation reaction of Sb(III) with O₂ was first order with respect to the Sb(III) concentration and inversely proportional to the H⁺ concentrations at a constant O₂ content of 0.22 × 10⁻³ M. Pseudo-first-order rate coefficients, k_obs, ranged from 3.5 × 10⁻⁸ s⁻¹ to 2.5 × 10⁻⁶ s⁻¹ at pH values between 10.9 and 12.9. The relationship between k_obs and pH was:

     

Modeling the effects of bond environment on equilibrium iron isotope fractionation in ferric aquo-chloro complexes

Four or five sets of ab initio models, including Unrestricted Hartree Fock (UHF) and hybrid Density Functional Theory (DFT) are calculated for each species in a series of aqueous ferric aquo-chloro complexes: , , , FeCl₃(H₂O)₃, FeCl₃(H₂O)₂, , FeCl₅H₂O²⁻, , ) in order to determine the relative isotopic fractionation among the complexes, to compare the results of different models for the same complexes, to examine factors that influence the magnitude of the isotopic fractionation, and to compare bond-partner-driven fractionation with redox-driven fractionation.Relative to , all models show a nearly linear decrease in ⁵⁶Fe/⁵⁴Fe as the number of Cl⁻ ions per Fe³⁺ ion increases, with slopes of −0.8‰ to −1.0‰ per Cl⁻ at 20 °C. At 20 °C, 1000 ln β (β = ⁵⁶Fe/⁵⁴Fe reduced partition function ratio relative to a dissociated Fe atom) values range from 8.93‰ to 9.73‰ for , 8.04-9.12‰ for , 7.61-8.73‰ for , 7.14-8.25‰ for , and 3.09-4.41‰ for . The fractionation between and ranges from 1.5‰ to 2.6‰, depending on the model; this is comparable in magnitude to fractionation effects due to Fe³⁺/Fe²⁺ redox reactions. β values from the UHF models are consistently higher than those from the hybrid DFT models.Isotopic fractionation is shown to be sensitive to differences in ligand bond stiffness (above), coordination number, bond length, and the frequency of the asymmetric Fe-X stretching vibrational mode, as predicted by previous theoretical studies. Complexes with smaller coordination numbers have higher 1000 ln β (7.46‰, 5.25‰, and 3.48‰ for , ,, respectively, from the B3LYP/6-31G(d) model). Species with the same number of chlorides but fewer waters also show the effect of coordination number on 1000 ln β: (7.46‰ vs. 7.05‰ for FeCl₃(H₂O)₂ vs. FeCl₃(H₂O)₃ and 5.25‰ vs. 4.94‰ for vs. FeCl₅H₂O²⁻ with the B3LYP/6-31G(d) model). As more Fe-Cl bonds substitute for Fe-OH₂ bonds (with a resulting decrease in β), the lengths of the Fe-Cl bonds and the Fe-O bonds increase.Preliminary modeling of shows an Fe³⁺/Fe²⁺ fractionation of 3.2‰ for the B3LYP/6-31G(d) model, in agreement with previous studies. The addition of an explicit outer hydration sphere of 12 H₂O molecules to models of improves agreement with measured vibrational frequencies and bond lengths; 1000 ln β increases by 0.8-1.0‰. An additional hydration sphere around increases 1000 ln β by only 0.1‰.Isotopic fractionations predicted for this simple system imply that ligands present in an aqueous iron environment are potentially important drivers of fractionation, and suggest that significant fractionation effects are likely in other aqueous systems containing sulfides or organic ligands. Fractionation effects due to both speciation and redox must be considered when interpreting iron isotope fractionations in the geological record.     

Solubility of NaCl in CO2 at high pressure and temperature: First experimental measurements

NaCl solubility in gaseous carbon dioxide has been measured in the pressure range from 30 to 70 MPa at 623 and 673 K. Our originally-designed high pressure apparatus allows in situ sampling of a portion of the fluid phase for chemical analysis. The results indicate that the solubility of NaCl increases with both temperature and pressure, and is about 4-5 orders of magnitude higher than saturated NaCl pressure values at the same temperature conditions (6.02 × 10⁻¹² at 623 K and 1.51 × 10⁻¹⁰ at 673 K). It is also 1-2 orders of magnitude greater than predictions according to the Equation of State of the ternary H₂O-CO₂-NaCl system by Duan, Moeller and Weare [Duan, Z., Moller, N., and Weare, J. H. (1995) Equation of state for the NaCl-H₂O-CO₂ system: prediction of phase equilibria and volumetric properties. Geochim. Cosmochim. Acta59, 2869] and has the opposite pressure dependence. The activity values of NaCl in the vapor phase, calculated from the experiments (with pure molten NaCl as a standard state in the vapor), have been fitted to the Darken Quadratic Formalism: , where, x_NaCl,v is mole the fraction of NaCl in the vapor phase, , , where P is the pressure in MPa and T the absolute temperature. Caution should be exerted while extrapolating this empirical equation far beyond the experimental P-T-compositional range.     

Solid solutions of trace Eu(III) in calcite: Thermodynamic evaluation of experimental data over a wide range of pH and pCO2

The thermodynamics of dilute Eu-calcite solid solutions formed under widely different pH-pCO₂ conditions at T = 25°C and p = 1 bar were investigated using three sets of Eu(III) uptake experiments, two of which were taken from the literature: (a) recrystallization in synthetic cement pore water at pH ∼ 13 and pCO₂ ∼ 10⁻¹³ bar (this work); (b) coprecipitation in 0.1 M NaClO₄ at pH ∼ 6 and pCO₂ ∼ 1 bar; (c) coprecipitation in synthetic seawater at pH ∼ 8 and pCO₂ ranging from 3 × 10⁻⁴ to 0.3 bar.Solid solution formation was modeled using the Gibbs energy minimization (GEM) method. In a first step (“forward” modeling), we tested ideal binary solid solution models between calcite and the Eu end-members Eu₂(CO₃)_3, EuNa(CO₃)₂, Eu(OH)CO₃ or Eu(OH)₃, for which solids with independently measured solubility products exist. None of these four binary solid solutions was capable of reproducing all three experimental datasets simultaneously. In a second step (“inverse” modeling), ideal binary solid solutions were constructed between calcite and the candidate Eu end-members EuO(OH), EuH(CO₃)₂ and EuO(CO₃)_0.5, for which no independent solubility products are available. For each single data point and each of these end-members, a free energy of formation with inherent activity coefficient term ( = G_α^o + RT lnγ_α) was estimated from “dual thermodynamic” GEM calculations. The statistical mean of was then calculated for each of the three datasets. A specific end-member was considered to be acceptable if a standard deviation of ± 2 kJ mol⁻¹ or less resulted for each single dataset, and if the mean -values calculated for the three datasets coincided. No binary solid solution with any of the seven above mentioned end-members proved to satisfy these criteria.The third step in our analysis involved consideration of ternary solid solutions with CaCO₃ as the major end-member and any two of the seven considered Eu trace end-members. It was found that the three datasets can only be reproduced simultaneously with the ternary ideal solid solution EuH(CO₃)₂ - EuO(OH) - CaCO₃, setting = −1773 kJ mol⁻¹ and = −955 kJ mol⁻¹, whereas all other end-member combinations failed. Our results are consistent with time-resolved laser fluorescence data for Cm(III) and Eu(III) indicating that two distinct species are incorporated in calcite: one partially hydrated, the other completely dehydrated. In conclusion, our study shows that substitution of trivalent for divalent cations in carbonate crystal structures is a more complex process than the classical isomorphic divalent-divalent substitution and may need consideration of multicomponent solid solution models.     

The solubility of platinum and gold in NaCl brines at 1.5 kbar, 600 to 800°C: A laser ablation ICP-MS pilot study of synthetic fluid inclusions

The concentration and distribution of Pt and Au in a fluid-melt system has been investigated by reacting the metals with S-free, single-phase aqueous brines (20, 50, 70 wt% eq. NaCl) ± peraluminous melt at a confining pressure of 1.5 kbar and temperatures of 600 to 800 °C, trapping the fluid in synthetic fluid inclusions (quartz-hosted) and vesicles (silicate melt-hosted), and quantifying the metal content of the trapped fluid and glass by laser ablation ICP-MS. HCl concentration was buffered using the assemblage albite-andalusite-quartz and fO₂ was buffered using the assemblage Ni-NiO. Over the range of experimental conditions, measured concentrations of Pt and Au in the brines (, ) are on on the order of 1-10³ ppm. Concentrations of Pt and Au in the melt (, ) are ∼35-100 ppb and ∼400-1200 ppb, respectively. Nernst partition coefficients (, ) are on the order of 10²-10³ and vary as a function of (non-Henry’s Law behavior). Trapped fluids show a significant range of metal concentrations within populations of inclusions from single experiments (∼ 1 log unit variability for Au; ∼2-3 log unit variability for Pt). Variability in metal concentration within single inclusion groups is attributed to premature brine entrapment (prior to metal-fluid-melt equilibrium being reached); this allows us to make only minimum estimates of metal solubility using metal concentrations from primary inclusions. The data show two trends: (i) maximum and average values of and in inclusions decrease ∼2 orders of magnitude as fluid salinity () increases from ∼4 to 40 molal (20 to 70 wt % eq. NaCl) at a constant temperature; (ii) maximum and average values of increase approximately 1 order of magnitude for every 100°C increase temperature at a fixed . The observed behavior may be described by the general expression:

     

Microbial ammonia oxidation and enhanced nitrogen cycling in the Endeavour hydrothermal plume

Ammonium was injected from the subseafloor hydrothermal system at the Endeavour Segment, Juan de Fuca Ridge, into the deep-sea water column resulting in an -rich (?177 nM) neutrally buoyant hydrothermal plume. This was quickly removed by both autotrophic ammonia oxidation and assimilation. The former accounted for at least 93% of total net removal, with its maximum rate in the neutrally buoyant plume (?53 nM d⁻¹) up to 10-fold that in background deep water. Ammonia oxidation in this plume potentially added 26-130 mg into the deep-sea water column. This oxidation process was heavily influenced by the presence of organic-rich particles, with which ammonia-oxidizing bacteria (AOB) were often associated (40-68%). AOB contributed up to 10.8% of the total microbial communities within the plume, and might constitute a novel lineage of β-proteobacterial AOB based on 16S rRNA and amoA phylogenetic analyses. Meanwhile, assimilation rates were also substantially enhanced within the neutrally buoyant plume (?26.4 nM d⁻¹) and accounted for at least 47% of total net removal rates. The combined oxidation and assimilation rates always exceeded total net removal rates, suggesting active in situregeneration rates of at least an order of magnitude greater than the particulate nitrogen flux from the euphotic zone. Ammonia oxidation is responsible for turnover of 0.7-13 days and is probably the predominant in situ organic carbon production process (0.6-13 mg C m⁻² d⁻¹) at early stages of Endeavour neutrally buoyant plumes.     

Thermodynamic properties of anhydrous smectite MX-80, illite IMt-2 and mixed-layer illite-smectite ISCz-1 as determined by calorimetric methods. Part I: Heat capacities, heat contents and entropies

The heat capacities of the anhydrous international reference clay minerals, smectite MX-80, illite IMt-2 and mixed-layer illite-smectite ISCz-1, were measured by low temperature adiabatic calorimetry and differential scanning calorimetry, from 6 to 520 K (at 1 bar). The samples were chemically purified and Na-saturated. Dehydrated clay fractions <2 μm were studied. The structural formulae of the corresponding clay minerals, obtained after subtracting the remaining impurities, are K_0.026Na_0.435Ca_0.010(Si_3.612Al_0.388) (Al_1.593Mg_0.228Ti_0.011)O₁₀(OH)₂ for smectite MX-80, K_0.762Na_0.044(Si_3.387Al_0.613) (Al_1.427Mg_0.241O₁₀(OH)₂ for illite IMt-2 and K_0.530Na_0.135(Si_3.565Al_0.435)(Al_1.709Mg_0.218Ti_0.005)O₁₀(OH)₂for mixed-layer ISCz-1. From the heat capacity values, we determined the molar entropies, standard entropies of formation and heat contents of these minerals. The following values were obtained at 298.15 K and 1 bar:

	(J mol−1 K−1)	S0 (J mol−1 K−1)
Smectite MX-80	326.13 ± 0.10	280.56 ± 0.16
Illite IMt-2	328.21 ± 0.10	295.05 ± 0.17
Mixed-layer ISCz-1	320.79 ± 0.10	281.62 ± 0.15

Full-size table

     

Sorption of yttrium and rare earth elements by amorphous ferric hydroxide: Influence of solution complexation with carbonate

The influence of solution complexation on the sorption of yttrium and the rare earth elements (YREEs) by amorphous ferric hydroxide was investigated at 25 °C over a range of pH (4.0-7.1) and carbonate concentrations . Distribution coefficients, defined as , where [MS_i]_T is the total concentration of sorbed YREE, M_T is the total YREE concentration in solution, and [S_i] is the concentration of amorphous ferric hydroxide, initially increased in magnitude with increasing carbonate concentration, and then decreased. The initial increase of is due to sorption of YREE carbonate complexes , in addition to sorption of free YREE ions (M³⁺). The subsequent decrease of , which is more extensive for the heavy REEs, is due to the increasing intensity of YREE solution complexation by carbonate ions. The competition for YREEs between solution complexation and surface complexation was modeled via the equation:

     

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:

     

Acoustic velocity measurements on Na2O-TiO2-SiO2 liquids: Evidence for a highly compressible TiO2 component related to five-coordinated Ti

Longitudinal acoustic velocities were measured at 1 bar in 10 Na₂O-TiO₂-SiO₂ (NTS) liquids for which previous density and thermal expansion data are reported in the literature. Data were collected with a frequency-sweep acoustic interferometer at centered frequencies of 4.5, 5, and 6 MHz between 1233 and 1896 K; in all cases, the sound speeds decrease with increasing temperature. Six of the liquids have a similar TiO₂ concentration (∼25 mol %), so that the effect of varying Na/Si ratio on the partial molar compressibility of the TiO₂ component can be evaluated. Theoretically based models for β_T and (∂V/∂P)_T as a function of composition and temperature are presented. As found previously for the partial molar volume of TiO₂ in sodium silicate melts, values of (13.7-18.8 × 10⁻²/GPa) vary systematically with the Na/Si and Na/(Si + Ti) ratio in the liquid. In contrast values of for the SiO₂ and Na₂O components (6.6 and 8.0 × 10⁻²/GPa, respectively, at 1573 K) are independent of composition. Na₂O is the only component that contributes to the temperature dependence of the compressibility of NTS liquids (1.13 ± 0.04 × 10⁻⁴/GPa K). The results further indicate that the TiO₂ component is twice as compressible as the Na₂O and SiO₂ components. The enhanced compressibility of TiO₂ appears to be related to the abundance of five-coordinated Ti (^[5]Ti) in these liquids, but not with a change in Ti coordination. Instead, it is proposed that the asymmetric geometry of ^[5]Ti in a square pyramidal site promotes different topological rearrangements in alkali titanosilicate liquids, which lead to the enhanced compressibility of TiO₂.