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 crosssection 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 sizedistributed ensembles.
Keywords:
Nanocrystals; Silicon nanocrystal ensemble; PhotoluminescenceBackground
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, electronhole 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 sizedependent absorption crosssection 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 siliconrich silicon oxide film with a thickness of 350 nm was deposited on a quartz substrate by a radio frequency cosputtering 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 SiO_{2} 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 S101411108s, 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 D^{−b} law, where 1 < b < 1.5 and D is the NC diameter. The relation between emission energy and size of NCs determined by highresolution 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, E_{g}, is given by the following equation:
where E_{Si} 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 sizedependent emission efficiency. For example, a single dangling bond on the surface will quench PL completely [10], and the chance of that will be strongly sizedependent. Other effects like exciton hopping [11,12] and nonradiative radiation via defects [13] can change the wavelengthdependent 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:
The NCs prepared by the procedure used in this study commonly have sizes whose logarithm is normally distributed [1416], described by Equation 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.
Figure 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 multiphoton absorption, undergo a fast nonradiative 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:
where N is the number of NCs in the illuminated volume that can emit a photon, σ_{abs} the absorption crosssection, and J_{pump} 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 crosssections for Si NCs of different sizes [18]. We can now write , where γ is a proportionality factor. In this case, Equation 4 can be written as follows:
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:
By converting E_{g} to wavelength with the relation , 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 crosssection, 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 powerdependent 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.
Figure 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).
Figure 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 sizedependent absorption crosssection, 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:
where β is the stretch parameter; τ_{PL}(E_{g}), the emission energydependent 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:
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 timedependent behavior, it is, thus, essential to remove contributions originating from the above nature.
Figure 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}(E_{g}), 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 [2224], and the importance of the absorption crosssection 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 sizedependent decay time could substantially alter the excitation powerdependent 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 nonradiative 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 nonlinear 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 crosssection 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 crosssection 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 timedependent 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 sizedependent 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 cowrote the paper. Both authors read and approved the final manuscript.
Authors’ information
DT and TG are postdoctoral 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

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

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 SiO_{2}.

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

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

Takeoka S, Fujii M, Hayashi S: Sizedependent photoluminescence from surfaceoxidized Si nanocrystals in a weak confinement regime.

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

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

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

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

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

Linnros J, Lalic N, Galeckas A, Grivickas V: Analysis of the stretched exponential photoluminescence decay from nanometersized silicon nanocrystals in SiO_{2}.

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

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

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

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.

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

Timmerman D, Izeddin I, Gregorkiewicz T: Saturation of luminescence from Si nanocrystals embedded in SiO_{2}.

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

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

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

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

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

Liu C, Kwon YK, Heo J: Laserinduced blueshift of the photoluminescence from PbS quantum dots in glasses.

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

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

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