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.     

Americium interaction with calcite and aragonite surfaces in seawater

Observations of the distribution of ²⁴¹Am in the marine environment indicate that Am has a high affinity for solid surfaces. The adsorption of Am onto calcite and aragonite surfaces from seawater and related solutions has been studied, in order to establish the interaction of Am with a major component of many marine sediments. Results indicate that Am is rapidly and strongly adsorbed. This occurs even when both dissolved Am concentrations and solid to solution ratios are low. The minimum value for $\text{K}{}_{\mathrm{D}}$ determined is 2 × 10⁵. Measurements of reaction kinetics established that Am is adsorbed from seawater at 40 times the rate per unit surface area on synthetic aragonite that it is on synthetic calcite. Approximately 15% of the difference is attributable to epitaxial influences, with the remainder being due to enhanced site competition by Mg on calcite relative to aragonite. The adsorption rate is first order with respect to Am concentration, but follows approximately the square root of the solid surface area to solution volume ratio.Adsorption rate of Am on biogenic aragonite and Mg-calcites are, within a given particle size range, close to equal. It is not possible to normalize these adsorption rates to surface area due to the differing microporous structure of biogenic carbonates. The Am adsorption rates on a shallow water calcium carbonate-rich sediment gave results which were predicted from, its mineralogie mixture of components.     

Interaction of calcium dioleate collector colloids with calcite and fluorite surfaces as revealed by AFM force measurements and molecular dynamics simulation

Spherical calcium dioleate particles (∼ 10 μm in diameter) were used as AFM (atomic force microscope) probes to measure interaction forces of the collector colloid with calcite and fluorite surfaces. The attractive AFM force between the calcium dioleate sphere and the fluorite surface is strong and has a longer range than the DLVO (Derjaguin–Landau–Verwey–Overbeek) prediction. The repulsive AFM force between the calcium dioleate sphere and the fluorite surface does not agree with the DLVO prediction. Consideration of non-DLVO forces, including the attractive hydrophobic force, was necessary to explain the experimental results. The non-DLVO interactions considered were justified by the different interfacial water structures at fluorite– and calcite–water interfaces as revealed by the numerical computation experiments using molecular dynamics simulation. The density of interfacial water at the fluorite surface is low and the fluorite surface is not strongly wetted by water molecules. In contrast to the water at the fluorite surface, water molecules at the calcite surface form tightly packed monolayer structures and the calcite surface is extensively hydrated by water molecules. The interfacial water structure agrees with the AFM force measurements and the flotation recovery data. The strong attraction between the calcium dioleate colloid and the fluorite surface, and the moderately wetted fluorite surface by water molecules explain the better flotation response of fluorite with the oleate collector colloid.     

Effects of magnesium ions on near-equilibrium calcite dissolution: Step kinetics and morphology

Dissolution kinetics at the aqueous solution-calcite interface at 50 °C were investigated using in situ atomic force microscopy (AFM) to reveal the influence of magnesium concentration and solution saturation state on calcite dissolution kinetics and surface morphology. Under near-equilibrium conditions, dissolved Mg²⁺ displayed negligible inhibitory effects on calcite dissolution even at concentrations of . Upon the introduction of , the solution saturation state with respect to calcite, , acted as a “switch” for magnesium inhibition whereby no significant changes in step kinetics were observed at Ω_calcite<0.2, whereas a sudden inhibition from Mg²⁺ was activated at Ω_calcite?0.2. The presence of the Ω-switch in dissolution kinetics indicates the presence of critical undersaturation in accordance with thermodynamic principles. The etch pits formed in solutions with exhibited a unique distorted rhombic profile, different from those formed in Mg-free solutions and in de-ionized water. Such unique etch pit morphology may be associated with the anisotropy in net detachment rates of counter-propagating kink sites upon the addition of Mg²⁺.     

Deposition rates of polysilicic acid with up to 10 M calcium ions

Cementitious materials used for radioactive waste repository construction complicate the performance assessment of radioactive waste systems because the use of cement may greatly alter the pH (8–13) of groundwater and release constituents such as calcium ions. Under such conditions, it is important to clarify also the dynamic behavior of silica (silicic acid), in order to evaluate the alteration in the chemical and physical properties of the fractured layer or the host rock surrounding the repository. Since silica undergoes polymerization, precipitation or dissolution depending on the pH and/or temperature, the behavior of silica would be greatly complicated in the presence of other ions. This study is focused on the deposition rates of polysilicic acid and soluble silicic acid with up to 10⁻³ M Ca ions. In the experiment, Na₂SiO₃ solution (250 mL, pH > 10, 298 K) was poured into a polyethylene vessel containing amorphous silica powder (0.5 g), and a buffer solution, HNO₃, and CaNO₃ as Ca ions were sequentially added into the vessel. The pH of the solution was set to 8. The silica, initially in a soluble form at pH > 10 (1.4 × 10⁻² M), became supersaturated and either deposited on the solid surface or changed into the polymeric form. Then the concentrations of both poly- and soluble silicic acid were monitored over a 40-day period. The decrease of polysilicic acid became slow with an increase in the concentration of Ca ions in the range of up to 10⁻³ M. In general, the addition of electrolytes to a supersaturated solution accelerates the aggregation and precipitation of polymeric species. However, the experimental result showed that polysilicic acid in the presence of Ca ions is apparently stable in solution, compared with that under a Ca-free condition. On the other hand, the concentration of soluble silicic acid in the presence of Ca ions immediately became metastable, that is, slightly higher than the solubility of soluble silicic acid. Its dynamic behavior was similar to that in the Ca-free condition.     

Adsorption of palmitic acid on calcite

The adsorption of palmitic acid by calcite was determined as a function of palmitic acid concentration and salinity. Adsorption isotherms were generally of the S-type, indicating probable hydrophobic interaction of adsorbed species. The effect of increased salinity was to increase adsorption from 0–25 ppt, then to decrease adsorption from 25–35 parts per thousand. The adsorption increase results from decreased solubility, whereas the adsorption decrease may be related to micelle formation. Dissolved Mg²⁺ was found to strongly inhibit fatty acid adsorption. Desorption of surface-bound phosphate was associated with palmitic acid adsorption. The relationship between adsorbed Mg²⁺, phosphate and palmitic acid suggests investigation of adsorption behavior of natural phospholipid compounds.     

Interactions of biopolymers with silica surfaces: Force measurements and electronic structure calculation studies

Pull-off forces were measured between a silica colloid attached to an atomic force microscope (AFM) cantilever and three homopolymer surfaces representing constituents of extracellular polymeric substances (EPS). The pull-off forces were −0.84 (±0.16), −0.68 (±0.15), and −2.37 (±0.31) nN as measured in water for dextran, phosphorylated dextran, and poly-l-lysine, respectively. Molecular orbital and density functional theory methods (DFT) were applied to analyze the measured pull-off forces using dimer clusters representing interactions between the three polymers and silica surfaces. Binding energies for each dimer were calculated with basis set superposition error (BSSE) and interpolated using corrections for silica surface hydroxyl density and silica charge density. The binding energies were compared with the normalized pull-off forces with the effective silica surface area contacting the polymer surfaces. The predicted binding energies at a −0.064 C/m² silica surface charge density corresponding to circum-neutral pH were −0.055, −0.029, and −0.338 × 10⁻¹⁸ J/nm² for the dimers corresponding to the silica surface with dextran, phosphorylated dextran, and poly-l-lysine, respectively. Polarizable continuum model (PCM) calculations with different solvents, silanol vibrational frequency calculations, and orbital interaction analysis based on natural bonding orbital (NBO) showed that phosphate groups formed stronger H-bonds with neutral silanols than hydroxyl and amino functional groups of polymers, implying that phosphate containing polymers would play important roles in EPS binding to silica surfaces.     

Coprecipitation of alkali metal ions with calcium carbonate

The coprecipitation of alkali metal ions (Li⁺, Na⁺, K⁺ and Rb⁺) with calcium carbonate has been studied experimentally and the following results have been obtained:

• 1.(1) Alkali metal ions are more easily coprecipitated with aragonite than with calcite.
• 2.(2) The relationship between the amounts of alkali metal ions coprecipitated with aragonite and their ionic radii shows a parabolic curve with a peak located at Na⁺ which has approximately the same ionic radius as Ca²⁺.
• 3.(3) However, the amounts of alkali metal ions coprecipitated with calcite decrease with increasing ionic radius of alkali metals.
• 4.(4) Our results support the hypothesis that
  • 4.1.(a) alkali metals are in interstitial positions in the crystal structure of calcite and do not substitute for Ca²⁺ in the lattice, but
  • 4.2.(b) in aragonite, alkali metals substitute for Ca²⁺ in the crystal structure.
• 5.(5) Magnesium ions in the parent solution increase the amounts of alkali metal ions (Li⁺, Na⁺, K⁺ and Rb⁺) coprecipitated with calcite but decrease those with aragonite.
• 6.(6) Sodium-bearing aragonite decreases the incorporation of other alkali metal ions (Li⁺, K⁺ and Rb⁺) into the aragonite.

     

Effects of phosphate ions on intergranular pressure solution in calcite: An experimental study

Intergranular pressure solution (IPS) is a coupled chemical-mechanical process of widespread importance that occurs during diagenesis and low-temperature deformation of sedimentary rocks. Laboratory experiments on IPS in halite, quartz, and calcite have largely concentrated on the mechanical aspects of the process. In this study, we report the effects of pore fluid chemistry, specifically varying phosphate ion concentration, on the mechanical compaction by IPS of fine-grained calcite powders at room temperature and 1 to 4 MPa applied effective stress. Phosphate was investigated because of its importance as a biogenic constituent of sea and pore waters. Increasing the pore fluid phosphate concentration from 0 to 10⁻³ mol/L systematically reduced compaction strain rates by up to two orders of magnitude. The sensitivity of the compaction strain rate to phosphate concentration was the same as the sensitivity of calcite precipitation rates to the addition of phosphate ions reported in the literature, suggesting that the rate of IPS in phosphate-bearing samples was controlled by calcite precipitation on pore walls. The results imply that IPS and associated porosity/permeability reduction rates in calcite sediments may be strongly reduced when pore fluids are enriched in phosphates, for example, through high biologic productivity or a seawater origin. Future modeling of IPS-related processes in carbonates must therefore take into account the effects of pore fluid chemistry, specifically the inhibition of interfacial reactions.     

Thickness and structure of the water film deposited from vapour on calcite surfaces

Synchrotron X-ray reflectivity (SXR) was used to measure the thickness of the water film that adsorbs on a {} cleavage surface of calcite (CaCO₃) in a sample chamber where relative humidity could be controlled within the range from <4% to 90%. Gases used to carry water vapour were initially either 100% N₂ or 100% CO₂. The product water film was remarkably constant in thickness at 15.5 Å (±1 Å) and independent of humidity. When N₂ was used as the carrier gas, this film displayed a gap in its electron density at between 0.6 and 2 Å distance from the calcite surface, depending on humidity. This implies that a change in the arrangement of water molecules occurs in direct proximity to the surface. This electron density discontinuity was measurably further from the calcite surface, at 3.4 Å, when CO₂ was used as the carrier gas. Except for this thin low density region proximate to the calcite surface, the density of the adsorbed water layer was 0.9 g cm⁻³, therefore suggesting a significant degree of ordering. Atomic force microscopy (AFM) images were completed in conjunction with the SXR measurements on similarly prepared calcite cleavage surfaces. AFM showed that terraces may be atomically flat over 1 μm or more. SXR corroborated this observation, with results showing that carefully cleaved surfaces had a starting root mean square (r.m.s.) roughness of ∼1.2 Å. Diffuse scatter measurements constrained the correlation lengths of these surfaces to be at least 18,000 Å. For comparison with the cleaved samples, a surface was also prepared by chemo-mechanical Syton polishing. This surface gave an r.m.s roughness by SXR that was an order of magnitude higher, equal to 12.1 Å. In this case, diffuse scatter resolved a correlation length of 950 Å, and revealed a fractal dimension that was higher than for the cleaved surface. On Syton polished samples, the water film determined by SXR was about twice as thick as for freshly cleaved surfaces, with a density of 1.0 g cm⁻³, equal to that of bulk water. However, surface roughness was too large to allow resolution of any gap in the electron density within the water layer proximate to the solid surface. Our AFM observations also confirm previous reports of calcite surface recrystallization. The electron density of the solid surface determined by SXR is indistinguishable from that of calcite, indicating that any material recrystallized within the adsorbed water film is compositionally indistinguishable from the calcite substrate.     

Volume and grain boundary diffusion of calcium in natural and hot-pressed calcite aggregates

 Calcium self-diffusion rates in natural calcite single crystals were experimentally determined at 700 to 900° C and 0.1 MPa in a stream of CO₂. Diffusion coefficients (D) were determined from ⁴²Ca concentration profiles measured with an ion microprobe. The Arrhenius parameters yield an activation energy (Q)=382±37 kJ/mol and pre-exponential factor (D₀)=0.13 m²/s, and there is no measurable anisotropy. Calcium grain boundary diffusion rates were experimentally determined in natural (Solnhofen) limestone and hot-pressed calcite aggregates at 650° to 850° C and 0.1 to 100 MPa pressure. The Solnhofen limestone was first pre-annealed for 24 h at 700° C and 100 MPa confining pressure under anhydrous conditions to produce an equilibrium microstructure for the diffusion experiments. Values for the product of the grain boundary diffusion coefficient (D′) and the effective grain boundary diffusion width (δ) were determined from ⁴²Ca concentration profiles measured with an ion microprobe. The results show that there is no measurable difference between D′δ values obtained for pre-annealed Solnhofen samples at 0.1 and 100 MPa or between hot-pressed calcite aggregates and pre-annealed Solnhofen samples. The temperature dependence for calcium grain boundary diffusion in Solnhofen samples annealed at 0.1 MPa is described by the Arrhenius parameters D^′ ₀δ=1.5×10⁻⁹ m³/s and Q=267±47 kJ/mol. Comparison of the results of this study with previously published data show that calcium is the slowest volume diffusing species in calcite. The calcium diffusivities measured in this study place constraints on several geological processes that involve diffusive mass transfer including diffusion-accommodated mechanisms in the deformation of calcite rocks. Received: 19 December 1994/Accepted: 30 June 1995     

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.     

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.     

Friction characteristics of Cd-rich carbonate films on calcite surfaces: implications for compositional differentiation at the nanometer scale

Lateral Force Microscopy (LFM) studies were carried out on cleaved calcite sections in contact with solutions supersaturated with respect to otavite (CdCO₃) or calcite-otavite solid solutions (SS) as a means to examine the potential for future application of LFM as a nanometer-scale mineral surface composition mapping technique. Layer-by-layer growth of surface films took place either by step advancement or by a surface nucleation and step advancement mechanisms. Friction vs. applied load data acquired on the films and the calcite substrate were successfully fitted to the Johnson Kendall Roberts (JKR) model for single asperity contacts. Following this model, friction differences between film and substrate at low loads were dictated by differences in adhesion, whereas at higher load they reflect differences in contact shear strength. In most experiments at fixed load, the film showed higher friction than the calcite surface, but the friction-load dependence for the different surfaces revealed that at low loads (0–40 nN), a calcian otavite film has lower friction than calcite; a result that is contrary to earlier LFM reports of the same system. Multilayer films of calcian-otavite displayed increasing friction with film thickness, consistent with the expectation that the film surface composition will become increasingly Cd-rich with increasing thickness. Both load- and thickness-dependence trends support the hypothesis that the contact shear strength correlates with the hydration enthalpy of the surface ions, thereby imparting friction sensitivity in the LFM to mineral-water interface composition.     

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.     

Elasticity of calcite: thermal evolution

The elastic properties of calcite have been determined by Brillouin spectroscopy for temperatures up to 600 °C. The results reveal that the variations of the aggregate bulk (K _VRH) and shear (G _VRH) moduli of calcite with respect to temperature can be approximately expressed as follows: $$\begin{aligned} K_{{{\text{VRH}}}} ({\text{GPa}}) & = 79.57-0.0230\;T\, (T\;{\text{in}}\;^{^\circ } {\text{C}}) \\G_{{{\text{VRH}}}} ({\text{GPa}}) & = 32.23 - 0.0097\;T. \\\end{aligned}$$ This indicates a nearly constant Poisson’s ratio (0.322) for calcite from 22 to 600 °C. A further analysis shows that the compressibility along the c axis (β _||) and that perpendicular to the c axis have comparable contributions to the volume compressibility of calcite, although the contribution of β _|| decreases with an increase in the temperature.     

Nanoscale phenomena during the growth of solid solutions on calcite {101¯4} surfaces

This work deals with the growth behaviour of calcite {101¯4} surfaces in contact with multicomponent aqueous solutions containing divalent cations (Ba²⁺, Sr²⁺, Mn²⁺, Cd²⁺, or Mg²⁺). The result is the formation of solid solutions, with calcite or aragonite as one of the end-members. In situ atomic force microscopy has revealed a wide variety of surface phenomena occurring during the formation of these solid solutions. Among them are: (1) the thickening of growth steps and the subsequent dissolution of surfaces followed by the nucleation of secondary three-dimensional nuclei on calcite surfaces, (2) the transition between growth mechanisms, (3) the formation of an epitaxial layer that armours the substrate from further dissolution and (4) the inhibitory effect of the newly formed surface on the subsequent growth (template effect). The two last phenomena can considerably limit coprecipitation as an effective mechanism for divalent metal uptake. All the phenomena described are a consequence of the interplay between thermodynamics, supersaturation of the aqueous solution with respect to the possible solid solutions and the crystallographic control of the surfaces on the cation incorporation, and indicates that there are many differences between the crystal growth of solid solutions and phases with fixed composition.