Abstract
The growth of inclined GaAs nanowires (NWs) during molecular beam epitaxy (MBE) on the rotating substrates is studied. The growth model provides explicitly the NW length as a function of radius, supersaturations, diffusion lengths and the tilt angle. Growth experiments are carried out on the GaAs(211)A and GaAs(111)B substrates. It is found that 20° inclined NWs are two times longer in average, which is explained by a larger impingement rate on their sidewalls. We find that the effective diffusion length at 550°C amounts to 12 nm for the surface adatoms and is more than 5,000 nm for the sidewall adatoms. Supersaturations of surface and sidewall adatoms are also estimated. The obtained results show the importance of sidewall adatoms in the MBE growth of NWs, neglected in a number of earlier studies.
Keywords:
Inclined GaAs nanowires; Molecular beam epitaxy; Surface diffusionIntroduction
A rapidly growing interest in selfstanding NWs of III–V compound semiconductors is explained by an interesting physics of their growth [16], crystal structure [6,7] transport [8] and optical [9] phenomena as well as a variety of promising applications in nanoelectronics [8] and nanophotonics [9,10]. III–V NWs with radii of the order of 10 nm and length up to ten micrometers are usually fabricated by metal organic chemical vapor deposition (MOCVD) [1,2] or MBE [37] via the socalled vapor–liquid–solid (VLS) mechanism [11] on the substrates activated by a metal (Au) growth catalyst. Due to their ability to accommodate strain in two dimensions, NWs are ideal candidates for monolithic integration of dissimilar semiconductor materials, e.g., of III–V compounds on Si [12,13]. For the controlled production of NWs with the desired morphology and crystal structure, it is important to understand major kinetic processes driving NW growth at given set of deposition conditions. Furthermore, theoretical and experimental investigations into NW formation mechanisms can provide important information on many kinetic characteristics (e.g., supersaturations, diffusion lengths and surface energies [15,1422]) that are otherwise not easy to measure or even define theoretically.
Since the dominant growth direction of NWs is , most growth experiments [16] are carried out on the (111) oriented surfaces, with NWs being perpendicular to the substrate. As regards the growth mechanisms of such NWs, semiconductor material is transferred to the drop (seated at the NW top) by different kinetic pathways: the direct impingement onto the drop surface and the surface diffusion of adatoms that first impinge the sidewalls and substrate [16,1422]. The diffusioninduced contribution into the overall growth rate is always dominant in MBE [36,14,22]. For very thin NWs, it is important to consider the Gibbs–Thomson (GT) effect of elevation of chemical potential caused by the curvature of the drop surface [2,15]. In MBE case, the flux directly impinging the sidewalls increases with the incident angle of the beam. The use of substrate orientation other than (111), resulting in the formation of inclined NWs with varying tilt angle (and consequently the incident angle of the beam), can therefore provide an additional parameter to alter the NW morphology. The incident angle of the flux impinging the tilted NW is, however, changing in time if the substrate is rotating, so a proper averaging should be introduced to generalize the existing models. The use of highindex substrates, implemented earlier, e.g., for the growth of InGaAs/GaAs quantum wells [23], can be of particular importance in connection with the phase perfection: as demonstrated in Refs. [7,24], the use of particular highindex GaAs substrates for the Auseeded MBE growth of GaAs NWs produces stacking faultfree zincblende structure.
In this work, we report on theoretical and experimental investigation into the Auassisted MBE of inclined GaAs NWs. A theoretical model of Refs. [6,15] is developed to describe the growth of inclined NWs and to find explicitly the dependence of NW growth rate as a function of tilt angle at given set of deposition conditions. We then carry out MBE growth experiments on the GaAs(211)A and GaAs(111)B substrates. From the analysis of scanning electron microscopy (SEM) images of different samples, we plot the lengthradius curves and fit them by theoretical dependences. This enables to deduce some important kinetic parameters of NW growth, in particular, the effective supersaturations and diffusion lengths on different GaAs surfaces.
Theoretical Model
The model of inclined NW is sketched in Fig. 1. We consider a single NW growing in a stationary mode in z direction making the tilt angle φ to the substrate normal. Neither lateral growth nor shadowing effect is taken into account. The model parameters include the impingement flux J, the incident angle of molecular beam α (in the case of III–V compounds, the growth rate is assumed as being limited by the incorporation of group III element so that α relates to the group III beam), the contact angle of the drop β and the drop surface energy γ. Under the standard assumption of a low concentration of group V element (As) in the drop [2,4,25], the value of γ must be between the surface energy of pure liquid group III element (Ga) and Au at the growth temperature T. Effective diffusion lengths on the substrate and sidewalls, limited either by desorption or incorporation to a growing surface layer [25], are denoted as λ_{s} and λ_{f}, respectively. The quantities θ_{v}θ_{s}θ_{f} and θ_{l} denote the activities (the effective supersaturations) of group III element in the vapor (v), surface adatom (s), sidewall adatom (f) and liquid (l) phases, with the usual definition of θ_{i} = exp(μ_{i}/k_{B}T), where μ_{i} is the chemical potential in phase i (measured, e.g., with respect to the solid phase) and k_{B} is the Boltzmann constant [6].
Figure 1. Illustration of the growth model with the parameters described in the text
The sidewall impingement rate must be corrected for the effective incident angle of the beam to the inclined NW. Below we consider the general case of rotating substrate, schematized in Fig. 2. The NW growth direction is given by the radius vector while the direction of the flux is parallel to the vector For the incident angle of the beam to the NWξ, this yields:
Figure 2. Direction parallel to the flux and the momentum growth direction of NW at time t (with respect to stationary coordinates x and y) defined by the tilt angle ϕ and the angle ψ = ωt, with ω as the angular velocity of substrate rotation
The activities of adatom phases [6,15] can be now put in the form
Here, σ_{s}, σ_{f} are the elementary areas of substrate and sidewall surfaces, and τ_{s}, τ_{f} are the corresponding adatom lifetimes. The quantity denotes the mean value of sinξ averaged over the substrate revolution. Obviously,
where cosξ is given by Eq. 1. Our numerical analysis shows that Eq. 3 can be well approximated (with less than 7% error) by the simplified formula where, in view of Eq. 1, With this approximation, we get
the expression used hereinafter for the effective incident angle to the inclined NWs. Obviously, Eq. 4 [as well as general Eqs. 1 and 3] is reduced to trivial identity at φ = 0, i.e., the growth of straight NWs is not influenced by the substrate rotation.
The detailed analysis of Refs. [6,15] shows that the exact solution for the stationary NW growth rate has the form
Here, L is the NW length, H = Vt is the effective deposition thickness at time t, V = JΩ_{s} cosα is the deposition rate (with Ω_{s} as the elementary volume in the solid phase) and is the effective diffusion length on the sidewalls (with D_{f} as the corresponding diffusion coefficient). The A term describes the direct impingement onto the drop surface [26], the desorption from the drop and the growth of surface layer. Since in our experiments the inequality β > π/2 + ξ [26] holds for most ξ during the substrate rotation, the corresponding averaging is almost independent on the incident angle ξ:
The desorption term contains the standard GT modification of liquid activity θ_{l}[2,6], with R_{GT} = (2γΩ_{l} sinβ)/(k_{B}T) as the characteristic GT radius and Ω_{l} as the elementary liquid volume. The desorption rate is inversely proportional to the ratio of activities in the vapor and infinitely large liquid alloy, The quantity ɛ = V_{s}/V accounts for the substrate growth, with V_{s} as the growth rate on the nonactivated surface [3].
The second, Ldependent term in the right hand side of Eq. 5 gives two diffusioninduced contributions, one originating from the adatoms impinging the sidewalls and migrating to the drop and another caused by the adatoms that first impinge the substrate and then diffuse to the drop along the sidewalls. The coefficients B (describing the sidewall adatoms) and C (describing the substrate adatoms) in the case of MBE are given by [6]
Here, and are the effective supersaturations of sidewall and surface adatoms with respect to the infinite liquid alloy, δ(R/λ_{s}) = K_{1}(R/λ_{s})/K_{0}(R/λ_{s}), is the effective diffusion length on the substrate surface (with D_{s} as the corresponding diffusion coefficient) and K_{n}(x) are the modified Bessel function of the second kind of order n. Further, the functions U(L/λ_{f}) and U′(L/λ_{f}) in Eq. 5 are defined as follows
with
At given MBE growth conditions (T, V, α, V/III fluxes ratio), the vapor–solid chemical potential and consequently the vapor supersaturation is well defined [25], but it generally tells nothing about the quantity θ_{vl} entering Eq. 6 for A. If, however, the growth temperature is low enough, we can safely neglect the desorption term in Eq. 6 (perhaps excluding very thin NWs) and eliminate unknown θ_{vl}. In this case, with the known material constants (providing the GT radius R_{GT}), contact angle β and tilt angle ϕ, the measured L(H) or L(R) curves of inclined NWs can be fitted by four parameters, two diffusion lengths λ_{s}λ_{f} and two supersaturations θ_{fl}θ_{sl}. At L/λ_{f} ≪ 1 and λ_{s}/λ_{f} ≪ 1, the dependence on λ_{f} disappears. The nonvanishing terms at L/λ_{f} → 0 and λ_{s}/λ_{f} → 0 (yielding also τ_{s}/τ_{f} → 0) reduce Eqs. 5, 8, 9 to the nonlinear growth equation of the form [27]
with coefficients
Thus, our model provides the exact solution for the NW growth rate with the tilt angle as the control parameter. This solution contains, however, a number of unidentified quantities that can be found only from the direct comparison with experiment.
Growth Experiments
Our growth experiments were carried out with a Riber 32 MBE setup equipped with solid sources supplying monomers of Ga and tetramers of As_{4}. Growth was performed on the GaAs (111)B and GaAs(211)A substrates. During the growth, the substrates rotation was applied. Incident angle of Ga flux α amounted to 16.7° to the normal. Substrate surfaces were deoxidized at 620°C, and a 30 nm GaAs buffer layers were grown to achieve atomically flat surface. Then, the substrate temperature was decreased and stabilized to T = 550°C for the NW growth on the both types of substrates. For catalyst deposition, Au source installed directly into the growth chamber in a regular effusion cell was used. This configuration enables to deposit the catalyst on chemically clean surface, and at the same time to control the substrate temperature and monitor the deposition process with reflection highenergy electron diffraction. An amount of Au equivalent to a uniform layer of ~1 nm was deposited on the substrate surface to promote the NW growth. This procedure resulted in the formation of Au droplets alloyed with the substrate constituents which could activate the NW growth. Since the growth temperature of 550°C is much higher than the lowest eutectic temperatures of Au–Ga alloy (339.4°C and 348.9°C [28]), the drop must remain liquid during the growth. For all samples, the nominal growth rate V was fixed to 0.2 nm/s, the V/III flux ratio was equal 3 and the effective deposition thickness was fixed to 360 nm with the corresponding deposition time of 30 min.
The morphology of asgrown GaAs NW ensembles was investigated using fieldemission SEM technique. Figures 3 and 4 present typical SEM images of GaAs NWs grown on the GaAs(111)B and GaAs(211)A substrates, respectively, under identical deposition conditions described hereinabove. It is seen that all NWs have almost uniform radius from base to top, which agrees with the results of Ref. [29] where the pronounced lateral growth [18,19] was observed only below 500°C. This justifies the assumption of NW elongation at R = const. As expected, the NWs grown on the (111)B substrate are perpendicular to the surface. From the analysis of plan and top view SEM images, the tilt angle of NWs on the (211)A substrate varies from 15 to 30°, with the statistical average around 20°, which equals the angle between the (211) plane and crystallographic direction. This confirms the dominant growth direction of inclined NWs. The measured lengths of selected NWs and the average thickness of surface layer on the nonactivated substrates are shown in the figures. Statistical analysis of SEM images presented in Figs. 3 and 4 enables to plot out the lengthradius dependences shown by points in Fig. 5.
Results and Discussion
From Fig. 5, the length of all NWs is noticeably larger than the deposition thickness (360 nm), and is almost 6 times larger for longest 2 μm NWs, which proves the diffusioninduced character of NW growth discussed previously in Ref. [16,1422]. The solid lines in Fig. 5 are obtained from general expressions given by Eqs. 2, 4, 5–9 with the parameters summarized in Table 1. The dotted lines correspond to simplified Eqs. 10 and 11 with the same fitting parameters. According to the data of Ref. [25], the vapor supersaturation θ_{v} = exp (Δμ_{v}/k_{B}T) with respect to the GaAs(111)B substrate equals approximately 183 at T = 550°C,, so that the assumption of θ_{vl} → ∞ (i.e., negligible reevaporation from the drop) looks reasonable. At T = 550°C, the average contact angle of the drops β = 120° and γ = 1.0 J/m^{2} (corresponding to approximately 40% Ga concentration in the liquid Au–Ga alloy during the growth [15]), the GT radius R_{GT} equals 5.8 nm. With neglect of desorption, the value of ɛ in Eq. 6 is estimated as ɛ ≅ H_{s}/H = 327/360 = 0.91, i.e. only 9% of material is distributed in the NWs and 91% remains in a 2D surface layer growing concomitantly with the NWs. As follows from Fig. 5, the simplified growth equation at L/λ_{f} → 0 is a good approximation to the general expressions for the parameters considered. A small difference which can be seen for the longest NWs with L > 2 μm is explained by the reevaporation of some Ga adatoms from the sidewalls. Due to this desorption, the actual length becomes smaller than that predicted by the simplified formula where the reevaporation is neglected. As in Refs. [2,15,25], theoretical model predicts the nonmonotonous behavior of the L(R) curves, reaching their maxima due to the balance of the GT and the diffusioninduced contributions into the overall growth rate. The GT effect suppresses completely the growth of straight GaAs with R < 9 nm and inclined GaAs NWs with R < 15 nm. The obtained estimate of minimum radius for the straight GaAs NWs is consistent with the result of Ref. [2] in the case of Aucatalyzed MOCVD of straight InAs NWs (~8 nm at T = 425°C).
Table 1. Growth conditions and fitting parameters for different GaAs NWs
As follows from the results summarized in Table 1, the effective diffusion length on the substrate surfaces (limited by the incorporation into a growing surface layer) appears to be only 12 nm for the both substrates studied. This estimate is noticeably smaller than the previously obtained results of 25 nm (Ref. [22]) and 35 nm (Ref. [14]) at 560°C. Such difference is most probably explained by the simplified growth equation used in Ref. [14,22], i.e., dL/dH = a_{0} instead of Eq. 10, resulting in the neglect of sidewall adatoms. It is noteworthy that the nonlinear terms in L/λ_{s} in Eqs. 10, 11 contain the contributions from sidewall adatoms through the λ_{f}independent coefficients (λ_{s}/λ_{f})B in the corresponding Eqs. 11 for a_{1} and a_{2}. These contributions cancel exactly only at B = 0, i.e., for straight NWs (ϕ = 0) and the beam being strictly perpendicular to the substrate (α = 0). Otherwise, the diffusion of adatoms directly impinging the NW sidewalls plays an important role even at the initial stage of NW growth with L/λ_{f} ≪ 1, the effect overlooked in the number of recent studies [2,14,22]. Since our NWs are relatively short, the fit obtained from general expressions given by Eqs. 5–9 becomes independent on λ_{f} at λ_{f} ≥ 5,000 nm. We can therefore conclude that the effective diffusion length of Ga atoms on the NW sidewalls (which should be constructed from six equivalent facets in the case of zincblende NWs or their wurtzite counterparts [25]) is more than 5,000 nm. This result is consistent with previous estimates, e.g., 3,000 nm at 590°C in Ref. [4].
As regards the obtained estimates for the effective supersaturations, the first obvious conclusion is that the inequalities θ_{sl} > 1 and θ_{fl} > 1 yield positive (i.e., directed from base to top) diffusion fluxes at R_{GT}/R ≪ 1 for the surface and sidewall adatoms, because the adatom chemical potentials are larger than the chemical potential of infinite liquid alloy [6]. The corresponding flux becomes negative only for sufficiently thin NWs due to the GT effect. For the both cases considered, the supersaturation of sidewall adatoms is several times larger than that of surface adatoms, which is qualitatively consistent with the strong inequality λ_{f}/λ_{s} ≫ 1. The supersaturation of sidewall adatoms is noticeably larger for the inclined NWs (16.3 against 11.2 for the straight NWs), which is again explained by a larger impingement onto the tilted sidewalls. The corresponding differences in chemical potentials in the adatom and infinite liquid phases, obtained from the relationships θ_{sl} = exp (Δμ_{sl}/k_{B}T); θ_{fl} = exp (Δμ_{fl}/k_{B}T), equal 24.0–36.5 meV for the surface and 172–199 meV for the sidewall adatoms.
To sum up, our results show that the diffusion of adatoms that first impinge the sidewalls has a tremendous effect on the growth rate. First, since the diffusivity of surface adatoms during MBE is fundamentally limited by the growing surface layer, and their supersaturation is much lower than that of the sidewall adatoms, the coefficient C in Eq. 5 is usually much smaller than B and even much smaller than A. Therefore, the initial growth stage should be controlled by the direct impingement onto the drop surface (A), while the contribution from the sidewall adatoms (B) rapidly increases as the NW elongates. Second, our experimental data and theoretical fits demonstrate that the inclined NWs grow much faster than the straight ones: from Fig. 4, the GaAs NWs on the GaAs(211)A substrate are more than 2 times larger at the same radii and otherwise identical deposition conditions. Better analysis could be performed with the experimental length–time dependences, where different contributions at different growth stages would be more easily distinguished [30,31]. As yet, however, we do not have in hand such experimental data for the inclined GaAs NWs. As the NWs grow, the shadowing effect might become important at a certain critical length which can be easily estimated with given incident angle of the beam and the NW density. We now plan to consider these effects from the viewpoint of the obtained results.
Acknowledgments
This work was partially supported by the 111 Project (No. B07005), Program for Changjiang Scholars and Innovative Research Team in University (No. IRT0609), National Basic Research Program of China (No. 2010CB327600), Russian Federal Agency for Science and Innovation (Contract No. 02.740.11.0383), the scientific program of Russian Academy of Sciences “Fundamental aspects of nanotechnologies and nanomaterials” and few grants of Russian Foundation for Basic Research.
Open Access
This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
References

Seifert W, Borgstrom M, Deppert K, Dick KA, Johansson J, Larsson MW, Martensson T, Skold N, Svensson CPT, Wacaser BA, Wallenberg LR, Samuelson L:
J. Cryst. Growth. 2004, 272:211.
COI number [1:CAS:528:DC%2BD2cXhtVamtbnN]; Bibcode number [2004JCrGr.272..211S]
Publisher Full Text 
Fröberg LE, Seifert W, Johansson J:
Phys. Rev. B. 2007, 76:153401.
Bibcode number [2007PhRvB..76o3401F]
Publisher Full Text 
Dubrovskii VG, Cirlin GE, Soshnikov IP, Tonkikh AA, Sibirev NV, Samsonenko YuB, Ustinov VM:
Phys. Rev. B. 2005, 71:205325.
Bibcode number [2005PhRvB..71t5325D]
Publisher Full Text 
Harmand JC, Patriarche G, PéréLaperne N, MératCombes MN, Travers L, Glas F:
Appl. Phys. Lett.. 2005, 87:203101.
Bibcode number [2005ApPhL..87t3101H]
Publisher Full Text 
J. Cryst. Growth. 2006, 286:394.
COI number [1:CAS:528:DC%2BD28XjsFan]; Bibcode number [2006JCrGr.286..394P]
Publisher Full Text 
Dubrovskii VG, Sibirev NV, Cirlin GE, Bouravleuv AD, Samsonenko YuB, Dheeraj DL, Zhou HL, Sartel C, Harmand JC, Patriarche G, Glas F:
Phys. Rev. B. 2009, 80:066940. Publisher Full Text

Shtrikman H, PopovitzBiro R, Kretinin A, Heiblum M:
Nano Lett.. 2009, 9:215.
COI number [1:CAS:528:DC%2BD1cXhsFWqsrzO]; Bibcode number [2009NanoL...9..215S]
PubMed Abstract  Publisher Full Text 
Bjork MT, Ohlsson BJ, Sass T, Persson AI, Thelander C, Magnusson MH, Deppert K, Wallenberg LR, Samuelson L:
Appl. Phys. Lett.. 2002, 80:1058.
COI number [1:CAS:528:DC%2BD38Xpt1ChsA%3D%3D]; Bibcode number [2002ApPhL..80.1058B]
Publisher Full Text 
Xiang J, Lu W, Hu YJ, Wu Y, Yan H, Lieber CM:
Nature. 2006, 441:489.
COI number [1:CAS:528:DC%2BD28XkvVyku7Y%3D]; Bibcode number [2006Natur.441..489X]
PubMed Abstract  Publisher Full Text 
Cirlin GE, Bouravleuv AD, Soshnikov IP, Samsonenko YuB, Dubrovskii VG, Arakcheeva EM, Tanklevskaya EM, Werner P:
Nanoscale Res. Lett.. 2010, 5:360.
COI number [1:CAS:528:DC%2BC3cXhsV2ls74%3D]; Bibcode number [2010NRL.....5..360C]
Publisher Full Text 
Appl. Phys. Lett.. 1964, 4:89.
COI number [1:CAS:528:DyaF2cXls1yhug%3D%3D]; Bibcode number [1964ApPhL...4...89W]
Publisher Full Text 
Chuang LC, Moewe M, Crankshaw S, Chase C, Kobayashi NP, ChangHasnain C:
Appl. Phys. Lett.. 2007, 90:043115.
Bibcode number [2007ApPhL..90d3115C]
Publisher Full Text 
Cirlin GE, Dubrovskii VG, Soshnikov IP, Sibirev NV, Samsonenko YuB, Bouravleuv AD, Harmand JC, Glas F:
Phys. Stat. Sol. RRL. 2009, 3:112.
COI number [1:CAS:528:DC%2BD1MXmvFKku78%3D]
Publisher Full Text 
Dubrovskii VG, Sibirev NV, Suris RA, Cirlin GE, Harmand JC, Ustinov VM:
Surf. Sci.. 2007, 601:4395.
COI number [1:CAS:528:DC%2BD2sXhtVOrt7rN]; Bibcode number [2007SurSc.601.4395D]
Publisher Full Text 
Dubrovskii VG, Sibirev NV, Cirlin GE, Soshnikov IP, Chen WH, Larde R, Cadel E, Pareige P, Xu T, Grandidier B, Nys JP, Stievenard D, Moewe M, Chuang LC, ChangHasnain C:
Phys. Rev. B. 2009, 79:205316.
Bibcode number [2009PhRvB..79t5316D]
Publisher Full Text 
Johansson J, Svensson CPT, Martensson T, Samuelson L, Seifert W:
J. Phys. Chem. B. 2005, 109:13567.
COI number [1:CAS:528:DC%2BD2MXlsFKhtbY%3D]
PubMed Abstract  Publisher Full Text 
Nanoscale Res. Lett.. 2007, 2:235.
COI number [1:CAS:528:DC%2BD2sXmvFylt70%3D]; Bibcode number [2007NRL.....2..235A]
Publisher Full Text 
J. Appl. Phys.. 2009, 105:114304.
Bibcode number [2009JAP...105k4304P]
Publisher Full Text 
Dubrovskii VG, Sibirev NV, Cirlin GE, Tchernycheva M, Harmand JC, Ustinov VM:
Phys. Rev. E. 2008, 77:031606.
COI number [1:STN:280:DC%2BD1czlvVWnsg%3D%3D]; Bibcode number [2008PhRvE..77c1606D]
Publisher Full Text 
J. Vac. Sci. Technol.. 2010, 28:329.
COI number [1:CAS:528:DC%2BC3cXkt1KmsLk%3D]

De Jong E, LaPierre RR, Wen JZ:
Nanotechnology. 2010, 21:045602.
Bibcode number [2010Nanot..21d5602D]
PubMed Abstract  Publisher Full Text 
Dubrovskii VG, Sibirev NV, Suris RA, Cirlin GE, Ustinov VM, Tchernysheva M, Harmand JC:
Semiconductors. 2006, 40:1075.
COI number [1:CAS:528:DC%2BD28Xpsl2nu78%3D]; Bibcode number [2006Semic..40.1075D]
Publisher Full Text 
Li Z, Wu J, Wang ZM, Fan D, Guo A, Li S, Yu SQ, Manasreh O, Salamo GJ:
Nanoscale Res. Lett.. 2010, 5:1079.
COI number [1:CAS:528:DC%2BC3cXmvFeksr8%3D]; Bibcode number [2010NRL.....5.1079L]
PubMed Abstract  Publisher Full Text  PubMed Central Full Text 
Shtrikman H, PopovitzBiro R, Kretinin A, Houben L, Heiblum M, Bukala M, Galicka M, Buczko R, Kacman P:
Nano Lett.. 2009, 9:1506.
COI number [1:CAS:528:DC%2BD1MXis1ars7s%3D]; Bibcode number [2009NanoL...9.1506S]
PubMed Abstract  Publisher Full Text 
Dubrovskii VG, Sibirev NV, Harmand JC, Glas F:
Phys. Rev. B. 2008, 78:235301.
Bibcode number [2008PhRvB..78w5301D]
Publisher Full Text 
Phys. Stat. Sol. B. 2010, 247:254.
COI number [1:CAS:528:DC%2BC3cXhtVOntL0%3D]; Bibcode number [2010PSSBR.247..254G]
Publisher Full Text 
Dubrovskii VG, Sibirev NV, Timofeeva MA:
Semiconductors. 2009, 43:1226.
COI number [1:CAS:528:DC%2BD1MXhtFaqtr7F]; Bibcode number [2009Semic..43.1226D]
Publisher Full Text 
Massalski TB (ed.) (Ed): Binary alloy phase diagrams. 1st edition. American Society for Metals, Metals Park, OH; 1986:260.
vol. 1

Harmand JC, Tchernycheva M, Patriarche G, Travers L, Glas F, Cirlin G:
J. Cryst. Growth. 2007, 301–302:853. Publisher Full Text

Tchernycheva M, Travers L, Patriarche G, Glas F, Harmand JC, Cirlin GE, Dubrovskii VG:
J. Appl. Phys.. 2007, 102:094313.
Bibcode number [2007JAP...102i4313T]
Publisher Full Text 
Glas F, Harmand JC, Patriarche G:
Phys. Rev. Lett.. 2010, 104:135501.
Bibcode number [2010PhRvL.104m5501G]
PubMed Abstract  Publisher Full Text