Xin
Ou
*ab,
Karl-Heinz
Heinig
a,
René
Hübner
a,
Jörg
Grenzer
a,
Xi
Wang
b,
Manfred
Helm
ac,
Jürgen
Fassbender
ac and
Stefan
Facsko
*a
aHelmholtz-Zentrum Dresden-Rossendorf, Institute of Ion Beam Physics and Materials Research, Bautzner Landstr. 400, 01328 Dresden, Germany. E-mail: s.facsko@hzdr.de
bState Key Laboratory of Functional Material for Informatics, Shanghai Institute of Microsystem and Information Technology, Chinese Academy of Sciences, Shanghai 200050, China. E-mail: ouxin@mail.sim.ac.cn
cTechnische Universität Dresden, 01062 Dresden, Germany
First published on 25th August 2015
Semiconductor quantum dots and wires are important building blocks for future electronic and optoelectronic devices. The common way of producing semiconductor nanostructures is by molecular beam epitaxy (MBE). In this additive growth process atoms are deposited onto crystalline surfaces and self-assemble into 3D structures. Here we present a subtractive process, in which surface vacancies are created by ion impacts. On terraces of crystalline surfaces their nucleation forms depressions which coarsen and finally lead to a self-organized 3D morphology. It is shown that this kind of spontaneous pattern formation is inherent to the ion induced erosion process on crystalline surfaces and is analogous to 3D growth by MBE. However, novel facets are found due to slightly different energetics and kinetics of ad-atoms and surface vacancies, especially at Ehrlich–Schwoebel step-edge barriers. Depending on the crystal orientation, three-fold, four-fold, six-fold symmetry, as well as extremely regular periodic nanogrooves can be produced on different orientations of group IV (Si, Ge) and III–V (GaAs, InAs) semiconductors.
Ion irradiation, on the other hand, is frequently used to etch materials by sputtering7 and can be easily scaled up to 300 mm wafer size. It can also be employed for nanopatterning surfaces by a self-organized, bottom-up approach.8 The resulting patterns are periodic ripple structures with periodicities down to 15 nm9 and can achieve quite high order as in the case of Si and Ge irradiated with 26 keV Au− ions.10 Such ripple patterns can be applied as templates for growing nanostructured thin films with strong anisotropic magnetic or optical properties.11 Therefore, a broad interest exists for this simple, high-throughput and inexpensive technology for nanopatterning surfaces. The main drawback of this method is that semiconductors are amorphized when ion irradiation is performed at room temperature.12 At elevated temperature, however, ion induced interstitials and vacancies are mobile enough to recombine or to reach the surface before a new ion hits the same region. Thus ion irradiation performed at temperatures higher than the recrystallization temperature prevents amorphization and leads to an additional surface instability due to the ES barrier for surface vacancies.13–20 Accordingly, in analogy to the case under growth conditions, this kind of surface instability leads to the formation of 3D nanostructures (Fig. 1a). Vacancies which are created by the ion impacts during irradiation, nucleate and coarsen to form pits that grow inside the surface.13 The formation of such structures has firstly been observed on metal surfaces, which remain crystalline even at room temperature, due to their high diffusion and non-covalent bonding.14–16
Recently, we have observed regular, dense patterns of inverse pyramids on Ge(001) surfaces irradiated by high fluence 1 keV Ar+ at temperature above 250 °C.21 These patterns exhibit {105} facets and strongly resemble mound patterns grown by MBE if they are reversed. Here, we present the formation of nanogrooves on the (001) surface of group III–V semiconductors (GaAs, InAs) with a zinc-blende structure. Although the symmetry of the surface is four-fold, like the (001) surface of Si and Ge, a symmetry-breaking due to the stacking of alternating planes of Ga (In) and As leads to an almost perfect two-fold symmetry of the resulting patterns even under normal incidence. In combination with enhanced ion-induced surface diffusion this symmetry-breaking driving force leads to faceted periodic nanogroove patterns oriented along the [10] direction with the highest degree of regularity observed so far. Similar ripple patterns have been observed recently by irradiation of GaAs with hyperthermal ion beams (30 eV Ar+), however with much lower regularity.22 Furthermore, on the (111) surfaces of Si and Ge with diamond structure irradiated with high fluence we found peculiar {123} and {356} facets that have not been observed so far. All of these ion induced patterns can be described by a universal generalized continuum equation21,23 by defining a non-equilibrium surface free energy with the required symmetry, giving the formation dynamics with the correct roughening and coarsening behavior.
Fig. 2a–f present typical patterns on different semiconductor surfaces after ion irradiation performed at normal incidence above the recrystallization temperature. Faceted nanostructures with different orientations and shapes develop, which are, with respect to the surface, reversed to similar structures found by MBE.24 The periodicity and the shape of the faceted structures can be tuned by the irradiation conditions (temperature and fluence) as well as by using different orientations of the crystalline surfaces, i.e. low-index crystal planes or different miscut angles from these planes. The formation of facets can be also concluded from the corresponding two-dimensional (2D) angle distributions25 (left) and three-dimensional (3D) zoom-in images (right) shown above the AFM images. The symmetry of the patterns reflects the crystal symmetry of the surfaces. For the (001) surface of GaAs (Fig. 2a) and InAs (Fig. 2d) the evolving patterns are periodic, faceted grooves oriented along the [10] direction. This is quite surprising, as the irradiation is performed at normal incidence, thus no asymmetry is imposed by the ion beam. As shown in Fig. 1b, the direction, where dimer rows of Ga (In) are formed due to a (2 × 1) surface reconstruction parallel to the steps,26 is energetically preferred over the formation of As dimer rows. Furthermore, the [10] direction is also the easy diffusion direction of the Ga (In) atoms.27 The 2D angle distribution shows two peaks at ±16° and ±19°, respectively. These polar angles can be identified with the {115} (15.79°) and {114} (19.47°) crystal planes, respectively (Fig. 2g).
For Si(001) and Ge(001) a 4-fold symmetry can be identified in Fig. 2b and e. For single-element semiconductors the dimer rows in the [10] and [110] directions consist of the same element, i.e. there is no longer a difference in their surface energy, and, consequently, the symmetry breaking disappears. The angle distribution exhibits peaks at 12° for Ge and 16° for Si, which are close to the {105} (11.31°) and {115} (15.79°) planes of the diamond lattice, respectively. The polar plot of the orientation of these facets with respect to the (001) surface is shown in Fig. 2h. The azimuthal orientation of the structures is different for Si(001) and Ge(001). On Si(001) the structures are oriented along the [110] direction, whereas on Ge(001) the structures align along the [100] direction. This is quite peculiar, because Si and Ge have the same crystal structure and exhibit usually the same (2 × 1) surface reconstruction of the (001) surface. The same facets and the rotation of the pattern by 45° between Si and Ge appear in homoepitaxial growth of Si and Ge as well.28–31 The emergence of these edges is attributed to an additional barrier at kink sites, similar to the barrier at step edges. However, the origin of the different orientation of the Si and Ge mound facets has not yet been conclusively clarified.
The patterns on the (111) surface of Si (Fig. 2c) and Ge (Fig. 2f) exhibit a different symmetry. On the Ge(111) surface, inverse pyramidal structures with an isotropic distribution are formed. The corresponding 2D angle distribution reveals a three-fold symmetry in the formation of the facets with a polar angle θ of 13°–16°. Noticeable is the development of a double peak structure in the 2D angle distribution separated by an azimuthal angle of ∼38°. Inspection of the geometrical projection of the possible facets of the diamond lattice with respect to the (111) surface reveals that the measured angle distribution coincides with the {356} facets. As can be seen in the polar plot of the crystal planes in Fig. 2i, six of these facets have the right polar angle of 15°. They appear as three pairs of two poles with a difference in the azimuthal angle of 38.2° and rotated by 120° to each other, in perfect agreement with the experimentally determined facet angles. These facets are typically not observed in homoepitaxial or heteroepitaxial growth of self-assembled Ge nanostructures. Thus, we can conclude that they are ion induced non-equilibrium stable crystal facets.
On the Si(111) surface inverse hexagonal pyramids are formed during ion irradiation. The structure size of these patterns exhibits a much broader distribution than for the patterns on Ge(111). Thus for Si(111) the wavelength selection at the beginning of the pattern formation is much weaker than for Ge(111). Furthermore, the 2D angle distribution shows six peaks with a polar angle of 21°–23° with a hexagonal symmetry. The facets, which correspond to these values are the {123} crystal planes (22.21°) of the diamond lattice, as can be seen in the geometrical projection of the planes with respect to the (111) surface in Fig. 2i.
Due to the fact that the facets are crystal planes, asymmetric facets can also be formed on surfaces with predefined miscut angles. Fig. 3 shows nanogroove patterns on vicinal GaAs(001) surfaces with miscut angles of 6° and 10° in the [110] direction, respectively. The wavelength of the formed nanogrooves increases from 46 nm on GaAs(001) surfaces without miscut (Fig. 1c) to 65 nm (Fig. 3a) and 90 nm (Fig. 3b) for the surface with 6° and 10° miscut angles, respectively. Fig. 3c and d show the cross-sectional TEM of GaAs(001) with a 10° miscut angle. The shape of these nanogrooves becomes asymmetrical with absolute facet angles to the (100) planes of 20° and 21°, respectively. To keep the same angle of the {114} facet with respect to the (001) surface on both sides one sidewall of the groove is elongated, thus increasing the periodicity of the pattern.32 These kinds of structures are particularly interesting to be applied as optical grating devices for ultraviolet light.33
The formation of 3D nanostructures on ion irradiated surfaces at temperatures above the recrystallization temperature can be regarded as the reverse mechanisms to 3D growth in homoepitaxy.16,21,24,34,35 The continuous ion irradiation creates a surface vacancy density, which is much higher than the vacancy density under equilibrium conditions, thus driving additional non-equilibrium surface currents. Due to the step edge barriers, the vacancy diffusion currents are “downhill”, leading to a surface instability.6,23,36 The orientation of the patterns along the crystalline directions and the selection of specific facets on the ion irradiated surfaces indicate that the diffusion currents are dominant and sputtering induced effects, like slope/curvature dependent erosion rates and mass redistribution, can be neglected. Thus the continuum equations, which are typically used to describe the morphology evolution on amorphized surfaces under ion irradiation, do not account for the crystal structure and cannot be used here. On binary materials an additional instability is expected37 due to preferential sputtering of one element changing the stoichiometry of the surface. The coupling of this altered surface layer and the topography can lead to the formation of periodic patterns as well. However, the crystalline GaAs (InAs) are Ga (In) terminated and no preferential sputtering or segregation of one element has been observed.38 Furthermore, even at temperatures below the recrystallization temperature we did not observe any pattern formation on GaAs or InAs, thus on these amorphized surfaces ion induced smoothing dominates at normal incidence.
(1) |
The surface current j includes three contributions: jHM = κ∇(∇2h) is the Herring–Mullins surface diffusion, jKPZ = σ∇(∇h)2 is a nonlinear current (conserved Kardar–Parisi–Zhang term) that breaks the up–down symmetry of the resulting surface, and jNE is the non-equilibrium surface current that includes the instability due to the ES step edge barrier and the formation of facets. Furthermore, a white noise term η(x, t) is included to account for the stochastic nature of the vacancy formation.
Similar to the chemical potential in equilibrium, which determines the vacancy and ad-atom concentrations in thermal equilibrium, a non-equilibrium potential UNE can be defined corresponding to the non-thermal ion induced concentrations of vacancies and ad-atoms. UNE(m) is proportional to the surface density of atoms and is a function of the local slope m = ∇h. The minimization of the non-equilibrium potential with respect to the surface slope induces a non-equilibrium surface current . Due to the high ion induced density of vacancies and their restricted interlayer transport due to the ES barrier the original surface becomes unstable, i.e. UNE(m) has a local maximum at zero slope m = (0,0). The interplay of this instability and the surface diffusion leads to a wavelength selection at the beginning of ion irradiation. At the minima of UNE(m) the surface current jNE vanishes leading to a selection of the preferred slopes (facets) at later times.
The different surface patterns in Fig. 4 result from appropriate choices of the non-equilibrium potential UNE(m), reflecting the symmetry of the crystalline surface and possessing minima at the respective, experimentally observed facet angles. In Fig. 4 the chosen non-equilibrium potential and results of the numerical integration of the continuum equation with the corresponding 2D angle distributions are presented. The good agreement between the experimentally observed pattern in Fig. 2 and results of the continuum equation in Fig. 4 indicates that the facet formation is strongly linked to surface currents driven by the minimization of the non-equilibrium effective free energy. In analogy to the determination of the chemical potential by observing the shape of crystals grown under equilibrium conditions, we can thus assess the shape and the first minima of the non-equilibrium potential UNE(m) in the vicinity of the original surface exposed to ion irradiation. It is known that under continuous vacancy creation crystal facets can become unstable and non-equilibrium steady state facets of different orientation can be formed.39
Materials | Ion species | Ion energy (eV) | Temperature window (°C) |
---|---|---|---|
Si(100) | Xe+ | 500 | 480–580 |
Ge(100) | Ar+ | 1000 | 250–430 |
Ge(111) | Ar+ | 1000 | 250–380 |
InAs(100) | Ar+ | 1000 | 180–430 |
GaAs(100) | Ar+ | 1000 | 200–480 |
In every integration step white noise with an amplitude of 0.05 has been added. The coefficients used for the numerical integrations for the different cases in Fig. 4 are as follows: (a) ε = 1, κ = 1, σ = −1, (b) ε = 1, κ = 4, σ = −1, (c) ε = 2, κ = 4, σ = −2, and (d) ε = 2, κ = 4, σ = −2. The coefficients for the numerical integration have been chosen solely by the intention to produce patterns with a structure density similar to the AFM images in Fig. 2.
This journal is © The Royal Society of Chemistry 2015 |