Molecular models of a hydrated calcite mineral surface

Hydrated mineral surfaces play an important role in many processes in biological, geological, and industrial applications. An energy force field was developed for molecular mechanics and molecular dynamics simulations of hydrated carbonate minerals and was applied to investigate the behavior of water on the calcite surface. The force field is a significant development for large-scale molecular simulations of these systems, and provides good agreement with experimental and previous modeling results. Simulations indicate that water molecules are significantly ordered near the calcite surface. The predominant surface configuration (75-80%) results from coordination of a water molecule with a single calcium cation-carbonate anion pair, while the less common situation involves water coordination with two ion pairs. Surface restructuring and variation in coordination in the water layers results in distinct distances for water oxygens above the calcite surface—a two-component first monolayer (2.3 and 3.0 Å) and a secondary monolayer (5.0 Å). The different coordinations also alter lateral displacement, hydrogen bonding, and surface-normal orientation of the water molecules. The ordering of water molecules and the formation of a unique hydrogen bonding network at the calcite surface influence the physical properties of the interfacial water. Surface exchange of water molecules is observed by molecular dynamics simulation to occur at a rate of one exchange per 10 ps. Diffusion coefficients derived from mean square displacement analysis of atomic trajectories indicate a dependence of water transport based on the distance of the water molecules from the calcite surface.     

An atomic force microscopy study of calcite dissolution in saline solutions: The role of magnesium ions

In situ Atomic Force Microscopy, AFM, experiments have been carried out using calcite cleavage surfaces in contact with solutions of MgSO₄, MgCl₂, Na₂SO₄ and NaCl in order to attempt to understand the role of Mg²⁺ during calcite dissolution. Although previous work has indicated that magnesium inhibits calcite dissolution, quantitative AFM analyses show that despite the fact that Mg²⁺ inhibits etch pit spreading, it increases the density and depth of etch pits nucleated on calcite surfaces and, subsequently, the overall dissolution rates: i.e., from 10^−11.75 mol cm⁻² s⁻¹ (in deionized water) up to 10^−10.54 mol cm⁻² s⁻¹ (in 2.8 M MgSO₄). Such an effect is concentration-dependent and it is most evident in concentrated solutions ([Mg²⁺] >> 50 mM). These results show that common soluble salts (especially Mg sulfates) may play a critical role in the chemical weathering of carbonate rocks in nature as well as in the decay of carbonate stone in buildings and statuary.     

Climatic dependence of stable carbon and oxygen isotope signals recorded in speleothems: From soil water to speleothem calcite

Understanding the relationship between stable isotope signals recorded in speleothems (δ¹³C and δ¹⁸O) and the isotopic composition of the carbonate species in the soil water is of great importance for their interpretation in terms of past climate variability. Here the evolution of the carbon isotope composition of soil water on its way down to the cave during dissolution of limestone is studied for both closed and open-closed conditions with respect to CO₂.The water entering the cave flows as a thin film towards the drip site. CO₂ degasses from this film within approx. 10 s by molecular diffusion. Subsequently, chemical and isotopic equilibrium is established on a time scale of several 10-100 s. The δ¹³C value of the drip water is mainly determined by the isotopic composition of soil CO₂. The evolution of the δ¹⁸O value of the carbonate species is determined by the long exchange time T_ex, between oxygen in carbonate and water of several 10,000 s. Even if the oxygen of the CO₂ in soil water is in isotopic equilibrium with that of the water, dissolution of limestone delivers oxygen with a different isotopic composition changing the δ¹⁸O value of the carbonate species. Consequently, the δ¹⁸O value of the rainwater will only be reflected in the drip water if it has stayed in the rock for a sufficiently long time.After the water has entered the cave, the carbon and oxygen isotope composition of the drip water may be altered by CO₂-exchange with the cave air. Exchange times, , of about 3000 s are derived. Thus, only drip water, which drips in less than 3000 s onto the stalagmite surface, is suitable to imprint climatic signals into speleothem calcite deposited from it.Precipitation of calcite proceeds with time constants, τ_p, of several 100 s. Different rate constants and equilibrium concentrations for the heavy and light isotopes, respectively, result in isotope fractionation during calcite precipitation. Since T_ex ? τ_p, exchange with the oxygen in the water can be neglected, and the isotopic evolution of carbon and oxygen proceed analogously. For drip intervals T_d < 0.1τ_p the isotopic compositions of both carbon and oxygen in the solution evolve linearly in time. The calcite precipitated at the apex of the stalagmite reflects the isotopic signal of the drip water.For long drip intervals, when calcite is deposited from a stagnant water film, long drip intervals may have a significant effect on the isotopic composition of the DIC. In this case, the isotopic composition of the calcite deposited at the apex must be determined by averaging over the drip interval. Such processes must be considered when speleothems are used as proxies of past climate variability.     

The role of background electrolytes on the kinetics and mechanism of calcite dissolution

The influence of background electrolytes on the mechanism and kinetics of calcite dissolution was investigated using in situ Atomic Force Microscopy (AFM). Experiments were carried out far from equilibrium by passing alkali halide salt (NaCl, NaF, NaI, KCl and LiCl) solutions over calcite cleavage surfaces. This AFM study shows that all the electrolytes tested enhance the calcite dissolution rate. The effect and its magnitude is determined by the nature and concentration of the electrolyte solution. Changes in morphology of dissolution etch pits and dissolution rates are interpreted in terms of modification in water structure dynamics (i.e. in the activation energy barrier of breaking water-water interactions), as well as solute and surface hydration induced by the presence of different ions in solution. At low ionic strength, stabilization of water hydration shells of calcium ions by non-paired electrolytes leads to a reduction in the calcite dissolution rate compared to pure water. At high ionic strength, salts with a common anion yield similar dissolution rates, increasing in the order Cl⁻ < I⁻ < F⁻ for salts with a common cation due to an increasing mobility of water around the calcium ion. Changes in etch pit morphology observed in the presence of F⁻ and Li⁺ are explained by stabilization of etch pit edges bonded by like-charged ions and ion incorporation, respectively. As previously reported and confirmed here for the case of F⁻, highly hydrated ions increased the etch pit nucleation density on calcite surfaces compared to pure water. This may be related to a reduction in the energy barrier for etch pit nucleation due to disruption of the surface hydration layer.     

Arsenate uptake by calcite: Macroscopic and spectroscopic characterization of adsorption and incorporation mechanisms

Batch uptake experiments and X-ray element mapping and spectroscopic techniques were used to investigate As(V) (arsenate) uptake mechanisms by calcite, including adsorption and coprecipitation. Batch sorption experiments in calcite-equilibrated suspensions (pH 8.3; PCO₂ = 10^−3.5 atm) reveal rapid initial sorption to calcite, with sorption rate gradually decreasing with time as available sorption sites decrease. An As(V)-calcite sorption isotherm determined after 24 h equilibration exhibits Langmuir-like behavior up to As concentrations of 300 μM. Maximum distribution coefficient values (K_d), derived from a best fit to a Langmuir model, are ∼190 L kg⁻¹.Calcite single crystals grown in the presence of As(V) show well-developed rhombohedral morphology with characteristic growth hillocks on surfaces at low As(V) concentrations (?5 μM), but habit modification is evident at As(V) concentrations ?30 μM in the form of macrostep development preferentially on the − vicinal surfaces of growth hillocks. Micro-X-ray fluorescence element mapping of surfaces shows preferential incorporation of As in the − vicinal faces relative to + vicinals. EXAFS fit results for both adsorption and coprecipitation samples confirm that As occurs in the 5+ oxidation state in tetrahedral coordination with oxygen, i.e., as arsenate. For adsorption samples, As(V) forms inner-sphere surface complexes via corner-sharing with Ca octahedra. As(V) coprecipitated with calcite substitutes in carbonate sites but with As off-centered, as indicated by two Ca shells, and with likely disruption of local structure. The results indicate that As(V) interacts strongly with the calcite surface, similar to often-cited analog phosphate, and uptake can occur via both adsorption and coprecipitation reactions. Therefore, calcite may be effective for partial removal of dissolved arsenate from aquatic and soil systems.     

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.     

Calcite dissolution kinetics in the system CaCO3-H2O-CO2 at high undersaturation

Dissolution rates of limestone covered by a water film open to a CO₂-containing atmosphere are controlled by the chemical composition of the CaCO₃-H₂O-CO₂ solution at the water-mineral interface. This composition is determined by the Ca²⁺-concentration at this boundary, conversion of CO₂ into H⁺ and in the solution, and by diffusional mass transport of the dissolved species from and towards the water-limestone interface. A system of coupled diffusion-reaction equations for Ca²⁺, , and CO₂ is derived. The Ca²⁺ flux rates at the surface of the mineral are defined by the PWP-empirical rate law. These flux rates by the rules of stoichiometry must be equal to the flux rates of CO₂ across the air-water interface. In the solution, CO₂ is converted into H⁺ and . At low water-film thickness this reaction becomes rate limiting. The time dependent diffusion-reaction equations are solved for free drift dissolution by a finite-difference scheme, to obtain the dissolution rate of calcite as a function of the average calcium concentration in the water film. Dissolution rates are obtained for high undersaturation. The results reveal two regimes of linear dissolution kinetics, which can be described by a rate law F = α_i(m_ic_eq − c), where c is the calcium concentration in the water film, c_eq the equilibrium concentration with respect to calcite. For index i = 0, a fast rate law, which here is reported for the first time, is found with α₀ = 3 × 10⁻⁶ m s⁻¹ and m₀ = 0.3. For c > m₀c_eq, a slow rate law is valid with α₁ = 3 × 10⁻⁷ m  s⁻¹ and m₁ = 1, which confirms earlier work. The numbers given above are valid for film thickness of several tenths of a millimetre and at 20 °C. These rates are proven experimentally, using a flat inclined limestone plate covered by a laminar flowing water film injected at an input point with known flow rate Q and calcium concentration. From the concentration measured after flow distance x the dissolution rates are determined. These experiments have been performed at a carbon-dioxide pressure of 0.00035 atm and also of 0.01 atm. The results are in good agreement to the theoretical predictions.     

Activity-composition relations in the system CaCO3-MgCO3 predicted from static structure energy calculations and Monte Carlo simulations

Thermodynamic mixing properties and subsolidus phase relations of the rhombohedral carbonate system, (1 − x) · CaCO₃ − x · MgCO₃, were modelled in the temperature range of 623-2023 K with static structure energy calculations based on well-parameterised empirical interatomic potentials. Relaxed static structure energies of a large set of randomly varied structures in a 4 × 4 × 1 supercell of calcite (a = 19.952 Å, c = 17.061 Å) were calculated with the General Utility Lattice Program (GULP). These energies were cluster expanded in a basis set of 12 pair-wise effective interactions. Temperature-dependent enthalpies of mixing were calculated by the Monte Carlo method. Free energies of mixing were obtained by thermodynamic integration of the Monte Carlo results. The calculated phase diagram is in good agreement with experimental phase boundaries.     

Specific effects of background electrolytes on the kinetics of step propagation during calcite growth

The mechanisms by which background electrolytes modify the kinetics of non-equivalent step propagation during calcite growth were investigated using Atomic Force Microscopy (AFM), at constant driving force and solution stoichiometry. Our results suggest that the acute step spreading rate is controlled by kink-site nucleation and, ultimately, by the dehydration of surface sites, while the velocity of obtuse step advancement is mainly determined by hydration of calcium ions in solution. According to our results, kink nucleation at acute steps could be promoted by carbonate-assisted calcium attachment. The different sensitivity of obtuse and acute step propagation kinetics to cation and surface hydration could be the origin of the reversed geometries of calcite growth hillocks (i.e., rate of obtuse step spreading < rate of acute step spreading) observed in concentrated (ionic strength, IS = 0.1) KCl and CsCl solutions. At low IS (0.02), ion-specific effects seem to be mainly associated with changes in the solvation environment of calcium ions in solution. With increasing electrolyte concentration, the stabilization of surface water by weakly paired salts appears to become increasingly important in determining step spreading rate. At high ionic strength (IS = 0.1), overall calcite growth rates increased with increasing hydration of calcium in solution (i.e., decreasing ion pairing of background electrolytes for sodium-bearing salts) and with decreasing hydration of the carbonate surface site (i.e., increasing ion pairing for chloride-bearing salts). Changes in growth hillock morphology were observed in the presence of Li⁺, F⁻ and , and can be interpreted as the result of the stabilization of polar surfaces due to increased ion hydration. These results increase our ability to predict crystal reactivity in natural fluids which contain significant amounts of solutes.     

Sorption and catalytic oxidation of Fe(II) at the surface of calcite

The effect of sorption and coprecipitation of Fe(II) with calcite on the kinetics of Fe(II) oxidation was investigated. The interaction of Fe(II) with calcite was studied experimentally in the absence and presence of oxygen. The sorption of Fe(II) on calcite occurred in two distinguishable steps: (a) a rapid adsorption step (seconds-minutes) was followed by (b) a slower incorporation (hours-weeks). The incorporated Fe(II) could not be remobilized by a strong complexing agent (phenanthroline or ferrozine) but the dissolution of the outmost calcite layers with carbonic acid allowed its recovery. Based on results of the latter dissolution experiments, a stoichiometry of 0.4 mol% Fe:Ca and a mixed carbonate layer thickness of 25 nm (after 168 h equilibration) were estimated. Fe(II) sorption on calcite could be successfully described by a surface adsorption and precipitation model (Comans & Middelburg, GCA51 (1987), 2587) and surface complexation modeling (Van Cappellen et al., GCA57 (1993), 3505; Pokrovsky et al., Langmuir16 (2000), 2677). The surface complex model required the consideration of two adsorbed Fe(II) surface species, >CO₃Fe⁺ and >CO₃FeCO₃H⁰. For the formation of the latter species, a stability constant is being suggested. The oxidation kinetics of Fe(II) in the presence of calcite depended on the equilibration time of aqueous Fe(II) with the mineral prior to the introduction of oxygen. If pre-equilibrated for >15 h, the oxidation kinetics was comparable to a calcite-free system (t_1/2 = 145 ± 15 min). Conversely, if Fe(II) was added to an aerated calcite suspension, the rate of oxidation was higher than in the absence of calcite (t_1/2 = 41 ± 1 min and t_1/2 = 100 ± 15 min, respectively). This catalysis was due to the greater reactivity of the adsorbed Fe(II) species, >CO₃FeCO₃H⁰, for which the species specific rate constant was estimated.     

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.     

The energetics of nanophase calcite

Calcium carbonate (CaCO₃) is an important component of the near-surface environment. Understanding the nature of its precipitation is important for a variety of environmental processes, as well as for the geologic sequestration of anthropogenic carbon dioxide. Calcite is the most thermodynamically stable bulk polymorph, but energy crossovers may exist that could favor the precipitation of vaterite or aragonite with decreasing particle size. The purpose of this study is to determine the surface energy of calcite, which is the first step towards understanding the effect of particle size on thermodynamic stability in the calcium carbonate system. The enthalpies of five well-characterized calcite samples (four nanophase and one bulk) were measured by acid solution isothermal and water adsorption calorimetric techniques. From the calorimetric data, the surface energies of calcite were determined to be 1.48 ± 0.21 and 1.87 ± 0.16 J/m² for hydrous and anhydrous surfaces. These values are similar to those measured for many oxides but larger than predicted from computational models for idealized calcite surfaces. The surfaces of synthetic CaCO₃ particles contain a range of planes and defect structures, which may give rise to the difference between the experimental and modeled values.     

Coprecipitation of chromate with calcite: Batch experiments and X-ray absorption spectroscopy

Batch experiments, combined with in situ spectroscopic methods, are used to examine the coprecipitation of Cr(VI) with calcite, including partitioning behavior, site-specific distribution of Cr on the surface of calcite single crystals, and local coordination of Cr(VI) in the calcite structure. It is found that the concentration of Cr incorporated in calcite increases with increasing Cr concentration in solution. The calculated apparent partition coefficient, , is highest at low Cr solution concentration, and decreases to a constant value with increasing Cr solution concentration. DIC images of the surface of calcite single crystals grown in the presence of exhibit well-defined growth hillocks composed of two pairs of symmetrically nonequivalent vicinal faces, denoted as + and −, which reflect the orientation of structurally nonequivalent growth steps. Micro-XRF mapping of the Cr distribution over a growth hillock shows preferential incorporation of Cr into the—steps, which is considered to result from differences in surface structure geometry. XANES spectra confirm that incorporated Cr is hexavalent, and no reduction of Cr(VI) in the X-ray beam was observed up to 2 days at room temperature. EXAFS fit results show the incorporated Cr(VI) has the expected first shell of 4 O at ∼1.64 ± 0.01 Å, consistent with . Best fit results show that the second shell is split with ∼2.5 Ca at ∼3.33 ± 0.05 and ∼2.2 Ca at ∼3.55 ± 0.05 Å, which confirms the incorporation of chromate into calcite. Consideration of possible local coordination indicates that significant distortion or disruption is required to accommodate in the calcite structure.     

Dissolution kinetics of calcite at 50-70 °C: An atomic force microscopic study under near-equilibrium conditions

Direct measurements of calcite faces were performed using in situ atomic force microscopy (AFM) to reveal the dissolution processes as a function of solution saturation state and temperature. Time-sequential AFM images demonstrated that step velocities at constant temperature increased with increasing undersaturation. The anisotropy of obtuse and acute step velocities appeared to become more significant as solutions approached equilibrium and temperature increased. At saturation state Ω > 0.02, a curvilinear boundary was formed at the intersection of two acute steps and the initially rhombohedral etch pit exhibited a nearly triangular shape. This suggests that the and steps may not belong to the calcite-aqueous solution equilibrium system. Further increase in the saturation state (Ω ? 0.3) led to a lack of etch pit formation and dissolution primarily occurred at existing steps, in accordance with Teng (2004). Analysis of step kinetics at different temperatures yielded activation energies of 25 ± 6 kJ/mol and 14 ± 13 kJ/mol for obtuse and acute steps, respectively. The inconsistencies in etch pit morphology, step anisotropy, and step activation energies from the present study with those of studies far-from-equilibrium can be explained by increased influence of the backward reaction, or growth, near-equilibrium. We propose that the backward reaction occurs preferentially at the acute-acute kink sites. The kinetics and effective activation energies of near-equilibrium calcite dissolution presented in this work provide accurate experimental data under likely CO₂ sequestration conditions, and thus are crucial to the development of robust geochemical models that predict the long-term performance of mineral-trapped CO₂.     

Sorption of phosphate onto calcite; results from batch experiments and surface complexation modeling

The adsorption of phosphate onto calcite was studied in a series of batch experiments. To avoid the precipitation of phosphate-containing minerals the experiments were conducted using a short reaction time (3 h) and low concentrations of phosphate (?50 μM). Sorption of phosphate on calcite was studied in 11 different calcite-equilibrated solutions that varied in pH, P_CO2, ionic strength and activity of Ca²⁺, and . Our results show strong sorption of phosphate onto calcite. The kinetics of phosphate sorption onto calcite are fast; adsorption is complete within 2-3 h while desorption is complete in less than 0.5 h. The reversibility of the sorption process indicates that phosphate is not incorporated into the calcite crystal lattice under our experimental conditions. Precipitation of phosphate-containing phases does not seem to take place in systems with ?50 μM total phosphate, in spite of a high degree of super-saturation with respect to hydroxyapatite (SI_HAP ? 7.83). The amount of phosphate adsorbed varied with the solution composition, in particular, adsorption increases as the activity decreases (at constant pH) and as pH increases (at constant activity). The primary effect of ionic strength on phosphate sorption onto calcite is its influence on the activity of the different aqueous phosphate species. The experimental results were modeled satisfactorily using the constant capacitance model with >CaPO₄Ca⁰ and either >CaHPO₄Ca⁺ or > as the adsorbed surface species. Generally the model captures the variation in phosphate adsorption onto calcite as a function of solution composition, though it was necessary to include two types of sorption sites (strong and weak) in the model to reproduce the convex shape of the sorption isotherms.     

The impact of salinity on the Mg/Ca and Sr/Ca ratio in the benthic foraminifera Ammonia tepida: Results from culture experiments

Over the last decade, sea surface temperature (SST) reconstructed from the Mg/Ca ratio of foraminiferal calcite has increasingly been used, in combination with the δ¹⁸O signal measured on the same material, to calculate the δ¹⁸O_w, a proxy for sea surface salinity (SSS). A number of studies, however, have shown that the Mg/Ca ratio is also sensitive to other parameters, such as pH or , and salinity. To increase the reliability of foraminiferal Mg/Ca ratios as temperature proxies, these effects should be quantified in isolation. Individuals of the benthic foraminifera Ammonia tepida were cultured at three different salinities (20, 33 and 40 psu) and two temperatures (10-15 °C). The Mg/Ca and Sr/Ca ratios of newly formed calcite were analyzed by Laser Ablation ICP-MS and demonstrate that the Mg concentration in A. tepida is overall relatively low (mean value per experimental condition between 0.5 and 1.3 mmol/mol) when compared to other foraminiferal species, Sr being similar to other foraminiferal species. The Mg and Sr incorporation are both enhanced with increasing temperatures. However, the temperature dependency for Sr disappears when the distribution factor D_Sr is plotted as a function of calcite saturation state (Ω). This suggests that a kinetic process related to Ω is responsible for the observed dependency of Sr incorporation on sea water temperature. The inferred relative increase in D_Mg per unit salinity is 2.8% at 10 °C and 3.3% at 15 °C, for the salinity interval 20-40 psu. This implies that a salinity increase of 2 psu results in enhanced Mg incorporation equivalent to 1 °C temperature increase. The D_Sr increase per unit salinity is 0.8% at 10 °C and 1.3% at 15 °C, for the salinity interval 20-40 psu.     

In situ AFM study of the dissolution and recrystallization behaviour of polished and stressed calcite surfaces

We investigated the dissolution behaviour of polished calcite surfaces in situ using a fluid-cell atomic force microscope. Polished calcite surfaces enabled us to study the effects of applied surface stress and crystallographic orientation on calcite dissolution pattern formation. Thin-sections of Iceland spar single-crystals polished either parallel or with a 5° miscut angle to cleavage planes were studied. Compressive surface stresses of up to 50 MPa were applied to some of the thin-section samples by means of elastic concave bending. Experiments were carried out in semi-stagnant deionized water under mainly transport limited dissolution conditions. Samples polished parallel to cleavage planes dissolved by the formation of etch-pits originating from polishing defects. The dissolution behaviour of 5° miscut surfaces was relatively unaffected by polishing defects, since no etch-pits developed in these samples. Dissolution of the miscut samples led to stepped or rippled surface patterns on the nanometer scale that coarsened during the first 30-40 min of the experiments. Possible reasons for the pattern-coarsening were: (i) progressive bunching of retreating dissolution steps and (ii) surface energy driven recrystallization (Ostwald ripening) under transport limited dissolution conditions. A flat polished miscut surface in calcite may recrystallize into a hill-and-valley structure in a (near-)saturated solution so as to lower its total surface free energy in spite of a larger surface area. No clear effect of applied stress on dissolution pattern formation has been observed.     

Experimental study of germanium adsorption on goethite and germanium coprecipitation with iron hydroxide: X-ray absorption fine structure and macroscopic characterization

Adsorption of germanium on goethite was studied at 25 °C in batch reactors as a function of pH (1-12), germanium concentration in solution (10⁻⁷ to 0.002 M) and solid/solution ratio (1.8-17 g/L). The maximal surface site density determined via Ge adsorption experiments at pH from 6 to 10 is equal to 2.5 ± 0.1 μmol/m². The percentage of adsorbed Ge increases with pH at pH < 9, reaches a maximum at pH ∼ 9 and slightly decreases when pH is further increased to 11. These results allowed generation of a 2-pK Surface Complexation Model (SCM) which implies a constant capacitance of the electric double layer and postulates the presence of two Ge complexes, and , at the goethite-solution interface. Coprecipitation of Ge with iron oxy(hydr)oxides formed during Fe(II) oxidation by atmospheric oxygen or by Fe(III) hydrolysis in neutral solutions led to high Ge incorporations in solid with maximal Ge/Fe molar ratio close to 0.5. The molar Ge/Fe ratio in precipitated solid is proportional to that in the initial solution according to the equation (Ge/Fe)_solid = k × (Ge/Fe)_solution with 0.7 ? k ? 1.0. The structure of adsorbed and coprecipitated Ge complexes was further characterized using XAFS spectroscopy. In agreement with previous data on oxyanions adsorption on goethite, bi-dentate bi-nuclear surface complexes composed of tetrahedrally coordinated Ge attached to the corners of two adjacent Fe octahedra represent the dominant contribution to the EXAFS signal. Coprecipitated samples with Ge/Fe molar ratios >0.1, and samples not aged in solution (<1 day) having intermediate Ge/Fe ratios (0.01-0.1) show 4 ± 0.3 oxygen atoms at 1.76 ± 0.01 Å around Ge. Samples less concentrated in Ge (0.001 < Ge/Fe < 0.10) and aged longer times in solution (up to 280 days) exhibit a splitting of the first atomic shell with Ge in both tetrahedral (R = 1.77 ± 0.02 Å) and octahedral (R = 1.92 ± 0.03 Å) coordination with oxygen. In these samples, octahedrally coordinated Ge accounts for up to ∼20% of the total Ge. For the least concentrated samples (Ge/Fe < 0.001-0.0001) containing lepidocrocite, 30-50% of total co-precipitated germanium substitutes for Fe in octahedral sites with the next-nearest environment dominated by edge-sharing GeO₆-FeO₆ linkages (R_Ge-Fe ∼ 3.06 Å). It follows from the results of our study that the largest structural change of Ge (from tetrahedral to octahedral environment) occurs during its coprecipitation with Fe hydroxide at Ge/Fe molar ratio ?0.0001. These conditions are likely to be met in many superficial aquatic environments at the contact of anoxic groundwaters with surficial oxygenated solutions. Adsorption and coprecipitation of Ge with solid Fe oxy(hydr)oxides and organo-mineral colloids and its consequence for Ge/Si fractionation and Ge geochemical cycle are discussed.     

Adsorption of Sr and Ba at the cleaved mica-water interface: Free energy profiles and interfacial structure

Monte Carlo simulations show that the adsorption position of the Sr²⁺ or Ba²⁺ ion on the cleaved muscovite surface is determined by the radius of the ion’s first hydration shell, hydrogen bonding of the first shell water molecules with the basal oxygens of muscovite as well as a requirement of minimization of the number of muscovite’s lattice cations in the ion’s first coordination shell. Accordingly, Sr²⁺ or Ba²⁺ adsorbs in ditrigonal cavities at a distance of 1.12 Å or 1.35 Å, respectively, from the basal surface on dehydrated mica and above tetrahedral substitutions at a height of 1.93 ± 0.02 Å or 2.15 ± 0.03 Å, respectively, at the highest simulated water coverage of 28 H₂O per ion. The ion’s displacement from a ditrigonal cavity occurs upon adsorption of 2 H₂O per ion for Sr²⁺ and 3 H₂O per ion for Ba²⁺. At a coverage of 28 H₂O per ion, outer-sphere adsorption of Sr²⁺ or Ba²⁺ at a height of 3.9 ± 0.2 Å or 4.17 ± 0.07 Å, respectively, is possible albeit unfavorable on the free energy scale by 107 ± 7 kJ/mol or 89 ± 13 kJ/mol, respectively, as compared to inner-sphere adsorption. Activation energies for the transformation between inner-sphere and outer-sphere adsorptions are calculated to be 121 ± 3 kJ/mol for Sr²⁺ and 99 ± 10 kJ/mol for Ba²⁺. A comparison of these values with those reported recently for Mg²⁺ and Ca²⁺ results in an adsorption affinity sequence Mg²⁺ > Ca²⁺ > Sr²⁺ > Ba²⁺ in agreement with the sequence predicted recently for low dielectric constant solids (which include mica) (Sverjensky, 2006). A recent resonant anomalous X-ray reflectivity study of Sr²⁺ adsorption on muscovite (Park et al., 2006) has questioned the common assumption (Stumm, 1992), which is supported by the present simulation results, that inner-sphere adsorption is stronger than outer-sphere adsorption. A modification of the cleaved muscovite surface as a result of Sr²⁺ adsorption in muscovite’s ditrigonal cavities and related destabilization of muscovite’s hydroxyl groups is proposed as a possible reason for this controversy.     

Effects of substrate structure and composition on the structure, dynamics, and energetics of water at mineral surfaces: A molecular dynamics modeling study

Molecular dynamics computer simulations of the molecular structure, diffusive dynamics and hydration energetics of water adsorbed on (0 0 1) surfaces of brucite Mg(OH)₂, gibbsite Al(OH)₃, hydrotalcite Mg₂Al(OH)₆Cl · 2H₂O, muscovite KAl₂(Si₃Al)O₁₀(OH)₂, and talc Mg₃Si₄O₁₀(OH)₂ provide new insight into the relationships between the substrate structure and composition and the molecular-scale structure and properties of the interfacial water. For the three hydroxide phases studied here, the differences in the structural charge on the octahedral sheet, cation occupancies and distributions, and the orientations of OH groups all affect the surface water structure. The density profiles of water molecules perpendicular to the surface are very similar, due to the prevalent importance of H-bonding between the surface and the water and to their similar layered crystal structures. However, the predominant orientations of the surface water molecules and the detailed two-dimensional near-surface structure are quite different. The atomic density profiles and other structural characteristics of water at the two sheet silicate surfaces are very different, because the talc (0 0 1) surface is hydrophobic whereas the muscovite (0 0 1) surface is hydrophilic. At the hydrophilic and electrostatically neutral brucite and gibbsite (0 0 1) surfaces, both donating and accepting H-bonds from the H₂O molecules are important for the development of a continuous hydrogen bonding network across the interfacial region. For the hydrophilic but charged hydrotalcite and muscovite (0 0 1) surfaces, only accepting or donating H-bonds from the water molecules contribute to the formation of the H-bonding network at the negatively and positively charged interfaces, respectively. For the hydrophobic talc (0 0 1) surface, H-bonds between water molecules and the surface sites are very weak, and the H-bonds among H₂O molecules dominate the interfacial H-bonding network. For all the systems studied, the orientation of the interfacial water molecules in the first few layers is influenced by both the substrate surface charge and the ability by the surfaces to facilitate H-bond formation. The first layer of water molecules at all surfaces is well ordered in the xy plane (parallel to the surface) and the atomic density distributions reflect the substrate crystal structure. The enhanced ordering of water molecules at the interfaces indicates reduced orientational and translational entropy. In thin films, water molecules are more mobile parallel to the surface than perpendicular to it due to spatial constraints. At neutral, hydrophilic substrates, single-monolayer surface coverage stabilizes the adsorbed water molecules and results in a minimum of the surface hydration energy. In contrast, at the charged and hydrophilic muscovite surface, the hydration energy increases monotonically with increasing water coverage over the range of coverages studied. At the neutral and hydrophobic talc surface, the adsorption of H₂O is unfavorable at all surface coverages, and the hydration energy decreases monotonically with increasing coverage.