Advanced
A Compact Quantum Model for Cylindrical Surrounding Gate MOSFETs using High-k Dielectrics
A Compact Quantum Model for Cylindrical Surrounding Gate MOSFETs using High-k Dielectrics
Journal of Electrical Engineering and Technology. 2014. Mar, 9(2): 649-654
Copyright © 2014, The Korean Institute of Electrical Engineers
  • Received : September 27, 2013
  • Accepted : October 31, 2013
  • Published : March 01, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
P. Vimala
Corresponding Author: Dept. of Electronics and Communication Engineering, Thiagarajar College of Engineering, Madurai, India. (vimalapalanichamy@tce.edu)
N. B. Balamurugan
Dept. of Electronics and Communication Engineering, Thiagarajar College of Engineering, Madurai, India. (nbbalamurugan@tce.edu)

Abstract
In this paper, an analytical model for Surrounding Gate (SG) metal-oxide- semiconductor field effect transistors (MOSFETs) considering quantum effects is presented. To achieve this goal, we have used variational approach for solving the Poission and Schrodinger equations. This model is developed to provide an analytical expression for inversion charge distribution function for all regions of device operation. This expression is used to calculate the other important parameters like inversion charge density, threshold voltage, drain current and gate capacitance. The calculated expressions for the above parameters are simple and accurate. This paper also focuses on the gate tunneling issue associated with high dielectric constant. The validity of this model was checked for the devices with different dimensions and bias voltages. The calculated results are compared with the simulation results and they show good agreement.
Keywords
1. Introduction
Surrounding gate (SG) MOSFETs have been proposed to extend the scaling limit of conventional MOSFETs to 10nm gate length and beyond, owing their excellent control of the short-channel effects and high current driving ability [1 - 5] . Several classical models have been developed for SG MOSFETs. A compact model should include the quantum effects too.
One the quantum effect is known as energy quantization. The effect of quantization of electronic energy becomes significant for 10 nm and thus it is extremely important to consider quantum effects in MOSFET models. In nanoscale regime, the SG MOSFETs show a different behaviour in comparison to their bulk MOSFETs counter parts is the inversion charge because it is strongly confined in the channel of these devices. As consequence, one of the most interesting parameters that can be modeled is the average inversion layer penetration (defined as the inversion layer centroid [6 - 8] ), which can be used to quantify the influence of quantum effects on the inversion charge of the MOSFETs [7] as well as in the channel capacitance. In this respect, the accurate calculation of these parameters is essential to develop good compact models to be used in circuit simulators. The second effect is that the direct tunneling of carriers between the gate electrode and silicon substrate known as gate tunneling. A solution to this gate leakage problem is the replacement of the silicon di-oxide with a material of higher dielectric constant (a concept abbreviated as high-k) [9] , [10] .The high dielectric constant means that the insulator can be thicker for a given gate capacitance, while presumably allowing exponentially less tunneling current than the oxide. These two phenomena should be incorporated into the classical model to make it complete.
Some models have been established for QM effects of Si Surrounding gate MOSFETs. The exact modeling of SG MOSFETs, which determines the charge distribution and current in the device is a complicated task. The reason is that the complete solution of the problem has to be obtained by solving the coupled Poisson and Schrodinger equation. Yuan [11] developed an analytical model for threshold voltage shift in the surrounding gate MOSFETs considering quantum confinement. Yu-Sheng [12] described the quantum effects in surrounding gate MOSFET and Moreno [13] has studied the quantum effects in the square GAA MOSFETs. However, extensive studies are not done to analyse the quantum effects in surrounding gate MOSFETs. Mathematical complexities of the above models are impractical for compact modeling. Roldan [14] have developed a model for centroid and inversion charge of Surrounding Gate MOSFETs including quantum effects. The wave functions as a function of different device parameters were obtained from the numerical simulator used by them. A good analytical model which simplifies the above problems is still to be developed. Our model is based on the analytical solution of Poisson’s-Schrodinger equation using variational approach. This approach involves guessing a reasonable, parametric form for a trial wave function ( ψ 0 ( r )).
Therefore, the purpose of this paper is to derive an accurate analytical quantum model for cylindrical surrounding gate MOSFETs .The rest of the paper is organized as follows: Section 2 presents the development of quantum analytical model for SG MOSFET by solving the coupled Poisson and Schrodinger equations. Section 3 gives the influence of quantization on device terminal parameters namely centroid, inversion charge density and I-V, C-V characteristics are extracted as closed form functions of radius, inversion charge density and bias voltage for different high-k dielectrics. Model verification and descriptions of main results are presented in section 4. Finally the conclusions are given in section 5.
2. Quantum Modeling
The structure and schematic and cross section view of the SG MOSFET is shown in Fig. 1 with radius of semiconductor cylinder R , dielectric gate insulator thickness thigh−k . For this structure, we solved the Poisson and Schrodinger equations (P-S) using variational approach [15] to obtain the quantum effects.
Modeling the quantum effects with a very high accuracy involves solving the coupled P-S self-consistently. The variational technique represents a completely different way of getting approximate energies and wave functions for quantum mechanical systems. Therefore, to model quantum effects, the method used in this paper refers to the variational approach which involves an analytical approximation of ground state energy based itself initially on the choice of wave function. Moreover, in order to simplify the different equations, we can focus exclusively on the lowest-energy subband wavefunction. By using several trial wave functions, we found the optimum wave function for SG MOSFETs as given below:
PPT Slide
Lager Image
where a 0 and b 0 are the unknown parameters. Normalization of (1) meaning gives an expression of a 0 in terms of b 0 .
PPT Slide
Lager Image
Schematic and Cross section of the silicon SG MOSFETs
Normalization of (1), which means
PPT Slide
Lager Image
gives,
PPT Slide
Lager Image
Considering the eigenfunction of (1), we use the variational approach to solve the Poisson equation as follows,
PPT Slide
Lager Image
where q is the electronic charge, n ( r ) is the electron density, ε si is the silicon permittivity and ϕ ( r ) is the electric potential of silicon region. Solving (3) analytically, we obtained the electric potential in terms of b 0 .
The 1D Schrodinger equation is written as,
PPT Slide
Lager Image
where, h is the Plank’s constant, me = 0.916 m 0 is the effective mass of electrons in lower energy level with (100) orientation, E 0 is the lowest subband energy. The value of kinetic and potential energies in the lowest subband is given by,
PPT Slide
Lager Image
From (5), the lowest subband total energy is written as,
PPT Slide
Lager Image
In the quantum variational approach, b 0 is evaluated by minimizing E 0 , i.e., dE 0 / db 0 = 0. A good approximation for b 0 is derived as,
PPT Slide
Lager Image
Note that the variational parameter ( b 0 ) is dependent on the inversion charge density ( Qinv ) and silicon radius R .
3. Compact Surrounding Gate MOSFET Model
For the compact quantum modeling purpose the parameters like inversion charge centroid, threshold voltage, inversion charge density, gate capacitance and drain current are extracted from the Poisson-Schrodinger equations. In this section, we develop an analytical model for the above parameters of SG MOSFETs which includes the quantum effects.
PPT Slide
Lager Image
In the classical case, the density maximum is always located directly at the interface. But in the quantum mechanical case, the electrons are distributed more toward the center of the silicon body. Based on [14] , the expression of inversion-charge density is given in (8).
The total capacitance is given by,
PPT Slide
Lager Image
where, Cox high-k and Csemi are calculated as,
PPT Slide
Lager Image
PPT Slide
Lager Image
where εhigh-k is the permittivity of high-k dielectric material. xI is the inversion charge centroid is calculated using wavefunction given in Eq. (1).
The inversion layer centroid is obtained by,
PPT Slide
Lager Image
The final expression for SG centroid is given in Eq. (13). The centroid model is validated for various values of R and Ninv .
PPT Slide
Lager Image
The drain current in SG MOSFET is obtained with the assumption that the mobility is independent of the position in the channel. It is expressed as,
PPT Slide
Lager Image
To calculate the drain current we have to solve (8) for V = 0 at the source end, Qinv = Qs and V = Vd at the drain end, Qinv = Qd . Based on [16] , the drain cuttent expression is obtained in terms of inversion charge density as given in Eq. (15).
PPT Slide
Lager Image
where μ is the mobility constant.
4. Results and Discussions
In this section, the accuracy of the proposed model is examined with the simulation results. The accuracy of the proposed model is verified using the commercially available TCAD Sentaurus device simulator. The set of parameters used for simulation are shown in Table-1 . In this section, we will also review the influence of various high-k gate dielectrics in cylindrical surrounding gate MOSFETs. The analysis was performed for a wide range of proposed gate dielectric high-k values such as 3.9, 7.5, 10, 15, 25, and 40. Table 2 shows the different gate dielectric materials and their corresponding high-k constant values.
Nomenclature of symbols
PPT Slide
Lager Image
Nomenclature of symbols
Different gate dielectric materials and their high-k constant values
PPT Slide
Lager Image
Different gate dielectric materials and their high-k constant values
The threshold voltage is examined in Fig. 3 with different silicon radius at different high-k dielectric constant. The threshold voltage value is decreases with the high-k dielectric constant increases. This reduction does not correspond to the decreased limit of inaccuracy of the classical case. But demonstrate the lowering of threshold voltage values in both classical and quantum approach and thus, resulting in small difference in threshold voltage for a given silicon radius value. The reduced threshold voltage promises a choice of low power consuming device together with low switching action, thereby improving the characteristics appreciably for fast switching applications.
Fig. 2 shows the inversion charge density variation with gate voltage at different high-k configurations together with different parameters. The trend of inversion charge density with gate voltage is in accordance with device physics. The replacement of high-k material shows an improved device performance with increased inversion charge density as compared with SiO 2 dielectric. This can be explained as; the high-k device offers high surface potential in comparison to SiO 2 device. This high potential results in inversion of carriers at reduced gate voltage and thus, offering more carriers in channel region. Thus, the insulating layer material is significant for obtaining device characteristics and is required to be monitored for optimum performance.
PPT Slide
Lager Image
Inversion charge density versus gate voltage at different High-k dielectric.
PPT Slide
Lager Image
Threshold voltage variations with silicon radius at different High-k dielectric
The drain current variations with drain voltage and gate source voltage are plotted in Figs. 4 and Fig. 5 respectively with different gate insulator materials. High-k materials are found to be promising to increase the drive current performance of MOSFETs. This enhanced performance is due to the fact that replacing the conventional SiO 2 layer with high-k layer results in high value of surface potential, thus increasing the inversion charge density. This increased surface potential thereby reduces the device threshold voltage, and increased carrier density leads to rise in device field effect mobility. All these effects i.e. increased surface potential, increased carrier density, reduced threshold voltage and increases mobility results in high drain current for circuit driver applications.
PPT Slide
Lager Image
Variation of drain current with gate voltage of cylindrical surrounding gate MOSFETs at different dielectric thickness
PPT Slide
Lager Image
I-V characteristics of cylindrical surrounding gate MOSFETs with R=8nm.
The gate capacitance of a trigate MOSFET may be obtained from dQinv/dVgs . The gate capacitance variation with gate source voltage is shown in Fig. 6 . It can be seen from the figure that introduction of high-k materials though offers high peak of gate capacitance compare with SiO 2 layer characteristics. This is due to the fact that for the initial gate voltages approximately 0v to 0.2v, the value of gate capacitance is abruptly zero. Then the capacitance value shoots up to its maximum due to high-k dielectric layer. The device characteristic with SiO 2 layer superimposes on device characteristics with TiO 2 . This shows an alternative to conventional oxide layer, which is facing the physical limit to scale down in order to meet the requirement of scaled devices.
PPT Slide
Lager Image
Gate capacitance characteristics of different high-k dielectric for different gate voltages
5. Conclusion
In this paper, we have presented an analytical energy quantization model for SG MOSFETs. The model is based on the variational approach, which takes into account of the silicon radius and gate to source voltage. The model has been tested by means of a comparison with simulation data for various device dimensions and bias values. We have developed the analytical model for other device parameters like inversion charge centroid, inversion charge, and drain current by considering quantum effects. The influences of high-k dielectric constants are also included in this model to reduce the gate tunneling current. A very good agreement between the modeled and simulated data was achieved for all the device parameters. Therefore the proposed approach is appropriate for the modeling of surrounding MOSFETs in the presence of quantum effects.
Acknowledgements
The authors would like to thank the financial support by WOS-A Scheme, Department of Science and Technology, New Delhi, Government of India through the grant SR/ WOS-A/ET-41/2011.
BIO
P.Vimala received her B.E and M.E degrees, both in electronics and communication engineering from Anna University, Chennai, Tamilnadu, India. She is currently pursuing the Ph.D degree in the Department of Electronics and Communication Engineering, Thiagarajar College of Engineering (TCE), Tamilnadu, India under Anna University, Chennai. She is awarded with Women Scientist Fellowship for three years for her research work from Department of Science and Technology (DST), New Delhi, Government of India. Her research interests focused on Quantum Modeling and Simulation of Multigate MOSFETs. She is currently engaged in the modeling of Energy quantization in tri-gate MOSFETs and analyzing its performance using the Synopsis TCAD Simulation tool.
N. B. Balamurugan received his B.E and M.E degrees, both in electro-nics and communication engineering from the Thiagarajar College of Engineering (TCE), Tamilnadu, India. He has obtained his Ph.D degree in nanoelectronics at Anna University, India. From 1998 to 2004, he worked as a lecturer in R.V.S.college of engineering and technology, Tamilnadu, India. He is currently as an Assistant Professor in the Department of Electronics and Communication Engineering, Thiagarajar College of Engineering (TCE), Tamilnadu, India. He has published more than 50 research papers as sole or joint author in the field of device modeling and simulation. His research interests include modeling and simulation of semiconductor device structures.
References
He J. , Liu F. , Bia W. , Feng J. , Zhang J. , Zhang X. 2007 “An approximate carrier-based compact model for fully depleted surrounding-gate MOSFETs with a finite doping body” Semiconductor Science and Technology 22 (6) 671 - 677    DOI : 10.1088/0268-1242/22/6/015
Son A. , Kim J. , Jeong N. , Choi J. , Shin H. 2009 “Improved Explicit current-Voltage Model for Long- Channel Undoped Surrounding-Gate Metal Oxide Semiconductor Field Effect Transistor” Journal of Applied Physics 48 412 - 413
Chen Y. , Luo J. 2001 “A Comparative Study of Double-Gate and Surrounding-Gate MOSFETs in Strong Inversion and Accumulation Using An Analytical Model” International Conference on Modeling 1 546 - 549
Jimenez D. , Inguiez B. 2004 “Continuous analytic IV model for surrounding-gate MOSFETs” IEEE Electron Device Letters 25 (8) 571 - 573    DOI : 10.1109/LED.2004.831902
Colinge J-P. 2004 “Multiple-gate SOI MOSFETs” Solid-State Electronics 48 (9-10) 897 - 905    DOI : 10.1016/j.sse.2003.12.020
López-Villanueva J. A. , Cartujo-Cassinello P. , Gámiz F. , Banqueri J. , Palma A.J. 2000 “Effects of the inversion layer centroid on the performance of double gate MOSFETs” IEEE Transaction on Electron Devices 47 (1) 141 - 146    DOI : 10.1109/16.817579
López-Villanueva J. A , Cartujo-Cassinello P. , Banqueri J. , Gámiz F. , Rodriguez S. 1997 “Effects of the inversion layer centroid on MOSFET behavior” IEEE Transaction on Electron Devices 44 (11) 1915 - 1922    DOI : 10.1109/16.641361
Rodríguez N. , Gámiz F. , Roldán J.B. 2007 “Modelling of the inversion layer centroid and polysilicon depletion effects on ultrathin gate oxide MOSFET behaviour: the influence of the crystallographic orientation” IEEE Transaction on Electron Devices 54 (4) 723 - 731    DOI : 10.1109/TED.2007.891854
Schenk A. , Heiser G. 2014 “Modeling and Simulation of tunneling through ultra-thin gate dielectrics” Journal of Applied Physics 81 (12) 7900 - 7908    DOI : 10.1063/1.365364
Govoreanu B. , Blomme P. , Henson K. , Houdt J.V. , Meyer K.D. 2004 “An effective model for analyzing tunneling gate leakage currents through ultrathin oxides and high-k gate stacks from Si inversion layers” Solid State Electronics 48 (4) 617 - 625    DOI : 10.1016/j.sse.2003.09.031
Yuan Yn. , Yu Bo. , Song Jooyoung , Taur Yuan 2009 “An Analytical model for threshold voltage shift due to quantum confinment in surrounding gate MOSFETs with anisotropic effective mass” Solid-State Electronics 53 (2) 140 - 144    DOI : 10.1016/j.sse.2008.10.010
Wu Yu-Sheng , Su Pin 2009 “Analytical quantum - confinement model for short-channel Gate-All- Around MOSFETs Under Subthreshold Region” IEEE Transactions on Electron Devices 56 2720 - 2725    DOI : 10.1109/TED.2009.2030714
Moreno E. , Roldan J.B. , Barrera D. , Godoy A. , Gamiz F. 2010 “An Analytical Model for Square GAA MOSFETs including Quantum effects” Solid-State Electronics 54 (11) 1463 - 1469    DOI : 10.1016/j.sse.2010.05.032
Roldan J.B. , Godoy Andres , Gamiz Francisco , Balaguer M. 2008 “Modeling the Centroid and the Inversion Charge in Cylindrical Surrounding Gate MOSFETs Including Quantum Effects” IEEE Transcations on Electron Devices 55 (1) 411 - 416    DOI : 10.1109/TED.2007.911096
Vimala P. , Balamurugan N.B. 2013 “Modelling the centroid and charge density in double-gate MOSFETs including quatum effcets” International journal of Electronics 100 (9) 1283 - 1295    DOI : 10.1080/00207217.2012.743081
Iniguez Benjamin , Jimenez David , Roig Jaume , Lluis Hamdy A. Hamid , Marsal F. , Pallares Josep 2005 “Explict continuous model for long-channel undoped surrounding gate MOSFETs” IEEE Transcations on Electron Devices 44 (8) 1868 - 1873