We demonstrate how to determine the size distribution of semiconductor nanocrystals in a quantitative manner using Raman spectroscopy employing an analytically defined multi-particle phonon confinement model. Results obtained are in excellent agreement with the other size analysis techniques like transmission electron microscopy and photoluminescence spectroscopy.
Analysis of the size distribution of nanocrystals is a critical requirement for the processing and optimization of their size-dependent properties. The common techniques used for the size analysis are transmission electron microscopy (TEM), X-ray diffraction (XRD) and photoluminescence spectroscopy (PL). These techniques, however, are not suitable for analyzing the nanocrystal size distribution in a fast, non-destructive and a reliable manner at the same time. Our aim in this work is to demonstrate that size distribution of semiconductor nanocrystals that are subject to size-dependent phonon confinement effects, can be quantitatively estimated in a non-destructive, fast and reliable manner using Raman spectroscopy. Moreover, mixed size distributions can be separately probed, and their respective volumetric ratios can be estimated using this technique. In order to analyze the size distribution, we have formulized an analytical expression of one-particle PCM and projected it onto a generic distribution function that will represent the size distribution of analyzed nanocrystal. As a model experiment, we have analyzed the size distribution of free-standing silicon nanocrystals (Si-NCs) with multi-modal size distributions. The estimated size distributions are in excellent agreement with TEM and PL results, revealing the reliability of our model.
Semiconductor nanocrystals draw attention as their electronic and optical properties can be tuned by simply changing their size in the range compared to their respective exciton-Bohr radii.1 These unique size-dependent features make these nanocrystals relevant for various technological applications. For instance, carrier multiplication effects, observed when a high energy photon is absorbed by the nanocrystals of CdSe, Si, and Ge, can be used in the concept of spectrum conversion in solar cell applications;2–4 or size-dependent optical emission from PbS-NCs and Si-NCs can be used in the light emitting diode (LED) applications.5,6 A precise knowledge and control on the nanocrystal size distribution will therefore play a determinant role on the reliability and the performance of these technological applications based on nanocrystals.
The commonly used techniques for the size distribution and morphology analysis of nanocrystals can be listed as X-ray diffraction (XRD), transmission electron microscopy (TEM), photoluminescence spectroscopy (PL), and Raman spectroscopy. XRD is a crystallographic technique that reveals morphological information of the analyzed material. From the broadening of the diffraction peak, estimation of the nanocrystal size is possible,7 however, obtaining a clear data is usually time consuming. Moreover, XRD can only enable the calculation of average of the nanocrystal size distribution. In the existence of multi-modal size distributions, size analysis with XRD can be misleading and result in wrong interpretations. TEM is a powerful technique that enables imaging of the nanocrystals.8 Although TEM is able to reveal the presence of individual distributions in a multi-modal size distribution, sample preparation issue is always an effort to be spent before the measurements. In addition, working on densely packed nanocrystal ensembles with different sizes is challenging because of the difficulty of individual nanocrystal imaging. Photoluminescence spectroscopy (PL) is an optical analysis technique, and optically active nanocrystals can be diagnosed. Nanocrystal size distribution is obtained from the size-dependent emission.9 Due to their poor optical properties of indirect band gap nanoparticles, large nanocrystals that are not subject to confinement effects, and defect-rich small nanocrystals cannot be detected by PL and the observed size distribution is only limited to nanocrystals with good optical properties. Although each of these abovementioned techniques has its own advantages, none of them have the capability of meeting the expectations (that is, being fast, non-destructive, and reliable) from and idealized size analysis technique.
Another means of size distribution analysis of nanocrystals is Raman spectroscopy. Raman spectroscopy is widely available in most of the labs, and it is a fast and non-destructive technique. In addition, in most cases, sample preparation is not required. Raman spectroscopy is a vibrational technique, which can be used to obtain information on different morphologies (crystalline or amorphous), and size-related information (from the size-dependent shift in the phonon modes that appear in the frequency spectrum) of the analyzed material.10 The unique feature of Raman spectroscopy is that, while size-dependent changes are observed as a shift in the frequency spectrum, the shape of the phonon peak (broadening, asymmetry) gives information on the shape of the nanocrystal size distribution. Therefore it is in principle possible to extract the necessary information, i.e., the mean size and the shape factor, from Raman spectrum to obtain the size distribution of nanocrystals analyzed. In the case of multi-modal size distributions sub-distributions can also be separately identified via deconvolution of the experimental Raman spectrum.
In the literature, two theories are commonly referred to model the effect of nanocrystal size distribution on the shape of the Raman spectrum. The bond polarizability model (BPM)11 describes the polarizability of a nanocrystal from the contributions of all the bonds within that size. The one-particle phonon confinement model (PCM)10 uses size-dependent physical variables, i.e., crystal momentum, phonon frequency and dispersion, and the degree of confinement, to define the Raman spectrum of a nanocrystal with a specific size. Since these physical variables depend on the size, an analytical representation of the PCM that can be explicitly formulized as a function of nanocrystal size can be defined. Projecting this expression on a generic size distribution function will therefore be able to account for the effect of size distribution within the PCM, which can be used to determine the nanocrystal size distribution from the experimental Raman spectrum.12
1. Planning of the Experiments
2. Raman Spectroscopy of the Nanocrystal of Interest
3. Size Distribution Determination of the Nanocrystal of Interest
For using Raman spectroscopy as a size analysis tool, a model to extract the size-related information from a measured Raman spectrum is needed. Figure 2 summarizes the analytical multi-particle phonon confinement model.12 All-size-dependent phonon confinement function (Figure 2c) is projected onto a generic size distribution function (Figure 2b), which is chosen as a lognormal distribution function. Given the amplitude (Figure 2d), full width at half maximum (Figure 2e), and the frequency shift (Figure 2f) value, this model can successfully be used to determine the size distribution.
Figure 3 envisages the use of multi-particle phonon confinement model to determine the size distribution of Si-NCs (details will follow). Si-NCs used in this analysis have a bimodal size distribution of small and large Si-NCs as shown in the TEM image.13 According to TEM size analysis (not shown here), small Si-NCs have a distribution in the range 2-10 nm, and large Si-NCs have a distribution in the range 40-120 nm. The analysis of the Raman spectrum at the left panel reveals that the size distribution of small Si-NCs are indeed in the range 2-10 nm. The distribution is lognormal with a mean size of 4.2 nm, and with a skewness (shape anisotropy factor) of 0.27.
Figure 4 represents a detailed comparative analysis of Si-NCs synthesized using different precursor flows in the plasma enhanced chemical vapor deposition (PECVD) system. For fitting the as-measured Raman data, we used two fitting functions knowing that we had two sub-distributions in the Si-NC mixture. Since Si does not show size-dependent peak-shift for the sizes larger than 20 nm, a bulk-like Lorentzian peak can be assigned for large Si-NCs, which are in the range 40-120 nm in this case (represented as “Large Si-NCs” in the plot). For small Si-NCs, we used the multi-particle phonon confinement model as the fitting function (represented as “Small Si-NCs” in the plot). The mean size, and the skewness of the size distribution are obtained from this fit, which are the necessary parameters to plot the size distribution demonstrated in Figure 2b. This fitting function can be integrated up to the size from which a peak-shift is not observed any more, i.e., 20 nm for Si-NCs. Results show that we can successfully determine the mean size, skewness, and the complete size distribution of Si-NCs (panel c and d) using Raman spectroscopy. Moreover, the volume fraction of small Si-NCs and large Si-NCs can be determined by the ratio of integrated peak areas. For Si-NCs synthesized using 3 sccs (standard cubic centimeters per second) of SiH4 flow, the volume fraction of small Si-NCs was 80%, while for the case of 10 sccs SiH4 flow, small Si-NC volume fraction is 88%.
Figure 5 demonstrates the comparison of determined average particle size of Si-NCs from various techniques. First, our analytical-PCM12 (stars) is in very good agreement with the PCM.10 Second, the results obtained from Raman spectroscopy are in good agreement with the results obtained from transmission electron microscopy (TEM) and photoluminescence spectroscopy (PL) (the PL size distribution is obtained using the model of Delerue et al.16). This proves the reliability of using Raman spectroscopy with the analytical-PCM for size analysis of Si-NCs. In addition, we also demonstrate the BPM,11 which is also used for size analysis of semiconductor nanocrystals. Figure 5 concludes also that the PCM predicts the size of a Si-NC from its size-dependent Raman shift better than the BPM does.
Figure 1. Representation of nanoparticles and the Raman spectrometer. a) Si-NCs deposited in an Ar/SiH4 gas mixture on plexiglass substrates using a PECVD tool. Si-NCs are in the form of a powder. The tone differences on the substrate are due to the differences in morphologies and thicknesses of Si-NC powder, which are exposed to different plasma regions during the synthesis13. As deposited morphologies of Si-NCs are ready for Raman spectroscopy measurements. The short side of the substrate is 2 cm. b) Reference sample, i.e., crystalline Si wafer, measurement in order to observe the Raman peak position of bulk Si. This information will be used as a reference point when determining the relative shift of Si-NCs from their bulk peak position. c) Image of the Raman spectrometer used for size determination studies. d) The screenshot of the software for performing and recording the data to be analyzed.
Figure 2. Formulas used in the analysis of size distribution of Si-NCs. a) Raman intensity of Si-NCs with size distribution. b) The generic size distribution function to determine Si-NC size distribution. c) Analytical representation of one-particle PCM for a Si-NC with a size D. d) The amplitude, e) the full width half maximum, and f) the vibration frequency representations of a Si-NC with a size D, which explicitly appear in c).
Figure 3. From Raman spectroscopy analysis to nanocrystal size distribution. As-measured data from Raman spectroscopy can be converted to a quantitative size distribution of nanocrystals using the multi-particle analytical-PCM.
Figure 4. Size and volume fraction analysis of Si-NCs. Raman spectrum of Si-NCs synthesized in a PECVD tool using a) 3 sccs and b) 10 sccs of SiH4(silane) gas flow, respectively. a) and b) demonstrates the deconvolution route for small and large Si-NCs. Deconvolution is done using a Lorentzian peak for bulk-like large Si-NCs and multi-particle analytical-PCM for small Si-NCs. Corresponding size distributions and volume fractions of small Si-NCs for 3 and 10 sccs SiH4 flow are demonstrated in c) and d), respectively. The average size of small Si-NCs is 4.2 nm with a skewness of 0.26 for panel c) and 3.7 nm with a skewness of 0.30 for panel d). The volume fractions are estimated as 80% and 88% for panel c) and d), respectively.
Figure 5. Comparison of size distribution of Si-NCs from various techniques. Size analysis of Si-NCs using various techniques (TEM and PL16) and size analysis using Raman spectroscopy resulted in excellent agreement. The results also demonstrate that the PCM leads to a more accurate size determination with respect to the BPM. This figure has been modified from Ref. 12 by permission of American Institute of Physics.
First discussion point is the critical steps within the protocol. In order not to have overlapping peaks with the material of interest, it is important to use another type of substrate material as mentioned in step 1.2. For instance, if Si-NCs are of interest, do not use silicon substrate for the Raman measurements. In Figure 1a, for instance, Si-NCs were synthesized on plexiglass substrates, which has completely flat signal roughly around the range of interest, i.e., 480-530 cm–1. In addition to measuring a bulk reference signal to estimate the shift of nanocrystal-related peak as mentioned in step 1.4, it is also critical to locate the exact position of the peak of the bulk material, since it can be also shifted depending on the operating and environment conditions of Raman spectroscope. For the case of Si-NCs, the reference sample is a crystalline silicon wafer, which is known to have the transverse optical (TO) mode at 521 cm–1.12 However, this can be shifted as result of temperature of the operating laser, which is liked with the ambient conditions, and intensity. Therefore, it is important to record a reference data every time before the measurements, and correct them with respect to the known peak positions from the literature. A high laser power can heat up the nanocrystals, and change their sizes, which result in a laser-induced shift in the Raman spectrum. Therefore, it is critical to determine the maximum laser power that can be used safely following the instructions at step 1.5. If inserting the estimated average size and skewness in step 3.4.3 into the fitting function cannot not cover the shape of the measured Raman spectrum, it means that the material analyzed is composed of sub-distributions. Then repeat the steps 3.4.2 and 3.4.3 again to add additional peaks for the fitting. The determination of the integral boundaries (step 3.4.7) for the size distribution is another critical point. The boundaries of the integral in the distribution function represents the smallest and the largest sizes in the nanocrystal size distribution. 1 nm is the smallest stable size for most of the nanocrystal systems.17 As the amount of size-dependent shift decreases with an increase of the nanocrystal size, set the upper limit of the integration to any size that does not exhibit any shift in the Raman spectrum (step 3.4.8). For instance, Si-NCs with sizes larger than 20 nm does not exhibit any shift in the Raman spectrum, they resemble bulk-like behavior.10,12 Therefore, for Si-NCs setting the upper boundary of the integral with any size greater than 20 nm will not change the result.
The PCM describes the size-dependent Raman peaks of nanocrystals through a complicated expression. In fact, the PCM depends on the phonon dispersion, the confinement function, and the vibration frequency, which all implicitly depend on the size, D. Furthermore, the PCM is for a specific size, and to determine the size distribution, it has to be projected on to a generic distribution function and get integrated over a size range. This procedure is complicated and up to now, experimentalists were using Raman spectroscopy mostly to determine the average size of Si-NCs from the shift of the as-measured Raman peak. On the other hand, the multi-particle analytical-PCM that we have formulized contains the size, D, as an explicit parameter and makes it possible to determine not only the average size, but also the complete distribution and it shape in a simple way using Raman spectroscopy.
As already noticed, Raman spectroscopy is able to determine the size distribution of nanocrystals within the confinement limit, which is about 20 nm for Si. Larger sizes cannot be analyzed for their size distribution as they do not exhibit a size-dependent feature in Raman spectrum, i.e., they have similar peak shapes and positions as bulk crystalline Si. This limitation holds for any type of nanocrystal system that exhibits size-dependent peak-shift in the Raman spectrum. However, the confinement limit can vary depending on the nanocrystal system. For instance, the confinement limit for Ge-NCs is about 15 nm.18
A concern during the size analysis of nanocrystals is the shift of the Raman peaks of nanocrystals because of the secondary reasons, which could result in wrong size interpretations using Raman spectroscopy. These secondary reasons are structural modification (grain growth or shape change) of nanocrystals under excessive laser power during Raman measurement, and the stress induced by the matrix, in which nanocrystals are embedded (if any). In order to avoid excessive laser heating, it is advised to start the measurement with the lowest possible laser power and increase it gradually to establish a clear signal. As long as the peak shape and position stays stable independent of the laser power used, it can be considered that the laser power used is in the safe limit.13 In addition, the limit beyond which, a heating related broadening of the Raman peaks is observed is defined in the literature as the Fano broadening.12,19 As long as the Fano broadening limit is not reached, nanoparticles are not subject to heating related modifications. If the stress-related peak-shift is present, it is unavoidable, and has to be accounted for before determining the size-dependent peak-shift. The amount of stress can be determined using X-ray diffraction (XRD), where a shift in the diffraction peak position is a measure of stress. The estimated stress can then be accounted in the Raman spectrum with an addition of the stress-dependent peak-shift term.18 If no stress is observed from XRD, direct analysis of size distribution can be performed from the as-measured Raman spectrum. Since all the nanoparticles have stable structures, as observed previously,14 the crystallinity is well established, and concerns related with poor crystalline structures are ruled out for the Raman analysis.
The multi-particle analytical-PCM demonstrated in Figure 2a is flexible in terms of the distribution function and the confinement function used. For instance, any type of generic distribution function can be replaced with lognormal distribution function, i.e. it could be a lognormal, normal, or a logistic function, without affecting the phonon confinement function. In addition, the analytical-PCM demonstrated in Figure 2c can be re-defined depending on the type of the material to be used. Some nanocrystal systems that exhibit size-dependent Raman peak-shifts (so that their size distributions can also be determined using Raman spectroscopy) are Ge-NCs,20 SnO2-NCs,21 TiO2-NCs,22 and diamond-NCs.23
The authors have nothing to disclose.
This work was part of the research programme of the Foundation for Fundamental Research on Matter (FOM), which is part of the Netherlands Organisation for Scientific Research (NWO). Authors of this work thank M. J. F. van de Sande for skillful technical assistance, M. A. Verheijen for TEM images, and the group of Tom Gregorkiewicz for PL measurements.
Raman Spectroscopy | Renishaw | In Via | Equipped with 514 nm Ar ion laser |
Wire 3.0 | Renishaw | Raman spectroscopy record tool | |
Mathematica | Wolfram | For fitting function and size determination | |
Substrate | Plexiglass (to avoid signal coincidence with Si-NCs) | ||
Si wafer | Reference to Si-NC peak position | ||
Photoluminescence Spectroscopy | 334 nm Ar laser. For optical size distribution. | ||
Transmission Electron Microscopy | Beam intensity 300 kV. For nanocrystal size and morphology determination. |