Structure and dynamics of the water films confined between edges of pyrophyllite: A first principle study

Edge sites of clay minerals play a key role for pH dependent sorption of ions from solutions of electrolytes. Pyrophyllite, Al₂[Si₄O₁₀](OH)₂, is an important structural prototype for a variety of 2:1 dioctahedral phyllosilicates but in contrast to the other clays has no permanent structural charge. The structure of thin water films confined between most common edges of 1Tc pyrophyllite: (0 1 0), (1 1 0) and (1 0 0), was analyzed by means of ab initio molecular dynamic simulations. The system setup allowed for a full flexibility of the interfaces and a proton exchange between the edges of pyrophyllite and water molecules in solution. The structure of hydrated surfaces is compared with the recent predictions of static geometry optimizations for edge-vacuum interfaces. All surfaces studied reveal a strong hydrophilic character of edge similar to the hydrated silica surface and the facets of simple layered hydroxides. Spontaneous proton transfer between different surface sites were observed in molecular dynamics simulations of the (0 1 0) interface. The proton bound to the SiOH site was found to exchange with the AlOH group by the mechanism . The direction of the proton transfer agrees with the scale of relative proton affinities for surface sites obtained from the static calculations. Alternatively, the proton attached to the AlOH₂ site exchanges with the AlOH group. In both reactions, the protons are transferred through the chains of hydrogen bonds formed between water molecules in the solution and the surface sites. The observed mechanisms might be one of the basic schemes for the surface proton diffusion in compacted clays. Kinetics of the proton transfer at edge sites is limited by the rate of rearrangements of the water molecules near interface.     

Weak interfacial water ordering on isostructural hematite and corundum (0 0 1) surfaces

Ordering of interfacial water at the hematite and corundum (0 0 1)-water interfaces has been characterized using in situ high resolution specular X-ray reflectivity measurements. The hematite (0 0 1) surface was prepared through an annealing process to produce a surface isostructural with corundum (0 0 1), facilitating direct comparison. Interfacial water was found to display a similar structure on this pair of isostructural surfaces. A single layer of adsorbed water having a large vibrational amplitude was present on each surface and additional ordering of water extended at least 1 nm into the bulk fluid, with the degree of ordering decreasing with increasing distance from the surfaces. Consistent with prior studies of the (0 1 2) and (1 1 0) surfaces of hematite and corundum, the configuration of water above the (0 0 1) surfaces is primarily controlled by the surface structure, specifically the arrangement of surface functional groups. However, interfacial water at the (0 0 1) surfaces displayed significantly larger vibrational amplitudes throughout the interfacial region than at other isostructural sets of hematite and corundum surfaces, indicating weaker ordering. Comparison of the vibrational amplitudes of adsorbed water on a series of oxide, silicate, and phosphate mineral surfaces suggests that the presence or absence of a substantial interfacial electrostatic field is the primary control on water ordering and not the surface structure itself. On surfaces for which charge originates dominantly through protonation-deprotonation reactions the controlling factor appears to be whether conditions exist where most functional groups are uncharged as opposed to the net surface charge. The doubly coordinated functional groups on hematite and corundum (0 0 1) surfaces are largely uncharged under slightly acidic to circumneutral pH conditions, leading to weak ordering, whereas singly coordinated groups on (0 1 2) and (1 1 0) surfaces of these phases are always charged, even when the net surface charge is zero, and induce strong water ordering. Surfaces lacking structural charge can thus be divided into two distinct classes that induce either strong or weak ordering of interfacial water. Surface functional group coordination is the ultimate control on this division as it determines the charge state of such groups under different protonation configurations. Ion adsorption and electron transfer processes may differ between these classes of surfaces because of the effect of water ordering strength on interfacial capacitances and hydrogen bonding.     

Hydrated goethite (α-FeOOH) (1 0 0) interface structure: Ordered water and surface functional groups

Goethite(α-FeOOH), an abundant and highly reactive iron oxyhydroxide mineral, has been the subject of numerous studies of environmental interface reactivity. However, such studies have been hampered by the lack of experimental constraints on aqueous interface structure, and especially of the surface water molecular arrangements. Structural information of this type is crucial because reactivity is dictated by the nature of the surface functional groups and the structure or distribution of water and electrolyte at the solid-solution interface. In this study we have investigated the goethite (1 0 0) surface using surface diffraction techniques, and have determined the relaxed surface structure, the surface functional groups, and the three dimensional nature of two distinct sorbed water layers. The crystal truncation rod (CTR) results show that the interface structure consists of a double hydroxyl, double water terminated interface with significant atom relaxations. Further, the double hydroxyl terminated surface dominates with an 89% contribution having a chiral subdomain structure on the (1 0 0) cleavage faces. The proposed interface stoichiometry is ((H₂O)(H₂O)OH₂OHFeOOFeR) with two types of terminal hydroxyls; a bidentate (B-type) hydroxo group and a monodentate (A-type) aquo group. Using the bond-valence approach the protonation states of the terminal hydroxyls are predicted to be OH type (bidentate hydroxyl with oxygen coupled to two Fe³⁺ ions) and OH₂ type (monodentate hydroxyl with oxygen tied to only one Fe³⁺). A double layer three dimensional ordered water structure at the interface was determined from refinement of fits to the experimental data. Application of bond-valence constraints to the terminal hydroxyls with appropriate rotation of the water dipole moments allowed a plausible dipole orientation model as predicted. The structural results are discussed in terms of protonation and H-bonding at the interface, and the results provide an ideal basis for testing theoretical predictions of characteristic surface properties such as pK_a , sorption equilibria, and surface water permittivity.     

XPS spectra of uranyl minerals and synthetic uranyl compounds. II: The O 1s spectrum

The O 1s spectrum is examined for 19 uranyl minerals of different composition and structure. Spectra from single crystals were measured with a Kratos Axis Ultra X-ray Photoelectron Spectrometer with a magnetic-confinement charge-compensation system. Well-resolved spectra with distinct maxima, shoulders and inflections points, in combination with reported and measured binding energies for specific O²⁻ species and structural data of the uranyl minerals are used to resolve the fine structure of the O 1s envelope. The resolution of the O 1s spectra includes, for the first time, different O²⁻ bands, which are assigned to O atoms linking uranyl with uranyl polyhedra (UOU) and O atoms of uranyl groups (OUO). The resolved bands in the O 1s spectrum occur at distinct ranges in binding energy: bands for (UOU) occur at 529.6-530.4 eV, bands for (OUO) occur at 530.6-531.4 eV, bands for O²⁻ in the equatorial plane of the uranyl polyhedra linking uranyl polyhedra with (TO_n) groups (T = Si, S, C, P, Se) (TO) occur at 530.9-532.2 eV; bands for (OH) groups in the equatorial plane of the uranyl polyhedra (OH) occur at 532.0-532.5 eV, bands of (H₂O) groups in the interstitial complex of the uranyl minerals (H₂O_interst) occur at 533.0-533.8 eV and bands of physisorbed (H₂O) groups on the surface of uranyl minerals (H₂O_adsorb) occur at 534.8-535.2 eV. Treatment of uranyl minerals with acidic solutions results in a decrease in UOU and an increase in OH. Differences in the ratio of OH : OUO between the surface and bulk structure is larger for uranyl minerals with a high number of UOU and TO species in the bulk structure which is explained by protonation of underbonded UO, UOU and TO terminations on the surface. The difference in the ratio of H₂O_interst : OUO between the bulk and surface structures is larger for uranyl minerals with higher percentages of H₂O_interst as well as, with a higher number of interstitial H₂O groups that are not bonded to interstitial cations, resulting in easier dehydration of interstitial H₂O groups in uranyl minerals during exposure to a vacuum.     

Adsorption of Rb and Sr at the orthoclase (0 0 1)-solution interface

Adsorption of Rb⁺ and Sr²⁺ at the orthoclase (0 0 1)-solution interface is probed with high-resolution X-ray reflectivity and resonant anomalous X-ray reflectivity. Specular X-ray reflectivity data for orthoclase in contact with 0.01 m RbCl solution at pH 5.5 reveal a systematic increase in electron density adjacent to the mineral surface with respect to that observed in contact with de-ionized water (DIW). Quantitative analysis indicates that Rb⁺ adsorbs at a height of 0.83 ± 0.03 Å with respect to the bulk K⁺ site with a nominal coverage of 0.72 ± 0.10 ions per surface unit mesh (55.7 Å²). These results are consistent with an ion-exchange reaction in which Rb⁺ occupies an inner-sphere adsorption (IS) site. In contrast, X-ray reflectivity data for orthoclase in contact with 0.01 m Sr(NO₃)₂ solution at pH 5.3 reveal few significant changes with respect to DIW. Resonant anomalous X-ray reflectivity was used to probe Sr²⁺ adsorption and to image its vertical distribution. This element-specific measurement reveals that Sr²⁺ adsorbs with a total coverage of 0.37 ± 0.02 ions per surface unit mesh, at a substantially larger height (3.28 ± 0.05 Å) than found for Rb⁺, and with a relatively broad density distribution (having a root-mean-square width of 1.88 ± 0.08 Å for a single-peak model), implying that Sr²⁺ adsorbs primarily as a fully-hydrated outer-sphere (OS), species. Comparison to a two-height model suggests that 13 ± 5% of the adsorbed Sr²⁺ may be present as an IS species. This partitioning implies a ∼5 kJ/mol difference in free energy between the IS and OS Sr²⁺ on orthoclase. Differences in the partitioning of Sr²⁺ between IS and OS species for orthoclase (0 0 1) and muscovite (0 0 1) suggest control by the geometry of the IS adsorption site. Results for the OS distribution are compared to predictions of the Poisson-Boltzmann equation in the strong coupling regime, which predicts an intrinsically narrow vertical diffuse ion distribution; the OS distribution might thus be thought of as the diffuse ion profile in the limit of high surface charge.     

Interfacial water structure on the (0 1 2) surface of hematite: Ordering and reactivity in comparison with corundum

Many geochemical reactions that control the composition of natural waters, contaminant fate and transport, and biogeochemical element cycling take place at the interface between minerals and aqueous solutions. A fundamental understanding of these important processes requires knowledge of the structure of mineral-water interfaces. High-resolution specular X-ray reflectivity was used to determine the structure of the hematite (0 1 2)-water interface. Relaxation of the surface was observed to be minor, and water was found to order near the hematite surface. Two sites of adsorbed water are inferred to be ordered laterally, one bridging between triply coordinated functional groups and the other bridging between the singly coordinated functional groups on the surface, as steric constraints limit the possible arrangements of water molecules occurring at the observed heights above the hematite surface. Relaxations of the hematite and corundum (0 1 2) surfaces, which are isostructural, are similar and limited primarily to the top most layer of the structures. No significant changes to the interfacial stoichiometry (i.e., partial occupancy of surface species) are observed in either case. The structure of interfacial water is similar on the hematite and corundum (0 1 2) surfaces as well, although water appeared to be less well ordered on the hematite surface. This may be due to expected differences in the oxygen exchange rates from surface functional groups or the apparent better matching of the corundum oxygen lattice to the natural structural ordering in water, and suggests that the dielectric constant gradients of interfacial water may differ on the two surfaces. Similar charging behavior is expected for these surfaces as similar types of surface functional groups are exposed. Although generally similar, subtle differences in the reactivity of hematite and corundum (0 1 2) surfaces to arsenate adsorption, and possibly the adsorption of other species, may be related to the difference in ordering of interfacial water observed in this study.     

Water ordering and surface relaxations at the hematite (1 1 0)-water interface

Structural characterization of iron oxide-water interfaces provides insight into the mechanisms through which these minerals control contaminant fate and element cycling in soil, sedimentary, and groundwater systems. Ordering of interfacial water and structural relaxations at the hematite (1 1 0) surface have been investigated in situ using high-resolution specular X-ray reflectivity. These measurements demonstrate that relaxations are constrained to primarily the top ∼5 Å of the surface. Near-surface iron atoms do not relax substantially, although the uppermost layer displays an increased distribution width, while the undercoordinated oxygens on the surface uniformly relaxed outward. Two sites of adsorbed water and additional layering of water farther from the surface were observed. Water fully covers the (1 1 0) surface and appears to form a continuous network extending into bulk solution, with positional order decreasing to that of a disordered bulk fluid within 1 nm. The arrangement of water is similar to that on the hematite (0 1 2) surface, which has a similar surface topography, although these surfaces display different vibrational amplitudes or positional disorder of adsorbed water molecules and average spacings of near-surface layered water. Comparison between these surfaces suggests that interfacial water ordering on hematite is controlled primarily by surface structure and steric constraints and that highly ordered water is likely common to most hematite-water interfaces.     

Solution behavior of reduced COH volatiles in silicate melts at high pressure and temperature

The solubility and solution mechanisms of reduced COH volatiles in Na₂OSiO₂ melts in equilibrium with a (H₂ + CH₄) fluid at the hydrogen fugacity defined by the iron-wüstite + H₂O buffer [f_H2(IW)] have been determined as a function of pressure (1-2.5 GPa) and silicate melt polymerization (NBO/Si: nonbridging oxygen per silicon) at 1400 °C. The solubility, calculated as CH₄, increases from ∼0.2 wt% to ∼0.5 wt% in the melt NBO/Si-range ∼0.4 to ∼1.0. The solubility is not significantly pressure-dependent, probably because f_H2(IW) in the 1-2.5 GPa range does not vary greatly with pressure. Carbon isotope fractionation between methane-saturated melts and (H₂ + CH₄) fluid varied by ∼14‰ in the NBO/Si-range of these melts.The (C..H) and (O..H) speciation in the quenched melts was determined with Raman and ¹H MAS NMR spectroscopy. The dominant (C..H)-bearing complexes are molecular methane, CH₄, and a complex or functional group that includes entities with CCH bonding. Minor abundance of complexes that include SiOCH₃ bonding is tentatively identified in some melts. There is no spectroscopic evidence for SiC or SiCH₃. Raman spectra indicate silicate melt depolymerization (increasing NBO/Si). The [CH₄/CCH]^melt abundance ratio is positively correlated with NBO/Si, which is interpreted to suggest that the (CCH)-containing structural entity is bonded to the silicate melt network structure via its nonbridging oxygen. The ∼14‰ carbon isotope fractionation change between fluid and melt is because of the speciation changes of carbon in the melt.     

Fe(II) adsorption on hematite (0 0 0 1)

The surface structure of α-Fe₂O₃(0 0 0 1) was studied using crystal truncation rod (CTR) X-ray diffraction before and after reaction with aqueous Fe(II) at pH 5. The CTR results show the unreacted α-Fe₂O₃(0 0 0 1) surface consists of two chemically distinct structural domains: an O-layer terminated domain and a hydroxylated Fe-layer terminated domain. After exposing the α-Fe₂O₃(0 0 0 1) surface to aqueous Fe(II), the surface structure of both co-existing structural domains was modified due to adsorption of Fe at crystallographic lattice sites of the substrate, resulting in six-coordinated adsorbed Fe at the surface. The average Fe-O bond lengths of the adsorbed Fe are consistent with typical Fe(III)-O bond lengths (in octahedral coordination), providing evidence for the oxidation of Fe(II) to Fe(III) upon adsorption. These results highlight the important role of substrate surface structure in controlling Fe(II) adsorption. Furthermore, the molecular scale structural characterization of adsorbed Fe provides insight into the process of Fe(II) induced structural modification of hematite surfaces, which in turn aids in assessing the effective reactivity of hematite surfaces in Fe(II) rich environments.     

A combined FTIR and TPD study on the bulk and surface dehydroxylation and decarbonation of synthetic goethite

The thermal dehydroxylation of a goethite-carbonate solid solution was studied with combined Fourier-transform infrared (FTIR)-Temperature programmed desorption (TPD) experiments. The TPD data revealed dehydroxylation peaks involving the intrinsic dehydroxylation of goethite at 560 K and a low temperature peak at 485 K which was shown to be associated to the release of non-stoichiometric water from the goethite bulk and surface. The FTIR and the TPD data of goethite in the absence of adsorbed carbonate species revealed the presence of adventitious carbonate mostly sequestered in the goethite bulk. The release of carbonate was however not only related to the dehydration of goethite but also from the crystallization of hematite at temperatures exceeding 600 K. The relative abundance of surface hydroxyls was shown to change systematically upon goethite dehydroxylation with a preferential stripping of singly-coordinated OH sites followed by a dramatic change in the dominance of the different surface hydroxyls upon the formation of hematite.     

Computer simulations of the interactions of the (0 1 2) and (0 0 1) surfaces of jarosite with Al, Cd, Cu and Zn

Jarosite is an important mineral on Earth, and possibly on Mars, where it controls the mobility of iron, sulfate and potentially toxic metals. Atomistic simulations have been used to study the incorporation of Al³⁺, and the M²⁺ impurities Cd, Cu and Zn, in the (0 1 2) and (0 0 1) surfaces of jarosite. The calculations show that the incorporation of Al on an Fe site is favorable on all surfaces in which terminal Fe ions are exposed, and especially on the (0 0 1) [Fe₃(OH)₃]⁶⁺ surface. Incorporation of Cd, Cu or Zn on a K site balanced by a K vacancy is predicted to stabilize the surfaces, but calculated endothermic solution energies and the high degree of distortion of the surfaces following incorporation suggest that these substitutions will be limited. The calculations also suggest that incorporation of Cd, Cu and Zn on an Fe site balanced by an OH vacancy, or by coupled substitution on both K and Fe sites, is unfavorable, although this might be compensated for by growth of a new layer of jarosite or goethite, as predicted for bulk jarosite. The results of the simulations show that surface structure will exert an influence on uptake of impurities in the order Cu > Cd > Zn, with the most favorable surfaces for incorporation being (0 1 2) [KFe(OH)₄]⁰ and (0 0 1) [Fe₃(OH)₃]⁶⁺.     

Bridging arsenate surface complexes on the hematite (0 1 2) surface

The fate of the oxoanion arsenate in diverse systems is strongly affected by its adsorption on the surfaces of iron (oxyhydr)oxide minerals. Predicting this behavior in the environment requires an understanding of the mechanisms of arsenate adsorption. In this study, the binding site and adsorption geometry of arsenate on the hematite (0 1 2) surface is investigated. The structure and termination of the hematite (0 1 2)-water interface were determined by high resolution X-ray reflectivity, revealing that two distinct terminations exist in a roughly 3:1 proportion. The occurrence of multiple terminations appears to be a result of sample preparation, and is not intrinsic to the hematite (0 1 2) surface. X-ray standing wave (XSW) measurements were used to determine the registry of adsorbed arsenate to the hematite structure, and thus the binding site and geometry of the resulting surface complex. Arsenate forms a bridging bidentate complex on two adjacent singly coordinated oxygen groups on each of the two distinct terminations present at the hematite surface. Although this geometry is consistent with that seen in past studies, the derived As-Fe distances are longer, the result of the topology of the FeO₆ octahedra on the (0 1 2) surface. As EXAFS-derived As-Fe distances are often used to determine the adsorption mechanism in environmental samples (e.g., mine tailings, contaminated sediments), this demonstrates the importance of considering the possible sorbent surface structures and arrangements of adsorbates when interpreting such data.As multiple functional groups are present and multiple binding geometries are possible on the hematite (0 1 2) surface, the XSW data suggest that formation of bridging bidentate surface complexes on singly coordinated oxygen sites is the preferred adsorption mechanism on this and most other hematite surfaces (provided those surfaces contain adjacent singly coordinated oxygen groups). These measurements also constrain the likely reaction stoichiometry, with only the protonation state of the surface complex undetermined. Although bridging bidentate inner-sphere surface complexes comprised the majority of the adsorbed arsenate present on the hematite (0 1 2) surface, there is an additional population of sorbed arsenate species that could not be characterized by the XSW measurements. These species are likely more disordered, and thus more weakly bound, than the bridging bidentate complexes, and may play a role in determining the fate, transport, and bioavailability of arsenate in the environment. Finally, the possibility of obtaining species-specific XSW measurements by tuning the incident beam energy to specific features in a XANES spectrum is described.     

In situ Atomic Force Microscopy study of dissolution of the barite (0 0 1) surface in water at 30 °C

The dissolution behavior of the barite (0 0 1) surface in pure water at 30 °C was investigated using in situ Atomic Force Microscopy (AFM), to better understand the dissolution mechanism and the microtopographical changes that occur during the dissolution, such as steps and etch pits. The dissolution of the barite (0 0 1) surface started with the slow retreat of steps, after which, about 60 min later, the <hk0> steps of one unit cell layer or multi-layers became two-step fronts (fast “f” and slow “s” steps) with a half-unit cell layer showing different retreat rates. The “f” step had a fast retreat rate (≈(14 ± 1) × 10⁻² nm/s) and tended to have a jagged step edge, whereas the “s” step (≈(1.8 ± 0.1) × 10⁻² nm/s) had a relatively straight front. The formation of the “f” steps led to the formation of a new one-layer step, where the front of the “s” step was overtaken by that of the immediate underlying “f” step. The “f” steps also led to the decrease of the <hk0> steps and the increase in the percentage of stable steps parallel to the [0 1 0] direction during the dissolution.Etch pits, which could be observed after about 90 min, were of three types: triangular etch pits with a depth of a half-unit cell, shallow etch pits, and deep etch pits. The triangular etch pits were bounded by the step edges parallel to [0 1 0], [1 2 0], and [] and had opposite orientations in the upper half and lower half layers. Shallow etch pits that had a depth of two or more half-unit cell layers had any two consecutive pits pointing in the opposite direction of each other. The triangular etch pit appeared to grow by simultaneously removal of a row of ions parallel to each direction from the three step edges. At first, deep etch pits were elongated in the [0 1 0] direction with a curved outline and then gradually developed to an angular form bounded by the {1 0 0}, {3 1 0}, and (0 0 1) faces. The retreat rate of the (0 0 1) face was much slower than those of the {1 0 0} and {3 1 0} and tended to separate into two rates ((0.13 ± 0.01) × 10⁻² nm/s for the deep etch pits derived from a screw dislocation and (0.07 ± 0.01) × 10⁻² nm/s for those from other line defects).The changes in the dissolution rate of a barite (0 0 1) surface during the dissolution were estimated using the retreat rates and densities of the various steps as well as the growth rates, density, and areas of the lateral faces of the deep etch pits that were obtained from this AFM analysis. Our results showed that the dissolution rate of the barite (0 0 1) surface gradually increased and approached the bulk dissolution rate because of the change in the main factor determining the dissolution rate from the density of the steps to the growth and the density of the deep etch pits on the surface.     

Interaction of muscovite (0 0 1) with Pu bearing solutions at pH 3 through ex-situ observations

The interaction of Pu³⁺ bearing solutions with the muscovite (0 0 1) basal plane is explored using a combination of ex-situ approaches including alpha-counting, to determine the Pu³⁺ adsorption isotherm, and X-ray reflectivity (XR) and resonant anomalous X-ray reflectivity (RAXR), to probe the interfacial structure and Pu-specific distribution, respectively. Pu uptake to the muscovite (0 0 1) surface from Pu³⁺ solutions in a 0.1 M NaClO₄ background electrolyte at pH 3 follows an approximate Langmuir isotherm with an apparent adsorption constant, K_app = 5 × 10⁴ M⁻¹, and with a maximum coverage that is consistent with the amount needed to fully compensate the surface charge by trivalent Pu. The XR results show that the muscovite surface reacted with a 10⁻³ M Pu³⁺ solution (at pH 3 with 0.1 M NaClO₄) and dried in the ambient environment, maintains a 30-40 Å thick layer, indicating the presence of a residual hydration layer (possibly including adventitious carbon). The RAXR results indicate that Pu sorbs on the muscovite surface with an intrinsically broad distribution with an average height of 18 Å, substantially larger than heights expected for any specifically adsorbed inner- or outer-sphere complexes. These results are discussed in the context of recent studies of cation adsorption trends on muscovite and the possible roles of Pu hydrolysis species in controlling the Pu-muscovite interactions.     

Structure, charge distribution, and electron hopping dynamics in magnetite (Fe3O4) (1 0 0) surfaces from first principles

For the purpose of improving fundamental understanding of the redox reactivity of magnetite, quantum-mechanical calculations were applied to predict Fe²⁺ availability and electron hopping rates at magnetite (1 0 0) surfaces, with and without the presence of adsorbed water. Using a low free energy surface reconstruction (½-Fe_tet layer relaxed into the Fe_oct (1 0 0) plane), the relaxed outermost layer of both the hydrated and vacuum-terminated surfaces were found to be predominantly enriched in Fe²⁺ within the octahedral sublattice, irrespective of the presence of adsorbed water. At room temperature, mobile electrons move through the octahedral sublattice by Fe²⁺-Fe³⁺ valence interchange small polaron hopping, calculated at 10¹⁰-10¹² hops/s for bulk and bulk-like (i.e., near-surface) environments. This process is envisioned to control sustainable overall rates of interfacial redox reactions. These rates decrease by up to three orders of magnitude (10⁹ hops/s) at the (1 0 0) surface, and no significant difference is observed for vacuum-terminated versus hydrated cases. Slower hopping rates at the surface appear to arise primarily from larger reorganization energies associated with octahedral Fe²⁺-Fe³⁺ valence interchange in relaxed surface configurations, and secondarily on local charge distribution patterns surrounding Fe²⁺-Fe³⁺ valence interchange pairs. These results suggest that, with respect to the possibility that the rate and extent of surface redox reactions depend on Fe²⁺ availability and its replenishment rate, bulk electron hopping mobility is an upper-limit for magnetite and slower surface rates may need to be considered as potentially rate-limiting. They also suggest that slower hopping mobilities calculated for surface environments may be amenable to Fe²⁺-Fe³⁺ site discrimination by conventional spectroscopic probes.     

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].     

Structure of the orthoclase (001)- and (010)-water interfaces by high-resolution X-ray reflectivity

High-resolution in situ X-ray specular reflectivity was used to measure the structures of orthoclase (001) and (010) cleavage surfaces in contact with deionized water at 25°C. X-ray reflectivity data demonstrate a high degree of structural similarity between these two orthoclase-water interfaces. Both interfacial structures include cleavage along the plane of minimal bond breakage resulting in surfaces terminated by non-bridging oxygens; structured water within 5 Å of the orthoclase surface (consisting of adsorbed species at the surface and layered water above the surface), with a featureless water profile beyond 5 Å; substitution of outermost K⁺ ions by an oxygen containing species (presumably H₃O⁺); and small structural displacements of the near surface atoms. The interfacial water structure, in comparison with recent results for other mineral-water interfaces, is intermediate between the minimal structure found at calcite-, barite-, and quartz-water interfaces and the more extensive structure found at the muscovite-water interface.     

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.     

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

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