Monte Carlo studies of liquid semiconductor surfaces: Si and Ge

We study the liquid-vapor interface of Si and Ge, and of Si doped with impurities, using the empirical Stillinger-Weber two- and three-body potentials. The surface tension of the pure elements is calculated by use of a direct evaluation of the free energy required to create the surface. For Si, both the surface tension and its temperature derivative are in good agreement with experiment. To within numerical accuracy, the free surfaces of both Si and Ge have a monotonically decreasing density, with a 10%-90% surface width of ∼2.2 Å in both cases. When large and small model impurities are introduced into pure liquid Si, they are found in the simulations to migrate, respectively, towards the surface and away from the surface. This behavior is consistent with the interpretation that impurities with atoms larger than Si tend to lower the surface tension.