A Monte Carlo simulation method for {111} surfaces of silicon and other diamond-cubic materials

We introduce a modified, solid-on-solid model for Monte Carlo computer simulation studies of equilibrium, growth and etching phenomena on the {111} surfaces of silicon and other diamond-cubic materials. Incorporation of the actual diamond-cubic crystal structure allows for the inclusion of possible effects of the structure\u27s rotational symmetry, coordination and the layer-stacking periodicity (normal to the surface) on the surface properties and dynamics. Results of extensive equilibrium simulations, including surface energies, surface specific heats, step energies and surface roughness, are presented to characterize the basic model. A broad peak in the surface specific heat and the universal behavior of the height-difference correlation function indicate a roughening transition at temperature TR. A preroughening transition is indicated by a sharp peak in the surface specific heat and a dramatic decrease in the step energy at a temperature TPR ≈ 0.43TR. The energies of stable and unstable steps differ significantly below TPR but become isotropic above TPR. Initially straight, unstable steps are observed to spontaneously microfacet below TPR. Preliminary simulations of growth and etching of surfaces, both with and without intersecting dislocations, are also presented