Reconstruction of silicon surfaces: a stochastic optimization problem

TL;DR

This paper introduces a stochastic optimization method using parallel-tempering Monte Carlo simulations to accurately predict silicon surface reconstructions, overcoming limitations of heuristic approaches and discovering lower-energy structures.

Contribution

The authors develop a novel stochastic optimization technique for silicon surface reconstruction prediction, demonstrating its effectiveness on Si(105) and identifying more stable structures.

Findings

Correct single-step rebonded structure obtained from various initial models.

Discovered double-step reconstructions with lower surface energies than previous models.

Method successfully explores complex surface energy landscapes.

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.