Linear free energy relationships for metal-ligand complexation: Monodentate binding to negatively-charged oxygen donor atoms

The analysis presented below suggests that the following equation models monodentate binding to negatively-charged oxygen donor atoms with no distinction between phenolic, carboxylic or inorganic hydroxide functional groups: log K_ML = α_O log K_HL + β_O; where K_ML is the metal-ligand formation constant, K_HL is the corresponding proton-ligand formation constant, and α_O and β_O are termed the Irving-Rossotti slope and intercept, respectively. Linear free energy relationships (LFERs) of this type are presented for 24 different metal ions complexing to negatively-charged oxygen donor atoms. Ligands selected for use in LFERs meet the following criteria: (i) they contain negatively-charged oxygen donor atoms (e.g., carboxylic acids and phenols), (ii) they are capable of only monodentate binding to metal ions, (iii) steric hindrances are not expected to influence the extent of metal-ligand binding, and (iv) the negatively-charged oxygen donor atom is the only functional group that imparts charge. The intercept of all LFERs was nearly zero for all metal ions investigated (β_O ≈ 0). The magnitude of α_O indicates the relative preference of metal binding to negatively-charged oxygen donor atoms and to the proton. Values of α_O can be used in QSARs (quantitative structure activity relationships) to estimate metal-NOM (natural organic matter) binding constants employed in the Windermere Humic Aqueous Model (WHAM) version V.     

Semi-empirical proton binding constants for natural organic matter

Average proton binding constants (K_H,i) for structure models of humic (HA) and fulvic (FA) acids were estimated semi-empirically by breaking down the macromolecules into reactive structural units (RSUs), and calculating K_H,i values of the RSUs using linear free energy relationships (LFER) of Hammett. Predicted log K_H,COOH and log K_H,Ph-OH are 3.73 ± 0.13 and 9.83 ± 0.23 for HA, and 3.80 ± 0.20 and 9.87 ± 0.31 for FA. The predicted constants for phenolic-type sites (Ph-OH) are generally higher than those derived from potentiometric titrations, but the difference may not be significant in view of the considerable uncertainty of the acidity constants determined from acid-base measurements at high pH. The predicted constants for carboxylic-type sites agree well with titration data analyzed with Model VI (4.10 ± 0.16 for HA, 3.20 ± 0.13 for FA; Tipping, 1998), the Impermeable Sphere model (3.50-4.50 for HA; Avena et al., 1999), and the Stockholm Humic Model (4.10 ± 0.20 for HA, 3.50 ± 0.40 for FA; Gustafsson, 2001), but differ by about one log unit from those obtained by Milne et al. (2001) with the NICA-Donnan model (3.09 ± 0.51 for HA, 2.65 ± 0.43 for FA), and used to derive recommended generic values. To clarify this ambiguity, 10 high-quality titration data from Milne et al. (2001) were re-analyzed with the new predicted equilibrium constants. The data are described equally well with the previous and new sets of values (R² ? 0.98), not necessarily because the NICA-Donnan model is overparametrized, but because titration lacks the sensitivity needed to quantify the full binding properties of humic substances. Correlations between NICA-Donnan parameters are discussed, but general progress is impeded by the unknown number of independent parameters that can be varied during regression of a model fit to titration data. The high consistency between predicted and experimental K_H,COOH values, excluding those of Milne et al. (2001), gives faith in the proposed semi-empirical structural approach, and its usefulness to assess the plausibility of proton stability constants derived from simulations of titration data.     

Sorption modelling on illite. Part II: Actinide sorption and linear free energy relationships

Sorption edge data for Ni(II), Co(II), Eu(III) and Sn(IV) [Bradbury M. H. and Baeyens B. (2009) Sorption modelling on illite. Part I: titration measurements and sorption of Ni(II), Co(II), Eu(III) and Sn(IV), Part I] on purified Na-Illite du Puy are available from some previous work, and some new measurements for Am(III), Th(IV), Pa(V) and U(VI) are presented here. All of these sorption edge measurements have been modelled with a 2 site protolysis non-electrostatic surface complexation and cation exchange (2SPNE SC/CE) sorption model for which the site types, site capacities and protolysis constants were fixed [Bradbury M. H. and Baeyens B. (2009), Part I]. In addition, two further data sets for the sorption of Am(III) and Np(V) on Illite du Puy, obtained from the literature, were also modelled in this work. Thus, surface complexation constants for the strong sites in the 2SPNE SC/CE sorption model for nine metals with valence states from II to VI have been obtained. A linear relationship between the logarithm of strong site metal binding constants, ^SK_x−1, and the logarithm of the corresponding aqueous hydrolysis stability constant, ^OHK_x, extending over nearly 35 orders of magnitude is established here for illite for these nine metals. Such correlations are often termed linear free energy relationships (LFER), and although they are quite common in aqueous phase chemistry, they are much less so in surface chemistry, especially over this large range. The LFER for illite could be described by the equation: where, “x” is an integer. A similar relationship has been previously obtained for montmorillonite, thus LFERs relating to the sorption on two of the most important clay minerals present in natural systems have been established. Such an LFER approach is an extremely useful tool for estimating surface complexation constants for metals in a chemically consistent manner. It provides a means of obtaining sorption values for radionuclides for which there are no measured values and thus allows gaps in missing sorption data to be filled. An ultimate goal of this approach is to develop a thermodynamic sorption database. This could then be used in radioactive waste management performance assessment studies to calculate sorption in natural systems, and thereby replace the current usage of single solid liquid distribution coefficients (K_d values) to describe radionuclide uptake. Finally, with the data now available, the 2SPNE SC/CE sorption model can be ported into reactive transport models allowing radionuclide migration to be calculated under spatially and temporally changing conditions.     

Lanthanide-humic substances complexation. I. Experimental evidence for a lanthanide contraction effect

The interaction of the lanthanides (Ln) with humic substances (HS) was investigated with a novel chemical speciation tool, Capillary Electrophoresis-Inductively Coupled Plasma Mass Spectrometry (CE-ICP-MS). By using an EDTA-ligand competition method, a bi-modal species distribution of LnEDTA and LnHS is attained, separated by CE, and detected online by sector field ICP-MS. We quantified the binding of all 14 rare earth elements (REEs), Sc and Y with Suwannee river fulvic acid, Leonardite coal humic acid, and Elliot soil humic acid under environmental conditions (pH 6-9, 0.001-0.1 mol L⁻¹ NaNO₃, 1-1000 nmol L⁻¹ Ln, 10-20 mg L⁻¹ HS). Conditional binding constants for REE-HS interaction (K_c) ranged from 8.9 < log K_c < 16.5 under all experimental conditions, and display a lanthanide contraction effect, Δ_LK_c: a gradual increase in K_c from La to Lu by 2-3 orders of magnitude as a function of decreasing ionic radius. HS polyelectrolyte effects cause K_c to increase with increasing pH and decreasing ionic strength. Δ_LK_c increases significantly with increasing pH, and likely with decreasing ionic strength. Based on a strong correlation between Δ_LK_c values and denticity for organic acids, we suggest that HS form a range of tri- to tetra-dentate complexes under environmental conditions. These results confirm HS to be a strong complexing agent for Ln, and show rigorous experimental evidence for potential REE fractionation by HS complexation.     

A critical review of three methods used for the measurement of mercury (Hg)-dissolved organic matter stability constants

Three experimental techniques – ion exchange, liquid–liquid extraction with competitive ligand exchange, and solid-phase extraction with competitive ligand exchange (CLE–SPE) – were evaluated as methods for determining conditional stability constants (K) for the binding of mercury (Hg²⁺) to dissolved organic matter (DOM). To determine the utility of a given method to measure stability constants at environmentally relevant experimental conditions, experimental results should meet three criteria: (1) the data must be experimentally valid, in that they were acquired under conditions that meet all the requirements of the experimental method, (2) the Hg:DOM ratio should be determined and it should fall within levels that are consistent with environmental conditions, and (3) the stability constants must fall within the detection window of the method. The ion exchange method was found to be limited by its detection window, which constrains the method to stability constants with log K values less than about 14. The liquid–liquid extraction method was found to be complicated by the ability of Hg–DOM complexes to partition into the organic phase. The CLE–SPE method was found to be the most suitable of these methods for the measurement of Hg–DOM stability constants. Stability constants for DOM isolates measured using the CLE–SPE method at environmentally relevant Hg:DOM ratios were log K = 25–30 (M⁻¹). These values are consistent with the strong Hg²⁺ binding expected for reduced S-containing binding sites.     

Metal-silicate partitioning of cesium: Implications for core formation

Here we present the first set of metal-silicate partitioning data for Cs, which we use to examine whether the primitive mantle depletion of Cs can be attributed to core segregation. Our experiments independently varied pressure from 5 to 15 GPa, temperature from 1900 to 2400 °C, metallic sulfur content from pure Fe to pure FeS, silicate melt polymerization, expressed as a ratio of non-bridging oxygens to tetrahedrally coordinated cations (nbo/t) from 1.26 to 3.1, and fO₂ from two to four log units below the iron-wüstite buffer. The most important controls on the partitioning behavior of alkalis were the metallic sulfur content, expressed as X_S, and the nbo/t of the silicate liquid. Normalization of X_S to 0.5 yielded the following expressions for D-values as a function of nbo/t: log D_Na = −2.0 + 0.44 × (nbo/t), log D_K = −2.4 + 0.67 × ( nbo/t), and log D_Cs = −3.2 + 1.17 × (nbo/t). Normalization of nbo/t to 2.7 resulted in the following equations for D-values as a function of S content: log D_Na = −4.1 + 6.4 × X_S, log D_K = −7.7 + 13.9 × X_S, and log D_Cs = −12.1 + 23.3 × X_S.There appears to be a negative pressure effect up to 15 GPa, but it should be noted that this trend was not present before normalization, and is based on only two measurements. There is a positive trend in cesium’s metal-silicate partition coefficient with increasing temperature. D_Cs exhibits the largest change and increased by a factor of three over 500 °C. The effect of oxygen fugacity has not been precisely determined but in general, lowering fO₂ by two log units resulted in a rise in all D-values of approximately an order of magnitude. In general, the sensitivity of partition coefficients to changing parameters increased with atomic number.The highest D-value for Cs observed in this study is 0.345, which was obtained at nbo/t of 2.7 and a metal phase of pure FeS. This metallic composition has far more S than has been suggested for any credible core-forming metal. We therefore conclude that the depletion of Cs in Earth’s mantle is either caused by radically different behavior of Cs at pressures higher than 15 GPa or is not related to core formation. Even so, we have shown that a planet with a sufficient S inventory may incorporate significant amounts of alkali elements into its core.     

Surface complexation modelling of Cd(II), Cu(II), Ni(II), Pb(II) and Zn(II) adsorption onto kaolinite

Proton binding constants for the edge and basal surface sites of kaolinite were determined by batch titration experiments at 25 °C in the presence of 0.1 M, 0.01 M and 0.001 M solutions of NaNO₃ and in the pH range 3-9. By optimizing the results of the titration experiments, the ratio of the edge sites to the basal surface sites was found to be 6:1. The adsorption of Cd(II), Cu(II), Ni(II), Zn(II) and Pb(II) onto kaolinite suspensions was investigated using batch adsorption experiments and results suggested that in the lower pH range the metallic cations were bound through non-specific ion exchange reactions on the permanently charged basal surface sites (X⁻). Adsorption on these sites was greatly affected by ionic strength. With increasing pH, the variable charged edge sites (SOH) became the major adsorption sites and inner-sphere specifically adsorbed monodentate complexes were believed to be formed. The effect of ionic strength on the extent of adsorption of the metals on the variable charged edge sites was much less than those on the permanently charged sites. Two binding constants, log K_(X2Me) and log K_(SOMe), were calculated by optimizing these constants in the computer program FITEQL. A model combining non-specific ion exchange reactions and inner-sphere specific surface complexations was developed to predict the adsorption of heavy metals onto kaolinite in the studied pH range. Linear free energy relationships were found between the edge site binding constants and the first hydrolysis constants of the metals.     

Solubility of plutonium(VI) carbonate in saline solutions

Among the plutonium oxidation states found to form in the environment, mobile plutonium(VI) can exist under oxidizing conditions and in waters with high chloride content due to radiolysis effects. We are investigating the solubility and speciation of plutonium(VI) carbonate under conditions relevant to natural waters and brines such as those found near some geologic radioactive waste repositories. The solid Pu(VI) phase PuO₂CO₃(s) was prepared and its solubility was measured in NaCl and NaClO₄ solutions in a CO₂ atmosphere as a function of pH and ionic strength (0.1-5.6 m). The concentration of soluble plutonium in solution was calculated from spectroscopic data and liquid scintillation counting. Spectroscopic measurements also revealed the plutonium oxidation state. The apparent solubility product of PuO₂CO₃(s) was determined at selected electrolyte concentrations to be, log K_s,0 = −13.95 ± 0.07 (0.1 m NaCl), log K_s,0 = −14.07 ± 0.13 (5.6 m NaCl), and log K_s,0 = −15.26 ± 0.11 (5.6 m NaClO₄). Specific ion interaction theory was used to calculate the solubility product at zero ionic strength, .     

Adsorption and surface complexation modeling of palladium,rhodium and platinum in surficial semi-arid soils and sediment

Sorption isotherms, time-dependent adsorption and surface complexation modeling studies were used to investigate the post-depositional mobility of three of the platinum group-elements (Pd, Rh, and Pt) in semi-arid soil and sediment samples with varying surface properties. The acidity constants (Log K_a1 and Log K_a2), optimized from batch titration data, ranged from 4.69 to 5.34 for Log K_a1 and from −6.51 to −7.61 for Log K_a2, suggesting the occurrence of both protonation and deprotonation reactions on the solid surfaces. Partition coefficients and removal rates of the metals had a general trend of Pd > Pt > Rh. The sediment sample, with the highest clay content and exchangeable cation concentrations, also had the highest affinity for the metals. The times required for sediment to adsorb 63% of the metals were 2.63 h, 4.08 h and 10.64 h for Pd, Pt and Rh, respectively. The FITEQL program successfully optimized the conditional binding constants of the metals on the solids from batch adsorption data. The constants decreased in the order of Pd > Rh > Pt, which was consistent with the observed high affinity of the solids for Pd. The modeling results also showed that aqueous Pd was the least sensitive to pH followed by Rh and Pt. However, metal adsorption below the points of zero net proton charges (ca. pH 6.7) is attributable to the involvement of permanent negatively charged binding sites in the adsorption process. Notably, partition coefficients, removal rates and conditional binding constants all showed a high affinity of Pd for the solids. A similarity between the model outputs and the batch adsorption data indicates the suitability of the model for describing the mobility and retention of the three metals in semi-arid soils and sediments.     

Stability of YREE complexes with the trihydroxamate siderophore desferrioxamine B at seawater ionic strength

Organic complexation of yttrium and the rare earth elements (YREEs), although generally believed to be important, is an understudied aspect of YREE solution speciation in the open ocean. We report the first series of stability constants for complexes of YREEs (except Ce and Pm) with the trihydroxamate siderophore desferrioxamine B (DFOB), representing a class of small organic ligands that have an extraordinary selectivity for Fe(III) and are found in surface seawater at low-picomolar concentrations. Constants were measured by potentiometric titration of DFOB (pH 3-10) in the presence of single YREEs, in simple media at seawater ionic strength (NaClO₄ or NaCl, I = 0.7 M). Under these circumstances, the terminal amine of DFOB does not deprotonate. The four acid dissociation constants of the siderophore were determined separately by potentiometric titration of DFOB alone.Values for the bidentate (log β₁), tetradentate (log β₂), and hexadentate (log β₃) complexes of La-Lu range from 4.88 to 6.53, 7.70 to 11.27, and 10.09 to 15.19, respectively, while Y falls between Gd and Tb in each case. Linear free-energy relations of the three stability constants with the first YREE hydrolysis constant, log , yield regression coefficients of >0.97. On the other hand, plots of the constants vs. the radius of the inner hydration sphere display an increasing deviation from linearity for the lightest REEs (La > Pr > Nd). This may signify steric constraints in DFOB folding around bulkier cations, a larger mismatch in coordination number, or a substantial degree of covalence in the YREE-hydroxamate bond.Complexes of the YREEs with DFOB are many orders of magnitude more stable than those with carbonate, the dominant inorganic YREE ligand in seawater. Speciation modeling with MINEQL indicates that, for an average seawater composition, the hexadentate complex could constitute as much as 28% of dissolved Lu at free DFOB concentrations as low as 10⁻¹³ M. Such conditions might occur when DFOB or other siderophores are present in excess over metals for which they have high affinity, like Fe(III) and Co(III), for example during plankton blooms. Even if it turns out that trihydroxamate siderophores are not the dominant organic YREE ligand in seawater, our results establish a benchmark for producing effects on YREE solution speciation comparable to that of DFOB: the free concentration of any weaker organic ligand L must exceed that of DFOB by a factor β₃/_Lβ₁, assuming its first order complex is formed in greatest abundance.     

Empirical prediction of the Pitzer’s interaction parameters for cationic Al species with both SiO2(aq) and CO2(aq): Implications for the geochemical modelling of very saline solutions

The Pitzer’s interaction parameters, λ_N–M, involving the Mth cationic Al species Al³⁺ or AlOH²⁺ or AlO⁺ and the Nth neutral species SiO_2(aq) (at temperatures of 25–300 °C) or CO_2(aq) (at temperatures of 25–150 °C), have been evaluated through empirical linear relationships between λ_N–M and the surface electrostatic field of the ionic species of interest. These relationships have been obtained starting from the known λ_N–M for both SiO_2(aq) and CO_2(aq) with the main dissolved cations. The Pitzer’s interaction parameter thus estimated for the pair CO_2(aq)–Al³⁺ at 25 °C, 0.327, is 20–40% higher than the corresponding values obtained from CO₂ solubilities in concentrated solutions of AlCl₃, 0.272 ± 0.010 (2σ), and Al₂(SO₄)₃, 0.232 ± 0.002 (2σ), partly corroborating the empirical approach adopted in this study. To test the Pitzer’s interaction coefficients for cationic Al species with aqueous SiO₂, the log K values of the kaolinite dissolution reaction have been computed starting from available experimental data at 23–25 °C and ionic strengths of 0.0001–0.12 mol/kg adopting, alternatively, the Pitzer’s equations and the Debye–Hückel equation. A satisfactory agreement has been found between the log K values obtained through these two approaches, with maximum deviations of 0.11–0.12 log units. This good convergence of results is encouraging as it represents a necessary condition to prove the reliability of the Pitzer’s interaction coefficients estimated in this work. These results are a first step to take into account specific interactions among solutes in concentrated electrolyte solutions, such as those hosted in sedimentary basins or geothermal waters, for instance through the Pitzer’s equations. However, experimental or field data at higher ionic strengths are absolutely necessary to validate the reliability of the Pitzer’s interaction coefficients determined in this study.     

In situ adsorption of mercury,methylmercury and other elements by iron oxyhydroxides and organic matter in lake sediments

Samples of authigenic material, sediment overlying water and oxic surface sediment (0–0.5-cm depth) from a perennially oxygenated lacustrine basin were analysed to investigate which solid phases are important for binding a suite of trace elements (Ag, As, Ca, Cd, Cu, Hg, In, methylmercury (MeHg), Mg, Mo, Pb, Sb and Zn). The authigenic material, which was collected with inert Teflon sheets deployed for several years across the sediment–water interface, contained mainly poorly crystallized Fe oxyhydroxides and natural organic matter, presumably humic substances derived from the watershed. Manganese oxyhydroxides were not present in the collected authigenic material due to the slightly acidic condition (pH = 5.6) of the lake that prevents the formation and recycling of these compounds. Conditional equilibrium constants for the adsorption of cationic (K_Fe–M) and anionic (K_Fe–A) trace elements onto the authigenic Fe oxyhydroxides were estimated from their concentrations in the authigenic material and in bottom water samples. These field-derived values of K_Fe–M and K_Fe–A were compared with those predicted by the surface complexation model, using laboratory-derived intrinsic adsorption constants and the water composition at the study site. Equilibrium constants (K_POM–M) were also calculated for the adsorption of the cationic trace elements onto the humic substances contained in the diagenetic material. The field-derived values of K_POM–M were compared to those predicted by the speciation code WHAM 6 for the complexation of the trace elements by dissolved humic substances in the lake. Combining the results of the present study with those on the distributions of trace elements in the porewater and solid-phase sediments reported in previous studies at the same site, it was determined whether the trace elements bind preferentially to Fe oxyhydroxides or natural organic matter in oxic sediments. The main inferences are that the anionic trace elements As, Mo and Sb, as well as the cationic metal Pb are preferentially bound to the authigenic Fe oxyhydroxides whereas the other trace elements, and especially Hg and MeHg, are preferentially bound to the humic substances.     

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:

     

Modelling the sorption of Mn(II), Co(II), Ni(II), Zn(II), Cd(II), Eu(III), Am(III), Sn(IV), Th(IV), Np(V) and U(VI) on montmorillonite: Linear free energy relationships and estimates of surface binding constants for some selected heavy metals and actinides

In solution thermodynamics, and more recently in surface chemistry, it is well established that relationships can be found between the free energies of formation of aqueous or surface metal complexes and thermodynamic properties of the metal ions or ligands. Such systematic dependencies are commonly termed linear free energy relationships (LFER). A 2 site protolysis non-electrostatic surface complexation and cation exchange (2SPNE SC/CE) model has been used to model “in house” and literature sorption edge data for eleven elements: Mn(II), Co(II), Ni(II), Zn(II), Cd(II), Eu(III), Am(III), Sn(IV), Th(IV), Np(V) and U(VI) to provide surface complexation constants for the strong sites on montmorillonite. Modelling a further 4 sets of sorption isotherms for Ni(II), Zn(II), Eu(III) and U(VI) provided complexation constants for the weak sites. The protolysis constants and site capacities derived for the 2SPNE SC/CE model in previous work were fixed in all of the calculations. Cation exchange was modelled simultaneously to provide selectivity coefficients. Good correlations between the logarithms of strong ^SK_x−1 and weak ^W1K_x−1 site binding constants on montmorillonite and the logarithm of the aqueous hydrolysis constants ^OHK_x were found which could be described by the following equations: Strong (≡S^SOH) sites:

SlogKX−1=8.1±0.3+(0.90±0.02)log^OHKX     

Complexation of cadmium to sulfur and oxygen functional groups in an organic soil

Cadmium (Cd) is a toxic trace element and due to human activities soils and waters are contaminated by Cd both on a local and global scale. It is widely accepted that chemical interactions with functional groups of natural organic matter (NOM) is vital for the bioavailability and mobility of trace elements. In this study the binding strength of cadmium (Cd) to soil organic matter (SOM) was determined in an organic (49% organic C) soil as a function of reaction time, pH and Cd concentration. In experiments conducted at native Cd concentrations in soil (0.23 μg g⁻¹ dry soil), halides (Cl, Br) were used as competing ligands to functional groups in SOM. The concentration of Cd in the aqueous phase was determined by isotope-dilution (ID) inductively-coupled-plasma-mass-spectrometry (ICP-MS), and the activity of Cd²⁺ was calculated from the well-established Cd-halide constants. At higher Cd loading (500-54,000 μg g⁻¹), the Cd²⁺ activity was directly determined by an ion-selective electrode (ISE). On the basis of results from extended X-ray absorption fine structure (EXAFS) spectroscopy, a model with one thiolate group (RS⁻) was used to describe the complexation (Cd²⁺ + RS⁻ ? CdSR⁺; log K_CdSR) at native Cd concentrations. The concentration of thiols (RSH; 0.047 mol kg⁻¹ C) was independently determined by X-ray absorption near-edge structure (XANES) spectroscopy. Log K_CdSR values of 11.2-11.6 (pK_a for RSH = 9.96), determined in the pH range 3.1-4.6, compare favorably with stability constants for the association between Cd and well-defined thiolates like glutathione. In the concentration range 500-54,000 μg Cd g⁻¹, a model consisting of one thiolate and one carboxylate (RCOO⁻) gave the best fit to data, indicating an increasing role for RCOOH groups as RSH groups become saturated. The determined log K_CdOOCR of 3.2 (Cd²⁺ +  RCOO⁻ ? CdOOCR⁺; log K_CdOOCR; pK_a for RCOOH = 4.5) is in accordance with stability constants determined for the association between Cd and well-defined carboxylates. Given a concentration of reduced sulfur groups of 0.2% or higher in NOM, we conclude that the complexation to organic RSH groups may control the speciation of Cd in soils, and most likely also in surface waters, with a total concentration less than 5 mg Cd g⁻¹ organic C.     

Temperature dependence of sulfide and sulfate solubility in olivine-saturated basaltic magmas

The sulfur concentration at pyrrhotite- and anhydrite-saturation in primitive hydrous basaltic melt of the 2001-2002 eruption of Mt. Etna was determined at 200 MPa, T = 1050-1250 °C and at log fO₂ from FMQ to FMQ+2.2 (FMQ is Fayalite-Magnetite-Quartz oxygen buffer). At 1050 °C Au sample containers were used. A double-capsule technique, using a single crystal olivine sample container closed with an olivine piston, embedded in a sealed Au₈₀Pd₂₀ capsule, was developed to perform experiments in S-bearing hydrous basaltic systems at T > 1050 °C. Pyrrhotite is found to be a stable phase coexisting with melt at FMQ-FMQ+0.3, whereas anhydrite is stable at FMQ+1.4-FMQ+2.2. The S concentration in the melt increases almost linearly from 0.12 ± 0.01 to 0.39 ± 0.02 wt.% S at FeS-saturation and from 0.74 ± 0.01 to 1.08 ± 0.04 wt.% S at anhydrite-saturation with T ranging from 1050-1250 °C. The relationships between S concentration at pyrrhotite and/or anhydrite saturation, MgO content of the olivine-saturated melt, T, and log fO₂ observed in this study and from previous data are used to develop an empirical model for estimating the magmatic T and fO₂ from the S and MgO concentrations of H₂O-bearing olivine-saturated basaltic melts. The model can also be used to determine maximum S concentrations, if fO₂ and MgO content of the melt are known. The application of the model to compositions of melt inclusions in olivines from Mt. Etna indicates that the most primitive magmas trapped in inclusions might have been stored at log fO₂ slightly higher than FMQ+1 and at T = 1100-1150 °C, whereas more evolved melts could have been trapped at T ? 1100 °C. These values are in a good agreement with the estimates obtained by other independent methods reported in the literature.     

Vapor pressures and evaporation coefficients for melts of ferromagnesian chondrule-like compositions

To determine evaporation coefficients for the major gaseous species that evaporate from silicate melts, the Hertz-Knudsen equation was used to model the compositions of residues of chondrule analogs produced by evaporation in vacuum by Hashimoto [Hashimoto A. (1983) Evaporation metamorphism in the early solar nebula-evaporation experiments on the melt FeO-MgO-SiO₂-CaO-Al₂O₃ and chemical fractionations of primitive materials. Geochem. J. 17, 111-145] and Wang et al. [Wang J., Davis A. M., Clayton R. N., Mayeda T. K., Hashimoto A. (2001) Chemical and isotopic fractionation during the evaporation of the FeO-MgO-SiO₂-CaO-Al₂O₃-TiO₂ rare earth element melt system. Geochim. Cosmochim. Acta 65, 479-494], in vacuum and in H₂ by Yu et al. [Yu Y., Hewins R. H., Alexander C. M. O’D., Wang J. (2003) Experimental study of evaporation and isotopic mass fractionation of potassium in silicate melts. Geochim. Cosmochim. Acta 67, 773-786], and in H₂ by Cohen et al. [Cohen B. A., Hewins R. H., Alexander C. M. O’D. (2004) The formation of chondrules by open-system melting of nebular condensates. Geochim. Cosmochim. Acta 68, 1661-1675]. Vapor pressures were calculated using the thermodynamic model of Ghiorso and Sack [Ghiorso M. S., Sack R. O. (1995) Chemical mass transfer in magmatic processes IV. A revised and internally consistent thermodynamic model for the interpolation and extrapolation of liquid-solid equilibria in magmatic systems at elevated temperatures and pressures. Contrib. Mineral. Petrol. 119, 197-212], except for the late, FeO-free stages of the Wang et al. (2001) and Cohen et al. (2004) experiments, where the CMAS activity model of Berman [Berman R. G. (1983) A thermodynamic model for multicomponent melts, with application to the system CaO-MgO-Al₂O₃-SiO₂. Ph.D. thesis, University of British Columbia] was used. From these vapor pressures, evaporation coefficients (α) were obtained that give the best fits to the time variation of the residue compositions. Evaporation coefficients derived for Fe_(g), Mg_(g), and SiO_(g) from the Hashimoto (1983) experiments are similar to those found by Alexander [Alexander C. M. O’D. (2004) Erratum. Meteoritics Planet. Sci. 39, 163] in his EQR treatment of the same data and also adequately describe the FeO-bearing stages of the Wang et al. (2001) experiments. From the Yu et al. (2003) experiments at 1723 K, α_Na = 0.26 ± 0.05, and α_K = 0.13 ± 0.02 in vacuum, and α_Na = 0.042 ± 0.020, andα_K = 0.017 ± 0.002 in 9 × 10⁻⁵ bar H₂. In the FeO-free stages of the Wang et al. (2001) experiments, α_Mg and α_SiO are significantly different from their respective values in the FeO-bearing portions of the same experiments and from the vacuum values obtained at the same temperature by Richter [Richter F. M., Davis A. M., Ebel D. S., Hashimoto A. (2002) Elemental and isotopic fractionation of Type B calcium-, aluminum-rich inclusions: experiments, theoretical considerations, and constraints on their thermal evolution. Geochim. Cosmochim. Acta 66, 521-540] for CMAS compositions much lower in MgO. When corrected for temperature, the values of α_Mg and α_SiO that best describe the FeO-free stages of the Wang et al. (2001) experiments also adequately describe the FeO-free stage of the Cohen et al. (2004) H₂ experiments, but α_Fe that best describes the FeO-bearing stage of the latter experiment differs significantly from the temperature-corrected value derived from the Hashimoto (1983) vacuum data.     

Groundwater chemistry of Al under Dutch sandy soils: Effects of land use and depth

Aluminium has received great attention in the second half of the 20th century, mainly in the context of the acid rain problem mostly in forest soils. In this research the effect of land use and depth of the groundwater on Al, pH and DOC concentration in groundwater under Dutch sandy soils has been studied. Both pH and DOC concentration play a major role in the speciation of Al in solution. Furthermore, the equilibrium with mineral phases like gibbsite, amorphous Al(OH)₃ and imogolite, has been considered. Agricultural and natural land use were expected to have different effects on the pH and DOC concentration, which in turn could influence the total Al concentration and the speciation of Al in groundwater at different depths (phreatic, shallow and deep). An extensive dataset (n = 2181) from the national and some provincial monitoring networks on soil and groundwater quality was used. Land use type and groundwater depth did influence the pH, and Al and DOC concentrations in groundwater samples. The Al concentration ranged from <0.4 μmol L⁻¹ at pH > 7 to 1941 μmol L⁻¹ at pH < 4; highest Al concentrations were found for natural-phreatic groundwater. The DOC concentration decreased and the median pH increased with depth of the groundwater. Natural-phreatic groundwater showed lower pH than the agricultural-phreatic groundwater. Highest DOC concentrations were found for the agricultural-phreatic groundwater, induced by the application of organic fertilizers. Besides inorganic complexation, the NICA-Donnan model was used to calculate Al³⁺ concentrations for complexation with DOC. Below pH 4.5 groundwater samples were mainly in disequilibrium with a mineral phase. This disequilibrium is considered to be the result of kinetic constraints or equilibrium with organic matter. Log K values were derived by linear regression and were close to theoretical values for Al(OH)₃ minerals (e.g. gibbsite or amorphous Al(OH)₃), except for natural-phreatic groundwater for which lower log K values were found. Complexation of Al with DOC is shown to be an important factor for the Al concentrations, especially at high DOC concentrations as was found for agricultural-phreatic groundwater.     

The effect of species diversity on metal adsorption onto bacteria

In this study, we measure proton, Pb, and Cd adsorption onto the bacteria Deinococcus radiodurans, Thermus thermophilus, Acidiphlium angustum, Flavobacterium aquatile, and Flavobacterium hibernum, and we calculate the thermodynamic stability constants for the important surface complexes. These bacterial species represent a wide genetic diversity of bacteria, and they occupy a wide range of habitats. All of the species, except for A. angustum, exhibit similar proton and metal uptake. The only species tested that exhibits significantly different protonation behavior is A. angustum, an acidophile that grows at significantly lower pH than the other species of this study. We demonstrate that a single, metal-specific, surface complexation model can be used to reasonably account for the acid/base and metal adsorption behaviors of each species. We use a four discrete site non-electrostatic model to describe the protonation of the bacterial functional groups, with averaged pK_a values of 3.1 ± 0.3, 4.8 ± 0.2, 6.7 ± 0.1, and 9.2 ± 0.3, and site concentrations of (1.0 ± 0.17) × 10⁻⁴, (9.0 ± 3.0) × 10⁻⁵, (4.6 ± 1.8) × 10⁻⁵, and (6.1 ± 2.3) × 10⁻⁵ mol of sites per gram wet mass of bacteria, respectively. Adsorption of Cd and Pb onto the bacteria can be accounted for by the formation of complexes with each of the bacterial surface sites. The average log stability constants for Cd complexes with Sites 1-4 are 2.4 ± 0.4, 3.2 ± 0.1, 4.4 ± 0.1, and 5.3 ± 0.1, respectively. The average log stability constants for Pb complexes with Sites 1-4 are 3.3 ± 0.2, 4.5 ± 0.3, 6.5 ± 0.1, and 7.9 ± 0.5, respectively. This study demonstrates that a wide range of bacteria exhibit similar proton and metal adsorption behaviors, and that a single set of averaged acidity constants, site concentrations, and stability constants for metal-bacterial surface complexes yields a reasonable model for the adsorption behavior of many of these species. The differences in adsorption behavior that we observed for A. angustum demonstrate that genetic differences do exist between the cell wall functional group chemistries of some bacterial species, and that significant exceptions to the typical bacterial adsorption behavior do exist.     

Oxygen fugacity dependence of Os solubility in haplobasaltic melt

Os equilibrium solubilities were determined at 1350 °C over a wide range of oxygen fugacities (−12 < log fO₂ < −7) applying the mechanically assisted equilibration technique (MAE) at 10⁵ Pa (= 1 bar). Os concentrations in the glass samples were analysed using ID-NTIMS. Additional LA-ICP-MS and SEM analyses were performed to detect, visualize and analyse the nature and chemistry of “nanonuggets.” Os solubilities determined range at a constant temperature of 1350 °C from 0.63 ± 0.04 to 37.4 ± 1.16 ppb depending on oxygen fugacity. At the highest oxygen fugacities, Os³⁺ can be confirmed as the main oxidation state of Os. At low oxygen fugacities (below log fO₂ = −8), samples are contaminated by nanonuggets which, despite the MAE technique, were still not removed entirely from the melt. However, the present results indicate that applying MAE technology does reduce the amount of nanonuggets present significantly, resulting in the lowest Os solubility results reported to date under these experimental conditions, and extending the experimentally accessible range of fO₂ for these studies to lower values. Calculated metal/silicate melt partition coefficients are therefore higher compared to previous studies, making Os more siderophile. Neglecting the as yet unknown temperature dependence of the Os metal/silicate melt partition coefficient, extrapolation of the obtained Os solubilities to conditions for core-mantle equilibrium, results in a , while metallic alloy/silicate melt partition coefficients range from 1.4 × 10⁶ to 8.6 × 10⁷, in agreement with earlier findings. Therefore remains too high by 2-4 orders of magnitude to explain the Os abundance in the Earth’s mantle as result of core-mantle equilibrium during core formation.