Sulfate and Molybdate Incorporation at the Calcite Water Interface: Insights from Ab Initio Molecular Dynamics

Sulfur and molybdenum trace impurities in speleothems (stalagmites and stalactites) can provide long and continuous records of volcanic activity, which are important for past climatic and environmental reconstructions. However, the chemistry governing the incorporation of the trace element-bearing species into the calcium carbonate phases forming speleothems is not well understood. Our previous work has shown that substitution of tetrahedral oxyanions [X04]2- (X = S and Mo) replacing [CO3]- in CaCO3 bulk phases (except perhaps for vaterite) is thermodynamically unfavorable with respect to the formation of competing phases, due to the larger size and different shape of the [X04]2'- tetrahedral anions in comparison with the flat [CO3]2- calcite anions, which implied that most of the incorporation would happen at the surface rather than at the bulk of the mineral. Here, we present an ab initio molecular dynamics study, exploring the incorporation of these impurities at the mineral water interface. We show that the oxyanion substitution at the aqueous calcite (10.4) surface is clearly favored over bulk incorporation, due to the lower structural strain on the calcium carbonate solid. Incorporation at surface step sites is even more favorable for both oxyanions, thanks to the additional interface space afforded by the surface line defect to accommodate the tetrahedral anion. Differences between sulfate and molybdate substitutions can be mostly explained by the size of the anions. The molybdate oxyanion is more difficult to incorporate in the calcite bulk than the smaller sulfate oxyanion. However, when molybdate is substituted at the surface, the elastic cost is avoided because the oxyanion protrudes out of the surface and gains stability via the interaction with water at the interface, which in balance results in more favorable surface substitution for molybdate than for sulfate. The detailed molecular-level insights provided by our calculations will be useful to understand the chemical basis of S- and Mo-based speleothem records