Open Access Article
T.
Dursap
a,
M.
Vettori
a,
A.
Danescu
a,
C.
Botella
a,
P.
Regreny
a,
G.
Patriarche
b,
M.
Gendry
a and
J.
Penuelas
*a
aInstitut des Nanotechnologies de Lyon-INL, UMR 5270 CNRS, Université de Lyon, École Centrale de Lyon, 36 avenue Guy de Collongue, F-69134 Ecully cedex, France. E-mail: jose.penuelas@ec-lyon.fr
bUniversité Paris-Saclay, CNRS, Centre de Nanosciences et de Nanotechnologies, 91120, Palaiseau, France
First published on 13th April 2020
It is well known that the crystalline structure of the III–V nanowires (NWs) is mainly controlled by the wetting contact angle of the catalyst droplet which can be tuned by the III and V flux. In this work we present a method to control the wurtzite (WZ) or zinc-blende (ZB) structure in self-catalyzed GaAs NWs grown by molecular beam epitaxy, using in situ reflection high energy electron diffraction (RHEED) diagram analysis. Since the diffraction patterns of the ZB and WZ structures differ according to the azimuth [1
0], it is possible to follow the evolution of the intensity of specific ZB and WZ diffraction spots during NW growth as a function of the growth parameters such as the Ga flux. By analyzing the evolution of the WZ and ZB spot intensities during NW growth with specific changes of the Ga flux, it is then possible to control the crystal structure of the NWs. ZB GaAs NWs with a controlled WZ segment have thus been realized. Using a semi-empirical model for the NW growth and our in situ RHEED measurements, the critical wetting angle of the Ga catalyst droplet for the structural transition is deduced.
The WZ phase is often obtained for gold-catalyzed III–V NWs, while the ZB phase is mostly obtained for self-catalyzed ones. In the latter case, the WZ phase is, however, often observed at the NW top near the Ga droplet and can be ascribed to the end of NW growth when the Ga flux is stopped and the Ga droplet is consumed under the As flux.18–24 The control of the ZB and WZ phases in self-catalyzed GaAs NWs has thus become a major challenge. Based on the models of F. Glas7 and V. Dubrovskii,25 occurrence of the WZ or ZB phase in self-catalyzed GaAs NWs has been explained by the position of the nucleation for a new atomic layer either at the triple phase line (TPL) or inside the droplet, respectively.19,22,26–28 It was recently shown that the nucleation position mainly depends on the droplet wetting angle and therefore on the catalyst droplet volume, for both Au-catalyzed and self-catalyzed GaAs NWs.29–32 A critical wetting angle βc1 in the 121°–124° range for Au catalyzed NWs29,30 and of about 125°–127° for self-catalyzed NWs31,32 has been experimentally observed for a transition from the WZ to the ZB crystal phase above this critical angle. In addition, a second critical angle βc2 in the 85°–100° range below which a transition from the WZ to the ZB crystal phase is also expected to occur has been experimentally observed.32 From the in situ transmission electron microscopy (TEM) observations by Jacobsson et al., above this angle a truncated facet is present at the NW top facet, thus determining that the nucleation site is inside the droplet leading to the ZB crystal phase.29 In self-catalyzed growth, the Ga droplet volume is mainly dependent on the growth conditions, in particular on the Ga and As fluxes. Many previous studies have thus reported the control of the GaAs NW crystal phase by tuning the Ga and/or As fluxes.22,27,31,33–37 In these studies the crystal phase characterization has been mainly performed ex situ by TEM. However, an in situ and real-time characterization tool also appeared to be very useful in order to characterize and possibly tune the crystal phase of the self-catalyzed GaAs NWs during growth. Compared to in situ X-ray diffraction35 and in situ TEM,29,30 the reflection high energy electron diffraction (RHEED) technique is commonly coupled with molecular beam epitaxy (MBE) to follow the structural properties of growing layers and nanostructures. Despite this, there are relatively few studies reporting RHEED observations during the self-catalyzed GaAs NW growth. Scarpellini et al.24 reported on the consumption of the Ga droplets at the end of the growth, while Rudolph et al.21 reported on the influence of the As flux on the InAs NW structural properties and Bastiman et al. reported on the incubation time of GaAs NWs.38 Only recently, in situ RHEED characterization of the NW growth coupled with ex situ TEM measurements was reported by Jo et al.39
In this work, we focus on the characterization of the growth of self-catalyzed GaAs NWs on a Si(111) substrate using in situ RHEED. In particular, we aim to control the formation of the ZB or WZ crystal phase of the GaAs NWs as a function of the Ga flux amounts by using the RHEED pattern. TEM measurements of some NWs were performed to check the obtained crystal phases.
0] azimuth of self-catalyzed GaAs NWs when both the WZ and ZB phases are present is shown in Fig. 1. The position of the spots is in agreement with epitaxial growth of NWs on Si(111): the GaAs [111] and [1
0]axes are parallel to the Si [111] and [1
0] axes, respectively. The corresponding spots of the WZ phase and twinned ZB phase are indexed. The (10–12) WZ and (002)t ZB spots whose intensities will be measured during growth are indicated by the green and red arrows, respectively. The SEM picture in Fig. 2(a) shows the NW morphology after 10 min of growth: their length is about 1.1 μm and their diameter is in the 40–60 nm range, with a density close to 1 NW μm−2. Some parasitic nanocrystals among the NWs are evidenced; however the substrate surface is still visible since the growth temperature was optimized in order to minimize this parasitic growth. Fig. 2(b) shows the typical RHEED pattern obtained along the [1
0]azimuth during the NW growth. Only ZB spots are observed indicating the growth of pure ZB NWs with twin planes. Fig. 2(c) shows the RHEED pattern at the end of the NW growth after closing the Ga shutter and cooling the sample under As4 flux. We can observe the presence of low intensity WZ spots (two of them are indicated by red arrows). Fig. 3(a–d) show TEM and HRTEM images (with the [1–10] zone axis) of a typical self-catalyzed GaAs NW. Under the applied growth conditions (TG = 600 °C and V/III flux ratio = 2.3) the NW exhibits a pure-ZB phase almost over its entire length. Both ZB variants can be observed due to the presence of some twin planes (marked in green in Fig. 3(a) and shown in Fig. 3(d)). Then, near the NW top, we observe a sequence with a transition zone about 60 nm in length with a high density of twin planes and some WZ segments (marked in blue in Fig. 3(a) and shown in Fig. 3(c)), and finally a pure WZ segment about 100 nm in length followed by a thin ZB segment about 10 nm in length (marked in red in Fig. 3(a) and shown in Fig. 3(b)).
This final sequence and the associated growth mechanism are well known thanks to previous studies.18–24,26 Indeed, the GaAs NW growth continues with the consumption of the Ga droplets leading to a decrease of its volume. Recent in situ TEM results obtained on Au- or Ga-catalyzed GaAs NWs29–32 have confirmed that nucleation of the ZB or WZ structure will occur depending on the droplet shape and more precisely on the wetting angle β which determines the location of the atomic layer nucleation. From these recent results, and as illustrated in Fig. 4(a), the final structural sequence of the NWs can be explained as follows depending on β: first, the NW growth is presumed to be with a large wetting angle βmax,31 typically greater than βc1 = 125°–127°, leading to nucleation inside the droplet, thereby giving the growth of a ZB phase according to the F. Glas model7 (Fig. 3(d) and 4(a) stage (i3)). Once the Ga flux is closed, the droplet starts to be consumed leading to a decrease of its volume and so of the wetting angle β, hence leading to nucleation at the TPL thereby giving the growth of a WZ phase, also according to the F. Glas model.7 This step leads first to the same ZB phase, as long as β is greater than βc1 (Fig. 4(a) stage (ii)), and then to a ZB/WZ phase mixing segment followed by a pure-WZ segment, when β becomes smaller than βc1 (Fig. 3(c and b) and Fig. 4(a) stage (iii)). Once β becomes typically lower than βc2 = 85°–100° a final ZB segment is formed (Fig. 3(b) and 4(a) stage (iv)) up to βmin where the droplet is unpinned from the top edge of the NW and moves freely on the (111) top facet maintaining a wetting angle equal to βmin and thus also leading to a ZB phase, according to previous studies,24,26 until total consumption of the droplet (Fig. 4(a) stage (v)).
The first purpose of this work was to follow this structural evolution via the RHEED pattern during NW growth and at the growth end when the Ga flux is closed. To proceed, RHEED videos were recorded along the [1
0] azimuth during the NW growth and the time-dependent intensities of the ZB and WZ spots were analyzed. The analyzed spots giving the ZB and WZ phase “intensities” are indicated in Fig. 1. Fig. 4(b) shows the evolution of the
(or
) intensity ratio (IR), where IZB and IWZ are the intensities of the ZB and WZ spots, respectively, as a function of the growth time. The green area corresponds to the NW growth with Ga and As fluxes (zones (i)). At time t = 0 seconds, WZ and ZB IRs are equal to 0.5, due to the absence of NWs (the measured intensities are those measured on the diffraction line of the Si(111) substrate surface). Once the Ga and As fluxes are opened, an increase of the ZB IR
(correspondingly, a decrease of the WZ IR
) is observed due to the growth of the parasitic ZB GaAs nanocrystals during the first few seconds, and then due to the nucleation and the growth of the NWs after a few tens of seconds38 (zone (i1)). After 90 seconds, we observe a short and slight decrease of the ZB IR (correspondingly, a short and slight increase of the WZ IR) for about 40 seconds (zone (i2)). Subsequently, we again notice an increase of the ZB IR (correspondingly, a decrease of the WZ IR), which tends to become constant with the growth time with a value close to 1 for the ZB IR (correspodingly close to 0 for the WZ IR) (zone (i3) and the corresponding stage (i3) in Fig. 4(a)), meaning that the NWs are entirely or quasi-entirely ZB at this moment (see Fig. 3(a)).
After 600 seconds, when the Ga flux is closed, we can assume that the Ga droplet volume starts to decrease thus leading to the sequence schematized in stages (ii) to (v) of Fig. 4(a). A delay is observable (zone (ii) in Fig. 4(b)) between the closing of the Ga flux and the beginning of the WZ IR increase (correspondingly, the ZB IR decrease). This delay measured at around 20 seconds is interpreted to be due to the consumption time of the droplets which is necessary to reach the wetting angle βc1 = 125°–127° leading to the WZ phase (corresponding to stage (ii) in Fig. 4(a)). Then, we observed an increase of the WZ IR during about 70 seconds (zone (iii) in Fig. 4(b) corresponding to stage (iii) in Fig. 4(a)), associated with the growth of the ZB/WZ phase mixing segment (zone t in Fig. 3(c)) and then with the pure-WZ segment (see Fig. 3(b)). A weak decrease of the WZ IR (zone (iv) in Fig. 4(b) corresponding to stage (iv) in Fig. 4(a)) is then visible and associated with the growth of the final ZB segment (see Fig. 3(b)) when β becomes typically lower than βc2 = 85–100°. It can be noted that this WZ-to-ZB phase transition is abrupt, not giving rise to a zone with structural defects. Finally, a stabilization of both ZB and WZ IRs is observed (zone (v) in Fig. 4(b) corresponding to stage (v) in Fig. 4(a)) and is associated with the end of NW growth when β = βmin, with the unpinning of the droplets from the top edge of the NWs24,26 and the final consumption of the Ga droplets. These RHEED measurements are perfectly in line with the structural evolution of the NWs observed by TEM measurements and we can therefore associate them with the size and wetting angle modifications of the Ga droplets, as illustrated in Fig. 4(a).
It should be noted that the observed RHEED diagram can be affected by the density and the length of the NWs, as well as by the incidence angle of the electron beam due to shadowing effects. All of these parameters can modify the RHEED probed depth. Nevertheless, in our specific configuration we estimated this probed depth to be around 400–500 nm. This estimation was obtained as explained in the ESI (Fig. 1 S.I.†). However, the purpose of this work is to monitor the evolution of the RHEED diagram during the growth, and in particular to relate the time dependence of the RHEED spots to the growing crystal phase. These shadowing effects have hence only a slight influence on our data.
intensity ratio. After the next Ga flux opening, it can be observed that the increase of the WZ IR did not stop immediately. This delay of about 26 seconds (see after) corresponds to Ga droplet refeeding which had not yet reached the wetting angle necessary for the transition to the ZB phase. In order to determine these durations more precisely, the first derivative of the WZ IR is plotted for the growth between 200 and 400 seconds (Fig. 6). The period without the Ga flux is represented by the red area. The ZB-to-WZ and WZ-to-ZB phase transitions are indicated by the vertical black lines, where the first derivative of the WZ IR is equal to zero. These phase transitions are used to estimate the lengths of the bottom ZB segment and of the WZ segment (including the ZB/WZ and WZ/ZB phase mixing segments) (see Table 1). In order to calculate the lengths of the bottom ZB and WZ segments, an average axial growth rate ν = 1.7 nm s−1 is used (from the average length of the bottom ZB + WZ segments measured on around twenty NWs). From Fig. 6, we approximate the growth duration of the WZ segment (including the ZB/WZ mixed segments) to be tWZ = 69 seconds, thus corresponding to a WZ segment length of ∼117 nm. A comparison between the calculated bottom ZB and WZ segment lengths and the measured bottom ZB and WZ segment lengths obtained by TEM measurements on around twenty NWs (a typical WZ segment is shown in Fig. 5(c)) is reported in Table 1.
| Growth time (seconds) | Calculated length (nm) (from RHEED) | Measured length (nm) (from TEM) | Simulated length (nm) | |
|---|---|---|---|---|
| Bottom ZB segment | 257 | ∼437 | 440 ± 70 | 450 |
| WZ segment (including the two ZB/WZ phase mixing segments) | 69 | ∼117 | 120 ± 20 | 106 |
By extending the semi-empirical growth model for Au-catalyzed GaAs25,33,43 and InAs44 NWs, we proposed in an earlier report by Vettori et al.45 a model for self-catalyzed GaAs NWs able to account for the Ga droplet evolution as a function of (a) the Ga and As atoms entering the droplets by direct impingement of the Ga and As fluxes and (b) the Ga atoms entering the droplets by diffusion on the SiO2-terminated Si substrate (for short NWs) and on the NW facets (more details can be found in the ESI section†). Motivated by the droplet stability at a solid angle, an original feature of the model revealed by Vettori et al.45 is the existence of an upper-limit wetting angle for the droplet (defined as βmax in the previous sections). As a consequence, besides the classical axial NW growth the model is able to predict conditions that trigger both the axial and the lateral NW growth depending on droplet evolution. In our case, we expect the numerical simulations to confirm the unique (critical) value of the wetting angle that characterizes the ZB-to-WZ and WZ-to-ZB phase transitions. Fig. 7 shows the time-evolution for: (a) the Ga and As amounts of atoms feeding the droplets, (b) the droplet and NW radii, (c) the wetting angle of the droplets and (d) the NW length. The Ga flux closing between 240 and 300 seconds induces a decrease of the droplet radius from 25 nm to 18 nm (shown in (b)) and of the wetting angle from 134° to 98° (shown in (c)). Moreover, the decrease of the wetting angle (or droplet volume) induces a decrease of the droplet capture surface for As atoms and, as a consequence, a slightly lower axial growth rate (shown in (d)). When the Ga flux is opened at 300 seconds the droplet radius and the wetting angle increase again. At 370 seconds the maximum wetting angle βmax is attained (∼135°) so that the NW radius increases (shown in (b)) in order to accommodate a higher droplet volume.
![]() | ||
| Fig. 7 Time evolution for (a) the number of Ga atoms (continuous blue line) and As atoms (dashed red line) feeding the droplet, (b) the droplet (continuous blue line) and NW (dashed red line) radii, (c) the wetting angle and (d) the NW length. The numerical results were obtained using the following numerical values: the Ga and As sources located at incidence angles equal to 27.9° and 41°, respectively, with respect to the normal to the substrate, the nominal Ga and As fluxes given by FGa = 3.53 atoms/sec.nm2 and FAs = 27.3 atoms per sec per nm2. The best fit was obtained using 38 nm and 1400 nm for the diffusion lengths on the SiO2-terminated Si substrate and on the NW facets, respectively, and the As concentration threshold in the droplet is fixed at 1%. From the experimental data, typical values for the initial conditions for the droplet are r = 18 nm and wetting angle = 90°. | ||
The vertical dotted lines in Fig. 7(c) and (d) correspond to t = 257 seconds and t = 326 seconds when the ZB-to-WZ and WZ-to-ZB phase transitions occur as measured in Fig. 6. The corresponding ZB and WZ segment lengths predicted by the numerical simulation are 450 nm and 106 nm, respectively, in good agreement with the TEM measurements of 440 ± 70 nm and 120 ± 20 nm, respectively. The numerical method predicts values of the wetting angle equal to 126° for the ZB-to-WZ transition and to 124° for the WZ-to-ZB transition, in agreement with the results recently reported by Kim et al.31 and measured by ex situ TEM.32 It should also be noted that this critical wetting angle is close to the values reported for Au-catalyzed GaAs NWs.29,30
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/d0na00273a |
| This journal is © The Royal Society of Chemistry 2020 |