We report low-energy locally stable structures for the clusters Si 20 and Si 21 . The structures were obtained by performing geometry optimizations within the local density approximation. Our calculated binding energies for these clusters are larger than any previously reported for this size regime. To aid in the experimental identification of the structures, we have computed the full vibrational spectra of the clusters, along with the Raman and IR activities of the various modes using a recently developed first-principles technique. These represent, to our knowledge, the first calculations of Raman and IR spectra for Si clusters of this size.
͓S0163-1829͑96͒06227-3͔Interest in the properties of Si clusters has grown in recent years, both as a result of the drive toward design and miniaturization in the electronics industry, and due to the many interesting and surprising properties of atomic clusters that have emerged from early experimental work.1 A crucial problem in this area involves the structure of the intermediate-sized clusters. The lowest-energy structures for small Si clusters (Nр13͒ appear well established by calculations 2-4 and for Nр10, experimental Raman spectra have confirmed these geometries. 5 At the other end of the size spectrum, large clusters are expected to take on bulklike structures. Between these limits, the structures of the intermediate-sized clusters are not known. The basic problem for theory is the complexity of the energy surfaces for these clusters and the vast number of local minima they are likely to contain. Various search strategies, coupled with a number of theoretical approaches, have been used to explore clusters in this size range, resulting in a wide variety of cluster models. [6][7][8][9][10][11][12] As was true for the clusters with Nр10, definitive identification of the structures is likely to come only through a combination of theory and experiment.In this paper we describe low-energy models for Si 20 and Si 21 , obtained using the local density approximation ͑LDA͒. These two structures have larger binding energies than alternative models that have appeared in the literature, making them leading candidates for minimum-energy structures. To further investigate the nature of these models and to encourage experimental work on clusters in this size range, we have computed the full vibrational spectrum for the clusters and the Raman and IR activities of the various modes. These calculations provide signatures that could be used to identify the clusters in experiments.The structures we report here were obtained using a standard conjugate-gradient algorithm with atomic forces calculated within the LDA. Our computational approach has been presented elsewhere. 15 Gaussian orbitals are used to represent the electronic states of the clusters and a numerical scheme is used to obtain accurate LDA forces and total energies. 15 The cluster geometries were relaxed until the largest force on every atom dropped below 0.001 hartree/ Bohr. Calculations performed here utilized a basis s...