Mariëtte
Wolthers
*ab,
Devis
Di Tommaso
a,
Zhimei
Du
a and
Nora H.
de Leeuw
ab
aDepartment of Chemistry, University College London, 20 Gordon Street, London WC1 H0AJ, U.K. E-mail: m.wolthers@ucl.ac.uk
bDepartment of Earth Sciences-Geochemistry, Utrecht University, P.O. Box 80021, 3508 TA Utrecht, the Netherlands
First published on 30th May 2013
It is generally accepted that cation dehydration is the rate-limiting step to crystal growth from aqueous solution. Here we employ classical molecular dynamics simulations to show that the water exchange frequency at structurally distinct calcium sites in the calcite surface varies by about two orders of magnitude. The decrease in water exchange frequency with progressive embedding of surface calcium ions is thought to be rate limiting to subsequent attachment of carbonate ions during calcite growth. Therefore, a process-based calcite growth kinetics model, reparameterized using the water exchange frequencies computed from molecular dynamics simulations, is used to illustrate the impact of these variations on kink-formation rate, step velocities and bulk growth rate. The calculated frequencies of kink formation show a strong variation with surface structures, which can be amplified depending on the saturation state and calcium to carbonate ratio of the solution. Modelled and measured step velocities and bulk growth rates are generally in agreement, showing that variations in calcite growth rates and step velocities observed experimentally might be at least partially induced by surface topography.
Calcite is the most stable long-term sink for carbon,11 and as such, precipitation of calcite directly from solution and through mineral carbonation is currently one of the most viable routes for carbon sequestration.12 The key issue in efficient and directed calcite precipitation and in determining the long-term stability of calcite-trapped CO2 is a fundamental understanding of the processes controlling crystal growth13–15 and dissolution.16
When a calcite crystal is cleaved, growing and/or dissolving, the dominant face exposed is the (104) surface,15,16 usually featuring acute and obtuse step edges around growth spirals and etch pits17 (Fig. 1). Direct observations of this surface have demonstrated that these structurally different step edges vary in spreading velocity during dissolution and growth.13,16,18,19 Many simulation studies have focused on the interaction of water with the (10
4) surface,20–25 but with the notable exception of Cooke and Elliott,24 no water exchange frequencies of surface sites have thus far been reported. Cooke and Elliott reported average residence times for water molecules in the hydration shell of surface calcium in the (10
4) surface, and in the step edges and corners of nano-particles of different sizes. They observed that for all of their calcite nano-particles, except the smallest one, the calculated water residence times varied significantly between corners, edges and faces.
![]() | ||
Fig. 1 Illustration of the various surface sites (a) top view of a growth island and (b) side view of an etch pit. No water molecules are shown for visibility and one pit edge was cut away as well. |
Simulations of the approach and catchment of ions during mineral growth have shown the importance of desolvation of both the surface and the approaching ion.5,22 Furthermore, while desolvation of an approaching calcium ion can be assisted by additives in solution such as small carboxylated molecules26 or phosphonate molecules,27 the structure of water at CaCO3 cluster–water interfaces has been shown to be the key in determining the CaCO3 nucleation pathway.14 However, no study has investigated the kinetics of progressive desolvation of surface calcium ions during calcite growth in relation to the surface structure. In a recent calcite growth kinetics model,7 a different approach to describing the calcite growth mechanism was chosen, partly following the general approach of Zhang and Nancollas6 for ionic crystal growth. It was assumed that the attachment frequency of calcium is controlled by its dehydration or water exchange frequency (∼10+8 s−1).3,22 The attachment frequency of carbonate and bicarbonate to the calcite surface was thought to be limited by the dehydration frequency of surface calcium sites, i.e. the progressive embedding of calcium ions, in agreement with previous experimental considerations.8 The resulting process-based growth rate equation7 agrees with measured calcite growth-step velocities and bulk growth rates over broad ranges of ionic strength, pH, solution stoichiometry and degree of supersaturation, since it takes into account surface speciation as well as solution composition changes. However, in calculating step velocities, it does not discriminate between acute and obtuse steps.
Given the crucial role played by kink sites and step edge structures of spirals and etch pits during calcite growth and dissolution, and the rate limitation imposed by water exchange events between solution and surface calcium, we have determined the water exchange frequency for a range of such sites at the calcite surface, differentiating between structurally distinct sites (face versus acute and obtuse edge and corner sites, Fig. 1). The results were used to refine a process-based calcite growth kinetics model.7 The refined growth model includes distinct step velocities for acute versus obtuse growth steps and is able to reproduce experimentally observed variability in bulk growth rates in relation to topographical variations.
The computer code used for the MD simulations was DL_POLY 2.20.34 In the DL_POLY code, the integration algorithms are based around the Verlet leap-frog scheme35 and we have used the Nosé–Hoover algorithm36 for the thermostat. The Nosé–Hoover parameters were set at 0.5 ps for both the thermostat and barostat relaxation times.
We have simulated a repeating calcite slab, containing 840 CaCO3 units, with a growth island of 16 CaCO3 units on one side and an etch pit of the same size in the other side of the slab (Fig. 1). The 22 Å gap between the repeating slabs was filled with 2048 water molecules, resulting in a simulation cell containing 14912 species including shells. The simulation cell was equilibrated for 140 ps NPT (constant number of particles, pressure and temperature), at P = 1 atm and 300 K, leading to a water density of ρ = 1.21 g cm−3.22 A similar procedure was carried out at 320 K and 340 K after which the statistics were collected for approximately three nanoseconds at 340 K and one nanosecond at the other temperatures. The simulation cell was tested successfully for absence of gap size effects and strain.
We have used the parameters for the short-range interactions in calcite derived empirically by Pavese et al.,30 who reported very good agreement between their simulated and experimental thermal dependence of structural and elastic properties of calcite. It has been shown in a previous study25 of the surface structures and stabilities of three calcium carbonate polymorphs, namely calcite, aragonite and vaterite, that the potential model derived by Pavese et al.30 for calcite is directly transferable to different calcium carbonate phases, accurately reproducing the experimental morphologies of all three polymorphs.
The potential parameters used for the intra- and inter-molecular water interactions are those initially described in a paper of MD simulations on MgO surfaces.37 For the interactions between water molecules and calcite surfaces, we have used the potential parameters previously fitted to calcite25 and successfully used in MD simulations of water adsorption at point defects and crystal dissolution from calcite steps.38 These potential parameters reproduce the experimental heat of formation of calcite from its aqueous ions to an acceptable degree of accuracy (within 20 kJ mol−1), even though the parameters were not fitted to this process. Moreover, de Leeuw and Parker25 have verified these potential parameters by simulating the structure of ikaite, a calcium carbonate hexahydrate, and found good agreement between calculated and experimental structural data. Calculating the change in enthalpy at 298 K for ikaite, calcite, and water, the change in interaction energy for the dissociation of ikaite per water molecule is 47 kJ mol−1. This compares to experimental values of 47–50 kJ mol−1,39 and suggests that the energies of interactions between calcium, carbonate and water are accurate and realistic. Recent evaluation of the water potential parameters used has shown that previously reported non-physical behaviour22,40 is an artefact of the size of the simulation cell.28 The cell dimensions used here are large enough to prevent such non-physical behaviour. Moreover, assessment of our potential with respect to the calcium–oxygen distances of small calcium ion–water clusters showed excellent agreement with density functional theory calculations.28 For an overview of all potential parameters and detailed re-evaluation of the force field see Wolthers et al.28
Our calcite (104) surfaces feature an etch pit and a growth island (Fig. 1), terminated by steps that have either an acute or obtuse angle to the surface. We refer to sites positioned in step edges as either acute or obtuse edge sites. Calcium corners are either acute or obtuse in our cell, while CO3 corner sites terminate one acute and one obtuse edge. Calcium and carbonate ions embedded within the atomically flat faces surrounding these low-coordinated surface features are referred to here as face sites.
For carbonate oxygen, Oc, an additional label indicates whether the oxygen points out of the surface towards the water molecules (‘type A’), if it is approximately level with carbon in the crystal truncation plane (‘type X’), or pointing in towards the bulk crystal (‘type B’). During equilibration and, to a lesser extent, over the course of the production runs, some carbonate surface groups rotated, for example turning type B oxygen into an A or X position. This rotation was observed to be most significant for growth island corners, and decreases for edge sites and face sites, with face sites showing only minor rotation during equilibration and none during production. This rotation caused counter-intuitive similar average water exchange frequencies observed for the various types of carbonate oxygen atoms.
![]() |
t sim (ps) |
|
logkwb (s−1) | logknormwc (s−1) | MRT (ps) |
---|---|---|---|---|---|
a Defined by the first minimum of the Ca–OW radial distribution functions, which lies at 3.5 Å for all calcium sites (Fig. S1a, ESI). b Water exchange frequency, i.e. the number of exchange events per second, average and standard deviation in ∼3 ns production. c Normalized to experimental value for calcium in water.45 d This study, shell model potential. e DFT (Car–Parrinello) MD simulations.43 f Hybrid ab initio (Hartree–Fock) quantum mechanical/molecular mechanical simulations.41 g Classical MD simulation at 300 K conducted in this study using a calcium–oxygen water Lennard-Jones potential parameterized by Dang and Smith42 and the SPC/E water model. h Classical MD simulation at 300 K conducted in this study using a calcium–oxygen water Lennard-Jones potential parameterized by Raiteri et al.21 and the SPC/E water model. i Classical MD at 300 K.24 j Experimental value.45 k Using Impey methods and 2 ps duration minimum for exchange events. | |||||
Flat face | 3071 | 78 | 10.41 ± 0.08 | 7.06 | 42.3 |
Pit corner | 3071 | 283 | 10.93 ± 0.16 | 7.43 | 11.0 |
Obtuse pit edge | 3071 | 942 | 11.48 ± 0.05 | 7.80 | 6.1 |
Acute pit edge | 3071 | 996 | 11.51 ± 0.06 | 7.82 | 6.2 |
Obtuse island edge | 3071 | 676 | 11.34 ± 0.07 | 7.70 | 8.5 |
Acute island edge | 3071 | 1338 | 11.65 ± 0.10 | 7.92 | 4.6 |
Obtuse island corner | 3071 | 1048 | 11.51 ± 0.15 | 7.82 | 9.0 |
Acute island corner | 3071 | 1791 | 11.76 ± 0.13 | 7.99 | 5.2 |
Ca2+ in 2027H2Od | 2988 | 3499 | 12.07 ± 0.02 | 8.20j | 7.0 |
Ca2+ in 53H2Oe | 19.3 | 5 | 11.41 | 23.3 | |
Ca2+ in 199H2Of | 21.8 | 4 | 11.26 | 43.2 (22.0k) | |
Ca2+ in 999H2Og | 1199 | 80 | 217.1 | ||
Ca2+ in 999H2Oh | 1199 | 44 | 10.5 | ||
Ca2+ in 255H2Oi | 1000 | 39 |
![]() | ||
Fig. 2 Average number of water exchange events at structurally different surface calcium sites at 340 K for three consecutive nanoseconds of production. The error bars indicate the standard deviation of the number of exchange events per ∼80 ps of production within each nanosecond. Note that for face sites, error bars are either equal to or smaller than the squares. |
Table 1 shows that for calcium ions at the flat face sites the observed water exchange frequency is nearly two orders of magnitude lower than for the hydrated calcium ion, which has been obtained from the simulation of one calcium ion in 2027 water molecules. The water exchange frequency at the surface significantly increases with decreasing surface coordination, i.e., going from face to edge to corner site, approaching the values for the hydrated calcium ion computed using our shell model potential, first principles Car–Parrinello MD simulations43 and hybrid ab initio (Hartree–Fock) quantum mechanical/molecular mechanical simulations.41 As expected, the water exchange events also increase with increasing temperature (Table S4, ESI†) and water exchange frequencies for carbonate oxygen atoms are several orders of magnitude higher than for surface calcium (Table S5, ESI†), thereby suggesting very strongly that calcium de-hydration is indeed the rate-determining step in the growth and dissolution of calcium carbonate.
At all temperatures studied, calcium water exchange frequencies increase in the same order: face < pit corners < obtuse edges = obtuse corners < acute edges < acute corners (Fig. 2 for 340 K). The same trend in water exchange frequencies was previously observed by Cooke and Elliott for calcite nanoparticles24 and for the stepped versus flat (104) surface, although their absolute frequencies at 300 K were slightly lower than observed here, which could be due to their different potential parameters to describe water–water interaction (Lennard-Jones22) and method employed to define exchange events (Impey method47). In fact, the value of the mean residence time of water in the first hydration shell of Ca2+ varies by approximately 20 ps if the same MD trajectory is analysed using either the Impey procedure or the “direct” method.41
Our results show that, if we were to follow the step-wise desolvation of a single fully hydrated calcium to a bulk lattice calcium during calcite growth, we would see its desolvation slow down: the transition of this calcium from a kink into an edge will be faster than its subsequent transition from an edge to a face position, while the final step from face to bulk calcium will be slower still. Moreover, because the water exchange frequency at the obtuse sites is always lower than the water exchange frequency at the acute sites (see Table 1), calcium desolvation will also be slower when calcium enters calcite via an obtuse step than through an acute step. Furthermore, if carbonate and bicarbonate attachment frequencies are indeed controlled by surface calcium dehydration,7,8 our results suggest that obtuse edge movement might be slowest during dissolution and spiral growth. This observation agrees well with atomic force microscopy (AFM) data under certain growth conditions (e.g. at Ca/CO3 activity ratios <1),13 but not consistently.9,10,13,19,48 To investigate this inconsistent behaviour of step velocities, and translate the MD simulations result to the macroscopic scale, we normalized the water exchange frequencies for different surface calcium sites (Table 1) and used them to refine a process-based calcite growth kinetics model.7
Given the mineralogical and biological significance of calcium, we have also reported in Table 1 the mean residence time (MRT) of water molecules in the first coordination sphere of Ca2+ computed from the potential model used in the present work, as well as the values obtained with other classical and first principles MD simulations. Although it is not the aim of this work to conduct an in-depth assessment of the available potential models for the simulation of the dynamical properties of calcium in water, it is interesting to note that the frequencies of water exchange obtained by the potential model used in this study are in better agreement with ab initio results than those obtained using a Lennard-Jones description for the calcium–oxygen interaction, which give exchange frequencies which are too slow.
In Fig. 3, the refined overall attachment frequencies for (bi-)carbonate and calcium ions at the acute and obtuse edges are compared with previously reported attachment frequencies. The previously reported frequencies all resulted from fitting AFM data by kinetic growth models for binary electrolyte crystals6 of varying levels of mechanistic detail. While Stack and Grantham10 applied the binary electrolyte growth model directly, Nielsen et al.48 extended it to include isotopic exchange mechanisms and Wolthers et al.7 extended the binary to ternary model to include bicarbonate into a consequently pH-dependent growth model. It must be noted that all models have a degree of interdependency amongst the parameters: a combination of lower attachment frequencies with lower detachment frequencies may lead to an equally good fit to step velocities as a set of higher frequencies. This explains, for example the large difference in frequencies of Stack and Grantham10 and Wolthers et al.,7 that were partly based on fitting the same data.
![]() | ||
Fig. 3 Comparison of carbonate, k[CO32−] (a), and calcium, k[Ca2+] (b), ion attachment frequencies at acute and obtuse edge sites derived from fitting AFM step velocities and constrained by MD simulations. [1] This study, plotting ![]() ![]() |
Because the MD simulation results have been used here to constrain the attachment frequencies, the refined process-based growth model parameters do not show this interdependency. The refined attachment frequencies show a clear but subtle variation with surface topography (in log scale, Fig. 3) and are in general agreement with the previously published values, although the interdependency amongst the latter makes a more detailed comparison ineffective.
We have calculated kink formation frequencies, step velocities and bulk growth rates for the specified sites at the (104) surface using the refined model (Tables S1–S3, ESI†). Fig. 4 shows the new results of kink formation at actual experimental conditions.46 Frequencies of kink formation show a variation with surface structure of up to two orders of magnitude, and this variation depends on solution composition (Fig. 4). The kink formation frequency differences are amplified by solution stoichiometry and degree of supersaturation, Ω, because solution composition strongly affects the kinetic ionic ratio (ri in Table S1, ESI†) and the anion detachment frequency (
B, Table S2, ESI†) in the model.
![]() | ||
Fig. 4 Kink formation rate variation calculated by the process-based growth model7 refined by assuming attachment frequencies of carbonate and bicarbonate ions that are identical to normalized water exchange frequencies for the various surface calcium sites (Table 1). Model conditions were identical to experimental conditions of Nehrke et al.46 Kink formation rates were calculated using the equation for ic in Table S1, ESI† and model parameter as listed in Table S2, ESI.† Note the different scales of the y-axes. |
Acute and obtuse step velocities, calculated using the refined model, together with measured AFM velocities10,13 are presented in Fig. 5. In the refined model, obtuse-step velocities agree well with measurements, in contrast to acute step velocities. Measured acute step velocities are generally about an order of magnitude lower than model acute step velocities (open squares and diamonds in Fig. 5). Some of the observed discrepancies may be related to small variations in experimental conditions during the experiments. The calcite surface speciation and process-based growth model calculations were performed on measured input solution chemistry data and initial calcite surface area. Small variations in solution chemistry over the course of the measurements, such as changes in solution composition due to gas exchange, will alter solution stoichiometry and Ω, and may lead to orders of magnitude differences in kink formation rate (Fig. 4) and therefore model step velocities. Another possible explanation for the remaining discrepancy between experimental and calculated acute step velocities might be related to the more dynamic environment at the acute edges (see values of log kw at acute sites in Table 1). Additives and electrolyte ions are known to change the solvation environment.43,50,51 The higher water exchange frequencies at acute edges compared to obtuse steps is likely to increase the probability of additive attachment,52 enhancing or inhibiting step propagation, thus narrowing the range of measured step velocities compared to modelled step velocities.
![]() | ||
Fig. 5 Absolute, experimentally measured and process-based growth model step velocities, calculated for the experimental conditions during measurements of (■,◆) obtuse and (□,◇) acute step velocities on growing calcite surfaces observed with AFM at {Ca2+}/{CO32−} = raq ≈ 1.0, pH 8.5 at 1 < Ω < 3.6 (squares13) and 7 < pH < 9, 10−3 < raq < 103, in NaCl background electrolyte and 2.2 < Ω < 2.7 (diamonds10). Solid line indicates x = y; RSS = 2.1 × 10−17 m s−1. Step velocities were calculated using the equation Vc= ρc·a·uc, where a is the mean ionic diameter for calcite (Table S3, ESI†) and equations for ρc and uc as listed in Table S1, ESI.† |
The variation in bulk growth rate for the various calcite (104) surface features with solution stoichiometry is illustrated in Fig. 6. Nehrke et al.46 observed an optimum in growth rate at a near one-to-one solution stoichiometry at constant pH of 10.2 and a constant Ω of 5 and 16 (open and closed triangles in Fig. 6, respectively). The different model lines in Fig. 6 represent bulk growth rates calculated using the attachment frequencies for the indicated surface sites and therefore show the observed variation in bulk growth rates that can be caused by variations in surface topography at identical solution composition. The model lines also indicate that, depending on surface topography, the dependency of bulk growth rate on solution stoichiometry (at constant pH and Ω) may vary strongly. As can be seen in Fig. 6, the model lines generally agree with the bulk growth rate data, except at Ω 5 near a 1
:
1 stoichiometry, where even the fastest model growth rates underestimate the observed growth rates. Some of the discrepancy may again originate from small variations between actual experimental solution composition and the input solution composition, which were used for the calcite surface speciation and process-based growth rate modelling, although these experiments were performed in a gas-tight closed system. Another explanation may be related to a change in actual calcite growth mechanism depending on solution stoichiometry.
![]() | ||
Fig. 6 Comparison of absolute bulk growth rates between measured in individual experiments46 and modelled data at experimental conditions for the various calcite (10![]() |
It is generally known that at low to moderate degrees of supersaturation, such as Ω 5 (Fig. 6a), calcite growth is dominated by spiral growth.13 Such a growth spiral consists of four step edge directions, two mirrored acute steps and two mirrored obtuse steps,17 and bulk growth rate at such conditions results from advancement of these steps, as well as kink formation at these steps.52 At higher degrees of supersaturation, such as Ω 16 (Fig. 6b), the formation of new growth islands (two-dimensional surface nucleation) becomes increasingly important.13 More recently, it has also been shown that two-dimensional surface nucleation may occur even at low to moderate degrees of supersaturation, depending on solution stoichiometry.50 Surface nucleation mechanisms are not implemented in the process-based calcite growth kinetics model and may explain some of the discrepancies between calculated and measured bulk growth data in Fig. 6. Nevertheless, the model rates in Fig. 6 generally agree with measured rates, suggesting that some of the variation in bulk growth rates obtained at identical solution conditions might have been caused by variations in seed-surface topography.
1) water exchange frequencies for surface calcium ions vary by up to one order of magnitude,
2) kink formation frequencies vary by up to two orders of magnitude, and
3) bulk growth rates can vary by more than two orders of magnitude at the same solution composition;
4) the dependency on solution stoichiometry of growth rate changes as well.
Our calculated bulk growth rates and obtuse step velocities generally compare favourably with measured obtuse step advancement rates. In the case of acute steps, disparity between process-based model and experiment may be explained by a higher sensitivity of the acute steps to additives or electrolyte ions due to their higher water exchange frequencies.
Since water exchange kinetics are key indicators of chemical reactivity,53 our results imply surface topography-dependent reactivity of calcite towards constituent and contaminant ions, which has important implications in, for example, the fields of biomineralization and mineral engineering. For example, when (nano-)calcite is synthesised with specific morphologies and sizes,54 our results may help aims to produce (nano-)particles with specific kinetic properties. Moreover, the effect of surface topography on the dynamics of interfacial water is likely to affect isotope fractionation48 and adsorption of foreign elements by calcite.52 For example, topography-induced variations in dynamics at the surface may cause discrepancies between experimental trace metal uptake studies if different calcite seed materials or different preparation techniques are used. This could explain, for example, the ongoing debate on Cd uptake by calcite, where Cd uptake mechanisms are related to calcite surface dynamics.55
This study shows that atomic scale information that can be obtained from MD calculations, such as water exchange kinetics at structurally different surface sites, can be “scaled up” to develop predictive models describing the reactivity of minerals from a macroscopic point of view.
Footnote |
† Electronic supplementary information (ESI) available: Tables listing key equations, parameters and nomenclature of the process-based calcite growth model. Tables listing water exchange frequencies at 300 and 320 K for surface calcium and 300, 320 and 340 K for surface oxygen in carbonate groups. See DOI: 10.1039/c3ce40249e |
This journal is © The Royal Society of Chemistry 2013 |