Europium retention onto clay minerals from 25 to 150 °C: Experimental measurements, spectroscopic features and sorption modelling

The sorption of Eu(III) onto kaolinite and montmorillonite was investigated up to 150 °C. The clays were purified samples, saturated with Na in the case of montmorillonite. Batch experiments were conducted at 25, 40, 80 and 150 °C in 0.5 M NaClO₄ solutions to measure the distribution coefficients (Kd) of Eu as a trace element (<10⁻⁶ mol/L) between the solution and kaolinite. For the Na-montmorillonite, we used Kd results from a previous study [Tertre, E., Berger, G., Castet, S., Loubet, M., Giffaut, E., 2005. Experimental study of adsorption of Ni²⁺, Cs⁺ and Ln³⁺ onto Na-montmorillonite up to 150 °C. Geochim. Cosmochim. Acta69, 4937-4948] obtained under exactly the same conditions. The number and nature of the Eu species sorbed onto both clay minerals were investigated by time resolved laser fluorescence spectroscopy (TRLFS) in specific experiments in the same temperature range. We identified a unique inner-sphere complex linked to the aluminol sites in both clays, assumed to be AlOEu²⁺ at the edge of the particles, and a second exchangeable outer-sphere complex for montmorillonite, probably in an interlayer position. The Kd values were used to adjust the parameters of a surface complexation model (DLM: diffuse layer model) from 25 to 150 °C. The number of Eu complexes and the stoichiometry of reactions were constrained by TRLFS. The acidity constants of the amphoteric aluminol sites were taken from another study [Tertre, E., Castet, S., Berger, G., Loubet, M., Giffaut, E. Acid/base surface chemistry of kaolinite and Na-montmorillonite at 25 and 60 °C: experimental study and modelling. Geochim. Cosmochim. Acta, in press], which integrates the influence of the negative structural charge of clays on the acid/base properties of edge sites as a function of temperature and ionic strength. The results of the modelling show that the observed shift of the sorption edge towards low pH with increasing temperature results solely from the contribution of the AlOEu²⁺ edge complexes. Finally, we successfully tested the performance of our model by confronting the predictions with experimental Kd data. We used our own data obtained at lower ionic strength (previous study) or higher suspension density and higher starting concentration (TRLFS runs, this study), as well as published data from other experimental studies [Bradbury, M.H., Baeyens, B., 2002. Sorption of Eu on Na and Ca-montmorillonite: experimental investigations and modeling with cation exchange and surface complexation. Geochim. Cosmochim. Acta66, 2325-2334; Kowal-Fouchard, A., 2002. Etude des mécanismes de rétention des ions U(IV) et Eu(III) sur les argiles: influence des silicates. Ph.D. Thesis, Université Paris Sud, France, 330p].     

Experimental study of neptunyl adsorption onto Bacillus subtilis

The subsurface mobility of Np is difficult to predict in part due to uncertainties associated with its sorption behavior in geologic systems. In this study, we measured Np adsorption onto a common gram-positive soil bacterium, Bacillus subtilis. We performed batch adsorption experiments with Np(V) solutions as a function of pH, from 2.5 to 8, as a function of total Np concentration from 1.29 × 10⁻⁵ M to 2.57 × 10⁻⁴ M, and as a function of ionic strength from 0.001 to 0.5 M NaClO₄. Under most pH conditions, Np adsorption is reversible and exhibits an inverse relationship with ionic strength, with adsorption increasing with increasing pH. At low pH in the 0.1 M ionic strength systems, we observed irreversible adsorption, which is consistent with reduction of Np(V) to Np(IV). We model the adsorption reaction using a nonelectrostatic surface complexation approach to yield ionic strength dependent NpO₂⁺-bacterial surface stability constants. The data require two bacterial surface complexation reactions to account for the observed adsorption behavior: R-L₁⁻ + NpO₂⁺ ↔ R-L₁-NpO₂^° and R-L₂⁻ + NpO₂⁺ ↔ R-L₂-NpO₂^°, where R represents the bacterium to which each functional group is attached, and L₁ and L₂ represent the first and second of four discrete site types on the bacterial surface. Stability constants (log K values) for the L₁ and L₂ reactions in the 0.001 M system are 2.3 ± 0.3 and 2.3 ± 0.2, and in the 0.1 M system the values are 1.7 ± 0.2 and 1.6 ± 0.2, respectively. The calculated neptunyl-bacterial surface stability constants are not consistent with values predicted using the linear free energy correlation approach from Fein et al. (2001), suggesting that possible unfavorable steric interactions and the low charge of NpO₂⁺ affects Np-bacterial adsorption.     

Kinetics and mechanisms of the leaching of low Fe sphalerite

The surface speciation and leaching kinetics of 38- to 75-μm sphalerite (0.45 wt.% Fe) particles reacted in O₂ purged perchloric acid (at pH 1.0) at 25, 40, 60, and 85 °C over a leach period of 144 h were investigated. In all cases, an initial rapid leach rate is observed followed by a slower leach rate. These two leach regimes can each be adequately modeled using straight-line interpolation, and thus two activation energies (E_a) have been derived. E_a for the fast and slow Zn dissolution rates were 33 ± 4 kJ mol⁻¹ and 34 ± 4 kJ mol⁻¹ respectively, suggesting the same rate-determining step.     

Effects of pH and temperature on dimerization rate of glycine: Evaluation of favorable environmental conditions for chemical evolution of life

To evaluate favorable environmental conditions for the chemical evolution of life, we studied the effects of pH and temperature on the dimerization rate of glycine (Gly: NH₂-CH₂-COOH), one of the simplest amino acids. Gly dimerizes to form glycylglycine (GlyGly), and GlyGly further reacts to form diketopiperazine (DKP). Gly solutions with pH ranging from 3.1 to 10.9 were heated for 1-14 days at 140 °C, and changes in concentrations of Gly, GlyGly, and DKP were evaluated. At pH 9.8, the experiments were conducted at 120, 140, 160, and 180 °C. The dimerization rate of Gly was nearly constant at pH 3-7 and increased with increasing pH from 7 to 9.8 and then decreased with further increases in pH. To elucidate the reason for this pH dependency, we evaluated the role of the three dissociation states of Gly (cationic state: Gly⁺, zwitterionic state: Gly^±, and anionic state: Gly⁻). For pH >6, the dominant forms are Gly^± and Gly⁻, and the molar fraction of Gly^± decreases and that of Gly⁻ increases with increasing pH. The dimerization rate was determined for each dissociation state. The reaction between Gly^± and Gly⁻ was found to be the fastest; the rate constant of the reaction between Gly^± and Gly⁻ was 10 times the size of that between Gly⁻ and Gly⁻ and 98 times that between Gly^± and Gly^±. The dimerization rate became greatest at pH 9.8 because the molar fractions of Gly^± and Gly⁻ are approximately equal at this pH. The dimerization rate increased with temperature, and an activation energy of 88 kJ mol⁻¹ was obtained. Based on these results and previous reports on the stability of amino acids under hydrothermal conditions, we determined that Gly dimerizes most efficiently under alkaline pH (∼9.8) at about 150 °C.     

Rare earth element sorption by basaltic rock: Experimental data and modeling results using the “Generalised Composite approach”

Sorption of the 14 rare earth elements (REE) by basaltic rock is investigated as a function of pH, ionic strength and aqueous REE concentrations. The rock sample, originating from a terrestrial basalt flow (Rio Grande do Sul State, Brazil), is composed of plagioclase, pyroxene and cryptocrystalline phases. Small amounts of clay minerals are present, due to rock weathering. Batch sorption experiments are carried out under controlled temperature conditions of 20 °C with the <125 μm fraction of the ground rock in solutions of 0.025 M and 0.5 M NaCl and at pH ranging from 2.7 to 8. All 14 REEs are investigated simultaneously with initial concentrations varying from 10⁻⁷ to 10⁻⁴ mol/L. Some experiments are repeated with only europium present to evaluate possible competitive effects between REE. Experimental results show the preferential retention of the heavy REEs at high ionic strength and circumneutral pH conditions. Moreover, results show that REE sorption increases strongly with decreasing ionic strength, indicating two types of sorption sites: exchange and specific sites. Sorption data are described by a Generalised Composite (GC) non-electrostatic model: two kinds of surface reactions are treated, i.e. cation exchange at >XNa sites, and surface complexation at >SOH sites. Total site density (>XNa + >SOH) is determined by measuring the cation exchange capacity (CEC = 52 μmol/m²). Specific concentrations of exchange sites and complexation sites are determined by fitting the Langmuir equation to sorption isotherms of REE and phosphate ions. Site densities of 22 ± 5 and 30 ± 5 μmol/m² are obtained for [>XNa] and [>SOH], respectively. The entire set of REE experimental data is modeled using a single exchange constant (log K_ex = 9.7) and a surface complexation constant that progressively increases from log K = −1.15 for La(III) to −0.4 for Lu(III).The model proves to be fairly robust in describing other aluminosilicate systems. Maintaining the same set of sorption constants and only adjusting the site densities, we obtain good agreement with the literature data on REE/kaolinite and REE/smectite sorption. The Generalised Composite non-electrostatic model appears as an easy and efficient tool for describing sorption by complex aluminosilicate mineral assemblages.     

Synthesis, characterization, and thermochemistry of K-Na-H3O jarosites

The thermochemistry of well-characterized synthetic K-H₃O, Na-H₃O and K-Na-H₃O jarosites was investigated. These phases are solid solutions that obey Vegard’s law. Electron probe microanalyses indicated lower alkali and iron contents than predicted from the theoretical end-member compositions, in agreement with thermal analyses, suggesting the presence of hydronium and “additional” water. The standard enthalpies of formation (ΔH°_f) of K-H₃O, Na-H₃O and K-Na-H₃O jarosites were determined by high-temperature oxide melt solution calorimetry. These enthalpies vary linearly with the K/H₃O, Na/H₃O and K/Na ratio, respectively. The enthalpy of formation of pure hydronium jarosite was also determined experimentally (ΔH°_f = −3741.6 ± 8.3 kJ.mol⁻¹), and it was used to evaluate ΔH°_f for the end-members KFe₃(SO₄)₂(OH)₆ (ΔH°_f = −3829.6 ± 8.3 kJ.mol⁻¹) and NaFe₃(SO₄)₂(OH)₆ (ΔH°_f = −3783.4 ± 8.3 kJ.mol⁻¹). Finally, enthalpies of dehydration (loss of the “additional” water) of some jarosites were determined and found to be near the enthalpy of vaporization of water, suggesting that the “additional” water is weakly bonded in the structure.     

Thermochemistry of yavapaiite KFe(SO4)2: Formation and decomposition

Yavapaiite, KFe(SO₄)₂, is a rare mineral in nature, but its structure is considered as a reference for many synthetic compounds in the alum supergroup. Several authors mention the formation of yavapaiite by heating potassium jarosite above ca. 400°C. To understand the thermal decomposition of jarosite, thermodynamic data for phases in the K-Fe-S-O-(H) system, including yavapaiite, are needed. A synthetic sample of yavapaiite was characterized in this work by X-ray diffraction (XRD), scanning electron microscopy (SEM), Fourier transform infrared spectroscopy (FTIR), and thermal analysis. Based on X-ray diffraction pattern refinement, the unit cell dimensions for this sample were found to be a = 8.152 ± 0.001 Å, b = 5.151 ± 0.001 Å, c = 7.875 ± 0.001 Å, and β = 94.80°. Thermal decomposition indicates that the final breakdown of the yavapaiite structure takes place at 700°C (first major endothermic peak), but the decomposition starts earlier, around 500°C. The enthalpy of formation from the elements of yavapaiite, KFe(SO₄)₂, ΔH°_f = −2042.8 ± 6.2 kJ/mol, was determined by high-temperature oxide melt solution calorimetry. Using literature data for hematite, corundum, and Fe/Al sulfates, the standard entropy and Gibbs free energy of formation of yavapaiite at 25°C (298 K) were calculated as S°(yavapaiite) = 224.7 ± 2.0 J.mol⁻¹.K⁻¹ and ΔG°_f = −1818.8 ± 6.4 kJ/mol. The equilibrium decomposition curve for the reaction jarosite = yavapaiite + Fe₂O₃ + H₂O has been calculated, at pH₂O = 1 atm, the phase boundary lies at 219 ± 2°C.     

Thermodynamics of iron oxides: Part III. Enthalpies of formation and stability of ferrihydrite (∼Fe(OH)3), schwertmannite (∼FeO(OH)3/4(SO4)1/8), and ε-Fe2O3

Enthalpies of formation of ferrihydrite and schwertmannite were measured by acid solution calorimetry in 5 N HCl at 298 K. The published thermodynamic data for these two phases and ε-Fe₂O₃ were evaluated, and the best thermodynamic data for the studied compounds were selected.Ferrihydrite is metastable in enthalpy with respect to α-Fe₂O₃ and liquid water by 11.5 to 14.7 kJ•mol⁻¹ at 298.15 K. The less positive enthalpy corresponds to 6-line ferrihydrite, and the higher one, indicating lesser stability, to 2-line ferrihydrite. In other words, ferrihydrite samples become more stable with increasing crystallinity. The best thermodynamic data set for ferrihydrite of composition Fe(OH)₃ was selected by using the measured enthalpies and (1) requiring ferrihydrite to be metastable with respect to fine-grained lepidocrocite; (2) requiring ferrihydrite to have entropy higher than the entropy of hypothetical, well-crystalline Fe(OH)₃; and (3) considering published estimates of solubility products of ferrihydrite. The ΔG°_f for 2-line ferrihydrite is best described by a range of −708.5±2.0 to −705.2±2.0 kJ•mol⁻¹, and ΔG°_f for 6-line ferrihydrite by −711.0±2.0 to −708.5±2.0 kJ•mol⁻¹.A published enthalpy measurement by acid calorimetry of ε-Fe₂O₃ was re-evaluated, arriving at ΔH°_f (ε-Fe₂O₃) = −798.0±6.6 kJ•mol⁻¹. The standard entropy (S°) of ε-Fe₂O₃ was considered to be equal to S° (γ-Fe₂O₃) (93.0±0.2 J•K⁻¹•mol⁻¹), giving ΔG°_f (ε-Fe₂O₃) = −717.8±6.6 kJ•mol⁻¹. ε-Fe₂O₃ thus appears to have no stability field, and it is metastable with respect to most phases in the Fe₂O₃-H₂O system which is probably the reason why this phase is rare in nature.Enthalpies of formation of two schwertmannite samples are: ΔH°_f (FeO(OH)_0.686(SO₄)_0.157•0.972H₂O) = −884.0±1.3 kJ•mol⁻¹, ΔH°_f (FeO(OH)_0.664(SO₄)_0.168•1.226H₂O) = −960.7±1.2 kJ•mol⁻¹. When combined with an entropy estimate, these data give Gibbs free energies of formation of −761.3 ± 1.3 and −823.3 ± 1.2 kJ•mol⁻¹ for the two samples, respectively. These ΔG_f° values imply that schwertmannite is thermodynamically favored over ferrihydrite over a wide range of pH (2-8) when the system contains even small concentration of sulfate. The stability relations of the two investigated samples can be replicated by schwertmannite of the “ideal” composition FeO(OH)_3/4(SO₄)_1/8 with ΔG°_f = −518.0±2.0 kJ•mol⁻¹.     

Activation energy for diffusion of helium in water-saturated, compacted Na-montmorillonite

It is important to clarify the migration behavior of hydrogen gas dissolved in water-saturated, compacted bentonite, which is a promising material for geologic disposal of high-level waste and TRU waste disposal. The diffusion coefficients of helium, which can be detected under extremely low background conditions, in water-saturated, compacted Na-montmorillonite were determined as a function of temperature by a transient diffusion method. The activation energies for diffusion of helium were then obtained. The activation energies were from 6.9 ± 4.8 to 19 ± 2.8 kJ mol^− 1 and were regarded to be independent of dry density. The activation energies of helium in water-saturated Na-montmorillonite were roughly equal to those in bulk water, 14.9 kJ mol^− 1, and in ice, from 11 to 13 kJ mol^− 1. It is possible that helium diffuses not only in pore water but also in interlayer water.     

Fe sorption onto nontronite (NAu-2)

The sorption of ferrous iron to a clay mineral, nontronite (NAu-2, a ferruginous smectite), was investigated under strictly anoxic conditions as a function of pH (3-10), Fe²⁺ concentration (0.01-50 mM), equilibration time (1-35 days), and ionic strength (0.01-0.5 M NaClO₄). The surface properties of NAu-2 were independently characterized to determine its fixed charge and amphoteric site density in order to interpret the Fe²⁺ sorption data. Fe²⁺ sorption to NAu-2 was strongly dependent on pH and ionic strength, reflecting the coupled effects of Fe²⁺ sorption through ion exchange and surface complexation reactions. Fe²⁺ sorption to NAu-2 increased with increasing pH from pH 2.5 to 4.5, remained constant from pH 4.5 to 7.0, increased again with further increase of pH from pH 7.0 to 8.5, and reached a maximum above pH 8.5. The Fe²⁺ sorption below pH 7.0 increased with decreasing ionic strength. The differences of Fe²⁺ sorption at different ionic strengths, however, diminished with increasing equilibration time. The Fe²⁺ sorption from pH 4.5 to 7.0 increased with increasing equilibration time up to 35 days and showed stronger kinetic behavior in higher ionic strength solutions. The kinetic uptake of Fe²⁺ onto NAu-2 is consistent with a surface precipitation mechanism although our measurements were not able to identify secondary precipitates. An equilibrium model that integrates ion exchange, surface complexation and aqueous speciation reactions reasonably well describes the Fe²⁺ sorption data as a function of pH, ionic strength, and Fe²⁺ concentration measured at 24 h of equilibration. Model calculations show that the species Fe(OH)⁺ was required to describe Fe²⁺ sorption above pH 8.0 satisfactorily. Overall, this study demonstrated that Fe²⁺ sorption to NAu-2 is affected by complex equilibrium and kinetic processes, likely caused by surface precipitation reactions.     

Impact of highly basic solutions on sorption of Cs to subsurface sediments from the Hanford site, USA

The effect of caustic NaNO₃ solutions on the sorption of ¹³⁷Cs to a Hanford site micaceous subsurface sediment was investigated as a function of base exposure time (up to 168 d), temperature (10°C or 50°C), and NaOH concentration (0.1 mol/L to 3 mol/L). At 10°C and 0.1 M NaOH, the slow evolution of [Al]_aq was in stark contrast to the rapid increase and subsequent loss of [Al]_aq observed at 50°C (regardless of base concentration). Exposure to 0.1 M NaOH at 10°C for up to 168 d exhibited little if any measurable effect on sediment mineralogy, Cs⁺ sorption, or Cs⁺ selectivity; sorption was well described with a two-site ion exchange model modified to include enthalpy effects. At 50°C, dissolution of phyllosilicate minerals increased with [OH]. A zeolite (tetranatrolite; Na₂Al₂Si₃O₁₀·2H₂O) precipitated in 0.1 M NaOH after about 7 days, while an unnamed mineral phase (Na₁₄Al₁₂Si₁₃O₅₁·6H₂O) precipitated after 4 and 2 days of exposure to 1 M and 3 M NaOH solutions, respectively. Short-term (16 h) Cs⁺ sorption isotherms (10⁻⁹-10⁻² mol/L) were measured on sediment after exposure to 0.1 M NaOH for 56, 112, and 168 days at 50°C. There was a trend toward slightly lower conditional equilibrium exchange constants (Δlog _Na^CsK_c ∼ 0.25) over the entire range of surface coverage, and a slight loss of high affinity sites (15%) after 168 days of pretreatment with 0.1 M base solution. Cs⁺ sorption to sediment over longer times was also measured at 50°C in the presence of NaOH (0.1 M, 1 M, and 3 M NaOH) at Cs⁺ concentrations selected to probe a range of adsorption densities. Model simulations of Cs⁺ sorption to the sediment in the presence of 0.1 M NaOH for 112 days slightly under-predicted sorption at the lower Cs⁺ adsorption densities. At the higher adsorption densities, model simulations under-predicted sorption by 57%. This under-prediction was surmised to be the result of tetranatrolite precipitation, and subsequent slow Na → Cs exchange. At higher OH concentrations, Cs⁺ sorption in the presence of base for 112 days was unexpectedly equal to, or greater than that expected for pristine sediment. The precipitation of secondary phases, coupled with the fairly unique mica distribution and quantity across all size-fractions in the Hanford sediment, appears to mitigate the impact of base dissolution on Cs⁺ sorption.     

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

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

Solubility of cyclooctasulfur in pure water and sea water at different temperatures

The solubility of cyclooctasulfur in water and sea water at various temperatures in the range between 4 and 80 °C was determined. Cyclooctasulfur in equilibrium with rhombic sulfur reacted with hot acidic aqueous potassium cyanide to form thiocyanate anion which was measured by anion chromatography. Sulfur solubility in pure water was found to increase with temperature by more than 78 times: from 6.1 nM S₈ at 4 °C to 478 nM S₈ at 80 °C. The following thermodynamic values for solubilisation of S₈ in water were calculated from the experimental data: K° = 3.01 ± 1.04 × 10⁻⁸, ΔG_r° = 42.93 ± 0.73 kJ mol⁻¹, ΔH_r° = 47.4 ± 3.6 kJmol⁻¹, ΔS_r° = 15.0 ± 11.7 J mol⁻¹ K⁻¹). Solubility of cyclooctasulfur in sea water was found to be 61 ± 13% of the solubility in pure water regardless of the temperature.     

Energetics of multicomponent diffusion in molten CaO-Al2O3-SiO2

The energetics of multicomponent diffusion in molten CaO-Al₂O₃-SiO₂ (CAS) were examined experimentally at 1440 to 1650°C and 0.5 to 2 GPa. Two melt compositions were investigated: a haplodacitic melt (25 wt.% CaO, 15% Al₂O₃, and 60% SiO₂) and a haplobasaltic melt (35% CaO, 20% Al₂O₃, and 45% SiO₂). Diffusion matrices were measured in a mass-fixed frame of reference with simple oxides as end-member components and Al₂O₃ as a dependent variable. Chemical diffusion in molten CAS shows clear evidence of diffusive coupling among the components. The diffusive flux of SiO₂ is significantly enhanced whenever there is a large CaO gradient that is oriented in a direction opposite to the SiO₂ gradient. This coupling effect is more pronounced in the haplodacitic melt and is likely to be significant in natural magmas of rhyolitic to andesitic compositions. The relative magnitude of coupled chemical diffusion is not very sensitive to changes in temperature and pressure.To a good approximation, the measured diffusion matrices follow well-defined Arrhenius relationships with pressure and reciprocal temperature. Typically, a change in temperature of 100°C results in a relative change in the elements of diffusion matrix of 50 to 100%, whereas a change in pressure of 1 GPa introduces a relative change in elements of diffusion matrix of 4 to 6% for the haplobasalt, and less than 5% for the haplodacite. At a pressure of 1 GPa, the ratios between the major and minor eigenvalues of the diffusion matrix λ₁/λ₂ are not very sensitive to temperature variations, with an average of 5.5 ± 0.2 for the haplobasalt and 3.7 ± 0.6 for the haplodacite. The activation energies for the major and minor eigenvalues of the diffusion matrix are 215 ± 12 and 240 ± 21 kJ mol⁻¹, respectively, for the haplodacite and 192 ± 8 and 217 ± 14 kJ mol⁻¹ for the haplobasalt. These values are comparable to the activation energies for self-diffusion of calcium and silicon at the same melt compositions and pressure. At a fixed temperature of 1500°C, the ratios λ₁/λ₂ increase with the increase of pressure, with λ₁/λ₂ varying from 2.5 to 4.1 (0.5 to 1.3 GPa) for the haplodacite and 4 to 6.5 (0.5 to 2.0 GPa) for the haplobasalt. The activation volumes for the major and minor eigenvalues of the diffusion matrix are 0.31 ± 0.44 and 2.3 ± 0.8 cm³ mol⁻¹, respectively, for the haplodacite and −1.48 ± 0.18 and −0.42 ± 0.24 cm³ mol⁻¹ for the haplobasalt. These values are quite different from the activation volumes for self-diffusion of calcium and silicon at the same melt compositions and temperature. These differences in activation volumes between the two melts likely result from a difference in the structure and thermodynamic properties of the melt between the two compositions (e.g., partial molar volume).Applications of the measured diffusion matrices to quartz crystal dissolution in molten CAS reveal that the activation energy and activation volume for quartz dissolution are almost identical to the activation energy and activation volume for diffusion of the minor or slower eigencomponent of the diffusion matrix. This suggests that the diffusion rate of slow eigencomponent is the rate-limiting factor in isothermal crystal dissolution, a conclusion that is likely to be valid for crystal growth and dissolution in natural magmas when diffusion in liquid is the rate-limiting factor.     

A Donnan potential model for metal sorption onto Bacillus subtilis

In this study, we conducted electrophoretic mobility, potentiometric titration, and metal sorption experiments to investigate the surface charge characteristics of Bacillus subtilis and the electrostatic interactions between metal cations and the cell surface electric field. Electrophoretic mobility experiments performed as a function of pH and ionic strength show an isoelectric point of pH 2.4, with the magnitude of the electrokinetic potential increasing with increasing pH, and decreasing with increasing ionic strength. Potentiometric titration experiments conducted from pH 2.4 to 9 yield an average surface charge excess of 1.6 μmol/mg (dry mass). Corresponding cell wall charge density values were used to calculate the Donnan potential (Ψ_DON) as function of pH and ionic strength. Metal sorption experiments conducted with Ca(II), Sr(II), and Ba(II) exhibit strong ionic strength dependence, suggesting that the metal ions are bound to the bacterial cell wall via an outer-sphere complexation mechanism. Intrinsic metal sorption constants for the sorption reactions were determined by correcting the apparent sorption constant with the Boltzmann factor. A 1:2 metal-ligand stoichiometry provides the best fit to the experimental data with log K₂^int values of 5.9 ± 0.3, 6.0 ± 0.2, 6.2 ± 0.2 for Ca(II), Sr(II), and Ba(II) respectively. Electrophoretic mobility measurements of cells sorbed with Ca(II), Sr(II), and Ba(II) support the 1:2 sorption stoichiometry. These results indicate that electrical potential parameters derived from the Donnan model can be applied to predict metal binding onto bacterial surfaces over a wide range of pH and ionic strength conditions.     

Diffusion and sorption of Cs, I and HTO in samples of the argillaceous Wakkanai Formation from the Horonobe URL, Japan: Clay-based modeling approach

Diffusion and sorption behaviors of cationic Cs⁺, anionic I⁻ and neutral HTO in samples of the Wakkanai Formation from the Horonobe underground research laboratory (URL), Japan, were investigated as a function of ionic strength (I) of synthetic groundwater by through-diffusion and batch sorption experiments and mechanistic modeling. The effective diffusivities (D_e) measured by through-diffusion experiments showed cation excess and anion exclusion effects, which were strongly dependent on I; D_e for Cs⁺ decreased as I increased, D_e for I⁻ showed the opposite dependency and D_e for HTO showed no dependence. The sorption of Cs⁺ measured by through-diffusion and batch sorption experiments were described by Freundlich isotherms with consistent parameters and decreased with I as a result of competitive ion exchange.Diffusion and sorption behaviors were interpreted by assuming the clay components of illite and smectite control diffusion and sorption mechanisms. The component additive (CA) sorption model, which includes illite and smectite contents and their ion exchange constants, provided a reasonable account of the Cs⁺ sorption trends measured as functions of I and Cs concentration. The diffusion model was developed by coupling the electrical double layer (EDL) model, describing the change of ionic concentrations (cation excess and anion deficit) and viscoelectric effects caused by electrostatic interaction at negatively charged clay surfaces, and a simplified pore model assuming one type of pore shape and includes their size distribution. When averaging the electrostatic effects by using the pore surface area distribution, the model could predict the cation excess and anion exclusion effects, and its dependence on I reasonably well. This result implies the nanoscale pores dominating the pore surface area can strongly impact on ionic diffusion in argillaceous rocks. The clay-based modeling approach described here provides a useful tool to predict ionic diffusion and sorption in argillaceous rocks.     

Sorption of Ni on Na-rectorite studied by batch and spectroscopy methods

Sorption of Ni²⁺ on Na-rectorite as a function of contact time, temperature, pH and fulvic acid (FA)/humic acid (HA) was studied under ambient conditions. A pseudo-second-order rate equation was used to simulate the kinetic sorption. The removal of Ni²⁺ increased with increasing pH. The presence of FA/HA enhanced the sorption of Ni²⁺ at low pH values, whereas no drastic effect of FA/HA on Ni²⁺ uptake to rectorite was found at high pH values. The diffuse layer model (DLM) fitted the experimental data of Ni²⁺ sorption in the absence and presence of FA/HA very well with the aid of FITEQL 3.2. The Langmuir, Freundlich and Dubinin–Radushkevich (D–R) models were used to simulate the sorption isotherms of Ni²⁺ at different temperatures. The thermodynamic data (ΔH⁰, ΔS⁰, ΔG⁰) were calculated from the temperature dependent sorption isotherms and the results suggested that the sorption process of Ni²⁺ on rectorite was spontaneous and endothermic. The sorption and species of Ni²⁺ on rectorite in the presence and absence of FA/HA was also investigated and characterized by XPS. The spectroscopic analysis indicated no drastic structural changes of Na-rectorite and the sorption of Ni²⁺ mainly occurred on the surface and at the edge position of Na-rectorite.     

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

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

Uranyl adsorption onto montmorillonite: Evaluation of binding sites and carbonate complexation

The fate and transport of uranium in contaminated soils and sediments may be affected by adsorption onto the surface of minerals such as montmorillonite. Extended X-ray absorption fine structure (EXAFS) spectroscopy has been used to investigate the adsorption of uranyl (UO₂²⁺) onto Wyoming montmorillonite. At low pH (∼4) and low ionic strength (10⁻³ M), uranyl has an EXAFS spectrum indistinguishable from the aqueous uranyl cation, indicating binding via cation exchange. At near-neutral pH (∼7) and high ionic strength (1 M), the equatorial oxygen shell of uranyl is split, indicating inner-sphere binding to edge sites. Linear-combination fitting of the spectra of samples reacted under conditions where both types of binding are possible reveals that cation exchange at low ionic strengths on SWy-2 may be more important than predicted by past surface complexation models of U(VI) adsorption on related montmorillonites. Analysis of the binding site on the edges of montmorillonite suggests that U(VI) sorbs preferentially to [Fe(O,OH)₆] octahedral sites over [Al(O,OH)₆] sites. When bound to edge sites, U(VI) occurs as uranyl-carbonato ternary surface complexes in systems equilibrated with atmospheric CO₂. Polymeric surface complexes were not observed under any of the conditions studied. Current surface complexation models of uranyl sorption on clay minerals may need to be reevaluated to account for the possible increased importance of cation exchange reactions at low ionic strengths, the presence of reactive octahedral iron surface sites, and the formation of uranyl-carbonato ternary surface complexes. Considering the adsorption mechanisms observed in this study, future studies of U(VI) transport in the environment should consider how uranium retardation will be affected by changes in key solution parameters, such as pH, ionic strength, exchangeable cation composition, and the presence or absence of CO₂.     

Structural evidence for the sorption of Ni(II) atoms on the edges of montmorillonite clay minerals: a polarized X-ray absorption fine structure study

The nature of surface complexes formed on Ni uptake onto montmorillonite (a dioctahedral smectite) has been investigated over an extended time period by polarized extended X-ray absorption fine structure (P-EXAFS) spectroscopy. Self-supporting films of Ni-sorbed montmorillonite were prepared by contacting Ni and montmorillonite at pH 7.2, high ionic strength (0.3 M NaClO₄), and low Ni concentration ([Ni]_initial = 19.9 μM) for 14- and 360-d reaction time. The resulting Ni concentration on the clay varied from 4 to 7 μmol/g. Quantitative texture analysis indicates that the montmorillonite particles were well orientated with respect to the plane of the film. The full width at half maximum of the orientation distribution of the c* axes of individual clay platelets about the normal to the film plane was 44.3° (14-d reaction time) and 47.1° (360-d reaction time). These values were used to correct the coordination numbers determined by P-EXAFS for texture effects. Ni K-edge P-EXAFS spectra were recorded at angles between the incident beam and the film normal equal to 10, 35, 55, and 80°. Spectral analysis led to the identification of three nearest cationic subshells containing 2.0 ± 0.5 Al at 3.0 Å and 2.0 ± 0.5 Si at 3.12 Å and 4.0 ± 0.5 Si at 3.26 Å. These distances are characteristic of edge-sharing linkages between Al and Ni octahedra and of corner-sharing linkages between Ni octahedra and Si tetrahedra, as in clay structures. The angular dependence of the Ni-Al and Ni-Si contributions indicates that Ni-Al pairs are oriented parallel to the film plane, whereas Ni-Si pairs are not. The study reveals the formation of Ni inner-sphere mononuclear surface complexes located at the edges of montmorillonite platelets and thus that heavy metals binding to edge sites is a possible sorption mechanism for dioctahedral smectites. Data analysis further suggests that either the number of neighboring Al atoms slightly increases from 1.6 to 2 or that the structural order of the observed surface complexes increases from 0.01 Ų to 0.005 Ų with increasing reaction time. On the basis of the low Ni-Al coordination numbers, it appears that over an extended reaction time period of 1 yr the diffusion of Ni atoms in the octahedral layer is not the major uptake mechanism of Ni onto montmorillonite.