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

Interfacial processes underlying the temperature-dependence of friction and wear of calcite single crystals

HYPOTHESIS

This study posits that thermal effects play a substantial role in influencing interfacial processes on calcite, and consequently impacting its mechanochemical properties.

EXPERIMENTS

This work interrogates the temperature-dependence of friction and wear at nanoscale contacts with calcite single crystals at low air humidity (≤ 3-10 % RH) by AFM.

FINDINGS

Three logarithmic regimes for the velocity-dependence of friction are identified. BelowTc ∼ 70 °C, where friction increases with T, there is a transition from velocity-weakening (W1) to velocity-strengthening friction (S1). AboveTc ∼ 70 °C, where friction decreases with T, a second velocity-strengthening friction regime (S0) precedes velocity-weakening friction (W1). The low humidity is sufficient to induce atomic scale changes of the calcite cleavage plane due to dissolution-reprecipitation, and more so at higher temperature and 10 % RH. Meanwhile, the surface softens above Tc -likely owing to lattice dilation, hydration and amorphization. These interfacial changes influence the wear mechanism, which transitions from pit formation to plowing with increase in temperature. Furthermore, the softening of the surface justifies the appearance of the second velocity-strengthening friction regime (S0). These findings advance our understanding of the influence of temperature on the interfacial and mechanochemical processes involving calcite, with implications in natural processes and industrial manufacturing.

The dynamic nature of natural and fatty acid modified calcite surfaces

Calcium carbonate, particularly in the form of calcite, is an abundant mineral widely used in both human-made products and biological systems. The calcite surface possesses a high surface energy, making it susceptible to the adsorption of organic contaminants. Moreover, the surface is also reactive towards a range of chemicals, including water. Consequently, studying and maintaining a clean and stable calcite surface is only possible under ultrahigh vacuum conditions and for limited amounts of time. When exposed to air or solution, the calcite surface undergoes rapid transformations, demanding a comprehensive understanding of the properties of calcite surfaces in different environments. Similarly, attention must also be directed towards the kinetics of changes, whether induced by fluctuating environments or at constant condition. All these aspects are encompassed in the expression "dynamic nature", and are of crucial importance in the context of the diverse applications of calcite. In many instances, the calcite surface is modified by adsorption of fatty acids to impart a desired nonpolar character. Although the binding between carboxylic acid groups and calcite surfaces is strong, the fatty acid layer used for surface modification undergoes significant alterations when exposed to water vapour and liquid water droplets. Therefore, it is also crucial to understand the dynamic nature of the adsorbed layer. This review article provides a comprehensive overview of the current understanding of both the dynamics of the calcite surface as well as when modified by fatty acid surface treatments.

Effect of the Water Content on the Adsorption of CO2 and CH4 in Calcite Slit Nanopores: Insights from GCMC, MD, and DFT

It is significant to understand the adsorption mechanisms of shale gas (CH4) and CO2 in shale formations to enhance CH4 recovery rates and enable geological CO2 storage. This study provides a comprehensive investigation into the adsorption behaviors of CO2 and CH4 within dry and hydrous calcite nanopores, utilizing a combination of grand canonical Monte Carlo simulations, molecular dynamics simulations, and density functional theory calculations. In dry calcite slits, the calculated results for the adsorption capacity, density profile, and isosteric heat of CO2 and CH4 reveal that CO2 possesses a stronger adsorption affinity, making it preferentially adsorb on the pore surface compared to CH4. In hydrous calcite slits, calculating the adsorption capacity and density profile of CO2 and CH4, the results show that the gas adsorption sites become progressively occupied by H2O molecules, leading to a substantial decrease in the adsorption capacity of CO2 and CH4. Furthermore, by analysis of the adsorption energy and electronic structure, the reason for the reduction of gas adsorption capacity caused by H2O is further revealed. This work has a deep understanding of the adsorption mechanisms of shale gas and CO2 in calcite and can offer valuable theoretical insights for the development of a CO2-enhanced shale gas recovery technology.

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

HYPOTHESIS

The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). We showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers.

EXPERIMENTS

First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances.

FINDINGS

Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

Structural water in amorphous carbonate minerals: ab initio molecular dynamics simulations of X-ray pair distribution experiments

Water is known to play a controlling role in directing mineralization pathways and stabilizing metastable amorphous intermediates in hydrous carbonate mineral MCO₃·nH₂O systems, where M²⁺ is a divalent metal cation. Despite this recognition, the nature of the controls on crystallization is poorly understood, largely owing to the difficulty in characterizing the dynamically disordered structures of amorphous intermediates at the atomic scale. Here, we present a series of atomistic models, derived from ab initio molecular dynamics simulation, across a range of experimentally relevant cations (M = Ca, Mg, Sr) and hydration levels (0 ≤ n ≤ 2). Theoretical simulations of the dependence of the X-ray pair distribution function on the hydration level n show good agreement with available experimental data and thus provide further evidence for a lack of significant nanoscale structure in amorphous carbonates. Upon dehydration, the metal coordination number does not change significantly, but the relative extent of water dissociation increases, indicating that a thermodynamic driving force exists for water dissociation to accompany dehydration. Mg strongly favors monodentate conformation of carbonate ligands and shows a marked preference to exchange monodentate carbonate O for water O upon hydration, whereas Ca and Sr exchange mono- and bidentate carbonate ligands with comparable frequency. Water forms an extensive hydrogen bond network among both water and carbonate groups that exhibits frequent proton transfers for all three cations considered suggesting that proton mobility is likely predominantly due to water dissociation and proton transfer reactions rather than molecular water diffusion.

Obtaining Consistent Free Energies for Ion Binding at Surfaces from Solution: Pathways versus Alchemy for Determining Kink Site Stability

Ion incorporation or removal from a solid at the interface with solution is a fundamental part of crystal growth. Despite this, there have been few quantitative determinations of the thermodynamics for such processes from atomistic molecular dynamics due to the associated technical challenges. In this study, we compute the free energies for ion removal from kink sites at the interface between NaCl and water as an illustrative example. To examine the influence of the free energy technique used, we compare methods that follow an explicit pathway for dissolution with those that focus on the thermodynamics of the initial and final states using metadynamics and free energy perturbation, respectively. While the initial results of the two approaches are found to be completely different, it is demonstrated that the thermodynamics can be reconciled with appropriate corrections for the standard states, thus illustrating the need for caution in interpreting raw free energy curves for ion binding as widely found in the literature. In addition, a new efficient approach is introduced to correct for the system size dependence of kink site energies both due to the periodic interaction of charges in an inhomogeneous dielectric system and due to the dipolar interactions between pairs of kinks along a row. Ultimately, it is shown that with suitable care, both classes of free energy techniques are capable of producing kink site stabilities that are consistent with the solubility of the underlying bulk solid. However, the precise values for individual kink sites exhibit a small systematic offset, which can be ascribed to the contribution of the interfacial potential to the pathway-based results. For the case of NaCl, the free energies of the kink sites relative to a 1 M aqueous solution for Na⁺ and Cl^- are found to be surprisingly different and of opposite sign, despite the ions having very similar hydration free energies.

The effect of a mixture of an ionic liquid and organic solvent on oxygen reduction reaction kinetics

Li-O₂ batteries attract great attention due to their promising theoretical energy density. One of the main obstacles on the way to achieving high energy density and good cyclability is positive electrode passivation by the Li₂O₂ discharge product as well as the presence of parasitic reactions that degrade electrode and electrolyte materials. To overcome these issues new electrolytes are being extensively searched for to ensure the bulk-mediated mechanism of the oxygen reduction reaction and inhibition of parasitic reactions. Different additives to organic solvents can significantly change the properties of electrolytes. This work is devoted to the effect of ionic liquids (ILs), which are proposed as an additive to the solvent due to their excellent solvation properties, high stability, low volatility and flammability. Using molecular dynamics simulations we investigate mixtures of the Pyr₁₄TFSI ionic liquid and dimethoxyethane (DME) with different volume fractions of the IL. Our calculations show that the presence of the ionic liquid in the electrolyte stabilises solvation shells around the ions, both involved in the oxygen reduction and parasitic reactions, slowing down the kinetics of Li⁺ and O₂^- association. This makes the usage of such mixtures promising for electrolyte design for Li-O₂ batteries.

A General Method for Calculating Solid/Liquid Interfacial Free Energies from Atomistic Simulations: Application to CaSO4.xH2O

We present a general method for computing interfacial free energies from atomistic simulations, which is particularly suitable for solid/liquid interfaces. Our method uses an Einstein crystal as a universal reference state and is more flexible than previous approaches. Surfaces with dipoles, complex reconstructions, and partially dissolved species are all easily accommodated within the framework. It may also be extended to calculating the relative free energies of different phases and other types of defect. We have applied our method to interfaces of bassanite and gypsum with water and obtained interfacial free energies of the order of 0.15 J/m2, of which approximately 50 % is due to entropic contributions. Our calculations of the interfacial free energy of NaCl with water obtained a value of 0.13 J/m2 of which only 19 % is from entropic contributions. We have also predicted equilibrium morphologies for bassanite and gypsum that compare well with experiments and previous calculations.

Electrokinetic behavior of artificial and natural calcites: A review of experimental measurements and surface complexation models

The surface charge of calcite in aqueous environments is essential to many industrial and environmental applications. Electrokinetic measurements are usually used to assess the calcite charging behavior and characterize its electrical double layer (EDL). Numerous surface complexation models (SCMs) have been proposed to interpret the effect of different surface interactions on the zeta potential. Because of their versatility, SCMs have also become important tools in reactive transport modeling. The research on enhanced oil recovery within the last decade has led to an increased number of publications reporting both zeta potential measurements and SCMs for calcite. Nonetheless, the measurements are often inconsistent and the reasons for choosing one model over another are unclear. In this work, we review the models proposed for calcite and address their main differences. We first collect a large number of published zeta potential measurements and then we fit a Diffuse Layer, Basic Stern, and Charge-Distribution Multi-Site Complexation models to a selected reliable dataset. For each model, we maintain a similar number of adjustable parameters. After optimizing the parameters of the models, we systematically compare their prediction capabilities against data obtained in monovalent and divalent electrolyte systems containing calcium, magnesium, sulfate, or carbonate. We show that, often, the discrepancies between the models and the experimental data can be explained by different levels of disequilibrium. Nonetheless, assumptions used in the development of the models may significantly reduce their extrapolability to variable chemical conditions. The poor agreement between the models tuned to electrokinetic data with surface charge measurements and dynamic retention from single-phase flowthrough tests show that zeta potential may not be the best type of data to characterize ion binding at the calcite surface. Including the effect of mineral impurities and temperature on the calcite surface speciation and electrokinetic behavior prevail as main challenges for reactive transport modeling.

Morphological characteristics of calcium carbonate crystallization in CO2 pre-cured aerated concrete

Early-stage CO₂ curing technology for alkaline construction materials (such as cement concrete) has gained increasing interest owing to the advantages of material properties improvement and high potential of CO₂ sinking. Less attention, however, has been paid to morphological characteristics of CaCO₃ in carbonated cement concrete. The crystal structure and micromorphology of CaCO₃ in an early-age aerated concrete (AC) cured under CO₂ gas pressures of 0.1, 1, and 2 bar were investigated. The fabricated AC has a high CO₂ sorption capacity (∼35 g CO₂ per 100 g cement in a 100 mm cube). The morphological characteristics of CaCO₃ were statistically analyzed in terms of long-axis length (b), short-axis length (a), and aspect ratio (K = b/a). As CO₂ pressure increases, b is almost unchanged from 0.8–1.8 μm, a decreases from 0.7 to 0.4 μm, and, consequently, K increases from 1.3 to 2.5. The different CaCO₃ crystal morphologies in AC are ascribed to the CO₂ pressure-associated crystal transformation processes: low gas pressure induces a symmetric CaCO₃ growth, while high gas pressure causes a faster calcite growth at the crystal tip ends. The findings would deepen the understanding of CaCO₃ crystal formation under different CO₂ curing pressures for tuning the microstructure of CO₂-cured cement concrete.

The work reports different morphological characteristics of CaCO₃ formed in an early-age aerated concrete (AC) under different CO₂ pressures, uncovering the physicochemical mechanisms of carbonation of cement-based materials affected by CO₂ curing.

Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

Reactive force fields for aqueous and interfacial magnesium carbonate formation

We develop Mg/C/O/H ReaxFF parameter sets for two environments: an aqueous force field for magnesium ions in solution and an interfacial force field for minerals and mineral-water interfaces. Since magnesium is highly ionic, we choose to fix the magnesium charge and model its interaction with C/O/H through Coulomb, Lennard-Jones, and Buckingham potentials. We parameterize the forcefields against several crystal structures, including brucite, magnesite, magnesia, magnesium hydride, and magnesium carbide, as well as Mg²⁺ water binding energies for the aqueous forcefield. Then, we test the forcefield for other magnesium-containing crystals, solvent separated and contact ion-pairs and single-molecule/multilayer water adsorption energies on mineral surfaces. We also apply the forcefield to the forsterite-water and brucite-water interface that contains a bicarbonate ion. We observe that a long-range proton transfer mechanism deprotonates the bicarbonate ion to carbonate at the interface. Free energy calculations show that carbonate can attach to the magnesium surface with an energy barrier of about 0.22 eV, consistent with the free energy required for aqueous Mg-CO₃ ion pairing. Also, the diffusion constant of the hydroxide ions in the water layers formed on the forsterite surface are shown to be anisotropic and heterogeneous. These findings can help explain the experimentally observed fast nucleation and growth of magnesite at low temperature at the mineral-water-CO₂ interface in water-poor conditions.

Adsorptive reduction of water hardness by a highly porous and regenerative geopolymer fabricated from coal fly ash waste with low-temperature calcination

In this research paper, potassium-activated geopolymer cubes (GeoC) fabricated from waste coal fly ash with low-temperature calcination were investigated as a water softening agent. The GeoC reduced water hardness contents by adsorbing calcium (Ca²⁺) and magnesium (Mg²⁺) ions from aqueous solutions. Batch experiments were conducted to investigate the adsorption performance for Ca²⁺ and Mg²⁺, including contact time, initial concentration of cations, and interference with competitive cations. The best performance for water hardness adsorption was found on GeoC-35, fabricated with the highest silicate ratio to hydroxide. The adsorption process reached equilibrium after a contact time of 6 h for Ca²⁺ and 24 h for Mg²⁺. The maximum adsorption capacity for Ca²⁺ and Mg²⁺ was 52.0 and 17.3 mg/g, respectively. Langmuir and pseudo-second-order models fitted the experimental data well, indicating that chemical reactions occurred on a homogeneous surface. The GeoC can also be reused for removing hardness. Furthermore, the increase in potassium and silicon concentration in solution varied directly with removal efficiency, suggesting that the aluminosilicate framework played a role in reducing water hardness via cationic exchange. The presence of competitive cations decreased adsorption ability, albeit it still exhibited an appreciable removal performance.

Structure and Surface Complexation at the Calcite(104)-Water Interface

Calcite is the most stable polymorph of calcium carbonate (CaCO₃) under ambient conditions and is ubiquitous in natural systems. It plays a major role in controlling pH in environmental settings. Electrostatic phenomena at the calcite-water interface and the surface reactivity of calcite in general have important environmental implications. They may strongly impact nutrient and contaminant mobility in soils and other subsurface environments, they control oil recovery from limestone reservoirs, and they may impact the safety of nuclear waste disposal sites. Besides the environmental relevance, the topic is significant for industrial applications and cultural heritage preservation. In this study, the structure of the calcite(104)-water interface is investigated on the basis of a new extensive set of crystal truncation rod data. The results agree with recently reported structures and resolve previous ambiguities with respect to the coordination sphere of surface Ca ions. These structural features are introduced into an electrostatic three-plane surface complexation model, describing ion adsorption and charging at the calcite-water interface. Inner surface potential data for calcite, as measured with a calcite single-crystal electrode, are used as constraints for the model in addition to zeta potential data. Ion adsorption parameters are compared with molecular dynamics simulations. All model parameters, including protonation constants, ion-binding parameters, and Helmholtz capacitances, are within physically and chemically plausible ranges. A PhreeqC version of the model is presented, which we hope will foster application of the model in environmental studies.

Effects of surface contamination on the interfacial properties of CO2/water/calcite systems

Understanding the wetting properties of reservoir rocks can be of great benefit for advanced applications such as the effective trapping and geological storage of CO₂. Despite their importance, not all mechanisms responsible for wetting mineral surfaces in subsurface environments are well understood. Factors such as temperature, pressure and salinity are often studied, achieving results with little unanimity; other possible factors are left somewhat unexplored. One such factor is the effect of contamination. In the present study, the effects of adding a non-aqueous organic contaminant, ethanol, on the CO₂-water interfacial tension (IFT) and the CO₂/water/calcite contact angle were investigated using molecular dynamics simulations. Within the conditions studied, relatively small amounts of ethanol cause a significant decrease in the CO₂-water IFTs, as well as a pronounced increase in the water-calcite-CO₂ three phase contact angle. The latter result is due to the decrease of the IFT between CO₂ and water and the strong adsorption of ethanol on the solid substrate. These findings could be helpful for explaining how impurities can affect experimental data and could lead to effective carbon sequestration strategies.

Molecular Dynamics Simulation of Ion Adsorption and Ligand Exchange on an Orthoclase Surface

Orthoclase (K-feldspar) is one of the natural inorganic materials, which shows remarkable potential toward removing heavy metal ions from aqueous solutions. Understanding the interactions of the orthoclase and metal ions is important in the treatment of saline wastewater. In this paper, molecular dynamics simulations were used to prove the adsorption of different ions onto orthoclase. The adsorption isotherms show that orthoclase has remarkable efficiency in the removal of cations at low ion concentrations. Aluminol groups are the preferential adsorption sites of cations due to higher negative charges. The adsorption types and adsorption sites are influenced by the valence, radius, and hydration stability of ions. Monovalent cations can be adsorbed in the cavities, whereas divalent cations cannot. The hydrated cation may form an outer-sphere complex or an inner-sphere complex in association with the loss of hydration water. Na+, K+, and Ca2+ ions mainly undergo inner-sphere adsorption and Mg2+ ions prefer outer-sphere adsorption. On the basis of simulation results, the mechanism of ion removal in the presence of orthoclase is demonstrated at a molecular level.

Atomistic insight into salinity dependent preferential binding of polar aromatics to calcite/brine interface: implications to low salinity waterflooding

This paper resolve the salinity-dependent interactions of polar components of crude oil at calcite-brine interface in atomic resolution. Molecular dynamics simulations carried out on the present study showed that ordered water monolayers develop immediate to a calcite substrate in contact with a saline solution. Carboxylic compounds, herein represented by benzoic acid (BA), penetrate into those hydration layers and directly linking to the calcite surface. Through a mechanism termed screening effect, development of hydrogen bonding between –COOH functional groups of BA and carbonate groups is inhibited by formation of a positively-charged Na⁺ layer over CaCO₃ surface. Contrary to the common perception, a sodium-depleted solution potentially intensifies surface adsorption of polar hydrocarbons onto carbonate substrates; thus, shifting wetting characteristic to hydrophobic condition. In the context of enhanced oil recovery, an ion-engineered waterflooding would be more effective than injecting a solely diluted saltwater.

Nanoscale insight into the relation between pressure solution of calcite and interfacial friction

Pressure solution of carbonate-based rocks participates in many geophysical and geochemical processes, but fundamental knowledge of the interfacial processes is still lacking. By concurrently pressing and sliding two single calcite crystals past each other, the pressure solution rate and the friction force between the crystals were concurrently measured in calcium-carbonate saturated water with an extended surface forces apparatus. These studies reveal that both a decrease and an increase in frictional strength can originate from the pressure-solution of calcite single crystals. By conducting nanoscale force measurements with an atomic force microscope, ion specific effects were unveiled at the level of a single asperity. Pressure solution is promoted when the interfacial water layers of calcite remain undisturbed under stress (e.g. with Ca(II)) and the dissolved ions and water lubricate the interface - a phenomenon called pressure-solution facilitated slip. The mechanically induced disruption of the hydration layers of the calcite surface (e.g. with Mg(II) and low Ni(II) concentration) correlates with the more fluid-like and lubricious behavior of the confined fluid in the absence of pressure solution. Charge neutralization of the calcite surface leads to an abrupt change of calcite's hydration layers, which promotes pressure-solution facilitated slip. This work advances the fundamental understanding of physicochemical interactions occurring at confined surfaces of stressed calcite.

First Steps towards Understanding the Non-Linear Impact of Mg on Calcite Solubility: A Molecular Dynamics Study

Magnesium (Mg2+) is one of the most common impurities in calcite and is known to have a non-linear impact on the solubility of magnesian calcites. Using molecular dynamics (MD), we observed that Mg2+ impacts overall surface energies, local free energy profiles, interfacial water density, structure and dynamics and, at higher concentrations, it also causes crystal surface deformation. Low Mg concentrations did not alter the overall crystal structure, but stabilised Ca2+ locally and tended to increase the etch pit nucleation energy. As a result, Ca-extraction energies over a wide range of 39 kJ/mol were observed. Calcite surfaces with an island were less stable compared to flat surfaces, and the incorporation of Mg2+ destabilised the island surface further, increasing the surface energy and the calcium extraction energies. In general, Ca2+ is less stable in islands of high Mg2+ concentrations. The local variation in free energies depends on the amount and distance to nearest Mg in addition to local disruption of interfacial water and the flexibility of surface carbonate ions to rotate. The result is a complex interplay of these characteristics that cause variability in local dissolution energies. Taken together, these results illustrate molecular scale processes behind the non-linear impact of Mg2+ concentration on the solubility of magnesium-bearing calcites.

Enthalpy of Uranium Adsorption onto Hematite

The influence of temperature on the adsorption of metal ions at the solid-water interface is often overlooked, despite the important role that adsorption plays in metal-ion fate and transport in the natural environment where temperatures vary widely. Herein, we examine the temperature-dependent adsorption of uranium, a widespread radioactive contaminant, onto the ubiquitous iron oxide, hematite. The multitemperature batch adsorption data and surface complexation models indicate that the adsorption of uranium, as the hexavalent uranyl (UO₂²⁺) ion, increases significantly with increasing temperature, with an adsorption enthalpy (ΔH_ads) of +71 kJ mol^-1. We suggest that this endothermic, entropically driven adsorption behavior is linked to reorganization of the uranyl-ion hydration and interfacial water structures upon U^VI adsorption at the hematite surface. Overall, this work provides fundamental insight into the thermodynamics driving metal-ion adsorption reactions and provides the specific enthalpy value necessary for improved predictive geochemical modeling of U^VI adsorption in the environment.

Carbon dioxide and water incorporation mechanisms in SrFeO3-δ phases: a computational study

With a higher propensity for low temperature synthesis routes along with a move toward lower solid oxide fuel cell operating temperatures, water and carbon dioxide incorporation in strontium ferrite is of importance. Despite this, the mechanisms are not well understood. In this work, classical-potential-based computational techniques are used to determine the favourability of water and CO₂ incorporation mechanisms in both SrFeO_3-δ and SrFeO_2.5. Our studies suggest that intrinsic Frenkel and Schottky type defects are unlikely to form, but that water and carbon dioxide incorporation are favourable in both phases. Water incorporation is likely for both the cubic and brownmillerite phases, with hydroxyl ions preferring to sit on octahedral oxygen sites in both structures, causing slight tilting of the shared octahedra. Interstitial hydroxyl ions are only likely for the brownmillerite phase, where the hydroxyl ions are most stable between adjacent FeO₄ tetrahedral chains. Carbon dioxide incorporation via carbonate defects is most favourable when a carbonate molecule exists on an iron site, preferring the iron site with lower oxygen coordination. This involves formation of multiple oxygen vacancies surrounding the iron site, and thus we conclude that carbonate can trap oxygen vacancies.

CrystalGrower: a generic computer program for Monte Carlo modelling of crystal growth

A Monte Carlo crystal growth simulation tool, CrystalGrower, is described which is able to simultaneously model both the crystal habit and nanoscopic surface topography of any crystal structure under conditions of variable supersaturation or at equilibrium. This tool has been developed in order to permit the rapid simulation of crystal surface maps generated by scanning probe microscopies in combination with overall crystal habit. As the simulation is based upon a coarse graining at the nanoscopic level features such as crystal rounding at low supersaturation or undersaturation conditions are also faithfully reproduced. CrystalGrower permits the incorporation of screw dislocations with arbitrary Burgers vectors and also the investigation of internal point defects in crystals. The effect of growth modifiers can be addressed by selective poisoning of specific growth sites. The tool is designed for those interested in understanding and controlling the outcome of crystal growth through a deeper comprehension of the key controlling experimental parameters.

Generic in silico methodology – CrystalGrower – for simulating crystal habit and nanoscopic surface topology to determine crystallisation free energies.

Hydrogen-Bond Structure and Low-Frequency Dynamics of Electrolyte Solutions: Hydration Numbers from ab Initio Water Reorientation Dynamics and Dielectric Relaxation Spectroscopy

We present an atomistic simulation scheme for the determination of the hydration number (h) of aqueous electrolyte solutions based on the calculation of the water dipole reorientation dynamics. In this methodology, the time evolution of an aqueous electrolyte solution generated from ab initio molecular dynamics simulations is used to compute the reorientation time of different water subpopulations. The value of h is determined by considering whether the reorientation time of the water subpopulations is retarded with respect to bulk-like behavior. The application of this computational protocol to magnesium chloride (MgCl2 ) solutions at different concentrations (0.6-2.8 mol kg-1 ) gives h values in excellent agreement with experimental hydration numbers obtained using GHz-to-THz dielectric relaxation spectroscopy. This methodology is attractive because it is based on a well-defined criterion for the definition of hydration number and provides a link with the molecular-level processes responsible for affecting bulk solution behavior. Analysis of the ab initio molecular dynamics trajectories using radial distribution functions, hydrogen bonding statistics, vibrational density of states, water-water hydrogen bonding lifetimes, and water dipole reorientation reveals that MgCl2 has a considerable influence on the hydrogen bond network compared with bulk water. These effects have been assigned to the specific strong Mg-water interaction rather than the Cl-water interaction.

Entropy-Enthalpy Compensation in Peptide Adsorption on Solid Surfaces: Dependence on Surface Hydration

Although protein adsorption at the solid-water interface is of immense importance, understanding the crucial role of the water phase in mediating protein-surface interactions is lacking, particularly due to the lack of fundamental thermodynamic data. Herein, we have performed complicated free energy calculations and successfully extracted the entropy and enthalpy changes of molecular adsorption on solids. Using the gold and graphene as the surface models with distinct affinities to the water phase, we successfully unravel the sharply opposite manners of entropy-enthalpy compensation in driving water and tripeptide adsorptions on two surfaces. Though the thermodynamic features of water adsorption on surface are enthalpically dominated based on the positions of free energy barriers and minima, the favorable entropy term significantly decreases the free energy barrier and further stabilizes the adsorbate at the adsorption site on the graphene surface. For the peptide, the shape of the adsorption free energy profile is jointly determined by the enthalpy and entropy changes, which, however, alternatively act the driving force to promote the peptide adsorption on the Au surface and graphene surface. The distinct structural and dynamic properties of solid-liquid interfaces account for the special role of the interfacial water phase in regulating the competitive relationship between the entropy and enthalpy variations.

Ion Distribution and Hydration Structure at Solid-Liquid Interface between NaCl Crystal and Its Solution

The interface structure between NaCl crystal and its solution has been investigated at the saturated concentration of 298 K by molecular dynamics simulations. We have found that there are many fine structures at this complex interface. Near the surface of crystal, most of Na⁺ only coordinate with water molecules, while almost all Cl^- coordinate with Na⁺ in addition to water molecules. An ion coordinating with more water molecules is farther away from the epitaxial position of lattice. As approaching to the interface, the first hydration shell of ions has the tendency of being ordered, while the orientation of dipole of water molecules in the first hydration shell becomes more disordered than that in the solution. Generally, the first hydration shell of Na⁺ is less affected by nearest Cl^-, whereas the first hydration shell of Cl^- is significantly affected by nearest Na⁺.

Reconsidering Calcium Dehydration as the Rate-Determining Step in Calcium Mineral Growth

The dehydration of cations is generally accepted as the rate-limiting step in many processes. Molecular dynamics (MD) can be used to investigate the dynamics of water molecules around cations, and two different methods exist to obtain trajectory-based water dehydration frequencies. Here, these two different post-processing methods (direct method versus survival function) have been implemented to obtain calcium dehydration frequencies from a series of trajectories obtained using a range of accepted force fields. None of the method combinations reproduced the commonly accepted experimental water exchange frequency of 10-8.2 s-1. Instead, our results suggest much faster water dynamics, comparable with more accurate ab initio MD simulations and with experimental values obtained using neutron scattering techniques. We obtained the best agreement using the survival function method to characterize the water dynamics, and we show that different method combinations significantly affect the outcome. Our work strongly suggests that the fast water exchange kinetics around the calcium ions is not rate-limiting for reactions involving dissolved/solvated calcium. Our results further suggest that, for alkali and most of the earth alkali metals, mechanistic rate laws for growth, dissolution, and adsorption, which are based on the principle of rate-limiting cation dehydration, need careful reconsideration.

Oxalic Acid Adsorption on Rutile: Molecular Dynamics and ab Initio Calculations

Detailed analysis of the adsorption of oxalic acid ions, that is, oxalate and hydrogenoxalate, on the rutile (110) surface was carried out using molecular dynamics augmented by free energy calculations and supported by ab initio calculations. The predicted adsorption on perfect nonhydroxylated and hydroxylated surfaces with surface charge density from neutral to +0.208 C/m² corresponding to pH values of about 6 and 3.7, respectively, agrees with experimental adsorption data and charge-distribution multisite ion complexation model predictions obtained using the most favorable surface complexes identified in our simulations. We found that outer-sphere complexes are the most favorable, owing to strong hydrogen binding of oxalic acid ions with surface hydroxyls and physisorbed water. The monodentate complex, the most stable among inner-sphere complexes, was about 15 kJ/mol higher in energy, but separated by a large energy barrier. Other inner-sphere complexes, including some previously suggested in the literature as likely adsorption structures such as bidentate and chelate complexes, were found to be unstable both by classical and by ab initio modeling. Both the surfaces and (hydrogen)oxalate ions were modeled using charges scaled to 75% of the nominal values in accord with the electronic continuum theory and our earlier parameterization of (hydrogen)oxalate ions, which showed that nominal charges exaggerate ion-water interactions.

Calcium Phosphate Deposition on Planar and Stepped (101) Surfaces of Anatase TiO2: Introducing an Interatomic Potential for the TiO2/Ca-PO4/Water Interface

Titanium is commonly employed in orthopaedic and dental surgery, owing to its good mechanical properties. The titanium metal is usually passivated by a thin layer of its oxide, and in order to promote its integration with the biological tissue, it is covered by a bioactive material such as calcium phosphate (CaP). Here, we have investigated the deposition of calcium and phosphate species on the anatase phase of titanium dioxide (TiO₂) using interatomic potential-based molecular dynamics simulations. We have combined different force fields developed for CaP, TiO₂, and water, benchmarking the results against density functional theory calculations. On the basis of our study, we consider that the new parameters can be used successfully to study the nucleation of CaP on realistic anatase and rutile TiO₂ nanoparticles, including surface defects.

Ion dissolution mechanism and kinetics at kink sites on NaCl surfaces

Desolvation barriers are present for solute-solvent exchange events, such as ligand binding to an enzyme active site, during protein folding, and at battery electrodes. For solution-grown crystals, desolvation at kink sites can be the rate-limiting step for growth. However, desolvation and the associated kinetic barriers are poorly understood. In this work, we use rare-event simulation techniques to investigate attachment/detachment events at kink sites of a NaCl crystal in water. We elucidate the desolvation mechanism and present an optimized reaction coordinate, which involves one solute collective variable and one solvent collective variable. The attachment/detachment pathways for Na⁺ and Cl^- are qualitatively similar, with quantitative differences that we attribute to different ion sizes and solvent coordination. The attachment barriers primarily result from kink site desolvation, while detachment barriers largely result from breaking ion-crystal bonds. We compute ion detachment rates from kink sites and compare with results from an independent study. We anticipate that the reaction coordinate and desolvation mechanism identified in this work may be applicable to other alkali halides.

The effect of solvation and temperature on the adsorption of small organic molecules on calcite

The effect of multicomponent mixture on small organic molecule adsorption on calcite at nonzero temperature was investigated.

Bonding, structural and thermodynamic analysis of dissociative adsorption of H3O+ ion onto calcite (101⁻4) surface: CPMD and DFT calculations

We used density functional theory (DFT) and Car-Parrinello molecular dynamics (CPMD) simulation to investigate the adsorption and bond formation of hydronium ion (H₃O⁺) onto a [Formula: see text] calcite surface. For surface coverage of 25% to 100%, the nature of H₃O⁺ interaction was explored through electron density and energetics in the context of bond critical points. The adsorbate-adsorbent structure was studied by simulation of pair correlation function. The results revealed that dissociation into water molecule(s) and proton(s) complements H₃O⁺ ion(s) adsorbtion. The H₂O molecule adsorbs onto the surface via its O atom, and interacts with surface calcium in a closed-shell mode; the H⁺ ion makes a covalent bond to the surface oxygen while maintaining H-bonding with water. Adsorption energies were diminished by 70-90 kJ mol^-1 when O_bridge-bonded H⁺ ions transferred to the O_terminal manually. While dissociative adsorption of H₃O⁺ ions is spontaneous at all surface coverages tested, the free energy was lowest at 75% coverage. Also, protonation of a completely pre-hydrated calcite surface leads to stronger interaction of water molecules with the surface. This unique outlook on hydrating calcite provides specific insights into biomineralization of this mineral, and helps depict further pH consequences in the field of biomaterial adsorption. Graphical abstract Dissociative adsorption of hydronium ion onto the surface of calcite.

Ambient atomic resolution atomic force microscopy with qPlus sensors: Part 1

Atomic force microscopy (AFM) is an enormous tool to observe nature in highest resolution and understand fundamental processes like friction and tribology on the nanoscale. Atomic resolution in highest quality was possible only in well-controlled environments like ultrahigh vacuum (UHV) or controlled buffer environments (liquid conditions) and more specified for long-term high-resolution analysis at low temperatures (∼4 K) in UHV where drift is nearly completely absent. Atomic resolution in these environments is possible and is widely used. However, in uncontrolled environments like air, with all its pollutants and aerosols, unspecified thin liquid films as thin as a single molecular water-layer of 200 pm or thicker condensation films with thicknesses up to hundred nanometer, have been a problem for highest resolution since the invention of the AFM. The goal of true atomic resolution on hydrophilic as well as hydrophobic samples was reached recently. In this manuscript we want to review the concept of ambient AFM with atomic resolution. The reader will be introduced to the phenomenology in ambient conditions and the problems will be explained and analyzed while a method for scan parameter optimization will be explained. Recently developed concepts and techniques how to reach atomic resolution in air and ultra-thin liquid films will be shown and explained in detail, using several examples. Microsc. Res. Tech. 80:50-65, 2017. © 2016 Wiley Periodicals, Inc.

Detection of Posner's clusters during calcium phosphate nucleation: a molecular dynamics study

Molecular dynamics simulations of calcium and phosphate ions in water show that Posner-like clusters originate during the aggregation process.

Molecular insight into the nanoconfined calcite-solution interface

Little is known about the influence of nanoconfinement on calcium carbonate mineralization. Here, colloidal probe atomic force microscopy is used to confine the calcite-solution interface with a silica microsphere and to measure Derjaguin-Landau-Verwey-Overbeek (DLVO) and non-DLVO forces as a function of the calcium concentration, also after charge reversal of both surfaces occurs. Through the statistical analysis of the oscillatory component of a strong hydration force, the subnanometer interfacial structure of the confined atomically flat calcite is resolved in aqueous solution. By applying a mechanical work, both water and hydrated counterions are squeezed out from the nanoconfined solution, leaving the calcite surface more negatively charged than the analogous unconfined surfaces. Layer size and applied work allow a distinction between the hydration states of the counterions in the Stern layer; we propose counterions to be inner- and outer-sphere calcium ions, with a population of inner-sphere calcium ions larger than on unconfined calcite surfaces. It is also shown that the composition of the nanoconfined solution can be tuned by varying calcium concentration. This is a fundamental study of DLVO and hydration forces, and of their connection, on atomically flat calcite. More broadly, our work scrutinizes the greatly unexplored relation between surface science and confined mineralization, with implications for diverse areas of inquiry, such as nanoconfined biomineralization, CO₂ sequestration in porous aquifers, and pressure solution and crystallization in confined hydrosystems.

Experimental-Theoretical Approach to the Adsorption Mechanisms for Anionic, Cationic, and Zwitterionic Surfactants at the Calcite-Water Interface

The adsorption of surfactants (DTAB, SDS, and CAPB) at the calcite-water interface was studied through surface zeta potential measurements and multiscale molecular dynamics. The ground-state polarization of surfactants proved to be a key factor for the observed behavior; correlation was found between adsorption and the hard or soft charge distribution of the amphiphile. SDS exhibits a steep aggregation profile, reaching saturation and showing classic ionic-surfactant behavior. In contrast, DTAB and CAPB featured diversified adsorption profiles, suggesting interplay between supramolecular aggregation and desorption from the solid surface and alleviating charge buildup at the carbonate surface when bulk concentration approaches CMC. This manifests as an adsorption profile with a fast initial step, followed by a metastable plateau and finalizing with a sharp decrease and stabilization of surface charge. Suggesting this competition of equilibria, elicited at the CaCO3 surface, this study provides atomistic insight into the adsorption mechanism for ionic surfactants on calcite, which is in accordance with experimental evidence and which is a relevant criterion for developing enhanced oil recovery processes.

Effects of carbon surface topography on the electrode/electrolyte interface structure and relevance to Li–air batteries

Carbon surface topography influences the solvent structure at the interface, concentration distribution of reactants (Li⁺, O₂), and their absorption kinetics.

Precise determination of water exchanges on a mineral surface

Solvent exchanges on solid surfaces and dissolved ions are a fundamental property important for understanding chemical reactions, but the rates of fast exchanges are poorly constrained.

Mechanistic insight into the displacement of CH4 by CO2 in calcite slit nanopores: the effect of competitive adsorption

Displacement of residual CH₄ by CO₂ in calcite slit nanopores owing to the competitive adsorption of CO₂ over CH₄.

Investigation of the Role of Polysaccharide in the Dolomite Growth at Low Temperature by Using Atomistic Simulations

Dehydration of water from surface Mg(2+) is most likely the rate-limiting step in the dolomite growth at low temperature. Here, we investigate the role of polysaccharide in this step using classical molecular dynamics (MD) calculations. Free energy (potential of mean force, PMF) calculations have been performed for water molecules leaving the first two hydration layers above the dolomite (104) surface under the following three conditions: without catalyst, with monosaccharide (mannose), and with oligosaccharide (three units of mannose). MD simulations reveal that there is no obvious effect of monosaccharide in lowering the dehydration barrier for surface Mg(2+). However, we found that there are metastable configurations of oligosaccharide, which can decrease the dehydration barrier of surface Mg(2+) by about 0.7-1.1 kcal/mol. In these configurations, the molecule lies relatively flat on the surface and forms a bridge shape. The hydrophobic space near the surface created by the nonpolar -CH groups of the oligosaccharide in the bridge conformation is the reason for the observed reduction of dehydration barrier.