Open Access Nano Express

Power-dependent spectral shift of photoluminescence from ensembles of silicon nanocrystals

Dolf Timmerman* and Tom Gregorkiewicz

Author Affiliations

Van der Waals-Zeeman Institute, University of Amsterdam, Science Park 904, Amsterdam, NL-1098 XH, The Netherlands

For all author emails, please log on.

Nanoscale Research Letters 2012, 7:389  doi:10.1186/1556-276X-7-389


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/7/1/389


Received:18 April 2012
Accepted:22 June 2012
Published:12 July 2012

© 2012 Timmerman and Gregorkiewicz; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Nanocrystals are widely studied for their tunable optical properties, most importantly increased luminescence efficiency and emission energy. Quantum confinement effects are found for many different types of nanocrystals, and these introduce a relation between the emission wavelength and the size of nanocrystals. When ensembles of nanocrystals with a distribution of sizes are studied, this can have profound effects on their luminescence spectra. Here, we show how photoluminescence spectra of ensembles of silicon nanocrystals can shift under different excitation conditions, resulting from differences in absorption cross-section of the individual nanocrystal sizes. This effect, together with the fact that after a pulsed excitation a silicon nanocrystal can only emit a single photon, determines how the distribution of excited nanocrystals changes and leads to the spectral shift for different excitation powers. Next to this effect, the influence of different radiative rates in such ensembles is also addressed. These notions are important for the interpretation of photoluminescence data for silicon nanocrystals but can be extended to any nanoparticle system comprising size-distributed ensembles.

Keywords:
Nanocrystals; Silicon nanocrystal ensemble; Photoluminescence

Background

Nanoparticles of many different semiconductor materials are intensively studied for their interesting optical properties. Quantum confinement in these nanostructures leads to strong altering of radiative rates and energy structure. In general, stronger confinement in smaller particles gives rise to shorter emission wavelengths [1]. Due to limitations in size selectivity on this scale, ensembles of nanoparticles always show a distribution of sizes and shapes. This results in inhomogeneously broadened photoluminescence (PL) spectra. For indirect band gap semiconductors, there is a special interest in the use of nanocrystals (NCs). In the bulk equivalent of these material, electron-hole radiative recombination has very low probability, since the difference in crystal momentum between the top of the valence band and bottom of conduction band can not be compensated by photon emission. In this case, optical transitions need to be accompanied by phonons, lowering their probability. For example, the radiative lifetime of excitons in bulk silicon is in the order of milliseconds. In Si NCs, the momentum conservation rule is relaxed, as k is not a good quantum number any more, and radiative phononless transitions become more probable. This leads to a substantial decrease of the radiative recombination time, which can be observed by photoluminescence lifetime measurements (from milliseconds to microseconds for the smallest Si NCs) [2,3].

PL spectroscopy is a powerful tool that is often used to determine specific properties of NCs. Most importantly, once a relationship has been determined between the size and the associated photon energy, it can be used as a relatively simple experimental way to determine sizes and distributions hereof.

In this work, we study PL spectra of silicon NCs and their change for different excitation fluences. A combination of size-dependent absorption cross-section and saturation effects in these NCs has a profound effect on the measured PL spectrum. We show simulations of this effect and compare them with experimental data. The effect of excitation conditions on the spectral profile will be discussed. Next to this, the influence of different PL lifetimes present in an ensemble of NCs is also addressed.

Methods

A silicon-rich silicon oxide film with a thickness of 350 nm was deposited on a quartz substrate by a radio frequency co-sputtering technique. Subsequent annealing of the sample at a temperature of 1,150°C in nitrogen atmosphere for 30 min induced formation of Si NCs. The layer of Si NCs in amorphous SiO2 created in this way was characterized by an excess Si amount of 5 vol.% and contained NCs with an average diameter of approximately 3.5 nm.

PL experiments were performed under pulsed excitation provided by an optical parameter oscillator, pumped by the 3rd harmonic of a Nd:YAG laser, with a wavelength of 416 nm (3 eV) and pulse width of approximately 9 ns. Luminescence was dispersed by a monochromator (Solar M266, SOLAR Laser Systems, Minsk, Republic of Belarus) and detected by a linear detector (Hamamatsu S10141-1108s, Hamamatsu Photonics K.K., Shizuoka Pref., Japan). The entire system was calibrated for its spectral response with a calibration source. All measurements were taken at room temperature.

Results

There is a lot of literature available on the dependence of the band gap of the Si NCs on the diameter. For an overview hereof, see [4]. For oxidized NCs, it is commonly reported that this dependence follows a Db law, where 1 < b < 1.5 and D is the NC diameter. The relation between emission energy and size of NCs determined by high-resolution transmission electron microscopy for the sample used in this study has been established earlier [5]. It was found that a good fit to the experimental data for the optical gap, Eg, is given by the following equation:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/389/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/389/mathml/M1">View MathML</a>

(1)

where ESi is the band gap of bulk Si, similar as the theoretical dependence determined in [6,7]. When an intrinsic yield of PL of 100% is assumed for all sizes [8,9], we can obtain a size distribution directly from PL spectrum, as the emission wavelength is related to the band gap. It is important to note that this is only valid for the optically active NCs in the entire distribution. Different effects can quench or enhance luminescence from Si NCs, which could introduce a size-dependent emission efficiency. For example, a single dangling bond on the surface will quench PL completely [10], and the chance of that will be strongly size-dependent. Other effects like exciton hopping [11,12] and non-radiative radiation via defects [13] can change the wavelength-dependent quantum efficiencies and, thus, the final emission spectrum. If such effects are quantified for a certain system, this could be used in order to get information about the entire distribution from the determined PL spectra. For the rest of the analysis done here, however, we assume that all NCs are optically active.

Rewriting Equation 1 gives the relation between diameter and wavelength:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/389/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/389/mathml/M2">View MathML</a>

(2)

The NCs prepared by the procedure used in this study commonly have sizes whose logarithm is normally distributed [14-16], described by Equation 3.

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/389/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/389/mathml/M3">View MathML</a>

(3)

where μ and σ are the mean and standard deviation, respectively, for the natural logarithm of the diameter D. By fitting the determined size distribution with Equation 3, we can determine these important parameters. Figure 1 shows the result of this analysis, where values of μ and σ of 1.496 and 0.296, respectively, show a good fit with the experimentally determined PL spectrum. We note here that the fit actually deviates for larger NCs and underestimates its contribution. This is a result of the fixed limit that is imposed on the minimum optical gap (Equation 1), whereas in realistic systems, this will be less strict by participation of phonons in emission.

thumbnailFigure 1. Spectrum and size distribution. PL spectrum obtained under high fluence excitation. Inset: The black line indicates the size distribution of the ensemble by converting the wavelengths from the measured PL spectrum with Equation 1. The dashed red line shows a fit of this distribution with Equation 3 with the fit parameters as shown in the figure.

In reference [17], it was shown that, under the condition that the excitation pulse is much shorter than the PL lifetime, every NC can emit a single photon only. Multiple excitons localized in the same NC, that have been excited by multi-photon absorption, undergo a fast non-radiative Auger recombination. Finally, only a single exciton remains for radiative recombination. The number of excited NCs N after an excitation pulse can, in this case, be described by Equation 4:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/389/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/389/mathml/M4">View MathML</a>

(4)

where N is the number of NCs in the illuminated volume that can emit a photon, σabs the absorption cross-section, and Jpump the pump fluence. From geometrical considerations, it is expected that the σabs scales with the square of the diameter of the NC, which is theoretically supported by calculations of optical absorption cross-sections for Si NCs of different sizes [18]. We can now write <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/389/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/389/mathml/M5">View MathML</a>, where γ is a proportionality factor. In this case, Equation 4 can be written as follows:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/389/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/389/mathml/M6">View MathML</a>

(5)

We can use now the expression for the distribution of NC sizes from Equation 3 and combine this with Equation 5 to arrive to an expression for the number of excited NCs as a function of excitation pump fluence and their associated band gap energy:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/389/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/389/mathml/M7">View MathML</a>

(6)

By converting Eg to wavelength with the relation <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/389/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/389/mathml/M8">View MathML</a>, the associated PL spectrum can be obtained. Qualitatively, this equation describes how the ratio of excited NCs for different sizes changes upon illumination with different excitation powers. So, under assumption that the PL quantum efficiency is identical for all sizes, it shows how the PL spectrum changes with laser power. Larger NCs have a larger absorption cross-section, so their emission will be saturated at smaller pump fluences than for the smaller NCs. This means that, when the pump fluence is increased even more, the smaller NCs can still increase their contribution to the total (ensemble) emission. As a result, the PL spectrum will shift to shorter wavelengths. Figure 2 presents PL spectra obtained with Equation 6 for six different values of pump fluence of 0.01, 0.03, 0.07, 0.3, 0.65, and 1. While the intensity increases, saturation sets in at longer wavelengths, and the peak shifts to shorter wavelengths. A comparison of experimental data and the simulations is shown in Figure 3. We note here that carrier multiplication, known to occur in similar systems [19], has, in this case, been eliminated by using an excitation photon energy below the threshold of two times the NC band gap. For clarity, the normalized spectra are shown. In the simulations, the same ratios of excitation fluence are used as in the experiment. The value of γ was determined from the saturation behavior at a single wavelength by the method presented in references [17,20]. In this case, Equation 4 is fitted to the power-dependent PL intensity to obtain σabs for a specific wavelength, which is subsequently used to deduce a value for γ. The simulation models the experimental spectra well; the peak shifts to shorter wavelengths, and the shift on the short wavelength side is more pronounced than that of longer wavelengths. The position of the peaks determined from the simulation and of the experiments shows a linear dependence. The small difference in the absolute numbers here is a result of the underestimation of contribution of larger NCs and, thus, longer wavelengths, which has been discussed earlier.

thumbnailFigure 2. Simulated power dependent spectra. Spectra obtained from Equation 6 for pump fluence of 0.01 (red), 0.03 (orange), 0.07 (green), 0.3 (light blue), 0.65 (blue), and 1 (pink).

thumbnailFigure 3. Experimental and simulated power dependent spectra (normalized). Left: Experimentally obtained normalized PL spectra for different excitation powers. The ratio of these powers is the same as the simulations in the left panel. Right: Simulated evolution of the normalized PL spectra for ratio’s of excitation power identical to that of Figure 2. Inset: Comparison of experimental and simulated peak wavelengths for the different excitation powers.

Next to a size-dependent absorption cross-section, there is also a dependence of the decay time [2]. Typically, smaller Si NCs have PL lifetimes in the order of 10 μs, while for large NCs, this increases to a few hundreds μs. This effect leads to red shift of PL spectra for longer detection times, which is often observed experimentally [21]. For Si NCs, most commonly a stretched exponential decay of PL intensity at a single wavelength is observed:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/389/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/389/mathml/M9">View MathML</a>

(7)

where β is the stretch parameter; τPL(Eg), the emission energy-dependent decay time. In the regime where all NCs are excited, we can take the initial distribution described by Equation (3) and combine it with Equation (7) to arrive to an expression that describes the evolution of PL spectrum in time:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/389/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/389/mathml/M10">View MathML</a>

(8)

For the wavelengths in the range of 600 to 850 nm, we have taken PL decays, and the traces have been fitted with the stretched exponential function to determine τPL and β. It was found that β = 0.55 over the entire range, and the PL lifetime increases exponentially with detection wavelength in this region (see inset to Figure 4). For longer wavelengths, the τPL was determined by extrapolation in order to be able to simulate spectra at different times after excitation. Figure 4 shows these results for t = 0, 10, 30, 50, 80, 100 and 150 μs. A clear shift of the maximum PL intensity is accompanying the change of delay, very similar as observed in reference [21]. In the interpretation of such a time-dependent behavior, it is, thus, essential to remove contributions originating from the above nature.

thumbnailFigure 4. Spectral evolution in time. Simulated spectra for different times after a pulsed excitation of the entire ensemble of NCs, acquired from Equation 8. The spectral dependence of PL decay time, τPL(Eg), has been obtained experimentally. Left inset: PL intensity decay for detection wavelength of 800 nm. The red line indicates a fit with a stretched exponential function. Right inset: Experimentally determined PL lifetimes for different wavelengths (black squares). The red line shows an extrapolation for longer wavelengths.

Discussion

The shift of PL spectra to shorter wavelengths when increasing the excitation power has been shown to occur in different types of NC ensembles. Different explanations have been suggested in order to account for this behavior [22-24], and the importance of the absorption cross-section has been noted as well [16,25]. Although the above proposed framework does not necessarily fully explain the experimentally observed shifts of PL spectra, it is important to account for the strong influence of the discussed effects. It should be emphasized that the analysis and simulations done here are based on the condition that the system is excited by a short pulse, with respect to the PL lifetime, and that the emerging luminescence is decaying fully in between the subsequent pulses. In experimental practice, this means that the pulse repetition rate should be small enough to allow all excited NCs to return to the ground state. Under continuous excitation, typical for the majority of experimental reports, the PL characteristics change considerably, as equilibrium conditions are important in this case. The size-dependent decay time could substantially alter the excitation power-dependent emission, and an even larger blueshift can be expected as the shorter PL lifetime of small NCs allows for more radiative recombinations. Now, the question arises under which experimental conditions PL spectra should be collected in order to get reliable information about the physical parameters of the system. For NCs showing long decay times, like Si, it is possible to excite the entire distribution by a pulsed excitation. After this, all multiple carriers simultaneously present in the same NC undergo a fast Auger non-radiative recombination, and every NC can, in principle, emit a photon. This would then be the best way to obtain information about the entire ensemble. However, high pump fluences might be necessary to arrive to this situation, and non-linear and heating effects can appear and influence PL properties [26]. Alternatively, at low powers where the average amount of absorbed photons per NC <<1, the shape of the PL spectrum will not change for varying power. This would be the best region to do comparative experiments, but the low intensity of PL signal might complicate the practical usability of this approach. It is also important to note that, in this case, PL will be dominated by contribution from large NCs. For more accurate information about smaller particles, a scaling based on absorption cross-section should be employed.

Conclusions

To sum up, it has been shown that PL spectra of typical Si NC ensembles show a relatively strong dependence on the excitation fluence. A shift of 30 nm in the maximum of PL intensity appears between spectra taken at a low power excitation and under saturation conditions. It has been demonstrated that this shift can be modeled fairly well by taking into account the differences in absorption cross-section for different sizes of NCs. This, and the notion that each Si NC can only emit a single photon after a pulsed excitation, determines the distribution of the ensemble that is excited and contributes to PL. The evolution of the excited state distribution with respect to excitation pump fluence can be modeled fairly well by these simulations. Additional simulations, taking into account also differences in PL lifetime in the ensemble, show spectral narrowing and red shift for peak intensity in time-dependent spectra. While the presented modeling demonstrates the importance of excitation conditions on the PL spectra of ensembles of Si NCs, this approach can be extended to all semiconductor NCs, where ensembles and size-dependent absorption are imminent.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

DT conceived the idea, developed theoretical descriptions and performed simulations. TG supervised the project. DT and TG co-wrote the paper. Both authors read and approved the final manuscript.

Authors’ information

DT and TG are post-doctoral researcher and full professor, respectively, at the van der Waals - Zeeman Institute at the University of Amsterdam.

Acknowledgements

The authors would like to thank The Netherlands Organisation for Scientific Research (NWO) for financial support. They also want to acknowledge Brent Huisman for the contribution in spectroscopic measurements.

References

  1. Alivisatos A: Semiconductor clusters, nanocrystals, and quantum dots.

    Science 1996, 271:933-937. OpenURL

  2. Garcia C, Garrido B, Pellegrino P, Ferre R, Moreno JA, Morante JR, Pavesi L, Cazzanelli M: Size dependence of lifetime and absorption cross section of Si nanocrystals embedded in SiO2.

    Appl Phys Lett 2003, 82:1595-1597. OpenURL

  3. de Boer WDAM, Timmerman D, Dohnalová K, Yassievich IN, Zhang H, Buma JW, Gregorkiewicz T: Red spectral shift and enhanced quantum efficiency in phonon-free photoluminescence from silicon nanocrystals.

    Nat Nano 2010, 5:878-884. OpenURL

  4. Belyakov V, Burdov VA, Lockwood R, Meldrum A: Silicon nanocrystals: fundamental theory and implications for stimulated emission.

    Adv Opt Technol 2008, 2008:279502. OpenURL

  5. Takeoka S, Fujii M, Hayashi S: Size-dependent photoluminescence from surface-oxidized Si nanocrystals in a weak confinement regime.

    Phys Rev B 2000, 62:24. OpenURL

  6. Delerue C, Allan G, Lannoo M: Theoretical aspects of the luminescence of porous silicon.

    Phys Rev B 1993, 48:11024. OpenURL

  7. Wang LW, Zunger A: Electronic structure pseudopotential calculations of large (∼1000 atoms) Si quantum dots.

    J Phys Chem 1994, 98:2158-2165. OpenURL

  8. Miura S, Nakamura T, Fujii M, Inui M, Hayashi S: Size dependence of photoluminescence quantum efficiency of Si nanocrystals.

    Phys Rev B 2006, 73:245333. OpenURL

  9. Walters RJ, Kalkman J, Polman A, Atwater HA, de Dood MJA: Photoluminescence quantum efficiency of dense silicon nanocrystal ensembles in SiO2.

    Phys Rev B 2006, 73:132302. OpenURL

  10. Lannoo M, Delerue C, Allan G: Theory of radiative and nonradiative transitions for semiconductor nanocrystals.

    J Lumin 1996, 70:170-184. OpenURL

  11. Linnros J, Lalic N, Galeckas A, Grivickas V: Analysis of the stretched exponential photoluminescence decay from nanometer-sized silicon nanocrystals in SiO2.

    J Appl Phys 1999, 86:6128. OpenURL

  12. Pavesi L: Influence of dispersive exciton motion on the recombination dynamics in porous silicon.

    J Appl Phys 1996, 80:216. OpenURL

  13. Zatryb G, Podhorodecki A, Misiewicz J, Cardin J, Gourbilleau F: On the nature of the stretched exponential photoluminescence decay for silicon nanocrystals.

    Nanoscale Res Lett 2011, 6:106. OpenURL

  14. Yorikawa H, Muramatsu S: Logarithmic normal distribution of particle size from a luminescence line-shape analysis in porous silicon.

    Appl Phys Lett 1997, 71:644-646. OpenURL

  15. Ledoux G, Guillois O, Porterat D, Reynaud C, Huisken F, Kohn B, Paillard V: Photoluminescence properties of silicon nanocrystals as a function of their size.

    Phys Rev B 2000, 62:15942-15951. OpenURL

  16. Vinciguerra V, Franzò G, Priolo F, Iacona F, Spinella C: Quantum confinement and recombination dynamics in silicon nanocrystals embedded in Si/SiO2 superlattices.

    J of Appl Phys 2000, 87:8165-8173. OpenURL

  17. Timmerman D, Izeddin I, Gregorkiewicz T: Saturation of luminescence from Si nanocrystals embedded in SiO2.

    Phys Stat Sol 2010, 207:183-187. OpenURL

  18. Poddubny AN, Prokofiev AA, Yassievich IN: Optical transitions and energy relaxation of hot carriers in Si nanocrystals.

    Appl Phys Lett 2010, 97:231116. OpenURL

  19. Timmerman D, Valenta J, Dohnalova K, de Boer WDAM, Gregorkiewicz T: Step-like enhancement of luminescence quantum yield of silicon nanocrystals.

    Nat Nanotech 2011, 6:710-713. OpenURL

  20. Timmerman D, Izeddin I, Stallinga P, Yassievich IN, Gregorkiewicz T: Space-separated quantum cutting with Si nanocrystals for photovoltaic applications.

    Nat Photonics 2008, 2:105-109. OpenURL

  21. Gusev OB, Prokofiev AA, Maslova OA, Terukov EI, Yassievich IN: Energy transfer between silicon nanocrystals.

    JETP Lett 2011, 93:147-150. OpenURL

  22. Kanemitsu Y, Ando M, Matsuura D, Kushida T, White CW: Photoluminescence from CdS nanocrystals fabricated by sequential ion implantation.

    J Lumin 2001, 94-95:235-238. OpenURL

  23. Liu C, Kwon YK, Heo J: Laser-induced blue-shift of the photoluminescence from PbS quantum dots in glasses.

    Chem Phys Lett 2008, 452:281-284. OpenURL

  24. Wen X, Dao LV, Hannaford P, Cho E, Cho YH, Green MA: Excitation dependence of photoluminescence in silicon quantum dots.

    New J Phys 2008, 9:337. OpenURL

  25. Jayatilleka H, Diamare D, Wojdak M, Kenyon AJ, Mokry CR: Probing energy transfer in an ensemble of silicon nanocrystals.

    J Appl Phys 2011, 110:033522. OpenURL

  26. Gardelis S, Nassiopoulou AG, Vouroutzis N, Frangis N: Effect of exciton migration on the light emission properties in silicon nanocrystal ensembles.

    J Appl Phys 2009, 105:113509. OpenURL