The degree of crystallinity of samples are measured by the Raman spectrometer. The Raman spectra in the range of 200–1000 cm− 1 of the P-B co-doped SiNx:H thin films prepared at different H2 flow are presented Fig. 1a. It can be seen that for all the samples, a significant broad band locating at ~ 480–550 cm− 1 is identified. A sharp peak centered approximately at ~ 520 cm− 1 is also observed. It is the Raman spectral signature of crystalline Si component, and can be found for the samples deposited at H2 flow = 300, 400 and 500 sccm. At H2 flow = 200 sccm, it can be concluded that there is no c-Si component in the thin film due to the absence of sharp peak. When the H2 flow increases from 300 sccm to 500 sccm, the sharp peak emerges and gradually grows, indicating that the crystallization of P-B co-doped SiNx:H thin film is promoted as more H2 is introduced in dilution. This is in good agreement with previous reports [20].
Quantitative analysis of crystallinity degree relies on the crystalline volume fraction (Fc), which help us understand the structural evolvement of the deposited thin films. For calculating Fc, all the spectra of P-B co-doped SiNx:H thin films in the wavenumber range of 384–552 cm− 1 are fitted through Gaussian deconvolution method, and the typical best-fitting spectrum is shown in the Fig. 1b for the sample prepared at H2 flow = 300 sccm. There are four Gaussian-shaped sub-peaks can be divided from the spectrum: (1) a broad peak assigned to the asymmetric Si-N bond stretching mode at ~ 465 cm− 1 [21]; (2) a broad peak of the amorphous silicon (a-Si) phase at ~ 480 cm− 1; (3) a narrow peak denoting the intermediate phase containing grain boundaries and Si ultra-nanocrystalline at ~ 510 cm− 1 [22–23]; (4) a narrow Lorentz-shaped peak corresponding to the component of crystalline Si at ~ 520 cm− 1. Considering the peaks of intermediate phase and c-Si co-exist in spectra, it can be deduced that the c-Si content in the thin films consists of Si QDs, and the amorphous content surrounding the Si QDs is composed of amorphous silicon nitride and amorphous silicon since Si-N bond and a-Si co-exist in the thin films. Considering the intermediate phase as a crystalline content, Fc can be estimated by the equation [24]:
F c= (Ia + Ii)/(βIa + Ii + Ic)
where Ia, Ii and Ic denote the integrated intensities of the amorphous, intermediate and crystalline component respectively. Ia is defined as the sum of integrated intensities of peaks corresponding to amorphous silicon nitride and a-Si. The cross-section ratio β depends on the size of the Si QDs and the excitation wavelength, and can be expressed as [24]:
β(D) = 0.1 + exp (-D/250)
where D is the QD size in nm, and can be assumed as a unity in the case of the tiny nanocrystallites of a few nanometers. The Fc value varying with H2 is summarized in Table.1. At H2 flow = 200 sccm, the thin film is completely amorphous. When H2 flow achieves at 300 sccm, c-Si content forms with the Fc of 31.3%. Then as H2 flow increases, the Fc gradually grows. The maximum Fc value of 47.1% is obtained at H2 flow = 500 sccm.
In order to directly view the detailed internal microstructure of P-B co-doped SiNx:H thin films, TEM micrographs are shown in Fig. 2. It can be seen in Fig. 2a that the thin film deposited at H2 flow = 200 sccm has a cluster-like structure, where no Si QDs are founded. This phenomenon coincides well with the Raman spectra as the thin film is completely amorphous. Figure 2b shows the typical morphology of Si QDs which are formed in the P-B co-doped SiNx:H thin film deposited at H2 flow = 400 sccm, there can be observed numerous tiny spherical-shaped Si QDs with dense distribution, the QD size is estimated to be in the range of ~ 5–8 nm. The inset of Fig. 2b reveals the magnified details of microstructure by presenting the HRTEM image of Si QDs. The clearly visible crystallographic planes are observed for the Si QDs. The inter planer spacing of the well-defined crystallographic planes can be estimated to be ~ 3.0 Å, corresponding to the Si (111) lattice planes. Therefore, it is believable that the Si QDs are highly crystallized. The average QD size estimated from the TEM images has been presented in Table.1. As H2 flow grows from 300 to 500 sccm, the average QD size increases from 4.76 to 8.58 nm. It is reasonable since the Fc value has the same varying tendency.
Table I. Morphological, compositional and optical properties of P-B co-doped SiNx:H thin films
H2 Flow (sccm) | Fc (%) | Average QD size (nm) | Optical band gap (eV) |
200 | N/A | N/A | 2.78 |
300 | 31.3 | 4.76 | 2.24 |
400 | 34.4 | 6.17 | 1.85 |
500 | 47.1 | 8.58 | 1.68 |
FTIR spectroscopy is carried out to investigate the chemical bonding configurations of P-B co-doped SiNx:H thin films. Figure 3 shows the FTIR absorbance spectra of the samples deposited with various H2 flow. Different absorption bands corresponding to the typical vibration modes can be identified for all the samples, namely the Si–H wagging mode of a-Si at ~ 630 cm− 1, the Si–N stretching mode at ~ 840–870 cm− 1, the N–H rocking mode at ~ 1150 cm− 1, the Si–H stretching mode at ~ 2100–2200 cm− 1, and the N–H stretching mode at ~ 3300–3350 cm− 1 [25]. As H2 flow increases, the intensity of Si-N absorption band shrinks significantly, indicating the decreasing volume fraction of SiNx:H. No P- or B-related absorption bands are found, it is because the P and B concentrations in the in-situ deposited thin films are rather low.
The elemental composition and bonding states of the deposited thin films are investigated by XPS measurement. Figure 4a shows the Si 2p core level XPS spectra for all the deposited thin films. While a typical best-fitted spectrum of the sample deposited at H2 flow = 400 sccm is presented in Fig. 5a. As seen, five Gauss-shaped sub-peaks are extracted from the spectrum. According to previous works, the Si 2p peak in the range of ~ 99.8 eV to ~ 104 eV are determined by five different Si binding states Sin+ (n = 0–4), and n represents the number of nitrogen atoms bonding with one silicon atoms [26]. Thereinto, the Si0 state at ~ 99.8 eV denotes the chemical structure of pure Si, meaning the presence of c-Si QDs as for our samples. When combining the Raman results with the Si 2p core level XPS spectra, it can be concluded that the SiNx:H network is a mixture including four kinds of Si binding states. At H2 flow = 200 sccm, a peak is observed centering near the position of Si4+ state. Since no Si0 peak can be found in the spectrum, the thin film is believed to be completely amorphous, which coincides well with the Raman results. As H2 flow increases, the peak gradually shifts towards the low energy side, implying the increment in the proportion of c-Si component.
Figure 4b shows the XPS spectra of P 2p core level for the deposited thin films. At ~ 129.4 eV and ~ 133.3 eV, two peaks can be found for all the samples. The deconvolution of spectra is carried out, and the peak at ~ 129.4 eV can be separated to two sub-peaks. As shown in Fig. 5b, the spectrum of the sample deposited at H2 flow = 400 sccm is presented. The two sub-peaks are respectively assigned to P-Si (~ 129.3 eV) and P-P (130.4 eV) chemical bonds, while the peak at ~ 133.3 eV denotes the P-N bonds in the matrix or at the QD interface [27]. With the increasing H2 flow, the intensity of P-Si/P-P peak gradually enhanced while that of P-N peak shrinks, implying that the P-Si and P-P bonds tend to form in the case of large H2 flow, whereas large Ar flow favors the formation of P-N bonds. According to the XPS spectra, it is possible that P-doping is promoted as more H2 flow are introduced. It is notable that P-doping efficiency should not be directly related to P-Si peak, because successful P-doping depends on the P atoms penetrating into c-Si QDs and occupying Si sites, whereas the P-Si bonds formed at the QD interface or inside SiNx:H matrix make no contribution to effective P-doping. The conclusion is the same for B-doping. Additionally, no peak corresponding to P-B bonds is visible in the spectra, meaning that there are seldom P-B compounds formed in the thin films.
Figure 4c shows the XPS spectra of B 1s core level for all the thin films. It can be observed that there is a broad peak at ~ 191 eV. This peak is attributed to the B-N chemical bonds, and it can be found for the thin film deposited at H2 flow = 200 sccm. Since B-N bonds only form in the matrix or at the QD interface, the B-N peak correlates with the B atoms locating outside c-Si QDs. As H2 flow increases, the B-N peak shrinks while the B-Si peak at ~ 187.eV emerges and grows. Similarly, the spectra of P 2p core level are the same in varying tendency, meaning that large H2 flow is favorable for both P and B to be bonded to Si. Figure 5c shows the deconvolution of spectrum for the sample deposited at H2 flow = 400 sccm. Besides B-N and B-Si peaks, there is another broad sub-peak at ~ 185 eV. According to previous studies, the peak is considered to be caused by the plasmon loss peak of Si [28]. Similar as P 2p spectra, no peak assigned to P-B bonds is found.
The atomic concentration ratios of compositional elements are shown in Fig. 6. All the atomic concentrations are calculated based on the XPS spectra. Figure 6a shows the atomic concentrations of Si and N varying with H2 flow. As seen, the Si concentration increases significantly with the increasing H2 flow whereas N concentration exhibits a opposite varying tendency. We calculate the N/Si concentration ratio, and find it decreases from 73.79–18.14%. Therefore, we conclude that N atoms incorporate into P-B co-doped SiNx:H thin films more easily in the case of Ar-riched dilution. Correspondingly, a great amount of Si-N bonds perform as the barrier which interrupts the orderly growth of c-Si crystallographic planes, so the crystallization is suppressed. Figure 6b shows the atomic concentrations of P and B varying with H2 flow. It is easily found that both the P and B concentrations gradually decrease with the increasing H2 flow. As the H2 flow increases from 200 sccm to 500 sccm, the calculated P/Si concentration ratio decreases from 2.56–0.49%, while the calculated B/Si concentration ratio decreases from 1.80–0.43%. Similar to N, it is more easily for P and B atoms to incorporate in thin films in the case of Ar-riched dilution. However, these P and B tend to bond with N other than Si according to the weak P-Si and B-Si peaks shown in Fig. 4. Correspondingly, there is no improvement in doping efficiency in spite of the larger impurity concentration.
The optical band gap energy (Eopt) is estimated by measuring the absorbance and reflectance using the Tauc equation which is defined as αhν = B(hν-Eopt)2 [29]. α is the absorption coefficient, h is the Planck’s constant, ν is the frequency of the radiation, B is the edge width parameter. The Eopt value as a function of H2 flow is shown in Table.1. With the H2 flow increasing from 200 to 500 sccm, the Eopt value decreases from 2.78 eV to 1.68 eV. The Eopt value is determined by both the SiNx:H matrix and Si QDs for our samples. So it is reasonable that Eopt grows wider under the combined effects of more SiNx:H content and smaller c-Si QDs in the case of Ar-riched dilution.
Figure 7 represents the PL spectra of the P-B co-doped SiNx:H thin films deposited with different H2 flow. All the PL spectra are within a wide photo energy range of ~ 1.25–3.50 eV. Two peaks can be found respectively centering at ~ (P1) and (P2), it can be deduced from constant peak position that the two PL peaks are defect-related. According to previous studied, P1 at ~ 3.01 eV is attributed to the recombination from the conduction band to N4+ level, while P2 at ~ 2.41 eV originates from the radiative recombination at Si-dangling bond (K0) center [8]. Additionally, a P3 peak with varying center in the range of 1.45 eV ~ 2.11 eV is observed for the thin films deposited at H2 flow = 300, 400 and 500 sccm. It is prominent that the PL mechanism of P3 is different from that of P1 and P2. Based on the varying tendency of QD size, P3 can be attributed to QCE since the red-shift takes place with the increasing QD size as the H2 flow increases from 300 to 500 sccm. Note that the intensity of P3 gradually decreases with the increasing QD size. This probably results from the variation of P- and B-doping concentration in c-Si QDs and will be illustrated later.
The comparison of PL performance for the co-doped and un-doped SiNx:H thin films deposited at H2 flow = 400 sccm is presented in Fig. 8. The three peaks P1, P2 and P3 can be observed for both spectra in which all the peaks are larger in intensity for the un-doped sample. For the co-doped sample, a red-shift takes place for the QCE-related P3 peak, suggesting the impurity levels are involved in the optical transitions. When P and B are simultaneously doped in Si QDs, the transition between donor and acceptor levels which locates within band gap leads to lower PL energy. Therefore, the red-shift is a sign of the realization of P-B co-doping, so it is verified that the P-B co-doped c-Si QD/SiNx:H thin films have been successfully prepared by PECVD in-situ deposition method.