Modeling crystal growth from solution with molecular dynamics simulations: approaches to transition rate constants

The feasibility of using the molecular dynamics (MD) simulation technique to study crystal growth from solution quantitatively, as well as to obtain transition rate constants, has been studied. The dynamics of an interface between a solution of Lennard-Jones particles and the (100) face of an fcc lattice comprised of solute particles have been studied using MD simulations, showing that MD is, in principle, capable of following growth behavior over large supersaturation and temperature ranges. Using transition state theory, and a nearest-neighbor approximation growth and dissolution rate constants have been extracted from equilibrium MD simulations at a variety of temperatures. The temperature dependence of the rates agrees well with the expected transition state theory behavior.