Home | About | Journals | Submit | Contact Us | Français |

**|**Nanoscale Res Lett**|**v.5(12); 2010**|**PMC2991238

Formats

Article sections

Authors

Related links

Nanoscale Res Lett. 2010; 5(12): 1873–1877.

Published online 2010 August 6. doi: 10.1007/s11671-010-9723-x

PMCID: PMC2991238

R Bergamaschini: ti.biminu.retam@inihcsamagreb.otrebor

Special Section(pp.1863-1951): 8th International Workshop on Epitaxial Semiconductors on Patterned Substrates and Novel Index Surfaces;Guest Editor: Stefano Sanguinetti

Received 2010 June 29; Accepted 2010 July 26.

Copyright ©2010 The Author(s)

This article has been cited by other articles in PMC.

We investigate ordered nucleation of Ge islands on pit-patterned Si(001) using an original hybrid Kinetic Monte Carlo model. The method allows us to explore long time-scale evolution while using large simulation cells. We analyze the possibility to achieve selective nucleation and island homogeneity as a function of the various parameters (flux, temperature, pit period) able to influence the growth process. The presence of an optimal condition where the atomic diffusivity is sufficient to guarantee nucleation only within pits, but not so large to induce significant Ostwald ripening, is clearly demonstrated.

Ge/Si(001) is often considered as the prototypical example of a system following the Stranski-Krastanow (SK) growth modality. The appearance of nanometric-sized, coherent Ge islands [1,2] following the formation of a thin wetting layer (WL), attracted widespread attention in view of the possible role that islands could play in developing future-generation devices. Two main problems were identified soon after the first experimental evidences of islands formation were reported in the literature. On a flat Si(001) substrate, indeed, islands tend to nucleate randomly (see, e.g., the discussion in Montalenti et al. [3]), precluding the possibility of obtaining ordered arrays. Secondly, under a wide range of deposition conditions, a bimodal distribution of islands (shallow {105} pyramids and steeper, multifaceted domes [4,5]) is obtained. Very interestingly, better lateral ordering and size homogeneity can be achieved by growing multistacked layers of Ge islands separated by Si spacer layers (SL) [6], the key role being played by the strain field originated by buried Ge islands at the surface of the outermost SL [7]. Recently, the possibility of enhancing lateral ordering directly from the first layer of deposited islands was also explored [8]. While in all the above reported references, some degree of lateral ordering was clearly demonstrated, much better results, both in terms of islands positioning and of size distribution, have been achieved by well-tuned patterning of the Si substrate. In particular, a suitable pit-patterning produced remarkable results under appropriate deposition conditions [9]. Highly ordered, unimodal distributions of dome islands are visible in the AFM images reported in Zhong and Bauer [9]. Further evidence of nice lateral ordering can be found also in Refs. [10-12].

In this work, we further explore ordering and homogeneity of Ge islands grown on pit-patterned Si(001). More specifically, we introduce an hybrid Kinetic Monte Carlo (h-KMC) method, developed in order to expand the typical length and time scales treatable in standard KMC, and we analyze the dependence of ordering and homogeneity on key parameters such as temperature, deposition flux and pit spacing.

On very general grounds, it is clear that reliable atomistic heteroepitaxial-growth simulations could strongly help in restricting the parameter-space to be experimentally sampled to obtain the desired result, e.g. in terms of islands distribution, positioning, etc. This is however a daunting task if one is willing to capture the whole atomic-scale complexity of Ge/Si(001) growth (with this respect, we notice that *prototypical* does not mean simple), including changes in surface reconstruction [13], long-ranged elastic fields [14], and huge spatial (experimentally determined island volumes easily reaching 10^{5} nm^{3} or more, pit extension and spacing being of the orders of hundreds of nm [9]) and time scales (experiments being performed at *human* time scales of several seconds). Despite the use of strong simplifications, such as solid-on-solid geometries [15] (which we shall also exploit), the gap between experiments and theory is still large. For example, standard KMC simulations [16-18] for growth on patterned substrates were reported for typical simulation-cell dimensions of 400 × 400 atomic sites or lower [19,20] and were restricted to very initial stages of growth, so that a direct comparison with experiments under realistic conditions was not attempted. Nevertheless, these simulations nicely demonstrated the existence of optimal ranges of temperature and deposition flux for obtaining homogeneous and laterally ordered distributions of islands. We therefore tried to keep the essential ingredients of the above approaches, using however a faster simulation method allowing for a reduction in the gap between simulation and experiments. Although details are different, our approach is conceptually similar to the one recently proposed by Mixa et al. [21] for simulating PbSe/PbEuTe, since it also mixes atomistic and continuum descriptions.

In our model, the growth substrate is defined as a two-dimensional squared lattice, with length parameter *a* ~4 Å, corresponding to the nearest neighboring distance on (001) planes in Si. Periodic boundary conditions are applied. On this grid, we reproduce the dynamics of single atoms and islands, assuming as mobile species only isolated atoms, in order to limit the number of possible events.

Adatoms are described as in the standard KMC: they diffuse between nearest neighboring sites through thermally activated hops, as shown in Fig. Fig.1.1. The hop-diffusion rate for an isolated adatom is determined by the Arrhenius relation

Schematic representation of the potential energy profile for adatom diffusion and detachment from islands on flat regions (**a**) and inside pits (**b**)

where ν_{0} is the hopping attempt frequency (here set to the standard value ν_{0} = 10^{13}s^{−1}), *E*_{S} is the diffusion barrier, *k* the Boltzmann constant and *T* the temperature. We use *E*_{S} = 1.1 eV, a value compatible with similar simulations in literature [22,23].

Adatoms aggregates of any dimension are defined as islands, *macroscopic* entities with volume-dependent shape and geometry: all the atomic motion inside an island is ignored assuming quick rearrangement processes to the mean equilibrium shape. An island nucleates once two adatoms reach the same lattice site forming a dimer and grows capturing adatoms diffused at its perimeter or directly deposited over its area. Coalescence between neighboring islands is possible too. Atomic detachment from islands is supposed to involve only atoms placed along the perimeter (as in Ross et al. [24]) through a thermally activated process whose barrier is assumed to be equal to the energy of an atom inside the island. Such quantity includes a chemical bonding contribution and an elastic and surface energy term. The former is taken equal to the interaction of an atom with the substrate *E*_{S} and with its *m* nearest neighbors *E*_{nn}. The latter should require in principle to solve the full elastic problem at each volume and configuration but, for the sake of simplicity, it is included in our model in an effective way as a chemical potential term μ(*V*). The effective detachment rate *R*_{det} for an island of volume *V* is then expressed as

where *N*_{p} is the number of atoms along the island perimeter. More specifically, we assume that each atom on the perimeter is bonded to three neighbors (*m* = 3, except for dimers and trimers) and for each bond we define an energy of *E*_{nn} = 0.25 eV (again, a common value from literature [22,23]). Atoms removed from an island are placed in a randomly chosen adjacent site and then treated as free adatoms. Since all events associated with islands only involve their base, the whole dynamics can be reduced to the growth plane, limiting the description to two dimensions. As a first approximation, we consider the *mean* island base as circular, with radius related to its volume. We distinguish between three typologies of structures: small 2D aggregates, {105} pyramids and domes. In the case of small atomic aggregates, for which it is not proper to consider a three-dimensional structure, we set μ(*V*) = 0 and assimilate them to cylinders one layer high, and perimeter determined by imposing conservation of the total atomic volume (sum of the atomic volumes of the atoms composing the 2D aggregate). For size greater than 3 nm^{3}, instead, we consider 3D islands. This magic value is simply determined by geometrical constraints (a {105} pyramid higher than a single (001) layer cannot occupy a smaller volume). Following the theoretical results reported in Brehm et al. [25], pyramids are transformed into domes when a critical volume of 2,400 nm^{3} is reached, so that the transformation into a steeper morphology is energetically favored. In the simulations, pyramids and domes are treated differently both from a geometrical point of view (different base corresponding to the same volume) and for what concerns the energetics. A different island chemical potential μ(*V*) (μ_{dome} < μ_{pyr} at sufficiently large *V*, reflecting the increased elastic-energy relaxation), also taken from Brehm et al. [25], is indeed attributed. Notice that in Ref. [25], the energy of the island as a function of the volume is reported for different WL thicknesses, the aim being to predict at what critical deposition islands start nucleating. Here however, for sake of simplicity, we shall assume a critical WL (~4 ML [25]) to be already formed, using only the *thick-film limit* energy values reported in the quoted reference. Figure Figure1a1a schematically summarizes our description of the dynamics on the flat surface.

The pit pattern on the surface is included in the model defining regularly spaced square areas on the reference lattice, each corresponding to a pit, associated with an extra barrier term *E*_{P}(*x*, *y*) with gaussian profile, as shown in Fig. Fig.1b.1b. In such a way, we obtain two effects: an increase in the nucleation probability into the pit, thanks to the slower adatoms diffusion, and a strong stabilization of islands grown there due to the reduced detachment rates. This is clearly an extremely simplified way for favoring nucleation in the pits which, as shown in Refs. [10,26], is driven not only by capillarity effects but also by pit-induced enhanced strain relaxation. In the simulations here below described, we considered pits 20 nm wide, using gaussians with maximum amplitude of 0.2 eV and a FWHM equal to 5 nm, and pit periods ranging from 40 to 80 nm.

In order to look for growth parameters ensuring positional and dimensional ordering for the islands, we have applied our h-KMC approach to a wide range of conditions.

In accordance with previous studies based on standard KMC approaches [19,20], our simulations show the existence of an optimal range of parameters enabling both positional and size ordering in the islands grown on the patterned substrate. Figure Figure22 shows some snapshots taken from our simulations at different temperature, deposition flux and pit spacing. From a more quantitative point of view, the size uniformity of islands inside the pits as a function of the growth temperature can be established through the distributions shown in Fig. Fig.3.3. The shown results are referred to the deposition of 1 ML of Ge. As already stressed in the previous Section, we assume that a critical (~4 ML [25]) WL is already present, so that our results are representative of a true coverage of ~5 ML.

Snapshots for the simulation after deposition of 1 ML of Ge on a pit-patterned substrate. *Red circles* represents 2D islands and *blue* ones are for pyramids. The surface shown corresponds to 480 × 480 nm^{2} (i.e., 1,200 × 1,200 lattice sites). **...**

Islands volume distributions from simulations at different temperatures (750, 850 and 950 K for the simulation parameters) at fixed flux (0.02 ML/s) and pit spacing (40 nm) after deposition of 1 ML of Ge. *Curves* are normalized and volumes are scaled with **...**

The existence of three distinct growth regimes is rather clear from a simple visual inspection of Fig. Fig.2.2. At low temperature (750 K) or high deposition flux (0.04 ML/s), there is no positional order: islands nucleate both in pits and in between, and they are non-uniform in size. Increasing the temperature or decreasing the flux, selective nucleation is achieved: mobility is now sufficient to nucleate in the energetically most favorable sites (i.e., at the center of the pits), forming stable islands and ensuring the desired positional order. If the temperature is further raised (950 K), however, communication between islands in different pits becomes important, so that some of the largers islands quickly grow at the expense of their neighbors, suppressing their growth (Ostwald ripening). As it is visible in the snapshot at high temperature in Fig. Fig.2a,2a, some pits remain almost empty so that the size control is poor. Consequently, the optimal regime for growth can be identified at intermediate temperature and deposition flux where it is possible to achieve good control both on island positioning and size: this condition corresponds to the case at intermediate temperature (850 K) shown in Figs. Figs.2a2a and and33.

Temperature and deposition flux permit to change the evolution from a regime to the other one defining the adatom effective diffusion length: increasing temperature we exponentially enhance the hopping rate, while reducing the flux we increase the time for the adatom motion ( 1/*F*). Thus, while the effect of small variations in the temperature is abrupt, changing gradually the flux, as shown in panel b) of Fig. Fig.2,2, we can slowly move from one regime to the other. In particular, from the figure we can notice that around each island there is a depleted region, the radius of which increases lowering the flux. Such an area roughly represents the island capture zone, and its extension with respect to the pit distance determines the growth modality. The higher is the temperature or the lower is the flux, the wider is the region from which a pit can capture adatoms and correspondingly the smaller is the space in between where other islands can nucleate. Approximately, when capture regions of islands in adjacent pits touch all adatoms can move to pits without nucleating outside so that the positional order is achieved. Moreover, in such a condition the region from which a pit acquires material is substantially the same for each one thus the islands grow very similar in size. Only when the growth condition gives raise to relevant overlapping in the capture areas around adjacent pits, there are competitive effects that produce ripening, as observed in the simulations at high temperature. Because the different growth regimes are determined by the ratio between the effective adatom diffusion length and pit spacing, the actual values of temperature and flux to achieve the optimal growth conditions are specific for the pattern geometry considered. This is evident comparing the two images in Fig. Fig.22 at intermediate temperature and flux but with different pit distance.

Despite the enhanced realistic conditions allowed for by our model, in terms of growth temperature and deposition flux, our simulation results still suffer from some limitations which do not allow for a direct comparison with experiments, e.g. with the results of Refs. [9,27]. The principal problem stems in the very large pit dimensions (e.g. lateral size 300 nm and depth ~50 nm), in the large period of the patterned area (~500 nm in Refs. [9,27], see [28] for experiments performed on smaller periods), in the typical volume of the islands (10^{5} nm^{3}) and in the presence of phenomena occurring at large enough volumes (Si/Ge intermixing [29], onset of plasticity [30]) which are not yet included in the presented approach. Nevertheless, the present simulations offer a clear qualitative picture of the three different growth regimes [19] characterizing growth on patterned substrates. We notice that very recent, still unpublished experimental results [27] confirmed the reported transition between random nucleation, ordered nucleation and Ostwald-ripening dominated patterns.

Interestingly, at variance with the observation of Ref. [9,27] reproduced by our model, recent experiments by Pascale et al. [12] revealed preferential nucleation sites at the pit border and not at its interior. We believe this is a strong indication of the role played by the detailed pit morphology, overlooked so far in the literature, and surely demanding for further theoretical investigations.

This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.

- Mo YW, Savage DE, Swartzentruber BS, Lagally MG. Phys. 1990. p. 1020. COI number [1:CAS:528:DyaK3cXlvVWksbs%3D]; Bibcode number [1990PhRvL..65.1020M] [PubMed] [Cross Ref]
- Eaglesham DJ, Cerullo M. Phys. 1990. p. 1943. COI number [1:CAS:528:DyaK3cXktFCqtLc%3D]; Bibcode number [1990PhRvL..64.1943E] [PubMed] [Cross Ref]
- Montalenti F, Marzegalli A, Capellini G, De Seta M, Miglio L. J. Phys.: Condens. Matter. 2007. p. 225001. Bibcode number [2007JPCM...19v5001M] [Cross Ref]
- Medeiros-Ribeiro G, Bratkovski AM, Kamins TI, Ohlberg DAA, Stanley Williams R. Science. 1998. p. 353. COI number [1:CAS:528:DyaK1cXmtlOiuw%3D%3D]; Bibcode number [1998Sci...279..353M] [PubMed] [Cross Ref]
- Ross FM, Tromp RM, Reuter MC. Science. 1999. p. 1931. COI number [1:CAS:528:DyaK1MXnvFSitLs%3D] [PubMed] [Cross Ref]
- Tersoff J, Teichert C, Lagally MG. Phys. 1996. p. 1675. COI number [1:CAS:528:DyaK28XhtlGnt7k%3D]; Bibcode number [1996PhRvL..76.1675T] [PubMed] [Cross Ref]
- Marchetti R, Montalenti F, Miglio L, Capellini G, De Seta M, Evangelisti F. Appl. 2005. p. 261919. Bibcode number [2005ApPhL..87z1919M] [Cross Ref]
- Capellini G, De Seta M, Evangelisti F, Zinovyev VA, Vastola G, Montalenti F, Miglio L. Phys. 2006. p. 106102. COI number [1:STN:280:DC%2BD283gtVGgtw%3D%3D]; Bibcode number [2006PhRvL..96j6102C] [PubMed] [Cross Ref]
- Zhong Z, Bauer G. Appl. 2004. p. 1922. COI number [1:CAS:528:DC%2BD2cXitFSnsLY%3D]; Bibcode number [2004ApPhL..84.1922Z] [Cross Ref]
- Zhong Z, Schwinger W, Schffäler F, Bauer G, Vastola G, Montalenti F, Miglio L. Phys. 2007. p. 176102. Bibcode number [2007PhRvL..98q6102Z] [Cross Ref]
- Zhang JJ, Stoffel M, Rastelli A, Schmidt OG, Jovanović V, Nanver LK, Bauer G. Appl. 2007. p. 173115. Bibcode number [2007ApPhL..91q3115Z] [Cross Ref]
- Pascale A, Berbezier I, Ronda A, Kelires PC. Phys. Rev. B. 2008. p. 075311. Bibcode number [2008PhRvB..77g5311P] [Cross Ref]
- Rastelli A, von Känel H, Albini G, Raiteri P, Migas DB, Miglio L. Phys. 2003. p. 216104. COI number [1:STN:280:DC%2BD3s3mtlalsA%3D%3D]; Bibcode number [2003PhRvL..90u6104R] [PubMed] [Cross Ref]
- Holý V, Stangl J, Zerlauth S, Bauer G, Darowski N, Lübbert D, Pietsch U. J. 1999. p. A234. Bibcode number [1999JPhD...32A.234H] [Cross Ref]
- Weeks JD, Gilmer GH. Adv. 1979. p. 157. COI number [1:CAS:528:DyaL3cXos1SgtQ%3D%3D] [Cross Ref]
- Bortz AB, Kalos MH, Lebowitz JL. J. 1975. p. 10. Bibcode number [1975JCoPh..17...10B] [Cross Ref]
- Voter A. Phys. Rev. B. 1986. p. 6819. COI number [1:CAS:528:DyaL2sXitlOlsQ%3D%3D]; Bibcode number [1986PhRvB..34.6819V] [PubMed] [Cross Ref]
- Fichthorn KA, Weinberg WH. J. 1991. p. 1090. COI number [1:CAS:528:DyaK3MXlt1ShtLk%3D]; Bibcode number [1991JChPh..95.1090F] [Cross Ref]
- Lee C, Barabasi A-L. Appl. 1998. p. 2651. COI number [1:CAS:528:DyaK1cXmvVSitbc%3D]; Bibcode number [1998ApPhL..73.2651L] [Cross Ref]
- Nurminen L, Kuronen A, Kaski K. Phys. Rev. B. 2000. p. 035407. Bibcode number [2001PhRvB..63c5407N] [Cross Ref]
- Mixa M, Holý V, Springholz G, Bauer G. Phys. Rev. B. 2009. p. 045325. Bibcode number [2009PhRvB..80d5325M] [Cross Ref]
- Mysliveček J, Schelling C, Schäffler F, Springholz G, Šmilauer P, Krug J, Voigtländer B. Surf. 2002. p. 193. Bibcode number [2002SurSc.520..193M] [Cross Ref]
- Kawamura T, Natori T. Surf. 1999. p. 148. COI number [1:CAS:528:DyaK1MXmsFWluro%3D]; Bibcode number [1999SurSc.438..148K] [Cross Ref]
- Ross FM, Tersoff J, Tromp RM. Phys. 1998. p. 984. COI number [1:CAS:528:DyaK1cXpsFaqtw%3D%3D]; Bibcode number [1998PhRvL..80..984R] [Cross Ref]
- Brehm M, Montalenti F, Grydlik M, Vastola G, Lichtenberger H, Hrauda N, Beck MJ, Fromherz T, Schäffler F, Miglio L, Bauer G. Phys. Rev. B. 2009. p. 205321. Bibcode number [2009PhRvB..80t5321B] [Cross Ref]
- Schülli TU, Vastola G, Richard M-I, Malachias A, Renaud G, Uhlík F, Montalenti F, Chen G, Miglio L, Schäffler F, Bauer G. Phys. 2009. p. 025502. Bibcode number [2009PhRvL.102b5502S] [PubMed] [Cross Ref]
- Grydlik M, Brehm M, Fromherz T, Bauer G. private communications.
- Dais C, Mussler G, Sigg H, Fromherz T, Auzelyte V, Solak HH, Grützmacher D. Europhys. 2008. p. 67017. Bibcode number [2008EL.....8467017D] [Cross Ref]
- De Seta M, Capellini G, Evangelisti F. J. 2002. p. 614. COI number [1:CAS:528:DC%2BD38Xkslersbs%3D]; Bibcode number [2002JAP....92..614D] [Cross Ref]
- Marzegalli A, Zinovyev VA, Montalenti F, Rastelli A, Stoffel M, Merdzhanova T, Schmidt OG, Miglio L. Phys. 2007. p. 235505. COI number [1:STN:280:DC%2BD1c%2FnslSitQ%3D%3D]; Bibcode number [2007PhRvL..99w5505M] [PubMed] [Cross Ref]

Articles from Nanoscale Research Letters are provided here courtesy of **Springer**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |