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.

Density functional theory based molecular dynamics study of solution composition effects on the solvation shell of metal ions

We present an ab initio molecular dynamics study of the alkali metal ions Li+, Na+, K+ and Cs+, and of the alkaline earth metal ions Mg2+ and Ca2+ in both pure water and electrolyte solutions containing the counterions Cl- and SO42-. Simulations were conducted using different density functional theory methods (PBE, BLYP and revPBE), with and without the inclusion of dispersion interactions (-D3). Analysis of the ion-water structure and interaction strength, water exchange between the first and second hydration shell, and hydrogen bond network and low-frequency reorientation dynamics around the metal ions have been used to characterise the influence of solution composition on the ionic solvation shell. Counterions affect the properties of the hydration shell not only when they are directly coordinated to the metal ion, but also when they are at the second coordination shell. Chloride ions reduce the sodium hydration shell and expand the calcium hydration shell by stabilizing under-coordinated hydrated Na(H2O)5+ complexes and over-coordinated Ca(H2O)72+. The same behaviour is observed in CaSO4(aq), where Ca2+ and SO42- form almost exclusively solvent-shared ion pairs. Water exchange between the first and second hydration shell around Ca2+ in CaSO4(aq) is drastically decelerated compared with the simulations of the hydrated metal ion (single Ca2+, no counterions). Velocity autocorrelation function analysis, used to probe the strength of the local ion-water interaction, shows a smoother decay of Mg2+ in MgCl2(aq), which is a clear indication of a looser inter-hexahedral vibration in the presence of chloride ions located in the second coordination shell of Mg2+. The hydrogen bond statistics and orientational dynamics in the ionic solvation shell show that the influence on the water-water network cannot only be ascribed to the specific cation-water interaction, but also to the subtle interplay between the level of hydration of the ions, and the interactions between ions, especially those of opposite charge. As many reactive processes involving solvated metal ions occur in environments that are far from pure water but rich in ions, this computational study shows how the solution composition can result in significant differences in behaviour and function of the ionic solvation shell.

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.

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.