The evolution of crystal growth morphologies was explored using a Monte Carlo model. The model combines diffusive transport in the nutrient with thermally activated and local configuration-dependent steps for attachment, surface diffusion and detachment. The solid pair interaction energy, φ, temperature, T, and chemical potential difference between nutrient and surface, Δμ, were used as input parameters. For a given set of simulation parameters, it was found that, due to the competing effects of bulk diffusion and interface kinetics, there was a critical size beyond which a crystal could not retain its macroscopically faceted shape. This critical size scales linearly with the mean free path in the nutrient, a. The stabilization of the growth shape by anisotropies in the kinetics was reduced through thermal and kinetic roughening. Thus, further increases either in T or Δμ, and/or decreased in φ, cause successive transitions from faceted to compact dendritic and side-branched dendritic morphologies. On interfaces with emerging screw dislocations the spiral growth mechanism dominates at low Δμ. When a was comparable to the lattice constant, the combination of bulk and surface diffusion reduces the terrace width near the centre of a growth spiral. At elevated T and Δμ, 2D nucleation-controlled growth could dominate in the more readily supplied corner and edge regions of a facet, while spiral growth prevails in its centre.

Monte Carlo Studies of the Interdependence of Crystal Growth Morphology, Surface Kinetics and Bulk Transport. R.Xiao, J.I.D.Alexander, F.Rosenberger: Faraday Discussions, 1993, 95, 85-95