Hostname: page-component-848d4c4894-nmvwc Total loading time: 0 Render date: 2024-06-18T16:42:47.674Z Has data issue: false hasContentIssue false

Theory of large-scale electronic structure calculation and nanostructures formed in silicon cleavage simulation: surface reconstruction, step and bending

Published online by Cambridge University Press:  01 February 2011

Takeo Hoshi
Department of Applied Physics, University of Tokyo, Hongo, Bunkyo-ku, Tokyo, Japan. CREST, JST, Motomachi, Kawaguchi-shi, Saitama, Japan.
Yusuke Iguchi
Department of Applied Physics, University of Tokyo, Hongo, Bunkyo-ku, Tokyo, Japan.
Takeo Fujiwara
Department of Applied Physics, University of Tokyo, Hongo, Bunkyo-ku, Tokyo, Japan. CREST, JST, Motomachi, Kawaguchi-shi, Saitama, Japan.
Get access


Several methodologies are developed for large-scale atomistic simulations with fully quantum mechanical description of electron systems. The important methodological concepts are (i) generalized Wannier state and (ii) Krylov subspace. Test calculations are done with upto 106 atoms using a standard workstation. As a practical nanoscale calculation, the 10 nm scale structure in cleavage process of silicon is investigated. Discussions are focused mainly on the stable (experimentally observed) (111)-(2 × 1) cleavage process. Surface reconstruction and step formation are observed and compared with experiments. These processes are analyzed by the quantum mechanical freedoms of electron system, such as the local density of states (LDOS). We also observed the bending of cleavage plane into the (111) plane from other planes, which is a direct evidence of the stability of the (111) cleavage mode. Finally, several common aspects between cleavage and other phenomena are discussed from a general viewpoint of 10 nm scale structure.

Research Article
Copyright © Materials Research Society 2005

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)



[1] Griffith, A. A., Philos. Trans. R. Soc. London, Ser. A 221, 163 (1920).Google Scholar
[2] Mott, N. F., Engineering 165, 16 (1948).Google Scholar
[3] As general REFERENCES of brittle fracture, Freund, L. B., ‘Dynamic fracture mechanics’, Cambridge university press (1989);Google Scholar
Lawn, B., ‘Fracture of brittle solids’, 2nd ed., Cambridge university press (1993).Google Scholar
[4] As a review, Neddermeyer, H., Rep. Prog. Phys, 59, 701 (1996).Google Scholar
[5] Pandey, K. C., Phys. Rev. Lett. 47, 1913 (1981).Google Scholar
[6] Northrup, J. E. and Cohen, M. L., Phys. Rev. Lett. 49, 1349 (1982).Google Scholar
[7] Feenstra, R. M., Thompson, W. A. and Fein, A. P. Phys. Rev. Lett. 56, 608 (1986);Google Scholar
Stroscio, J. A., Feenstra, R. M. and Fein, A. P. Phys. Rev. Lett. 57, 2579 (1986).Google Scholar
[8] Ancilotto, F., Andreoni, W., Selloni, A., Car, R., and Parrinello, M., Phys. Rev. Lett. 65, 3148 (1990).Google Scholar
[9] Huang, Y. M., Spence, J. C. H., Sankey, O. F., and Adams, G. B., Surf. Sci. 256, 344 (1991).Google Scholar
[10] Spence, J. C. H., Huang, Y. M., and Sankey, O., Acta Metall. Mater. 41, 2815 (1993).Google Scholar
[11] Rohlfing, M. and Louie, S. G., Phys. Rev. Lett. 83, 856 (1999).Google Scholar
[12] See the following recent work or REFERENCES therein; Stekolnikov, A. A., Furthmüller, J., and Bechstedt, F., Phys. Rev. B 65, 115318 (2002).Google Scholar
[13] Henzler, M., Surf. Sci. 36, 109 (1973).Google Scholar
[14] As well as the (111) plane, the (110) cleavage plane is also experimentally observed but is less favorable than the (111) cleavage plane. Experimental STM images are found in the following paper; Lutz, M. A., Feenstra, R. M., and Chu, J. O., Surf. Sci. 328, 215 (1995). See also the recent theoretical work and REFERENCES therein;Google Scholar
Pérez, R. and Gumbsch, P., Phys. Rev. Lett. 84, 5347 (2000).Google Scholar
[15] The cleaved Ge(111) surface also shows the same (2 × 1) structure as the Si surface, while the ground state surface structure is not the (7 × 7) structure but the c(2 × 8) structure. [4, 12]Google Scholar
[16] Hoshi, T. and Fujiwara, T., Phys, J.. Soc. Jpn. 69, 3773 (2000).Google Scholar
[17] Hoshi, T. and Fujiwara, T., Surf. Sci. 493, 659 (2001).Google Scholar
[18] Hoshi, T. and Fujiwara, T., Phys, J.. Soc. Jpn. 72, 2429 (2003).Google Scholar
[19] Geshi, M., Hoshi, T., and Fujiwara, T., J. Phys. Soc. Jpn. 72, 2880 (2003).Google Scholar
[20] Takayama, R., Hoshi, T., and Fujiwara, T., J. Phys. Soc. Jpn. 73, 1519 (2004).Google Scholar
[21] Hoshi, T., Iguchi, Y. and Fujiwara, T., preprint (cond-mat/0409142).Google Scholar
[22] The Physical Society of Japan 2004 Autumn Meeting, 15aTB-3; Takayama, R., Hoshi, T., Sogabe, T., Zhang, S.-L. and Fujiwara, T., in preparation.Google Scholar
[23] Kohn, W., Phys. Rev. B7, 4388 (1973).Google Scholar
[24] Kohn, W., Chem. Phys. Lett. 208, 167 (1993).Google Scholar
[25] In several papers, the Wannier states are constructed in a post process of the standard electronic structure calculations with eigen states. For silicon crystal, the ab initio calculation was carried out using plane wave bases (Marzari, N. and Vanderbilt, D., Phys. Rev. B56, 12847 (1997)). or using muffin-tin orbital basesGoogle Scholar
(Andersen, O. K., Saha-Dasgupta, T., and Ezhov, S., Bull. Mater. Sci. 26, 19 (2003)). Though their methodologies in constructing Wannier states are quite different, the resultant wave functions show the character of sp 3 bonding orbital with ‘tail’. The ‘tail’ part of the wave function has node structures on neighbor bond sites, as a consequence of the orthogonality between the Wannier states.Google Scholar
[26] As textbooks, Golub, G. H. and Van Loan, C. F., Matrix Computations, third ed., Johns Hopkins University Press, London (1996);Google Scholar
van der Vorst, H. A., Iterative Krylov methods for large linear systems, Cambridge University Press (2003).Google Scholar
[27] Feenstra, R. M. and Stroscio, J. A., Phys. Rev. Lett. 59, 2173 (1987).Google Scholar
[28] Tokumoto, H., Wakiyama, S., Miki, K., Murakami, H., Okayama, S., and Kajimura, K., J. Vac. Sci. Technol. B 9 695 (1991);Google Scholar
Komeda, T., Gwo, S., and Tokumoto, H., Jpn. J. Appl. Phys. 35, 3724 (1996).Google Scholar
[29] Mera, Y., Hashizume, T., Maeda, K., and Sakurai, T., Ultramicroscopy 42–44, 915 (1992).Google Scholar
[30] Chadi, D. J. and Chelikowsky, J. R., Phys. Rev. B 24, 4892 (1981).Google Scholar
[31] Shklyaev, A. A. and Ichikawa, M., Phys. Rev. B 65 045307 (2002).Google Scholar