Early stages of carbonate mineralization revealed from molecular simulations: Implications for biomineral formation

The carbonate mineral constituents of many biomineralized products, formed both in and ex vivo, grow by a multi-stage crystallization process that involves the nucleation and structural reorganization of transient amorphous phases. The existence of transient phases and cluster species has significant implications for carbonate nucleation and growth in natural and engineered environments, both modern and ancient. The structure of these intermediate phases remains elusive, as does the nature of the disorder to order transition, however, these process details may strongly influence the interpretation of elemental and isotopic climate proxy data obtained from authigenic and biogenic carbonates. While molecular simulations have been applied to certain aspects of crystal growth, studies of metal carbonate nucleation are strongly inhibited by the presence of kinetic traps that prevent adequate sampling of the potential landscape upon which the growing clusters reside within timescales accessible by simulation. This research addresses this challenge by marrying the recent Kawska-Zahn (KZ) approach to simulation of crystal nucleation and growth from solution with replica-exchange molecular dynamics (REMD) techniques. REMD has been used previously to enhance sampling of protein conformations that occupy energy wells that are separated by sizable thermodynamic and kinetic barriers, and is used here to probe the initial formation and onset of order within hydrated calcium and iron carbonate cluster species during nucleation. Results to date suggest that growing clusters initiate as short linear ion chains that evolve into two- and three-dimensional structures with continued growth. The planar structures exhibit an obvious 2d lattice, while establishment of a 3d lattice is hindered by incomplete ion desolvation. The formation of a dehydrated core consisting of a single carbonate ion is observed when the clusters are ~0.75 nm. At the same size a distorted, but discernible calcite-type lattice is also apparent. Continued growth results in expansion of the dehydrated core, however, complete desolvation and incorporation of cations into the growing carbonate phase is not achieved until the cluster grows to ~1.2 nm. Exploration of the system free energy along the crystallization path reveals "special" cluster sizes that correlate with ion desolvation milestones. The formation of these species comprise critical bottlenecks on the energy landscape and for the establishment of order within the growing clusters.