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.     

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.     

Uranyl adsorption at (0 1 0) edge surfaces of kaolinite: A density functional study

We report density functional investigations of kaolinite edge surfaces and uranyl adsorption thereon. Applying periodic slab models, we studied the (0 1 0) surface of kaolinite as an example of kaolinite edge facets which are expected to be highly reactive and to adsorb preferentially metal ions. Among the four terminations of the (0 1 0) surface, we selected the two most likely ones and determined their structures to be affected by solvation. On these modified surfaces, we explored bidentate inner-sphere adsorption complexes of uranyl, at single metal center sites, Al(O,OH), and sites of mixed type, AlOH-SiO. On one of the terminations hydrolysis of uranyl was found to occur. Comparison of key calculated structure parameters with available experimental data suggests an extension of the prevailing interpretation and implies that a set of uranyl complexes may coexist on edge surfaces.     

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.     

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.     

Molecular dynamics simulations of the orthoclase (0 0 1)- and (0 1 0)-water interfaces

Molecular dynamics simulations of water in contact with the (0 0 1) and (0 1 0) surfaces of orthoclase (KAlSi₃O₈) were carried out to investigate the structure and dynamics of the feldspar-water interface, contrast the intrinsic structural properties of the two surfaces, and provide a basis for future work on the diffusion of ions and molecules in microscopic mineral fractures. Electron density profiles were computed from the molecular dynamics trajectories and compared with those derived experimentally from high-resolution X-ray reflectivity measurements by Fenter and co-workers [Fenter P., Cheng L., Park C., Zhang H. and Sturchio N. C. (2003a) Structure of the orthoclase (0 0 1)- and (0 1 0)-water interfaces by high-resolution X-ray reflectivity. Geochim. Cosmochim. Acta67, 4267-4275]. For each surface, three scenarios were considered whereby the interfacial species is potassium, water, or a hydronium ion. Excellent agreement was obtained for the (0 0 1) surface when potassium is the predominant interfacial species; however, some discrepancies in the position of the interfacial peaks were obtained for the (0 1 0) surface. The two surfaces showed similarities in the extent of water ordering at the interface, the activation energies for water and potassium desorption, and the adsorption localization of interfacial species. However, there are also important differences between the two surfaces in the coordination of a given adsorbed species, adsorption site densities, and the propensity for water molecules in surface cavities and those in the first hydration layer to coordinate to surface bridging oxygen atoms. These differences may have implications for the extent of dissolution in the low-pH regime since hydrolysis of Si(Al)OSi(Al) bonds is a major dissolution mechanism.     

Structural characterization of poorly-crystalline scorodite, iron(III)-arsenate co-precipitates and uranium mill neutralized raffinate solids using X-ray absorption fine structure spectroscopy

X-ray absorption fine structure (XAFS) is used to characterize the mineralogy of the iron(III)-arsenate(V) precipitates produced during the raffinate (aqueous effluent) neutralization process at the McClean Lake uranium mill in northern Saskatchewan, Canada. To facilitate the structural characterization of the precipitated solids derived from the neutralized raffinate, a set of reference compounds were synthesized and analyzed. The reference compounds include crystalline scorodite, poorly-crystalline scorodite, iron(III)-arsenate co-precipitates obtained under different pH conditions, and arsenate-adsorbed on goethite. The poorly-crystalline scorodite (prepared at pH 4 with Fe/As = 1) has similar As local structure as that of crystalline scorodite. Both As and Fe K-edge XAFS of poorly-crystalline scorodite yield consistent results on As-Fe (or Fe-As) shell. From As K-edge analysis the As-Fe shell has an inter-atomic distance of 3.33 ± 0.02 Å and coordination number of 3.2; while from Fe K-edge analysis the Fe-As distance and coordination number are 3.31 ± 0.02 Å and 3.8, respectively. These are in contrast with the typical arsenate adsorption on bidentate binuclear sites on goethite surfaces, where the As-Fe distance is 3.26 ± 0.03 Å and coordination number is close to 2. A similar local structure identified in the poorly-crystalline scorodite is also found in co-precipitation solids (Fe(III)/As(V) = 3) when precipitated at the same pH (pH = 4): As-Fe distance 3.30 ± 0.03 Å and coordination number 3.9; while at pH = 8 the co-precipitate has As-Fe distance of 3.27 ± 0.03 Å and coordination number about 2, resembling more closely the adsorption case. The As local structure in the two neutralized raffinate solid series (precipitated at pH values up to 7) closely resembles that in the poorly-crystalline scorodite. All of the raffinate solids have the same As-Fe inter-atomic distance as that in the poorly-crystalline scorodite, and a systematic decrease in the As-Fe coordination is observed when pH is progressively increased; the basic poorly-crystalline scorodite structural feature remains in the raffinate solid up to pH 7.     

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.     

Removal of arsenic from aqueous solution by natural siderite and hematite

Batch and column experiments were conducted to examine the capability of naturally formed hematite and siderite to remove As from drinking water. Results show that both minerals were able to remove As from aqueous solutions, but with different efficiencies. In general, each material removed arsenate much more efficiently than As–DMA (dimethylarsinic acid), with the lowest adsorption efficiency for arsenite. The best removal efficiency for As species was obtained using a hematite, with a grain size range between 0.25 and 0.50 mm. The adsorption capacity for inorganic As(V) reached 202 μg/g. The pH generally had a great impact on the arsenate removal by the Fe minerals studied, while arsenite removal was slightly dependent on the initial pH of between 3 and 10. The presence of phosphate always had a negative effect on arsenate adsorption, due to competitive adsorption between them. A column packed with hematite in the upper half and siderite in the lower half with a grain size range of 0.25–0.5 mm proved to be an efficient reactive filter for the removal of all As species, causing a decrease in As concentration from 500 μg/L (including 200 μg/L As(V) as arsenate, 200 μg/L As(III) as arsenite and 100 μg/L As(V) as DMA) to less than 10 μg/L after 1055 pore volumes of water were filtered at a flow rate of 0.51 mL/min. After 2340 pore volumes passed through the column filter, the total inorganic As in the effluent was less than 5 μg/L. The total As load in the column filter was estimated to be 0.164 mg/g. Results of μ-synchrotron X-ray fluorescence analysis (μ-XRFA) suggest that coatings of fresh Fe(III) oxides, formed on the surface of the siderite grains after two weeks of operation, greatly increased the adsorption capacity of the filling material towards As.     

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.     

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.     

Structure and oxidation state of hematite surfaces reacted with aqueous Fe(II) at acidic and neutral pH

Structural changes and surface oxidation state were examined following the reaction of hematite (0 0 1), (0 1 2), and (1 1 0) with aqueous Fe(II). X-ray reflectivity measurements indicated that Fe(II) induces changes in the structure of all three surfaces under both acidic (pH 3) and neutral (pH 7) conditions. The structural changes were generally independent of pH although the extent of surface transformation varied slightly between acidic and neutral conditions; no systematic trends with pH were observed. Induced changes on the (1 1 0) and (0 1 2) surfaces include the addition or removal of partial surface layers consistent with either growth or dissolution. In contrast, a <1 nm thick, discontinuous film formed on the (0 0 1) surface that appears to be epitaxial yet is not a perfect extension of the underlying hematite lattice, being either structurally defective, compositionally distinct, or nanoscale in size and highly relaxed. Resonant anomalous X-ray reflectivity measurements determined that the surface concentration of Fe(II) present after reaction at pH 7 was below the detection limit of approximately 0.5-1 μmol/m² on all surfaces. These observations are consistent with Fe(II) oxidative adsorption, whereby adsorbed Fe(II) is oxidized by structural Fe(III) in the hematite lattice, with the extent of this reaction controlled by surface structure at the atomic scale. The observed surface transformations at pH 3 show that Fe(II) oxidatively adsorbs on hematite surfaces at pH values where little net adsorption occurs, based on historical macroscopic Fe(II) adsorption behavior on fine-grained hematite powders. This suggests that Fe(II) plays a catalytic role, in which an electron from an adsorbed Fe(II) migrates to and reduces a lattice Fe(III) cation elsewhere, which subsequently desorbs in a scenario with zero net reduction and zero net adsorption. Given the general pH-independence and substantial mass transfer involved, this electron and atom exchange process appears to be a significant subsystem within macroscopic pH-dependent Fe(II) adsorption.     

Simultaneous inner- and outer-sphere arsenate adsorption on corundum and hematite

The ability to predict the fate and transport of arsenic in aquatic environments, its impact on water quality and human health, and the performance and cost-effectiveness of water treatment systems relies on understanding how it interacts with solid surfaces. In situ resonant surface X-ray scattering measurements of arsenate adsorption at pH 5 in 0.01 M NaCl on corundum and hematite (012) surfaces demonstrate that arsenate surface complexation is unexpectedly bimodal, adsorbing simultaneously as inner- and outer-sphere species. In addition, this bimodal behavior is found to be independent of the total arsenate solution concentration, and thus surface coverage, over the range of 10⁻⁶ to 10⁻³ M. Alternative mechanisms to produce the observed As distributions, such as arsenate dimerization or surface precipitation of an aluminum or ferric arsenate, are inconsistent with the experimentally-determined total and As-specific density profiles. Based on the location of the outer-sphere arsenate in relation to the surfaces studied, possible binding mechanisms include electrostatic attraction, hydrogen bonding to surface oxygen functional group, and configurational stabilization by interfacial water. Although the observation of outer-sphere arsenate surface complexes on a metal oxide surface is unprecedented, it is unclear if such species were absent in previous molecular-scale studies, as it is difficult for methods commonly used to investigate the mechanisms of arsenate adsorption to conclusively identify or rule out the presence of outer-sphere species when inner-sphere species are also present.     

A simple model for the effect of hydration on the distribution of ferrous iron at reduced hematite (0 1 2) surfaces

Using a simple ionic model with polarizable oxygen ions and dissociating water molecules, we have calculated the energetics governing the distribution of Fe(II)/Fe(III) ions at the reduced (2 × 1) surface of α-Fe₂O₃ (hematite) (0 1 2) under dry and hydrated conditions. The results show that systems with Fe(II) ions located in the near-surface region have lower potential energy for both dry and hydrated surfaces. The distribution is governed by coupling of the ferrous iron centers to positive charge associated with missing oxygen atoms on the dry reduced (2 × 1) (0 2 1) surface. As the surface is hydroxylated, the missing oxygen rows are filled and protons from dissociated water molecules become the positive charge centers, which couple more weakly to the ferrous iron centers. At the same time, the first-layer iron centers change from fourfold or fivefold coordination to sixfold coordination lowering the potential energy of ferric iron in the first layer and favoring migration of ferrous iron from the immediate surface sites. This effect can also be understood as reflecting stronger solvation of Fe(III) by the adsorbed water molecules and by hydrolysis reactions favoring Fe(III) ions at the immediate surface. The balance between these two driving forces, which changes as a function of hydration, provides a compelling explanation for the anomalous coverage dependence of water desorption in ultra-high vacuum experiments.     

Cation sorption on the muscovite (0 0 1) surface in chloride solutions using high-resolution X-ray reflectivity

The structure and mechanism of cation sorption at the (0 0 1) muscovite-water interface were investigated in 0.01 and 0.5 m KCl, CsCl, and CaCl₂ and 0.01 m BaCl₂ solutions at slightly acidic pH by high-resolution X-ray reflectivity. Structural relaxations of atom positions in the 2M₁ muscovite were small (?0.07 Å) and occurred over a distance of 30 to 40 Å perpendicular to the interface. Cations in all solutions were sorbed dominantly in the first and second solution layers adjacent to the mineral surface. The derived heights of the first solution layer in KCl and CsCl solutions, 1.67(6)-1.77(7) and 2.15(9)-2.16(2) Å, respectively, differ in magnitude by the approximate difference in crystallographic radii between K and Cs, and correspond closely to the interlayer cation positions in bulk K- and Cs-mica structures. The first solution layer heights in CaCl₂ and BaCl₂ solutions, 2.46(5)-2.56(11) and 2.02(5) Å, respectively, differ in a sense opposite to that expected based on crystallographic or hydrated radii of the divalent cations. The derived ion heights in all solutions imply that there is no intercalated water layer between the first solution layer and the muscovite surface. Molecular compositions were assigned to the first two solution layers in the electron density profiles using models that constrain the number density of sorbed cations, water molecules, and anions by considering the permanent negative charge of the muscovite and average solution density. The models result in partial charge balance (at least 50%) by cations sorbed in the first two layers in the 0.01 m solutions and approximately full charge balance in the 0.5 m solutions. Damped oscillations of model water density away from the first two solution layers agree with previous X-ray reflectivity results on the muscovite (0 0 1) surface in pure water.     

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)₃]⁶⁺.     

Solubility products of amorphous ferric arsenate and crystalline scorodite (FeAsO4 · 2H2O) and their application to arsenic behavior in buried mine tailings

Published solubility data for amorphous ferric arsenate and scorodite have been reevaluated using the geochemical code PHREEQC with a modified thermodynamic database for the arsenic species. Solubility product calculations have emphasized measurements obtained under conditions of congruent dissolution of ferric arsenate (pH < 3), and have taken into account ion activity coefficients, and ferric hydroxide, ferric sulfate, and ferric arsenate complexes which have association constants of 10^4.04 (FeH₂AsO₄²⁺), 10^9.86 (FeHAsO₄⁺), and 10^18.9 (FeAsO₄). Derived solubility products of amorphous ferric arsenate and crystalline scorodite (as log K_sp) are −23.0 ± 0.3 and −25.83 ± 0.07, respectively, at 25 °C and 1 bar pressure. In an application of the solubility results, acid raffinate solutions (molar Fe/As = 3.6) from the JEB uranium mill at McClean Lake in northern Saskatchewan were neutralized with lime to pH 2-8. Poorly crystalline scorodite precipitated below pH 3, removing perhaps 98% of the As(V) from solution, with ferric oxyhydroxide (FO) phases precipitated starting between pH 2 and 3. Between pH 2.18 and 7.37, the apparent log K_sp of ferric arsenate decreased from −22.80 to −24.67, while that of FO (as Fe(OH)₃) increased from −39.49 to −33.5. Adsorption of As(V) by FO can also explain the decrease in the small amounts of As(V)(aq) that remain in solution above pH 2-3. The same general As(V) behavior is observed in the pore waters of neutralized tailings buried for 5 yr at depths of up to 32 m in the JEB tailings management facility (TMF), where arsenic in the pore water decreases to 1-2 mg/L with increasing age and depth. In the TMF, average apparent log K_sp values for ferric arsenate and ferric hydroxide are −25.74 ± 0.88 and −37.03 ± 0.58, respectively. In the laboratory tests and in the TMF, the increasing crystallinity of scorodite and the amorphous character of the coexisting FO phase increases the stability field of scorodite relative to that of the FO to near-neutral pH values. The kinetic inability of amorphous FO to crystallize probably results from the presence of high concentrations of sulfate and arsenate.     

Uranyl adsorption and surface speciation at the imogolite-water interface: Self-consistent spectroscopic and surface complexation models

Macro- and molecular-scale knowledge of uranyl (U(VI)) partitioning reactions with soil/sediment mineral components is important in predicting U(VI) transport processes in the vadose zone and aquifers. In this study, U(VI) reactivity and surface speciation on a poorly crystalline aluminosilicate mineral, synthetic imogolite, were investigated using batch adsorption experiments, X-ray absorption spectroscopy (XAS), and surface complexation modeling. U(VI) uptake on imogolite surfaces was greatest at pH ∼7-8 (I = 0.1 M NaNO₃ solution, suspension density = 0.4 g/L [U(VI)]_i = 0.01-30 μM, equilibration with air). Uranyl uptake decreased with increasing sodium nitrate concentration in the range from 0.02 to 0.5 M. XAS analyses show that two U(VI) inner-sphere (bidentate mononuclear coordination on outer-wall aluminol groups) and one outer-sphere surface species are present on the imogolite surface, and the distribution of the surface species is pH dependent. At pH 8.8, bis-carbonato inner-sphere and tris-carbonato outer-sphere surface species are present. At pH 7, bis- and non-carbonato inner-sphere surface species co-exist, and the fraction of bis-carbonato species increases slightly with increasing I (0.1-0.5 M). At pH 5.3, U(VI) non-carbonato bidentate mononuclear surface species predominate (69%). A triple layer surface complexation model was developed with surface species that are consistent with the XAS analyses and macroscopic adsorption data. The proton stoichiometry of surface reactions was determined from both the pH dependence of U(VI) adsorption data in pH regions of surface species predominance and from bond-valence calculations. The bis-carbonato species required a distribution of surface charge between the surface and β charge planes in order to be consistent with both the spectroscopic and macroscopic adsorption data. This research indicates that U(VI)-carbonato ternary species on poorly crystalline aluminosilicate mineral surfaces may be important in controlling U(VI) mobility in low-temperature geochemical environments over a wide pH range (∼5-9), even at the partial pressure of carbon dioxide of ambient air (p_CO2 = 10^−3.45 atm).     

Surface simulation studies of the hydration of white rust Fe(OH)2, goethite α-FeO(OH) and hematite α-Fe2O3

Computer modelling techniques were used to elucidate the hydration behaviour of three iron (hydr)oxide minerals at the atomic level: white rust, goethite and hematite. A potential model was first adapted and tested against the bulk structures and properties of eight different iron oxides, oxyhydroxides and hydroxides, followed by surface simulations of Fe(OH)₂, α-FeO(OH) and α-Fe₂O₃. The major interaction between the adsorbing water molecules and the surface is through interaction of their oxygen ions with surface iron ions, followed by hydrogen-bonding to surface oxygen ions. The energies released upon the associative adsorption of water range from 1 to 17 kJ mol⁻¹ for Fe(OH)₂, 26 to 80 kJ mol⁻¹ for goethite and 40 to 85 kJ mol⁻¹ for hematite, reflecting the increasing oxidation of the iron mineral. Dissociative adsorption at goethite and hematite surfaces releases larger hydration energies, ranging from 120 to 208 kJ mol⁻¹ for goethite and 76 to 190 kJ mol⁻¹ for hematite.The thermodynamic morphologies of the minerals, based on the calculated surface energies, agree well with experimental morphologies, where these are available. When the partial pressures required for adsorption of water from the gas phase are plotted against temperature for the goethite and hematite surfaces, taking into account experimental entropies for water, it appears that these minerals may well be instrumental in the retention of water during the cyclic variations in the atmosphere of Mars.     

Structure of hydrated Zn at the rutile TiO2 (110)-aqueous solution interface: Comparison of X-ray standing wave, X-ray absorption spectroscopy, and density functional theory results

Adsorption of Zn²⁺ at the rutile TiO₂ (110)-aqueous interface was studied with Bragg-reflection X-ray standing waves (XSW), polarization-dependent surface extended X-ray absorption fine structure (EXAFS) spectroscopy, and density functional theory (DFT) calculations to understand the interrelated issues of adsorption site, its occupancy, ion-oxygen coordination and hydrolysis. At pH 8, Zn²⁺ was found to adsorb as an inner-sphere complex at two different sites, i.e., monodentate above the bridging O site and bidentate between two neighboring terminal O sites. EXAFS results directly revealed a four or fivefold first shell coordination environment for adsorbed Zn²⁺ instead of the sixfold coordination found for aqueous species at this pH. DFT calculations confirmed the energetic stability of a lower coordination environment for the adsorbed species and revealed that the change to this coordination environment is correlated with the hydrolysis of adsorbed Zn²⁺. In addition, the derived adsorption locations and the occupancy factors of both sites from three methods agree well, with some quantitative discrepancies in the minor site location among the XSW, EXAFS, and DFT methods. Additional XSW measurements showed that the adsorption sites of Zn²⁺ were unchanged at pH 6. However, the Zn²⁺ partitioning between the two sites changed substantially, with an almost equal distribution between the two types of sites at pH 6 compared to predominantly monodentate occupation at pH 8.