Prerpints.org logo
Preprint
Article

Can Nanowires Coalesce?

Altmetrics

Downloads

123

Views

21

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

21 September 2023

Posted:

26 September 2023

You are already at the latest version

Alerts
Abstract
Coalescence of nanowires and other three-dimensional structures into continuous film is desirable for growing low dislocation density III-nitride and III-V materials on lattice-mismatched substrates, and interesting from the fundamental viewpoint. Here, we develop a growth model for vertical nanowires which, under rather general assumptions on the solid-like coalescence process within the Kolmogorov crystallization theory, results in a morphological diagram for the asymptotic coverage of a substrate surface. The coverage is presented as a function of two variables, the material collection efficiency on the top nanowire facet a, and the normalized surface diffusion flux of adatoms from the NW sidewalls b. The full coalescence of nanowires is possible only when a=1 regardless of b. At a>1, which is usual for vapor-liquid-solid growth with a catalyst droplet, nanowires can only partly merge but never coalesce into continuous film. These results can be used for predicting and controlling the degree of surface coverage by nanowires and three-dimensional islands by tuning the surface density, droplet size, adatoms diffusivity, and geometry of the initial structures in the vapor-liquid-solid, selective area or self-induced growth processes in different epitaxy techniques.
Keywords: 
Subject: Physical Sciences  -   Condensed Matter Physics

1. Introduction

Semiconductor nanowires (NWs) offer almost unlimited possibilities for the bottom-up bandgap engineering and design of Si-integrated III-V optoelectronic device structures [1]. Vertical NWs can be grown by the vapor-liquid-solid (VLS) method [2-4] with a foreign catalyst metal, often Au [2], which can be replaced with a group III metal in the self-catalyzed VLS approach for III-V NWs [3]; selective area epitaxy (SAE) without any droplets [5]; or using the self-induced (SI) nucleation mechanism typical for III-nitride NWs on Si(111) substrates with different interlayers [6,7]. High aspect ratio NWs and other elongated structures are able to relax elastic stress induced by lattice mismatch on their strain-free sidewalls [8-10]. Therefore, they can be grown on dissimilar substrates such as Si without forming misfit dislocations [9]. If vertical NWs grow in both vertical and radial directions, they may merge together or even fully coalesce into continuous film. The coalescence process was studied in much detail for high density III-nitride NWs [11-16], and showed to yield nearly dislocation-free planar layers of GaN and AlGaN on Si [13,15,16]. The NW coalescence was also investigated for vertical InAs [17], ZnO [18], and Ge [19] NWs, where it can be desirable (for example, for obtaining continuous III-V films on Si with a possibly reduced dislocation density compared to planar epi-layers) or undesirable (for applications requiring isolated NWs). The coalescence process in dense ensembles of NWs and other elongated structures which are able to change their aspect ratio in different stages of epitaxial growth is also interesting from the fundamental viewpoint.
Here, we develop the existing models for the vertical NW coalescence [11,20,21] by considering an ensemble of vertical NWs in the coalescence stage and describing the large time asymptote of the surface coverage as a function of the NW geometry (surface density, droplets on the NW tops or flat tops without any droplets) and the diffusion length of adatoms on the NW sidewalls. We also consider the NW growth kinetics and the evolution of the surface coverage in the pre-coalescence stage. Our approach can also be used for studying the coalescence of three-dimensional (3D) islands or planar NWs in SAE [22,23], where the full coalescence of the initially separated structures is always required to obtain 100% coverage of a growth template [24-26]. Such a treatment was not performed so far to our knowledge, and should be useful for understanding and controlling the late growth stages and coalescence process of different 3D structures in a wide range of material systems.

2. Model

We consider an ensemble of vertical NWs (or 3D islands) which grow in both vertical and horizontal directions with respect to the substrate surface. They nucleate with random positions and a surface density N = P − 2 , where P is the average distance between the NW centers. For VLS or SAE growths, the NW density is determined by the pitch separating the pre-existing droplets or patterned pinholes in a dielectric mask, while in the SI growth of GaN NWs it depends on the nucleation mechanism of 3D islands which subsequently become NWs [7,11]. In the SAE process, the randomness in the NW positions is less than on unpatterned substrates. All NWs before coalescence have the same shape such that their cross-sectional area is given by c R 2 , with R  as the NW radius and c  as a shape constant ( c = π  for cylindrical and c = 3 3 / 2 for hexahedral NWs). Different stages of NW growth are illustrated in Figure 1. We are interested in the coalescence process, which should usually occur in a late stage of growth where the NW ensemble consumes the entire flux v of semiconductor material (group III atoms for III-V NWs) sent from vapor [27-29]. In molecular beam epitaxy (MBE), this corresponds to the full shadowing of the substrate surface which either re-emit material or contribute into the total diffusion flux coming to the NW tops [28]. In vapor phase epitaxy (VPE) techniques, this corresponds to the saturation of the total flux per NW which includes the growth species desorbed from different surfaces [29]. Introducing the average NW length L , we can write
d d t θ L = v ,
where θ is the surface coverage. In stage (a) of isolated growth where θ = c R 2 N ≪ 1 , with R as the average NW radius, this is equivalent to d ( c R 2 N ) / d t = v , as in Refs. [28,29]. The growth regime described by Equation (1) starts at a certain NW length L * and coverage θ * (most probably before coalescence, corresponding to θ * = c R * 2 N , with R * as the NW radius at the average length L * [28]). The pre-coalescence stage of NW growth will be discussed in the next section. Integrating Equation (1), we find
θ L = h , h = θ * L * + v ( t − t * ) ,
where t * is the moment of time after which the NW growth kinetics is governed by Equation (1). This equation shows that the volume of the NW ensemble per unit surface area θ L scales linearly with time, and is proportional to the effective deposition thickness h .
For the vertical NW growth rate, we use the standard model of NW growth assuming that tall enough NWs elongate due to direct impingement of material onto their tops and surface diffusion of adatoms (group III adatoms for III-V or III-nitride NWs) from the upper part of NWs to their tops [4,28,30-32], re-written in terms of θ :
θ d L d t = v a θ + Λ u .
Here, a is the collection efficiency on the NW top, with a > 1 for most VLS NWs having droplets on their tops due to a larger surface area of the droplet compared to c R 2   [4,33], and a = 1 for catalyst-free SAE or SI NWs. Λ = Λ i n c denotes the effective collection length of adatoms on the NW sidewalls, limited by surface incorporation to allow for the radial NW growth, and u is the total NW perimeter per unit surface area. At   θ = c R 2 N ≪ 1 and u = 2 c R N , this is reduced to the usual growth law d L / d t = v ( a + 2 Λ / R )   for isolated NWs [4,28,30-32]. NWs can grow infinitely long only when Λ ≥ 0 , which is the interesting case considered in what follows. At Λ < 0 , corresponding to “negative” diffusion of material from the NW tops to their sidewalls, NWs can only grow to a finite length [34] and should then fully coalesce due to the radial growth. The saturation of length at large t for Λ < 0 is seen directly from Equation (3). The growth law given by Equation (3) is equivalent to
d L d h = a + Λ u θ   ,
where u / θ is the average perimeter per unit cross-sectional area of NWs. This elongation law is valid at any θ , starting from θ ≪ 1 to unity. In particular, the NW perimeter u tends to zero at θ → 1 . Therefore, we have L = a h at h → ∞ , which is consistent with Equation (2) at   θ → 1 only when a = 1 . This gives a very simple answer to the main question in the title of this work. Vertical NWs can fully coalesce into continuous film only without any magnifying effect on the material collection at the NW top. In particular, VLS NWs can partly merge, but will never coalesce into continuous film due to a magnifying effect of a catalyst droplet. This conclusion was drawn earlier in Ref. [28]; however, it was based on modeling the growth kinetics of isolated NWs. Any catalyst-free NWs with a = 1 will finally coalesce into continuous film regardless of the value of Λ , simply because surface diffusion becomes ineffective at u → 0 . For example, catalyst-free III-N NWs can fully coalesce [16].
To describe the degree of incomplete coalescence of NWs using Equations (2) and (4), we require a model for the NW perimeter per unit surface area u as a function of θ . According to the Kolmogorov crystallization model, the coverage θ is given by θ = 1 − e x p ⁡ ( − θ 0 ) , where θ 0 is the so-called extended volume, or the coverage in the absence of coalescence [35]. The Kolmogorov model accounts for Poissonian nucleation with random position on a substrate surface, solid-like coalescence process which usually occurs for solid-state structures, and with neglect of the boundary effect. Applying this to the two-dimensional (2D) projection of the NW ensemble in the substrate plane, with θ 0 = c R 2 N , we have
θ = 1 − e − c R 2 N ,     u = d θ d R = 2 c R N e − c R 2 N .
These expressions apply when the nucleation step is short compared to the duration of the whole NW growth process, which is guaranteed for most VLS [4] and SAE [5,36] NWs and is a good approximation for SI GaN NWs [7,11]. It is interesting to note that due to Equation (2) for the total NW volume, we do not need to introduce any model for the NW radial growth rate, whereas the function V t = d R / d t is essential in the Kolmogorov model [35]. Using Equations (5), we obtain
u θ = 2 c N 1 − θ − ln ⁡ 1 − θ .
At small θ = c R 2 N ≪ 1 , this is reduced to u = 2 c N θ = 2 c N R and u / θ = 2 c N / θ = 2 / R . Figure 2 shows the normalized perimeters u and u / θ (in the units of 2 c N ) given by Equation (6), compared to the approximations u / ( 2 c N ) = θ and [ u / θ ] / ( 2 c N ) = 1 / θ for isolated NWs at θ ≪ 1 . The curves for isolated NWs become inaccurate after θ ≅ 0.1 . The perimeter of NWs per unit area of the substrate surface u reaches its maximum at θ m = 1 − 1 / e ≅ 0.39 , corresponding to the beginning of the coalescence process. The perimeter of NWs per unit area of their top facets   u / θ decreases monotonically with θ . This decrease is sharper than the simple 1 / θ   dependence for isolated NWs.
Using Equation (6) in Equation (4), we obtain the elongation law which accounts for the NW coalescence process:
d L d h = a + b L h − 1 − l n 1 − L h , b = 2 Λ c N ,
where h / L = θ according to Equation (2). Hence, the right hand side can equivalently be presented as a function of θ . Numerical soluion of this equation gives the mean NW length L , and simultaneously the surface coverage θ , as a function of h or the growth time. The generalized growth model given by Equation (7), together with a model for the surface coverage θ , is the main result of this work. In the next section, we analyze the asymptotic behavior of the coverage at h → ∞ and L / h → c o n s t depending on the two control parameters a ≥ 1 and b > 0 , and consider the NW growth kinetics before the asymptotic stage.

3. Results and discussion

First, we note that at θ = h / L = c R 2 N ≪ 1 , corresponing to low coverage and high NW aspect ratio, Equations (7) and (2) are reduced to the result of Ref. [28]:
d L d h ≅ a + b L h ,   R ≅ h c N L
This special type of the Chini equation for L ( h ) can be resolved in the analytic form [28]. At h → ∞ , the NW length scales linearly with h : L → ( b 2 + 4 a + b ) 2 / 4 . Without coalescence, the coverage is given by θ = c R 2 N , and equals h / L according to Equation (8) for R . Therefore, the asymptotic coverage at h → ∞ is given by [28]
θ ∞ ≅ 4 b + b 2 + 4 a 2 .
This result, obtained within a simplified model without coalescence, gives correctly the asymptotic coverage θ ∞ = 1   only at a = 1 and b = 0 .
The asymptotic coverage in the general case with coalescence is obtained from Equation (7) at L / h → c o n s t . Using L = c h in Equation (7) and finding the unknown c , the result is given by
a θ ∞ = 1 − b 1 − θ ∞ − l n 1 − θ ∞ .
This transcendent equation for θ ∞ cannot be analytically resoled, but enables one to plot the asymptotic coverage as a function of a and b . Figure 3 shows θ ∞ versus a for different b . As discussed above, the full coalescence of NWs occurs only at a = 1 and regardless of b . For any a > 1 , NWs can only partly merge but will never fully coalesce. The asymptotic coverage monotonically decreases with a for any b , showing that a more efficient material collection at the NW top (enhanced by the presence of a catalyst droplet for VLS NWs) reduces the asymptotic surface coverage. For a given a > 1 , θ ∞ decreases with b . According to Equation (7) for b , the coverage is decreased for higher diffusivities of adatoms on the NW sidewalls having larger diffusion lengths Λ , and for larger surface densities N , or smaller pitches P in regular NW arrays. Decrease of the surface coverage and the corresponding increase of the NW aspect ratio for denser NW ensembles may look counter-intuitive at the first glance. However, this effect is well known, for example, for InAs NWs grown by MBE on patterned SiOx/Si(111) substrates [36], and is explained by the enhanced shadowing of denser NWs and lower material flux per NW contributing to the radial growth [28,29,37]. It does not contradict the fact that merging of NWs occurs earlier in denser NW ensembles, simply because the late stage of growth corresponding to the full shadowing of the substrate surface starts earlier than in sparse ensembles (corresponding to shorter times t * in Equation (1)). This effect will be considered in detail in what follows.
The approximation for isolated NWs given by Equation (9) predicts correctly some qualitative trends, in particular, the possibility of the full coalescence at a = 1 (but only at b = 0 , which is wrong), and a decrease of the asymptotic coverage with b . However, it underestimates the coverage for any a at any positive b , that is, in the growth regimes with surface diffusion of adatoms. These inaccuracies, seen in Figure 3, are related to the overestimated perimeter per NW cross-sectional surface, as shown in Figure 2. In the growth model for isolated NWs, their total perimeter never decreases and never shrinks to zero. This gives a non-vanishing contribution of surface diffusion into the NW elongation rate, resulting in a higher NW aspect ratio and lower coverage of a substrate surface. In the refined model with the coalescence process included, the total NW perimeter reaches a maximum and then decreases to zero, leading to an almost negligible surface diffusion at large θ and zero diffusion flux of adatoms at θ → 1 .
Let us now consider the growth kinetics of a NW ensemble starting from the very beginning of growth. Re-formulating the results of Refs. [28,29] in terms of the coverage θ rather than c R 2 N , which gives the coverage only at c R 2 N ≪ 1 , we have
d θ L d H = a θ + L u 2 − a θ − L u ,   a θ + L u ≤ 1 ,
d θ L d H = 1 ,     a θ + L u > 1 ,
d L d H = a + Λ u θ 2 − a θ − L u ,   a θ + L u ≤ 1 ,
d L d H = a + Λ u θ ,   a θ + L u > 1 .
Here, H = v t is the 2D equivalent deposition thickness, and the parameters a and Λ have the same meaning as before. According to Equations (11), the volume of NWs per unit surface area increases due to the direct impingement of material onto the droplet surface (for VLS NWs with a > 1 ) or flat tops (for catalyst-free SAE or SI NWs with a = 1 ), and material collection on the sidewall surfaces of NWs having the total surface area L u . The multiplying factor 2 − a θ − L u describes the contribution from the re-emitted species. This model applies when desorption of material (group III adatoms for III-V or III-nitride NWs) from the NW sidewalls is negligible, in which case all adatoms either diffuse to the NW top and contribute into the axial growth or get incorporated into the radial shell growing around each NW. Such growth occurs until a NW ensemble start to collect the entire flux of the arriving material. After this moment of time, corresponding to a θ + L u = 1 , the growth kinetics is governed by Equation (1). Equations (12), which actually apply for L ≥ Λ , with Λ as the incorporation-limited diffusion length of adatoms on the NW sidewalls, show that the axial NW growth rate decreases with its length due to a lower fraction of re-emitted species landing on each NW and becomes a + Λ u / θ   after saturation, as given by Equation (4). These growth equations are written in the case of VPE [29], where
a = 2 1 + c o s β ,
with β as the contact angle of a catalyst droplet resting on the NW top ( β = 0 for non-VLS NWs). In the directional MBE technique, the growth equations are modified [28]. In particular, the parameter a depends on the droplet contact angle β and the group III beam angle with respect to the substrate normal α [33].
For both VPE and MBE NW growths, Equations (11) and (12) yield the evolution of the coverage with the average NW length in the form
d θ d L = u 1 − Λ / L a + Λ u / θ ,   a θ + L u ≤ 1 .
We start the analysis with the simplest NW growth with negligible surface diffusion of adatoms, corresponding to Λ = b = 0 . This case gives the lowest axial NW growth rate and hence the highest coverage compared to the general case with surface diffusion at b > 0 . Using d θ / d L = ( d θ / d R ) d R / d L = u d R / d L , Equation (14) yields d R / d L = 1 / a for the NW radius. Therefore, R scales linearly with the NW length:
R = R 0 + L a ,
where R 0 is the initial NW radius at L = 0 . Using this in Equation (5), the surface coverage is obtained in the form
θ = 1 − e x p − c R 0 + L / a 2 N , a θ + L u ≤ 1
The increase of coverage is described by this equation until the beginning of the late growth stage, which is governed by Equations (2) and (4). At Λ = b = 0 , these equations give the linear asymptotes θ L = θ * L * + H − H * and L = L * + a ( H − H * ) , where H * is the deposition thickness at which the total flux per NW saturates at its maximum level. Therefore, further evolution of the surface coverage with the deposition thickness or NW length is given by
θ = θ * L * + H − H * L * + a ( H − H * ) = θ * L * + ( L − L * ) / a L ,   a θ + L u > 1 ,
giving the asymptotic coverage θ ∞ = 1 / a at H → ∞ or L → ∞ .
Figure 4 shows the evolution of the coverage with the NW length, obtained from Equation (16) at a fixed initial NW radius of 30 nm and a fixed inter-NW distance of 500 nm, for cylindrical NWs with c = π . The asymptotic coverage   θ ∞ = 1 / a at Λ = 0 depends only on a , regardless of the NW density, initial radius and the details of growth kinetics. From Figure 4, only NWs with flat tops ( a = 1 ) will fully coalesce. VLS NWs with a small a = 2 will partly merge, but never coalesce into continuous film. NWs with larger a will remain separated at any time or for any length. The exponential increase of coverage given by Equation (16) to a much slower increase given by Equation (17) occurs at short NW lengths (from 190 to 270 nm in this example). Therefore, the slow increase of coverage in the late stage of NW growth provides a good approximation for nearly all growth times excluding the initial stage before the saturation of the total flux received by a NW ensemble. This observation also holds in the general case with surface diffusion of adatoms ( b > 0 ), because the diffusion-induced terms in the growth equations become negligible compared to the direct impingement at small enough u .
Most catalyst droplets on top of untapered VLS NWs have the droplet contact angles larger than 900 [4]. For example, the contact angles of Ga droplets on top of untapered self-catalyzed GaAs NWs typically range from 90o to 140o, depending on the Ga/As flux ratio [38,39]. Figure 5 shows the asymptotic coverage of the substrate surface as a function of the droplet contact angle β for VLS NWs grown by different epitaxy techniques. The dependence a ( β ) for VPE growth method is obtained from Equation (13). For the directional MBE growth, it is calculated using the expressions of Ref. [33] at a beam angle of 45o. In this example, merging of NWs is possible only in VPE growth for small enough contact angles β < 103o. As mentioned above, VLS growth without surface diffusion of adatoms yields the maximum possible coverage, because the axial NW growth rate increases and the radial NW growth rate decreases for larger diffusion lengths Λ . Therefore, VLS NWs with larger contact angles on their tops will remain separated at infinite length or growth time, and the asymptotic coverage will be less than in Figure 5 for VLS NWs with surface diffusion.
Presenting θ as a function of R , Equation (14) with arbitrary Λ > 0 takes the form
d R d L = 1 − Λ / L a + Λ u / θ .
Using Equations (5) for θ and u and integrating, we obtain the analytic relationship between the NW radius (or the coverage θ ) and length
L Λ − 1 − l n L Λ = a R − R 0 Λ + l n θ θ 0 .
Here, θ is given by Equation (5), R 0 is the NW radius at L = Λ , and θ 0 = 1 − e x p − c R 0 2 N is the initial coverage at R = R 0 . At Λ → 0 , this is reduced to the linear R ( L ) dependence given by Equation (15).
For high temperature growths with enhanced desorption from the NW sidewalls, the NW top collects material from the upper NW section of length Λ = Λ i n c , while the desorption-limited diffusion length equals Λ d e s [29,37]. The radial NW growth can occur only when Λ d e s > Λ i n c [29], In this case, Equation (14) rewrites [29]
d θ d L = u Λ d e s − Λ i n c L ( a + Λ u / θ ) ,   a θ + L u ≤ 1 .
Presenting θ as a function of R and repeating the above procedure, the analytic relationship between the NW length and coverage modifies to
L Λ d e s = θ θ 1 γ e x p γ a R − R 1 Λ i n c , γ = 1 Λ d e s / Λ i n c − 1 .
Here, R 1 is the NW radius at L = Λ d e s , and θ 1 = 1 − e x p − c R 1 2 N is the coverage at R = R 1 .
SI GaN NWs on different substrates, emerging from the Volmer-Weber 3D islands, usually have very high surface density, up to 1011 cm-2 [6,7,11,32,40], which is why the coalescence process starts quite early. The typical base radius of the initial 3D islands equals 5 nm [7,32,37]. Desorption of Ga adatoms from the NW sidewalls should be present in high-temperature MBE growth of GaN NWs at around 800 oC [32]. According to Equation (21), in the initial stage of growth of separated GaN NWs, corresponding to γ a ( R − R 1 ) / Λ i n c ≪ 1 and c R 2 N ≪ 1 , the NW length scales as a power law of its radius:
L Λ d e s ≅ R R 1 2 γ
with the power exponent γ = 1.23 according to Ref. [7].
Figure 6 shows the coverage versus L for a model system with R 0 = 5 nm, Λ d e s = 40 nm [7,32], γ = 1.23 [7] and c = 3 3 / 2 for hexahedral NWs, at three different surface densities N = 1010, 5 × 1010 and 1011 cm-2, corresponding to the NW separations P = 100, 45 and 32 nm. The curves are obtained from Equation (21) at a = 1 in the early stage and from numerical solution of Equations (2) and (4) after the flux saturation. These curves show the importance of the NW surface density (determined in the nucleation stage of 3D islands in the case of SI GaN NWs) in the kinetics of the coalescence process. The low surface density NWs ( N = 1010 cm-2, P = 100 nm) remain separated even at a large L of 2000 nm, although they will finally coalesce due to a = 1 . The higher surface density NWs start to merge in the early stage of growth. The highest surface density NWs ( N = 1011 cm-2, P = 32 nm) are almost fully coalesced at L = 2000 nm.

4. Conclusions

To summarize, the developed model predicts the following picture of the NW coalescence process. Any VLS NWs with the enhanced material collection efficiencies at their tops such that a > 1 cannot fully coalesce even without any surface diffusion of adatoms from the NW sidewalls (at b = 0 ), unless their droplets are removed ( a = 1 ) or the direction of the adatom diffusion flux is reversed ( b < 0 ). However, merging or partial coalescence of VLS NWs is entirely possible at low enough a and b according to Figure 3. Catalyst-free SAE or SI NWs with a = 1 can coalesce into continuous film regardless of the strength of surface diffusion, because the adatom diffusion flux is disabled for large enough coverages θ . NWs with larger a and b will reach a lower asymptotic coverage. This corresponds to a lower coverage and higher aspect ratio of denser NWs and material systems with longer diffusion lengths of adatoms on their sidewalls. The coalescence process usually occurs in the late stage of growth where a NW ensemble collects all material sent from vapor. However, our model is capable of describing the NW growth kinetics from the early stage where NWs grow separately at low surface coverage or only start to merge. The model is purely geometrical and contains the characteristics of a concrete material system and epitaxy technique only in the control parameters a (the droplet contact angle for VLS NWs and the beam direction in MBE) and b (the pitch of a NW array and the diffusion length of adatoms on the NW sidewalls). Therefore, it can be applied to a wide range of material systems and growth technologies, including VLS and catalyst-free NWs of III-V, III-nitride and oxide materials and elemental semiconductors. The model is not even restricted to vertical NWs and can treat any 3D surface islands which are able to change their aspect ratio in the course of growth. We now plan to consider the time-dependent coalescence process in different ensembles of NWs and other 3D structures and epitaxy techniques from the viewpoint of the obtained results. In particular, we plan to investigate the coalescence process in dense ensembles of SI GaN NWs grown by plasma-assisted MBE on Si substrates.

Funding

This research was supported by the research grant of St. Petersburg State University (ID 94033852).

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. McIntyre, P. C.; Fontcuberta i Morral, A. Semiconductor nanowires: to grow or not to grow? Mater. Today Nano 2020, 9, 100058. [Google Scholar] [CrossRef]
  2. Wagner, R. S.; Ellis, W. C. Vapor-liquid-solid mechanism of single crystal growth. Appl. Phys. Lett. 1964, 4, 89. [Google Scholar] [CrossRef]
  3. Colombo, C.; Spirkoska, D.; Frimmer, M.; Abstreiter, G.; Fontcuberta i Morral, A. Ga-assisted catalyst-free growth mechanism of GaAs nanowires by molecular beam epitaxy, Phys. Rev. B 2008, 77, 155326. [Google Scholar] [CrossRef]
  4. Dubrovskii, V. G.; Glas, F. Vapor–liquid–solid growth of semiconductor nanowires, In: N. Fukata, R. Rurali (eds) Fundamental Properties of Semiconductor Nanowires. Springer, Singapore, 2020. [CrossRef]
  5. Tomioka, K.; Ikejiri, K.; Tanaka, T.; Motohisa, J.; Hara, S.; Hiruma, K.; Fukui, T. Selective-area growth of III-V nanowires and their applications, Mater. Res. 2011, 26, 2127. [Google Scholar]
  6. Sanchez-Garcia, M. A.; Calleja, E.; Monroy, E.; Sanchez, F. J.; Calle, F.; Muñoz, E.; Beresford, R. J. The effect of the III/V ratio and substrate temperature on the morphology and properties of GaN- and AlN-layers grown by molecular beam epitaxy on Si(1 1 1). J. Cryst. Growth 1998, 183, 23. [Google Scholar] [CrossRef]
  7. Dubrovskii, V. G.; Consonni, V.; Geelhaar, L.; Trampert, A.; Riechert, H. Scaling growth kinetics of self-induced GaN nan-owires. Appl. Phys. Lett. 2012, 100, 153101. [Google Scholar] [CrossRef]
  8. Glas, F. Critical dimensions for the plastic relaxation of strained axial heterostructures in free-standing nanowires. Phys. Rev. B 2006, 74, 121302 (R). [Google Scholar] [CrossRef]
  9. Chuang, L. C.; Moewe, M.; Chase, C.; Kobayashi, N. P.; Chang-Hasnain, C.; Crankshaw, S. Critical diameter for III-V nanowires grown on lattice-mismatched substrates. Appl. Phys. Lett. 2007, 90, 043115. [Google Scholar] [CrossRef]
  10. Dubrovskii, V. G.; Sibirev, N. V.; Zhang, X.; Suris, R. A. Stress-driven nucleation of three-dimensional crystal islands: from quantum dots to nanoneedles. Cryst. Growth Des. 2010, 10, 3949. [Google Scholar] [CrossRef]
  11. Kaganer, V. M.; Fernández-Garrido, S.; Dogan, P.; Sabelfeld, K. K.; Brandt, O. Nucleation, growth, and bundling of GaN nanowires in molecular beam epitaxy: Disentangling the origin of nanowire coalescence. Nano Lett. 2016, 16, 3717. [Google Scholar] [CrossRef]
  12. Fernández-Garrido, S.; Kaganer, V. M.; Hauswald, C.; Jenichen, B.; Ramsteiner, M.; Consonni, V.; Geelhaar, L.; Brandt, O. Correlation between the structural and optical properties of spontaneously formed GaN nanowires: a quantitative evaluation of the impact of nanowire coalescence. Nanotechnology 2014, 25, 455702. [Google Scholar] [CrossRef] [PubMed]
  13. Li, Q.; Lin, Y.; J. Creighton, R.; Figiel, J. J.; Wang, G. T. Nanowire-templated lateral epitaxial growth of low-dislocation density nonpolar a-plane GaN on r-plane sapphire. Adv. Materials 2009, 21, 2416. [Google Scholar] [CrossRef]
  14. Conroy, M.; Zubialevich, V. Z.; Li, H.; Petkov, N.; Holmes, J. D.; Parbrook, P. J. Ultra-high-density arrays of defect-free AlN nanorods: A “space-filling” approach. ACS Nano 2016, 10, 1988. [Google Scholar] [CrossRef] [PubMed]
  15. Fan, S.; Zhao; S, Liu; X.; Mi, Z. Study on the coalescence of dislocation-free GaN nanowires on Si and SiOx. J. Vac. Sci. Tech. 2014, B 32, 02C114.
  16. Le, B. H.; Zhao, S.; Liu, X.; Woo, S. Y.; Botton, G. A.; Mi, Z. Controlled coalescence of AlGaN nanowire Arrays: An architecture for nearly dislocation-free planar ultraviolet photonic device applications. Adv. Materials 2016, 28, 8446. [Google Scholar] [CrossRef] [PubMed]
  17. Mandl, B.; Dey, A. W.; Stangl, J.; Cantoro, M.; Wernersson, L. E.; Bauer, G.; Samuelson, L.; Deppert,K. ; Thelander, C. Self-seeded, position-controlled InAs nanowire growth on Si: A growth parameter study. J. Cryst. Growth 2011, 334, 51. [Google Scholar] [CrossRef] [PubMed]
  18. Zhu, L.; Phillips, M. R.; Ton-That, C. Coalescence of ZnO nanowires grown from monodisperse Au nanoparticles. CrystEngComm 2015, 17, 4987. [Google Scholar] [CrossRef]
  19. Lotty, O.; Hobbs, R.; O’Regan, C.; Hlina, J.; Marschner, C.; O’Dwyer, C.; Petkov, N; Holmes, J. D. Self-seeded growth of germanium nanowires: Coalescence and Ostwald ripening. Chemistry of Materials 2013, 25, 215. [Google Scholar] [CrossRef]
  20. Sabelfeld, K.; Kablukova, E. Stochastic model and simulation of growth and coalescence of spontaneously formed GaN nanowires in molecular beam epitaxy. Comput. Mater. Sci. 2016, 125, 284. [Google Scholar] [CrossRef]
  21. Wu, Z.; Yang, X.; Wang, Z. Size effect on the spontaneous coalescence of nanowires. Nanotechnology 2019, 30, 245601. [Google Scholar] [CrossRef]
  22. Chi, C.-Y.; Chang, C.-C.; Hu, S.; Yeh, T.-W.; Cronin, S. B.; Dapkus, P. D. Twin-free GaAs nanosheets by selective area growth: Implications for defect-free nanostructures. Nano Lett. 2013, 13, 2506. [Google Scholar] [CrossRef]
  23. Yuan, X.; Pan, D.; Zhou, Y.; Zhang, X.; Peng, K.; Zhao, B.; Deng, M.; He, J.; Tan, J. H. H.; Jagadish, C. Selective area epitaxy of III–V nanostructure arrays and networks: Growth, applications, and future directions. Appl. Phys. Reviews 2021, 8, 021302. [Google Scholar] [CrossRef]
  24. Ramanandan, S. P.; Tomić, P.; Morgan, N. P.; Giunto, A.; Rudra, A.; Ensslin, K.; Ihn, T.; Fontcuberta i Morral, A. Coherent hole transport in selective area grown Ge nanowire networks. Nano Lett. 2022, 22, 4269. [Google Scholar] [PubMed]
  25. Morgan, N. P.; Dubrovskii, V. G.; Stief, A.-K.; Dede, D.; Sanglé-Ferriere, M.; Rudra, A.; Piazza, V.; Fontcuberta i Morral, A. From layer-by-layer growth to nanoridge formation: Selective area epitaxy of GaAs by MOVPE. Cryst. Growth Des. 2023, 23, 5083. [Google Scholar] [CrossRef] [PubMed]
  26. Li, Q.; Jiang, H.; Lau, K. M. Coalescence of planar GaAs nanowires into strain-free three-dimensional crystals on exact (001) silicon. J. Cryst. Growth 2016, 454, 19. [Google Scholar] [CrossRef]
  27. Oehler, F.; Cattoni, A.; Scaccabarozzi, A.; Patriarche, J.; Glas, F.; Harmand, J. C. Measuring and modeling the growth dynamics of self-catalyzed GaP nanowire arrays. Nano Lett. 2018, 18, 701. [Google Scholar] [CrossRef]
  28. Dubrovskii, V. G. Theory of MBE growth of nanowires on reflecting substrates. Nanomaterials 2022, 12, 253. [Google Scholar] [CrossRef]
  29. Dubrovskii, V. G. Theory of MOCVD growth of III-V nanowires on patterned substrates. Nanomaterials 2022, 12, 2632. [Google Scholar]
  30. Cirlin, G. E.; Dubrovskii, V. G.; Sibirev, N. V.; Soshnikov, I. P.; Samsonenko, Yu. B.; Tonkikh, A. A.; Ustinov, V. M. The diffusion mechanism in the formation of GaAs and AlGaAs nanowhiskers during the process of molecular-beam epitaxy. Semiconductors 2005, 39, 557. [Google Scholar] [CrossRef]
  31. Johansson, J.; Magnusson, M. H. From diffusion limited to incorporation limited growth of nanowires, J. Cryst. Growth 2019, 525, 125192. [Google Scholar]
  32. Consonni, V.; Dubrovskii, V. G.; Trampert, A.; Geelhaar, L.; Riechert, H. Quantitative description for the growth rate of self-induced GaN nanowires. Phys. Rev. B 2012, 85, 155313. [Google Scholar]
  33. Glas, F. Vapor fluxes on the apical droplet during nanowire growth by molecular beam epitaxy. Phys. Status Solidi B 2010, 247, 254. [Google Scholar] [CrossRef]
  34. Dubrovskii, V. G.; Sibirev, N. V.; Cirlin, G. E.; Bouravleuv, A. D.; Samsonenko, Yu. B.; Dheeraj, D. L.; Zhou, H. L.; Sartel, C.; Harmand, J. C.; Patriarche, G.; Glas, F. Role of non-linear effects in nanowire growth and crystal phase. Phys. Rev. B 2009, 80, 205305. [Google Scholar] [CrossRef]
  35. Kolmogorov, A. N. On the statistical theory of the crystallization of metals. Bull. Acad. Sci. URSS (Cl. Sci. Math. Nat.) 1937, 3, 355. [Google Scholar]
  36. Hertenberger, S.; Rudolph, D.; Bichler, M.; Finley, J. J.; Abstreiter, G.; Koblmüller, G. Growth kinetics in position-controlled.
  37. and catalyst-free InAs nanowire arrays on Si(111) grown by selective area molecular beam epitaxy. J. Appl. Phys. 2010, 108, 114316. 2010; 108, 114316.
  38. Dubrovskii, V. G. Modeling catalyst-free growth of III-V nanowires: empirical and rigorous approaches. Nanomaterials 2023, 13, 1253. [Google Scholar] [CrossRef]
  39. Panciera, F.; Baraissov, Z.; Patriarche, G.; Dubrovskii, V. G.; Glas, F.; Travers, L.; Mirsaidov, U.; Harmand, J. C. Phase selection in self-catalyzed GaAs nanowires. Nano Lett. 2020, 20, 1669. [Google Scholar] [CrossRef]
  40. Kim, W.; Dubrovskii, V. G.; Vukajlovic-Plestina, J.; Tütüncüoglu, G.; Francaviglia, L.; Güniat, L.; Potts, H.; Friedl, M.; Leran, J.-B.; Fontcuberta i Morral, A. Bistability of contact angle and its role in achieving quantum-thin self-assisted GaAs nanowires.
  41. Nano Lett. 2018, 18, 49.
  42. Olszewski, K.; Sobanska, M.; Dubrovskii, V. G.; Leshchenko, E. D.; Wierzbicka, A.; Zytkiewicz, Z. R. Geometrical selection of GaN nanowires grown by plasma-assisted MBE on polycrystalline ZrN layers. Nanomaterials 2023, 13, 2587. [Google Scholar] [CrossRef] [PubMed]
Preprints 85823 g001
Figure 1. Stages of the growth process resulting in (a) isolated NWs, (b) partly merged NWs, and (c) fully coalesced NWs with θ = 1 . (d) Mosaic structure of merged SI GaN NWs [7] in stage (b). Similar picture of the coalescence process applies to 3D islands in SAE, where the islands emerge with a small size and a given density (determined in the nucleation step) but then merge together to fill a growth template.
Figure 1. Stages of the growth process resulting in (a) isolated NWs, (b) partly merged NWs, and (c) fully coalesced NWs with θ = 1 . (d) Mosaic structure of merged SI GaN NWs [7] in stage (b). Similar picture of the coalescence process applies to 3D islands in SAE, where the islands emerge with a small size and a given density (determined in the nucleation step) but then merge together to fill a growth template.
Preprints 85823 g001
Preprints 85823 g002
Figure 2. Normalized perimeters per unit surface area u   and per surface area of the NW tops u / θ versus θ (the solid lines). The dashed lines show the same functions for isolated NWs.
Figure 2. Normalized perimeters per unit surface area u   and per surface area of the NW tops u / θ versus θ (the solid lines). The dashed lines show the same functions for isolated NWs.
Preprints 85823 g002
Preprints 85823 g003
Figure 3. Asymptotic coverage θ ∞ versus a at three different b shown in the legend (the solid lines), obtained from Equation (10). The dashed lines show the approximations given by Equation (9) for isolated NWs, which underestimate the coverage for any a at b > 0 . At b = 0 , both models give the same result θ ∞ = 1 / a . The full coalescence occurs only at a = 1 . The asymptotic coverage at a fixed a > 1 decreases with b . The shaded zone corresponds to partly merged NWs at θ ∞ > θ m = 0.39 .
Figure 3. Asymptotic coverage θ ∞ versus a at three different b shown in the legend (the solid lines), obtained from Equation (10). The dashed lines show the approximations given by Equation (9) for isolated NWs, which underestimate the coverage for any a at b > 0 . At b = 0 , both models give the same result θ ∞ = 1 / a . The full coalescence occurs only at a = 1 . The asymptotic coverage at a fixed a > 1 decreases with b . The shaded zone corresponds to partly merged NWs at θ ∞ > θ m = 0.39 .
Preprints 85823 g003
Preprints 85823 g004
Figure 4. Coverage of a substrate surface with NWs that start from a fixed radius of 30 nm and have a fixed surface density corresponding to an average separation of 500 nm, at different a shown in the legend. NWs will fully coalesce only at a = 1 , as in Figure 3. NWs will partly merge at a = 2 . Other NWs will remain separated for any length.
Figure 4. Coverage of a substrate surface with NWs that start from a fixed radius of 30 nm and have a fixed surface density corresponding to an average separation of 500 nm, at different a shown in the legend. NWs will fully coalesce only at a = 1 , as in Figure 3. NWs will partly merge at a = 2 . Other NWs will remain separated for any length.
Preprints 85823 g004
Preprints 85823 g005
Figure 5. Asymptotic coverage as a function of the droplet contact angle for VLS NWs grown by VPE and MBE at a beam angle of 45o. Only VPE NWs with small contact angles β < 103o will partly merge.
Figure 5. Asymptotic coverage as a function of the droplet contact angle for VLS NWs grown by VPE and MBE at a beam angle of 45o. Only VPE NWs with small contact angles β < 103o will partly merge.
Preprints 85823 g005
Preprints 85823 g006
Figure 6. Evolution of the coverage with the NW length at three different NW surface densities corresponding to the average separations P shown in the legend. At a low surface density of 1010 cm-2 ( P = 100 nm), NWs do not even start to merge at an average length of 2000 nm, while the higher surface density NWs with N = 5 × 1010 and 1011 cm-2 merge early. At the highest surface density of 1011 cm-2, the NW film almost reaches continuity at L = 2000 nm.
Figure 6. Evolution of the coverage with the NW length at three different NW surface densities corresponding to the average separations P shown in the legend. At a low surface density of 1010 cm-2 ( P = 100 nm), NWs do not even start to merge at an average length of 2000 nm, while the higher surface density NWs with N = 5 × 1010 and 1011 cm-2 merge early. At the highest surface density of 1011 cm-2, the NW film almost reaches continuity at L = 2000 nm.
Preprints 85823 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated