1. Introduction
Semiconductor NWs, particularly III-V NWs, are widely considered fundamental building blocks for nanoscience and nanotechnology and useful for applications in nanoelectronics and nanophotonics [
1,
2,
3,
4,
5]. Very efficient elastic stress relaxation on strain-free NW side facets allows for dislocation-free growth in material systems with high lattice mismatch [
6,
7,
8]. For example, planar growth of InAs on Si substrates (lattice mismatch = 11.6%) is difficult and leads to Volmer–Weber islands [
9], which often contain crystallographic defects, while fully coherent growth of InAs NWs on Si is possible, provided that the NW diameter is smaller than critical (~25 nm [
4]). A high aspect ratio (length over radius) is crucial for the applications and fundamental physical properties of NWs, such as one-dimensional transport of charge carriers, directional light emission, crystal purity within axial or radial heterostructures, abruptness of heterointerfaces, and crystal phase switching in III-V NWs [
10]. These features explain the importance of controlling the NW dimensions during growth.
III-V NWs are often fabricated by MBE via metal-catalyzed vapor–liquid–solid (VLS) growth (with either Au [
11,
12,
13,
14,
15,
16] or group III [
17,
18,
19] droplets) or catalyst-free SAG [
20]. At low temperatures, NW growth transitions to the vapor–solid–solid (VSS) mode controlled by solid-phase diffusion through a frozen nanoparticle [
21]. Catalyst-free SAG of NWs necessarily requires patterning of the substrate surface. SAG-MBE growth of VLS III-V NWs is achieved by preparation of Au seeds inside the pinholes in SiO
x [
16] or SiN
x [
22] mask layers or pre-deposition of a group III metal (usually, Ga or In) into the pinholes in SiO
x/Si(111) for the self-catalyzed VLS process [
18,
19]. SAG of NWs offers several important advantages over growth on unpatterned substrates [
12,
13,
14,
15], including the absence of parasitic layer between the NWs [
16] and, consequently, very low surface roughness, improved thermal conductivity in the substrate plane, regular positioning, and narrow size distributions of NWs in terms of both lengths and radii [
10].
Due to the known high volatility of group V species such as As and P, their desorption from a catalyst nanoparticle and different NW facets occurs even at low temperatures. Surface diffusion of group V atoms is negligible [
10]. Conversely, group III atoms may diffuse on different surfaces or incorporate with the NW sidewalls without desorption. The total balance of group III atoms is therefore more relevant for modeling NW growth [
18]. There are two main mechanisms for material transport of group III atoms from the substrate surface to NWs. The first is surface diffusion of adatoms, which stick to the substrate and then diffuse to the NW base and along the sidewalls [
23], as usually considered in modeling NW growth on unpatterned “sticky” substrates [
12,
13,
14,
16]. The second is re-emission of group III atoms from the reflecting substrate covered with an inert mask [
17,
18,
22,
24]. Several authors [
16,
22,
24] attempted to model the growth kinetics of NW arrays grown by MBE on reflecting substrates using a cosine law (Lambert scattering) for re-emitted flux without rigorous justification. Furthermore, shadowing of both direct and re-emitted [
16] or direct [
16,
22] fluxes was not taken into account. An important step forward was taken in Ref. [
18], where Ga-catalyzed VLS growth of GaP NWs in patterned arrays of pinholes in SiO
2/Si(111) was monitored using GaAsP markers. A model was developed that allowed for calculation of the total Ga flux into the NW (contributing to the axial and radial NW growths along with the droplet swelling), without any free parameters. One important conclusion was that the specular re-emission model gave the best fit to the data compared to the cosine law or random angular reflection. However, no analytical expression for the total flux was given, and no NW growth modeling was presented.
In this work, I try to fill the gap by developing an analytic model for the total group III flux influenced by specular re-emission and shadowing. This allows one to determine the shadowing length of NWs corresponding to the full shadowing of the substrate surface in MBE. Different stages of NW growth are considered, including purely axial growth, axial and radial growths after reaching the critical diffusion length, and the asymptotic growth stage, where each NW receives a maximum flux determined by the pitch of a regular array. Good fits with the data on the growth kinetics of Ga-catalyzed GaP NWs, Au-catalyzed InP NWs, and catalyst-free InAs NWs are obtained. The model is quite general and should work equally well for metal-catalyzed and catalyst-free NWs in different material systems, where MBE SAG is performed on a masked substrate.
2. Assumptions and Model Parameters
A full description of MBE growth of III-V NWs on reflecting substrates that accounts for all possible factors, influencing the growth process and NW morphology, is beyond reach to this end. Here, a simplified analytical model is developed using the following assumptions. First, an NW is approximated as a cylindrical rod with uniform radius
from base to top, neglecting all tapering effects. An approximately cylindrical shape is often observed in self-catalyzed VLS III-V NWs [
17,
18], some Au-catalyzed III-V NWs [
13], and catalyst-free SAE III-V NWs [
20]. However, more complex tapered or pencil-like shapes are also possible [
14]. Second, MBE growth on rotating substrates is considered, where the particular geometry of an array of pinholes (for example, square or hexagonal) should not critically influence the NW growth process. Therefore, only the NW surface density
enters the results. It can be related to the surface area per NW
as
, where
is the pitch of the square array. The results can easily be re-formulated for any geometry (for example, using
, where
is the pitch of the hexagonal array [
18]). Some works, for example, Refs. [
18,
24], explicitly take the array geometry into account, in which case deriving any analytical expression for the reflected flux is very difficult. Consequently, here we only study the dependence of the NW growth kinetics on
or
, assuming the influence of the array geometry as being a second-order effect. Third, an idealized ensemble of NWs having identical lengths
and radii
is considered. Fourth, the incorporation limited diffusion length of group III adatoms on the NW sidewalls,
, is introduced [
10]. With the neglect of evaporation of group III adatoms from the NW side facets, surface incorporation leads to radial growth.
In the absence of desorption of group III atoms from the NW sidewalls and top (with or without a catalyst nanoparticle), the total balance of group III atoms is given by
. Here,
is the direct group III flux onto the surface (nm/s),
is the group III flux reflected from the substrate,
is the NW volume including the nanoparticle resting on its top for VLS or VSS growth process, and
, with
as the direct group III atomic current onto the NW (nm
3/s). When the substrate surface is entirely blocked by NWs due to the shadowing effect in the directional MBE technique, the reflected flux becomes zero. Starting from this moment of time, the total balance of group III atoms gives
(at
), as in Ref. [
18], meaning that each NW receives the maximum flux
, which equals the group III flux times the surface area per one NW
. Thus, the variables considered are the time-dependent NW length
and radius
, with the control parameters
,
,
, and the initial NW radius
.
3. Total Material Flux into the Nanowire
MBE growth of NWs on a reflecting substrate is illustrated in
Figure 1. The group III flux onto the substrate surface equals
, with
as the total group III flux and
as the beam angle with respect to the substrate normal. The nominal thickness of planar material deposited by the moment of time
equals
. I take into account re-emission (reflection) of group III atoms from the mask surface and shadowing of the mask by the NW array. Re-emission from the NW sidewalls or catalyst nanoparticles is neglected. Hence, any material exchange between the NWs [
25,
26] is not considered. Surface growth on the mask is assumed negligible. As mentioned above, all NWs in the array are considered identical in terms of their length and radius and modeled as cylinders with length
L above the mask surface and uniform radius
R.
The total current into the NW,
, contains contributions from the direct (
) and reflected (
) fluxes. The maximum current is given by
, as discussed above. The direct current is given by
if
where
is the NW surface area (including the nanoparticle surface area for metal-catalyzed NWs or NW top facet for catalyst-free SAE NWs) exposed to the direct flux. The expression
is the common definition of the direct material flux onto the NW [
10,
18,
27,
28]. It shows that the volume of group III atoms impinging the NW per unit time equals the group III flux times the total surface area of this NW exposed to the flux. However, this expression does not account for the shadowing effect. When
becomes larger than
, the direct flux can no longer increase. Instead, it stabilizes at the maximum value
, meaning that each NW in the array of
NWs per unit surface area receives a part of the total flux
.
Similarly, the reflected current impinging the NW is
, where
is the NW surface area (including its top) exposed to the reflected flux. From the total balance of group III atoms, it follows that
and
if
. These considerations give the total flux in the form
For the total flux normalized to its maximum value
, Equation (1) takes the form
For VLS NWs, the total group III flux equals the derivative of the NW volume plus the droplet volume with respect to time
Here,
is the known geometrical function of the droplet contact angle
[
6]. For catalyst-free SAG NWs,
, and only the first term remains on the right side of Equation (3). This model for the total flux is purely geometrical, as in some previous works on NW growth modeling [
14,
27], and captures the two major effects: (i) re-emission from the mask and (ii) shadowing of the mask surface, which finally leads to a given maximum flux per NW determined by the array pitch [
28].
The above expressions are insensitive to the re-emission law. In the case of specular reflection, where group III atoms just bounce off the mask surface [
18], the surface areas exposed to the direct and reflected fluxes are given by
Here, the first terms stand for the surface area of NW sidewalls exposed to the flux and are identical for the direct and reflected fluxes in the case of specular re-emission. The droplet surface areas intercepted by the fluxes contain different geometrical coefficients
and
, which can be obtained using the approach of Ref. [
29] as functions of the two angles
and
. In the case of catalyst-free SAG NWs, one simply has
, meaning that the flux impinging the flat top surface is the same as for the substrate surface and no reflected atoms can impinge the NW top.
Introducing
, Equation (2) can be put in the dimensionless form
with
given by Equation (4). This dependence is shown in
Figure 2 for catalyst-free SAG at a constant NW radius
. The dashed curve shows the approximation
, which becomes more accurate for NWs with higher aspect ratios. It clearly shows that the flux impinging the NW side facets equals twice the direct flux at the beginning of growth due to re-emission and becomes
at
due to shadowing.
Full shadowing of the mask surface occurs at
. From Equation (4), one obtains the shadowing length in the form
Here,
is the NW radius and
is the geometrical function of the droplet angle
reached at the moment of time where
. This
depends on the NW radius because thicker NWs shadow a larger surface area, and hence the full shadowing of the mask occurs at a shorter NW length. The shadowing length is shown in
Figure 3 as a function of the pitch (a) and radius (b) for different beam angles, at a fixed
of 135°, corresponding to
2.0. As expected, the shadowing length increases with the pitch and decreases for thicker NWs and larger beam angles. The radial NW growth will never occur
, with
being the diffusion length of group III adatoms on the NW sidewalls, limited by surface incorporation in the absence of desorption, as discussed above. In this case, group III atoms are collected from the top part of the NW of height
, while the rest of the NW is shadowed. All these atoms will diffuse to the NW top due to
, and none of them will incorporate with the NW sidewalls. The published estimates for the Ga diffusion length on (110) side facets of <111>-oriented GaAs NWs equal 1500 nm [
14], 1800 nm [
30], and more than 2000 nm [
13]. The shaded zones in
Figure 3 correspond to the range of pitches and NW radii, where no radial growth occurs at
1500 nm.
It is well known that fabrication of thin III-V NWs by MBE is difficult, partly due to the radial growth [
10]. According to our results, the radial growth can be more easily suppressed in dense arrays of NWs corresponding to smaller pitches and in MBE systems with larger beam angles
for group III atoms. It is desirable to ensure the condition
, where the NW radius
should be uniform from base to top and defined by the size of the growth seeds.
4. Evolution of Nanowire Length at a Constant Radius
The governing equation for the NW growth kinetics is given by
with
and
determined by Equation (4). Below, I will use
Meaning that the droplet volume remains constant during growth. This result is exact for catalyst-free SAG NWs, where the droplet volume is zero. It may also be justified for Au-catalyzed VLS NWs [
10,
12,
13,
15], in which case there is always a steady-state solution for a time-independent droplet volume, regulated by chemical potential of group V atoms in an Au-III-V droplet. However, Equation (8) is only approximate for self-catalyzed VLS NWs, where a group III droplet may swell or shrink depending on the atomic V/III ratio [
18,
31,
32]. In the case of droplet swelling under excessive Ga flux, subsequent radial growth by step flow starting from the NW top leads to enlargement of the NW radius [
18]. In any case, the ratio of the droplet volume over the NW volume scales as
, and hence the approximation given by Equation (8) is justified for long enough NWs having high aspect ratios.
Using Equations (4), (7), and (8) at
, one arrives at
The exact solution for the NW length is given by
with coefficients
The NW elongation after the full shadowing of the substrate surface is given by
where
is the deposition thickness at which
.
Figure 4 shows the evolution of the NW length with deposition thickness at a fixed NW radius of 50 nm, a beam angle of 32.5°, and a droplet contact angle of 125° for three different pitches of 200 nm, 300 nm, and 400 nm. The corresponding growth parameters are summarized in
Table 1. As expected, the NWs elongate faster for larger pitches. It is noteworthy that
in most cases, corresponding to a pitch-independent characteristic thickness
, which does not change significantly with the pitch according to
Table 1. Overall, the NWs elongate faster at the beginning of growth, while the length evolution becomes linear in
. The initial growth stage is, however, more complex than the exponential increase in the NW length with time. The exponential growth stage was theoretically predicted and experimentally observed for NW growth on sticky substrates [
13,
14,
33,
34,
35] and actually regardless of the epitaxy technique (directional MBE method or vapor-phase epitaxy). This important difference stems from the fact that the number of group III atoms collected by the NW sidewalls in MBE on sticky substrates is proportional to
, while in our case, it scales as
at the beginning and gradually decreases to
in the course of growth. In both cases, the group III flux converges to
due to shadowing, but NWs on a reflecting substrate grow faster before that.
5. Radial Growth in the Intermediate Growth Stage
The case of
corresponds to the most complex scenario of the NW growth kinetics and morphology evolution, even in the case of cylindrical NW geometry without tapering. The radial growth starts before the full shadowing of the mask (when
reaches
), where the total influx of group III atoms continues to increase. The moment of time where
=
, with
. Further evolution of the NW volume is given by
assuming again a time-independent droplet volume. This equation is insufficient for finding
and
separately, which is why a second equation for
is needed. Using similar considerations as presented in
Section 2 for the upper part of an NW of height
, the evolution of the NW length is given by
Here, the first term gives the direct flux impinging the NW top section [
13,
14], while the second term describes the contribution from re-emitted group III atoms. From Equations (13) and (14), the NW radius evolves according to
These two equations can be resolved only numerically.
However, with neglect of the contribution from the flux impinging the droplet or the top NW facet, Equation (14) simplifies to
From Equations (15) and (16), the evolution of NW radius with its length is given by
This equation has the solution
Therefore, the NW volume scales exponentially with its length
This relationship contains only one fitting parameter and is independent of the re-emission and shadowing. The latter are described by the factors in Equations (15) and (16) but cancel in Equation (17) for . Exponential dependence is expected to be more accurate for higher aspect ratio NWs and at larger diffusion length , corresponding to an almost negligible contribution of the NW top into the total collection area for group III atoms impinging the NW.
6. Asymptotic Growth Stage
This late stage of NW growth occurs for
, where the substrate surface is entirely shadowed by the NW array and hence is exactly identical for MBE growth of NWs on any substrate. The condition of a time-independent droplet volume becomes almost exact because the aspect ratio
increases in the course of growth. The governing equations are reduced to
with the initial conditions
and
. Integrating Equation (20), the NW radius is obtained in the form
Substitution of Equation (22) into Equation (21) yields
This is a special type of the Chini equation considered recently in Ref. [
36], where it was used from the beginning of growth, assuming that all diffusive group III adatoms are collected by NWs from a “sticky” substrate surface. In our case, the solution is given by
and is controlled by the two parameters
and
. The latter is inversely proportional to the array pitch, and hence the NW length should finally decrease with the pitch, the trend, which is inversed with respect to the growth start (see
Figure 4).
In the large time limit, the
ratio and the radius saturate to the constants
Therefore, the NW length scales linearly with the deposition thickness and increases for larger and . The maximum possible NW radius corresponds to coalescence of the NW array into continuous film. According to Equation (25), the steady-state radius decreases for larger and . The NWs will coalesce only at (the absence of surface diffusion from the NW sidewalls to the top) and , as in catalyst-free SAG NWs. Whenever , surface diffusion of group III adatoms does not allow for coalescence of NWs into continuous film, that is, the NW length always remains larger than and the radius smaller than .
Figure 5 show the evolution of the NW lengths and radii, obtained from Equations (22) and (24) at a fixed
°,
° (
)
75 nm,
500 nm, and three different pitches
300 nm, 400 nm, and 500 nm. The corresponding growth parameters are summarized in
Table 2. For the smallest pitch of 300 nm, the radius stays almost constant, meaning that its steady-state limit has been reached in the earlier growth steps. At a constant radius, the NW length scales linearly with
. Increasing the pitch leads to a more pronounced radial growth such that the steady-state NW radius enlarges to 110 nm for
400 nm and 145 nm for
500 nm. As a consequence, the NW lengths increase more sublinearly with the deposition thickness. In the final step, the length of NWs grown in the 400 nm pitch array becomes smaller than that in the 300 nm pitch array. Such an inverse pitch dependence of the NW length was observed experimentally in Ref. [
37].
7. Theory and Experiment
Self-catalyzed MBE growth of GaP NWs of Ref. [
18] was performed at 600 °C using the SAG approach in a regular hexagonal array of patterned holes in a SiO
2 mask layer on Si(111) (
216,506 nm
2), with
32.5° and
0.135 nm/s. Therefore, 60 min growth corresponded to
486 nm. In total, 122 GaAsP markers were introduced at fixed time intervals to monitor ex situ the axial growth and radial extension of individual NWs. This linear dependence is used to convert the marker number to the Ga deposition thickness. The droplet contact angle stayed nearly constant at
135°, corresponding to the region of zincblende NWs, where the droplet volume increases by enlarging the base radius [
38,
39,
40]. An almost untapered NW geometry was kept by step flow radial growth starting from the NW top [
18]. In 60 min growth, the NW radius enlarged from 20 to 75 nm, approximately following the parabolic dependence
nm. The NW axial growth rate increased linearly with
from 0.3 nm/s to 1 nm/s at
105.3 nm and stayed constant at 1 nm/s until the end of growth. Therefore, I use
nm for
nm and
for
105.3 nm.
With these
and
, Equations (2) and (4) provide the fit to the measured Ga flux shown in
Figure 6, which works equally well compared to a dedicated model of Ref. [
18]. The discrepancy can be seen only for the highest Ga fluxes measured, where the model of Ref. [
18] is not perfect either. Overall, our model predicts the maximum flux being reached at
, similarly to the model of Ref. [
18], while the measured flux seems to decrease at the end of growth. The reason for this mismatch should be considered more carefully. One possible explanation is the partial evaporation of Ga atoms from the NW tip or the droplet in the late stage of growth, which is forbidden in our model.
Catalyst-free InAs NWs of Ref. [
20] were grown by SAG-MBE in patterned arrays of pinholes in SiO
2/Si(111) at 480 °C. Square arrays with pitches
ranging from 250 to 3000 nm were investigated by measuring the average values of the NW lengths and diameters grown under identical conditions for different times. The data for the average values of NW lengths and diameters are given in
Table 3, with the corresponding error bars shown in
Figure 7a. Both lengths and diameters increased sublinearly with time [
20]. However, if one plots the NW volume versus length, the exponential dependence fits all the datapoints within the error bars. The best fit is obtained using Equation (19) with
560 nm and
50 nm, as shown in
Figure 7a. Including the shadowing effect, which is important for the smallest pitches, leads to the curves shown by the dashed lines in
Figure 7a.
In Ref. [
16], SAG of Au-catalyzed InP NWs was performed by chemical beam epitaxy (another directional deposition method that is similar to MBE in this respect) in hexagonal arrays of patterned holes in SiO
2 on InP(111)B substrates at 420 °C. Arrays with different pitches were investigated ranging from 100 to 700 nm. The growth time was 15 min, corresponding to a nominal planar growth of 58 nm thick InP for the data shown in
Figure 7b. The In beam angle
was 45°. InAs markers were used to measure the growth kinetics of a thin InP NW having approximately uniform radius
12 nm from base to top. Although the authors claimed no significant pitch dependence of the NW length [
16], the exponential curve based on Lambert re-emission without the shadowing effect shows a discrepancy for the longest growth time (dashed line in
Figure 7b. The bold line in
Figure 7b shows the fit obtained from Equations (10) to (13) for a hexagonal array with a pitch of 200 nm and gives a better quantitative agreement with the data compared to the original fit without shadowing.