Molecular dynamics simulation of crystal dissolution from calcite steps

Molecular-dynamics simulations were used to model two stepped {10 (1) over bar 4} surfaces of the calcium carbonate polymorph calcite. The acute monat omic steps were found to be more stable than the obtuse monatomic steps. Th e initial stages of dissolution from the steps were considered in vacuo and in water. In vacuo CaCO3 was shown to dissolve preferentially from the obt use step. In aqueous environment both stepped surfaces are stabilized by th e presence of the water molecules although the relative stabilities remain similar. Using potential parameters that reproduce experimental enthalpies of the dissolution of calcite crystal, the formation of the double kinks on the obtuse step is shown to cost less energy than dissolution from the acu te step, probably due to the lower stability of the obtuse surface. The sim ulations suggest that formation of thr kink sites on the dissolving edge of the obtuse step of calcite is the rate determining step and this edge is p redicted to dissolve preferentially, which is in agreement with experimenta l findings of calcite dissolution under aqueous conditions. [S0163-1829(99) 12039-3].