This paper proposes an analytical model for the asymmetric double-gate (ADG) TFET operation. Two-dimensional Poisson equation is solved based on the depletion approximation [1,2]. The two-dimensional potential function is given by appropriate boundary conditions. Internal electric field is calculated from the potential model and tunnel current is numerically calculated using the Kane’s tunnel current generation model [3,4]. We used the commercial TCAD simulator in order to examine the analytical model [5]. Schematic of ADG-TFET device structure is shown in Fig. 1. It is assumed that p-type source and n-type drain regions are highly doped. Channel region is a p-type with the doping level of Na . tox denotes the gate oxide thickness and tsi denotes the semiconductor film thickness. Hereafter, we call the top gate “front gate” and the bottom gate “back gate”. Potential profile of the channel region is derived from the following two-dimensional Poisson equation based on the depletion approximation. Equation (1) where q denotes the elemental charge and εsi denotes the permittivity of semiconductor. Here, we assume the following potential function with the cubic function of coordinate y [1]. Equation (2) where four factors a(x), b(x), c(x), and d(x) are functions of coordinate x. We introduce the following boundary conditions for the potential in solving eq. (1). Equation (3) where ψf(x) and ψb(x) stand for front and back surface potentials at y=0 and y=t si, respectively. Vfg,fb and Vbg,fb are the front-gate and the back-gate flat-band voltages, respectively. Vfg and Vbg are the front- and back-gate voltages, respectively. εox is the permittivity of gate insulator. Drain current of ADG TFET is numerically calculated by integrating carriers tunneling from the valence band of the source region to the conduction band of the channel region. We can calculate the drain current (IDS ) using the carrier generation rate (G) due to the band-to-band tunnel process. Equation (4) where we use the Kane’s model for the carrier generation rate (G) [3,4]. Equation (5) where |E| (=(Ex 2+Ey 2)1/2 ) is the electric field strength, A and B are constants given byEquation (Si) Equation (Ge) We perform the TCAD simulation [5] in order to examine the validity of the analytical model. In device simulations for Si-based TFETs, we assume the gate material with the work function of 4.1 eV for the front- and the back-gate electrodes. Figure 2 shows the calculated potential profiles at y=0nm, 2.5nm, and 5.0nm for tsi =10nm. Using the potential function revealed in Fig. 2, local electric field values are calculated along with the x direction and the y direction separately. The potential function analytically given and resulting electric fields successfully reproduced those given by TCAD simulations (not shown here). We found a slight difference of potential profile near the surface around the source junction for a relatively low-doping source region because of the depletion in the source region [6]. We also calculated drain current characteristics. Figure 3 demonstrates ID-Vfg characteristics for tsi =10nm and 20nm, where the condition of Lg >2ts is satisfied. In Fig. 5, the model basically reproduces the TCAD result. On the other hand, the model doesn’t well reproduce TCAD-based ID-Vgf characteristics for tsi =30nm and 100nm (not shown here), where Lg <2tsi . This difference stems from the linear approximation of the in-depth electric field in the model. As the condition of Lg >2tsi is important in the future, this is not a crucial issue of the model. In addition, we confirmed the model proposed here is available for Ge-based TFET (not shown here).
Read full abstract