van Sebille et al. Nanoscale Research Letters (2016) 11:355 DOI 10.1186/s11671-016-1567-6
NANO EXPRESS
Open Access
Optimizing Silicon Oxide Embedded Silicon Nanocrystal Inter-particle Distances Martijn van Sebille1* , Jort Allebrandi1 , Jim Quik1 , René A. C. M. M. van Swaaij1 , Frans D. Tichelaar2 and Miro Zeman1 Abstract We demonstrate an analytical method to optimize the stoichiometry and thickness of multilayer silicon oxide films in order to achieve the highest density of non-touching and closely spaced silicon nanocrystals after annealing. The probability of a nanocrystal nearest-neighbor distance within a limited range is calculated using the stoichiometry of the as-deposited film and the crystallinity of the annealed film as input parameters. Multiplying this probability with the nanocrystal density results in the density of non-touching and closely spaced silicon nanocrystals. This method can be used to estimate the best as-deposited stoichiometry in order to achieve optimal nanocrystal density and spacing after a subsequent annealing step. Keywords: Silicon nanocrystal, Silicon oxide, Inter-particle distance, Stoichiometry, Spacing
Background Silicon nanocrystals embedded in a high band gap silicon alloy are interesting candidates for top cells of multijunction solar cells of which the band gap can be tuned by the nanocrystal size [1]. The ability to tune the material’s band gap allows us to minimize thermalization losses and thereby increase the solar cell efficiency. The mean nanocrystal size and size distribution are crucial parameters in determining the optical properties of the material [2, 3] and electronic transport properties in photovoltaic devices [4]. Embedded silicon nanocrystals can be made by annealing silicon-rich silicon alloy films, and this is typically performed using a tube furnace or rapid thermal annealing furnace [5]. Upon annealing at temperatures between 600 and 900 °C, phase separation of the excess silicon occurs, creating amorphous silicon nanoparticles surrounded by an amorphous silicon oxide matrix. Annealing at temperatures in excess of 900 °C leads to crystallization of these amorphous nanoparticles [6]. Using films containing alternating layers of stoichiometric and silicon-rich silicon alloys allows for the control over the nanocrystal size, limited by the silicon-rich layer *Correspondence: [email protected] 1 Photovoltaic Materials and Devices, Delft University of Technology, Mekelweg 4, 2628 CD Delft, The Netherlands Full list of author information is available at the end of the article
thickness [7, 8]. Various charge transport mechanisms for embedded silicon nanocrystal have been suggested, including direct tunneling [9], trap-assisted tunneling [10], and hopping [11]. No clear consensus exists concerning the exact mechanisms, especially concerning the role of defects in the matrix and at the nanocrystal interface [10–14]. Nonetheless, the total charge transport is expected to be highly dependent on the nanocrystal spacing and the choice of dielectric material [15]. For SiO2 films, inter-particle spacing up to 2 nm is acceptable, which provides a minimum mobility of 10−1 cm2 V−1 s, as calculated by Green et al. [1]. The nanocrystal density in the silicon-rich layers can be controlled by tuning the composition of these layers during deposition. A low silicon content leads to relatively few isolated nanocrystals, and increasing the excess silicon content will eventually lead to clustering of nanocrystals, shown schematically in Fig. 1. When the nanocrystal density is too low, the probability of nearest-neighbor nanocrystal within 2 nm is too low. In contrast, when the excess silicon content is too high, nanocrystals are so closely spaced that they start clustering, which reduces the quantum confinement in these crystals. This means there is an optimal composition to achieve a limited nanocrystal spacing, while limiting clustering. In this letter, we demonstrate an analytical method to optimize the composition and thickness of multilayer
© 2016 The Author(s). Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
van Sebille et al. Nanoscale Research Letters (2016) 11:355
a
Page 2 of 7
b
c
Fig. 1 Nanocrystal spacing and clustering. Nanocrystals formed in silicon-rich layers with relatively low (a), medium (b), and high (c) excess silicon, separated by stoichiometric buffer layers
silicon oxide films in order to achieve the highest density of non-touching and closely spaced silicon nanocrystals after annealing.
Methods We deposited approximately 100-nm-thick a-SiOx :H films on Spectrosil 2000 quartz substrates in a radio frequency plasma-enhanced chemical vapor deposition (PECVD) reactor, operating at 13.56 MHz. The following deposition parameters were used: a power density of 2.1 × 10−2 W cm−2 , a deposition pressure of 1.4 mbar, and a substrate temperature of 95 °C. The film composition was varied by changing the SiH4 over CO2 flow ratio from 0.07 to 0.37. A H2 flow rate of 200 sccm was used for all depositions. n- and p-type films were fabricated by including PH3 and B2 H6 flows, respectively. The dopant over SiH4 flow rate ratio was 2.0 × 10−3 . The atomic compositions of the silicon-rich and buffer layers were determined using x-ray photoelectron spectroscopy measurements, using a Thermo Scientific K-Alpha setup. The film surface was etched with an ion gun prior to measurements to remove surface contamination. Annealing was carried out using a Tempress horizontal tube stack or a Solaris 100 RTA furnace for 1 h and 3 min, respectively. All samples were annealed at 1000 °C, at atmospheric pressure and in pure nitrogen gas. The composition of the buffer layer used in these experiments is SiO1.3 . Measurements show that this stoichiometry is sufficiently high to prevent crystallization for the annealing conditions used (not shown here). Raman spectra were measured to determine the crystallinity, using a Renishaw inVia setup in backscattering geometry, with a 25-mW Ar laser as excitation source with a wavelength of 514 nm and focused on a spot of approximately 1 μm. The crystallinity XC is the ratio of the Si–Si bonds in crystalline phase over the Si–Si bonds in amorphous and crystalline phase [16] and is calculated as follows: XC =
ITO,c-Si , σ ITO,a-Si + ITO,c-Si
(1)
where ITO,c-Si and ITO,a-Si are the integrated TO phonon modes of crystalline and amorphous silicon, respectively.
σ is a factor to correct for the difference in scattering cross section between these modes and is set to 0.8 [16]. Imaging the silicon crystals in the amorphous silicon layer was done using a FEI Tecnai F20ST/STEM transmission electron microscope (TEM) operated at 200 kV. Thin samples for TEM were prepared in cross section following a standard procedure after gluing the two samples together face to face: a 500-μm-thick lamella was cut out using a diamond saw, subsequently thinned to approximately 15-μm thickness by mechanical polishing, glued on a copper support ring, and argon ion-milled to electron transparency. The silicon nanocrystals were marked using the freehand selection tool in ImageJ [17]. The surface area was then determined, and an effective diameter was recorded.
Results and Discussion Figure 2 illustrates nanocrystals with radius r in a multilayer structure, including their parameters needed to determine the inter-particle distance d.
Fig. 2 Nanocrystals in a multilayer structure shown schematically. Nanocrystals in a multilayer structure shown schematically, including the nanocrystal radius r, buffer layer thickness t, and inter-particle distance d. The enclosing box around a nanocrystal is shown for the right-hand nanocrystal
van Sebille et al. Nanoscale Research Letters (2016) 11:355
Page 3 of 7
We assume that the mean nanocrystal diameter equals the silicon-rich layer thickness. In order to validate this assumption, a multilayer sample with silicon-rich and buffer layer thicknesses of 3 and 1 nm, respectively, has been measured with high-resolution TEM, shown in Fig. 3a. The histogram of the obtained nanocrystal diameters is shown in Fig. 3b. The mean nanocrystal diameter obtained from TEM is 2.4 nm. Figure 4 shows the mean nanocrystal diameter as a function of its silicon-rich layer thickness of this sample, as well as data obtained by Gutsch et al. [18]. For thin silicon-rich layer thicknesses, the deviation between the mean nanocrystal diameter and the sample’s silicon-rich layer thickness is reasonable. For thicker layers, the deviation increases. However, we should note that
a Fig. 4 Nanocrystal diameter in a multilayer sample as a function of the silicon-rich layer thickness. The mean nanocrystal diameter D¯ (solid symbols) and the mean equivalent diameter D¯ equiv (open symbols) for samples with varying silicon-rich layer thicknesses. The black data points are obtained from Gutsch et al. [18]. The dashed line represents the equality between the nanocrystal diameter and the silicon-rich layer thickness
b
instead of being interested in nanocrystal diameters, we are interested in the volume these nanocrystals occupy, since that allows us to predict the nanocrystal density and their inter-particle distance for varying stoichiometries and crystallinities. The size distribution of such multilayer samples is log-normally shaped, as was observed by Gutsch et al. [18] and can be seen from Fig. 3b. Since the volume of the nanoparticles depends on the third power of their radius, the larger nanocrystals have a greater contribution to the mean volume V¯ . This is given by 4 3 3 πr , (2) V¯ = nNC where r and nNC are the nanocrystal radius and the number of nanocrystals obtained from TEM, respectively. The ¯ equiv can be equivalent diameter of the mean nanocrystal D expressed by V¯ ¯ equiv = 2 3 . (3) D 4 3π
Fig. 3 TEM image of a multilayer sample and its nanocrystal diameter histogram. a Cross-sectional high-resolution TEM image of an annealed multilayer sample with silicon-rich and buffer layer thicknesses of 3 and 1 nm, respectively. b The histogram of the sample’s nanocrystal diameters. Approximately 250 nanocrystals were measured. The histogram is fitted with a log-normal probability density function with μ = 0.83 nm and σ = 0.27 nm
Combining Eqs. (2) and (3) results in ¯ equiv = 2 3 r3 , D
(4)
where r3 represents the mean value of r3 . The equivalent diameter of the sample shown in Fig. 3 is 2.6 nm and is shown in Fig. 4 along with the equivalent diameters of
van Sebille et al. Nanoscale Research Letters (2016) 11:355
Page 4 of 7
the data obtained by Gutsch et al. [18]. Because of the asymmetrical, log-normally shaped nanocrystal size distributions, all equivalent diameters are greater than their corresponding mean diameters. In general, the equivalent diameters are very close to the assumed equality between the nanocrystal diameter and the silicon-rich layer thickness. This result implies that our assumption is reasonable, at least up to silicon-rich layer thicknesses up to 4.5 nm. Fortunately, this range is most interesting for photovoltaic purposes because of their increased confinement. Note that we do not include a core/shell structure in this approach. An amorphous sub-oxide shell is likely to form around silicon nanocrystals [7, 19]. Iacona et al. measured a shell to be approximately 1 nm thick [20]. This thickness corresponds with theoretical calculations and experimental measurements by Daldosso et al. [21]. However, Queeney et al. determined the shell thickness to be less than 6 Å [22]. Zimina et al. measured shell thicknesses of 2 to 5 Å and suggested a dependence on nanocrystal size [23]. These disagreeing results complicate the incorporation of a core/shell structure in our method. Since the aim of this publication is to provide a simple method to predict the inter-nanocrystal distance, we do not include the core/shell structure. In order to optimize the density of non-touching and closely spaced silicon nanocrystals, the nanocrystal density in the silicon-rich layers should be determined first. This depends on the excess silicon in these layers and can be calculated from its composition as follows: x x Si. (5) SiOx → SiO2 + 1 − 2 2 The excess silicon can be in amorphous or crystalline phase, so the total atomic density in the layer ρlayer is given by ρlayer = c-Si ρc-Si + a-Si ρa-Si + SiO2 ρSiO2 ,
n2D =
ρc-Si,layer 2r. NNC
(11)
Figure 5a shows the 2D nanocrystal density for a sample with silicon-rich and buffer layer thicknesses of 3 and 1 nm, respectively, using ρc-Si , ρa-Si , and ρSiO2 5.0 × 1028 , 5.0 × 1028 , and 2.2 × 1028 m−3 , respectively [24–27].
a
b
c
(6)
where ρc-Si and ρa-Si are the atomic densities of c-Si and a-Si, respectively; ρSiO2 is the molecular density of SiO2 ; and c-Si , a-Si , and SiO2 are their respective atomic and molecular percentages. Using the definition of crystallinity and Eq. (5), the atomic percentages of c-Si, a-Si, and SiO2 can be written as x · 100% (7) c-Si = XC 1 − 2 x a-Si = (1 − XC ) 1 − · 100% (8) 2 x (9) SiO2 = · 100% 2 The atomic density of c-Si in the layer can then be determined by ρc-Si,layer = c-Si ρlayer .
The number of atoms in a nanocrystal is NNC = 4 3 3 πr ρc-Si , which can be used to calculate the 2D nanocrystal density n2D in a silicon-rich layer with thickness 2r
(10)
Fig. 5 Nanocrystal density, NN probability, and NN density. The 2D nanocrystal (NC) density (a), the probability of finding a nearest neighbor (NN) within 2 nm (b), and the density of nanocrystals with a NN within 2 nm (c) as a function of the silicon-rich layer composition and crystallinity for a sample with silicon-rich and buffer layer thicknesses of 3 and 1 nm, respectively. The black diamonds represent tube furnace annealed intrinsic samples. The red, magenta, and blue squares show intrinsic, p-type and n-type samples annealed using RTA
van Sebille et al. Nanoscale Research Letters (2016) 11:355
Page 5 of 7
The 2D nanocrystal density is highest for silicon-rich layers with high crystallinity and low stoichiometry. However, in that case, the nanocrystal density can be so high that crystals cluster together as illustrated in Fig. 1c. To find the fraction of nanocrystals that are properly spaced, we use the probability density function F of finding a nearest neighbor at distance d for a nanocrystal in a multilayer sample. For randomly distributed point particles, this is given by [28]
F (d) = 4j + 2 n2D πd exp − 2j + 1 n2D πd2 ×
n2D πt 2 j j + 1 2j + 1 , (12) exp 3 where t is the buffer layer thickness and j = r/t. Note that this probability density function describes centerto-center distances. Integrating this equation from 2r to 2r+2 nm provides the edge-to-edge probability of finding a nearest neighbor within 2 nm for a multilayer structure with silicon-rich layer thickness of 3 nm. This is shown in Fig. 5b and confirms that the optimal probability of finding the nearest neighbor between 0 to 2 nm is very low in the range where the 2D nanocrystal density is highest, caused by clustering of the nanocrystals. We can find an optimum stoichiometry for a given crystallinity using the result shown in Fig. 5b, but aside from proper spacing of nanocrystals, we are also interested in a high nanocrystal density. Integrating Eq. (12) and multiplying with the 2D nanocrystal density provides the non-touching nanocrystal density with a nearest neighbor within d nNN (d) = n2D
2r+d
2r
F (d) dd.
(13)
The nanocrystal density with a nearest neighbor within 2 nm for a sample with silicon-rich layer thickness of 3 nm is shown in Fig. 5c. The highest non-touching nanocrystal density for this structure can be achieved by tuning the silicon-rich layer’s composition from pure Si to approximately SiO1.5 , with crystallinity values of 0.15 to 1, respectively. A too high crystallinity for layers with relatively low stoichiometry will lead to clustering, while too low crystallinity in layers with high stoichiometry will result in separated but too isolated nanocrystals. We deposited multilayer films with silicon-rich layer and buffer layer thicknesses of 3 and 1 nm, respectively, and varied the silicon-rich layer’s compositions and doping. The films were annealed in a tube furnace or rapid thermal annealing furnace. The crystallinity of these films after annealing are shown in Fig. 5.1 The crystallinity
in our films does not exceed 0.4. We expect that this is caused by incomplete phase separation prior to crystallization, as observed before for silicon oxide films deposited using PECVD [29]. This means that only part of all excess silicon (see Eq. (5)) clusters into silicon nanoparticles, which can subsequently crystallize. The rest of the excess silicon remains in the surrounding matrix, which will not be SiO2 , but has a lower stoichiometry. This in turn will lower its energy barrier, increasing the tunneling probability and possibly enlarge the inter-particle distance at which the mobility remains sufficiently high [30]. However, for simplicity, we will keep an inter-particle distance range from 0 to 2 nm. Assuming a maximum achievable crystallinity of 0.4 for PECVD films, an optimal stoichiometry to achieve the highest density of non-touching, closely spaced nanocrystals can be found. This optimal stoichiometry is SiO0.84 . In contrast, films deposited using magnetron sputtering are reported to lead to complete phase separation [29]. Assuming all silicon clusters crystallize upon annealing, this will lead to a crystallinity equal to unity. In reality, the sub-oxide shell around the nanocrystal core will limit complete crystallization [31], but for simplicity, we assume a crystallinity equal to unity. In that case, the optimal stoichiometry of the siliconrich layers is approximately SiO1.4 for this structure. Note that in both cases, the 2D nanocrystal density with nearest neighbor within 2 nm is 1.3 × 1012 cm−2 . This value corresponds well with results obtained experimentally by Laube et al. for single 4.5-nm-thick layers [32] and is slightly lower than experimental results obtained by Gutsch et al. for single 3.5-nm-thick layers [18]. However, we should note that these reported values are the total nanocrystal density, while we estimated the isolated nanocrystal density with nearest neighbor within 2 nm. This excludes clustered nanocrystals and too isolated nanocrystals, which inevitably leaves out a portion of the total nanocrystal density. Furthermore, the thickness of the silicon-rich layer affects the estimated nanocrystal density, with lower values for thicker layers. The optimal stoichiometry decreases for increasing silicon-rich layer thicknesses, as shown in Fig. 6. Note that the model’s accuracy decreases for greater silicon-rich layers thicknesses, as shown in Fig. 4. Nonetheless, from a theoretical perspective, the observed trend for thick silicon-rich layers is still interesting. This is caused by the differences in volume between a nanocrystal and its enclosing box (see Fig. 2). The volume of a nanocrystal is 43 πr3 , and its enclosing box is approximately 2r(2r + d)2 . The volume ratio of the enclosing box over the nanocrystal decreases with increasing siliconrich layer thickness, explaining the trend observed in Fig. 6 for relatively large silicon-rich layer thicknesses. However, for very small silicon-rich layer thicknesses, the optimal composition does not vary. To explain this, we
van Sebille et al. Nanoscale Research Letters (2016) 11:355
Page 6 of 7
thicknesses, the probability density functions broaden. The optimal composition depends on the probability of finding a nearest neighbor within a limited range (2 nm for SiOx ). For silicon-rich layer thicknesses up to approximately 2 nm, the probability of finding a nearest neighbor beyond 2 nm is negligible, as can be observed in Fig. 7. For these thicknesses, there is no reason to increase the nanocrystal density, since more closely packed nanocrystals will not increase the probability of finding a nearest neighbor within 2 nm. On the contrary, an increase in nanocrystal density will increase the probability of clustering.
Conclusions
Fig. 6 Optimal stoichiometry. The optimal stoichiometry for films with varying silicon-rich layer thicknesses for a sample with crystallinity of 1 (solid line) and 0.4 (dashed line). The buffer layer thickness is kept constant at 1 nm
look closer into their probability density functions, shown in Fig. 7. The peak at 1 nm is caused by the availability of nearest neighbors in the neighboring silicon-rich layers. The probability density function broadens for larger siliconrich layer thicknesses because its standard deviation is √ related to the 2D nanocrystal density by σ2D ∝ 1/ n2D [28]. Since n2D decreases for greater silicon-rich layer
We demonstrated an analytical method to optimize the composition of a silicon-rich layer for different crystallinities and thicknesses in order to achieve the highest density of non-touching and closely spaced silicon nanocrystals after annealing. The optimal stoichiometry depends on the crystallinity decreases for increasing silicon-rich layer thicknesses. However, for very small silicon-rich layer thicknesses, the optimal composition does not vary. This method can be used to find the best as-deposited composition in order to achieve optimal nanocrystal density and spacing after a subsequent annealing step.
Endnote 1 Note that the buffer layer composition is not stoichiometric and therefore will contain some excess silicon as well. This excess silicon will increase the amorphous Si– Si bond density, resulting in an underestimation of the crystallinity. However, since the buffer layer thickness is only 1 nm, compared to 3 nm for the silicon-rich layer, we expect this effect to be limited and assume it can be neglected. Acknowledgements We acknowledge financial support for this research from ADEM, A green Deal in Energy Materials of the Ministry of Economic Affairs of The Netherlands (http://www.adem-innovationlab.nl). We also thank ing. M.W.G.M. Verhoeven from the Eindhoven University of Technology for his assistance with the XPS measurements. Authors’ Contributions MS and JA initiated the research and performed the calculations. JA and JQ performed the sample fabrication and characterization. MS, RS, and MZ directed the experiments and the drafting of the paper. FT performed and analyzed the TEM measurements. MS wrote the manuscript. All authors read and approved the final manuscript.
Fig. 7 Nearest-neighbor distance probability density functions. The nearest-neighbor distance probability density functions for samples with silicon-rich layer thicknesses of 1.5, 2, and 3 nm and buffer layer thickness of 1 nm, calculated using their respective optimal compositions. The gray area depicts the range of desired nearest-neighbor distances. Shorter distances imply clustering, and greater distances lead to an insufficient tunneling probability
Competing Interests The authors declare that they have no competing interests. Author details 1 Photovoltaic Materials and Devices, Delft University of Technology, Mekelweg 4, 2628 CD Delft, The Netherlands. 2 Kavli Institute of Nanoscience, Quantum Nanoscience, Delft University of Technology, Lorentzweg 1, 2628 CJ Delft, The Netherlands.
van Sebille et al. Nanoscale Research Letters (2016) 11:355
Received: 28 May 2016 Accepted: 26 July 2016
References 1. Green MA, Cho EC, Cho Y, Huang Y, Pink E, Trupke T, Lin A, Fangsuwannarak T, Puzzer T, Conibeer G, Corkish R (2005) All-silicon tandem cells based on artificial semiconductor synthesised using silicon quantum dots in a dielectric matrix. In: Proceedings of the 20th European Photovoltaic Solar Energy Conference and Exhibition. pp 3–7 2. Heitmann J, Müller F, Zacharias M, Gösele U (2005) Silicon nanocrystals: size matters. Adv Mater 17(7):795–803. doi:10.1002/adma.200401126 3. Künle M, Kaltenbach T, Löper P, Hartel A, Janz S, Eibl O, Nickel KG (2010) Si-rich a-SiC:H thin films: structural and optical transformations during thermal annealing. Thin Solid Films 519(1):151–157. doi:10.1016/j.tsf.2010.07.085 4. Puthen-Veettil B, Patterson R, König D, Conibeer G, Green MA (2014) The impact of disorder on charge transport in three dimensional quantum dot resonant tunneling structures. J Appl Phys 116(16):163707. doi:10.1063/1.4899207 5. Schnabel M, Weiss C, Löper P, Wilshaw PR, Janz S (2015) Self-assembled silicon nanocrystal arrays for photovoltaics. Phys Stat Solidi A 212(8):1649–1661. doi:10.1002/pssa.201431764 6. Yi LX, Heitmann J, Scholz R, Zacharias M (2002) Si rings, Si clusters, and Si nanocrystals—different states of ultrathin SiOx layers. Appl Phys Lett 81(22):4248. doi:10.1063/1.1525051 7. López-Vidrier J, Hernández S, Hartel AM, Hiller D, Gutsch S, Löper P, López-Conesa L, Estradé S, Peiró F, Zacharias M, Garrido B (2011) Structural and optical characterization of size controlled silicon nanocrystals in SiO2 /SiOx Ny multilayers. Energy Procedia 10:43–48. doi:10.1016/j.egypro.2011.10.150 8. Conibeer G, Perez-Wurfl I, Hao X, Di D, Lin D (2012) Si solid-state quantum dot-based materials for tandem solar cells. Nanoscale Res Lett 7(1):193. doi:10.1186/1556-276x-7-193 9. Garrido B, Hernández S, Berencén Y, López-Vidrier J, Ramírez J, Blázquez O, Mundet B (2015) Electrical transport in Si-based nanostructured superlattices. In: Light Energy Harvesting with Group IV Nanostructures. Pan Stanford. pp 255–316. doi:10.1201/b18090-9. http://dx.doi.org/10. 1201/b18090-9 10. Mavilla NR, Chavan V, Solanki CS, Vasi J (2016) Study of temperature-dependent charge conduction in silicon-nanocrystal/SiO2 multilayers. Thin Solid Films 612:41–48. doi:10.1016/j.tsf.2016.05.029 11. Balberg I (2011) Electrical transport mechanisms in three dimensional ensembles of silicon quantum dots. J Appl Phys 110(6):061301. doi:10.1063/1.3637636 12. Jiang CW, Green MA (2006) Silicon quantum dot superlattices: modeling of energy bands, densities of states, and mobilities for silicon tandem solar cell applications. J Appl Phys 99(11):114902. doi:10.1063/1.2203394 13. Stegner AR, Pereira RN, Klein K, Lechner R, Dietmueller R, Brandt MS, Stutzmann M, Wiggers H (2008) Electronic transport in phosphorus-doped silicon nanocrystal networks. Phys Rev Lett 100(2):026803. doi:10.1103/physrevlett.100.026803 14. Lechner R, Stegner AR, Pereira RN, Dietmueller R, Brandt MS, Ebbers A, Trocha M, Wiggers H, Stutzmann M (2008) Electronic properties of doped silicon nanocrystal films. J Appl Phys 104(5):053701. doi:10.1063/1.2973399 15. Conibeer G, Green MA, König D, Perez-Wurfl I, Huang S, Hao X, Di D, Shi L, Shrestha S, Puthen-Veetil B, So Y, Zhang B, Wan Z (2010) Silicon quantum dot based solar cells: addressing the issues of doping, voltage and current transport. Prog Photovolt Res Appl 19(7):813–824. doi:10.1002/pip.1045 16. Smit C, van Swaaij RACMM, Donker H, Petit AMHN, Kessels WMM, van de Sanden MCM (2003) Determining the material structure of microcrystalline silicon from Raman spectra. J Appl Phys 94(5):3582. doi:10.1063/1.1596364 17. Schneider CA, Rasband WS, Eliceiri KW (2012) NIH image to ImageJ: 25 years of image analysis. Nat Methods 9(7):671–675. doi:10.1038/nmeth.2089 18. Gutsch S, Hiller D, Laube J, Zacharias M, Kübel C (2015) Observing the morphology of single-layered embedded silicon nanocrystals by using temperature-stable TEM membranes. Beilstein J Nanotechnol 6:964–970. doi:10.3762/bjnano.6.99
Page 7 of 7
19. Hiller D (2015) Silicon nanoclusters embedded in dielectric matrices: nucleation, growth, crystallization, and defects. In: Valenta J, Mirabella S (eds). Silicon nanoclusters embedded in dielectric matrices: nucleation, growth, crystallization, and defects. Pan Stanford Publishing. pp 99–144. doi:10.1201/b18090-5. http://dx.doi.org/10.1201/b18090-5 20. Iacona F, Bongiorno C, Spinella C, Boninelli S, Priolo F (2004) Formation and evolution of luminescent Si nanoclusters produced by thermal annealing of SiOx films. J Appl Phys 95(7):3723–3732. doi:10.1063/1.1664026 21. Daldosso N, Luppi M, Dalba G, Pavesi L, Rocca F, Priolo F, Franzò G, Iacona F, Degoli E, Magri R, Ossicini S (2003) Experimental and theoretical joint study on the electronic and structural properties of silicon nanocrystals embedded in SiO2 : active role of the interface region. In: MRS Proc Vol. 770. p 6. http://journals.cambridge.org/action/displayAbstract? fromPage=online&aid=8028079&fileId=S1946427400137005 22. Queeney KT, Weldon MK, Chang JP, Chabal YJ, Gurevich AB, Sapjeta J, Opila RL (2000) Infrared spectroscopic analysis of the Si/SiO2 interface structure of thermally oxidized silicon. J Appl Phys 87(3):1322. doi:10.1063/1.372017 23. Zimina A, Eisebitt S, Eberhardt W, Heitmann J, Zacharias M (2006) Electronic structure and chemical environment of silicon nanoclusters embedded in a silicon dioxide matrix. Appl Phys Lett 88(16):163103. doi:10.1063/1.2193810 24. Zhang XG (2001) Silicon/electrolyte interface. Springer. doi:10.1007/0-306-47921-4_2. http://dx.doi.org/10.1007/0-306-47921-4_ 2 25. Kessels WMM, Sanden MCMVD, Severens RJ, Ijzendoorn LJV, Schram DC (1998) Hydrogen in a-Si:H deposited by an expanding thermal plasma: a temperature, growth rate and isotope study. In: MRS Proceedings. Cambridge University Press Vol. 507. pp 529–534. doi:10.1557/proc-507-529. http://dx.doi.org/10.1557/proc-507-529 26. Sze SM, Ng KK (2006) Appendix G: properties of Si and GaAs. Wiley-Blackwell. doi:10.1002/9780470068328.app7. http://dx.doi.org/10. 1002/9780470068328.app7 27. Zhang XG (2001) Anodic oxide. Springer. doi:10.1007/0-306-47921-4_3. http://dx.doi.org/10.1007/0-306-47921-4_3 28. Tournus F (2011) Random nanoparticle deposition: inter-particle distances in 2d, 3d, and multilayer samples. J Nanopart Res 13(10):5211–5223. doi:10.1007/s11051-011-0506-9 29. Franzò G, Miritello M, Boninelli S, Lo Savio R, Grimaldi MG, Priolo F, Iacona F, Nicotra G, Spinella C, Coffa S (2008) Microstructural evolution of SiOx films and its effect on the luminescence of Si nanoclusters. J Appl Phys 104(9):094306. doi:10.1063/1.3006735 30. Conibeer G, Green MA, Corkish R, Cho Y, Cho EC, Jiang CW, Fangsuwannarak T, Pink E, Huang Y, Puzzer T, Trupke T, Richards B, Shalav A, Lin K-l (2006) Silicon nanostructures for third generation photovoltaic solar cells. Thin Solid Films 511–512:654–662. doi:10.1016/j.tsf.2005.12.119 31. López-Vidrier J, Hernández S, Hiller D, Gutsch S, López-Conesa L, Estradé S, Peiró F, Zacharias M, Garrido B (2014) Annealing temperature and barrier thickness effect on the structural and optical properties of silicon nanocrystals/SiO2 superlattices. J Appl Phys 116(13):133505. doi:10.1063/1.4896878 32. Laube J, Gutsch S, Wang D, Kübel C, Zacharias M, Hiller D (2016) Two-dimensional percolation threshold in confined Si nanoparticle networks. Appl Phys Lett 108(4):043106. doi:10.1063/1.4940971
NANO EXPRESS
Open Access
Optimizing Silicon Oxide Embedded Silicon Nanocrystal Inter-particle Distances Martijn van Sebille1* , Jort Allebrandi1 , Jim Quik1 , René A. C. M. M. van Swaaij1 , Frans D. Tichelaar2 and Miro Zeman1 Abstract We demonstrate an analytical method to optimize the stoichiometry and thickness of multilayer silicon oxide films in order to achieve the highest density of non-touching and closely spaced silicon nanocrystals after annealing. The probability of a nanocrystal nearest-neighbor distance within a limited range is calculated using the stoichiometry of the as-deposited film and the crystallinity of the annealed film as input parameters. Multiplying this probability with the nanocrystal density results in the density of non-touching and closely spaced silicon nanocrystals. This method can be used to estimate the best as-deposited stoichiometry in order to achieve optimal nanocrystal density and spacing after a subsequent annealing step. Keywords: Silicon nanocrystal, Silicon oxide, Inter-particle distance, Stoichiometry, Spacing
Background Silicon nanocrystals embedded in a high band gap silicon alloy are interesting candidates for top cells of multijunction solar cells of which the band gap can be tuned by the nanocrystal size [1]. The ability to tune the material’s band gap allows us to minimize thermalization losses and thereby increase the solar cell efficiency. The mean nanocrystal size and size distribution are crucial parameters in determining the optical properties of the material [2, 3] and electronic transport properties in photovoltaic devices [4]. Embedded silicon nanocrystals can be made by annealing silicon-rich silicon alloy films, and this is typically performed using a tube furnace or rapid thermal annealing furnace [5]. Upon annealing at temperatures between 600 and 900 °C, phase separation of the excess silicon occurs, creating amorphous silicon nanoparticles surrounded by an amorphous silicon oxide matrix. Annealing at temperatures in excess of 900 °C leads to crystallization of these amorphous nanoparticles [6]. Using films containing alternating layers of stoichiometric and silicon-rich silicon alloys allows for the control over the nanocrystal size, limited by the silicon-rich layer *Correspondence: [email protected] 1 Photovoltaic Materials and Devices, Delft University of Technology, Mekelweg 4, 2628 CD Delft, The Netherlands Full list of author information is available at the end of the article
thickness [7, 8]. Various charge transport mechanisms for embedded silicon nanocrystal have been suggested, including direct tunneling [9], trap-assisted tunneling [10], and hopping [11]. No clear consensus exists concerning the exact mechanisms, especially concerning the role of defects in the matrix and at the nanocrystal interface [10–14]. Nonetheless, the total charge transport is expected to be highly dependent on the nanocrystal spacing and the choice of dielectric material [15]. For SiO2 films, inter-particle spacing up to 2 nm is acceptable, which provides a minimum mobility of 10−1 cm2 V−1 s, as calculated by Green et al. [1]. The nanocrystal density in the silicon-rich layers can be controlled by tuning the composition of these layers during deposition. A low silicon content leads to relatively few isolated nanocrystals, and increasing the excess silicon content will eventually lead to clustering of nanocrystals, shown schematically in Fig. 1. When the nanocrystal density is too low, the probability of nearest-neighbor nanocrystal within 2 nm is too low. In contrast, when the excess silicon content is too high, nanocrystals are so closely spaced that they start clustering, which reduces the quantum confinement in these crystals. This means there is an optimal composition to achieve a limited nanocrystal spacing, while limiting clustering. In this letter, we demonstrate an analytical method to optimize the composition and thickness of multilayer
© 2016 The Author(s). Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
van Sebille et al. Nanoscale Research Letters (2016) 11:355
a
Page 2 of 7
b
c
Fig. 1 Nanocrystal spacing and clustering. Nanocrystals formed in silicon-rich layers with relatively low (a), medium (b), and high (c) excess silicon, separated by stoichiometric buffer layers
silicon oxide films in order to achieve the highest density of non-touching and closely spaced silicon nanocrystals after annealing.
Methods We deposited approximately 100-nm-thick a-SiOx :H films on Spectrosil 2000 quartz substrates in a radio frequency plasma-enhanced chemical vapor deposition (PECVD) reactor, operating at 13.56 MHz. The following deposition parameters were used: a power density of 2.1 × 10−2 W cm−2 , a deposition pressure of 1.4 mbar, and a substrate temperature of 95 °C. The film composition was varied by changing the SiH4 over CO2 flow ratio from 0.07 to 0.37. A H2 flow rate of 200 sccm was used for all depositions. n- and p-type films were fabricated by including PH3 and B2 H6 flows, respectively. The dopant over SiH4 flow rate ratio was 2.0 × 10−3 . The atomic compositions of the silicon-rich and buffer layers were determined using x-ray photoelectron spectroscopy measurements, using a Thermo Scientific K-Alpha setup. The film surface was etched with an ion gun prior to measurements to remove surface contamination. Annealing was carried out using a Tempress horizontal tube stack or a Solaris 100 RTA furnace for 1 h and 3 min, respectively. All samples were annealed at 1000 °C, at atmospheric pressure and in pure nitrogen gas. The composition of the buffer layer used in these experiments is SiO1.3 . Measurements show that this stoichiometry is sufficiently high to prevent crystallization for the annealing conditions used (not shown here). Raman spectra were measured to determine the crystallinity, using a Renishaw inVia setup in backscattering geometry, with a 25-mW Ar laser as excitation source with a wavelength of 514 nm and focused on a spot of approximately 1 μm. The crystallinity XC is the ratio of the Si–Si bonds in crystalline phase over the Si–Si bonds in amorphous and crystalline phase [16] and is calculated as follows: XC =
ITO,c-Si , σ ITO,a-Si + ITO,c-Si
(1)
where ITO,c-Si and ITO,a-Si are the integrated TO phonon modes of crystalline and amorphous silicon, respectively.
σ is a factor to correct for the difference in scattering cross section between these modes and is set to 0.8 [16]. Imaging the silicon crystals in the amorphous silicon layer was done using a FEI Tecnai F20ST/STEM transmission electron microscope (TEM) operated at 200 kV. Thin samples for TEM were prepared in cross section following a standard procedure after gluing the two samples together face to face: a 500-μm-thick lamella was cut out using a diamond saw, subsequently thinned to approximately 15-μm thickness by mechanical polishing, glued on a copper support ring, and argon ion-milled to electron transparency. The silicon nanocrystals were marked using the freehand selection tool in ImageJ [17]. The surface area was then determined, and an effective diameter was recorded.
Results and Discussion Figure 2 illustrates nanocrystals with radius r in a multilayer structure, including their parameters needed to determine the inter-particle distance d.
Fig. 2 Nanocrystals in a multilayer structure shown schematically. Nanocrystals in a multilayer structure shown schematically, including the nanocrystal radius r, buffer layer thickness t, and inter-particle distance d. The enclosing box around a nanocrystal is shown for the right-hand nanocrystal
van Sebille et al. Nanoscale Research Letters (2016) 11:355
Page 3 of 7
We assume that the mean nanocrystal diameter equals the silicon-rich layer thickness. In order to validate this assumption, a multilayer sample with silicon-rich and buffer layer thicknesses of 3 and 1 nm, respectively, has been measured with high-resolution TEM, shown in Fig. 3a. The histogram of the obtained nanocrystal diameters is shown in Fig. 3b. The mean nanocrystal diameter obtained from TEM is 2.4 nm. Figure 4 shows the mean nanocrystal diameter as a function of its silicon-rich layer thickness of this sample, as well as data obtained by Gutsch et al. [18]. For thin silicon-rich layer thicknesses, the deviation between the mean nanocrystal diameter and the sample’s silicon-rich layer thickness is reasonable. For thicker layers, the deviation increases. However, we should note that
a Fig. 4 Nanocrystal diameter in a multilayer sample as a function of the silicon-rich layer thickness. The mean nanocrystal diameter D¯ (solid symbols) and the mean equivalent diameter D¯ equiv (open symbols) for samples with varying silicon-rich layer thicknesses. The black data points are obtained from Gutsch et al. [18]. The dashed line represents the equality between the nanocrystal diameter and the silicon-rich layer thickness
b
instead of being interested in nanocrystal diameters, we are interested in the volume these nanocrystals occupy, since that allows us to predict the nanocrystal density and their inter-particle distance for varying stoichiometries and crystallinities. The size distribution of such multilayer samples is log-normally shaped, as was observed by Gutsch et al. [18] and can be seen from Fig. 3b. Since the volume of the nanoparticles depends on the third power of their radius, the larger nanocrystals have a greater contribution to the mean volume V¯ . This is given by 4 3 3 πr , (2) V¯ = nNC where r and nNC are the nanocrystal radius and the number of nanocrystals obtained from TEM, respectively. The ¯ equiv can be equivalent diameter of the mean nanocrystal D expressed by V¯ ¯ equiv = 2 3 . (3) D 4 3π
Fig. 3 TEM image of a multilayer sample and its nanocrystal diameter histogram. a Cross-sectional high-resolution TEM image of an annealed multilayer sample with silicon-rich and buffer layer thicknesses of 3 and 1 nm, respectively. b The histogram of the sample’s nanocrystal diameters. Approximately 250 nanocrystals were measured. The histogram is fitted with a log-normal probability density function with μ = 0.83 nm and σ = 0.27 nm
Combining Eqs. (2) and (3) results in ¯ equiv = 2 3 r3 , D
(4)
where r3 represents the mean value of r3 . The equivalent diameter of the sample shown in Fig. 3 is 2.6 nm and is shown in Fig. 4 along with the equivalent diameters of
van Sebille et al. Nanoscale Research Letters (2016) 11:355
Page 4 of 7
the data obtained by Gutsch et al. [18]. Because of the asymmetrical, log-normally shaped nanocrystal size distributions, all equivalent diameters are greater than their corresponding mean diameters. In general, the equivalent diameters are very close to the assumed equality between the nanocrystal diameter and the silicon-rich layer thickness. This result implies that our assumption is reasonable, at least up to silicon-rich layer thicknesses up to 4.5 nm. Fortunately, this range is most interesting for photovoltaic purposes because of their increased confinement. Note that we do not include a core/shell structure in this approach. An amorphous sub-oxide shell is likely to form around silicon nanocrystals [7, 19]. Iacona et al. measured a shell to be approximately 1 nm thick [20]. This thickness corresponds with theoretical calculations and experimental measurements by Daldosso et al. [21]. However, Queeney et al. determined the shell thickness to be less than 6 Å [22]. Zimina et al. measured shell thicknesses of 2 to 5 Å and suggested a dependence on nanocrystal size [23]. These disagreeing results complicate the incorporation of a core/shell structure in our method. Since the aim of this publication is to provide a simple method to predict the inter-nanocrystal distance, we do not include the core/shell structure. In order to optimize the density of non-touching and closely spaced silicon nanocrystals, the nanocrystal density in the silicon-rich layers should be determined first. This depends on the excess silicon in these layers and can be calculated from its composition as follows: x x Si. (5) SiOx → SiO2 + 1 − 2 2 The excess silicon can be in amorphous or crystalline phase, so the total atomic density in the layer ρlayer is given by ρlayer = c-Si ρc-Si + a-Si ρa-Si + SiO2 ρSiO2 ,
n2D =
ρc-Si,layer 2r. NNC
(11)
Figure 5a shows the 2D nanocrystal density for a sample with silicon-rich and buffer layer thicknesses of 3 and 1 nm, respectively, using ρc-Si , ρa-Si , and ρSiO2 5.0 × 1028 , 5.0 × 1028 , and 2.2 × 1028 m−3 , respectively [24–27].
a
b
c
(6)
where ρc-Si and ρa-Si are the atomic densities of c-Si and a-Si, respectively; ρSiO2 is the molecular density of SiO2 ; and c-Si , a-Si , and SiO2 are their respective atomic and molecular percentages. Using the definition of crystallinity and Eq. (5), the atomic percentages of c-Si, a-Si, and SiO2 can be written as x · 100% (7) c-Si = XC 1 − 2 x a-Si = (1 − XC ) 1 − · 100% (8) 2 x (9) SiO2 = · 100% 2 The atomic density of c-Si in the layer can then be determined by ρc-Si,layer = c-Si ρlayer .
The number of atoms in a nanocrystal is NNC = 4 3 3 πr ρc-Si , which can be used to calculate the 2D nanocrystal density n2D in a silicon-rich layer with thickness 2r
(10)
Fig. 5 Nanocrystal density, NN probability, and NN density. The 2D nanocrystal (NC) density (a), the probability of finding a nearest neighbor (NN) within 2 nm (b), and the density of nanocrystals with a NN within 2 nm (c) as a function of the silicon-rich layer composition and crystallinity for a sample with silicon-rich and buffer layer thicknesses of 3 and 1 nm, respectively. The black diamonds represent tube furnace annealed intrinsic samples. The red, magenta, and blue squares show intrinsic, p-type and n-type samples annealed using RTA
van Sebille et al. Nanoscale Research Letters (2016) 11:355
Page 5 of 7
The 2D nanocrystal density is highest for silicon-rich layers with high crystallinity and low stoichiometry. However, in that case, the nanocrystal density can be so high that crystals cluster together as illustrated in Fig. 1c. To find the fraction of nanocrystals that are properly spaced, we use the probability density function F of finding a nearest neighbor at distance d for a nanocrystal in a multilayer sample. For randomly distributed point particles, this is given by [28]
F (d) = 4j + 2 n2D πd exp − 2j + 1 n2D πd2 ×
n2D πt 2 j j + 1 2j + 1 , (12) exp 3 where t is the buffer layer thickness and j = r/t. Note that this probability density function describes centerto-center distances. Integrating this equation from 2r to 2r+2 nm provides the edge-to-edge probability of finding a nearest neighbor within 2 nm for a multilayer structure with silicon-rich layer thickness of 3 nm. This is shown in Fig. 5b and confirms that the optimal probability of finding the nearest neighbor between 0 to 2 nm is very low in the range where the 2D nanocrystal density is highest, caused by clustering of the nanocrystals. We can find an optimum stoichiometry for a given crystallinity using the result shown in Fig. 5b, but aside from proper spacing of nanocrystals, we are also interested in a high nanocrystal density. Integrating Eq. (12) and multiplying with the 2D nanocrystal density provides the non-touching nanocrystal density with a nearest neighbor within d nNN (d) = n2D
2r+d
2r
F (d) dd.
(13)
The nanocrystal density with a nearest neighbor within 2 nm for a sample with silicon-rich layer thickness of 3 nm is shown in Fig. 5c. The highest non-touching nanocrystal density for this structure can be achieved by tuning the silicon-rich layer’s composition from pure Si to approximately SiO1.5 , with crystallinity values of 0.15 to 1, respectively. A too high crystallinity for layers with relatively low stoichiometry will lead to clustering, while too low crystallinity in layers with high stoichiometry will result in separated but too isolated nanocrystals. We deposited multilayer films with silicon-rich layer and buffer layer thicknesses of 3 and 1 nm, respectively, and varied the silicon-rich layer’s compositions and doping. The films were annealed in a tube furnace or rapid thermal annealing furnace. The crystallinity of these films after annealing are shown in Fig. 5.1 The crystallinity
in our films does not exceed 0.4. We expect that this is caused by incomplete phase separation prior to crystallization, as observed before for silicon oxide films deposited using PECVD [29]. This means that only part of all excess silicon (see Eq. (5)) clusters into silicon nanoparticles, which can subsequently crystallize. The rest of the excess silicon remains in the surrounding matrix, which will not be SiO2 , but has a lower stoichiometry. This in turn will lower its energy barrier, increasing the tunneling probability and possibly enlarge the inter-particle distance at which the mobility remains sufficiently high [30]. However, for simplicity, we will keep an inter-particle distance range from 0 to 2 nm. Assuming a maximum achievable crystallinity of 0.4 for PECVD films, an optimal stoichiometry to achieve the highest density of non-touching, closely spaced nanocrystals can be found. This optimal stoichiometry is SiO0.84 . In contrast, films deposited using magnetron sputtering are reported to lead to complete phase separation [29]. Assuming all silicon clusters crystallize upon annealing, this will lead to a crystallinity equal to unity. In reality, the sub-oxide shell around the nanocrystal core will limit complete crystallization [31], but for simplicity, we assume a crystallinity equal to unity. In that case, the optimal stoichiometry of the siliconrich layers is approximately SiO1.4 for this structure. Note that in both cases, the 2D nanocrystal density with nearest neighbor within 2 nm is 1.3 × 1012 cm−2 . This value corresponds well with results obtained experimentally by Laube et al. for single 4.5-nm-thick layers [32] and is slightly lower than experimental results obtained by Gutsch et al. for single 3.5-nm-thick layers [18]. However, we should note that these reported values are the total nanocrystal density, while we estimated the isolated nanocrystal density with nearest neighbor within 2 nm. This excludes clustered nanocrystals and too isolated nanocrystals, which inevitably leaves out a portion of the total nanocrystal density. Furthermore, the thickness of the silicon-rich layer affects the estimated nanocrystal density, with lower values for thicker layers. The optimal stoichiometry decreases for increasing silicon-rich layer thicknesses, as shown in Fig. 6. Note that the model’s accuracy decreases for greater silicon-rich layers thicknesses, as shown in Fig. 4. Nonetheless, from a theoretical perspective, the observed trend for thick silicon-rich layers is still interesting. This is caused by the differences in volume between a nanocrystal and its enclosing box (see Fig. 2). The volume of a nanocrystal is 43 πr3 , and its enclosing box is approximately 2r(2r + d)2 . The volume ratio of the enclosing box over the nanocrystal decreases with increasing siliconrich layer thickness, explaining the trend observed in Fig. 6 for relatively large silicon-rich layer thicknesses. However, for very small silicon-rich layer thicknesses, the optimal composition does not vary. To explain this, we
van Sebille et al. Nanoscale Research Letters (2016) 11:355
Page 6 of 7
thicknesses, the probability density functions broaden. The optimal composition depends on the probability of finding a nearest neighbor within a limited range (2 nm for SiOx ). For silicon-rich layer thicknesses up to approximately 2 nm, the probability of finding a nearest neighbor beyond 2 nm is negligible, as can be observed in Fig. 7. For these thicknesses, there is no reason to increase the nanocrystal density, since more closely packed nanocrystals will not increase the probability of finding a nearest neighbor within 2 nm. On the contrary, an increase in nanocrystal density will increase the probability of clustering.
Conclusions
Fig. 6 Optimal stoichiometry. The optimal stoichiometry for films with varying silicon-rich layer thicknesses for a sample with crystallinity of 1 (solid line) and 0.4 (dashed line). The buffer layer thickness is kept constant at 1 nm
look closer into their probability density functions, shown in Fig. 7. The peak at 1 nm is caused by the availability of nearest neighbors in the neighboring silicon-rich layers. The probability density function broadens for larger siliconrich layer thicknesses because its standard deviation is √ related to the 2D nanocrystal density by σ2D ∝ 1/ n2D [28]. Since n2D decreases for greater silicon-rich layer
We demonstrated an analytical method to optimize the composition of a silicon-rich layer for different crystallinities and thicknesses in order to achieve the highest density of non-touching and closely spaced silicon nanocrystals after annealing. The optimal stoichiometry depends on the crystallinity decreases for increasing silicon-rich layer thicknesses. However, for very small silicon-rich layer thicknesses, the optimal composition does not vary. This method can be used to find the best as-deposited composition in order to achieve optimal nanocrystal density and spacing after a subsequent annealing step.
Endnote 1 Note that the buffer layer composition is not stoichiometric and therefore will contain some excess silicon as well. This excess silicon will increase the amorphous Si– Si bond density, resulting in an underestimation of the crystallinity. However, since the buffer layer thickness is only 1 nm, compared to 3 nm for the silicon-rich layer, we expect this effect to be limited and assume it can be neglected. Acknowledgements We acknowledge financial support for this research from ADEM, A green Deal in Energy Materials of the Ministry of Economic Affairs of The Netherlands (http://www.adem-innovationlab.nl). We also thank ing. M.W.G.M. Verhoeven from the Eindhoven University of Technology for his assistance with the XPS measurements. Authors’ Contributions MS and JA initiated the research and performed the calculations. JA and JQ performed the sample fabrication and characterization. MS, RS, and MZ directed the experiments and the drafting of the paper. FT performed and analyzed the TEM measurements. MS wrote the manuscript. All authors read and approved the final manuscript.
Fig. 7 Nearest-neighbor distance probability density functions. The nearest-neighbor distance probability density functions for samples with silicon-rich layer thicknesses of 1.5, 2, and 3 nm and buffer layer thickness of 1 nm, calculated using their respective optimal compositions. The gray area depicts the range of desired nearest-neighbor distances. Shorter distances imply clustering, and greater distances lead to an insufficient tunneling probability
Competing Interests The authors declare that they have no competing interests. Author details 1 Photovoltaic Materials and Devices, Delft University of Technology, Mekelweg 4, 2628 CD Delft, The Netherlands. 2 Kavli Institute of Nanoscience, Quantum Nanoscience, Delft University of Technology, Lorentzweg 1, 2628 CJ Delft, The Netherlands.
van Sebille et al. Nanoscale Research Letters (2016) 11:355
Received: 28 May 2016 Accepted: 26 July 2016
References 1. Green MA, Cho EC, Cho Y, Huang Y, Pink E, Trupke T, Lin A, Fangsuwannarak T, Puzzer T, Conibeer G, Corkish R (2005) All-silicon tandem cells based on artificial semiconductor synthesised using silicon quantum dots in a dielectric matrix. In: Proceedings of the 20th European Photovoltaic Solar Energy Conference and Exhibition. pp 3–7 2. Heitmann J, Müller F, Zacharias M, Gösele U (2005) Silicon nanocrystals: size matters. Adv Mater 17(7):795–803. doi:10.1002/adma.200401126 3. Künle M, Kaltenbach T, Löper P, Hartel A, Janz S, Eibl O, Nickel KG (2010) Si-rich a-SiC:H thin films: structural and optical transformations during thermal annealing. Thin Solid Films 519(1):151–157. doi:10.1016/j.tsf.2010.07.085 4. Puthen-Veettil B, Patterson R, König D, Conibeer G, Green MA (2014) The impact of disorder on charge transport in three dimensional quantum dot resonant tunneling structures. J Appl Phys 116(16):163707. doi:10.1063/1.4899207 5. Schnabel M, Weiss C, Löper P, Wilshaw PR, Janz S (2015) Self-assembled silicon nanocrystal arrays for photovoltaics. Phys Stat Solidi A 212(8):1649–1661. doi:10.1002/pssa.201431764 6. Yi LX, Heitmann J, Scholz R, Zacharias M (2002) Si rings, Si clusters, and Si nanocrystals—different states of ultrathin SiOx layers. Appl Phys Lett 81(22):4248. doi:10.1063/1.1525051 7. López-Vidrier J, Hernández S, Hartel AM, Hiller D, Gutsch S, Löper P, López-Conesa L, Estradé S, Peiró F, Zacharias M, Garrido B (2011) Structural and optical characterization of size controlled silicon nanocrystals in SiO2 /SiOx Ny multilayers. Energy Procedia 10:43–48. doi:10.1016/j.egypro.2011.10.150 8. Conibeer G, Perez-Wurfl I, Hao X, Di D, Lin D (2012) Si solid-state quantum dot-based materials for tandem solar cells. Nanoscale Res Lett 7(1):193. doi:10.1186/1556-276x-7-193 9. Garrido B, Hernández S, Berencén Y, López-Vidrier J, Ramírez J, Blázquez O, Mundet B (2015) Electrical transport in Si-based nanostructured superlattices. In: Light Energy Harvesting with Group IV Nanostructures. Pan Stanford. pp 255–316. doi:10.1201/b18090-9. http://dx.doi.org/10. 1201/b18090-9 10. Mavilla NR, Chavan V, Solanki CS, Vasi J (2016) Study of temperature-dependent charge conduction in silicon-nanocrystal/SiO2 multilayers. Thin Solid Films 612:41–48. doi:10.1016/j.tsf.2016.05.029 11. Balberg I (2011) Electrical transport mechanisms in three dimensional ensembles of silicon quantum dots. J Appl Phys 110(6):061301. doi:10.1063/1.3637636 12. Jiang CW, Green MA (2006) Silicon quantum dot superlattices: modeling of energy bands, densities of states, and mobilities for silicon tandem solar cell applications. J Appl Phys 99(11):114902. doi:10.1063/1.2203394 13. Stegner AR, Pereira RN, Klein K, Lechner R, Dietmueller R, Brandt MS, Stutzmann M, Wiggers H (2008) Electronic transport in phosphorus-doped silicon nanocrystal networks. Phys Rev Lett 100(2):026803. doi:10.1103/physrevlett.100.026803 14. Lechner R, Stegner AR, Pereira RN, Dietmueller R, Brandt MS, Ebbers A, Trocha M, Wiggers H, Stutzmann M (2008) Electronic properties of doped silicon nanocrystal films. J Appl Phys 104(5):053701. doi:10.1063/1.2973399 15. Conibeer G, Green MA, König D, Perez-Wurfl I, Huang S, Hao X, Di D, Shi L, Shrestha S, Puthen-Veetil B, So Y, Zhang B, Wan Z (2010) Silicon quantum dot based solar cells: addressing the issues of doping, voltage and current transport. Prog Photovolt Res Appl 19(7):813–824. doi:10.1002/pip.1045 16. Smit C, van Swaaij RACMM, Donker H, Petit AMHN, Kessels WMM, van de Sanden MCM (2003) Determining the material structure of microcrystalline silicon from Raman spectra. J Appl Phys 94(5):3582. doi:10.1063/1.1596364 17. Schneider CA, Rasband WS, Eliceiri KW (2012) NIH image to ImageJ: 25 years of image analysis. Nat Methods 9(7):671–675. doi:10.1038/nmeth.2089 18. Gutsch S, Hiller D, Laube J, Zacharias M, Kübel C (2015) Observing the morphology of single-layered embedded silicon nanocrystals by using temperature-stable TEM membranes. Beilstein J Nanotechnol 6:964–970. doi:10.3762/bjnano.6.99
Page 7 of 7
19. Hiller D (2015) Silicon nanoclusters embedded in dielectric matrices: nucleation, growth, crystallization, and defects. In: Valenta J, Mirabella S (eds). Silicon nanoclusters embedded in dielectric matrices: nucleation, growth, crystallization, and defects. Pan Stanford Publishing. pp 99–144. doi:10.1201/b18090-5. http://dx.doi.org/10.1201/b18090-5 20. Iacona F, Bongiorno C, Spinella C, Boninelli S, Priolo F (2004) Formation and evolution of luminescent Si nanoclusters produced by thermal annealing of SiOx films. J Appl Phys 95(7):3723–3732. doi:10.1063/1.1664026 21. Daldosso N, Luppi M, Dalba G, Pavesi L, Rocca F, Priolo F, Franzò G, Iacona F, Degoli E, Magri R, Ossicini S (2003) Experimental and theoretical joint study on the electronic and structural properties of silicon nanocrystals embedded in SiO2 : active role of the interface region. In: MRS Proc Vol. 770. p 6. http://journals.cambridge.org/action/displayAbstract? fromPage=online&aid=8028079&fileId=S1946427400137005 22. Queeney KT, Weldon MK, Chang JP, Chabal YJ, Gurevich AB, Sapjeta J, Opila RL (2000) Infrared spectroscopic analysis of the Si/SiO2 interface structure of thermally oxidized silicon. J Appl Phys 87(3):1322. doi:10.1063/1.372017 23. Zimina A, Eisebitt S, Eberhardt W, Heitmann J, Zacharias M (2006) Electronic structure and chemical environment of silicon nanoclusters embedded in a silicon dioxide matrix. Appl Phys Lett 88(16):163103. doi:10.1063/1.2193810 24. Zhang XG (2001) Silicon/electrolyte interface. Springer. doi:10.1007/0-306-47921-4_2. http://dx.doi.org/10.1007/0-306-47921-4_ 2 25. Kessels WMM, Sanden MCMVD, Severens RJ, Ijzendoorn LJV, Schram DC (1998) Hydrogen in a-Si:H deposited by an expanding thermal plasma: a temperature, growth rate and isotope study. In: MRS Proceedings. Cambridge University Press Vol. 507. pp 529–534. doi:10.1557/proc-507-529. http://dx.doi.org/10.1557/proc-507-529 26. Sze SM, Ng KK (2006) Appendix G: properties of Si and GaAs. Wiley-Blackwell. doi:10.1002/9780470068328.app7. http://dx.doi.org/10. 1002/9780470068328.app7 27. Zhang XG (2001) Anodic oxide. Springer. doi:10.1007/0-306-47921-4_3. http://dx.doi.org/10.1007/0-306-47921-4_3 28. Tournus F (2011) Random nanoparticle deposition: inter-particle distances in 2d, 3d, and multilayer samples. J Nanopart Res 13(10):5211–5223. doi:10.1007/s11051-011-0506-9 29. Franzò G, Miritello M, Boninelli S, Lo Savio R, Grimaldi MG, Priolo F, Iacona F, Nicotra G, Spinella C, Coffa S (2008) Microstructural evolution of SiOx films and its effect on the luminescence of Si nanoclusters. J Appl Phys 104(9):094306. doi:10.1063/1.3006735 30. Conibeer G, Green MA, Corkish R, Cho Y, Cho EC, Jiang CW, Fangsuwannarak T, Pink E, Huang Y, Puzzer T, Trupke T, Richards B, Shalav A, Lin K-l (2006) Silicon nanostructures for third generation photovoltaic solar cells. Thin Solid Films 511–512:654–662. doi:10.1016/j.tsf.2005.12.119 31. López-Vidrier J, Hernández S, Hiller D, Gutsch S, López-Conesa L, Estradé S, Peiró F, Zacharias M, Garrido B (2014) Annealing temperature and barrier thickness effect on the structural and optical properties of silicon nanocrystals/SiO2 superlattices. J Appl Phys 116(13):133505. doi:10.1063/1.4896878 32. Laube J, Gutsch S, Wang D, Kübel C, Zacharias M, Hiller D (2016) Two-dimensional percolation threshold in confined Si nanoparticle networks. Appl Phys Lett 108(4):043106. doi:10.1063/1.4940971
