Blog

Equivalent circuit of a silicon–lithium p–i–n nuclear radiation detector | Scientific Reports

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  13, Article number: 12525 (2023 ) Cite this article Water Jet Disrupter

Equivalent circuit of a silicon–lithium p–i–n nuclear radiation detector | Scientific Reports

Nuclear radiation detectors are indispensable for research in the field of nuclear radiation, X-ray spectroscopy and other areas. Interest in silicon p–i–n detectors of nuclear radiation is increasing today due to the possibility of their operation under normal conditions. In this paper, an equivalent circuit of a silicon–lithium p–i–n nuclear radiation detector is proposed. The proposed circuit is obtained using the classical Shockley equation for silicon semiconductors and the telegraph equations. The parameters of the equivalent circuit were determined using the multiple regression method. As a result of simulation of the model in the MATLAB Simulink graphical development environment, the amplitude-frequency and phase-frequency characteristics of the proposed model were obtained. Using the Monte Carlo method, the alpha-decay of the uranium isotope \({}_{92}{}^{233}\mathrm{U}\) , thorium isotope \({}_{90}{}^{227}\mathrm{Th}\) and americium isotope \({}_{95}{}^{241}\mathrm{Am}\) the alpha-decay spectrum was obtained. Obtained alpha-decay spectra coincides with the experimental data, presented in previous works of other authors.

Felix Pino, Jessica Carolina Delgado, … Sandra Moretto

Hamideh Daneshvar, Kavoos Ghordoei Milan, … Armin Mosayebi

Sergey I. Maximenko, Jim E. Moore, … Phillip P. Jenkins

Semiconductor p–i–n structured detectors are used in many fields of research as precision instruments1, especially in high-energy physics experiments2. The appearance of detectors with a larger detection area generated great interest in them because they significantly improved the efficiency of the detectors and made it possible to register weakly intense charged particles3. However, today, despite the fact that the physical processes in p–i–n diodes and their characteristics have been well studied, scientists are still working on developing large-sized semiconductor detectors based on p–i–n structures4,5,6. Large sized Si(Li) detectors are used in medical imaging, high-energy astrophysics, Compton polarimetry, monitoring of nuclear waste7 The main problems in improving p–i–n detectors of large sizes are related to its development technology8,9 and the development of optimal readout electronics for these detectors10,11. In12,13 authors showed application of silicon p–i–n diodes for spectroscopy. Equivalent circuit of a p–i–n diode was presented and preamplification noise was investigated.

Dementyev et al.14 in their work broadly studied the readout electronics of p–i–n detectors. In their work, authors brought valuable evidence about the pros and cons of p–i–n diodes as X–ray detectors. As an advantage of p–i–n detectors, they emphasize the following characteristics: resistance to the magnetic field; compact size; low operating voltage; inherent stability, and long lie time. As disadvantages of the p–i–n detectors, the authors mentioned the following characteristics: domain energy resolutions of p–i–n detectors are at low energies therefore they need a high gain preamplifier system relatively poor timing resolution and problems related to accepting high counting rate. A number of these problems were solved by some groups of authors, for instance, Muminov et al.15,16 proposed a unique technology for the fabrication of large-sized Si(Li) p–i–n detectors with help of double-sided diffusion and drift of Li ions into monocrystalline silicon. Applying this technology authors could obtain large-sized Si(Li) p–i–n detectors, where they could increase the counting rate of the detector due to its size and did increase its efficiency due to the uniform distribution of Li ions in the i- region. Mostly used technology for increasing the count rate and resolution of detectors is to use various cooling technologies17,18 during detector operation. In order to achieve a high-speed counting rate Gontard et al.19, have designed a high‐bandwidth circuit at the expense of electronic noise and used a prototype of the electronic circuit connected to a detector, aiming to detect single electron events of 200 keV.

Elshennawy and Sunil20 in their recent work described the architecture and device of p–i–n diode detectors for gamma radiation. In their work, they use a model of the simple diode in order to simulate signals from gamma rays. They divided the array of p–i–n diodes into clusters, where each cluster had ten p–i–n diodes. These clusters were used to create the sensitive area of the detector. As a result of their work, the authors show that their diode had a very good discriminating ability, for a 20 keV in particle energy the dynamic differential value for the diode was as high as 12 mV. Authors in21,22 demonstrated evaluation of dielectric properties of semiconductors. In23 authors considered changes of semiconductor’s dielectric properties after high gamma irradiation doses.

The above articles show a high interest in this topic and a large number of studies are devoted to Si(Li) p–i–n detectors of the gamma and X-ray range. Although the scientific interest in this area is quite large9,24, its amount is not sufficient to meet practical purposes.

One of the most important characteristics of any electrical devices is their accurate equivalent circuit. There are various equivalent models have been developed for p–i–n photodiodes, close relative to Si(Li) p–i–n detectors, that allow predicting the behavior of a semiconductor device under different conditions25,26,27. The equivalent circuit of a semiconductor detector was proposed in the works28,29. In these works, a charge separation method is proposed for determining the coordinates of a charged particle hit. In the work30 modeling p–i–n detector of nuclear radiation in SILVACON development environment. A simulation of the manufacture of a p–i–n structure was made and a simulation of the operation of the detector was carried out. The simulation is based on the classic Shockley equation for an ideal diode. In the work31 equivalent circuit of p–i–n structure under reverse bias was developed. In this work Poisson’s continuity equation to determine charge carrier concentration. The paper also presents an equivalent circuit of a p–i–n structure in the form of an RC chain of n links and a study of the dependence of resistance and capacitance on the magnitude of the reverse voltage. Authors of the work32 reviewed application of p–i–n diodes for gamma and X-ray detection. In this paper, the main attention is paid to the scheme of charge collection and preamplification of the detector signal.

In our recent papers15,33, we proposed a new method of obtaining large-sized Si(Li) p–i–n detectors and investigated the physical processes during the formation of the i-region. In order to deeply explain the processes of the newly obtained detector here, in the current work, we proposed modeling and designing a signal formation procedure in these detectors using the classical Shockley equation for silicon semiconductors and a system of telegraph equations. This paper shows the use of the multiple regression method to determine the values of the equivalent circuit elements using the telegraph equation. The resulting model was simulated and the spectra of alpha particles were obtained during the decay of some isotopes.

When designing automatic readout electronics for detectors, it is necessary to take into account the behavior of the system in various operating conditions. For this purpose, this paper shows the simulation of p–i–n detectors of nuclear radiation by the equivalent substitution method. The silicon–lithium nuclear radiation detector is a semiconductor with a p–i–n structure. Works34,35,36 show equivalent circuits of p–i–n diodes, which are the baseline for our research. In33 we showed the distribution of lithium ions in a silicon crystal under the action of a homogeneous electric field during the creation of the detector. Here our task is to simulate the response of the p–i–n structure to external excitation using an equivalent transformation. In this paper, the modeling process can be divided into three stages: modeling the reverse current through a semiconductor detector at the moment of radiation detection, determining the parameters of an equivalent circuit based on a telegraph equation, and simulating the operation of an equivalent circuit. The first stage consists of an approximate representation of the detector in the form of an ideal diode under reverse bias conditions. Then such a diode is described by the Shockley equation, and the generation of charge carriers is described by the continuity equation taking into account recombination and ionization. At the second stage, the found function will be equated to the classical telegraph partial differential equation for current. At the third stage, a simulation of the resulting equivalent circuit will be made.

Semiconductor silicon Si(Li) detectors are used to obtain alpha, beta and gamma radiation spectra. Let's consider a model of a semiconductor diode with different conductivity of p, n and i-regions. When a semiconductor device is connected to reverse bias, a reverse current occurs in it, as shown in Fig. 1.

Electric circuit of p–i–n detector.

Figure 1 shows the electrical circuit of a semiconductor nuclear radiation detector in operating mode. The voltage V is shifted in the opposite direction to expand the sensitive area. The voltage incident on the detector at the reverse bias can be calculated from the expression (1). The reverse current is obtained from Ohm’s law in differential form (2). The current is created by minority charge carriers in the p, n and i-regions

where \({q}_{e}\) —elementary charge, \(x\) —positional coordinate,\(t\) – time \(, {n}_{p}\left(x\right)\) —the electron concentration in the p-region, \({p}_{n}\left(x\right)\) – the hole concentration in the n-region, \({n}_{i}\left(x\right)\) —the intrinsic concentration of charge carriers in the i-\(region, {v}_{n}\) , \({v}_{p}\) , \({v}_{i}\) —diffusion velocities of electrons, holes n, p and i-regions. The velocity of charge carriers in the corresponding areas \({v}_{n}=\frac{{L}_{n}}{{\tau }_{n}}\) , \({v}_{p}=\frac{{L}_{p}}{{\tau }_{p}}\) , \({v}_{i}=\frac{{L}_{i}}{{\tau }_{i}}\) , where \(L_{n} , L_{p} , L_{i}\) —diffusion length of charge carriers in the corresponding region, the diffusion length for the corresponding region is defined as \(L=\) \(\sqrt{D\tau }\) , \({\tau }_{n}\) , \({\tau }_{p}\) , \({\tau }_{i}\) —lifetime of charge carriers in the corresponding region. The distribution of the concentration of charge carriers in the i-region obeys an exponential law, and can be expressed by the following expression (10)28,29.

where \({\mathrm{n}}_{\mathrm{i}0}\) —the concentration of charge carriers in the i—region after drift, \({\mathrm{W}}_{\mathrm{i}}\) —thickness of the i—region. Let’s apply Shockley's formula for Eq. (2). To do this, we divide the Shockley equation by the cross-sectional area of the detector and apply the definition of current density (4).

where \(k\) —Boltzmann constant, \(T\) —absolute temperature. This current density will be established before the detection of radioactive decay particles. At the moment of detection, there is a sharp jump in current for a short period of time. Since the sensitive area of the detector is the i–region, electron–hole pairs are generated at the time of detection. We write down the Poisson continuity equation for electrons and holes at the time of detection, neglecting recombination (5, 6).

where \(v = \mu {E}_{f}\) drift velocity of charge carriers, \(\mu\) —mobility of the corresponding charges, \({E}_{f}\) —electric field. The solution of these partial differential equations will be the following general solutions:

where \({C}_{1}, {C}_{2}\) – integration constants. The generation function can be expressed in terms of the energy of a radioactive particle and the energy required to generate one electron–hole pair as follows (8):

where E – the energy of a radioactive particle, ε – the energy required to create one electron–hole pair, \(\varepsilon =\frac{14}{5}{E}_{g}+{\varphi E}_{R}\) ,34,35. Eg—band gap of semiconductor material, φER—optical phonon losses, V0 = Sx—unit of volume in which electron–hole pairs were generated, \(\delta (t)\) —Dirac delta function, R(x)—recombination function. Generation function has concentration units. Lost charge for recombination can be expressed as follows.

where q0—charge of generated electrons when a nuclear particle falls, \(\lambda = \frac{{t}_{tr}}{\tau }\) 37—relative charge loss, ttr—time of plasma track, τ—lifetime of charge carriers. During time of plasma track high density of charged particles shields an external electric field. When the plasma track time interval ends all charge carriers are dispersed in the i-region. From charge value can be derived number of charged particles and their concentration, which depends on energy of nuclear particle and material of detector. Recombination function R(x), which also has units of concentration, can be expressed taking into account Eq. (9) as follows (10):

where V0 – is initial volume where at the moment of generation all electron hole pairs located before recombination, \({V}_{0}=\pi {{r}_{tr}}^{2}{l}_{tr}\) 36, \({r}_{tr}\) —track radius, \({l}_{tr}\) —length of track of α-particle in detector, S(x)—function of flared base cylinder, x—coordinate.

We consider the initial conditions to be \(x=d\) , where d—coordinates of the i—region where electron hole pair generated, and t = 0. Then the initial concentration at the time of particle detection:

From here, integrating (8) and substituting the initial conditions (11) into the system (7), we obtain:

And finally, for the concentration of charge carriers in the i-region, we obtain:

Then the total current through the detector consists of the current flowing through it before detection and the current through the detector at the time of detection (14).

\({n}_{p}\) = \({p}_{n}\) = \({n}_{i0}\)

\({\tau }_{n}\) = \({\tau }_{p}\) = \({\tau }_{i}\)

\(S=\pi {r}^{2}\) , r = 55 mm

Figure 2 shows the dependence of the current flowing through a semiconductor detector at different energies of α-particles. Figure 3 shows the dependence of the current flowing through a semiconductor detector at different energies of α-particles on detector thickness.

Dependence of current on radioactive particle energy.

Dependence of current density on radioactive particle energy and detector thickness.

Detector resolution FWHM (full width of the distribution at one half of its maximum height) is determined by Fano factor which is the ratio of the standard deviation of the created electron–hole pairs to the average value of electron–hole pairs (15)38.

where F—Fano factor, N—generated electron–hole pairs number, \(\overline{\mathrm{N} }\) —average generated electron–hole pairs number. Theoretical calculations of the Fano factor are given in the works39,40,41. For silicon factor Fano can be calculated using the following Eq. (16)40 :

In the works42,43 Fano factor is equal 0,117–0,118. In the work44 it is shown that Fano factor is decreasing by 5 times from 0.5 to 0.1 during last 40 years.

FWHM is determined from the Eq. (17):38

To simulate p–i–n detectors of nuclear radiation under various conditions, an equivalent electrical circuit was obtained, which is based on the simulation results obtained above using the Shockley model. Let's imagine the p–i–n detector as a two-port network, as shown in Fig. 4.

Equivalent circuit of the p–i–n detector.

Let's make a system of telegraph equations for this circuit (18).

where u—the voltage on the detector at the moment of the fall of the charged particle, i—the current at the moment of the fall of the charged particle, R—the resistance of the detector, L—the equivalent inductance, C—the capacitance of the detector, G—the conductivity. Since in Eq. (12) we obtained the current, then in the future from system (18) we will use the second equation. We bring the second equation from the system (18) to the form (19).

where \(\mathrm{a}=\frac{1}{GR}\) , \(b=\frac{LC}{GR}\) , \(c=\frac{RC+GL}{GR}\) , \(X=\frac{{\partial }^{2}i}{\partial {x}^{2}}\) , \(Y=\frac{{\partial }^{2}i}{\partial {t}^{2}}\) , \(Z=\frac{\partial i}{\partial t}\) , d—an arbitrary constant.

The left part of Eq. (14) is described by Eq. (12) and can be calculated for various parameters of this equation. In the right part there are partial derivatives of the second order of the current in the positional coordinate and in time, as well as a partial derivative of the first order in time. Consequently, having the values of the current depending on time and on the positional coordinate, partial derivatives of the first and second order can be numerically obtained with respect to the corresponding variables. Figure 5 shows the partial derivatives of the first (5a) and second order (5b) in time, as well as the second order in the positional coordinate for the α-particle energy of 4 MeV and 7 MeV.

(a) First order partial derivatives of current in time, (b) second order partial derivatives of current in time and coordinate.

The task of determining arbitrary variables a, b, c and d can be determined using the multiple regression method. At the same time, the coefficient of determination R2 of multiple regression is 0.98124. Substituting the coefficients obtained from the model, we obtain the expression (20).

To return to the replaced variables and find the parameters of the equivalent circuit, we solve the following system (21), which follows from the replacement of variables for (19).

The system consists of three equations, but contains four variables. The capacitance of the p–i–n structure as shown in28 has an order of tens of pF, and depends on the inversely biased voltage. Then the values of R, L, G will depend on the capacity. Solving system of equations we obtain R, L, G for different capacitances of detector. Figure 6 shows the dependencies of the parameters R, L, G respectively of the equivalent circuit on the capacitance of the p–i–n diode.

Parameters of the equivalent p–i–n detector circuit.

The resulting equivalent circuit of a semiconductor p–i–n detector of nuclear radiation is useful in modeling the response of the detector to external influences in the construction of readout electronics.

Figure 7 shows an equivalent detector circuit made in MATLAB Simulink. Figure 8 shows the results of modeling the frequency response of an equivalent circuit. As can be seen from the graph, in the high frequency range up to several tens of gigahertz, the reaction is a linear function. At the same time, a phase shift of -90° is observed in this range.

Modeling frequency response using MATLAB Simulink.

Frequency response of equivalent circuit.

In the resulting equivalent circuit, it is necessary to add readout electronics to simulate the response of the detector at the time of particle detection. Figure 9 shows the equivalent circuit of the detector and readout electronics, made in MATLAB Simulink. The source of the signals is a block of rectangular pulses, followed by a block of a first-order derivative to create a Dirac function. Next, the signal passes through the equivalent detector circuit and the readout electronic unit, which is an integrating circuit and an amplifier. An integrating circuit is necessary to obtain the charge function as an integral of current and voltage as a function of charge. Thus, the voltage amplitude is proportional to the energy of the particles, and the number of particles is equal to the number of pulses detected by the circuit. The number of particles is determined using the pulse counter block.

Equivalent circuit of the detector and readout electronics.

To simulate the detector operation, we use the particle energy distribution during the alpha-decay of the uranium isotope \({}_{92}{}^{233}\mathrm{U}\) , thorium isotope \({}_{90}{}^{227}\mathrm{Th}\) and americium isotope \({}_{95}{}^{241}\mathrm{Am}\) . To apply the Monte Carlo method, we present the energy distribution function as a step function, adding up the probabilities of each event sequentially. We will generate random numbers from 0 to 1 and determine which probability interval this random number belongs to. Having determined the probability interval, we determine the corresponding energy level and, accordingly, the amplitude of the signal generator in the equivalent circuit.

Performing this operation about 107 times and counting the number of random numbers that fell on a particular probability interval, we will plot the dependence of the number of particles on the energy of the particles. Figure 10 shows alpha-decay spectra of a) uranium isotope \({}_{92}{}^{233}\mathrm{U}\) b) thorium isotope \({}_{90}{}^{227}\mathrm{Th}\) and c) americium isotope \({}_{95}{}^{241}\mathrm{Am}\) . Alpha decay spectrum of uranium isotope \({}_{92}{}^{233}\mathrm{U}\) has one peak at energy 4.824 MeV with FWHM 8 keV, other peaks are at 4.783 MeV and 4.729 MeV. As can be seen from the alpha-decay spectrum of thorium, particles with an energy of about 5.98 MeV have the greatest contribution with FWHM = 9 keV, another peak is observed in the energy region of 5.76 MeV. Alpha decay spectrum of americium \({}_{95}{}^{241}\mathrm{Am}\) has three peaks at energies 5.486 MeV, 5.443 MeV and 5.389 MeV with FWHM = 6 keV.

Alpha-decay spectra of (a) uranium isotope \({}_{92}{}^{233}\mathrm{U}\) (b) thorium isotope \({}_{90}{}^{227} \mathrm{Th}\) and (c) americium isotope \({}_{95}{}^{241}\mathrm{Am}\).

In the works28,29 generation of current totally determined by telegraph equation. Presented analytical solution obtained differential equations gives opportunity to judge about accumulated charge, however, the processes of recombination and ionization and the physical properties of the semiconductor structure are not taken into account, which is one of the most important problems in our work. It was reflected in Eqs. (1)–(4) and (8)–(10). The results of the simulation of the detector in the work30, in particular, the current strength of the order of microamperes is consistent with the results obtained in our model in Fig. 3. In this case, the value of the linear current depends on the energy of the charged particle, as shown in Fig. 2.

The simulation results of the obtained model using the Monte Carlo method show satisfactory results. The obtained alpha-decay spectra coincide with the experimental data shown in45,46,47. Theoretical calculations of FWHM closely related with calculation of Fano factor. Equation (16) obtained in the work40. In the work48 Fano factor obtained for silicon is equal 0.07. In our model Fano factor, calculated using Eq. (16) is equal 0.0895.

In our previous works5,15 the capacity of the developed detectors is of the order of tens of pF, the resistance of the detectors is of the order of tens of kΩ. Detector diameter 110 mm, thickness 8–10 mm. In this paper, the conductivity values G were calculated taking into account the approximation using multiple regression, as shown in Fig. 6, and have values of the order of 10–15 Sm, then the resistance of the i-region is of the order of 1015 Ω. The series resistance R and inductance L are also of the high order of 1011. This may be due to the use of the telegraph equation to describe the flow of electric current. The rapid signal decay that is observed when using the Shockley equation results in a rapid voltage drop. The rapid decrease in current in turn leads to high resistance, which is what we are seeing.

As a result of the work, an equivalent circuit of a p–i–n nuclear radiation detector was obtained. The model is based on the Shockley diode equation for a semiconductor, without taking into account the effects associated with the interaction of the crystal lattice of the p–i–n structure, the effects of lithium ions in the i-region on the generation and movement of charge carriers. This approach makes it possible to simplify the model for linearization and the possibility of using the multiple regression method to obtain an equivalent detector circuit by identifying general trends in changes in physical quantities and without detailing the effects occurring in the nodes of the semiconductor crystal lattice. Using experimental data on particle energy distribution during the alpha-decay of the uranium isotope \({}_{92}{}^{233}\mathrm{U}\) , thorium isotope \({}_{90}{}^{227}\mathrm{Th}\) and americium isotope \({}_{95}{}^{241}\mathrm{Am}\) the alpha-decay spectrum was obtained using the Monte Carlo method in relation to the proposed equivalent circuit. Obtained alpha-decay spectra coincides with the experimental data, presented in previous works of other authors.

All data generated or analyzed during this study are included in this published article.

Ramírez‐Jiménez, F. J., Mondragón‐Contreras, L., & Cruz‐Estrada, P. Application of PIN diodes in physics research. In AIP Conference Proceedings, Vol. 857, No. 2 (American Institute of Physics, 2006).

Loupilov, A., Sokolov, A. & Gostilo, V. X-ray Peltier cooled detectors for X-ray fluorescence analysis. Radiat. Phys. Chem. 61(3–6), 463–464 (2001).

Article  ADS  CAS  Google Scholar 

Bosco, G. L. Development and application of portable, hand-held X-ray fluorescence spectrometers. TrAC Trends Anal. Chem. 45, 121–134 (2013).

Perez, K. et al. Fabrication of low-cost, large-area prototype Si (Li) detectors for the GAPS experiment. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 905, 12–21 (2018).

Article  ADS  CAS  Google Scholar 

Muminov, R. A., Radzhapov, S. A. & Saimbetov, A. K. Developing Si (Li) nuclear radiation detectors by pulsed electric field treatment. Tech. Phys. Lett. 35, 768–769 (2009).

Article  ADS  CAS  Google Scholar 

Muminov, R. A., Saymbetov, A. K. & Toshmurodov, Y. K. Special features of formation of high-performance semiconductor detectors based on αSi–Si (Li) heterostructures. Instrum. Exp. Tech. 56, 32–33 (2013).

Protic, D. et al. Large-volume Si (Li) orthogonal-strip detectors for Compton-effect-based instruments. IEEE Trans. Nucl. Sci. 52(6), 3181–3185 (2005).

Article  ADS  CAS  Google Scholar 

Muminov, R. A., Radzhapov, S. A. & Saimbetov, A. K. silicon–lithium telescopic detector in one crystal. At. Energy 106(2), 141 (2009).

Rogers, F. et al. Large-area Si (Li) detectors for X-ray spectrometry and particle tracking in the GAPS experiment. J. Instrum. 14(10), P10009 (2019).

Gao, W. et al. Design and characterization of a low-noise front-end readout ASIC in 0.18-μm CMOS technology for CZT/Si–PIN detectors. IEEE Trans. Nucl. Sci. 65(5), 1203–1211 (2018).

Article  ADS  CAS  Google Scholar 

Qian, Yi. et al. Development of front-end readout electronics for silicon strip detectors. Chin. Phys. C 37(1), 016101 (2013).

Article  ADS  CAS  Google Scholar 

Gramsch, E. et al. Silicon PIN photodetectors in high-resolution nuclear spectroscopy. Nucl. Instrum. Methods Phys. Res. Sect. A 311(3), 529–538 (1992).

Fann, S.-S., Jiang, Y.-L. & Hwang, H.-L. Operating principles and performance of a novel A-Si: H pin-based x-ray detector for medical image applications. IEEE Trans. Electron Devices 50(2), 341–346 (2003).

Article  ADS  CAS  Google Scholar 

Dementyev, D. V. et al. CsI (Tl) photon detector with PIN photodiode readout for a Kμ3 T-violation experiment. Nucl. Instrum. Methods Phys. Res. Sect. A 440(1), 151–171 (2000).

Article  ADS  CAS  Google Scholar 

Muminov, R. A. et al. Double sided diffusion and drift of lithium ions on large volume silicon detector structure. J. Semicond. Technol. Sci. 17(5), 591–596 (2017).

Muminov, R. A. et al. Physical features of double sided diffusion of lithium into silicon for large size detectors. J. Nano Electron. Phys. 11(2), 02031–1 (2019).

Murty, V. R. K. & Devan, K. R. S. On the suitability of Peltier cooled Si–PIN detectors in transmission experiments. Radiat. Phys. Chem. 61(3–6), 495–496 (2001).

Article  ADS  CAS  Google Scholar 

Cesareo, R., Gigante, G. E. & Castellano, A. Thermoelectrically cooled semiconductor detectors for non-destructive analysis of works of art by means of energy dispersive X-ray fluorescence. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 428(1), 171–181 (1999).

Article  ADS  CAS  Google Scholar 

Gontard, L. C. et al. Detecting single-electron events in TEM using low-cost electronics and a silicon strip sensor. Microscopy 63(2), 119–130 (2014).

Article  CAS  PubMed  Google Scholar 

Elshennawy, A., Marianno, C. M. & Khatri, S. P., Architecture and 3D device simulation of a PIN diode-based Gamma radiation detector. In Proceedings of the 23rd ACM International Conference on Great Lakes Symposium on VLSI (2013).

Akay, D., Gökmen, U. & Ocak, S. B. An evaluation of dielectric qualities by using frequency dependence in superbenzene-ring based organic polymer-semiconductors. Mater. Chem. Phys. 245, 122708 (2020).

Akay, D., Gokmen, U. & Ocak, S. B. Ionizing radiation influence on rubrene-based metal polymer semiconductors: direct information of intrinsic electrical properties. JOM 72(6), 2391–2397 (2020).

Article  ADS  CAS  Google Scholar 

Akay, D., Gokmen, U. & Ocak, S. B. Radiation-induced changes on poly(methyl methacrylate) (PMMA)/lead oxide (PbO) composite nanostructure. Phys. Scr. 94(11), 115302 (2019).

Article  ADS  CAS  Google Scholar 

Kozai, M. et al. Developing a mass-production model of large-area Si (Li) detectors with high operating temperatures. Nucl. Instrum. Methods Phys. Res. Sect. A 947, 162695 (2019).

Steinbach, A. H., et al. Equivalent circuit modelling of pin photodiodes for 40 Gb/s receivers. In The 15th Annual Meeting of the IEEE Lasers and Electro-Optics Society, Vol. 2 (IEEE, 2002).

Wang, G. et al. A time-delay equivalent-circuit model of ultrafast pin photodiodes. IEEE Trans. Microw. Theory Tech. 51(4), 1227–1233 (2003).

Malyshev, S. A. & Chizh, A. L. PIN photodiodes for optical control of microwave circuits. IEEE J. Sel. Top. Quantum Electron. 10(4), 679–685 (2004).

Article  ADS  CAS  Google Scholar 

Kalbitzer, S. & Melzer, W. On the charge dividing mechanism in position sensitive detectors. Nucl. Inst. Methods 56(2), 301–304 (1967).

Laegsgaard, E. Position-sensitve semiconductor detectors. Nucl. Inst. Methods 162(1–3), 93–111 (1979).

Article  ADS  CAS  Google Scholar 

Shrivastava, S., & Rabinder, H. Process simulation and analysis for PIN detector. In 2009 IEEE Circuits and Systems International Conference on Testing and Diagnosis (IEEE, 2009).

Senhouse, L. S. Reverse biased PIN diode equivalent circuit parameters at microwave frequencies. IEEE Trans. Electron Devices 3, 314–322 (1966).

Ramírez‐Jiménez, F. J. PIN diode detectors. In: AIP Conference Proceedings, Vol. 1026, No. 1 (American Institute of Physics, 2008).

Saymbetov, A. et al. Physical processes during the formation of silicon–lithium pin structures using double-sided diffusion and drift methods. Materials 14(18), 5174 (2021).

Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

Bryant, A. T. et al. Two-step parameter extraction procedure with formal optimization for physics-based circuit simulator IGBT and pin diode models. IEEE Trans. Power Electron. 21(2), 295–309 (2006).

Zhang, Ao. & Gao, J. Comprehensive analysis of linear and nonlinear equivalent circuit model for GaAs-PIN diode. IEEE Trans. Ind. Electron. 69(11), 11541–11548 (2021).

Mili, S., Aguili, C. L. & Aguili, T. Study of fractal-shaped structures with pin diodes using the multi-scale method combined to the generalized equivalent circuit modeling. Prog. Electromagn. Res. B 27, 213–233 (2011).

Miller, G. L., Gidson, W. M. & Donovan, P. F. Semiconductor particle detectors. Annu. Rev. Nucl. Sci. 12(1), 189–220 (1962).

Article  ADS  CAS  Google Scholar 

Eberhardt, J. E. Fano factor in semiconductor radiation detectors. Diss. UNSW Sydney (1970).

Klein, C. A. Bandgap dependence and related features of radiation ionization energies in semiconductors. J. Appl. Phys. 39(4), 2029–2038 (1968).

Article  ADS  CAS  Google Scholar 

Pell, E. M. Ion drift in an n-p Junction. J. Appl. Phys. 31(2), 291–302 (1960).

Article  ADS  CAS  Google Scholar 

Alkhazov, G. D., Komar, A. P. & Vorob’ev, A. A. Ionization fluctuations and resolution of ionization chambers and semiconductor detectors. Nucl. Instrum. Methods 48(1), 1–12 (1967).

Article  ADS  CAS  Google Scholar 

Mazziotta, M. N. Electron–hole pair creation energy and Fano factor temperature dependence in silicon. Nucl. Instrum. Methods Phys. Res. Sect. A 584(2–3), 436–439 (2008).

Article  ADS  CAS  Google Scholar 

Lowe, B. G. & Sareen, R. A. A measurement of the electron–hole pair creation energy and the Fano factor in silicon for 5.9 keV X-rays and their temperature dependence in the range 80–270 K. Nucl. Instrum. Methods Phys. Res. Sect. A 576(2–3), 367–370 (2007).

Article  ADS  CAS  Google Scholar 

Devanathan, R. et al. Signal variance in gamma-ray detectors—A review. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 565(2), 637–649 (2006).

Article  ADS  CAS  Google Scholar 

Haas, R. et al. Alpha spectrometric characterization of thin 233U sources for 229(m)Th production. Radiochim. Acta 108(12), 923–941 (2020).

Ishhanov, BS, Kapitonov, IM, Yudin, NP Particles and Atomic Nuclei, 672 (M.: LKI Publishing House, 2019).

Kandlakunta , P. & Cao , L. Gamma-ray rejection, or detection, with gadolinium as a converter.Radiation.Prot.Dosim.151(3), 586–590 (2012).

Zulliger, H. R. & Aitken, D. W. Fano factor fact and fallacy. IEEE Trans. Nucl. Sci. 17(3), 187–195 (1970).

Article  ADS  CAS  Google Scholar 

This research has been funded by the Science Committee of the Ministry of Science and Higher Education of the Republic of Kazakhstan (Grant No. AP09058014).

Al-Farabi Kazakh National University, Almaty, Kazakhstan

Ahmet Saymbetov, Zhang Jing, Madiyar Nurgaliyev, Nursultan Japashov, Nurzhigit Kuttybay, Ainur Kapparova, Batyrbek Zholamanov, Sayat Orynbassar & Nursultan Koshkarbay

Physical-Technical Institute, Uzbekistan Academy of Sciences, Tashkent, Uzbekistan

Tashkent Institute of Irrigation and Agricultural Mechanization Engineers, Tashkent, Uzbekistan

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

R.M., A.S. and N.J. contributed significantly to conceptualization and interpretation of results. Y.T., M.N. and N.K. contributed to modeling the response of p–i–n structure on nuclear radiation. Z.J., N.K. and B.Z. contributed to designing equivalent circuit of Si(Li) detector. A.K. and S.O. contributed to modeling p–i–n detector using Monte-Carlo method. All authors read and approved the final version of the manuscript.

The authors declare no competing interests.

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

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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/4.0/.

Saymbetov, A., Muminov, R., Jing, Z. et al. Equivalent circuit of a silicon–lithium p–i–n nuclear radiation detector. Sci Rep 13, 12525 (2023). https://doi.org/10.1038/s41598-023-39710-5

DOI: https://doi.org/10.1038/s41598-023-39710-5

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Equivalent circuit of a silicon–lithium p–i–n nuclear radiation detector | Scientific Reports

Explosive Devices Disruptor Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.