First-Principles Prediction of Surface Wetting

We have developed a method for predicting the solvation contribution to solid-liquid interfacial tension (IFT) based on density functional theory and the implicit solvent model COSMO-RS. Our method can be used to predict wetting behavior for a solid surface in contact with two liquids. We benchmarked our method against measurements of contact angle from water-in-oil on silica wafers and a range of self-assembled monolayers (SAMs) with different compositions, ranging from oil-wet to water-wet. We also compared our predictions to literature data for wetting of a polydimethylsilane surface. By explicitly including deprotonation for silica surfaces and carboxylic acid SAMs, very good agreement was obtained with experimental data for nearly all surfaces. Poor agreement was found for amine-terminated SAMs, which could be the result of both method and model insufficiencies and impurities known to be present for such surfaces. Solid-liquid IFT cannot be measured directly, making predictions such as from our method all the more important.

Surface Reactivity and Dissolution Properties of Alumina-Silica Glasses and Fibers

The ability of bulk glass and fibers to react in aqueous solution, with organic polymers and coupling agents, depends on the surface charge, reactivity, and adsorption properties of the glass surface, i.e. the character and density of surface -OH groups, whereas glass and fiber chemical stability and biosolubility depend on the resistance to dissolution. If glass dissolution products are accumulated in a media, they can change the surface properties by specific adsorption. We determined the -OH surface concentration, reactivity, adsorption, and dissolution properties of aluminosilicate glasses containing various modifiers and compared the results with the behavior of complex mineral wool fibers. Using proton consumption and element release from batch surface titration experiments, over the range 5 < pH < 10, surface -OH adsorption properties were modeled with the FITEQL program. During titration, network modifiers in the glass subsurface are preferentially replaced by protons, resulting in cation accumulation in the solution and formation of a leached layer enriched with Si on the solid. The behavior of Al was different. At 5 < pH < 9, only very small amounts of Al were found in the leachates, which can be explained by almost complete Al adsorption as stable surface complexes, i.e. >XOAl(OH)₂ (where X = Si or Al and > represents the surface). At pH > 9, divalent cations adsorbed specifically, as >XOMe⁺ complexes (Me = Ca or Mg). This deeper understanding of the surface behavior of glasses and fibers is important for the design of composite materials, for applications in biology and medicine and in materials production in general, as well as for understanding natural processes, such as global uptake estimates of CO₂ during rock weathering.

Mechanistic insight into biopolymer induced iron oxide mineralization through quantification of molecular bonding

Microbial production of iron (oxyhydr)oxides on polysaccharide rich biopolymers occurs on such a vast scale that it impacts the global iron cycle and has been responsible for major biogeochemical events. Yet the physiochemical controls these biopolymers exert on iron (oxyhydr)oxide formation are poorly understood. Here we used dynamic force spectroscopy to directly probe binding between complex, model and natural microbial polysaccharides and common iron (oxyhydr)oxides. Applying nucleation theory to our results demonstrates that if there is a strong attractive interaction between biopolymers and iron (oxyhydr)oxides, the biopolymers decrease the nucleation barriers, thus promoting mineral nucleation. These results are also supported by nucleation studies and density functional theory. Spectroscopic and thermogravimetric data provide insight into the subsequent growth dynamics and show that the degree and strength of water association with the polymers can explain the influence on iron (oxyhydr)oxide transformation rates. Combined, our results provide a mechanistic basis for understanding how polymer-mineral-water interactions alter iron (oxyhydr)oxides nucleation and growth dynamics and pave the way for an improved understanding of the consequences of polymer induced mineralization in natural systems.

Thermodynamic and Kinetic Parameters for Calcite Nucleation on Peptoid and Model Scaffolds: A Step toward Nacre Mimicry

The production of novel composite materials, assembled using biomimetic polymers known as peptoids (N-substituted glycines) to nucleate CaCO₃, can open new pathways for advanced material design. However, a better understanding of the heterogeneous CaCO₃ nucleation process is a necessary first step. We determined the thermodynamic and kinetic parameters for calcite nucleation on self-assembled monolayers (SAMs) of nanosheet-forming peptoid polymers and simpler, alkanethiol analogues. We used nucleation rate studies to determine the net interfacial free energy (γ _net) for the peptoid-calcite interface and for SAMs terminated with carboxyl headgroups, amine headgroups, or a mix of the two. We compared the results with γ _net determined from dynamic force spectroscopy (DFS) and from density functional theory (DFT), using COSMO-RS simulations. Calcite nucleation has a lower thermodynamic barrier on the peptoid surface than on carboxyl and amine SAMs. From the relationship between nucleation rate (J ₀) and saturation state, we found that under low-saturation conditions, i.e. <3.3 (pH 9.0), nucleation on the peptoid substrate was faster than that on all of the model surfaces, indicating a thermodynamic drive toward heterogeneous nucleation. When they are taken together, our results indicate that nanosheet-forming peptoid monolayers can serve as an organic template for CaCO₃ polymorph growth.

Achieving arsenic concentrations of <1 μg/L by Fe(0) electrolysis: The exceptional performance of magnetite

Consumption of drinking water containing arsenic at concentrations even below the World Health Organization provisional limit of 10 μg/L can still lead to unacceptable health risks. Consequently, the drinking water sector in the Netherlands has recently agreed to target 1 μg/L of arsenic in treated water. Unfortunately, in many poor, arsenic-affected countries, the costs and complexity of current methods that can achieve <1 μg/L are prohibitive, which highlights the need for innovative methods that can remove arsenic to <1 μg/L without costly support infrastructure and complicated supply chains. In this work, we used Fe(0) electrolysis, a low cost and scalable technology that is also known as Fe(0) electrocoagulation (EC), to achieve <1 μg/L residual dissolved arsenic. We compared the arsenic removal performance of green rust (GR), ferric (oxyhydr)oxides (Fe(III) oxides) and magnetite (Mag) generated by EC at different pH (7.5 and 9) in the presence of As(III) or As(V) (initial concentrations of 200-11,000 μg/L). Although GR and Fe(III) oxides removed up to 99% of initial arsenic, neither Fe phase could reliably meet the 1 μg/L target at both pH values. In contrast, EC-generated Mag consistently achieved <1 μg/L, regardless of the initial As(V) concentration and pH. Only solutions with initial As(III) concentrations ≥2200 μg/L resulted in residual arsenic >1 μg/L. As K-edge X-ray absorption spectroscopy showed that Mag also sorbed arsenic in a unique mode, consistent with partial arsenic incorporation near the particle surface. This sorption mode contrasts with the binuclear, corner sharing surface complex for GR and Fe(III) oxides, which could explain the difference in arsenic removal efficiency among the three Fe phases. Our results suggest that EC-generated Mag is an attractive method for achieving <1 μg/L particularly in decentralized water treatment.

Wetting Properties of the CO2-Water-Calcite System via Molecular Simulations: Shape and Size Effects

Assessment of the risks and environmental impacts of carbon geosequestration requires knowledge about the wetting behavior of mineral surfaces in the presence of CO₂ and the pore fluids. In this context, the interfacial tension (IFT) between CO₂ and the aqueous fluid and the contact angle, θ, with the pore mineral surfaces are the two key parameters that control the capillary pressure in the pores of the candidate host rock. Knowledge of these two parameters and their dependence on the local conditions of pressure, temperature, and salinity is essential for the correct prediction of structural and residual trapping. We have performed classical molecular dynamics simulations to predict the CO₂-water IFT and the CO₂-water-calcite contact angle. The IFT results are consistent with previous simulations, where simple point charge water models have been shown to underestimate the water surface tension, thus affecting the simulated IFT values. When combined with the EPM2 CO₂ model, the SPC/Fw water model indeed underestimates the IFT in the low-pressure region at all temperatures studied. On the other hand, at high pressure and low temperature, the IFT is overestimated by ∼5 mN/m. Literature data regarding the CO₂/water/calcite contact angle on calcite are contradictory. Using our new set of force field parameters, we performed NVT simulations at 323 K and 20 MPa to calculate the contact angle of a water droplet on the calcite {10.4} surface in a CO₂ atmosphere. We performed simulations for both spherical and cylindrical droplet configurations for different initial radii to study the size dependence of the water contact angle on calcite in the presence of CO₂. Our results suggest that the contact angle of a cylindrical droplet, is independent of droplet size, for droplets with a radius of 50 Å or more. On the contrary, spherical droplets make a contact angle that is strongly influenced by their size. At the largest size explored in this study, both spherical and cylindrical droplets converge to the same contact angle, 38°, indicating that calcite is strongly wetted by water.

Effect of Divalent Cations on the Interaction of Carboxylate Self-Assembled Monolayers

Interactions between organic molecules in aqueous environments, whether in the fluid phase or adsorbed on solids, are often affected by the cations present in the solution. We investigated, at nanometer scale, how surface carboxylate interactions are influenced by dissolved divalent cations: Mg²⁺, Ca²⁺, Sr²⁺, and Ba²⁺. Self-assembled monolayer (SAM) surfaces with exposed terminations of alkyl, -CH₃, carboxylate, -COO^- , or dicarboxylate, -DiCOO^-, were deposited on gold-coated tips and substrates. We used atomic force microscopy (AFM), in chemical force mapping (CFM) mode, to measure adhesion forces between various combinations of SAMs on the tip and substrate, in solutions of 0.5 M NaCl, that contained 0.012 M of one of the divalent cations. The type of cation, the number of carboxyl groups that interact, and their structure on the SAM influenced adhesion between the surfaces. The effect of the reference solution, which only contains Na⁺ cations, on adhesion force was mainly attributed to van der Waals and hydrophobic forces, explaining the lower force in systems that are more hydrophilic, i.e., -COO^--COO^-, and higher force for more hydrophobic systems. For charged surfaces, i.e., -COO^- and -DiCOO^-, in divalent cation solutions results were consistent with ion bridging. The inclusion of a hydrophobic surface, i.e., the -CH₃-COO^- or -CH₃-DiCOO^- system, decreased the possibility for strong cation bridging with the charged surface, resulting in lower adhesion. For systems including -COO^-, the adhesion force series followed the inverse cation hydrated radius trend (Na⁺ ≈ Mg²⁺ < Sr²⁺ < Ca²⁺ < Ba²⁺) whereas -DiCOO^- was responsible for lower adhesion force and modified trends, depending on the corresponding surface in the system. Differences in force magnitude between the monolayers were correlated with lower charge availability on the -DiCOO^- surface as a result of fewer active sites, probably because of the tendency of exposed malonate surface groups to interact between them, as well as high rigidity, resulting from the molecule structure. The characteristic response of the -DiCOO^- surface in solutions of Sr²⁺ and Ca²⁺ was correlated with possible malonate complexation modes. Comparison with previous studies suggested that the strong response of a -DiCOO^- surface to Sr²⁺ resulted from bidentate chelation, whereas Ca²⁺ response was attributed to alpha-mode association to malonate.

Numerical simulations of NMR relaxation in chalk using local Robin boundary conditions

The interpretation of nuclear magnetic resonance (NMR) data is of interest in a number of fields. In Ögren (2014) local boundary conditions for random walk simulations of NMR relaxation in digital domains were presented. Here, we have applied those boundary conditions to large, three-dimensional (3D) porous media samples. We compared the random walk results with known solutions and then applied them to highly structured 3D domains, from images derived using synchrotron radiation CT scanning of North Sea chalk samples. As expected, there were systematic errors caused by digitalization of the pore surfaces so we quantified those errors, and by using linear local boundary conditions, we were able to significantly improve the output. We also present a technique for treating numerical data prior to input into the ESPRIT algorithm for retrieving Laplace components of time series from NMR data (commonly called T-inversion).

Specific Ion Effects on the Interaction of Hydrophobic and Hydrophilic Self-Assembled Monolayers

Interactions between mineral surfaces and organic molecules are fundamental to life processes. The presence of cations in natural environments can change the behavior of organic compounds and thus alter the mineral-organic interfaces. We investigated the influence of Na⁺, Mg²⁺, Ca²⁺, Sr²⁺, and Ba²⁺ on the interaction between two models, self-assembled monolayers, that were tailored to have hydrophobic -CH₃ or hydrophilic -COO(H) terminations. Atomic force microscopy in chemical force mapping mode, where the tips were functionalized with the same terminations, was used to measure adhesion forces between the tip and substrate surfaces, to gather fundamental information about the role of these cations in the behavior of organic compounds and the surfaces where they adsorb. Adhesion force between hydrophobic surfaces in 0.5 M NaCl solutions that contained 0.012 M divalent cations did not change, regardless of the ionic potential, that is, the charge per unit radius, of the cation. For systems where one or the other surface was functionalized with carboxylate, -COO(H), mostly in its deprotonated form, -COO^-, a reproducible change in the adhesion force was observed for each of the ions. The trend of increasing adhesion force followed the pattern: Na⁺ ≈ Mg²⁺ < Sr²⁺ < Ca²⁺ < Ba²⁺, suggesting that ionic potential, thus hydrated radius, controls the interaction. The presence of a -CH₃ surface in the asymmetric system leads to lower adhesion forces than in the hydrophilic system, whereas the ionic trend remains the same. Although specific ion effects are felt in both systems, the lower adhesion force in the asymmetric system, compared with the hydrophilic system, implies that the -CH₃ surface plays an important role.

Adsorption of organic ligands on low surface charge clay minerals: the composition in the aqueous interface region

An understanding of the mechanisms that control the adsorption of organic molecules on clay minerals is of interest in several branches of science and industry. Oil production using low salinity injection fluids can increase yields by as much as 40% over standard injection with seawater or formation water. The mechanism responsible for the low salinity response is still debated, but one hypothesis is a change in pore surface wettability. Organic contamination in soil and drinking water aquifers is a challenge for municipal water suppliers and for agriculture. A better understanding is needed for how mineral species, solution composition and pH affect the desorption of low molecular weight organic ligands from clay minerals and consequently their wettability. We used X-ray photoelectron spectroscopy under cryogenic conditions to investigate the in situ composition in the mineral-solution interface region in a series of experiments with a range of pH and ion concentrations. We demonstrate that both chlorite and kaolinite release organic molecules under conditions relevant for low salinity water flooding. This release increases with a higher solution pH but is only slightly affected by the character of the organic ligand. This is consistent with the observation that low salinity enhanced oil recovery correlates with the presence of chlorite and kaolinite. Our results indicate that the pore surface charge and salinity of formation water and injection fluids are key parameters in determining the low salinity response. In general, our results imply that clay mineral surface charge influences the composition in the interface through an affinity for organic molecules.

Composition in the Interface between Clay Mineral Surfaces and Divalent Cation Electrolytes

The interfacial free energy of a solid, which determines its adsorption properties, depends on interactions between the surface and the fluid. A change in surface composition can completely change the behavior of the solid. Decades of work have explored adsorption and its effects at solid-fluid interfaces from the macroscopic perspective and using molecular modeling, so the concept of the electric double layer (EDL) is well established in the community. However, direct, molecular level, experimental observations of the composition within the interface region, and its change with time and conditions, are not abundant. We used cryogenic X-ray photoelectron spectroscopy (cryoXPS) to observe the composition in the clay mineral-solution interface region as a function of bulk solution composition, on illite and chlorite in MgCl₂ and CaCl₂ electrolytes, over a range of concentrations (1-125 mM), in situ, on vitrified samples. These samples were prepared from very thin smears of centrifuged wet paste that were instantaneously chilled to liquid N₂ temperature. They preserved the adsorbed solution in its amorphous state, maintaining the location of the ions and water with respect to the solid, without the disruption that occurs during drying or the rearrangement that results as water crystallizes during freezing. With decreasing ionic strength, we could directly monitor the loss of negative charge in the interface region, producing an anion deficiency, as predicted by theory. The Cl^-/Me²⁺ ratio dropped below 1 for chlorite at 12-25 mM MeCl₂ and for illite at 75-100 mM. In addition to better understanding of clay mineral behavior in solution, this work demonstrates that only those clay minerals where surface charge density is the same or lower than that for chlorite contribute to a low salinity enhanced oil recovery response (LS EOR). This explains many of the contradictory results from studies about the role of clay minerals in LS EOR.

Controlling biomineralisation with cations

The production of polymers for controlling calcite growth is a well-known approach in biomineralising organisms. Numerous studies have shown that polymers significantly influenced the growth rate and morphology of CaCO₃ but little is known about how the polymers are actually controlled by the organisms. Here we show that cations control the effect of polysaccharides and that these processes have been in place for at least 60 million years. We studied the interaction between cleaved samples of pure calcite and ancient coccolith associated polysaccharides (aPS) that we had extracted from the samples of Cretaceous chalk, in solutions that contained one of the common seawater cations, K⁺, Ca²⁺, Mg²⁺ and Sr²⁺. With atomic and chemical force microscopy (AFM and CFM), we showed that K⁺, Ca²⁺ and Sr²⁺ complex aPS through a weak, outer sphere bonding, giving the aPS affinity to sites on steps and terraces. In contrast, Mg²⁺ enhanced the formation of stronger and longer aPS complexes, resulting in low affinity to calcite terraces and strong affinity to steps. It is known that adsorption is influenced by ionic potential and ionic strength. Our results show that cation-polysaccharide complexing can modify the effectiveness of the polymer. Thus, creating organic molecules with cation complexing ability is an effective strategy for regulating mineral growth, both now and in the past.

Hydrolytic Stability of 3-Aminopropylsilane Coupling Agent on Silica and Silicate Surfaces at Elevated Temperatures

3-Aminopropylsilane (APS) coupling agent is widely used in industrial, biomaterial, and medical applications to improve adhesion of polymers to inorganic materials. However, during exposure to elevated humidity and temperature, the deposited APS layers can decompose, leading to reduction in coupling efficiency, thus decreasing the product quality and the mechanical strength of the polymer-inorganic material interface. Therefore, a better understanding of the chemical state and stability of APS on inorganic surfaces is needed. In this work, we investigated APS adhesion on silica wafers and compared its properties with those on complex silicate surfaces such as those used by industry (mineral fibers and fiber melt wafers). The APS was deposited from aqueous and organic (toluene) solutions and studied with surface sensitive techniques, including X-ray photoelectron spectroscopy (XPS), atomic force microscopy (AFM), streaming potential, contact angle, and spectroscopic ellipsometry. APS configuration on a model silica surface at a range of coverages was simulated using density functional theory (DFT). We also studied the stability of adsorbed APS during aging at high humidity and elevated temperature. Our results demonstrated that APS layer formation depends on the choice of solvent and substrate used for deposition. On silica surfaces in toluene, APS formed unstable multilayers, while from aqueous solutions, thinner and more stable APS layers were produced. The chemical composition and substrate roughness influence the amount of deposited APS. More APS was deposited and its layers were more stable on fiber melt than on silica wafers. The changes in the amount of adsorbed APS can be successfully monitored by streaming potential. These results will aid in improving industrial- and laboratory-scale APS deposition methods and increasing adhesion and stability, thus increasing the quality and effectiveness of materials where APS is used as a coupling agent.

A Microkinetic Model of Calcite Step Growth

In spite of decades of research, mineral growth models based on ion attachment and detachment rates fail to predict behavior beyond a narrow range of conditions. Here we present a microkinetic model that accurately reproduces calcite growth over a very wide range of published experimental data for solution composition, saturation index, pH and impurities. We demonstrate that polynuclear complexes play a central role in mineral growth at high supersaturation and that a classical complexation model is sufficient to reproduce measured rates. Dehydration of the attaching species, not the mineral surface, is rate limiting. Density functional theory supports our conclusions. The model provides new insights into the molecular mechanisms of mineral growth that control biomineralization, mineral scaling and industrial material synthesis.

Modelling how incorporation of divalent cations affects calcite wettability-implications for biomineralisation and oil recovery

Using density functional theory and geochemical speciation modelling, we predicted how solid-fluid interfacial energy is changed, when divalent cations substitute into a calcite surface. The effect on wettability can be dramatic. Trace metal uptake can impact organic compound adsorption, with effects for example, on the ability of organisms to control crystal growth and our ability to predict the wettability of pore surfaces. Wettability influences how easily an organic phase can be removed from a surface, either organic compounds from contaminated soil or crude oil from a reservoir. In our simulations, transition metals substituted exothermically into calcite and more favourably into sites at the surface than in the bulk, meaning that surface properties are more strongly affected than results from bulk experiments imply. As a result of divalent cation substitution, calcite-fluid interfacial energy is significantly altered, enough to change macroscopic contact angle by tens of degrees. Substitution of Sr, Ba and Pb makes surfaces more hydrophobic. With substitution of Mg and the transition metals, calcite becomes more hydrophilic, weakening organic compound adsorption. For biomineralisation, this provides a switch for turning on and off the activity of organic crystal growth inhibitors, thereby controlling the shape of the associated mineral phase.

Microstructure and micromechanics of the heart urchin test from X-ray tomography

The microstructure of many echinoid species has long fascinated scientists because of its high porosity and outstanding mechanical properties. We have used X-ray microtomography to examine the test of Echinocardium cordatum (heart urchin), a burrowing cousin of the more commonly known sea urchins. Three dimensional imaging demonstrates that the bulk of the test is composed of only two distinct, highly porous, fenestrated regions (stereom), in which the thickness of the struts is constant. Different degrees of porosity are achieved by varying the spacing of the struts. Drawing an analogy to vertebrate trabecular bone, where for example, human bone has a connectivity density of ≈1/mm(3), we measure up to 150,000 strut connections per mm(3). Simulations of mechanical loading using finite element calculations indicate that the test performs at very close to the optimum expected for foams, highlighting the functional link between structure and mechanical properties.

Adsorption of ethanol and water on calcite: dependence on surface geometry and effect on surface behavior

Molecular dynamics (MD) simulations were used to explore adsorption on calcite, from a 1:1 mixture of ethanol and water, on planar {10.4} and stepped, i.e. vicinal, surfaces. Varying the surface geometry resulted in different adsorption patterns, which would directly influence the ability of ethanol to control calcite crystal growth, dissolution, and adsorption/desorption of other ions and molecules. Ethanol forms a well-ordered adsorbed layer on planar faces and on larger terraces, such as between steps and defects, providing little chance for water, with its weaker attachment, to displace it. However, on surfaces with steps, adsorption affinity depends on the length of the terraces between the steps. Long terraces allow ethanol to form a well-ordered, hydrophobic layer, but when step density is high, ethanol adsorption is less ordered, allowing water to associate at and near the steps and even displacing pre-existing ethanol. Water adsorbed at steps forms mass transport pathways between the bulk solution and the solid surface. Our simulations confirm the growth inhibiting properties of ethanol, also explaining how certain crystal faces are more stabilized because of their surface geometry. The -O(H) functional group on ethanol forms tight bonds with calcite; the nonpolar, -CH3 ends, which point away from the surface, create a hydrophobic layer that changes surface charge, thus wettability, and partly protects calcite from precipitation and dissolution. These tricks could easily be adopted by biomineralizing organisms, allowing them to turn on and off crystal growth. They undoubtedly also play a role in the wetting properties of mineral surfaces in commercial CaCO3 manufacture, oil production, and contamination remediation.

Adhesion of alkane as a functional group on muscovite and quartz: dependence on pH and contact time

The interactions between mineral surfaces and organic molecules in water control many processes in nature and in the production of modern materials. To improve the understanding of fluid-surface interactions, we investigated the interface behavior of quartz and muscovite, a model for clay minerals, in aqueous solutions where the pH and composition were controlled. We used atomic force microscopy (AFM) in chemical force mapping (CFM) mode to measure adhesion using tips functionalized with alkyl, -CH3. By combining adhesion forces measured as a function of pH, with data from streaming potential experiments and DLVO calculations, we were able to determine the surface charge density. We observed increased adhesion between the mineral surface and the hydrophobic tips as the contact time increased from 7 ms to ∼2 s. The diffusion of dissolved ions takes time, and density functional theory (DFT) calculations did not indicate a strong hydration of the mineral surfaces. Therefore, we interpret that the loss of ions from the confined space between the tip and sample is a likely explanation of the correlation between the dwell time and adhesion. The maximum adhesion increase with dwell time for muscovite, i.e., 400 ± 77 pN, was considerably larger than for quartz, 84 ± 15 pN, which fits with the different surface structure and composition of the two minerals. We propose two mechanisms to explain these results: (1) cations that are structured in the solution and on the surface remain associated at the tip-sample interface initially but diffuse away during extended contact time and (2) adventitious carbon, the organic material that comes spontaneously from air and solution, can diffuse to the tip-sample interface during contact. This material decreases the surface energy by aggregating near the alkyl tip and increases adhesion between the tip and sample.

Nickel adsorption on chalk and calcite

Nickel uptake from solution by two types of chalk and calcite was investigated in batch sorption studies. The goal was to understand the difference in sorption behavior between synthetic and biogenic calcite. Experiments at atmospheric partial pressure of CO2, in solutions equilibrated with calcite and chalk and pH ranging from 7.7 to 8.8, explored the influence of initial concentration and the amount and type of sorbent on Ni uptake. Adsorption increases with increased surface area and pH. A surface complexation model describes the data well. Stability constants for the Ni surface complex are log KNi=-1.12 on calcite and log KNi=-0.43 and -0.50 on the two chalk samples. The study confirms that synthetic calcite and chalk both take up nickel, but Ni binds more strongly on the biogenic calcite than on inorganically precipitated, synthetic powder, because of the presence of trace amounts of polysaccharides and clay nanoparticles on the chalk surface.

Predicting hydration energies for multivalent ions

We have predicted the free energy of hydration for 40 monovalent and multivalent cations and anions using density functional theory and the implicit solvent model COnductor like Screening MOdel for Real Solvents (COSMO-RS) at the Becke-Perdew (BP)/Triple zeta valence with polarization functions (TZVP) level. Agreement with experimental data for monovalent and divalent ions is good and shows no significant systematic errors. Predictions are noticeably better than with standard COSMO. The agreement with experimental data for trivalent and tetravalent ions is slightly worse and shows systematic errors. Our results indicate that quantum chemical calculations combined with COSMO-RS solvent treatment is a reliable method for treating multivalent ions in solution, provided one hydration shell of explicit water molecules is included for metal cations. The accuracy is not high enough to allow absolute predictions of hydration energies but could be used to investigate trends for several ions, thanks to the low computational cost, in particular for ligand exchange reactions.

Incorporation of monovalent cations in sulfate green rust

Green rust is a naturally occurring layered mixed-valent ferrous-ferric hydroxide, which can react with a range of redox-active compounds. Sulfate-bearing green rust is generally thought to have interlayers composed of sulfate and water. Here, we provide evidence that the interlayers also contain monovalent cations, using X-ray photoelectron spectroscopy and synchrotron X-ray scattering. For material synthesized with Na(+), K(+), Rb(+), or Cs(+), interlayer thickness derived from basal plane spacings correlates with the radius of the monovalent cation. In addition, sequential washing of the materials with water showed that Na(+) and K(+) were structurally fixed in the interlayer, whereas Rb(+) and Cs(+) could be removed, resulting in a decrease in the basal layer spacing. The incorporation of cations in the interlayer opens up new possibilities for the use of sulfate green rust for exchange reactions with both anions and cations: e.g., radioactive Cs.

Interactions of Na+, K+, Mg2+, and Ca2+ with benzene self-assembled monolayers

Interactions between cations and organic molecules are found throughout nature, from the functionality and structure of proteins in humans and animals to the exchange of ions in minerals in soil and oil reservoirs with the fluid phases. We have explored the behavior of the s-block elements that are most common in the natural world, namely, Na(+), K(+), Mg(2+), and Ca(2+). Specifically, we investigated how these ions affect the interactions between surfaces covered by self-assembled monolayers (SAMs) terminated with benzene molecules. We used a flat oxidized silicon substrate and an atomic force microscopy (AFM) tip that were both functionalized with 11-phenoxyundecane-1-thiol and measured the adhesion force between them in solutions of each of the four chloride salts. We observed that the adhesion increased in the order of the Hofmeister series: K(+) < Na(+) ≈ Mg(2+) < Ca(2+). Supplementary evidence from X-ray photoelectron spectroscopy (XPS) allowed us to conclude that K(+) binds in the benzene layers, creating a positive surface charge on the benzene-covered surfaces, thus leading to lower adhesion in KCl solutions than in pure water. Evidence suggested that Ca(2+) does not bind to the surfaces but forms bridges between the layers, leading to higher adhesion than in pure water. In Na(+) and Mg(2+) solutions, adhesion is quite similar to that in pure water, indicating a lack of interaction between these two ions and the surfaces, or at least that the interaction is too weak to be detected by our measurements. The results of our studies clearly show that even a nonpolar, hydrophobic molecule, such as benzene, has a role to play in the behavior of aqueous solutions and that it interacts differently depending on which ions are present. Even ions from the same column in the periodic table behave differently.

Change in organic molecule adhesion on α-alumina (sapphire) with change in NaCl and CaCl2 solution salinity

We investigated the adhesion of two functional groups to α-alumina as a model for the adsorption of organic molecules on clay minerals. Interactions between organic compounds and clay minerals play an important role in processes such as drinking water treatment, remediation of contaminated soil, oil recovery, and fabricating complicated nanomaterials, and there have been claims that organic compound-clay mineral interaction created the ordering that is necessary for the genesis of life. In many organisms, interaction between organic molecules and biominerals makes it possible to control the growth of bones, teeth, and shells. Adhesion of carboxylic acid, -COO(H), and pyridine, -C5H5N(H(+)), on the {0001} plane of α-alumina wafers has been investigated with atomic force microscopy (AFM) in chemical force mapping (CFM) mode. Both functional groups adhered to α-alumina in deionized water at pH < 5, and adhesion decreased as NaCl or CaCl2 concentration increased. X-ray photoelectron spectroscopy (XPS) showed that Na(+) and Ca(2+) adsorbed to the α-alumina surface at pH < 5, decreasing surface interaction with the carboxylic acid and pyridine groups. We interpret the results as evidence that the tips adhere to alumina through hydrogen bonding when only water is present. In solutions containing NaCl and CaCl2, cations are adsorbed but Cl(-) is not. When NaCl solutions are replaced by CaCl2, Ca(2+) replaces Na(+), but rinsing with ultrapure deionized water (pH 5.6) could not restore the original protonated surface. The results demonstrate that the alumina surface at pH 3 has a higher affinity for inorganic cations than for -COO(H) or -C5H5N(H(+)), in spite of the known positive surface charge of α-alumina {0001} wafers. These results demonstrate that solution salinity plays an important role in surface properties, controlling surface tension (i.e., contact angle) and adsorption affinity on α-alumina and, by analogy, on clay minerals.

Predicting the pKa and stability of organic acids and bases at an oil-water interface

We have used density functional theory and the implicit solvent model, COSMO-RS, to investigate how the acidity constant, pKa, of organic acids and bases adsorbed at the organic compound-aqueous solution interface changes, compared to its value in the aqueous phase. The pKa determine the surface charge density of the molecules that accumulate at the fluid-fluid interface. We have estimated the pKa by comparing the stability of the protonated and unprotonated forms of a series of molecules in the bulk aqueous solution and at an interface where parts of each molecule reside in the hydrophobic phase and the rest remains in the hydrophilic phase. We found that the pKa for acids is shifted by ∼1 pH unit to higher values compared to the bulk water pKa, whereas they are shifted to lower values by a similar amount for bases. Because this pKa shift is similar in magnitude for each of the molecules studied, we propose that the pKa for molecules at a water-organic compound interface can easily be predicted by adding a small shift to the aqueous pKa. This shift is general and correlates with the functional group. We also found that the relative composition of molecules at the fluid-fluid interface is not the same as in the bulk. For example, species such as carboxylic acids are enriched at the interface, where they can dominate surface properties, even when they are a modest component in the bulk fluid. For high surface concentrations of carboxylic acid groups at an interface, such as a self-assembled monolayer, we have demonstrated that the pKa depends on the degree of deprotonation through direct hydrogen bonding between protonated and deprotonated acidic headgroups.

Strontium, nickel, cadmium, and lead substitution into calcite, studied by density functional theory

We have used density functional theory to predict the ion exchange energies for divalent cations Ni(2+), Sr(2+), Cd(2+), and Pb(2+) into a calcite {10.4} surface in equilibrium with water. Exchange energies were calculated for substitution into the topmost surface layer, at the mineral-fluid interface, and into the second layer of the solid. This information can be used as an indicator for cation substitution in the bulk phase, such as for the uptake of toxic metals from the environment and the growth of secondary phases. In both the surface and in the second layer, Ni(2+), Cd(2+), and Pb(2+) substitute exothermically and Sr(2+) substitutes endothermically. Our results agree with published experimental data that demonstrate trace metal coprecipitation with calcite as a sink for Ni(2+), Cd(2+), and Pb(2+), whereas Sr(2+) has a distribution constant significantly smaller than 1. Ni(2+) substitution is favored at the mineral-fluid interface compared with bulk substitution, which also agrees with experimental data. Our results predict that Ni(2+), Cd(2+), and Pb(2+) form a stable solid solution with calcite. Successful prediction of the experimental results gives us confidence in our ability to predict the divalent cation preference for surfaces rather than for sites within the bulk crystal structure, which cannot be directly derived from experiment.

Energies of the adsorption of functional groups to calcium carbonate polymorphs: the importance of -OH and -COOH groups

The adsorption behavior of calcium carbonate is an important factor in many processes in nature, industry, and biological systems. We determined and compared the adsorption energies for a series of small molecules of different sizes and polarities (i.e., water, several alcohols, and acetic acid) on three synthetic CaCO3 polymorphs (calcite, aragonite, and vaterite). We measured isosteric heats of adsorption from vapor adsorption isotherms for 273 < T < 293 K, and we used XRD and SEM to confirm that samples did not change phase during the experiments. Density functional calculations and molecular dynamics simulations complemented the experimental results and aided interpretation. Alcohols with molecular mass greater than that of methanol bind more strongly to the calcium carbonate polymorphs than water and acetic acid. The adsorption energies for the alcohols are typical of chemisorption and indicate alcohol displacement of water from calcium carbonate surfaces. This explains why organisms favor biomolecules that contain alcohol functional groups (-OH) to control which polymorph they use, the crystal face and orientation, and the particle shape and size in biomineralization processes. This new insight is also very useful in understanding organic molecule adsorption mechanisms in soils, sediments, and rocks, which is important for predicting the behavior of mineral-fluid interactions when the challenge is to remediate contaminated groundwater aquifers or to produce oil and gas from reservoirs.

Atomic modifications by synchrotron radiation at the calcite-ethanol interface

This article reports on studies of the chemical alterations induced by synchrotron radiation at the calcite-ethanol interface, a simple model system for interfaces between minerals and more complex organic molecules containing OH groups. A combination of X-ray reflectivity and X-ray photoelectron spectroscopy of natural calcite, cleaved in distilled ethanol to obtain new clean interfaces, indicated that, during a 5 h period, the two top atomic layers of calcite, CaCO(3), transform into calcium oxide, CaO, by releasing CO(2). Also, the occupation of the first ordered layer of ethanol attached to calcite by hydrogen bonds almost doubles. Comparison between radiated and non-radiated areas of the same samples demonstrate that these effects are induced only by radiation and not caused by aging. These observations contribute to establishing a time limit for synchrotron experiments involving fluid-mineral interfaces where the polar OH group, as present in ethanol, plays a key role in their molecular structure and bonding. Also, the chemical evolution observed in the interface provides new insight into the behavior of some complex organic molecules involved in biomineralization processes.

Molecular ordering of ethanol at the calcite surface

To produce biominerals, such as shells, bones, and teeth, living beings create organic compounds that control the growth of the solid phase. Investigating the atomic scale behavior of individual functional groups at the mineral-fluid interface provides fundamental information that is useful for constructing accurate predictive models for natural systems. Previous investigations of the activity of coccolith-associated polysaccharides (CAP) on calcite, using atomic force microscopy (AFM) [Henriksen, K., Young, J. R., Bown, P. R., and Stipp, S. L. S. Palentology 2004, 43 (Part 3), 725-743] and molecular dynamics (MD) modeling [Yang, M., Stipp, S. L. S., and Harding, J. H. Cryst. Growth Des. 2008, 8 (11), 4066-4074], have suggested that OH functional groups control polysaccharide attachment. The purpose of this work was to characterize, using X-ray reflectivity (XR) combined with molecular dynamics (MD) simulations, the structuring on calcite of a layer of the simplest carbon chain molecule that contains an OH group, ethanol (CH(3)-CH(2)-OH). We found evidence that EtOH forms a highly ordered structure at the calcite surface, where the first layer molecules bond with calcite. The ethanol molecules stand up perpendicularly at the interface or nearly so. As a consequence, the fatty, CH(3) ends form a new surface, about 6 Å from the termination of the bulk calcite, and beyond that, there is a thin gap where ethanol density is low. Following is a more disordered layer that is two to three ethanol molecules thick, about 14 Å, where density more resembles that of bulk liquid ethanol. The good agreement between theory and experiment gives confidence that a theoretical approach can offer information about behavior in more complex systems.

Binding of ethanol on calcite: the role of the OH bond and its relevance to biomineralization

The interaction of OH-containing compounds with calcite, CaCO(3), such as is required for the processes that control biomineralization, has been investigated in a low-water solution. We used ethanol (EtOH) as a simple, model, OH-containing organic compound, and observed the strength of its adsorption on calcite relative to OH from water and the consequences of the differences in interaction on crystal growth and dissolution. A combination of atomic force microscopy (AFM) and molecular dynamics (MD) simulations showed that EtOH attachment on calcite is stronger than HOH binding and that the first adsorbed layer of ethanol is highly ordered. The strong ordering of the ethanol molecules has important implications for mineral growth and dissolution because it produces a hydrophobic layer. Ethanol ordering is disturbed along steps and at defect sites, providing a bridge from the bulk solution to the surface. The strong influence of calcite in structuring ethanol extends further into the liquid than expected from electrical double-layer theory. This suggests that in fluids where water activity is low, such as in biological systems optimized for biomineralization, organic molecules can control ion transport to and from the mineral surface, confining it to specific locations, thus providing the organism with control for biomineral morphology.

Interaction of ethanol and water with the {1014} surface of calcite

Molecular dynamics simulations have been used to model the interaction between ethanol, water, and the {1014} surface of calcite. Our results demonstrate that a single ethanol molecule is able to form two interactions with the mineral surface (both Ca-O and O-H), resulting in a highly ordered, stable adsorption layer. In contrast, a single water molecule can only form one or other of these interactions and is thus less well bound, resulting in a more unstable adsorption layer. Consequently, when competitive adsorption is considered, ethanol dominates the adsorption layer that forms even when the starting configuration consists of a complete monolayer of water at the surface. The computational results are in good agreement with the results from atomic force microscopy experiments where it is observed that a layer of ethanol remains attached to the calcite surface, decreasing its ability to interact with water and for growth at the {1014} surface to occur. This observation, and its corresponding molecular explanation, may give some insight into the ability to control crystal form using mixtures of different organic solvents.

Uptake and release of cerium during Fe-oxide formation and transformation in Fe(II) solutions

Fe-oxides are ubiquitous in soils and sediments and form during Fe(0) corrosion. Depending on redox conditions and solution composition, Fe-oxides such as ferrihydrite, goethite, magnetite, and green rust (GR) may form. These phases typically have high surface area and large affinity for adsorption of trace components. Further, Fe(II)-Fe(III) (hydr)oxides are redox active. Cerium, a member of the lanthanide family, can be used as an analogue for the tri- and tetra-valent actinides found in radioactive waste, expected to be stored in subsurface repositories. In experiments with ferrihydrite, Ce(III) was effectively scavenged from Fe(II)-bearing solutions within 5 min at pH 7. During transformation of ferrihydrite to green rust, however, all Ce(III) was released to solution. By varying initial solution Fe(II):Fe(III) ratio, magnetite and goethite formed together with GR(Na,SO(4)), resulting in decreased Ce(III) release. X-ray photoelectron spectroscopy revealed Ce(III) adsorbed on magnetite. When Fe-oxides were synthesized by air oxidation of Fe(II) solutions at pH 7, GR(Na,SO(4)) played a catalytic role in the oxidation of Ce(III) to Ce(IV) by O(2), removing more than 90% of the dissolved Ce. Transmission electron microscopy revealed that it formed discrete nanocrystals of CeO(2(s)). These results demonstrate that Fe-oxide interaction with radionuclides is likely to depend strongly on the local redox conditions. By analogy with Ce, the trivalent actinides are not expected to be sequestered by preformed GR in anoxic environments. Our results also suggest that trivalent actinides and lanthanides are released when dissimilatory iron reduction of Fe(III)-oxides leads to GR formation However, under oxidizing conditions, GR may influence radionuclide mobility by catalyzing their transformation to a higher oxidation state.

Identification of green rust in groundwater

Green rust, a family of Fe(II), Fe(III) layered double hydroxides, is believed to be present in environments close to the Fe(II)/ Fe(III) transition zone. Attempts to identify members of this family in nature have proven difficult because the material is oxidized after only a few minutes exposure to air. In this paper, we present a sampling method for capturing green rust so it is not oxidized. We then we used the method to identify the compound in a groundwater sample taken below the water table from fractures in granite. X-ray diffraction patterns were weak, but clearly identical to those of synthetic GR(CO3), the green rustfamily memberwhere carbonate and water occupy the interlayer between the iron-hydroxide layers. The method was then tested on samples taken from an artesian well and a deep underground experimental station, both within the Fe(II)/ Fe(III) redox zone. In both cases, GR(CO3), could be identified. Currently, transport models for predicting the behavior of contaminants in groundwater do not include parameters for green rust This work demonstrates they should.

Fluoride removal by calcite: evidence for fluorite precipitation and surface adsorption

Fluoride contamination of groundwater, both anthropogenic and natural, is a major problem worldwide. In this study, fluoride removal by crushed limestone (99% pure calcite) was investigated by batch studies and surface-sensitive techniques from solutions with fluoride concentrations from 150 micromol/L (3 mg/L) to 110 mM (approximately 2100 mg/L). Surface-sensitive techniques, including atomic force microscopy (AFM) and X-ray photoelectron spectroscopy (XPS) as well as zeta potential measurements, confirm that, in addition to precipitation reactions, adsorption of fluoride also occurs. Results indicate that fluoride adsorption occurs immediately over the entire calcite surface with fluorite precipitating at step edges and kinks, where dissolved Ca2+ concentration is highest. The PHREEQ geochemical model was applied to the observed data and indicates that existing models, especially at low fluoride concentrations and high pH (>7.5) are not equipped to describe this complex system, largely because the PHREEQ model includes only precipitation reactions, whereas a combination of adsorption and precipitation parameters are required.

Modeling steps and kinks on the surface of calcite

This work presents modeling results on the cleavage face of calcite as well as on steps and isolated kinks on this face. We used static lattice energy minimization and interatomic potentials fitted to bulk properties. The energy needed to cleave a bulk calcite crystal along the [1 0 (-)1 4] plane was calculated to be 0.59 J m(-2) in agreement with previous studies using the same potentials. The perfect surface reconstructs in the top few atomic layers, but its symmetry corresponds to the bulk termination. By contrast, the (1 0 (-)1 4) surface with cleavage steps present reconstructs to form a (2 x 1) super cell. This may help explain experimental observations of (2 x 1) symmetry on calcite surfaces. The energy required to form a monatomic obtuse step is calculated to be 1.3 x 10(-10) J m(-1) and for the acute step, 2.4 x 10(-10) J m(-1), suggesting that obtuse steps dominate on cleaved surfaces. Along the two types of steps, a total of 16 kink geometries exist. We calculated kink defect energy with two different approaches: one where kink pairs were added onto infinitely long steps and one where kinks were placed inside pits on a cleavage surface. Calculations on infinitely long steps show that for vacuum conditions, kink pairs possess roughly identical formation energy, about 1.2-2.2 eV, so based on energetics one cannot expect significant differences in kink site frequency

Eu3+ uptake by calcite: preliminary results from coprecipitation experiments and observations with surface-sensitive techniques

A lack of information in databases for contamination risk assessment about the transport behaviour of the trivalent f-orbital elements in groundwater systems where calcite is at equilibrium motivated this study of Eu(3+) uptake. The free drift technique was used to examine the effects of Eu(3+) concentration, presence of Na(+) or K(+) and temperature, as well as calcite nucleation and precipitation kinetics, on the partitioning of calcite. Changes in surface composition and morphology resulting from exposure of single crystals of Iceland spar to Eu(3+)-bearing solutions were observed with X-ray photoelectron spectroscopy (XPS) and atomic force microscopy (AFM). First results confirm that calcite has high affinity for Eu(3+). Rates of nucleation and precipitation strongly affect the extent of uptake but the presence of Na(+) and K(+) has no effect, suggesting formation of solid solution as CaCO(3)-EuOHCO(3). Surface-sensitive techniques prove that Eu(3+) is adsorbed to calcite even when the surface is dissolving and adsorption is not accompanied by precipitation of a separate Eu(3+)-solid phase. Adsorbed Eu modifies calcite's dissolution behaviour, roughening terraces and rounding step edges, and producing surface morphology where some surface sites appear blocked. Results imply that Eu(3+) concentrations in natural calcites are limited by Eu(3+) availability rather than by a lack of ability to fit into calcite's atomic structure. This behaviour can probably be expected for other trivalent rare Earth elements (REE), actinides and fission products whose behaviour is similar to that of Eu(3+). These elements are likely to be incorporated within the calcite bulk in systems where it is precipitating and the demonstrated strong partitioning ensures some uptake even where calcite is at or under saturation.