
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

72560
10.1038/s41598-024-72560-3
Article
Properties of phosphorus-boron co-doped c-Si quantum dots/SiNx:H thin film prepared by PECVD in-situ deposition
Gu Zhifeng luoyanggzf@163.com

1
Shan Feng 2
Liu Jia lj2000ppp933@163.com

2
1 https://ror.org/05d80kz58 grid.453074.1 0000 0000 9797 0900 School of Software, Henan University of Science and Technology, Luoyang, 471023 People’s Republic of China
2 https://ror.org/04nraex26 grid.459728.5 0000 0000 9694 8429 Department of Mathematics and Physics, Luoyang Institute of Science and Technology, Luoyang, 471023 People’s Republic of China
16 9 2024
16 9 2024
2024
14 216122 4 2024
9 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Co-doping of phosphorus and boron elements into crystalline silicon quantum dot (c-Si QD) is an effective approach for enhancing the photoluminescence (PL) performance. In this paper, we report on the preparation of hydrogenated silicon nitride (SiNx:H) thin films embedded with phosphorus-boron co-doped c-Si QDs via plasma enhanced chemical vapor deposition route. Mixed dilution including hydrogen (H2) and argon (Ar) is applied in the in-situ deposition process for optimizing the deposition process. The P-B co-doped c-Si QD/SiNx:H thin films exhibit a wide range of PL spectra. The emission is greatly improved especially for the short-wavelength light when compared to the SiOx:H thin film containing P-B co-doped c-Si QDs. The effects of H2/Ar flow ratio on the structural and optical characteristics of thin films are systematically investigated through a series of characterizations. Experimental results show that various properties, such as crystallinity, QD size, optical band gap and doping concentrations, are effectively controlled by tuning H2/Ar flow ratio. Based on the red-shift of QCE-related PL peak, the successful P-B co-doping into Si QDs are verified. Finally, a comprehensive discussion has been made to analyze the influence of H2-Ar mixed dilution on the film growth and impurity doping in detail in this paper.

Keywords

Si quantum dot
Silicon nitride
Co-doping
PECVD
Subject terms

Materials for optics
Quantum dots
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Crystalline silicon quantum dot (c-Si) has been proved to be a promising candidate for novel optoelectronics and photovoltaic devices by virtue of quantum confinement effect (QCE). Similar as other QD materials, the properties of c-Si QD are determined by QCE1–3. Accordingly, by simply controlling QD size, the emission/absorption spectra, radiative lifetimes, and the conductivity are able to be effectively modified for the intended application4–8. As benefit from the peculiarities, c-Si QD show great potential in manufacturing high-performance optoelectrical devices, in particular for photoluminescence (PL). When comparing with other QD materials (e.g. group II−VI QD, group III−V QD, perovskite QD), c-Si QD is competitive for commercial application owing to the nontoxicity and low cost for PL application.

As for c-Si QD, its bandgap energy is higher than that of bulk Si (~ 1.12 eV) because of QCE. It means a poor PL performance in near-infrared region. One approach for solving the problem is introducing impurities into Si QD, and PL energy can be effectively lowered. Nevertheless, single n-type or p-type doping in c-Si QDs leads to the PL quenching which results from the Auger interaction between photo-excited electron–hole pairs and carriers which are supplied by dopant atoms9. With the aim at eliminating Auger effect, researchers adopt co-doping method which enable P and B impurity atoms to be doped in c-Si QD simultaneously, and PL quenching can be inhibited through the charge compensation of n-type and p-type carriers10,11. Additionally, a further shift towards low-energy side for PL peak takes place as induced by the donor—acceptor (D−A) transitions deriving from P-B co-doping12, thus promoting the infrared emission.

Nowadays the most commonly used method for producing P-B co-doped Si QD is co-sputtering. Through this way, P-B co-doped c-Si QDs can be obtained as embedded in silicon oxide (SiOx) matrix11,13,14. With a large bandgap energy over ~ 9.0 eV, SiOx is considered as a suitable matrix material for enhancing QCE. However, the Si=O bonds at the c-Si QD/SiOx interface results in the interface states which locate within band gap. The interface states weaken the emission of short-wavelength light, thus acting as an obstacle for the realization of Si-based full-color PL devices15. In order to extend the PL wavelength region, the influence of Si=O bonds should be eliminated. Researchers attempt to obtain P-B co-doped colloidal Si QDs by etching away SiOx matrix16 or directly synthesis through non-thermal plasma method17. Then the independent P-B co-doped colloidal Si QDs can be further assembled to thin film by spin-coating for device application18.

Although the thin films based on P-B co-doped colloidal Si QDs show better substrate compatibility and emission ability in short wavelength range, current production techniques applied is relative complex and unsuitable for the high-efficiency film fabrication. In comparison, PECVD in-situ deposition is more efficient in preparing impurity-doped Si QD-based thin film under low temperature since there is just one step for fabricating the Si QD-based thin film with the structure of Si QDs/dielectric matrix19,20. Among various matrix materials, silicon nitride has prominent advantages in PL application by virtue of the relatively large band gap energy (~ 5.3 eV). Moreover, for Si QDs embedded in silicon nitride matrix, the PL emission in short wavelength range can be enhanced due to the absence of Si=O bonds. To our knowledge, there remains no reports concerning the investigations on the P-B co-doping of Si QD-based thin films with the structure of c-Si QD/SiNx.

With the aim at realizing wide PL spectrum, we have investigated the properties of P-B co-doped c-Si QD/SiNx:H thin films. PECVD in-situ deposition technic has been carried out for preparing samples in this work. The mixed dilution including H2 and Ar is adopted for promoting crystallization and doping since the effect of H2-Ar mixed dilution has been verified previously21. Various characterizations have been implemented for studying the properties of P-B co-doped c-Si QD/SiNx:H thin films, and the effect of dilution on crystallization and doping process has been illustrated in detail.

Experimental

In this research, P-B co-doped Si QD/SiNx:H thin films were deposited though PECVD route. A conventional capacitively coupled RF (13.56 MHZ) PECVD system has been employed with the electrode area of 144 cm2 and the parallel electrode distance of 2 cm. Before experiments, p-type silicon (100) wafer and quartz glass were clean by acetone, ethanol and deionized water in subsequence for the utilization as substrate. When the background vacuum was evacuated to 5 × 10–7 torr, the dilution gases H2 and Ar, as well as the precursors silane (SiH4) and ammonia (NH3) were introduced into reaction chamber, then followed by the doping source gases phosphine (PH3) and diborane (B2H6). The flow of SiH4, NH3, PH3 and B2H6 were fixed at 5, 5, 0.25 and 0.125 sccm respectively. The total flow of H2-Ar mixed dilution was fixed at 500 sccm while the H2 flow ration RH = [H2/(Ar + H2)] varied from 0.4 to 1. During deposition, the reaction pressure was maintained at 0.5 torr for all the samples, while the substrate temperature and RF power were set to 150 °C and 120 W respectively. The film thicknesses were controlled to be ~ 300 nm relying on the deposition rates. In addition, the thin films for TEM observation were directly deposited on carbon-coated copper grids with the thickness of 30 ~ 60 nm.

For confirming the existence of Si QDs in thin films, room-temperature Raman measurement was performed using an in Via-Reflex Raman spectrometer with a laser at 514 nm for the films on quartz. High-resolution transmission electron micrographs (HRTEM) images were obtained by means of a JEOL-JEM 2100 transmission electron microscope operating at 200 kV. A structural characterization of the films on quartz was carried out utilizing a fourier transform infrared spectrometer (Nicolet Is10) for studying FTIR spectra. An x-ray photoelectron spectroscopy (XPS) ESCALAB 250 was applied for the films on Si wafers to investigate the chemical bonding configurations of samples as well as atomic concentrations of critical elements. A Horiba Jobin–Yvon T64000 system with a 325 nm He-Cd laser as the excitation source was applied for studying the photoluminescence (PL) properties of the films on quartz. The optical band gap of the films on quartz was estimated from the transmittance spectra obtained by a UV-1601 UV–visible spectrophotometer, with the thickness measured by a Dektak 6 M step profiler.

Results

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 RH 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 RH = 0.6, 0.8 and 1. At RH = 0.4, it can be concluded that there is no c-Si component in the thin film due to the absence of sharp peak. When the RH increases from 0.6 to 1, 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 reports22.Fig. 1 Raman spectra of P-B co-doped SiNx:H thin films deposited under different RH (a). For calculating Fc, the spectra are deconvoluted to several independent sub-peaks. The deconvolution spectrum of sample deposited at At RH = 0.8 is shown in (b).

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 At RH = 0.8, and other spectra are shown in Fig. S1. 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−123; (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−124,25; (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 Equation26:Fc=(Ii+Ic)/(β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 as26:β(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 RH is summarized in Table.1. At RH = 0.4, the thin film is completely amorphous. When RH flow achieves at 0.6, c-Si content forms with the Fc of 31.3%. Then as RH increases, the Fc gradually grows. The maximum Fc value of 47.1% is obtained at RH = 1.Table 1 Morphological, compositional and optical properties of P-B co-doped SiNx:H thin films.

RH	Fc (%)	Average QD size (nm)	Optical band gap (eV)	
0.4	N/A	N/A	2.78	
0.6	31.3	4.76	2.24	
0.8	34.4	6.17	1.85	
1	47.1	8.58	1.68	

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 RH = 0.4 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 RH = 0.8, 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 RH grows from 0.6 to 1 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.Fig. 2 TEM micrographs of P-B co-doped SiNx:H thin films deposited at RH = 0.4 (a) and RH = 0.8 (b). Insets of (b) show the typical HRTEM images of Si QDs.

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)227. α 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 RH is shown in Table.1. With the RH increasing from 0.4 to 1, the Eopt value decreases from 2.78 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.

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 RH. 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−128. As RH increases, the intensity of Si–N absorption band shrinks significantly, indicating the decreasing volume fraction of SiNx:H. Similarly, the intensities of Si–H wagging mode, N–H rocking mode and Si–H stretching mode are gradually decreasing with RH, it implies that there are more Si–H, Si–N and N–H bonds formed in the film deposited with larger Ar flow. In contrast, the N–H stretching mode appears no obvious change. 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.Fig. 3 FTIR spectra of P-B co-doped SiNx:H thin films deposited under different RH.

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 RH = 0.8 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–104 eV is determined by five different Si binding states Sin+ (n = 0–4), and n represents the number of nitrogen atoms bonding with one silicon atoms29. 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 RH = 0.4 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 RH increases, the peak gradually shifts towards the low energy side, implying the increment in the proportion of c-Si component.Fig. 4 XPS spectra of P-B co-doped SiNx:H thin films: Si 2p (a), P 2p (b) and B 1 s (c).

Fig. 5 XPS spectra of the sample deposited at RH = 0.8 sccm: Si 2p (a), P 2p (b) and B 1 s (c).

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 RH = 0.8 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 interface30. With the increasing RH, 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 RH 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 1 s 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 RH = 0.4. 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 RH 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 RH = 0.8. 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 Si27. 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 RH. As seen, the Si concentration increases significantly with the increasing RH whereas N concentration exhibits a opposite varying tendency. We calculate the N/Si concentration ratio, and find it decreases from 73.79 to 18.14% when RH increases from 0.4 to 1. 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 RH. It is easily found that both the P and B concentrations gradually decrease with the increasing RH. As the RH increases from 200 to 500 sccm, the calculated P/Si concentration ratio decreases from 2.56 to 0.49%, while the calculated B/Si concentration ratio decreases from 1.80 to 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.Fig. 6 variation of atomic concentrations of Si and N with RH (a), variation of atomic concentrations of P and B with RH (b).

Figure 7 represents the PL spectra of the P-B co-doped SiNx:H thin films deposited with different RH. 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+ level8,31, while P2 at ~ 2.41 eV originates from the radiative recombination at Si-dangling bond (K0) center8,32. Additionally, a P3 peak with varying center in the range of 1.45–2.11 eV is observed for the thin films deposited at RH = 0.6, 0.8 and 1. 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 RH increases from 0.6 to 1. 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.Fig. 7 Photoluminescence spectra of P-B co-doped SiNx:H thin films deposited with different RH.

The comparison of PL performance for the co-doped and un-doped SiNx:H thin films deposited at RH = 0.8 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.Fig. 8 Photoluminescence spectra of co-doped and un-doped samples deposited at RH = 0.8

Discussion

According to the measurement results presented above, P-B co-doped c-Si QD/SiNx:H thin films have been successfully prepared through PECVD in-situ deposition method. The production consists of only one step, and the process temperature is limited as low as 150 °C. As for the co-sputtering method, the formation of c-Si QD depends on the high process temperature (generally > 1000 °C) which is beyond the limit of most substrate materials. On the other hand, impurity atoms can be driven towards QD interface and finally move out of c-Si QD in the case of high temperature. This phenomenon is known as self-purification effect33, and it plays a negative role in impurity doping. Unlike co-sputtering method, plasma plays the key role other than temperature for the in-situ deposition. In this process, the Si-based and N-based precursors are dissociated in plasma to be ionized Si-based and N-based radicals, these ionized radicals gradually agglomerate to fabricate amorphous network in which c-Si QDs are self-assembled at nucleation centers. Meanwhile, the ionized B-based and P-based radicals captured by the defects in the growing Si QDs are rapidly overlaid by the Si-based radicals which reach subsequently. As a result, these B and P atoms are trapped in Si QDs. Since the whole process lasts only a few milliseconds14, there is no chance for impurity atoms to escape, instead they permanently locate inside QDs and form displacement doping. Evidently, self-purification effect become rather weak in this plasma-dominated route which is a kinetic growth process. Impurity doping becomes more efficient due to the little heat created during deposition.

The growth mechanism of Si QD-based thin film prepared by the way of PECVD in-situ deposition has been intensively studied, and the effect of dilution gas on film growth have been clearly clarified34–36. As for this research, H2 and Ar are applied for creating plasma. It is well known that H plasma is able to induce crystallization for Si-based film, while the Ar plasma has the same function. What’s more, Ar are much lower in threshold energy comparing with H2 since there are two kinds of metastable states Ar*33. It means that Ar plasma possesses higher energy when excited at the same power, so it is more efficient in creating nucleation centers by breaking weak Si–Si bonds. Furthermore, dissociation takes place in Ar plasma via direct interaction or the electron-accelerating route37. Consequently, dissociation rates of various kinds of molecules can be significantly improved when Ar flow is increased, enabling the incorporation of more atoms into network. As for our samples, the N/Si concentration ratio is only 18.14% at RH = 1. Considering the same flow of SiH4 and NH3, obviously the dissociation of NH3 is insufficient in the case of H2-riched dilution. As more Ar flow are introduced, the N/Si concentration ratio grow larger, indicating that the dissociation rate of NH3 greatly improves, and it becomes close to that of SiH4 in Ar-riched dilution. As verified by previous research, larger N/Si concentration ratio leads to smaller Si QDs and lower Fc15. This is because that more incorporated N leads to larger volume of SiNx:H which plays a role as the barrier against the growth of c-Si QD by inhibiting the fabrication of crystallographic plane. Our results coincide well with the previous researches, and the N/Si concentration ratio can be easily tuned by varying the flow ratio of H2/Ar.

Similar as the N/Si concentration ratio, both the P/Si and B/Si concentration ratios increase when larger Ar flow is applied, as revealed in Fig. 6b. It is reasonable as the PH3 and B2H6 molecules are dissociated more efficiently in Ar-riched dilution. More P and B atoms are able to diffuse into the thin film and form n- or p-type doping. However, the impurity atoms outside c-Si QDs make no contribution to effective doping. Therefore, B/Si or P/Si concentration ratio should not be used for accurately measuring doping level. In consideration of this point, we have to identify the P-B co-doping effect indirectly through analyzing the PL performance. According to Fig. 8, the realization of P-B co-doping in Si QDs can be verified according to the red-shift of QCE-related P3 peak. Notably, the defect-related P1 and P2 peaks shrink for the P-B co-doped sample, it is possible that the impurity atoms passivate the defects, and correspondingly inhibit the defect-related radiative recombination. Moreover, P3 peak exhibits a significant shrinkage with the decreasing RH as shown in Fig. 7. It implies that when simultaneously doped by P and B, larger Si QDs possess lower PL intensity. Hori et al. has demonstrated that the large co-doped Si QDs are p-type semiconductor since donors and acceptors are not well compensated. As a result, the Fermi level is very close to the HOMO level38. As for our samples, it can be deduced that the Auger recombination in the uncompensated larger Si QDs induces the PL quenching, and this effect will enhance with the increasing QD size.

Since our aim is synthesizing a type of film material which possesses excellent PL properties while is suitable for efficient production. The P-B co-doped c-Si QDs/SiNx:H thin films exhibit a wide range of PL spectra from ~ 1.25 to ~ 3.50 eV. In contrast, the PL energy of P-B co-doped c-Si QDs embedded in SiOx can hardly surpass 2.1 eV because of the Si=O bonds39. It suggests the P–B co-doped c-Si QDs/SiNx:H thin films are more appropriate for the application on light-emitting devices by virtue of the wide range of PL spectra. What’s more, the production of P-B co-doped c-Si QDs/SiNx:H thin films is simply as there is only one step, superior to the multi-step process such as the co-sputtering technic, as well as the methods for obtaining co-doped colloidal Si QDs. Based on all the illustration above, the P–B co-doped Si QDs/SiNx:H thin films should be considered as a promising candidate for novel high-performance PL devices.

Conclusion

In conclusion, we have investigated the fabrication process as well as the properties of P-B co-doped Si QDs/SiNx:H thin films prepared by PECVD in-situ deposition method. H2-Ar mixed dilution has been introduced during deposition. Through various measurements, we have demonstrated the successful co-doping of P and B into the Si QDs embedded in SiNx:H matrix by applying H2-Ar mixed dilution. It is also found that by tuning H2/Ar flow ratio, various properties, such as QD size, optical band gap and PL peak, can be effectively controlled. As for the QD-related PL peak, PL quenching has been observed together with a red-shift of peak position as QD size increases. This has been attributed to the imperfect compensation of P and B atoms in large Si QDs. In contrast with the undoped sample, the shrinkage in the intensity of QD-related PL peak has proved the successful co-doping. All the spectra are in a wide PL energy range from infrared region to ultraviolet region. In summary, P-B co-doped Si QD/SiNx:H thin film can be fabricated by PECVD in-situ deposition method at relatively low substrate temperature. This work will be enlightening for the design of novel QD-based luminescence device.

Supplementary Information

Supplementary Figures.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-72560-3.

Author contributions

Zhifeng Gu wrote the main manuscript text. Feng Shan and Jia Liu reviewed the manuscript.

Data availability

The authors declare that the data supporting the findings of this study are available within the paper, its supplementary information files.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Park NM Choi CJ Seong TY Park SJ Quantum confinement in amorphous silicon quantum dots embedded in silicon nitride Phys. Rev. Lett. 2001 86 7 10.1103/PhysRevLett.86.1355
Park, N. M., Choi, C. J., Seong, T. Y. & Park, S. J. Quantum confinement in amorphous silicon quantum dots embedded in silicon nitride. Phys. Rev. Lett. 86, 7 (2001).10.1103/PhysRevLett.86.1355
2. Terada S Ueda H Ono T Saitow K Orange−red Si quantum dot LEDs from recycled rice husks ACS Sustain. Chem. Eng. 2022 10 1765 1776 10.1021/acssuschemeng.1c04985
Terada, S., Ueda, H., Ono, T. & Saitow, K. Orange−red Si quantum dot LEDs from recycled rice husks. ACS Sustain. Chem. Eng. 10, 1765–1776 (2022).10.1021/acssuschemeng.1c04985
3. Li D Chen J Sun T Zhang Y Xu J Li W Chen K Enhanced subband light emission from Si quantum dots/SiO2 multilayers via phosphorus and boron co-doping Opt. Express. 2022 30 12308 10.1364/OE.453086 35472868
Li, D. et al. Enhanced subband light emission from Si quantum dots/SiO2 multilayers via phosphorus and boron co-doping. Opt. Express. 30, 12308 (2022).35472868 10.1364/OE.453086
4. Masaadeh Q Kaplani E Chao Y Luminescent downshifting silicon quantum dots for performance enhancement of polycrystalline silicon solar cells Electronics 2022 11 2433 10.3390/electronics11152433
Masaadeh, Q., Kaplani, E. & Chao, Y. Luminescent downshifting silicon quantum dots for performance enhancement of polycrystalline silicon solar cells. Electronics 11, 2433 (2022).10.3390/electronics11152433
5. Yu X Ma Z Shen Z Li W Chen K Xu J Xu L 3D NAND flash memory based on double-layer NC-Si floating gate with high density of multilevel storage Nanomaterials 2022 12 2459 10.3390/nano12142459 35889681
Yu, X. et al. 3D NAND flash memory based on double-layer NC-Si floating gate with high density of multilevel storage. Nanomaterials 12, 2459 (2022).35889681 10.3390/nano12142459
6. Tsubata R Gotoh K Matsumi M Wilde M Inoue T Kurokawa Y Fukutani K Usami N Silicon nanocrystals embedded in nanolayered silicon oxide for crystalline silicon solar cells ACS Appl. Nano Mater. 2022 5 1820 1827 10.1021/acsanm.1c03355
Tsubata, R. et al. Silicon nanocrystals embedded in nanolayered silicon oxide for crystalline silicon solar cells. ACS Appl. Nano Mater. 5, 1820–1827 (2022).10.1021/acsanm.1c03355
7. Xu Q Cheong IT Meng L Veinot JGC Wang X Silicon surface passivation for silicon-colloidal quantum dot heterojunction photodetectors ACS. Nano. 2021 15 18429 18436 10.1021/acsnano.1c08002 34757719
Xu, Q., Cheong, I. T., Meng, L., Veinot, J. G. C. & Wang, X. Silicon surface passivation for silicon-colloidal quantum dot heterojunction photodetectors. ACS. Nano. 15, 18429–18436 (2021).34757719 10.1021/acsnano.1c08002
8. Sain B Das D Tunable photoluminescence from nc-Si/a-SiNx: H quantum dot thin films prepared by ICP-CVD Phys. Chem. Chem. Phys. 2013 15 3881 10.1039/c3cp43875a 23407687
Sain, B. & Das, D. Tunable photoluminescence from nc-Si/a-SiNx: H quantum dot thin films prepared by ICP-CVD. Phys. Chem. Chem. Phys. 15, 3881 (2013).23407687 10.1039/c3cp43875a
9. Joo BS Jang S Gu M Jung N Han M Effect of Auger recombination induced by donor and acceptor states on luminescence properties of silicon quantum dots/ SiO2 multilayers J. Alloy. Compd. 2019 801 568 572 10.1016/j.jallcom.2019.06.171
Joo, B. S., Jang, S., Gu, M., Jung, N. & Han, M. Effect of Auger recombination induced by donor and acceptor states on luminescence properties of silicon quantum dots/ SiO2 multilayers. J. Alloy. Compd. 801, 568–572 (2019).10.1016/j.jallcom.2019.06.171
10. Fujii M Toshikiyo K Takase Y Yamaguchi Y Hayashi S Effect of Auger recombination induced by donor and acceptor states on luminescence properties of silicon quantum Dots/SiO2 multilayers J. Appl. Phys. 2003 94 1990 10.1063/1.1590409
Fujii, M., Toshikiyo, K., Takase, Y., Yamaguchi, Y. & Hayashi, S. Effect of Auger recombination induced by donor and acceptor states on luminescence properties of silicon quantum Dots/SiO2 multilayers. J. Appl. Phys. 94, 1990 (2003).10.1063/1.1590409
11. Fujii M Yamaguchi Y Takase Y Ninomiya K Hayashi S Control of photoluminescence properties of Si nanocrystals by simultaneously doping n- and p-type impurities Appl. Phys. Lett. 2004 85 1158 10.1063/1.1779955
Fujii, M., Yamaguchi, Y., Takase, Y., Ninomiya, K. & Hayashi, S. Control of photoluminescence properties of Si nanocrystals by simultaneously doping n- and p-type impurities. Appl. Phys. Lett. 85, 1158 (2004).10.1063/1.1779955
12. Limpens R Fujii M Neale NR Gregorkiewicz T Negligible electronic interaction between photoexcited electron hole pairs and free electrons in phosphorus-boron co-doped silicon nanocrystals J. Phys. Chem. C. 2018 122 6397 6404 10.1021/acs.jpcc.7b12313
Limpens, R., Fujii, M., Neale, N. R. & Gregorkiewicz, T. Negligible electronic interaction between photoexcited electron hole pairs and free electrons in phosphorus-boron co-doped silicon nanocrystals. J. Phys. Chem. C. 122, 6397–6404 (2018).10.1021/acs.jpcc.7b12313
13. Nakamura T Adachi S Fujii M Miura K Yamamoto S Phosphorus and boron codoping of silicon nanocrystals by ion implantation: Photoluminescence properties Phys. Rev. B 2012 85 045441 10.1103/PhysRevB.85.045441
Nakamura, T., Adachi, S., Fujii, M., Miura, K. & Yamamoto, S. Phosphorus and boron codoping of silicon nanocrystals by ion implantation: Photoluminescence properties. Phys. Rev. B 85, 045441 (2012).10.1103/PhysRevB.85.045441
14. Fujii M Sugimoto H Kano S Silicon quantum dot with heavily boron and phosphorus codoped shell Chem. Commun. 2018 54 4375 4389 10.1039/C8CC01612G
Fujii, M., Sugimoto, H. & Kano, S. Silicon quantum dot with heavily boron and phosphorus codoped shell. Chem. Commun. 54, 4375–4389 (2018).10.1039/C8CC01612G
15. Kim TW Cho CH Kim BH Park SJ Quantum confinement effect in crystalline silicon quantum dots in silicon nitride grown using SiH4 and NH3 Appl. Phys. Lett. 2006 88 123102 10.1063/1.2187434
Kim, T. W., Cho, C. H., Kim, B. H. & Park, S. J. Quantum confinement effect in crystalline silicon quantum dots in silicon nitride grown using SiH4 and NH3. Appl. Phys. Lett. 88, 123102 (2006).10.1063/1.2187434
16. Fukuda M Fujii M Sugimoto H Imakita K Hayashi S Surfactant-free solution-dispersible Si nanocrystals surface modification by impurity control Opt. Lett. 2011 36 4026 4028 10.1364/OL.36.004026 22002374
Fukuda, M., Fujii, M., Sugimoto, H., Imakita, K. & Hayashi, S. Surfactant-free solution-dispersible Si nanocrystals surface modification by impurity control. Opt. Lett. 36, 4026–4028 (2011).22002374 10.1364/OL.36.004026
17. Limpens R Pach GF Neale NR Nonthermal Plasma-synthesized phosphorus-boron co-doped Si nanocrystals: A new approach to nontoxic NIR-emitters Chem. Mater. 2019 31 4426 4435 10.1021/acs.chemmater.9b00810
Limpens, R., Pach, G. F. & Neale, N. R. Nonthermal Plasma-synthesized phosphorus-boron co-doped Si nanocrystals: A new approach to nontoxic NIR-emitters. Chem. Mater. 31, 4426–4435 (2019).10.1021/acs.chemmater.9b00810
18. Sugimoto H Fujii M Imakita K Synthesis of boron and phosphorus codoped all-inorganic colloidal silicon nanocrystals from hydrogen silsesquioxane Nanoscale 2014 6 12354 10.1039/C4NR03857F 25192435
Sugimoto, H., Fujii, M. & Imakita, K. Synthesis of boron and phosphorus codoped all-inorganic colloidal silicon nanocrystals from hydrogen silsesquioxane. Nanoscale 6, 12354 (2014).25192435 10.1039/C4NR03857F
19. Das D Samanta S Reverse Meyer-Neldel rule prevailing in the hole transport of B-doped nc-SiOx: H thin films sustaining degeneracy and performing as suitable window of nc-Si solar cells Physica E 2021 128 114615 10.1016/j.physe.2020.114615
Das, D. & Samanta, S. Reverse Meyer-Neldel rule prevailing in the hole transport of B-doped nc-SiOx: H thin films sustaining degeneracy and performing as suitable window of nc-Si solar cells. Physica E 128, 114615 (2021).10.1016/j.physe.2020.114615
20. Das D Samanta S Advanced nanocrystallinity with widened optical gap realized via microstructural control in P-doped silicon oxide thin films used as window layer in nc-Si solar cells Mater. Chem. Phys. 2020 243 122628 10.1016/j.matchemphys.2020.122628
Das, D. & Samanta, S. Advanced nanocrystallinity with widened optical gap realized via microstructural control in P-doped silicon oxide thin films used as window layer in nc-Si solar cells. Mater. Chem. Phys. 243, 122628 (2020).10.1016/j.matchemphys.2020.122628
21. Zhang X Wu A Shi S Qin F Influence of Ar/H2 ratio on the characteristics of boron-doped nc-Si: H films prepared by electron cyclotron resonance plasma-enhanced chemical vapor deposition Surf. Coat. Technol. 2013 228 S412 S415 10.1016/j.surfcoat.2012.05.056
Zhang, X., Wu, A., Shi, S. & Qin, F. Influence of Ar/H2 ratio on the characteristics of boron-doped nc-Si: H films prepared by electron cyclotron resonance plasma-enhanced chemical vapor deposition. Surf. Coat. Technol. 228, S412–S415 (2013).10.1016/j.surfcoat.2012.05.056
22. Cheng Q Tam E Xu S Ostrikov K Si quantum dots embedded in an amorphous SiC matrix: Nanophase control by non-equilibrium plasma hydrogenation Nanoscale 2010 2 594 600 10.1039/b9nr00371a 20644764
Cheng, Q., Tam, E., Xu, S. & Ostrikov, K. Si quantum dots embedded in an amorphous SiC matrix: Nanophase control by non-equilibrium plasma hydrogenation. Nanoscale 2, 594–600 (2010).20644764 10.1039/b9nr00371a
23. Mercaldo LV Esposito EM Veneri PD Fameli G Mirabella S Nicotra G First and second-order Raman scattering in Si nanostructures within silicon nitride Appl. Phys. Lett. 2010 97 153112 10.1063/1.3501133
Mercaldo, L. V. et al. First and second-order Raman scattering in Si nanostructures within silicon nitride. Appl. Phys. Lett. 97, 153112 (2010).10.1063/1.3501133
24. Zhang P Li S Li D Ren L Qin Z Jiang L Xu J Quantum size-dependent luminescence and nonlinear optical properties of silicon quantum dots/SiO2 multilayer Opt. Laser. Technol. 2023 157 108706 10.1016/j.optlastec.2022.108706
Zhang, P. et al. Quantum size-dependent luminescence and nonlinear optical properties of silicon quantum dots/SiO2 multilayer. Opt. Laser. Technol. 157, 108706 (2023).10.1016/j.optlastec.2022.108706
25. Cheng Q Xu S Ostrikov K low-temperature synthesis of nc-Si in high-density, non-equilibrium plasmas: Enabling nanocrystallinity at very low hydrogen dilution J. Mater. Chem. 2009 19 5134 5140 10.1039/b904227j
Cheng, Q., Xu, S. & Ostrikov, K. low-temperature synthesis of nc-Si in high-density, non-equilibrium plasmas: Enabling nanocrystallinity at very low hydrogen dilution. J. Mater. Chem. 19, 5134–5140 (2009).10.1039/b904227j
26. Raha D Das D Controlling the growth of nanocrystalline silicon by tuning negative substrate bias Sol. Energy Mater. Sol. Cells 2011 95 3181 3188 10.1016/j.solmat.2011.06.048
Raha, D. & Das, D. Controlling the growth of nanocrystalline silicon by tuning negative substrate bias. Sol. Energy Mater. Sol. Cells 95, 3181–3188 (2011).10.1016/j.solmat.2011.06.048
27. Park S Cho E Song D Conibeer G Green MA n-Type silicon quantum dots and p-type crystalline silicon heteroface solar cells Sol. Energy Mater. Sol. Cells. 2009 93 684 690 10.1016/j.solmat.2008.09.032
Park, S., Cho, E., Song, D., Conibeer, G. & Green, M. A. n-Type silicon quantum dots and p-type crystalline silicon heteroface solar cells. Sol. Energy Mater. Sol. Cells. 93, 684–690 (2009).10.1016/j.solmat.2008.09.032
28. Kim BH Cho CH Kim TW Park NM Sung GY Photoluminescence of silicon quantum dots in silicon nitride grown by NH3 and SiH4 Appl. Phys. Lett. 2005 86 091908 10.1063/1.1872211
Kim, B. H., Cho, C. H., Kim, T. W., Park, N. M. & Sung, G. Y. Photoluminescence of silicon quantum dots in silicon nitride grown by NH3 and SiH4. Appl. Phys. Lett. 86, 091908 (2005).10.1063/1.1872211
29. Karcher R Ley L Johnson RL Electronic structure of hydrogenated and unhydrogenated amorphous SiNx (0≤x≤1.6): A photoemission study Phys. Rev. B 1984 30 1896 10.1103/PhysRevB.30.1896
Karcher, R., Ley, L. & Johnson, R. L. Electronic structure of hydrogenated and unhydrogenated amorphous SiNx (0≤x≤1.6): A photoemission study. Phys. Rev. B 30, 1896 (1984).10.1103/PhysRevB.30.1896
30. Wu PJ Wang YC Chen IC Influence of phosphorous doping on silicon nanocrystal formation in silicon-rich silicon nitride films J. Phys. D: Appl. Phys. 2013 46 125104 10.1088/0022-3727/46/12/125104
Wu, P. J., Wang, Y. C. & Chen, I. C. Influence of phosphorous doping on silicon nanocrystal formation in silicon-rich silicon nitride films. J. Phys. D: Appl. Phys. 46, 125104 (2013).10.1088/0022-3727/46/12/125104
31. Robertson J Electronic structure of silicon nitride Philos. Mag. B 1991 63 47 10.1080/01418639108224430
Robertson, J. Electronic structure of silicon nitride. Philos. Mag. B 63, 47 (1991).10.1080/01418639108224430
32. Wang M Li D Yuan Z Yang D Que D Photoluminescence of Si-rich silicon nitride: Defect-related states and silicon nanoclusters Appl. Phys. Lett. 2007 90 131903 10.1063/1.2717014
Wang, M., Li, D., Yuan, Z., Yang, D. & Que, D. Photoluminescence of Si-rich silicon nitride: Defect-related states and silicon nanoclusters. Appl. Phys. Lett. 90, 131903 (2007).10.1063/1.2717014
33. Dalpian GM Chelikowsky JR Self-purification in semiconductor nanocrystals Phys. Rev. Lett. 2006 96 226802 10.1103/PhysRevLett.96.226802 16803336
Dalpian, G. M. & Chelikowsky, J. R. Self-purification in semiconductor nanocrystals. Phys. Rev. Lett. 96, 226802 (2006).16803336 10.1103/PhysRevLett.96.226802
34. Chaibi F Jemai R Aguas H Khemakhem H Khirouni K The effects of argon and helium dilution in the growth of nc-Si: H thin films by plasma-enhanced chemical vapor deposition J. Mater. Sci. 2018 53 3672 3681 10.1007/s10853-017-1791-1
Chaibi, F., Jemai, R., Aguas, H., Khemakhem, H. & Khirouni, K. The effects of argon and helium dilution in the growth of nc-Si: H thin films by plasma-enhanced chemical vapor deposition. J. Mater. Sci. 53, 3672–3681 (2018).10.1007/s10853-017-1791-1
35. Das D Raha D Bhattacharya K Evolution of nc-Si network and the control of Its growth by He/H2 plasma assistance in SiH4 at PECVD J. Nanosci. Nanotechnol. 2009 9 5614 5621 10.1166/jnn.2009.1151 19928275
Das, D., Raha, D. & Bhattacharya, K. Evolution of nc-Si network and the control of Its growth by He/H2 plasma assistance in SiH4 at PECVD. J. Nanosci. Nanotechnol. 9, 5614–5621 (2009).19928275 10.1166/jnn.2009.1151
36. Veprek S Iqbal Z Brunner J Schärli M Preparation and properties of amorphous phosphorus nitride prepared in a low-pressure plasma Phil. Mag. B. 1981 43 527 547 10.1080/01418638108222114
Veprek, S., Iqbal, Z., Brunner, J. & Schärli, M. Preparation and properties of amorphous phosphorus nitride prepared in a low-pressure plasma. Phil. Mag. B. 43, 527–547 (1981).10.1080/01418638108222114
37. Kole A Chaudhuri P Nanocrystalline silicon and silicon quantum dots formation within amorphous silicon carbide by plasma enhanced chemical vapor deposition method controlling the Argon dilution of the process gases Thin Solid Films 2012 522 45 49 10.1016/j.tsf.2012.02.078
Kole, A. & Chaudhuri, P. Nanocrystalline silicon and silicon quantum dots formation within amorphous silicon carbide by plasma enhanced chemical vapor deposition method controlling the Argon dilution of the process gases. Thin Solid Films 522, 45–49 (2012).10.1016/j.tsf.2012.02.078
38. Hori Y Kano S Sugimoto H Imakita K Fujii M Size-dependence of acceptor and donor levels of boron and phosphorus codoped colloidal silicon nanocrystals Nano Lett. 2016 16 2615 2620 10.1021/acs.nanolett.6b00225 26998965
Hori, Y., Kano, S., Sugimoto, H., Imakita, K. & Fujii, M. Size-dependence of acceptor and donor levels of boron and phosphorus codoped colloidal silicon nanocrystals. Nano Lett. 16, 2615–2620 (2016).26998965 10.1021/acs.nanolett.6b00225
39. Wolkin MV Jorne J Fauchet PM Allan G Delerue C Electronic states and luminescence in porous silicon quantum dots: The role of oxygen Phys. Rev. Lett. 1999 82 197 10.1103/PhysRevLett.82.197
Wolkin, M. V., Jorne, J., Fauchet, P. M., Allan, G. & Delerue, C. Electronic states and luminescence in porous silicon quantum dots: The role of oxygen. Phys. Rev. Lett. 82, 197 (1999).10.1103/PhysRevLett.82.197
