Boundary-reaction-diffusion model for oscillatory zoning in binary crystals grown from solution

Oscillatory Zoning (OZ) is a phenomenon exhibited by many geologically formed crystals. It is characterized by quasiperiodic oscillations in the composition of a solid solution, caused by self-organization. We present a model for OZ. The growth mechanism applied includes species diffusion through the solution bulk, particle adsorption, surface diffusion, and subsequently desorption or incorporation into the crystal. This mechanism, in particular, can provide the synchronization effects necessary to reproduce the layered structure of experimentally obtained crystals, lacking in other models. We conduct a linear stability analysis combined with numerical simulations. Our results reproduce the experimental findings with respect to the patterns formed and a critical supersaturation necessary for OZ to occur.