Investigating the effect of some parameters of the channel on the characteristics of tunneling carbon nanotube field-effect transistor

This paper studies p-i-n tunneling carbon nanotube field-effect transistor to investigate the effect of various parameters of the channel on the characteristics of tunneling carbon nanotube field-effect transistor. Tunneling carbon nanotube field-effect transistor (T-CNTFET) has been simulated using non-equilibrium Green’s function (NEGF), and the transmission was conducted through inelastic scattering. Besides the evaluation of device performance, various parameters of the channel were also compared. One of the parameters is considered as the variable, while other parameters of the channel are constant. Then, improved characteristics were discussed by selection of some channel parameters. T-CNTFET with CNT (10, 0) with oxide thickness = 1 nm shows reduced sub-threshold swing (18 mV/decade).


Introduction
While technology-based silicon will reveal its limitations, some emerging devices are investigated to find appropriate substitution to silicon. Carbon nanotube FETs (CNTFETs) are illustrated to replace silicon in the future [1]. Carbon nanotubes were explored by Sumio Iijima in 1991 [2]. Carbon nanotubes have properties, such as caring high current, ballistic transport, mechanical stability, and dynamic load, which make them appropriate substances to replace silicon as the channel of field-effect transistors. Ballistic transport in p-type CNTFETs using quasi-planar gate geometries and higher on-currents to MOSFET has been proved by experiments [3,4]. Cylindrical carbon nanotube with gate-all-around (GAA) leads to an increase in miniaturization and also provides electrostatic control of the gate [5]. There is charge pile-up in the channel CNTFET, which leads to the increase in the off-current and lowers the I ON /I OFF ratio. To overcome this problem, a gate-controlled tunneling device (T-CNTFET) is offered. A steep inverse sub-threshold slope and the better controlled off-current were obtained by T-CNTFETs. In T-CNTFET, the increase in the potential difference between the gate and the drain contact leads to strong band-bending near the drain contact [6]. The T-CNTFET has been simulated by self-consistent solution between the Poisson and Schrodinger equations [3]. As can be found from [3], to suppress the ambipolar behavior of CNTFET and to improve the performance of this device, the asymmetric double-halo doping has been applied in CNTFET. In addition, for this purpose, CNTFETs with graded double-halo channel are presented as in [7]. Transmission is considered ballistic in both [3] and [7]. In this paper, we investigate the effect of changing a channel parameter on characteristics of T-CNTFET, provided that the other channel parameters remain constant. The investigated channel parameters include channel length (L), channel diameter, and gate oxide thickness. Finally, on-current (Ion), tunneling current, on/off-current ratio (I ON /I OFF ), and sub-threshold swing (SS) were the characteristics investigated in this research.
Device structure Figure 1 demonstrates the illustration of T-CNTFET structure. The CNT is made of a zigzag (16, 0) of 0.63-nm radius and length of L = 15 nm. A coaxial gate is placed around the channel region of the nanotube. The dielectric layer HfO 2 has 2 nm thickness, and the relative dielectric constant is e r = 16. HfO 2 separates the gate from CNT 2 nm. T-CNTFET included in p-doped source and n-doped drain regions with 30 nm length which have been doped by 2 nm -1 doping.

Computational method
The proposed device is simulated by the self-consistent solution of the poison and Schrodinger equations within the NEGF formalism. The nanoscale system is simulated by the NEGF method out of equilibrium. The electrostatic potential is obtained by the Poisson equation to compute the Hamiltonian of the system [3]. Carbon nanotube is coaxially gated, so symmetric characteristics of cylindrical coordinates create the potential, and the charge is constant around the nanotube. Consequently, the Poisson equation is essentially a 2D problem along the tube (z-direction) and the redial direction (r-direction). Eigen states around the tube circumferential direction (modes) are plane wave vectors satisfying the periodic boundary condition [7]: where U J (r, z), e, and pðr; Z j Þ are the electrostatic potential, the constant dielectric, and the net charge density distribution, which include dopant density as well, respectively. The retarded Green's function for the device in matrix form is computed as where P S and P D are the self-energies of the source and the drain, respectively. g þ , E, I, and H are an infinitesimal positive value, the energy, the identity matrix, and the Hamiltonian of the CNT, respectively: where b 2q = 2tcos (pq/n), t & 3 eV, and N are the nearest neighbor hopping parameter and the total number of carbon rings along the device, respectively. The diagonal elements U J are assumed as the no-site electrostatic potential along the tube surface as in [8]. The self-energy matrix can be interpreted as a boundary condition of the Schrodinger equation. In this paper, we investigated a self-energy for semi-infinite that leads to boundary conditions, which enables to consider the CNT as connected to infinitely long CNTs at its ends. The source self-energy function P S has all its entries zero except for the (1, 1) factor: In addition, P D has only its (N, N) element nonzero, and it is given by an equation similar to E q with U 1 displaced by U N as in [8]. The two-dimensional nanotube lattice of a (n, 0) zigzag CNT was transformed to n decoupled one dimensional modes by doing a basis transform from the real space to the mode space in the circumferential direction. Under typical bias conditions, the few modes that are relevant to electronic transport are treated. For this mode, the charge density is calculated by integrating the local density-of-states (LDOS) over energy: where e, sgn (E), E FS,D , D i S/D, and (E, z) are the electron charge, the sign function, the source (drain) Fermi level, the LDOS due to the source (drain) contact as computed by the NEGF method, respectively. In this paper, we have used single-pi model for all simulations; nanotube conduction and valence bands are considered symmetric, so the charge neutrality level, E N (z), lies at the middle of bandgap. The potentials at the source/drain and gate electrodes are assumed fixed as the boundary conditions. The iteration between the atomistic quantum transport equation and the electrostatic equation continues until self-consistency is obtained, and then, Landaur-Buttiker formula computes the source-drain current: In this formula, T(E) is transmission coefficient [7].
Results and discussion Figure 2 shows the behavior of T-CNTFET clearly. The carriers can tunnel from the source to the drain provided that conduction band (E C ) in the channel is higher than the valence band (E V ) of the source. For the equilibrium condition, the whole length of the channel region acts as a tunneling barrier. On the other hand, when E C in the channel is lower than the E V in the source, the height of the tunneling barrier will decrease at the source-channel side and as a result of band-to-band tunneling current (applied V DS [ 0). The phonon-induced states can create further tunneling paths of carriers that lead to higher tunneling current under scattering as in [6]. Thus, the transmission was assumed with inelastic scattering. Figure 3 shows that when the device is on, the drain current saturates at high V DS . But increasing V DS would finally result in tunneling at the drain-channel side. This ambipolar behavior limits the performance of the device. As V GS increases, the bending of the bands at the source side will increase more and consequently causes tunneling barrier thinner and rises the transmission r as well. It is effective less than 60 mV/ decade swing, and leads to increase on-current essentially (Fig. 4). On the other hand, applying a negative V GS leads to band-to-band tunneling at the drain-channel side, and results in ambipolar conduction. There is the current minimum at V GS = V DS /2 for symmetric source/drain doping. Figure 4 shows I DS -V GS characteristics for T-CNTFET at different V DS values. As can be observed from Fig. 4, the drain current saturates at high V DS for T-CNTFET. The ambipolar behavior leads to the high off-current, also the sub-threshold slope and the on-current increase. Due to the rise of off-current, the I ON /I OFF ratio reduces subsequently. Therefore, we investigate various channel parameters, to achieve improved characteristics. Then, with choosing some channel parameters, we improve some characteristics of T-CNTFET. The diameter of CNT is inversely proportional to bandgap energy. The bandgap energy is the minimum energy needed to break a covalent bond in the semiconductor crystal, so as to free an electron for conduction. This means that lower bandgap energy creates a better conduction material [9,10].
Diameter CNT is acquired from the following equation: That a 0 = 0.425 nm is carbon-to-carbon distance, n and m are described CNT diameter and zero, respectively. The threshold voltage (V th ) is described as the voltage required to switch on a transistor. A threshold voltage of the intrinsic CNT channel can be approximated as the halfbandgap that is an inverse function of the diameter: where a = 2.49 Å is the carbon-to-carbon distance, V p = 3.033 eV is the carbon p-p bond energy in the tight bonding model, e is the unit electron charge, and D CNT is the CNT diameter [10]. Table 1 shows that as CNT diameter increases, threshold voltage decreases and leads to more conduction. The threshold voltage is defined for various values of n, as shown in Table 1. Electronics industry is looking for more reductions in the size of MOSFETs. The MOSFET channel length is a measure for the size of its other parameters. The final proposal is to achieve the channel length 10 nm and the less. In this structure, coaxial gate covers all around channel thus are defined as one of the channel parameters. Figure 5 is investigated with various CNT diameter at a constant L = 15 nm and oxide thickness = 2 nm, different channel length at the constant channel diameter (n = 16) and oxide thickness = 2 nm, and also different oxide thickness at the fixed L = 15 nm and CNT (16, 0) for I DS -V GS of the T-CNTFET. As can be seen from Fig. 5, on-current and tunneling current increase due to increase in CNT diameter at the constant L = 15 nm and oxide thickness = 2 nm. The on-current and tunneling current increase due to the fact that the CNT diameter is increased, so the bandgap energy lowers which result in the higher conduction. In addition, Fig. 5 demonstrates, for constant CNT diameter and CNT length, reduction oxide thickness leads to the increase in I ON and tunneling current. Thinner gate oxide results in the reduction of gate-source fringing. Due to this reduction, the height of the barrier energy decreases at both source-channel and drain-channel sides which result in lower tunneling barriers and I ON . Gate-drain fringing can repress the ambipolar conduction but degenerates highfrequency operation by increased fringe capacitances. As CNT diameter and gate oxide thickness are constant, I ON and tunneling current do not vary by reducing channel length (for inelastic scattering transport). For further observation, we investigated I DS -V GS figure for the range   Fig. 6, separately. Figure 6 demonstrates that for L = 15, 20, and 30 nm, the reduction of channel length does not reduce I ON but I OFF increases, because direct tunneling from source to drain for shortchannel length (L G B 10 nm). As can be found from  Figure 7 shows the I ON /I OFF characteristics for T-CNTFET with all various parameters of the channel. As can be observed from Fig. 7 with decreasing CNT diameter, I ON /I OFF ratio increases. The reason of this increase is that with reduction of the CNT diameter, the off-current reduced, it leads to higher I ON /I OFF ratio, which causes an increased I ON /I OFF ratio for the T-CNTFET with smaller channel diameter. For CNTs (10, 0), (13, 0), (16, 0), and (19, 0), I ON /I OFF ratio is [10 6 , [10 4 , [10 3 , and [10 2 , respectively. Furthermore, Fig. 7 illustrates for the range of oxide thickness include in 1, 1.5, and 2 nm provided that other parameters are constant. ON/ OFF-current ratio is constant and nearly [10 3 . In addition, it demonstrates for L = 15, 20, and 30 nm, I ON /I OFF ratio is nearly constant, and [10 3 which is due to fixed I ON and I OFF for various channel lengths. For further evaluation, Fig. 8 was investigated. As can be seen from Fig. 8 due to large off-current L = 10 nm compared to other channel length, I ON /I OFF is decreased and its value is [10 2 . The largest I ON /I OFF ratio are observed by T-CNTFET with CNT (10, 0), L = 15 nm, and oxide thickness = 2 nm. The sub-threshold swing (SS) is defined as the change in gate bias needed to change the sub-threshold drain current by one decade, as given by: I DS and V GS are descried as drain-source current and gate-source voltage, respectively [5,11]. Figure 9 demonstrates sub-threshold swing as a function of the gatesource voltage for variation in CNT diameter, provided that other channel parameter is constant. As can be found from Fig. 9 as channel diameter decreases, the band gap energy will increase more, making the tunneling barrier increasingly ticker and leads to the decrease in the drain   transconductance also threshold voltage increases and results in on-current/sub-threshold swing decrease/ increase. This causes less than 60 mV/decade sub-threshold swing. For CNTs (10, 0) and (13, 0), sub-threshold swing is 53 and 56 mV/decade, respectively. In addition, for CNTs (16, 0) and (19, 0), sub-threshold swing is not achieved less than 60 mV/decade swing. Reduction of oxide thickness does not affect sub-threshold swing. Figure 10 shows that subthreshold swing is constant while channel length is varied. Swing sub-threshold is 60 mV/decade for the range of the channel lengths 30, 20, and 15 nm, respectively. For L = 10 nm, sub-threshold swing is not observed less than 60 mV/ decade which that is undesirable. Reason of this increase is high off-current for the channel length 10 nm. Figure 11 demonstrates that for various values of oxide thickness include in 2, 1.5, and 1 nm, SS achieves 62 mV/decade which more than 60 mV/decade conventional limits. To improve characteristics of tunneling carbon nanotube fieldeffect transistor, we select some channel parameters, and then investigate characteristics of the T-CNTFET regarding the selected parameters and compared to T-CNTFET with default parameters. As can be observed from Fig. 12, reduction channel diameter causes the reduction at off-current. For CNT with same diameter (10, 0), reduction oxide thickness leads to reduction at off-current that results in higher I ON /I OFF ratio, nearly [10 7 . In addition, sub-threshold swing of T-CNTFET with CNT (16, 0) and (10, 0) with various oxide thickness is compared in Fig. 13. It is evident    Fig. 13 that T-CNTFET with CNT (10, 0) with oxide thickness = 1 nm shows reduced sub-threshold swing (18 mV/decade) in comparison with T-CNTFET with CNT (10, 0) with oxide thickness = 2 nm and T-CNTFET with CNT (16, 0). Swing sub-threshold is decreased rather than [6]. This is attributed to control of the gate voltage of T-CNTFET (10, 0) with oxide thickness = 1 nm over the channel.

Conclusion
The NEGF formalism has been applied to simulate the electronic properties of the T-CNTFET structure with various channel parameters. The simulated characteristics were compared with those of the changing one of the channel parameter if other channel parameters are constant, the improved characteristics were discussed, finally. Simulations show that using CNT with larger bandgap leads to lower subthreshold swing and higher I ON /I OFF -current, while reducing on-current and tunneling current. In this paper, the highest value of the I ON /I OFF ratio ([10 7 ) and the least sub-threshold swing (18 mV/decade) observed by T-CNTFET with CNT (10, 0), oxide thickness = 2 nm and CNT (10, 0), and oxide thickness = 1 nm, respectively. Finally, simulation results demonstrate that changing the channel parameters alters characteristics of the devices, essentially.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://crea tivecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.