Loading [MathJax]/jax/output/HTML-CSS/jax.js
J. Semicond. > 2018, Volume 39 > Issue 9 > 094010

SEMICONDUCTOR DEVICES

Multivariate rational regression and its application in semiconductor device modeling

Yuxi Hong, Dongsheng Ma and Zuochang Ye

+ Author Affiliations

 Corresponding author: Zuochang Ye, zuochang@tsinghua.edu.cn

DOI: 10.1088/1674-4926/39/9/094010

PDF

Turn off MathJax

Abstract: Physics equation-based semiconductor device modeling is accurate but time and money consuming. The need for studying new material and devices is increasing so that there has to be an efficient and accurate device modeling method. In this paper, two methods based on multivariate rational regression (MRR) for device modeling are proposed. They are single-pole MRR and double-pole MRR. The two MRR methods are proved to be powerful in nonlinear curve fitting and have good numerical stability. Two methods are compared with OLS and LASSO by fitting the SMIC 40 nm MOS-FET I–V characteristic curve and the normalized mean square error of Single-pole MRR is 3.02×108 which is 4 magnitudes less than an ordinary least square. The I–V characteristics of CNT-FET and performance indicators (noise factor, gain, power) of a low noise amplifier are also modeled by using MRR methods. The results show MRR methods are very powerful methods for semiconductor device modeling and have a strong nonlinear curve fitting ability.

Key words: multivariate rational regressionMRRsemiconductor device modelingvector fitting

Accurate device modeling is essential for circuit simulation and design. In 1996, BSIM3 Version 3 (commonly abbreviated as BSIM3v3) was established by SEMATECH as the first industry-wide standard of its kind[1]. It has since been widely used by most semiconductor and IC design companies world-wide for device modeling and CMOS IC design. Though the BSIM model is accurate, it needs a long time to adjust for non-ideal effects. Meanwhile, Moore’s law has nearly come to an end and lots of new devices show up which need study[24]. It is unwise to invest a huge amount of resources to model a new device; we may just test its circuit characteristics rather than business use. Efficient modeling methods should be proposed for these purposes.

A lot of methods have been proposed for the similar usage. Gustanven proposed a lot of black box Macro-modeling methods for huge complex systems such as the high-voltage system and electronic magnetic system[5]. The whole large systems are viewed as a black box with input and output, and statistical regression models are built to describe these systems. By using this method, the complexity of systems is greatly reduced[68]. The same thought can be used in semiconductor device modeling since a semiconductor device can also be viewed as a black box. If an ideal statistical regression model is found, we can save a lot of time and money. For example, the carbon nanotubes field effect transistor (CNT-FET) is a promising candidate for MOS-FET, which generates much less heat and runs as fast as MOS-FET[9]. But CNT-FET shows different I–V characteristics with a different manufacturing process, which is hard for us to build a universal physical model for all of them. A black box statistical modeling is a choice for efficiently modeling different CNT-FET[10, 11].

The keys of statistical modeling are the choice of model and at what extent the model can fit the true data to. Numerous statistical models are proposed for this task. Ordinary least square (OLS) and its regularization method (Ridge and LASSO) are the most commonly used regression approaches[12]. But they are mainly used for linear regression and may not fit for semiconductor device modeling, since the device has a highly nonlinear property. In this paper, we propose two numerical methods for semiconductor device modeling, single-pole denominator numerator fitting (single-pole MRR) and double-pole MRR. They have a great nonlinear curve fitting ability and good numerical stability which are critical in semiconductor device modeling.

A lot of methods have been proposed for the curve fitting task. Suppose we have p dimension input attributes and one dimension output observed data points from the curve, (xi,yi),i=1,...,K,xiRP,yiR,PN+. We want to find a model ˆyi=g(xi) to best fit our original curve. OLS (Ordinary Least Square) has been widely used as a linear regression method, which supposes the input and output has a simple linear relationship. Firstly, using system input attributes x to construct a feature vector X , then, combining these features linearly to get one output ˆy=g(x;w)=XTw,wRP . OLS tries to minimize the least square error of true data and model output to find optimal parameters, this is also called objective function,

w=argminw||yXTw||2. (1)

If XTX is invertible, the optimal solution of coefficient w is w* = (XTX)–1XTy[1]. The Lasso is a shrinkage method for OLS, which adds the L1 norm to the objective function,

w=argminw||yXTw||2+λPj=1wj. (2)

It can shrink some of the parameters to be exactly zeros if positive penalty λ is large enough to achieve higher predict accuracy and interpretation[1]. Though a basis function such as polynomial basis and triangle basis can be used for nonlinear expansion, the OLS and Lasso suffer from its linear nature and behave not so well on a nonlinear curve fitting task.

Compared with the linear models, the rational model, y=N(x)/D(x)=XTwn/XTwd , has nonlinearity in nature, thus is more suitable for fitting a nonlinear curve. The optimal rational model for a system can be obtained through minimizing objective function,

w=argminw||yN(X)D(X)||2. (3)

However, Eq. (3) cannot be solved directly in closed form as there are unknown parameters in the denominator.

Since the 1950s, considerable effort has been devoted to the development of methods for parameters extraction of rational function. Levy, Sanathanan and Koerner, Lawrence and Rogers, and Stahl have presented various techniques by posing linear least squares problems. Pintelon and Guilaume analyze these and several other techniques[1316]. Vector-Fitting(VF), introduced by Gustavsen and Semlyen and using partial fraction basis, has been widely accepted as a robust modeling method for approximating frequency domain responses[1719].

In this work, we proposed single-pole MRR and double-pole MRR methods based on vector fitting, which are of great help in a nonlinear curve fitting task. The power of MRR is shown by numerical examples involving artificially created binary function, SMIC 40 nm NMOS DC characteristic, CNT-FET, and the LNA performance model.

Consider the rational function approximation,

y=N(X)D(X)=w0+XTwn1XTwd, (4)

where X is a feature vector formed by input attributes x (x is a column vector and its elements are input variable, like x1,x2,x3 and so on). We can use a different basis function to map input attributes x to feature vector X . For example, X=[x1,x2,x3,x21,x22,x23,...]T .

We solve Eq. (4) by transforming it into an iterating OLS problem. We knew that equation y=XTw can be solved by objective function equation 1. In rational function approximation Eq. (4), we have object function,

wn,wd=argminwn,wd||yw0+XTwn1XTwd||2. (5)

Suppose in our iteration process, we have a set of {wtd},t=0,1,...w0d=0 and wt1d is known in the ith iteration, where t is the iteration step of the whole process. We multiply a factor ||1XTwtd1XTwt1d||2 to Eq. (5), and let wtd=wt1d+w+d ,

||1XTwtd1XTwt1d||2argminwn,wd||yw0+XTwtn1XTwtd||2=argmin||wn,wd(1XTw+d1XTwt1d)(yw0+XTwtn1XTwtd)||2=argminwn,wd||yXTw+dy1XTwt1dw0+XTwtn1XTwt1d||2=argminwn,wd||yw0+XTwtn+XTyw+d1XTwt1d||2. (6)

Here in the last equation, a new OLS problem is made where the feature vector is [11XTwt1d,XT1XTwt1d,XTy1XTwt1d] and coefficient is [w0,wtn,w+d] . By solving this OLS problem we get w+d and wtd can be solved by wtd=wt1d+w+d . In the (t+1)th iteration, wtd is known, and we continue the above iteration until the convergence of wd . If limtwtd=C , we will have limtw+d=wtdwt1d=0 and the factor will be limt||1XTwtd1XTwt1d||2=1 . We will get the optimal solution wn,wd if wd converges.

Consider the rational function approximation,

y=N(X)D(X)=w01+XTwn11XTwd1+w02+XTwn21XTwd2. (7)

Though we can still multiply a factor ||1XTwtd11XTwt1d1 1XTwtd21XTwt1d2||2 to its objective function, the equation may have too many matrix products such as XTwd1XTwd2 and is hard to solve. The same as single-pole MRR, we transform Eq. (7) by adding two additional elements in the tth iteration,

wt01+XTwtn11XTwt1d1+wt02+XTwtn21XTwt1d2=(1XTw+d11XTwt1d1XTw+d21XTwt1d2)y, (8)
wt01+XTwtn11XTwt1d1+wt02+XTwtn21XTwt1d2+XTyw+d11XTwt1d1+XTyw+d21XTwt1d2=y. (9)

Here wt1d1+w+d1=wtd1 , wt1d2+w+d2=wtd2 . A new OLS problem is made where the feature vector is [11XTwt1d1, XT1XTwt1d1,11XTwt1d2,XT1XTwt1d2,XTy1XTwt1d1,XTy1XTwt1d2] and the coefficient is [wt01,wtn1,wt02,wtn2,w+d1,w+d2] . We continue the above iteration until the convergence of wd1,wd2 . Different from the single-pole MRR, if we set wd1,wd2 to zero as the initial value, the equation will suffer from rank deficiency, the solution is we get the value of wd1,wd2 from a random normal distribution with zero mean and small variance. In fact, if we get the initial value of wd from a random normal distribution rather than zero in single-pole MRR, the convergence value of wd will not change. This means the random initial value of wd will not affect the final convergence, and the same goes for double-pole MRR.

The key step of MRR is accurately solving the over-determined Eq. (7). However, because the condition of such a problem is poor, the solving of the normal equation is of reduced numerical stability and may result in large errors in the solution. Besides, if there are degrees of magnitude difference between y and x, the rank of matrix A might be rank deficient, which also lead to the worst solution. In order to circumvent these cases, normalizing input attributes x and target output y before the parameters extraction procedure is usually of great help.

Semiconductor device I–V characteristic value might be small in number and varies in a wide range, e.g. the Idof NMOS-FET may vary from 1×105to1×1013 . Log transformation is important in this circumstance. However, this will also require high precision of the model since inverse operation, exponentiation, will zoom in the prediction error dramatically.

To illustrate the validity of the proposed method, we firstly consider an artificially created function defined y=sin(2x21+x22)/(2x21+x22),x1,x2[3,3] . We adopt Single-pole MRR with the decic polynomial basis. Figs. 1(a) and 1(b) show the origin function and MRR fitting result. The normalized mean-squared-error (NMSE) is about 7.45 × 10−10. Figs. 1(c) and 1(d) show the step NMSEs of two MRR methods, and that single-pole MRR achieves better accuracy and stability. So in this fitting task, we adopt single-pole MRR as fitting model.

Figure  1.  (Color online) The original function and MRR fitting function. (a) Original artificial function. (b) Single-pole MRR fitting result. (c) Step NMSEs of single-pole MRR. (d) Step NMSEs of double-pole MRR.

DC characteristics of NMOS-FET are used to show the performance of MRR, compared with OLS and LASSO. As BSIM has been widely used in industry for years and, to some extent, it can represent the authentic physical property of NMOS-FET, so we use Cadence and SPICE to get DC simulation data of SMIC 40 nm NMOS-FET with a channel width of 1 μm and channel length of 40 nm. In this case, we choose two independent variables, Vd and Vg , and one dependent variable, Id .

Fig. 2 shows the fitting results of different algorithms for the dataset. We plot the IdVd curve and IdVg curves to show the fitting performance. Table 1 shows the parameters and NMSE of the fitting result. Single-pole MRR with sextic polynomial, double-pole MRR with quartic polynomial and OLS nonic polynomial are used for contrast. The number of parameters in each model is close. Because of the regularization property, we choose the nonic polynomial for LASSO to include as many meaningful features as possible. After 10-fold cross validation, the optimal regularization parameter λ in Eq. (2) is 1 and the number of non-zero coefficient is 36. Normalization and logarithm transformation are both used as preprocessing methods. Two MRR algorithms outperform OLS and LASSO, as shown in Fig. 2. For convenience, we zoom-in the IdVg curve and clip two ranges of Vg from each curve. The detail of the fitting result can then be seen clearly in Fig. 3.

Table  1.  NMSE and parameters number results of different algorithms.
Algorithm Single-pole MRR
(sextic polynomial)
Double-pole MRR
(quartic polynomial)
OLS
(nonic polynomial)
LASSO
(nonic polynomial)
Parameter 55 58 55 36
NMSE 3.0157 × 10–8 2.631 × 10–6 4.357 × 10–4 5.082 × 10–4
DownLoad: CSV  | Show Table
Figure  2.  (Color online) The fitting results of different algorithms. (a) IdVd curve by single-pole MRR. (b) IdVg curve by single-pole MRR. (c) IdVd curve by double-pole MRR. (d) IdVg curve by double-pole MRR. (e) IdVd curve by OLS. (f) IdVg curve by OLS. (g) IdVd curve by LASSO. (h) IdVg curve by LASSO.
Figure  3.  (Color online) Zoom-in pictures of IdVg curves fitting by different algorithms. (a) Zoom-in of IdVg curve by single-pole MRR, Vg ranges from 0.15 to 0.45 V. (b) Zoom-in of IdVg curve by single-pole MRR, Vg ranges from 0.65 to 1.1 V. (c) Zoom-in of IdVg curve by double-pole MRR, Vg ranges from 0.15 to 0.45 V. (d) Zoom-in of IdVg curve by double-pole MRR, Vg ranges from 0.65 to 1.1 V. (e) Zoom-in of IdVg curve by OLS, Vg ranges from 0.15 to 0.45 V. (f) Zoom-in of IdVg curve by single-pole MRR, Vg ranges from 0.65 to 1.1 V. (g) Zoom-in of IdVg curve by LASSO, Vg ranges from 0.15 to 0.45 V. (h) Zoom-in of IdVg curve by LASSO, Vg ranges from 0.65 to 1.1 V.

The Table 1 details the NMSE of four algorithms. As the data is very small, the Normalized Mean Square Error (NMSE) is used, as shown in Eq. 10. Single-pole MRR with sextic polynomial

1mKi=1[timax(t)yimax(t)]2. (10)

NMSE is 3.0157×108 and double-pole MRR with cubic polynomial’s NMSE is 2.631×106 . The result shows MRRs have a much better nonlinear curve fitting ability than OLS and LASSO. We also compare the computation complexity with OLS and LASSO in this paper. Double-pole MRR reach smaller NMSE than OLS and LASSO with less parameters. It is more efficient than the classics. The time complexity of MRR is iteration time plus the time period of a single iteration of MRR. The time of a single iteration of MRR is equal to other classics. The MRR will converge in less than 10 iterations. This will be discussed in Section 4.2.2. In all of the experiments in this paper, the iteration time is no longer than 1 s. The experiment environment is Intel i7-4790 CPU with 32G memory and MATLAB 2017b.

Although single-pole MRR reaches the minimum of NMSE in the results, double-pole MRR uses a lower polynomial (cubic polynomial) and the fewest parameters to reach the second best fitting solution. In fact, it is the best model with polynomial of low degree (1, 2, 3). Double-pole MRR with cubic polynomial’s NMSE is 1.593 × 10−5. This is very helpful in the circumstance where the number of input attributes is large and polynomial of a high degree cannot be achieved due to the limit of CPU memory, e.g. a sextic polynomial of 12 input attributes requires 1237 GB CPU memory.

The convergence of the algorithm is critical, for our parameters extraction method is based on iteration. As has been analyzed, the stability of every least squared problem would be improved after data normalization. Single-pole MRR has very good numerical stability with all polynomial while double-pole MRR shows good numerical stability for a polynomial of low degrees. When the polynomial degree is high, double-pole MRR may suffer from rank deficiency and the numerical stability will decrease. We get the sequential NMSE of two MRRs’ iteration with a quadratic polynomial. The result is plotted in Fig. 4 in log scale. Fig. 4(a) shows step NMSEs of the global optimal fitting results for 40 nm MOSFET. Single-pole MRR adopts a sextic polynomial and double-pole MRR adopts a quartic polynomial. Fig. 4(b) shows step NMSEs of two MRR methods when they adopt the quadratic polynomial. The best NMSE of single-pole MRR and double-pole MRR with the quadratic polynomial are 5.94×102 and 2.851×105 .

In this section, we apply single-pole MRR to CNT-FET DC behavior modeling. Fig. 4 visualizes the measured Id with respect to Vd and Vg for a CNT-FET ( w=3μm and l=15nm ). Vg ranges from −2 to 0 V with a step size 0.15 V and Vd ranges from −2 to 0.02 V with a step size 0.02 V. We choose a septic polynomial to generate feature vector X.

Figure  4.  (Color online) Step NMSEs of two MRR methods in different polynomial degrees. (a) Step NMSEs of MRR methods in the situation of a polynomial of high degree. (b) Step NMSEs of MRR methods in the situation of a polynomial of low degree.

Normalized mean square error (NMSE) of Idis 5.961×104 . Fig. 5 shows the original data surface and fitting data surface generated by single-pole MRR. Figs. 6(a)6(d) shows the fitting result of output characteristics (|Id|–|Vd|) and transfer characteristics (|Id|–|Vg|) respectively. The prediction fits the measured data well. Figs. 6(e) and 6(f) show the step NMSEs of two MRR methods, they all can converge but the single-pole MRR outperforms the double-pole MRR in the global optimal NMSE, so we choose single-pole MRR as our final fitting model. Considering the manufacturing process of CNT-FET is not mature, there is still great uncertainty in devices, so we think this fitting result is acceptable.

Figure  5.  (Color online) The observed CNT-FET data point surface and MRR fitting surface. (a) Original data of CNT-FET. (b) Single-pole MRR fitting result.
Figure  6.  (Color online) |Id|–|Vd|, |Id|–|Vg| curves and step NMSEs by MRR methods. (a) |Id|–|Vd| curve by single-pole MRR. (b) |Id|–|Vg| curve by single-pole MRR. (c) Zoom-in |Id|–|Vg| curve, Vg ranges from 0–0.8 V. (d) Zoom-in |Id|–|Vg| curve, Vg ranges from 1–2 V. (e) Step NMSEs of single-pole MRR. (f) Step NMSEs of double-pole MRR.

In this section we apply MRR methods to model LNA performance characteristics and demonstrate our algorithm in multi-variables (more than 2) regression. Generally, we get a certain performance of a circuit by means of solving the KCL and KVL equations. The solution is accurate but very time-consuming. If we want to do behavioral simulation and optimization for the RF circuit, directly mathematical mapping between design parameters and performance is highly useful.

We chose an LNA which worked at 4GHz and use Cadence to get simulation data for modeling. Figs. 7(a)7(c) show the fitting result of three LNA indicators, NF (noise factor), gain and power by single-pole MRR. Figs. 7(d)7(f) are step NMSEs of single-pole MRR and Figs. 7(g)7(i) are step NMSEs of double-pole MRR. The performance of two MRR methods are close. We sort the sample point in an increased order and record the index sequence, then we plot data points, predictions and errors according to the sorted index sequence in one figure. We can see that the training result is good. Normalized mean square error (NMSE) of NF, gain and power are 0.0031, 0.00710, and 3.5887×107 .

Figure  7.  (Color online) The original function and MRR fitting function. (a) Noise factor fitting by single-pole MRR. (b) Gain fitting by single-pole MRR. (c) Power fitting by single-pole MRR.(d) Step NMSEs of single-pole MRR for noise factor. (e) Step NMSEs of single-pole MRR for gain. (f) Step NMSEs of single-pole MRR for power. (g) Step NMSEs of double-pole MRR for noise factor. (h) Step NMSEs of double-pole MRR for gain. (i) Step NMSEs of double-pole MRR for power.

This paper proposes a family of numerical methods---MRR to approximate an unknown system and extract model parameters. We firstly use single-pole MRR to fit an artificial function and the result is extremely good. Then we compare the performance among single-pole MRR, double-pole MRR, OLS and LASSO on SMIC 40 nm DC characteristics dataset. The results show single-pole MRR has the highest fitting precision and double-pole performs better than single-pole MRR in circumstances of a low degree polynomial. The MRR methods have a more powerful nonlinear curve fitting ability than OLS and LASSO and are proved to be numerically stable. CNT-FET and LNA performance indicators are also modeled, of which the fitting results are good as well. But there are two key points for using MRR methods. Firstly, users have to pay close attention to the numerical stability of MRR methods. We have done one artificial function fitting task and three device model fitting tasks for the convergence analysis of two MRR methods. The results show single-pole MRR has better numerical stability than double-pole MRR. Secondly, the datasets used for the fitting curve should be well-distributed and not be sparse. MRR methods are powerful in fitting a highly-nonlinear function but can also lead to overfitting if the dataset is ill-distributed. Our paper shows the MRR methods are good choices for semiconductor devices statistical modeling as well as other highly-nonlinear curve fitting tasks.



[1]
Sheu B J, Scharfetter D L, Ko P K, et al. BSIM: Berkeley short-channel IGFET model for MOS transistors. IEEE J Solid-State Circuits, 1987, 22(4): 558 doi: 10.1109/JSSC.1987.1052773
[2]
Waldrop M M. The chips are down for Moore’s law. Nature, 2016, 530(7589): 144 doi: 10.1038/530144a
[3]
Kim S H, Song W, Jung M W, et al. Carbon nanotube and graphene hybrid thin film for transparent electrodes and field effect transistors. Adv Mater, 2014, 26(25): 4247 doi: 10.1002/adma.v26.25
[4]
Chen B, Zhang P, Ding L, et al. Highly uniform carbon nanotube field-effect transistors and medium scale integrated circuits. Nano Lett, 2016, 16(8): 5120 doi: 10.1021/acs.nanolett.6b02046
[5]
Grivet-Talocia S, Gustavsen B. Black-box macromodeling and curve fitting, in passive macromodeling: theory and applications. John Wiley & Sons, Inc, 2015
[6]
Bañuelos-Cabral E S, Gutiérrez-Robles J A, Gustavsen B. Rational fitting techniques for the modeling of electric power components and systems using MATLAB environment. InTech, 2017
[7]
Grivet-Talocia S, Gustavsen B. Black-box macromodeling and its EMC applications. IEEE Electromagn Compat Mag, 2016, 5(3): 71
[8]
Bañuelos-Cabral E S, Gutiérrez-Robles J A, Gustavsen B, et al. Enhancing the accuracy of rational function-based models using optimization. Electr Power Syst Res, 2015, 125: 83 doi: 10.1016/j.jpgr.2015.04.001
[9]
Avouris P, Chen Z H, Perebeinos V. Carbon-based electronics. Nat Nanotechnol, 2017, 2(10): 605
[10]
Che Y C, Chen H T, Gui H, et al. Review of carbon nanotube nanoeletronics and macroelectronics. Semicond Sci Technol, 2014, 29(7): 073001 doi: 10.1088/0268-1242/29/7/073001
[11]
Derenskyi V, Gomulya W, Rios J M S, et al. Carbon nanotube network ambipolar field‐effect transistors with 108 on/off ratio. Adv Mater, 2014, 26(34): 5969 doi: 10.1002/adma.201401395
[12]
Friedman J, Hastie T, Tibshirani R. The elements of statistical learning. New York: Springer, 2001
[13]
Levy E C. Complex-curve fitting. IRE Trans Autom Control, 1959, 4(1): 37
[14]
Sanathanan C, Koerner J. Transfer function synthesis as a ratio of two complex polynomials. IEEE Trans Autom Control, 1963, 8(1): 56 doi: 10.1109/TAC.1963.1105517
[15]
Lawrence P J, Rogers G J. Sequential transfer-function synthesis from measured data. Proceedings of the Institution of Electrical Engineers, 1979, 126(1): 104 doi: 10.1049/piee.1979.0020
[16]
Whitfield A H. Transfer function synthesis using frequency response data. Int J Control, 1986, 43(5): 1413 doi: 10.1080/00207178608933548
[17]
Pintelon R, Guillaume P, Rolain Y, et al. Parametric identification of transfer functions in the frequency domain-a survey. IEEE Trans Autom Control, 1994, 39(11): 2245 doi: 10.1109/9.333769
[18]
Gustavsen B, Semlyen A. Rational approximation of frequency domain responses by vector fitting. IEEE Trans Power Deliv, 1999, 14(3): 1052 doi: 10.1109/61.772353
[19]
Grivet-Talocia S, Gustavsen B. The vector fitting algorithm, in passive macromodeling: theory and applications. John Wiley & Sons, Inc, 2015
Fig. 1.  (Color online) The original function and MRR fitting function. (a) Original artificial function. (b) Single-pole MRR fitting result. (c) Step NMSEs of single-pole MRR. (d) Step NMSEs of double-pole MRR.

Fig. 2.  (Color online) The fitting results of different algorithms. (a) IdVd curve by single-pole MRR. (b) IdVg curve by single-pole MRR. (c) IdVd curve by double-pole MRR. (d) IdVg curve by double-pole MRR. (e) IdVd curve by OLS. (f) IdVg curve by OLS. (g) IdVd curve by LASSO. (h) IdVg curve by LASSO.

Fig. 3.  (Color online) Zoom-in pictures of IdVg curves fitting by different algorithms. (a) Zoom-in of IdVg curve by single-pole MRR, Vg ranges from 0.15 to 0.45 V. (b) Zoom-in of IdVg curve by single-pole MRR, Vg ranges from 0.65 to 1.1 V. (c) Zoom-in of IdVg curve by double-pole MRR, Vg ranges from 0.15 to 0.45 V. (d) Zoom-in of IdVg curve by double-pole MRR, Vg ranges from 0.65 to 1.1 V. (e) Zoom-in of IdVg curve by OLS, Vg ranges from 0.15 to 0.45 V. (f) Zoom-in of IdVg curve by single-pole MRR, Vg ranges from 0.65 to 1.1 V. (g) Zoom-in of IdVg curve by LASSO, Vg ranges from 0.15 to 0.45 V. (h) Zoom-in of IdVg curve by LASSO, Vg ranges from 0.65 to 1.1 V.

Fig. 4.  (Color online) Step NMSEs of two MRR methods in different polynomial degrees. (a) Step NMSEs of MRR methods in the situation of a polynomial of high degree. (b) Step NMSEs of MRR methods in the situation of a polynomial of low degree.

Fig. 5.  (Color online) The observed CNT-FET data point surface and MRR fitting surface. (a) Original data of CNT-FET. (b) Single-pole MRR fitting result.

Fig. 6.  (Color online) |Id|–|Vd|, |Id|–|Vg| curves and step NMSEs by MRR methods. (a) |Id|–|Vd| curve by single-pole MRR. (b) |Id|–|Vg| curve by single-pole MRR. (c) Zoom-in |Id|–|Vg| curve, Vg ranges from 0–0.8 V. (d) Zoom-in |Id|–|Vg| curve, Vg ranges from 1–2 V. (e) Step NMSEs of single-pole MRR. (f) Step NMSEs of double-pole MRR.

Fig. 7.  (Color online) The original function and MRR fitting function. (a) Noise factor fitting by single-pole MRR. (b) Gain fitting by single-pole MRR. (c) Power fitting by single-pole MRR.(d) Step NMSEs of single-pole MRR for noise factor. (e) Step NMSEs of single-pole MRR for gain. (f) Step NMSEs of single-pole MRR for power. (g) Step NMSEs of double-pole MRR for noise factor. (h) Step NMSEs of double-pole MRR for gain. (i) Step NMSEs of double-pole MRR for power.

Table 1.   NMSE and parameters number results of different algorithms.

Algorithm Single-pole MRR
(sextic polynomial)
Double-pole MRR
(quartic polynomial)
OLS
(nonic polynomial)
LASSO
(nonic polynomial)
Parameter 55 58 55 36
NMSE 3.0157 × 10–8 2.631 × 10–6 4.357 × 10–4 5.082 × 10–4
DownLoad: CSV
[1]
Sheu B J, Scharfetter D L, Ko P K, et al. BSIM: Berkeley short-channel IGFET model for MOS transistors. IEEE J Solid-State Circuits, 1987, 22(4): 558 doi: 10.1109/JSSC.1987.1052773
[2]
Waldrop M M. The chips are down for Moore’s law. Nature, 2016, 530(7589): 144 doi: 10.1038/530144a
[3]
Kim S H, Song W, Jung M W, et al. Carbon nanotube and graphene hybrid thin film for transparent electrodes and field effect transistors. Adv Mater, 2014, 26(25): 4247 doi: 10.1002/adma.v26.25
[4]
Chen B, Zhang P, Ding L, et al. Highly uniform carbon nanotube field-effect transistors and medium scale integrated circuits. Nano Lett, 2016, 16(8): 5120 doi: 10.1021/acs.nanolett.6b02046
[5]
Grivet-Talocia S, Gustavsen B. Black-box macromodeling and curve fitting, in passive macromodeling: theory and applications. John Wiley & Sons, Inc, 2015
[6]
Bañuelos-Cabral E S, Gutiérrez-Robles J A, Gustavsen B. Rational fitting techniques for the modeling of electric power components and systems using MATLAB environment. InTech, 2017
[7]
Grivet-Talocia S, Gustavsen B. Black-box macromodeling and its EMC applications. IEEE Electromagn Compat Mag, 2016, 5(3): 71
[8]
Bañuelos-Cabral E S, Gutiérrez-Robles J A, Gustavsen B, et al. Enhancing the accuracy of rational function-based models using optimization. Electr Power Syst Res, 2015, 125: 83 doi: 10.1016/j.jpgr.2015.04.001
[9]
Avouris P, Chen Z H, Perebeinos V. Carbon-based electronics. Nat Nanotechnol, 2017, 2(10): 605
[10]
Che Y C, Chen H T, Gui H, et al. Review of carbon nanotube nanoeletronics and macroelectronics. Semicond Sci Technol, 2014, 29(7): 073001 doi: 10.1088/0268-1242/29/7/073001
[11]
Derenskyi V, Gomulya W, Rios J M S, et al. Carbon nanotube network ambipolar field‐effect transistors with 108 on/off ratio. Adv Mater, 2014, 26(34): 5969 doi: 10.1002/adma.201401395
[12]
Friedman J, Hastie T, Tibshirani R. The elements of statistical learning. New York: Springer, 2001
[13]
Levy E C. Complex-curve fitting. IRE Trans Autom Control, 1959, 4(1): 37
[14]
Sanathanan C, Koerner J. Transfer function synthesis as a ratio of two complex polynomials. IEEE Trans Autom Control, 1963, 8(1): 56 doi: 10.1109/TAC.1963.1105517
[15]
Lawrence P J, Rogers G J. Sequential transfer-function synthesis from measured data. Proceedings of the Institution of Electrical Engineers, 1979, 126(1): 104 doi: 10.1049/piee.1979.0020
[16]
Whitfield A H. Transfer function synthesis using frequency response data. Int J Control, 1986, 43(5): 1413 doi: 10.1080/00207178608933548
[17]
Pintelon R, Guillaume P, Rolain Y, et al. Parametric identification of transfer functions in the frequency domain-a survey. IEEE Trans Autom Control, 1994, 39(11): 2245 doi: 10.1109/9.333769
[18]
Gustavsen B, Semlyen A. Rational approximation of frequency domain responses by vector fitting. IEEE Trans Power Deliv, 1999, 14(3): 1052 doi: 10.1109/61.772353
[19]
Grivet-Talocia S, Gustavsen B. The vector fitting algorithm, in passive macromodeling: theory and applications. John Wiley & Sons, Inc, 2015
1

Giant spin injection into semiconductor and THz pulse emission

Zheng Feng

Journal of Semiconductors, 2019, 40(7): 070201. doi: 10.1088/1674-4926/40/7/070201

2

Studies of Water V. Five Phonons in Protonic Semiconductor Lattice Model of Pure Liquid Water

Binbin Jie, Chihtang Sah

Journal of Semiconductors, 2017, 38(7): 071001. doi: 10.1088/1674-4926/38/7/071001

3

Semiconductor steady state defect effective Fermi level and deep level transient spectroscopy depth profiling

Ken K. Chin, Zimeng Cheng

Journal of Semiconductors, 2016, 37(9): 092003. doi: 10.1088/1674-4926/37/9/092003

4

Memristive SRAM cell of seven transistors and one memristor

Patrick W. C. Ho, Haider Abbas F. Almurib, T. Nandha Kumar

Journal of Semiconductors, 2016, 37(10): 104002. doi: 10.1088/1674-4926/37/10/104002

5

Progress in complementary metal-oxide-semiconductor silicon photonics and optoelectronic integrated circuits

Hongda Chen, Zan Zhang, Beiju Huang, Luhong Mao, Zanyun Zhang, et al.

Journal of Semiconductors, 2015, 36(12): 121001. doi: 10.1088/1674-4926/36/12/121001

6

Estimation of the optical loss in bent-waveguide superluminescent diodes by an analytical method

Qi An, Peng Jin, Zhanguo Wang

Journal of Semiconductors, 2015, 36(6): 064008. doi: 10.1088/1674-4926/36/6/064008

7

Large-signal characterizations of DDR IMPATT devices based on group Ⅲ-Ⅴ semiconductors at millimeter-wave and terahertz frequencies

Aritra Acharyya, Aliva Mallik, Debopriya Banerjee, Suman Ganguli, Arindam Das, et al.

Journal of Semiconductors, 2014, 35(8): 084003. doi: 10.1088/1674-4926/35/8/084003

8

Simulating and modeling the breakdown voltage in a semi-insulating GaAs P+N junction diode

A. Resfa, Brahimi.R. Menezla, M. Benchhima

Journal of Semiconductors, 2014, 35(8): 084002. doi: 10.1088/1674-4926/35/8/084002

9

The effects of electron irradiation on the optical properties of the organic semiconductor polypyrrole

J. V. Thombare, M. C. Rath, S. H. Han, V. J. Fulari

Journal of Semiconductors, 2013, 34(9): 093001. doi: 10.1088/1674-4926/34/9/093001

10

AC-electronic and dielectric properties of semiconducting phthalocyanine compounds:a comparative study

Safa'a M. Hraibat, Rushdi M-L. Kitaneh, Mohammad M. Abu-Samreh, Abdelkarim M. Saleh

Journal of Semiconductors, 2013, 34(11): 112001. doi: 10.1088/1674-4926/34/11/112001

11

Microstructural properties of over-doped GaN-based diluted magnetic semiconductors grown by MOCVD

Tao Zhikuo, Zhang Rong, Xiu Xiangqian, Cui Xugao, Li Li, et al.

Journal of Semiconductors, 2012, 33(7): 073002. doi: 10.1088/1674-4926/33/7/073002

12

A theoretical model of the femtosecond laser ablation of semiconductors considering inverse bremsstrahlung absorption

Lin Xiaohui, Zhang Chibin, Ren Weisong, Jiang Shuyun, Ouyang Quanhui, et al.

Journal of Semiconductors, 2012, 33(4): 046002. doi: 10.1088/1674-4926/33/4/046002

13

Electrical characterization of the organic semiconductor Ag/CuPc/Au Schottky diode

Mutabar Shah, M. H. Sayyad, Kh. S. Karimov

Journal of Semiconductors, 2011, 32(4): 044001. doi: 10.1088/1674-4926/32/4/044001

14

Neurocomputing van der Pauw function for the measurement of a semiconductor's resistivity without use of the learning rate of weight vector regulation

Li Hongli, Sun Yicai, Wang Wei, Harry Hutchinson

Journal of Semiconductors, 2011, 32(12): 122002. doi: 10.1088/1674-4926/32/12/122002

15

Modeling and experimental research on a removal mechanism during chemical mechanical polishing at the molecular scale

An Wei, Zhao Yongwu, Wang Yongguang

Journal of Semiconductors, 2010, 31(11): 116005. doi: 10.1088/1674-4926/31/11/116005

16

Comparison of electron transmittances and tunneling currents in an anisotropic TiNx/HfO2/SiO2/p-Si(100) metal–oxide–semiconductor (MOS) capacitor calculated using exponential- and Airy-wavefunction app

Fatimah A. Noor, Mikrajuddin Abdullah, Sukirno, Khairurrijal

Journal of Semiconductors, 2010, 31(12): 124002. doi: 10.1088/1674-4926/31/12/124002

17

Terahertz Semiconductor Quantum Well Devices

Liu H C, Luo H, Ban D, Wachter M, Song C Y, et al.

Chinese Journal of Semiconductors , 2006, 27(4): 627-634.

18

Electronic Structure of Semiconductor Nanocrystals

Li Jingbo, Wang Linwang, Wei Suhuai

Chinese Journal of Semiconductors , 2006, 27(2): 191-196.

19

Temperature Dependence of Vacuum Rabi Splitting in a Single Quantum Dot-Semiconductor Microcavity

Zhu Kadi, Li Waisang

Chinese Journal of Semiconductors , 2006, 27(3): 489-493.

20

Passively Mode Locked Diode-End-Pumped Yb∶YAB Laser with High Reflectivity Type Semiconductor Saturable Absorption Mirror

Wang Yonggang, Ma Xiaoyu, Xue Yinghong, Sun Hong,Zhang Zhigang,and Wang Qingyue

Chinese Journal of Semiconductors , 2005, 26(2): 250-253.

1. Kumar, C.V., Sheela, J.J.J. Investigating the analysis between the logistic regression compared over J 48 decision tree in predicting semiconductor device health. AIP Conference Proceedings, 2024, 3193(1): 020175. doi:10.1063/5.0232947
2. Di Caprio, U., Wu, M., Elmaz, F. et al. Hybrid modelling of a batch separation process. Computers and Chemical Engineering, 2023. doi:10.1016/j.compchemeng.2023.108319
  • Search

    Advanced Search >>

    GET CITATION

    Yuxi Hong, Dongsheng Ma, Zuochang Ye. Multivariate rational regression and its application in semiconductor device modeling[J]. Journal of Semiconductors, 2018, 39(9): 094010. doi: 10.1088/1674-4926/39/9/094010
    Y X Hong, D S Ma, Z C Ye, Multivariate rational regression and its application in semiconductor device modeling[J]. J. Semicond., 2018, 39(9): 094010. doi: 10.1088/1674-4926/39/9/094010.
    shu

    Export: BibTex EndNote

    Article Metrics

    Article views: 4602 Times PDF downloads: 62 Times Cited by: 2 Times

    History

    Received: 07 January 2018 Revised: 12 March 2018 Online: Uncorrected proof: 25 April 2018Accepted Manuscript: 26 April 2018Published: 01 September 2018

    Catalog

      Email This Article

      User name:
      Email:*请输入正确邮箱
      Code:*验证码错误
      Yuxi Hong, Dongsheng Ma, Zuochang Ye. Multivariate rational regression and its application in semiconductor device modeling[J]. Journal of Semiconductors, 2018, 39(9): 094010. doi: 10.1088/1674-4926/39/9/094010 ****Y X Hong, D S Ma, Z C Ye, Multivariate rational regression and its application in semiconductor device modeling[J]. J. Semicond., 2018, 39(9): 094010. doi: 10.1088/1674-4926/39/9/094010.
      Citation:
      Yuxi Hong, Dongsheng Ma, Zuochang Ye. Multivariate rational regression and its application in semiconductor device modeling[J]. Journal of Semiconductors, 2018, 39(9): 094010. doi: 10.1088/1674-4926/39/9/094010 ****
      Y X Hong, D S Ma, Z C Ye, Multivariate rational regression and its application in semiconductor device modeling[J]. J. Semicond., 2018, 39(9): 094010. doi: 10.1088/1674-4926/39/9/094010.

      Multivariate rational regression and its application in semiconductor device modeling

      DOI: 10.1088/1674-4926/39/9/094010
      More Information
      • Corresponding author: zuochang@tsinghua.edu.cn
      • Received Date: 2018-01-07
      • Revised Date: 2018-03-12
      • Published Date: 2018-09-01

      Catalog

        /

        DownLoad:  Full-Size Img  PowerPoint
        Return
        Return