Free Energy of Adsorption of Water and Metal Ions on the {101̄4} Calcite Surface

Adsorption Properties of Triethylene Glycol on a Hydrated {101̅4} Calcite Surface and Its Effect on Adsorbed Water

Molecular dynamics (MD) and Born-Oppenheimer MD (BOMD) simulations were employed to investigate adsorption of aqueous triethylene glycol (TEG) on a hydrated {101̅4} calcite surface at 298 K. We analyzed the orientation of TEG adsorbed on calcite, as well as the impact of TEG on the water density and adsorption free energy. The adsorption energies of TEG, free energy profiles for TEG, details of hydrogen bonding between water and adsorbed TEG, and dihedral angle distribution of adsorbed TEG were estimated. We found that while the first layer of water was mostly unaffected by the presence of adsorbed TEG, the density of the second water layer was decreased by 71% at 75% surface coverage of TEG. TEG primarily attached to the calcite surface via two adjacent adsorption sites. Hydrogen bonds between water and adsorbed TEG in the second layer almost exclusively involved the hydroxyl oxygen of TEG. The adsorption energy of TEG on calcite in a vacuum environment calculated by classical MD amounted to 217 kJ/mol, which agreed very well with estimates found by using BOMD. Adsorption on hydrated calcite yielded a drastically lower value of 33 kJ/mol, with the corresponding adsorption free energy of 55.3 kJ/mol, giving an entropy increase of 22.3 kJ/mol due to adsorption. We found that the presence of TEG resulted in a decreased magnitude of the adsorption free energy of water, thus decreasing the calcite wettability. This effect can have a profound effect on oil and gas reservoir properties and must be carefully considered when evaluating the risk of hydrate nucleation.

Infrared Signature of the Cation-π Interaction between Calcite and Aromatic Hydrocarbons

The cation-π interaction is proposed as an important mechanism for the adsorption of aromatic hydrocarbons having non-zero quadrupole moments by mineral surfaces. Direct evidence supporting such a mechanism is, however, limited. Using the model mineral calcite, we probe the cation-π interaction with adsorbed benzene, toluene, and ethylbenzene (BTE) molecules using attenuated total reflectance Fourier transform infrared spectroscopy. We show that the presence of calcite increases the energy required to excite the synchronized bending of aromatic C-H bonds of BTE molecules. The unique conformation of this vibrational mode indicates that the planar aromatic rings of BTE molecules are constrained in a tilted face-down position by the cation-π interaction, as further confirmed by density functional theory calculations. Our results suggest that the shift of the excitation energy of the aromatic C-H bending may be used as an infrared signature for the cation-π interaction occurring on mineral surfaces.

Modelling the effects of salt solutions on the hydration of calcium ions

Classical molecular dynamics simulations of several aqueous alkali halide salt solutions have been used to determine the effect of electrolytes on the structure of water and the hydration properties of calcium ions. Compared with the simulations of Ca(2+) ions in pure liquid water, the frequency of water exchange in the first hydration shell of calcium, which is a fundamental process in controlling the reactivity of calcium(ii) aqua-ions, is drastically reduced in the presence of other electrolytes in solution. The strong stabilization of the hydration shell of Ca(2+) occurs not only when the halide anions are directly coordinated to calcium, but also when the alkali and halide ions are placed at or outside the second coordination shell of Ca(2+), suggesting that the reactivity of the first solvation shell of the calcium ion can be influenced by the specific affinity of other ions in solution for the water molecules coordinated to Ca(2+). Analysis of the hydrogen-bonded structure of water in the vicinity of the calcium ion shows that the average number of hydrogen bonds per water molecules, which is 1.8 in pure liquid water, decreases as the concentration of alkali-halide salts in solution increases, and that the temporal fluctuations of hydrogen bonds are significantly larger than those obtained for Ca(2+) in pure liquid water. This effect has been explained in terms of the dynamics of reorganization of the O-H···X(-) (X = F, Cl and Br) hydrogen bond. This work shows the importance of solution composition in determining the hydrogen-bonding network and ligand-exchange dynamics around metal ions, both in solution and at the mineral-water interfaces, which in turn has implications for interactions occurring at the mineral-water interface, ultimately controlling the mobilization of ions in the environment as well as in industrial processes.

Atomic Resolution of Calcium and Oxygen Sublattices of Calcite in Ambient Conditions by Atomic Force Microscopy Using qPlus Sensors with Sapphire Tips

Characterization and imaging at the atomic scale with atomic force microscopy in biocompatible environments is an ongoing challenge. We demonstrate atomically resolved imaging of the calcite (101̅4) surface plane using stiff quartz cantilevers ("qPlus sensors", stiffness k = 1280 N/m) equipped with sapphire tips in ambient conditions without any surface preparation. With 10 atoms in one surface unit cell, calcite has a highly complex surface structure comprising three different chemical elements (Ca, C, and O). We obtain true atomic resolution of calcite in air at relative humidity ranging from 20% to 40%, imaging atomic steps and single atomic defects. We observe a great durability of sapphire tips with their Mohs hardness of 9, only one step below diamond. Depending on the state of the sapphire tip, we resolve either the calcium or the oxygen sublattice. We determine the tip termination by comparing the experimental images with simulations and discuss the possibility of chemical tip identification in air. The main challenges for imaging arise from the presence of water layers, which form on almost all surfaces and have the potential to dissolve the crystal surface. Frequency shift versus distance spectra show the presence of at least three ordered hydration layers. The measured height of the first hydration layer corresponds well to X-ray diffraction data and molecular dynamic simulations, namely, ∼220 pm. For the following hydration layers we measure ∼380 pm for the second and third layer, ending up in a total hydration layer thickness of at least 1 nm. Understanding the influence of water layers and their structure is important for surface segregation, surface reactions including reconstructions, healing of defects, and corrosion.

Atomistic molecular dynamics simulations of carbohydrate-calcite interactions in concentrated brine

We report atomistic molecular dynamics simulations to study the interactions of a model carbohydrate monomer (Glucopyranose) and calcite slabs in brine. We show that the interactions between the sugar molecules and the mineral decrease with increasing salinity. The decrease is due to the formation of salt layers on the calcite surfaces, which screen the carbohydrate-calcite hydrogen bonding. This screening effect depends on the affinities of calcite surface to specific ions as well as to the carbohydrate molecules.

Simulating preferential sorption of tartrate on prismatic calcite surfaces

Calculated binding energies of favorable adsorption configurations emphasize the importance of surface charge/energetics, structural match and water layers in mineral–organic interactions.

Comparison of select polarizable and non-polarizable water models in predicting solvation dynamics of water confined between MgO slabs

We present a molecular dynamics simulation study in which we compare and contrast the performance of a polarizable shell water potential model and non-polarizable water force field-extended simple point charge (SPC/EF) model in predicting the solvation dynamics of confined water molecules sandwiched between MgO(100) slabs. Structural features based on radial distribution functions, atomic density profiles, adsorption patterns, orientational ordering and dynamical correlations such as diffusional characteristics, hydrogen bonding lifetimes and residence probabilities are used as metrics for comparison. The simulations yield significant ordering of water molecules in the two layers adjacent to the oxide interface and the extent of ordering decreases with increasing distance from the oxide-water interface. These results elucidate that the dependence of local ordering and solvation dynamics on the molecular geometry and charge distribution, observed for typical three- and four-site water models, is generally lost for confined water if polarization is explicitly included. While the interfacial water structure predicted by the polarizable and non-polarizable models are similar, the confinement and interface proximity effects on the solvation dynamics are seen to be more pronounced for polarizable water models in comparison to non-polarizable ones. The study also shows that the polarizable water model over predicts the orientational order and under predicts the transport properties of confined water. In addition, analysis of the orientational preferences and hydrogen bonding characteristics of water near oxide interfaces suggests a higher degree of tetrahedral disorder in the polarizable shell compared to the non-polarizable SPC/E flexible model. The origin of the differences in solvation behavior of confined water between oxide slabs is analyzed based on the energetic contributions of the dispersive and electrostatic terms in the two force fields. Our findings suggest some new considerations regarding the role of polarization terms in predicting confinement and interface proximity effects that may guide future development of reliable polarizable water models for confined liquids.

Direct visualization of single ions in the Stern layer of calcite

Calcite is among the most abundant minerals on earth and plays a central role in many environmental and geochemical processes. Here we used amplitude modulation atomic force microscopy (AFM) operated in a particular regime to visualize single ions close to the (1014) surface of calcite in solution. The results were acquired at equilibrium, in aqueous solution containing different concentrations of NaCl, RbCl, and CaCl(2). The AFM images provide a direct and atomic-level picture of the different cations adsorbed preferentially at certain locations of the calcite-water interface. Highly ordered water layers at the calcite surface prevent the hydrated ions from directly interacting with calcite due to the energy penalty incurred by the necessary restructuring of the ions' solvation shells. Controlled removal of the adsorbed ions from the interface by the AFM tip provides indications about the stability of the adsorption site. The AFM results show the familiar "row pairing" of the carbonate oxygen atoms, with the adsorbed monovalent cations located adjacent to the most prominent oxygen atoms. The location of adsorbed cations near the surface appears better defined for monovalent ions than for Ca(2+), consistent with the idea that Ca(2+) ions remain further away from the surface of calcite due to their larger hydration shell. The precise distance between the different hydrated ions and the surface of calcite is quantified using MD simulation. The preferential adsorption sites found by MD as well as the ion residence times close to the surface support the AFM findings, with Na(+) ions dwelling substantially longer and closer to the calcite surface than Ca(2+). The results also bring new insights into the problem of the Stern and electrostatic double layer at the surface of calcite, showing that parameters such as the thickness of the Stern layer can be highly ion dependent.

Microsolvation and hydration enthalpies of CaC₂O₄(H₂O) n (n=0-16) and C₂O₄²⁻(H₂O) n (n=0-14): an ab initio study

We studied hydrated calcium oxalate and its ions at the restricted Hartree-Fock RHF/6-31G* level of theory. Performing a configurational search seems to improve the fit of the HF/6-31G* level to experimental data. The first solvation shell of calcium oxalate contains 13 water molecules, while the first solvation shell of oxalate ion is formed by 14 water molecules. The first solvation shell of Ca(II) is formed by six water molecules, while the second shell contains five. At 298.15 K, we estimate the asymptotic limits (infinite dilution) of the total standard enthalpies of hydration for Ca(II), oxalate ion and calcium oxalate as -480.78, -302.78 and -312.73 kcal mol(-1), resp. The dissociation of hydrated calcium oxalate is an endothermic process with an asymptotic limit of +470.84 kcal mol(-1).

Ion desolvation as a mechanism for kinetic isotope fractionation in aqueous systems

Molecular dynamics simulations show that the desolvation rates of isotopes of Li(+), K(+), Rb(+), Ca(2+), Sr(2+), and Ba(2+) may have a relatively strong dependence on the metal cation mass. This inference is based on the observation that the exchange rate constant, k(wex), for water molecules in the first hydration shell follows an inverse power-law mass dependence (k(wex) ∝ m(-γ)), where the coefficient γ is 0.05 ± 0.01 on average for all cations studied. Simulated water-exchange rates increase with temperature and decrease with increasing isotopic mass for each element. The magnitude of the water-exchange rate is different for simulations run using different water models [i.e., extended simple point charge (SPC/E) vs. four-site transferrable intermolecular potential (TIP4P)]; however, the value of the mass exponent γ is the same. Reaction rate theory calculations predict mass exponents consistent with those determined via molecular dynamics simulations. The simulation-derived mass dependences imply that solids precipitating from aqueous solution under kinetically controlled conditions should be enriched in the light isotopes of the metal cations relative to the solutions, consistent with measured isotopic signatures in natural materials and laboratory experiments. Desolvation effects are large enough that they may be a primary determinant of the observed isotopic fractionation during precipitation.

Atomistic simulations of calcium uranyl(VI) carbonate adsorption on calcite and stepped-calcite surfaces

Adsorption of actinyl ions onto mineral surfaces is one of the main mechanisms that control the migration of these ions in environmental systems. Here, we present computational classical molecular dynamics (MD) simulations to investigate the behavior of U(VI) in contact with different calcite surfaces. The calcium-uranyl-carbonate [Ca(2)UO(2)(CO(3))(3)] species is shown to display both inner- and outer-sphere adsorption to the flat {101̅4} and the stepped {314̅8} and {31̅2̅16} planes of calcite. Free energy calculations, using the umbrella sampling method, are employed to simulate adsorption paths of the same uranyl species on the different calcite surfaces under aqueous condition. Outer-sphere adsorption is found to dominate over inner-sphere adsorption because of the high free energy barrier of removing a uranyl-carbonate interaction and replacing it with a new uranyl-surface interaction. An important binding mode is proposed involving a single vicinal water monolayer between the surface and the sorbed complex. From the free energy profiles of the different calcite surfaces, the uranyl complex was also found to adsorb preferentially on the acute-stepped {314̅8} face of calcite, in agreement with experiment.

Diffusion and adsorption of uranyl carbonate species in nanosized mineral fractures

Atomistic simulations were performed to study the diffusion and adsorption of Ca(2)UO(2)(CO3)3 and of some of its constituent species, i.e., UO(2)2+, CO(3)2–, and UO(2)CO3, in feldspar nanosized fractures. Feldspar is important to uranium remediation efforts at the U.S. Department of Energy Hanford site as it has been found in recent studies to host contaminants within its intragrain fractures. In addition, uranyl carbonate species are known to dominate U(VI) speciation in conditions relevant to the Hanford site. Molecular dynamics (MD) simulations showed that the presence of the feldspar surface diminishes the diffusion coefficients of all of the species considered in this work and that the diffusion coefficients do not reach their bulk aqueous solution values in the center of a 2.5 nm fracture. Moreover, the MD simulations showed that the rate of decrease in the diffusion coefficients with decreasing distance from the surface is greater for larger adsorbing species. Free energy profiles of the same species adsorbing on the feldspar surface revealed a large favorable free energy of adsorption for UO(2)2+ and UO(2)CO3, which are able to adsorb to the surface with their uranium atom directly bonded to a surface hydroxyl oxygen, whereas adsorption of CO(3)2– and Ca(2)UO(2)(CO3)3, which attach to the surface via hydrogen bonding from a surface hydroxyl group to a carbonate oxygen, was calculated to be either only slightly favorable or unfavorable.

The structure and dynamics of hydrated and hydroxylated magnesium oxide nanoparticles

An understanding of the structure of water on metal oxide nanoparticles is important due to its involvement in a number of surface processes, such as in the modification of transport near surfaces and the resulting impact on crystal growth and dissolution. However, as direct experimental measurements probing the metal oxide-water interface of nanoparticles are not easily performed, we use atomistic simulations using experimentally derived potential parameters to determine the structure and dynamics of the interface between magnesium oxide nanoparticles and water. We use a simple strategy to generate mineral nanoparticles, which can be applied to any shape, size, or composition. Molecular dynamics simulations were then used to examine the structure of water around the nanoparticles, and highly ordered layers of water were found at the interface. The structure of water is strongly influenced by the crystal structure and morphology of the mineral and the extent of hydroxylation of the surface. Comparison of the structure and dynamics of water around the nanoparticles with their two-dimensional flat surface counterparts revealed that the size, shape, and surface composition also affects properties such as water residence times and coordination number.

Nanominerals and nanoparticles in feed coal and bottom ash: implications for human health effects

Environmental and human health risk assessments of nanoparticle effects from coal and bottom ash require thorough characterisation of nanoparticles and their aggregates. In this manuscript, we expand the study of human exposure to nanosized particles from coal combustion sources (typically <100 nm in size), characterising the complex micromineralogy of these airborne combustion-derived nanomaterials. Our study focuses on bottom ash generated in the Santa Catarina power station (Brazil) which uses coal enriched in ashes, many potential elements (e.g. Cr and Ni) and pyrite. Transmission electron microscope data reveal nanoscale C deposits juxtaposed with and overgrown by slightly larger aluminosilicate (Al-Si) glassy spheres, oxides, silicates, carbonated, phosphates and sulphates. Iron oxides (mainly hematite and magnetite) are the main bottom ash products of the oxidation of pyrite, sometimes via intermediate pyrrhotite formation. The presence of iron oxide nanocrystals mixed with silicate glass particles emphasises the complexity of coal and bottom ash micromineralogy. Given the potentially bioreactive nature of such transition metal-bearing materials, there is likely to be an increased health risk associated with their inhalation.

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.

Molecular simulations of water and ion diffusion in nanosized mineral fractures

Molecular dynamics simulations were carried out to investigate the effects of confinement and of the presence of the mineral surface on the diffusion of water and electrolyte ions in nanosized mineral fractures. Feldspar was used as a representative mineral because recent studies found that it is an important mineral that hosts contaminants within its intragrain fractures at the U.S. Department of Energy Hanford site. Several properties of the mineral-water interface were varied, such as the fracture width, the ionic strength of the contacting solution, and the surface charge,to provide atomic-level insights into the diffusion of ions and contaminants within intragrain regions. In each case, the self-diffusion coefficient of water and that of the electrolyte ions were computed as a function of distance from the mineral surface. Our calculations reveal a 2.0-2.5 nm interfacial region within which the self-diffusion coefficient of water and that of the electrolyte ions decrease asthe diffusing species approach the surface. As a result of the extent of the interfacial region, water and electrolyte ions are predicted to never reach bulk-like diffusion in fractures narrower than approximately 5 nm. The average diffusion coefficient along the mineral fracture was computed as a function of fracture width and indicated that the surface effects only become negligible for fractures several tens of nanometers wide. The molecular dynamics results improve our conceptual models of ion transport in nanoscale pore regions surrounded by mineral surfaces in porous media.

Structural models of bioactive glasses from molecular dynamics simulations

The bioactive mechanism, by which living tissues attach to and integrate with an artificial implant through stable chemical bonds, is at the core of many current medical applications of biomaterials, as well as of novel promising applications in tissue engineering. Having been employed in these applications for almost 40 years, soda-lime phosphosilicate glasses such as 45S5 represent today the paradigm of bioactive materials. Despite their strategical importance in the field, the relationship between the structure and the activity of a glass composition in a biological environment has not been studied in detail. This fundamental gap negatively affects further progress, for instance, to improve the chemical durability and tailor the biodegradability of these materials for specific applications. This paper reviews recent advances in computer modelling of bioactive glasses based on molecular dynamics simulations, which are starting to unveil key structural features of these materials, thus contributing to improve our fundamental understanding of how bioactive materials work.

Atomistic simulations of calcite nanoparticles and their interaction with water

Molecular dynamics (MD) simulations have been used to study the stability of calcite nanoparticles ranging in size from 18 to 324 f.u., both in vacuo and in the presence of explicit water molecules. In vacuo, the smallest particles become highly disordered during the MD simulation due to rotation and translation of the undercoordinated CO(3) (2-) anions at the edges of the particles. As the nanoparticle size increases, the influence of the fully coordinated bulk ions begins to dominate and long-range order is seen both in the Ca-C pair distribution functions and in the degree of rotational order of the CO(3) (2-) anions. However, when explicit water is added to the system, the molecules in the first hydration layer complete the coordination shell of the surface ions, preserving structural order even in the smallest of the nanoparticles. Close to particle surface, the structure of the water itself shows features similar to those seen close to planar periodic (1014) surfaces, although the molecules are far less tightly bound.

Adsorption and self-assembly of surfactant/supercritical CO2 systems in confined pores: a molecular dynamics simulation

A coarse-grained molecular dynamics simulation has been carried out to study the adsorption and self-organization for a model surfactant/supercritical CO2 system confined in the slit-shape nanopores with amorphous silica-like surfaces. The solid surfaces were designed to be CO2-philic and CO2-phobic, respectively. For the CO2-philic surface, obviously surface adsorption is observed for the surfactant molecules. The various energy profiles were used to monitor the lengthy dynamics process of the adsorption and self-assembly for surfactant micelles or monomers in the confined spaces. The equilibrium properties, including the morphologies and micelle-size distributions of absorbed surfactants, were evaluated based on the equilibrium trajectory data. The interaction between the surfactant and the surface produces an obvious effect on the dynamics rate of surfactant adsorption and aggregation, as well as the final self-assembly equilibrium structures of the adsorbed surfactants. However, for the CO2-phobic surfaces, there are scarcely adsorption layers of surfactant molecules, meaning that the CO2-phobic surface repels the surfactant molecules. It seems to conclude that the CO2 solvent depletion near the interfaces determines the surface repellence to the surfactant molecules. The effect of the CO2-phobic surface confinement on the surfactant micelle structure in the supercritical CO2 has also been discussed. In summary, this study on the microscopic behaviors of surfactant/Sc-CO2 in confined pores will help to shed light on the surfactant self-assembly from the Sc-CO2 fluid phase onto solid surfaces and nanoporous media.