Modeling of the growth of GaAs–AlGaAs core–shell nanowires

  1. 1 ,
  2. 1,2 and
  3. 1
1Department of Engineering Sciences and Applied Mathematics, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208-3125, USA
2Department of Materials Science and Engineering, Northwestern University, 2225 Campus Drive, Evanston, Illinois 60208-3030, USA
  1. Corresponding author email
Guest Editor: A. S. Sidorenko
Beilstein J. Nanotechnol. 2017, 8, 506–513. https://doi.org/10.3762/bjnano.8.54
Received 11 Aug 2016, Accepted 10 Feb 2017, Published 24 Feb 2017
Letter
cc by logo
Album

Abstract

Heterostructured GaAs–AlGaAs core–shell nanowires with have attracted much attention because of their significant advantages and great potential for creating high performance nanophotonics and nanoelectronics. The spontaneous formation of Al-rich stripes along certain crystallographic directions and quantum dots near the apexes of the shell are observed in AlGaAs shells. Controlling the formation of these core–shell heterostructures remains challenging. A two-dimensional model valid on the wire cross section, that accounts for capillarity in the faceted surface limit and deposition has been developed for the evolution of the shell morphology and concentration in AlxGa1−xAs alloys. The model includes a completely faceted shell–vapor interface. The objective is to understand the mechanisms of the formation of the radial heterostructures (Al-rich stripes and Al-poor quantum dots) in the nanowire shell. There are two issues that need to be understood. One is the mechanism responsible for the morphological evolution of the shells. Analysis and simulation results suggest that deposition introduces facets not present on the equilibrium Wulff shapes. A balance between diffusion and deposition yields the small facets with sizes varying slowly over time, which yield stripe structures, whereas deposition-dominated growth can lead to quantum-dot structures observed in experiments. There is no self-limiting facet size in this case. The other issue is the mechanism responsible for the segregation of Al atoms in the shells. It is found that the mobility difference of the atoms on the {112} and {110} facets together determine the non-uniform concentration of the atoms in the shell. In particular, even though the mobility of Al on {110} facets is smaller than that of Ga, Al-rich stripes are predicted to form along the {112} facets when the difference of the mobilities of Al and Ga atoms is sufficiently large on {112} facets. As the size of the shell increases, deposition becomes more important. The Al-poor dots are obtained at the apices of {112} facets, if the attachment rate of Al atoms is smaller there.

Findings

Core–shell nanowires with heterostructures hold great promise in photonic and electronic applications because of their high sensitivity to electronic and magnetic fields. However, controlling the formation of these heterostructures remains a challenge because they are typically embedded in 3D matrices. One of the potential solutions to this problem is to create heterostructures near the edges of the nanowires. Classical approaches to create this kind of heterostructure are to deposit materials with lattice parameters different from that of the wires, such as growing a Ge shell on a Si nanowire [1], so that the quantum dots form as a result of stresses generated by the lattice parameter mismatch. Recently, GaAs has been used as substrates to deposit AlGaAs on their sides. Heiss et al. [2] observed that the deposition of the AlGaAs layer leads to Al-poor quantum dots located within the fork-like Al-rich stripes along the [Graphic 1] directions (Figure 1).

[2190-4286-8-54-1]

Figure 1: Schematic of a core–shell nanowire (left) with possible configurations of the cross section (right) observed in experiments. These are the six-fold symmetry configuration observed in [3] with Al-rich {112} facets along the corners of the hexagonal core (top right); the three-fold symmetry configuration observed in [4] (middle right) and the Al-poor quantum dots (red) at the apices adjunct with the Al-rich stripes observed in [2] (bottom right).

Yet, there is a very small elastic mismatch between AlGaAs shell and GaAs core. So, the question is, how do these heterostructures form in this experiment? Moreover, many other experiments also observe radial heterostructures with the same materials [3]. Jiang et al. [5] also observed Al-rich stripes along the [Graphic 2] directions. Zhang et al. [4] obtained a three-fold symmetry structure with Al-rich stripes (Figure 1). Thus, an understanding of the mechanisms leading to these heterostructures could be very helpful in controlling the formation of these core–shell nanowires.

There are two issues that need to be understood. One is the mechanism responsible for the morphological evolution of the shells. The surface energy density on the {112} facets is at least three times larger than that on the {110} facets. Energy minimization thus tells us that there are no {112} facets in the equilibrium state. Thus, the question is, how do the {112} facets, which are not present on the Wulff shape, arise in the configuration of the shell of nanowires? The other issue is the mechanism responsible for the segregation of Al atoms in the shells. More specifically, under certain deposition conditions, Al atoms move slower than Ga atoms along the {110} facets [6]. If the atoms move from {110} to {112} facets, then the Al atoms should be left behind on the {110} facets rather than accumulate in the [Graphic 3] directions. Thus, the question is how the Al-rich stripes could form along {112} facets when the atoms move from {110} facets to {112} facets and the mobility of Al is lower than that of Ga on {110} facets [6]. Following this question is that of how the Al-poor dots form directly following the Al-rich stripes.

These problems are solved in [7] and [8], respectively. For simplicity, at first a pure material is deposited around the hexagonal core to investigate the mechanisms governing the morphological evolution. In [7], we proposed a mathematical model for the deposition of pure materials onto a faceted nanowire. In the model, the growth of each facet is due to the deposition of atoms from the vapor and surface–atom diffusion from its neighbouring facets (Figure 2).

[2190-4286-8-54-2]

Figure 2: Schematic of cross section of a core of a nanowire that is a six-fold symmetric structure with {112} facets along the corners of the hexagonal core. The lengths of the facets are l{112} and l{110}, respectively. In this structure, the angle θ between a {110} facet and its neighboring {112} facet is π/6.

It is suggested in experiments that bulk diffusion is negligible relative to surface diffusion for typical temperatures and length scales. Moreover, the motion of the atoms along the surface due to surface diffusion is driven by the difference of the chemical potentials between neighbouring facets. The average chemical potential on the i-th facet is calculated by

[2190-4286-8-54-i1]
(1)

where ni is the unit normal vector of the i-th facet and γ(ni) is the surface energy density on the i-th facet. Thus, if there is deposition, the mass conservation around the surface gives the following equations:

[2190-4286-8-54-i2]
(2)
[2190-4286-8-54-i3]
(3)
[2190-4286-8-54-i4]
(4)

where ui is the normal velocity of the i-th facet, [Graphic 4] is the surface diffusion flux along i-th facet and μi is the chemical potential on the i-th facet, Ω is the atomic volume and Fi is the total deposition rate of the material in the outer normal direction of the i-th facet, Mi is the mobility of the atoms on the i-th facet, and li is the length of the i-th facet. Equation 2Equation 4 together with the continuity conditions of surface diffusion flux and chemical potential at the corners give a complete equation system of the evolution of the configuration of a crystal. Thus, facet lengths evolve as

[2190-4286-8-54-i5]
(5)
[2190-4286-8-54-i6]
(6)

where ΔF{112} = −2F{112}cotθ + 2F{112}cscθ, ΔF{110} = 2F{112}cscθ − 2F{110}cotθ (θ is shown in Figure 2 and θ = π/6 in this six-fold symmetric configuration). Here, v{112} and v{110} are the normal velocities of the corresponding facets due to diffusion of the surface atoms. This fully faceted model without deposition flux was first proposed by Carter and co-workers [9]. The derivation for the continuum model can be found in [7,10-12]. Different from the focus of [9], in [7], we did a systematical theoretical analysis and numerical simulations for the case shown in Figure 2.

Here, in order to determine the influence of diffusion and deposition on the morphological evolution of the shell, consider first a surface evolution due only to surface diffusion, i.e., without deposition. Figure 3 shows the evolution of a dodecagonal crystal with initially l{112} = l{110} (red lines) to a hexagonal crystal due to surface diffusion. Under the condition γ{112} = 3γ{110}, Equation 1 implies that [Graphic 5]. Thus, surface diffusion tends to move the atoms from {110} to {112} facets and eventually the {112} facets disappear (black line). In addition, through this numerical example, it is seen that the dynamic simulation result yields the equilibrium Wulff shape.

[2190-4286-8-54-3]

Figure 3: Evolution of a dodecagonal crystal (bounded by red lines) to a hexagonal crystal (bounded by black lines) due to surface diffusion.

Now, consider deposition in the process of growth of the shell with the same conditions as the last numerical simulation, i.e., F{112} = 1.1F{110} with F{110} being large enough. These deposition rates can be different on facets of different orientations due to different sticking coefficients or exchange rates between the surface and bulk. Instead of getting a hexagonal shape, a dodecagonal shell is obtained (Figure 4a). Moreover, as the thickness of the shell gets large, the ratio between the lengths of the facets approaches a non-zero constant, i.e., l{112}/l{110} ≈ ΔF{112}F{110} = 0.202 (Figure 4b). This is because diffusion becomes negligible compared with deposition as the size of the shell gets large (see Equation 6). Moreover, when the thickness of the shell is large enough, this ratio does not change, even though the size of the facets continue to grow; this is a self-similar shape, the “kinetic Wulff shape” because it is obtained when kinetic effects dominate. Thus, it is shown that deposition allows for the possibility of {112} facets that are not present on the equilibrium Wulff shape.

[2190-4286-8-54-4]

Figure 4: (a) Evolution of a dodecagonal shell to its kinetic Wulff shape due to surface diffusion and deposition. (red: six {110} facets; blue: six {112} facets.) (b) Evolution of the ratio between the lengths of the facets shown in (a).

According to the discussions above, surface diffusion tends to remove the {112} facets whereas deposition tends to preserve them. Figure 5a shows the balance between the surface diffusion and deposition giving a slowly varying facet size at early times. The numerical simulation is compared with the experiment result in [3] (Figure 5b) giving almost quantitative agreement. When a certain thickness of the shell is reached, the size of the {112} facets increases rapidly (Figure 6). As shown in the inset in Figure 6b, the length of {112} facets evolves with the increasing of the thickness of the shell, even though the rate of change is much smaller than that of {110} facets at early time. This implies that the length of the {112} facets is never constant over time.

[2190-4286-8-54-5]

Figure 5: (a) Cross section of the core–shell nanowire with six {112} facets (blue color) along the corners of six {110} facets (red color) in the numerical simulation; the stripe structure along the [Graphic 6] directions within 25 nm of thickness of the shell. (b) The experimental result showing Al-rich stripes along the [Graphic 7] directions (dark lines) obtained in [3].

[2190-4286-8-54-6]

Figure 6: (a) Cross section of the core–shell nanowire with six {112} facets (blue color) along the corners of six {110} facets (red color) in the numerical simulation when the shell of the nanowire is large enough. (b) Evolution of the length of the facets with respect to the thickness of the shell.

This conclusion departs from the traditional understanding on the formation of the Al-rich stripes along the [Graphic 8] directions. In [2,3,13], the claim is that the Al-rich stripes form because the length of the {112} facets reaches its “self-limiting” size and does not change further upon growth of the shell. Such self-limiting facets are absent in our results. (In order to see that such a self-limiting facet size does not exist, we argue by contradiction. Suppose there was a self-limiting facet size for the {112} facets, as the size of the shell of the nanowire increases, the length of {110} facets must then increase. Thus, the ratio of the length of the {112} facets to the length of the {110} facets cannot be a constant. This contradicts the attainment a self-similar shape, i.e., the kinetic Wulff shape in which the ratio of the length of {112} facets to the length of {110} facets is a constant, when the size of the nanowire is large enough.) Rigorous proof of the non-existence of a self-limiting facet size can be found in [7].

Given that there is no self-limiting size for any facet, another set of deposition rates can yield the dot configuration that is formed spontaneously by a transient growth (Figure 7b). The configuration obtained in the numerical simulation is very close to that seen experimentally (Figure 7a).

[2190-4286-8-54-7]

Figure 7: (a) This is partial view of a cross section of the shell of core–shell nanowire with one {112} facet along the corner of two neighboring {110} facets. (b) Detail of a Al-poor quantum dot located within the fork-like Al-rich stripes obtained in [2].

Besides the dot configuration shown in Figure 7b, Heiss et al. [2] also observed the segregation of Al in the shell of the nanowire (see the Al concentration shown in Figure 7b). It is obvious that a model assuming pure materials cannot address the mechanisms leading to the non-homogeneous distribution of Al in the shell. A two-dimensional fully faceted model [8] was developed for the growth of the shell with an AxB1−x (AlxGa1−xAs) alloy on a hexagonal core. Surface morphology and composition on the i-th facet evolve as

[2190-4286-8-54-i7]
(7)
[2190-4286-8-54-i8]
(8)

where δ is the ratio of the density of lattice sites on the surface to that in the bulk. For each component m, ξm is the local composition fraction of atoms m of the surface layer and

[Graphic 9]

Moreover, ξA + ξB = 1. [Graphic 10] is the surface-diffusion flux of component m. μm is the chemical potential at the surface. Msm is the surface diffusivity of atoms m. Fn is the total deposition rate of atoms and is positive for growth, and xm is the fraction of component m on the surface due to the deposition flux. Cm is the local composition of the “bulk” material immediately beneath the surface layer. This model suggests that the changing rate of the concentration of A atoms is due to the surface diffusion of the atoms along the surface, the physical deposition of the atoms from the deposition fluxes, and the motion of the surface, which leads to the atoms leaving or entering the surface from below. Similar to the pure-material model, the chemical potentials can be related to the surface energy density (γ) and bulk energy density (gα) in its average sense, as follows:

[2190-4286-8-54-i9]
(9)
[2190-4286-8-54-i10]
(10)

where

[Graphic 11]

is the average surface energy density of the i-th facet. In addition, the advancing surface leaves behind material with composition CA determined by the equilibrium with the surface layer 1/δ. Equation 7Equation 10 cooperating with the continuity conditions at the joint points of two neighboring facets for the surface diffusion fluxes, the chemical potentials at the surface, and the surface concentrations give a complete system for the evolution of morphology and surface concentration. A derivation of the model on faceted surfaces with anisotropic surface energy density and numerical simulation details can be found in [8]. Similar models on smooth surfaces with isotropic surface energy density were proposed in [14,15].

The mobility of Ga atoms on the {110} facets is difficult to measure. However, the measurement in [6] gives a diffusion of Ga that is ten times larger than that of Al atoms, i.e., [Graphic 12]. The diffusivities of Al and Ga are likely functions of surface reconstructions and thus can be functions of the deposition conditions. For simplicity, take

[Graphic 13]

and

[Graphic 14]

where γ0(n) is the surface energy density for pure atoms A or atoms B, kB is Boltzmann’s constant, T is the absolute temperature, Γ0 and ρ0 are the lattice sites densities on the surface and in the bulk, respectively. In the case considered here un > 0, the advancing surface leaves behind material with composition CA determined by the equilibrium with the surface layer, (1/δ)∂γ/∂ξA = dgα/dCA[8,14,16,17]. In experiments, [Graphic 15] and γ0 ≈ 1 J/m2, Γ0 ≈ 1018/m2, ρ0 ≈ 1028/m3 and T ≈ 1000 K. Usually, the alloy used in the deposition flux in the experiments is Al0.33Ga0.67As, i.e., xA = 0.33 in Equation 7.

In the numerical simulation, suppose the Al atoms diffuse much faster than Ga atoms on the {112} facets. For example, [Graphic 16]. This is exaggerated, but it serves to illustrate the basic concepts that lead to stripes and dots. Under these conditions, it is shown first that both Al and Ga atoms tend to move from {110} to {112} facets because [Graphic 17] (Figure 8a). Second, Figure 8b shows that the portion of Al atoms subjected to the tangential diffusion fluxes from {112} to {110} facets along the surface at the corners of {112} and {110} facets (JsA/(JsA + JsB)) are much larger than the portion of Al atoms subjected to the deposition flux that introduces the accumulation of Al atoms along the {112} facets (Figure 9). Moreover, this effect is enhanced within a certain thickness of the shell because surface diffusion dominates over deposition at early times. As time proceeds, the thickness of the shell of the nanowire gets large, diffusion gets weaker with a decrease of [Graphic 18]. Thus, JsA and JsB decrease (Figure 8a,c). The dynamic growth process switches from diffusion-dominated to deposition-dominated. This is why the concentration of Al atoms along the surface tends to be a constant equal to the concentration of Al atoms in the deposition flux when the thickness of the shell is large enough (Figure 9b).

[2190-4286-8-54-8]

Figure 8: (a) Evolution of the difference of the average chemical potentials of the facets with the thickness of the shell(× 2.5 nm). (b) Evolution of the portion of Al atoms taken in the diffusion flux along the surface. (c) Evolution of the diffusion fluxes of Al atoms (denoted by A) and Ga atoms (denoted by B) from {110} facets to {112} facets.

[2190-4286-8-54-9]

Figure 9: (a) Distribution of the concentration of Al in the shell of the nanowire. (b) Evolution of the average concentration of Al on different facets versus the thickness of the shell (in the direction normal to {110} facets).

If it is harder for Al atoms to attach on the {112} than on the {110} facets from the deposition flux (for example, [Graphic 19] and [Graphic 20] = 0.297, [Graphic 21] = 0.330) then the concentration of Al atoms along each facet tends to be a constant that is equal to its attachment rate on the corresponding facets in the deposition-dominated regime (Figure 10). The behaviour of the segregation is similar to the results of the former numerical simulation in the surface-diffusion-dominated regime (in this case [Graphic 22]). This explains why the Al-poor quantum dots can grow adjacent to an Al-rich stripe in some experiments [2].

[2190-4286-8-54-10]

Figure 10: (a) Distribution of the concentration of Al in the shell of the nanowire. (b) Evolution of the average concentration of Al on different facets with the thickness of the shell (in the direction normal to {110} facets).

The presence of stripes and dots also strongly depends on the diffusivities of Al and Ga on the {110} facets. As the ratio between the diffusion coefficients on {110} facets approach unity, the ratio between the diffusion coefficients of Al between {110} and {112} facets necessary to get stripes and dots becomes much smaller than 1000. In particular, according to the results shown in Figure 11, a change from [Graphic 23] = 0.1 to [Graphic 24] = 0.5 will lead to a decrease of the critical value of [Graphic 25] to yield Al-rich stripes along the [Graphic 26] directions down to about a tenth.

[2190-4286-8-54-11]

Figure 11: Influence of the ratio between the mobilities of different atoms and deposition rates on the surface concentration segregation. “+” means [Graphic 27] whereas “−” means [Graphic 28].

It is worth mentioning that even though the morphological evolution is coupled with the mass transfer, as shown in Equation 7 and Equation 8, the surface segregation does not influence the morphological evolution significantly. This is because γ0 ≈ 1 J/m2 and kBΓ0T ≈ 0.01 J/m2. Hence, [Graphic 29], which implies that [Graphic 30] is essentially determined by the surface energy densities of the pure materials on the facets. Therefore, the mechanisms for the morphological evolution in pure materials [7] are still applicable in the case of alloys.

Acknowledgements

The authors would like to thank Jean-Noël Aqua, Lincoln J. Lauhon and Anna Fontcuberta i Morral for helpful discussions. This research was supported by the Office of Naval Research under Grant no. N00014-14-1-0697.

References

  1. Pan, L.; Lew, K.-K.; Redwing, J. M.; Dickey, E. C. Nano Lett. 2005, 5, 1081–1085. doi:10.1021/nl050605z
    Return to citation in text: [1]
  2. Heiss, M.; Fontana, Y.; Gustafsson, A.; Wüst, G.; Magen, C.; O’Regan, D. D.; Luo, J. W.; Ketterer, B.; Conesa-Boj, S.; Kuhlmann, A. V.; Houel, J.; Russo-Averchi, E.; Morante, J. R.; Cantoni, M.; Marzari, N.; Arbiol, J.; Zunger, A.; Warburton, R. J.; Fontcuberta i Morral, A. Nat. Mater. 2013, 12, 439–444. doi:10.1038/nmat3557
    Return to citation in text: [1] [2] [3] [4] [5] [6]
  3. Rudolph, D.; Funk, S.; Döblinger, M.; Morkötter, S.; Hertenberger, S.; Schweickert, L.; Becker, J.; Matich, S.; Bichler, M.; Spirkoska, D.; Zardo, I.; Finley, J. J.; Abstreiter, G.; Koblmüller, G. Nano Lett. 2013, 13, 1522–1527. doi:10.1021/nl3046816
    Return to citation in text: [1] [2] [3] [4] [5]
  4. Zheng, C.; Wong-Leung, J.; Gao, Q.; Tan, H. H.; Jagadish, C.; Etheridge, J. Nano Lett. 2013, 13, 3742–3748. doi:10.1021/nl401680k
    Return to citation in text: [1] [2]
  5. Jiang, N.; Gao, Q.; Parkinson, P.; Wong-Leung, J.; Mokkapati, S.; Breuer, S.; Tan, H. H.; Zheng, C. L.; Etheridge, J.; Jagadish, C. Nano Lett. 2013, 13, 5135–5140. doi:10.1021/nl4023385
    Return to citation in text: [1]
  6. Kasu, M.; Kobayashi, N. J. Cryst. Growth 1997, 170, 246–250. doi:10.1016/S0022-0248(96)00522-2
    Return to citation in text: [1] [2] [3]
  7. Zhang, Q.; Aqua, J.-N.; Voorhees, P. W.; Davis, S. H. J. Mech. Phys. Solids 2016, 91, 73–93. doi:10.1016/j.jmps.2016.02.033
    Return to citation in text: [1] [2] [3] [4] [5] [6]
  8. Zhang, Q.; Voorhees, P. W.; Davis, S. H. J. Mech. Phys. Solids 2017, 100, 21–44. doi:10.1016/j.jmps.2016.12.005
    Return to citation in text: [1] [2] [3] [4]
  9. Carter, W. C.; Roosen, A. R.; Cahn, J. W.; Taylor, J. E. Acta Metall. Mater. 1995, 43, 4309–4323. doi:10.1016/0956-7151(95)00134-H
    Return to citation in text: [1] [2]
  10. Cahn, J. W.; Carter, W. C. Metall. Mater. Trans. A 1996, 27, 1431–1440. doi:10.1007/BF02649804
    Return to citation in text: [1]
  11. Garcke, H. Jahresber. Dtsch. Math.-Ver. 2013, 115, 63–100. doi:10.1365/s13291-013-0066-2
    Return to citation in text: [1]
  12. Carter, W. C.; Taylor, J. E.; Cahn, J. W. JOM 1997, 49, 30–36. doi:10.1007/s11837-997-0027-2
    Return to citation in text: [1]
  13. Biasiol, G.; Gustafsson, A.; Leifer, K.; Kapon, E. Phys. Rev. B 2002, 65, 205306. doi:10.1103/PhysRevB.65.205306
    Return to citation in text: [1]
  14. Tersoff, J. Appl. Phys. Lett. 2003, 83, 353–355. doi:10.1063/1.1592304
    Return to citation in text: [1] [2]
  15. Spencer, B. J.; Voorhees, P. W.; Tersoff, J. Phys. Rev. B 2001, 64, 235318. doi:10.1103/PhysRevB.64.235318
    Return to citation in text: [1]
  16. Tu, Y.; Tersoff, J. Phys. Rev. Lett. 2004, 93, 216101. doi:10.1103/PhysRevLett.93.216101
    Return to citation in text: [1]
  17. Tu, Y.; Tersoff, J. Phys. Rev. Lett. 2007, 98, 096103. doi:10.1103/PhysRevLett.98.096103
    Return to citation in text: [1]
Other Beilstein-Institut Open Science Activities