删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

Comparative study of various methods for extraction of multi- quantum wells Schottky diode parameter

本站小编 Free考研考试/2022-01-01




1.
Introduction




The Schottky diode is a microelectronic component based on semiconductor materials. The metallization of this device can give an Ohmic or Schottky contact, according to work function of the used metal. Electrical characterizations of Schottky diode are considered as an important stage before the device final fabrication. Among this electrical characterization, the current voltage (I–V) characteristics are used. Indeed, the extraction of Schottky barrier diodes parameters like ideality factor n, barrier height фb, series resistance Rs and saturation current Is, has attracted the attention of many researchers[1-6].



In the literature, the parameters extraction is investigated by using analytical and heuristics technics. The analytical approaches such as Cheung, Werner, and Nord[1, 7, 8] are only useful in narrow regimes. This limitation put consequently the researchers to develop other techniques and methods of extraction. Heuristics algorithms mimic biological behavior to resolve mathematical optimization problem. This is the reason why many researchers use these algorithms to solve parameters extraction issue. The main advantage of this kind of extraction methods is the global convergence for all the desired search space[9, 10].



Among the proposed works in the literature, the modified ABC algorithm (MABC), vertical optimization, genetic algorithm (GA) and simulated annealing (SA) have been used to extract the SBD parameters[5, 11]. However, these works do not give a detailed comparative study with analytical methods and they have been tested at small search spaces.



In this paper, the main contribution is the elaboration of comparative study for different methods of extraction, tested at large range of temperature (from very low temperature 100 K to room temperature 300 K). More precisely, four technics have been developed such as Cheung[1], PSO, ABC and DE to deal with this concern. Concerning the used heuristic technics, the PSO algorithm originally mimics the social behavior and the movement of organisms, such as a flock of birds or a school of fish. Meanwhile, the ABC algorithm simulates the foraging behavior of bees in nature. While the last algorithm, which is DE is a population-based algorithm developed to solve real optimization problems. This algorithm requires four main steps, such as initialization, mutation, recombination, and selection. For more details about these algorithms, see Refs. [5, 11, 12].




2.
Problem statement





2.1
Mathematical model of a Schottky barrier diode




The mathematical model of a Schottky barrier diode (SBD), which is described by its forward bias current, is given by the following expression[13-15]:









$${{I =; }}{{{I}}_{
m{s}}}left{ {{
m{exp}}left[ {frac{{{{q}}left( {{{V}} - {{I}}{{{R}}_{
m{s}}}}
ight)}}{{{{nKT}}}}}
ight]- 1}
ight}.$$

(1)



In Eq. (1), Is is the saturation current. The following formula is used to compute its value[15]:









$${{{I}}_{
m{s}}} = AA^*{{{T}}^{
m{2}}}{
m{exp}}left( {frac{{- q{{
m{phi}}_{
m{b}}}}}{{{{KT}}}}}
ight).$$

(2)



In these equations, I and V are the diode current and bias voltage, respectively. K and A* denote Boltzmann and Richardson constants, respectively; q is the electronic charge; A defines the Schottky contact area; T is the operating temperature. n denotes the diode ideality factor, Rs is the series resistance; and фb is the Schottky barrier height.




2.2
Parameters extraction strategy




The accuracy of the SBD model (given by Eq. (1)) depends on the accurate value of its parameters (more precisely n, Rs and фb). For this purpose, an efficient strategy of SBD parameters extraction has been proposed in this work.



Here, the SBD parameters extraction is considered as an optimization problem, in which the cost criterion to be optimized (in this paper minimized) is the root mean square error (RMSE) between real experimental and estimated Schottky diode currents, as follows:









$${
m{RMSE}} = sqrt {frac{{{I}}}{{{L}}}sumlimits_{i = 1}^{{L}} {{{left[ {{f_{
m{i}}}left( {{{{I}}_{{
m{meas}}}},{{{V}}_{{
m{meas}}}}, eta }
ight)}
ight]}^2}} } ,$$

(3)



where









$$begin{array}{l}!!!!!! {{f}}({{{I}}_{{
m{meas}}}},{{{V}}_{{
m{meas}}}},eta) = {{{I}}_{{
m{meas}}}}; -qquadqquad;;;;;;quad left( {{{{I}}_{
m{s}}}left{ {{
m{exp}}left[ {dfrac{{{{q}}left( {{{{V}}_{{
m{meas}}}} - {{{I}}_{{
m{meas}}}}{{{R}}_{
m{s}}}}
ight)}}{{{{nkT}}}}}
ight] - 1}
ight}}
ight).end{array}$$

(4)



In Eqs. (3) and (4), $eta $ is a vector that includes the model parameters, and it is given by:








$$ eta = left[ {{{n}},{{{I}}_{
m{s}}},{{{R}}_{
m{s}}}}
ight].$$



L denotes the size of the experimental (I–V) data; Vmeas and Imeas are the measured data gotten from the recorded I–V characteristics.



In this paper, the SBD parameters are accurately extracted by using three different heuristics algorithms and an analytical method. Indeed, the heuristic algorithms will change the value of each parameter, in each cycle and within preselected range, until we obtain the lowest cost function (RMSE). Fig. 1 summarizes the main steps of Schottky diode parameters extraction.






onerror="this.onerror=null;this.src='http://www.jos.ac.cn/fileBDTXB/journal/article/jos/2020/10/PIC/19110022-1.jpg'"
class="figure_img" id="Figure1"/>



Download



Larger image


PowerPoint slide






Figure1.
(Color online) The Schottky diode parameters extraction strategy.





3.
Experiments




In this work, a multi-quantum wells (MQWs) GaAs/Al0.33Ga0.67As Schottky diode, which was fabricated at Nottingham University (School of Physics and Astronomy Applied)[16, 17], has been utilized to test the effectiveness of the developed methods.



The electrical characterization of this diode has been accomplished in the DLTS characterization room of Nottingham University. This room includes the characterization equipment needed to test and characterize electronic components at low/high frequencies and temperatures.



The experimental I–V data of this Schottky diode have been collected and measured by using Agilent 41501B pulse generator and Keithly 428 external I–V meter. The DLTS bench also includes a computer for the DLTS software. The different components of DLTS bench are shown in Fig. 2.






onerror="this.onerror=null;this.src='http://www.jos.ac.cn/fileBDTXB/journal/article/jos/2020/10/PIC/19110022-2.jpg'"
class="figure_img" id="Figure2"/>



Download



Larger image


PowerPoint slide






Figure2.
(Color online) The DLTS characterization bench at Nottingham University.




The experimental (I–V) data have been measured for different temperatures with a sampling rate of 10 ms.




4.
Results and discussions




Fig. 3 shows the experimental forward I–V characteristics of the GaAs/Al0.33Ga0.67As Schottky diode. These characteristics were measured at voltage range (0.2–0.7 V) and for low temperature starting from 100 to 300 K. This choice was taken to observe the mechanism of conduction of our sample (Schottky diode) under low temperature conditions (–173 °C).






onerror="this.onerror=null;this.src='http://www.jos.ac.cn/fileBDTXB/journal/article/jos/2020/10/PIC/19110022-3.jpg'"
class="figure_img" id="Figure3"/>



Download



Larger image


PowerPoint slide






Figure3.
(Color online) The I–V characteristics of n-type Al0.33Ga0.67As Schottky diode for different temperatures.




It can be seen from Fig. 3 that by rising temperature, the current value is also increased. This indicates a usual thermal activation phenomenon[13, 16, 18]. The obtained experimental (I–V) characteristics were used to extract the main Schottky diode parameters (n, Rs, Is, фb). The extraction was done by using different algorithms such as DE, PSO and ABC algorithms, additionally with an analytical method[1].



Finally, a comparative study was achieved to evaluate the accuracy of those methods. The setting parameters of those algorithms are given in Table 1. The variation range of the SBD parameters is given in Table 2.






DEPSOABC
Population size: 40
CR: 0.2
MF: 0.5
Swarm size: 100
Inertia weight: 0.4–0.7
C1, C2: 1.2,1.6
Colony size: 160
Limit: 3 × 160 = 480
Cycle: 5000





Table1.
The setting parameters of DE, PSO and ABC algorithms.



Table options
-->


Download as CSV





DEPSOABC
Population size: 40
CR: 0.2
MF: 0.5
Swarm size: 100
Inertia weight: 0.4–0.7
C1, C2: 1.2,1.6
Colony size: 160
Limit: 3 × 160 = 480
Cycle: 5000








ParameterIs (A)nRs (?)
Variation boundary[10–9, 10–6][0, 20][0, 104]





Table2.
The variation range of SBD parameters.



Table options
-->


Download as CSV





ParameterIs (A)nRs (?)
Variation boundary[10–9, 10–6][0, 20][0, 104]





In this simulation, each algorithm was executed 30 times to compute the Mean and STD metrics. Table 3 and Table 4 summarize the obtained parameters for the different competitor methods. These results are the key elements for the comparative study.






TemparatureParameterDEPSOABCCheung
T = 100 KIs1.069 × 10–71.352 × 10–74.816 × 10–75.00 × 10–8
n9.60210.07613.288.61
Rs1.717 × 1031.633 × 1031.286 × 1032217.09
фb0.1760.1740.1630.20
RMSE1.816 × 10–72.715 × 10–72.131 × 10–61.064 × 10–5
Mean1.816 × 10–71.086 × 10–72.693 × 10–6
STD8.050 × 10–221.487 × 10–74.180 × 10–7
T = 120 KIs1.105 × 10–71.083 × 10–71.784 × 10–75.00 × 10–8
n7.9077.8758.6647.12
Rs1.630 × 1031.637 × 1031.520 × 1032.114 × 103
фb0.2140.2150.200.24
RMSE1.679 × 10–71.690 × 10–71.498 × 10–61.305 × 10–5
Mean1.679 × 10–76.763 × 10–82.635 × 10–6
STD2.813 × 10–229.261 × 10–84.234 × 10–7
T = 160 KIs1.437 × 10–71.402 × 10–79.953 × 10–76.00 × 10–8
n5.9625.9379.9355.27
Rs1.483 × 1031.482 × 1034.58 × 1021.875 × 103
фb0.2900.2910.260.33
RMSE4.011 × 10–74.116 × 10–72.601 × 10–61.394 × 10–5
Mean4.011 × 10–71.646 × 10–72.831 × 10–6
STD1.123 × 10–212.254 × 10–74.306 × 10–7





Table3.
The obtained results for the temperature range (100–160 K).



Table options
-->


Download as CSV





TemparatureParameterDEPSOABCCheung
T = 100 KIs1.069 × 10–71.352 × 10–74.816 × 10–75.00 × 10–8
n9.60210.07613.288.61
Rs1.717 × 1031.633 × 1031.286 × 1032217.09
фb0.1760.1740.1630.20
RMSE1.816 × 10–72.715 × 10–72.131 × 10–61.064 × 10–5
Mean1.816 × 10–71.086 × 10–72.693 × 10–6
STD8.050 × 10–221.487 × 10–74.180 × 10–7
T = 120 KIs1.105 × 10–71.083 × 10–71.784 × 10–75.00 × 10–8
n7.9077.8758.6647.12
Rs1.630 × 1031.637 × 1031.520 × 1032.114 × 103
фb0.2140.2150.200.24
RMSE1.679 × 10–71.690 × 10–71.498 × 10–61.305 × 10–5
Mean1.679 × 10–76.763 × 10–82.635 × 10–6
STD2.813 × 10–229.261 × 10–84.234 × 10–7
T = 160 KIs1.437 × 10–71.402 × 10–79.953 × 10–76.00 × 10–8
n5.9625.9379.9355.27
Rs1.483 × 1031.482 × 1034.58 × 1021.875 × 103
фb0.2900.2910.260.33
RMSE4.011 × 10–74.116 × 10–72.601 × 10–61.394 × 10–5
Mean4.011 × 10–71.646 × 10–72.831 × 10–6
STD1.123 × 10–212.254 × 10–74.306 × 10–7








TemparatureParameterDEPSOABCCheung
T = 200 KIs2.347 × 10–72.304 × 10–75.560 × 10–77.00 × 10–8
n5.0645.0516.2044.19
Rs1.282 × 1031.280 × 1031.017 × 1031.758 × 103
фb0.3620.3630.3470.42
RMSE7.96 × 10–78.005 × 10–71.441 × 10–61.829 × 10–5
Mean7.968 × 10–73.202 × 10–72.970 × 10–6
STD5.989 × 10–224.384 × 10–78.012 × 10–7
T = 220 KIs9.245 × 10–79.165 × 10–71 × 10–68.00 × 10–8
n6.6236.6046.7993.79
Rs4.871 × 1024.92 × 1024.405 × 1021.736 × 103
фb0.3760.3760.3750.46
RMSE2.119 × 10–62.149 × 10–62.211 × 10–62.315 × 10–5
Mean2.119 × 10–68.596 × 10–73.134 × 10–6
STD9.470 × 10–221.177 × 10–62.519 × 10–7
T = 240 KIs1.212 × 10-71.364 × 10–78.229 × 10–71.00 × 10–7
n3.3723.4866.0783.46
Rs1.799 × 1031.622 × 103100.231.693 × 103
фb0.4560.4540.4160.51
RMSE1.768 × 10–85.763 × 10–84.810 × 10–64.147 × 10–5
Mean1.768 × 10–82.307 × 10–86.223 × 10–6
STD1.544 × 10–223.159 × 10–86.973 × 10–7
T = 300 KIs2.063 × 10–74.208 × 10–78.450 × 10–81.00 × 10–7
n2.7623.4682.8742.82
Rs1.560 × 1036.496 × 10217.6761.480 × 103
фb0.5680.5490.5910.63
RMSE2.539 × 10–83.959 × 10–79.856 × 10–61.971 × 10–4
Mean2.539 × 10–81.583 × 10–71.329 × 10–5
STD3.766 × 10–222.168 × 10–71.170 × 10–5





Table4.
The obtained results for the temperature range (200–300 K).



Table options
-->


Download as CSV





TemparatureParameterDEPSOABCCheung
T = 200 KIs2.347 × 10–72.304 × 10–75.560 × 10–77.00 × 10–8
n5.0645.0516.2044.19
Rs1.282 × 1031.280 × 1031.017 × 1031.758 × 103
фb0.3620.3630.3470.42
RMSE7.96 × 10–78.005 × 10–71.441 × 10–61.829 × 10–5
Mean7.968 × 10–73.202 × 10–72.970 × 10–6
STD5.989 × 10–224.384 × 10–78.012 × 10–7
T = 220 KIs9.245 × 10–79.165 × 10–71 × 10–68.00 × 10–8
n6.6236.6046.7993.79
Rs4.871 × 1024.92 × 1024.405 × 1021.736 × 103
фb0.3760.3760.3750.46
RMSE2.119 × 10–62.149 × 10–62.211 × 10–62.315 × 10–5
Mean2.119 × 10–68.596 × 10–73.134 × 10–6
STD9.470 × 10–221.177 × 10–62.519 × 10–7
T = 240 KIs1.212 × 10-71.364 × 10–78.229 × 10–71.00 × 10–7
n3.3723.4866.0783.46
Rs1.799 × 1031.622 × 103100.231.693 × 103
фb0.4560.4540.4160.51
RMSE1.768 × 10–85.763 × 10–84.810 × 10–64.147 × 10–5
Mean1.768 × 10–82.307 × 10–86.223 × 10–6
STD1.544 × 10–223.159 × 10–86.973 × 10–7
T = 300 KIs2.063 × 10–74.208 × 10–78.450 × 10–81.00 × 10–7
n2.7623.4682.8742.82
Rs1.560 × 1036.496 × 10217.6761.480 × 103
фb0.5680.5490.5910.63
RMSE2.539 × 10–83.959 × 10–79.856 × 10–61.971 × 10–4
Mean2.539 × 10–81.583 × 10–71.329 × 10–5
STD3.766 × 10–222.168 × 10–71.170 × 10–5





From the results of Table 3 and Table 4, the main observation is the large difference of the RMSE values between the heuristics algorithms and the analytical method. So, the analytical method is not a good candidate to solve this problem. Fig. 4 and Fig. 5 show the weakness of Cheung method to extract the SBD parameters against DE algorithm, at 100 and 300 K, respectively.






onerror="this.onerror=null;this.src='http://www.jos.ac.cn/fileBDTXB/journal/article/jos/2020/10/PIC/19110022-4.jpg'"
class="figure_img" id="Figure4"/>



Download



Larger image


PowerPoint slide






Figure4.
(Color online) Experimental and simulated I–V characteristics from 0 to 0.7 V at 100 K.






onerror="this.onerror=null;this.src='http://www.jos.ac.cn/fileBDTXB/journal/article/jos/2020/10/PIC/19110022-5.jpg'"
class="figure_img" id="Figure5"/>



Download



Larger image


PowerPoint slide






Figure5.
(Color online) Experimental and simulated I–V characteristics from 0 to 0.7 V at 300 K.




Our analysis of the results of the algorithms shows that all of these methods are adequate solutions for this problematic. Moreover, the DE algorithm gives accurate parameters (n, Rs, Is, фb) along the entire temperature range, compared with the others. This accuracy is due to the ability of the DE to converge to the lowest RMSE.



Additionally, the robustness of DE algorithm is verified by comparing the STD metric for the three heuristic algorithms. Indeed, STD values of DE are in scale of 10–22, while are in 10–7 and 10–8 for the others.



Fig. 6 shows the convergence rate of DE, PSO and ABC algorithms at T = 100 K. This figure indicates that ABC algorithm presents a rapid convergence till cycle #800. Up to this cycle, ABC algorithm shows slow convergence. The optimal solution was achieved at iteration #4700. In contrast to ABC algorithm, the others show an apparent rapid convergence speed. Where, it significantly increases in case of DE algorithm compared with PSO and ABC.






onerror="this.onerror=null;this.src='http://www.jos.ac.cn/fileBDTXB/journal/article/jos/2020/10/PIC/19110022-6.jpg'"
class="figure_img" id="Figure6"/>



Download



Larger image


PowerPoint slide






Figure6.
(Color online) Convergence characteristics of DE, PSO and ABC algorithms for T = 100 K.




The variation of RMSE at each running cycle for DE, PSO and ABC algorithm for T = 100 K is shown in Fig. 7. By analyzing the results in this figure, we can see that by increasing the number of running cycle (30) for the different algorithms, the DE algorithm shows the best robustness. In fact, its RMSE value does not change at almost all running cycles compared with PSO and ABC. Consequently, the SBD parameters do not change.






onerror="this.onerror=null;this.src='http://www.jos.ac.cn/fileBDTXB/journal/article/jos/2020/10/PIC/19110022-7.jpg'"
class="figure_img" id="Figure7"/>



Download



Larger image


PowerPoint slide






Figure7.
(Color online) Variation of RMSE at each running cycle for DE, PSO and ABC algorithm for T = 100 K.





5.
Conclusion




In this paper, a comparative study was elaborated for different extraction methods. This comparison is divided on two sections; first, the performance of heuristics algorithms was compared with those of analytical method (Cheung method); and second, the efficiency of DE, PSO and ABC algorithms was investigated.



The obtained results show that the performance metrics, convergence rate speed and the accuracy of the parameters determined by the DE algorithm are better than ABC, PSO and Cheung method. In addition, DE is more robust than the other algorithms. As a result, DE is an effective, robust and fast approach for the parameter extraction of the Schottky diode.



相关话题/Comparative study various