M. Wolthers

The surface chemistry of divalent metal carbonate minerals; a critical assessment of surface charge and potential data using the charge distribution multi-site ion complexation model

The Charge Distribution MUltiSite Ion Complexation or CD–MUSIC modeling approach is used to describe the chemical structure of carbonate mineral-aqueous solution interfaces. The new model extends existing surface complexation models of carbonate minerals, by including atomic scale information on the surface lattice and the adsorbed water layer. In principle, the model can account for variable proportions of face, edge and kink sites exposed at the mineral surface, and for the formation of inner- and outer-sphere surface complexes. The model is used to simulate the development of surface charges and surface potentials on divalent carbonate minerals as a function of the aqueous solution composition. A comparison of experimental data and model output indicates that the large variability in the observed pH trends of the surface potential for calcite may in part reflect variable degrees of thermodynamic disequilibrium between mineral, solution and, when present, gas phase during the experiments. Sample preparation and non-stoichiometric surfaces may introduce further artifacts that complicate the interpretation of electrokinetic and surface titration measurements carried out with carbonate mineral suspensions. The experimental artifacts, together with the high sensitivity of the model toward parameters describing hydrogen bridging and bond lengths at the mineral-water interface, currently limit the predictive application of the proposed CD–MUSIC model. The results of this study emphasize the need for internally consistent experimental data sets obtained with well-characterized mineral surfaces and in situ aqueous solution compositions (that is, determined during the charge or potential measurements), as well as for further molecular dynamic simulations of the carbonate mineral-water interface to better constrain the bond lengths and the number plus valence contribution of hydrogen bridges associated with different structural surface sites.

Coupled dynamics of iron and phosphorus in sediments of an oligotrophic coastal basin and the impact of anaerobic oxidation of methane.

Studies of phosphorus (P) dynamics in surface sediments of lakes and coastal seas typically emphasize the role of coupled iron (Fe), sulfur (S) and P cycling for sediment P burial and release. Here, we show that anaerobic oxidation of methane (AOM) also may impact sediment P cycling in such systems. Using porewater and sediment profiles for sites in an oligotrophic coastal basin (Bothnian Sea), we provide evidence for the formation of Fe-bound P (possibly vivianite; Fe3(PO4)2 .8H2O) below the zone of AOM with sulfate. Here, dissolved Fe2+ released from oxides is no longer scavenged by sulfide and high concentrations of both dissolved Fe2+ (>1 mM) and PO4 in the porewater allow supersaturation with respect to vivianite to be reached. Besides formation of Fe(II)-P, preservation of Fe-oxide bound P likely also contributes to permanent burial of P in Bothnian Sea sediments. Preliminary budget calculations suggest that the burial of Fe-bound P allows these sediments to act as a major sink for P from the adjacent eutrophic Baltic Proper.

Calcite surface structure and reactivity: molecular dynamics simulations and macroscopic surface modelling of the calcite-water interface.

Calcite-water interactions are important not only in carbon sequestration and the global carbon cycle, but also in contaminant behaviour in calcite-bearing host rock and in many industrial applications. Here we quantify the effect of variations in surface structure on calcite surface reactivity. Firstly, we employ classical Molecular Dynamics simulations of calcite surfaces containing an etch pit and a growth terrace, to show that the local environment in water around structurally different surface sites is distinct. In addition to observing the expected formation of more calcium-water interactions and hydrogen-bonds at lower-coordinated sites, we also observed subtle differences in hydrogen bonding around acute versus obtuse edges and corners. We subsequently used this information to refine the protonation constants for the calcite surface sites, according to the Charge Distribution MUltiSite Ion Complexation (CD-MUSIC) approach. The subtle differences in hydrogen bonding translate into markedly different charging behaviour versus pH, in particular for acute versus obtuse corner sites. The results show quantitatively that calcite surface reactivity is directly related to surface topography. The information obtained in this study is not only crucial for the improvement of existing macroscopic surface models of the reactivity of calcite towards contaminants, but also improves our atomic-level understanding of mineral-water interactions.

Variations in calcite growth kinetics with surface topography: molecular dynamics simulations and process-based growth kinetics modelling

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.

Mg/Ca Partitioning Between Aqueous Solution and Aragonite Mineral: A Molecular Dynamics Study

We have calculated the concentrations of Mg in the bulk and surfaces of aragonite CaCO(3) in equilibrium with aqueous solution, based on molecular dynamics simulations and grand-canonical statistical mechanics. Mg is incorporated in the surfaces, in particular in the (001) terraces, rather than in the bulk of aragonite particles. However, the total Mg content in the bulk and surface of aragonite particles was found to be too small to account for the measured Mg/Ca ratios in corals. We therefore argue that most Mg in corals is either highly metastable in the aragonite lattice, or is located outside the aragonite phase of the coral skeleton, and we discuss the implications of this finding for Mg/Ca paleothermometry.

How minerals dissolve

Local surface processes drive calcite dissolution [Also see Report by Laanait et al.] During mineral growth and dissolution, material (in the form of atoms, molecules, complexes, or clusters) is transported through water to or from the mineral surface. Even in a solution that is at equilibrium with respect to a given mineral, there is transport to and from the surface. Although this equilibrium transport is in balance over large scales, the rates of this transport can vary locally depending on the energy landscape of the mineral surface (1). On page 1330 of this issue, Laanait et al. (2) present an elegant way of measuring and visualizing how dissolution rates vary across the calcite surface with time and with dissolution mechanism.

In Situ XAS of the Solvothermal Decomposition of Dithiocarbamate Complexes

An in situ XAS study of the solvothermal decomposition of iron and nickel dithiocarbamate complexes was performed in order to gain understanding of the decomposition mechanisms. This work has given insight into the steps involved in the decomposition, showing variation in reaction pathways between the iron and nickel dithiocarbamates, and the non-innocent role of oleylamine as the solvent and capping agent in the reaction.

How minerals dissolve : Local surface processes drive calcite dissolution

Local surface processes drive calcite dissolution [Also see Report by Laanait et al.] During mineral growth and dissolution, material (in the form of atoms, molecules, complexes, or clusters) is transported through water to or from the mineral surface. Even in a solution that is at equilibrium with respect to a given mineral, there is transport to and from the surface. Although this equilibrium transport is in balance over large scales, the rates of this transport can vary locally depending on the energy landscape of the mineral surface (1). On page 1330 of this issue, Laanait et al. (2) present an elegant way of measuring and visualizing how dissolution rates vary across the calcite surface with time and with dissolution mechanism.