Gap states controlled transmission through 1D Metal-Nanotube junction

Gap states controlled transmission through 1D Metal-Nanotube junction

D. Talukdar Saha Institute of Nuclear Physics, I/AF Bidhannagar, Kolkata 700 064, India Nanoscience Center, Department of Physics, P.O.Box 35, FI-40014 University of Jyväskylä, Finland    P. Yotprayoonsak Nanoscience Center, Department of Physics, P.O.Box 35, FI-40014 University of Jyväskylä, Finland    C.D. Mukherjee Saha Institute of Nuclear Physics, I/AF Bidhannagar, Kolkata 700 064, India    K. K. Bardhan Kalpana Chawla Centre for Space and Nano Sciences, 3F Swamiji Nagar, Kolkata 700 030, India    B. Karmakar Saha Institute of Nuclear Physics, I/AF Bidhannagar, Kolkata 700 064, India
August 31, 2019

Understanding the nature of metal/1D-semiconductor contacts such as metal/carbon nanotubes is a fundamental scientific and technological challenge for realizing high performance transistorsFrancois ; Franklin . A Schottky Barrier(SB) is usually formed at the interface of the metal electrode with the semiconducting carbon nanotube. As yet, experimentalAppenzeller ; Chen ; Heinze ; Derycke and numerical Leonard ; Jimenez studies have generally failedSvensson to come up with any functional relationship among the relevant variables affecting carrier transport across the SB owing to their unique geometries and complicated electrostatics. Here, we show that localized states called the metal induced gap states (MIGS)Tersoff ; Leonard already present in the barrier determines the transistor drain characteristics. These states seem to have little or no influence near the ON-state of the transistor but starts to affect the drain characteristics strongly as the OFF-state is approached. The role of MIGS is characterized by tracking the dynamics of the onset bias, of non-linear conduction in the drain characteristics with gate voltage . We find that varies with the zero-bias conductance for a gate bias as a power-law: with an exponent . The origin of this power-law relationship is tentatively suggested as a result of power-law variation of effective barrier height with , corroborated by previous theoretical and experimental resultsAppenzeller . The influence of MIGS states on transport is further verified independently by temperature dependent measurements. The unexpected scaling behavior seem to be very generic for metal/CNT contact providing an experimental forecast for designing state of the art CNT devices.

In nanoscale FETs depending upon the band bending at the metal-nanotube interface a Schottky Barrier(SB) is formed at the source and drain contacts. The charge injection in such a device is mainly determined by tunneling across the SB, whose width can be modulated by the gate voltage, . The resulting current through the SB (a function of ) is usually given by the Landauer-Bttiker formulaDatta :


where is the transmission probability across the junction at bias for carriers having energy . Linear - output characteristics are expected for this formula at low drain bias for ballistic tubes or wires with no defects. For a fixed gate bias , however one expects deviations from linear behavior with increasing drain bias at some onset bias, due to lowering of barrier heights leading to increased tunneling currents. At the onset of non-linear conduction , the electronic energy will be almost equal to the effective barrier height resulting in increased tunneling probabilitySimmons :


where is the source drain onset voltage and is the effective SB height which is a function of the . Eq. 2 implies that there is a voltage scale for non-linearity in the system and for each there will be a corresponding onset drain bias, which can be used as a measure of the effective barrier height .”Effective barrier height” is assumed as the flatband conditions are not applicable at finite temperatures, gate and drain bias and also contains the effect of other details e.g., metal-nanotube coupling as a function of curvature. In literature the standard approach Chen ; Appenzeller so far has been to calculate the effective SB height by self-consistently solving the Poisson and Schrdinger equation and experimentally by studying the device ON-current, . As non-ohmic behavior is usually obtained when there is increased tunneling due to barrier modification by the source drain bias , it is quite evident that the dynamics of the barrier change is contained in the drain characteristics -. The threshold bias for non-linearity can thus be used to characterize the SB. The question then is: Is there any general functional relationship between the relevant parameters affecting this change which should in principle be valid for all 1-D nanotransitors irrespective of device geometry?

The question can be addressed by analyzing the non-linear conduction properties, in particular the onset bias, as a function of gate voltage, . To find the onset or the crossover bias, there are mainly two methods used in literature. In the first method a crossover point is chosen at which the conductance deviates significantly from the zero field conductance by some arbitrary fraction , such that the conductance . In the second methodTalukdar11 ; Talukdar12 based on the existence of a single voltage scale , a set of curves can be characterized by a number , called the nonlinearity exponent. In this scaling approach, the conductance is given by the scaling relation:


where Q is a scaling function and the bias scale at each variable value, M (gate voltage) is given by the phenomenological relation


where is the conductance at zero bias at gate voltage and is a constant. After evaluating from the first method the scaling can be applied for the data following Eq. 3. Thus, both the methods are essentially equivalent Talukdar11 ; Talukdar12 for determination of the nonlinearity exponent . For CNTFETs the barrier is thin near the ON-state gradually becoming thicker towards the OFF-state. Accordingly, the threshold bias for non-linearity in the characteristics is also expected to vary systematically with . In the OFF state the onset bias for non-linearity should start at a higher voltage, than near ON-state where the corresponding should start earlier. Consequently, from Eq. 4 we expect that the non-linearity exponent to be negative for CNTFETs.

Figure 1: Figure 1. The transfer and ouput characteristics of the CNT transistor (a)Schematic and AFM image of the bottom gated CNT transistor with 300 nm oxide thickness. (b)Transfer characteristics of a CNTFET for different temperatures taken at =10 mV. changes weakly as temperature is lowered signifying increased tunneling at lower temperatures. (c) Output - characteristics for the same device at 120 K at different values of gate voltage in the subthreshold region. The curves become increasingly non-linear while approaching the OFF-state.

A typical CNT device architecture with a representative AFM and SEM image is shown in Fig. 1(a)-(c). The CNT is lying on a 300nm layer with Si as backgate connected with Pd contacts (see Methods for details). All measurements were done on small diameter tubes with d1.5nm as the SB height is inversely proportional to the nanotube diameterChen . More than four devices were studied at different temperatures in this study. Note that usually the fabrication process incorporates some additional tunneling barrier in the devices in addition to the SB already presentSvensson . The thickness of the material and diameter of the CNT is a crucial parameter for the SB of the CNTFET as they influence the electric field at the contact which controls its switching action.

Figure 1 (d) presents the variation of with as a function of temperature in a backgated CNT transistor with a channel length of 300 nm and diameter 1.4 nm. Resistance more than 6.5 k (quantum limit h/44.5k) or low value of ON-current are indicative of the existence of SB at the contacts for p-type hole transportJavey . Note that tunneling currents increases as becomes more negative indicating the thinning of barrier for current injection switching the transistor from OFF to ON state. Here, the current modulation by directly translates into corresponding change in . Below 120 K the -current becomes independent of temperature indicating that the dominant contribution to the current injection is due to tunneling with negligible contribution from thermal assisted tunnelingAppenzeller ; Martel ; Svensson . At this temperature - characteristics in Fig. 1 (e) contains complete information about the transport mechanism through the SB with minimal interference from thermal excitations. In the ON-state when the SB is thin the current increases rapidly with and is linear while near the OFF-state the rate of current increase is relatively slower and non-linear with . Although, the - characteristics usually give a qualitative understanding of the situation, for developing a framework to analyze the non-linear transport data, the equivalent conductance - description is usually more useful as described below.

Figure 2: Figure 2. Equivalent characterization of CNT transitors (a) Variation of the conductance vs. drain bias of the same data of figure 1 (c) at different values of gate bias. The dotted line schematically shows the motion of onset field as the gate bias is changed. The saturation of conductance seems to have a threshold of 160 mV which is due to emission of optical phonons. (b) The scaling of the same data as in figure 2a to achieve data collapse. There seems to be an envelope for the scaling curve from where deviations are due to optical phonon scattering.

The systematics of the evolution of - drain characteristics with changing is better illustrated by replotting them as - shown in Fig. 2 (a). The - curves are qualitatively identical, initially linear and after a threshold onset voltage start becoming non-linear. This general behavior of conductance for all measured tubes can be understood in the following manner. With increasing barrier thickness i.e., increasingly positive the - curves tend to become non-ohmic at smaller drain bias, . The systematics can be put in perspective by tracking the onset-bias for non-linear conduction as a function of , shown by the dotted line in the figure. It is clear that onset voltage is intricately related to the zero-bias conductance modulated by the gate bias. However the direction of movement of with is opposite to what is expected from Eq. (4). A distinctive feature of the - characteristics is the low values of onset voltages , a measure of the effective barrier heights, varying between mV. This is much lower than the effective barrier values obtained usually in literature which varies between Appenzeller ; Chen ; Heinze . The cause of this discrepancy will be discussed in a later section.

It is clear from Figure 2(a) that the curves for different can be collapsed into a single curve using the scaling procedure discussed earlier. The results of this scaling are illustrated in Fig. 2(b) where is represented as a function of . The detailed scaling procedure is described in greater detail elsewhereTalukdar11 ; Talukdar12 . The remarkable data collapse underlines the fact that there is a field scale for non-linearity in CNT transistors. The existence of the scaling is very significant here as it implies that the curves for a CNT transistors are identical differing by only a constant. The deviations from the scaled curve at high-bias (0.16 eV) is due to the well known electron backscattering due to emission of zone-boundary or optical phonons Yao ; Back ; Dresselhaus in single walled carbon nanotubes. The scaling holds for all the studied tubes and different temperatures.

Figure 3: Figure 3. Determination of functional relation ships for characterization of the Schottky barriers for hole injection: (a) Plot of onset voltage with zero bias conductance which is a function of gate voltage in the subthreshold region. The linear fits represent a power law between the two with the slopes mentioned. (b) log-log plot of effective schottky barrier height vs. representing competition between thermoionic emission and quantum mechanical tunneling obtained from Ref. Appenzeller . The plot indicates that the barrier height changes as a power law with gate voltage in CNTs.

To examine the scaling behavior, we plot and on a log-log plot as shown in Fig. 2(a). The solid line indicates that there exists a power-law relationship between and i.e.,


with a positive non-linearity exponent of as suggested by the orientation of the dotted line. The existence of the power-law automatically implies at high fields the conductance goes as a power-law with bias (see appendix), which is however masked in this case as conductance starts to saturate much earlier due to optical phonon scattering (0.16eV)Yao ; Back . We extract the values of non-linearity exponent for this tube at different temperatures (supplementary information). Surprisingly, the sign of the exponent is positive for all the measured tubes at different temperatures. Positive exponent implies that for thinner barriers (ON) higher drain voltage is required for the onset of non-linear conduction than near the OFF-state (thicker barriers). Existence of positive exponents in these tubes is counterintuitive to the expected negative exponent as discussed earlier in Eq. (4). With regard to this, it is of prime importance to know the dynamics of SB change with , as depends on the nature of SB which determines , implicitly dependent on . To understand this dynamics, we note that Appenzeller et. al.Appenzeller simulated and experimentally determined the effective barrier heights as a function of using the non-equilibrium Green’s function (NEGF) formalism and a thermal activation model respectively. Strikingly, their data shown in Fig. 3(b) indicates that the effective barrier height of the metal-CNT interface also varies as a power-law with , implying the power-law behavior is strongly linked to the modulation of the barrier.

Figure 4: Figure 4. Schematic for band alignment for different gate voltages at metal/nanotube interface (a) ON-state: Qualitative response of conduction through the schottky barrier at the metal/nanotube interface at a gate voltage close to threshold voltage where the transistor ON. The barrier shown by the colored band is sufficiently thin assisting direct tunneling (b) Intermediate state: The band bends for a gate voltage below the threshold where the SB is sufficiently thick making the MIGS slowly relevant for transport (c) OFF state: The SB is the thickest and the transport through the localized MIGS states dominates.

To understand this anomalous behavior and account for the origin of functional relationships obtained we need to examine the different band bending scenarios in the CNT device for different gate bias conditions ranging from OFF to ON. In Fig. 4(a) the CNT-device is in an ON-state in 4(b) at an intermediate gate bias and for 4 (c) the device is OFF. As hole transport is considered we mainly concentrate on the the valence band. The shaded bands represent the SB on the source and drain sides respectively. In figure 4 (a) when the transistor is ON the barrier thickness is reduced leading to direct tunneling and subsequently lower values of for non-linear conduction. For close to OFF-state the barrier is considerably thicker resulting in lower tunneling currents in turn requiring larger for non-linear conduction. However, the experimental situation is completely opposite. The discrepancy in the non-linear behavior, manifested distinctly in the dependence of on the barrier thickness, and thus conductance demonstrates the role of the localized states for transport. These states located in the metal-CNT interface are the Metal Induced gap states (MIGS)Tersoff ; Tersoff1 ; Heine ; Louie whose role in carrier transport was not recognized in literature. These states arise as a result of electronic wave functions being abruptly terminated at the metal-CNT junction and extends only up to a few nanometer into the CNTs getting localized close to the interface. These interface states have little or no effect on electrical behavior near the ON-state when the barrier is thin as transport is through direct tunneling through the barrier leading to a particular onset bias for non-linearity . However, near the OFF-state when the barrier is sufficiently thick the dominant conduction mechanism crosses over from direct tunneling to transport through localized states in the barrier as function of barrier thickness. This can explain the lower thresholds for non-linear conduction with increasing barrier heights as shown earlier which is a characteristic of tunneling through localized statesXu . The localized states at the interface states thus play a crucial role for controlling the device characteristics in CNTFETs and sets the scale for non-linearity in these devices. The MIGS states present in the barrier can also explain the asymmetry of the characteristics commonly obtained for carbon nanotube transistors Li ; Khoo which is due to different configurations of MIGS states present in each metal-semiconductor barrier on the source and the drain side.

The origin of the power-law relationship between and might arise due to the variation of barrier height as a power-law with having an exponent ’y’ as obtained by Appenzeller et. al. Physically, the exponent ’y’ indicates a ratio of and the electric field at the contact. As depends on the the configuration of localized states at the interface varying with as a power-law it indicates the ratio of and the distribution of localized states at the interface.

Figure 5: Figure 5. Temperature dependence for a CNTFET (a) The experimental transfer chracteristics for different temperatures at = 10 mV. varies by an order of magnitude with temperature. The inset shows the same transfer characteristics with values in logarithmic scale. (b) Output characteristics plotted as at different temperatures. The onset voltage start decreasing with lowering temperatures shown by the dotted line. At 4.2K the is completely non-linear. (c) Scaling of the same data leading to data collapse as shown. and are the zero bias conductance and onset bias at some temperature T. The inset shows the log-log plots of and with linear fits and slopes mentioned.

To verify the physical reasoning described above we made additional temperature dependent measurements on our devices. In Figure 5 (a) the gate characteristics of a device in both linear- and log-scale is shown. The drain characteristics - at zero gate bias is shown in Figure 5 (b). Zero gate bias is chosen so as to keep the barrier without being affected by the gate which will enable us to see the intrinsic barrier characteristics. Also for illustration of this method it is advantageous to take the drain characteristics - at zero gate bias as they are more non-linear. It can be seen from the figure that the drain curves start to become progressively non-linear with decreasing temperature. To find the onset bias shown by the dashed line, we can use the same scaling procedure described above obtain a perfect data collapse for normalized conductance vs. the normalized bias . The bias scale appears to be a power-law function of the Ohmic conductance and decreases with decreasing temperature as shown in the inset of Fig. 5c (circles). This is similar to the situation where the conductance was modulated by . This can be understood qualitatively in the following manner. As the temperature decreases the contribution of thermally assisted tunneling for transport is reduced. Initially - curves are therefore ohmic at high temperature as the main contribution to carrier injection is by thermal activation. However, with decreasing temperature, thermal activation is reduced significantly and the effect of localized states in the barrier can no more be neglected. With progressive decrease of temperature more and more localized states in the barrier (MIGS) start taking part in conduction leading to lower onset bias with decrease of temperature. However, it is still not clear why varies with conductance as a power law. It must be mentioned here that the curve at is not used for the scaling as the conduction mechanism is different in this case due to the onset of coulomb blockadePostma ; Park effects.

In conclusion, our experiments provide first clear experimental evidence that even though the transistor action observed in nanoscale CNTFETs is based on transport across a SB barrier the - curves can only be explained if we consider the barriers to contain localized MIGS states. Using experiments and comparing with previous simulations and experimental data from literature we suggest that the effective SB heights and threshold bias for non-linearity in CNTFETs are modulated by the gate voltage/source-drain voltage in a power-law fashion with a nonlinearity exponent . Further, the formalism can explain the discrepancy commonly observed between theoretical and experimental SB values for any low dimensional material. Our method of analyzing nonlinear transport data can also be used to characterize the cut off working linear regime voltage of a CNTFET or nanowire FET device.

Methods In our studies the SWNT devices used here fabricated by spin depositing suspension of nanotubes with (1,2)-dichloroethane onto highly p-doped Silicon wafers covered with 300 nm thick thermally grown . AFM (atomic force microscope) was used for locating the tubes and standard e-beam lithography steps were used to electrically contact the nanotubes. Pd of thickness 25 nm was deposited over the CNTs in order to improve the contact resistance. Low temperature transport measurements were performed in a RF-shielded room using home made dipstick.

Appendix The field-dependent conductance as shown in Figure. 2a and 5b clearly show that at high fields conductance becomes independent of the variational parameter (Gate bias and Temperature respectively). If the conductance is continued to be described by Eq. (3) it is obvious that at large fields the conductance varies with bias as a power-law: where ’’ is the high field exponent and .

We acknowledge helpful discussions with Prof. M. Ahlskog.


  • (1) L. Francois, and A. A. Talin, Nat. Nanotech. 6, 773 (2011).
  • (2) A. D. Franklin, and Z. Chen, Nat. Nanotech. 5, 858 (2010).
  • (3) J. Appenzeller, M. Radosavljevic’, J. Knoch, and P. Avouris., Phys. Rev. lett. 92 048301 (2004).
  • (4) Z. Chen, J. Appenzeller, J. Knoch, Y.M. Lin, and P. Avouris, Nano Lett. 5, 1497 (2005).
  • (5) W. Kim, A. Javey, R. Tu, J. Cao, Q.Wang, and H. Dai., App. Phys. Lett. 87 173101 (2005).
  • (6) S. Heinze, J. Tersoff, R. Martel, V. Derycke, J. Appenzeller, and P. Avouris, Phys. Rev. Lett. 89, 106801 (2002).
  • (7) V. Derycke, R. Martel, J. Appenzeller, and P. Avouris, Appl. Phys. Lett. 80, 2773 (2002).
  • (8) J. Svensson, and E.B. Campbell, J. Appl. Phys. 110, 111101 (2011).
  • (9) J. Tersoff, Phys. Rev. lett. 52 465 (1984).
  • (10) F. L onard, and J. Tersoff. Phys. Rev. lett. 84 4693 (2000).
  • (11) D. Jimenez, X. Cartoixa, E. Miranda, J. Sune, F. A. Chaves, and S. Roche, Nanotechnology 18, 025201 (2007).
  • (12) S. Datta, Electronic transport in mesoscopic systems ( Cambridge university press, 1997).
  • (13) J. G. Simmons, J. Appl. Phys. 34, 1793-1803 (1963).
  • (14) S.M. Sze, and K. K. Ng , Physics of semiconductor devices ( Wiley. com, 2006).
  • (15) D. Talukdar, U. N. Nandi, K. K. Bardhan, C. C. Bufon, T. Heinzel, A. De, and C. D. Mukherjee, Phys. Rev. B 84,054205 (2011).
  • (16) K. K. Bardhan, D. Talukdar, U. N. Nandi, and C. D. Mukherjee, Accepted in Phys. Rev. B.
  • (17) Z. Yao, C. L. Kane, and C. Dekker, Phys. Rev. Lett. 84, 2941 (2000).
  • (18) J. H. Back, C. L. Tsai, S. Kim, S. Mohammadi, and M. Shim, Phys. Rev. Lett. 103, 215501 (2009).
  • (19) M. S. Dresselhaus, G. Dresselhaus, R. Saito, A. Jorio, Phys. Rep. 409, 47 (2005).
  • (20) F. leonard, and J. Tersoff, Phys. Rev. Lett. 84, 4693 (2000).
  • (21) V. Heine, , Phys. Rev. 138 A1689 (1965).
  • (22) S.G..Louie, and M.L.Cohen, Phys. Rev. B 13 2461 (1976).
  • (23) Y. Xu, D. Ephron, and M. R. Beasley, Phys. Rev. B 52 2843 (1995).
  • (24) H. Li, Q. Zhang, and N. Marzari, Nano letters 8 64 (2008).
  • (25) K.H. Khoo, J. B. Neaton, Y. W. Son, M. L. Cohen, and S. G. Louie, textNano letters 8 2900 (2008).
  • (26) A. Javey, J. Guo, Q. Wang, M. Lundstrom, and H. Dai., Nature 424, 654 (2003).
  • (27) R. Martel, V. Derycke, C. Lavoie, J. Appenzeller, K. K. Chan, J. Tersoff, and Ph. Avouris., Phys. Rev. lett. 87 256805 (2001).
  • (28) H. Postma, T. Teepen, Z. Yao, M. Grifoni, and C. Dekker., Science 293 76 (2001).
  • (29) J. W. Park, A. N. Pasupathy, J. I. Goldsmith, C. Chang, Y. Yaish, J. R. Petta, M Rinkoski., Nature 417 722 (2002).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description