Reconstruction of silicon surfaces: a stochastic optimization problem

Abstract

Over the last two decades, scanning tunnelling microscopy (STM) has become one of the most important ways to investigate the structure of crystal surfaces. STM has helped achieve remarkable successes in surface science such as finding the atomic structure of Si(111) and Si(001). For high-index Si surfaces the information about the local density of states obtained by scanning does not translate directly into knowledge about the positions of atoms at the surface. A commonly accepted strategy for identifying the atomic structure is to propose several possible models and analyze their corresponding {\em simulated} STM images for a match with the experimental ones. However, the number of good candidates for the lowest-energy structure is very large for high-index surfaces, and heuristic approaches are not likely to cover all the relevant structural models. In this article, we take the view that finding the atomic structure of a surface is a problem of stochastic optimization, and we address it as such. We design a general technique for predicting the reconstruction of silicon surfaces with arbitrary orientation, which is based on parallel-tempering Monte Carlo simulations combined with an exponential cooling. The advantages of the method are illustrated using the Si(105) surface as example, with two main results: (a) the correct single-step rebonded structure [e.g., Fujikawa {\em et al.}, Phys. Rev. Lett. 88, 176101 (2002)] is obtained even when starting from the paired-dimer model [Mo {\em et al.}, Phys. Rev. Lett. 65, 1020 (1990)] that was assumed to be correct for many years, and (b) we have found several double-step reconstructions that have lower surface energies than any previously proposed double-step models.

13 pages, 9 figures

References in corpus (9)

• Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms
• Exchange Monte Carlo Method and Application to Spin Glass Simulations
• Parallel Tempering Algorithm for Conformational Studies of Biological Molecules
• Multidimensional replica-exchange method for free-energy calculations
• The incomplete beta function law for parallel tempering sampling of classical canonical systems
• Structure and Stability of Si(114)-(2x1)
• Finding Low-Temperature States with Parallel Tempering, Simulated Annealing and Simple Monte Carlo
• Strain induced stabilization of stepped Si and Ge surfaces near (001)
• Structure and stability of the Si(105) surface

Cited by in corpus (12)

• Parallel Tempering: Theory, Applications, and New Perspectives
• Classical potential describes martensitic phase transformations between the $\alpha$, $\beta$ and $\omega$ titanium phases
• Finding the reconstructions of semiconductor surfaces via a genetic algorithm
• Thermodynamic Stability of Molybdenum Oxycarbides Formed from Orthorhombic Mo2C in Oxygen-rich Environments
• Structure of Si(114) determined by global optimization methods
• Systematic investigation of the structure of the Si(553)-Au surface from first principles
• Interfacial Energy of Copper Clusters in Fe-Si-B-Nb-Cu alloys
• Model reconstructions for the Si(337) orientation
• Evolutionary approach for finding the atomic structure of steps on stable crystal surfaces
• Determination of ground-state and low-energy structures of perovskite superlattices from first principles
• The stability of strained H:Si(105) and H:Ge(105) surfaces
• Adsorbate phase transitions on nanoclusters from nested sampling