Contact resistance in graphene channel transistors
Contact resistance in graphene channel transistors
Carbon letters. 2013. Jul, 14(3): 162-170
Copyright ©2013, Korean Carbon Society
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : May 24, 2013
  • Accepted : June 30, 2013
  • Published : July 31, 2013
Export by style
Cited by
About the Authors
Seung Min, Song
Byung Jin, Cho

The performance of graphene-based electronic devices is critically affected by the quality of the graphene-metal contact. The understanding of graphene-metal is therefore critical for the successful development of graphene-based electronic devices, especially field-effect-transistors. Here, we provide a review of the peculiar properties of graphene-metal contacts, including work function pinning, the charge transport mechanism, the impact of the process on the contract resistance, and other factors.
1. Introduction
Since the first demonstration of top-gated graphene transistor in 2007 [1] , there has been major progress in the development of graphene devices, especially radio-frequency transistors [2 - 18] . The most of the improvement was achieved by the development of the gate oxide deposition process, exemplified by inserting interfacial layers or thin Al film before the atomic layer deposition process [5 , 19 , 20] ; the formation of relatively low contact resistance with Pd or Ti [21 , 22] , or the implementation of short-channel devices by a self-alignment process [10 , 16 , 17 , 23 , 24] . Although significant advances in terms of performance have been made, the contact resistance is still poorly understood despite the fact that it is one of the most critical barriers preventing further improvements [21 , 25 - 33] . In this paper, a short review of the properties of graphene-metal contacts is given. The contents include work function pinning as induced by a metal contact, charge transport at the contact region, and a discussion of the impact of the process conditions on the characterization process. It should be emphasized that the dominant factors which significantly determine the intrinsic contact resistance have not been clearly identified. To exploit the superior electrical properties of graphene fully, a reliable and the reproducible metal contact process scheme must be developed in the future.
2. The work function pinning of graphene under metal
When a metal and a semiconductor come into contact, charges are transferred from higher energy states to lower energy states until the energy levels are balanced in an equilibrium state. Due to the small density of states in a semiconductor, there is energy band bending along the vertical direction, as shown in Fig. 1 a, in what is known as the depletion region with depletion width W d . In a similar manner, graphene can also be modeled as a semiconductor with a zero band gap, and band bending is necessary if graphene has a finite thickness as a bulk semiconductor materials. Because graphene has no band gap, instead of an ‘energy band,’ the ‘Dirac energy level’ (E Dr ) would be a proper term for graphene, as illustrated in Fig. 1 b. Dirac energy level bending in graphene along the lateral direction does not occur, but there is a work function shift in graphene ( Fig. 1 c). It was first theoretically estimated that graphene on various crystalline metals can be doped due to the presence of an interface
PPT Slide
Lager Image
Energy band diagram of (a) a metal/semiconductor, (b) metal/three-dimensional graphene, and (c) metal/graphene. EDr is the Dirac energy level, ΔΦb is the Schottky barrier height, ΔΦgr/m is the work function change due to the metal-induced doping, and Wd and LB are respectively the depletion width and the length of the charge transfer region. The red solid arrows indicate the path of the carriers. (d) Schematic of the energy band diagram in graphene at the contact region. Graphene underneath metal is doped in holes and the channel graphene is negatively doped by the back-gate voltage. The thicker EDr underneath the metal represents the broadening of the density of states near the Dirac energy level. (e) The model of the graphene-metal contact includes the contact resistivity (ρc) and sheet resistance (Rs). Transfer length (LT) is the effective contact length and is defined with the distance where the potential drops to 1/e times from the edge..
dipole layer induced by the charge transfer [34] . The metals were classified in two types, and the characteristics of graphene were found to be significantly altered by chemisorption on Co, Ni, and Pd. In contrast, graphene weakly bonds to Al, Cu, Ag, Au, and Pt and preserves its conical electronic structure. The weak bonding leads to a work function shift in graphene, and the amount of this shift is determined by the work function of the corresponding metal.
In addition to the shift of the work function of graphene underneath metal, there is a gradual bending of the Dirac energy level along the horizontal direction near the contact edge due to the work function difference between the exposed graphene and the graphene underneath the metal, as shown in Fig. 1 d [35] . This region, known as the charge transfer region, has been successfully measured by Kelvin probe force microscopy [36 , 37] , and photocurrent experiments [38 - 41] . Clear Dirac energy level shifts have been observed near the edge of the graphene–metal contact, and metal-induced work function pinning has also been clearly identified. The presence of a charge transfer region leads to the formation of a p-n junction in the graphene channel, which increases the resistance in graphene compared to a p-p’ or n-n’ junction [42] . This introduces important electrical characteristics in graphene, such as asymmetric conductance behaviour coupled with back–gate voltage dependence [21 , 25 , 42 - 44] . It also leads to errors during the determination of the sheet resistance and contact resistance during transmission length measurements (TLM) or four-point probe measurements (see the details in Section 4) [43 , 45] . Fig. 2 shows schematic illustrations of the energy band bending and electrical characteristics with the back gate voltage. If there is no metal-induced doping and if the work function at the contact region is independent of the back gate voltage, the V G -I D characteristic curve will show symmetry, as depicted in Fig. 2 d. In contrast, the curve is asymmetric with respect to the Dirac voltage and the hole current is usually higher than the electron conductance when a contact is formed by Pd or Ni [21 , 41 , 44 , 46] . This indicates that the work function of graphene under the corresponding metal is pinned to a higher value than its original work function (that of the exposed graphene), and higher resistance is observed when the external field lowers the Dirac energy level of the channel below Fermi level, thereby leading to the formation of a p-n junction, as illustrated in Figs. 2 b and e.
To understand the detailed carrier transport behavior of graphene field-effect transistors (FETs), the precise work function of graphene under metal is one type of basic information. Previously reported values measured by scanning photocurrent microscopy were evaluated from indirect measurements, e.g., by finding the gate voltage that yields a flat channel potential profile [38 - 41] . As the flat band voltage ( V FB ) of a metal-oxidesemiconductor (MOS) capacitor is directly related to the work function of the gate material, the work function of graphene under metal ( Ф gr/m ) can be measured by capacitance-voltage ( C-V ) measurements with good accuracy. A detailed and thorough analysis of C-V measurements using a metal-graphene-oxidesemiconductor (MGOS) capacitor structure was, therefore, conducted with various metal electrodes [47] . It is interesting that the work function of graphene under Cr/Au coincides with that of Cr/Au. This indicates that the work function of graphene follows the work function of the overlying metal. The same result was also observed from a Ni contact. In contrast, the work function of graphene under a noble metal such as Pd/Au or Au does not match that of the corresponding metal. The exact work
PPT Slide
Lager Image
Schematics of the modulation of EDr in graphene by the back gate voltage at three different conditions at the contact region: (a) work function pinning without doping, (b) work function pinning and doping, and (c) doping without work function pinning. The blue dotted lines are the results with a negative voltage and red dotted lines are those at a positive voltage. The channel and contact graphene are assumed to be doped with ΔEgr and ΔΦgr/m, respectively. (d-f ) Schematics of the electrical characteristics under each condition.
PPT Slide
Lager Image
(a) Capacitance-voltage (C-Vg) measurement results with various Cr/Au area of metal-oxide-semiconductor capacitor which is illustrated in the inset. (b) The differentiated curves in Fig. 3a clearly show two peaks. (c) The peak heights have strong dependence on the metal area, indicating the presence of two different work functions in graphene underneath the metal and the exposed graphene. (d) The extracted work function of graphene underneath the metal (Φgr/m) and that of the corresponding metals (Φm) (Reprinted with permission from earlier work [47]. Copyright 2012, American Chemical Society).
function of graphene under a noble metal was measured to be ~4.62 eV, as shown in Fig. 3 c, and this is in good agreement with the previously reported values of the potential step at the graphene/ Pd interface (approximately 95~120 m e V) [39 - 41] . This indicates that the charge transfer between graphene and metal depends wholly on the species of metal. The measured work functions are, however, not consistent with previous theoretical results showing that the estimated work function of graphene under Au is ~4.74 eV. The Ni contact is also expected to alter the electrical characteristics of graphene [34] . As similar discord between experimental and calculation results has also been reported for contacts between carbon nanotubes and metal [48] , and inconsistencies among the calculations have been reported [49] , it is believed that the previous theoretical prediction may have some limitations when used to explain the behavior of the work function of graphene under a metal electrode.
3. Charge transport in contact
Carrier transport in a graphene-metal contact system is understood to take place in two cascading steps. Charges in the metal are first injected into the underlying graphene and then flow towards the graphene channel [21] . The former step is significantly affected by the contact metal species and the metallization process. The interaction between the graphene and metal, the work function of the graphene, or residues at the interface are crucial factors determining the transmission of charges at the interface. The injected carriers then face a built-in potential barrier, the charge transfer region, formed by the pinning of the work function in graphene due to the metal contact as noted in the previous section. The interfacial potential barrier width, shape, and
PPT Slide
Lager Image
Extracted contact resistances with the back gate voltage (Rc-Vg) with various metals. The inset shows a microscopic image of a transmission length measurement device (Reprinted with permission from earlier work [47]. Copyright 2012, American Chemical Society).
height can crucially affect the carrier tunneling probability [21] . In addition to the two interfaces, the conduction modes of both the graphene at the channel and the contact region contribute to the contact resistance. Therefore, the lowest number of modes between them and the carrier transmission probability at the two interfaces determine the overall contact conductance.
As the conduction mode of graphene contributes to the properties of the contact, the contact resistance shows gate-dependent behavior, as shown in Fig. 4 . Similar to the bell-shaped resistance of graphene as a function of the gate voltage, the contact resistance also reaches its maximum value near the Dirac voltage due to the absence of a conduction mode in the channel. In contrast to the channel, it is not clearly understood if the conduction mode of graphene underneath a metal contact can be modulated by back gate bias. The number of peaks in the contact resistance (or the channel resistance) offer a clue about the modulation of the density-of-state (DOS) in graphene at the contact region, as shown in Figs. 2 c and f [21 , 41] . In a simple situation without charge puddles in the channel and the broadening of DOS in the contact region, if graphene at the contact and channel region can be modulated by the back gate bias and if each has a different energy difference between the Fermi level and Dirac point energy, the contact resistance will exhibit two peaks. One of them originates from the minimum DOS in the channel and the other peak corresponds to the minimum DOS in the graphene at the contact. The single peak in the contact resistance measured as a function of the back gate voltage indicates that graphene underneath metal cannot be modulated by the bias. Most reports on the contact resistance measured with the back gate voltage have shown a single peak, indicating that the metal-induced pinning of the work function of graphene is independent of the external field.
The determination of a single peak in the contact resistance is not easy in real situations. One of the reasons is the large standard deviation of the contact resistance near the Dirac voltage, leading to an unclear position of the peak or an ambiguous situation when attempting to count the number of peaks. Such larger error bars originates from the charge puddles in graphene at the channel or the contact, or for instance by the extraction method. More details about the larger distribution of the contact resistance are given in Section 4. The larger broadening energy in the contact region also makes the peaks vague and makes it difficult to identify the presence of peaks [21] . In addition, the broadening of DOS in the graphene channel results in the broadening of the peak, making the separation of the peaks difficult. There are, however, few reports on the modulation of DOS in graphene underneath metal. The asymmetric transfer characteristics of fully covered top-gated graphene FETs have been explained in terms of the impact of the back-gate bias on the Fermi level shift in the graphene under the metal [41 , 46 , 50] . Chen and Appenzeller [46] also suggested that the energy dispersion of graphene underneath metal can be altered due to the much lower dependence on the back-gate bias. The similar broadening of the DOS of graphene at the contact area has been introduced to explain the observation of a single minimum during asymmetric conductance [51] . From simulation results in the referenced study, it was estimated that the transfer characteristics will have two minima when the energy dispersion of graphene in the contact region is intrinsic and the Fermi level is not aligned with that of the channel graphene, as two the Dirac levels in both the channel and contact region can be modulated together with back-gate bias. Interestingly, two minima have been measured from Ni-contacted FETs, indicating the free of metal-induced pinning, and it has been suggested that the oxidation of the metal surface occurs at the contact interface [44] . It was noted that the amount of the modulation of DOS in the graphene at the contact region can differ from that in the channel due to the different coupling ratio contributed by the metal [50] . As the coupling ratio at the contact is wholly determined by the graphene-metal interaction, the distance between the graphene and the metal, the metal-induced broadening, and the inhomogeneous broadening may determine the amount of modulation at the contact region [21] . Because the influence of the back-gate voltage on the work function pinning of graphene underneath the metal can play a critical role in the operation of the device, clearer and direct evidence must be identified in the future.
At a high negative or positive back-gate bias with respect to the Dirac voltage, the conduction mode of the channel is large enough and the contact resistance is determined by the charge transmittance at the two interfaces and the DOS of the graphene at the contact point. A higher conduction mode or higher work function difference as regards the Fermi level of the graphene underneath the metal is expected to show lower contact resistance [21] . It has been found, however, that a higher work function of graphene does not guarantee lower contact resistance [47] . Ni contact results in the work function of graphene exceeding 5 eV, but the contact resistance is about 10 times higher compared to a Pd contact in a high-bias regime. This indicates that the work function of graphene underneath metal is not a critical factor affecting the contact resistance.
One of the possible reasons for the non-critical contribution of the work function may be the current crowding effect [25 , 52] . During graphene-metal contact, the charges flow within a finite length at the edge and the contact resistance is a linear function of the width of the metal, but not a function of the length of the metal on the graphene [25] . The transfer length, which is the length at the potential decreases to e -1 from the edge, as shown in Fig. 1 e, is known to be shorter than 1 μm [25 , 31] . The transfer length with a Pd contact has been found to range from 50 to 250 nm [21] , while that with Ti contact reaches only about 200 nm [53] . It was noted that such a short Lt leads to current crowding at the contact edge, which generates local heating and results in degraded contact performance or undesirable electrical properties such as an increase in 1/f noise [52 , 54] .
The carrier transport at the p-p’ and n-n’ interfacial built-in potential is known to be nearly transparent due to Klein tunneling [21 , 55 , 56] . In contrast, transmission through the p-n junction undergoes relatively high resistance, as explained previous section. The charge transfer through the interfacial barrier in graphene at the contact edge has been analyzed with various potential profiles [21] . In addition to exponential potential profiles, linear, 1/ x and x -1/2 characteristics have been evaluated in these calculations, with the result showing that the potential shape does not have a significant impact on the charge transmission probability. The width of the barrier has also been investigated, and the contact resistance was found to be insensitive to the width as long as the width is larger than 20 nm.
Carrier transport at metal-graphene junctions was thoroughly investigated in an earlier study [21] . The injection efficiency at the interface has been classified into two regimes: one characterized by the diffusive transport limit and the other by the ballistic transport limit. Regarding the ballistic transport limit, where the scattering mean-free path is longer than the graphene-metal coupling length, the injection efficiency is expected to be close to one and the Pd contact is estimated to exhibit transmission efficiency of around 75%, indicating the ballistic transport type. The graphene-metal coupling length can be strongly affected by the interaction between the π-orbitals in graphene and the metal d-states. Therefore, the graphene-metal distance, wettability, or the deformation potential due to the metal contact can play significant roles in the coupling length. In terms of the graphene- metal distance, Ti is expected to show the lowest contact resistance as compared to Pd, Pt, Cu, and Au [57] . On the other hand, utilizing higher quality graphene, which exhibits higher mobility and less scattering, leads to the longer mean-free path. In addition, the scattering mean-free path can also be affected by phonon interaction and by coulomb scattering with the supporting substrate. Among the various treatments on SiO 2 , the O 2 plasma treatment on SiO 2 results in stronger interaction with graphene compared to Ar or H 2 gas treatments, leading to degraded mobility. It has been reported that HF-dipped SiO 2 shows smaller hysteresis and higher mobility compared to other treatment processes such as that with O 2 plasma or re-oxidation. Different substrates also cause different doping states in graphene [58 , 59] . These results indicate that the interaction between graphene and substrate can strongly affect the contact properties and the graphene-metal interaction.
4. Process dependence on the contact resistance
At an early stage of study on contact resistance in graphene FETs, no gate voltage dependence was thought to exist [29 , 30] . However, strong dependence on the gate bias was noted, sometimes even showing negative values [19 , 21 , 26 , 43 , 47] . This inconsistency may stem from large errors in the extraction method or insufficient data sets [26 , 30] . Although the TLM and fourprobe measurement methods are the common ways to extract the contact resistance, they can contain errors originating from the basic assumptions of the methods.
The TLM method is based on the assumption of homogenous properties of the graphene at both the channel and contact regions, such as the Dirac voltage (V Dirac ), the sheet resistance with respect to the normalized voltage (V g -V Dirac ), and the contribution of the metal or the surrounding insulators [43] . These ideal conditions are not possible in the real situations, and errors mostly likely arise. This problem is especially more severe with two terminal devices. The dominant origins of the inhomogeneous properties are the charge puddles in graphene [60 , 61] and the presence of a charge transfer region due to the metal contact, as explained in the previous section. Charge puddles and/ or chemical doping can be induced by the absorption of ambient molecules or charge transfer from the substrate or residual materials such as the photoresist. Such inhomogeneity can lead to the broadening of the Dirac voltage and result in the suppression of the peak resistance at the Dirac voltage. Many have reported that the contact resistance measured as a function of the gate voltage contains relatively large error bars near the Dirac voltage, as shown in Fig. 4 [21 , 47] . The charge transfer region induced by the metal contact can play a significant role when the dimensions of the devices are compatible with the carrier mean free path [26] . If doping extends far into the graphene in short channel devices, it can lead to inconsistent sheet resistance in the channel region, provoking errors in measurement. Blake et al. clearly pointed out this issue, reporting that a macroscopic inhomogeneity induced by either chemical doping or metal contacts is an important factor when evaluating the properties of graphene near the Dirac voltage [26] . Fig. 5 illustrates the possible reason behind the errors in the extraction method for the four-point probe method and the TLM method [43] . The inhomogeneous sheet resistance can even lead to negative contact resistance.
There have been many reports on charge puddles or substrateinduced doping [61 - 65] , and it can be suppressed by choosing the optimum dielectric materials or by fabrication devices with well-developed process conditions [66 - 69] . It has been found that the intrinsic doping of graphene and low hysteresis are obtained by depositing polymers on top of the substrate, such as a hexamethyldisilazane or phenylsilane self-assembled monolayer [24 , 66] . The utilization of boron-nitride film as a substrate has also been shown to preserve the excellent properties of graphene and significantly suppress the charge inhomogeneity [70 - 73] . A photoresist layer or chemicals which remain on the graphene during the device fabrication process seriously affect the inhomogeneity of the graphene. To reduce the contamination on the contact region, a plasma treatment prior to metal deposition during the lift-off process has been introduced [74 , 75] .
PPT Slide
Lager Image
(a) A possible explanation for the extraction errors from the four–point probe measurement. The contact resistance is obtained from (Rtotal-Rch S·L/l)/2, where Rtotal is the total resistance between the source and drain, Rch S is the channel resistance between the two voltage probes, L is the length between the source and the drain, and l the length between the two voltage probes. When the Dirac voltages are not coincident due to inhomogeneous doping in the channel, negative values or undesired results may arise. Larger error bars near the Dirac voltage normally originate for this reason, but a value at a higher or a lower gate voltage is close to the actual value (Reprinted with permission from earlier research [43]. Copyright 2011, The Japan Society of Applied Physics). (b) In the transmission line model, the inhomogeneous also could lead the extraction errors. Although the data fit in the linear relation in an ideal situation, higher or lower contact resistance can result when the sheet resistances are not consistent with the length between the source and the drain (Lgr), as shown by red-triangle dots or the blue square dots.
The low-power plasma in such a process not only removes resist residue on graphene but also creates defects in the contact region of graphene, thereby enhancing the metal-graphene interaction. An O 2 plasma treatment process improves the contact resistance by nearly 6000 times compared to a conventional process [74] . It has also been revealed that this technique results in similar levels of contact resistivity regardless of the metal species. In addition to O 2 , various plasma types such as Ar, N 2 , and H 2 have been exposed to graphene, though O 2 plasma has been found to be more effective in making graphene chemically reactive, thereby enhancing the adhesion with the metal [75] , It was noted that the metal deposition process also can affect the performance of graphene devices. It has been found that sputtering Ti on graphene greatly degrades the device performance compared to the e-beam evaporation method [76] . Although the sputtering method also generates defects in graphene like the O 2 plasma treatment, the situation is somewhat different. Nonetheless, the exact mechanism causing this difference has not been well identified.
5. Contact module process
At the current stage of graphene technology, the immature process technology is one of the critical obstacles to overcome in order to obtain good contact. The preparation of graphene, the metal deposition process, the substrate condition, or even the measurement environment can all have an impact on the measured contact properties. This can be noticed from the face, where there are widely varied or inconsistent results regarding the contact resistance in previous reports. For example, the Ti contact was predicted to show poor results in one report [25] , but in other reports, it results in superior contact properties over other metals [22 , 57] . The conventional facial contact structure fabricated by the lift-off process is also a limitation when seeking to obtain good contact performance. The transport of carriers in the vertical direction of graphene is known to be different from the nature of the transport on the basal plane. The superior electric properties of graphene are attributed to the π-orbital electrons on the basal plane, while the conduction of the carriers along the c-axis basically originates from the thermal excitation or from impurity-assisted interlayer hopping [77 , 78] . The carrier transport in the conventional structure with which the metal is deposited on the graphene basal plane is therefore expected to limit how much the performance can be enhanced. To improve the transmission probability of carriers, the dual gate structure which graphene is sandwiched with top and buried metals has been introduced and it results improved contact resistance of 40% compared to the conventional top-contacted structure [79] . Although it is a new approach, the structure also explores the facial contact. An alternative approach to overcome the facial contact is to use the edge contact structure which is formed by a covalent bond between the metal atoms and the graphene; it is expected to show significantly lower contact resistance compared to the facial contact. Theoretical reports estimated that the contact resistivity could be as low as 10 to 10 3 times the facial contact based on a quantum tunneling process [57] . Very recently, contact modules considering the edge contact resistance were developed [80] . The proposed contact module results in an improvement of approximately 20~60% compared to the conventional method involving a lift-off process.
6. Summary
Graphene-metal contact regions exhibit unique characteristics compared to the conventional semiconductor-metal contact. Graphene can be doped by the contacting metal and the work function strongly depends on the metal species. Work function pinning can form a p-n junction in the graphene channel, leading to asymmetrical conductance behavior and contact resistance. The contact resistance is significantly affected by the density of states of the graphene in the contact and channel regions, engineercausing the gate voltage to modulate the contact resistance. The charge transport at the graphene-metal interface and at the edge of the metal contact determines the intrinsic contact resistance, which must be minimized. Although the exact mechanism has not been suitably explained thus far, the potential profile at the contact edge does not seem to contribute to the contact resistance significantly. As the scattering mean-free path of graphene at the contact point can affect the intrinsic contact resistance, the influence of the substrate can also play a major role in addition to the quality of the graphene. In order to achieve good contact properties, associated process technologies must mature as state-of-the-art CMOS technology. An advanced contact process scheme must also be developed, and the edge contact structure may be a solution in this area.
This work was supported by research grants from the NationalResearch Foundation of Korea (NRF) (2008-2002744, 2010-0029132, and 2011-0031638).
Lemme MC , Echtermeyer TJ , Baus M , Kurz H (2007) A graphene fieldeffectdevice. IEEE Electron Device Lett 28 282 -    DOI : 10.1109/LED.2007.891668
Meric I , Baklitskaya N , Kim P , Shepard KL (2008) RF performance oftop-gated, zero-bandgap graphene field-effect transistors. IEEE International Electron Devices Meeting San Francisco, CA 1 -    DOI : 10.1109/IEDM.2008.4796738
Lin YM , Jenkins KA , Valdes-Garcia A , Small JP , Farmer DB , Avouris P (2009) Operation of graphene transistors at gigahertz frequencies. Nano Lett 9 422 -    DOI : 10.1021/Nl803316h
Lin YM , Jenkins K , Farmer D , Valdes-Garcia A , Avouris P , Sung CY , Chiu HY , Ek B (2009) Development of graphene FETs forhigh frequency electronics. IEEE International Electron Devices Meeting Baltimore, MD 1 -    DOI : 10.1109/IEDM.2009.5424378
Farmer DB , Chiu HY , Lin YM , Jenkins KA , Xia FN , Avouris P (2009) Utilization of a buffered dielectric to achieve high field-effect carriermobility in graphene transistors. Nano Lett 9 4474 -    DOI : 10.1021/Nl902788u
Dimitrakopoulos C , Lin YM , Grill A , Farmer DB , Freitag M , Sun YN , Han SJ , Chen ZH , Jenkins KA , Zhu Y , Liu ZH , McArdle TJ , Ott JA , Wisnieff R , Avouris P (2010) Wafer-scale epitaxial graphene growth on the Si-face of hexagonal SiC (0001) for high frequency transistors. J Vac Sci Technol, B 28 985 -    DOI : 10.1116/1.3480961
Lin YM , Chiu HY , Jenkins KA , Farmer DB , Avouris P , Valdes-Garcia A (2010) Dual-gate graphene FETs with f(T) of 50 GHz. IEEE Electron Device Lett 31 68 -    DOI : 10.1109/LED.2009.2034876
Lin YM , Dimitrakopoulos C , Jenkins KA , Farmer DB , Chiu HY , Grill A , Avouris P (2010) 100-GHz transistors from wafer-scale epitaxial graphene. Science 327 662 -    DOI : 10.1126/science
Pince E , Kocabas C (2010) Investigation of high frequency performancelimit of graphene field effect transistors. Appl Phys Lett 97 173106 -    DOI : 10.1063/1.3506506
Liao L , Lin YC , Bao M , Cheng R , Bai J , Liu Y , Qu Y , Wang KL , Huang Y , Duan X (2010) High-speed graphene transistors with a selfaligned nanowire gate. Nature 467 305 -    DOI : 10.1038/nature09405
Chauhan J , Guo J (2011) Assessment of high-frequency performance limitsof graphene field-effect transistors. Nano Res 4 571 -    DOI : 10.1007/s12274-011-0113-1
Das S , Appenzeller J (2011) An all-graphene radio frequency lownoise amplifier. IEEE Radio Frequency Integrated Circuits Symposium Baltimore, MD 1 -    DOI : 10.1109/RFIC.2011.5940628
Koswatta SO , Valdes-Garcia A , Steiner MB , Lin YM , Avouris P (2011) Ultimate RF potential of carbon electronics. IEEE Trans Microwave Theory Tech 59 2739 -    DOI : 10.1109/tmtt.2011.2150241
Moon JS , Curtis D , Zehnder D , Kim S , Gaskill DK , Jernigan GG , Myers-Ward RL , Eddy CR , Campbell PM , Lee KM , Asbeck P (2011) Low-phase-noise graphene FETs in ambipolar RF applications. IEEE Electron Device Lett 32 270 -    DOI : 10.1109/led.2010.2100074
Wu Y , Lin Y , Bol AA , Jenkins KA , Xia F , Farmer DB , Zhu Y , Avouris P (2011) High-frequency, scaled graphene transistors ondiamond-like carbon. Nature 472 74 -    DOI : 10.1038/nature09979
Badmaev A , Che YC , Li Z , Wang C , Zhou CW (2012) Self-aligned fabricationof graphene RF transistors with T-shaped gate. ACS Nano 6 3371 -    DOI : 10.1021/nn300393c
Cheng R , Bai JW , Liao L , Zhou HL , Chen Y , Liu LX , Lin YC , Jiang S , Huang Y , Duan XF (2012) High-frequency self-aligned graphene transistors with transferred gate stacks. Proc Natl Acad Sci U S A 109 11588 -    DOI : 10.1073/pnas.1205696109
Wu YQ , Jenkins KA , Valdes-Garcia A , Farmer DB , Zhu Y , Bol AA , Dimitrakopoulos C , Zhu WJ , Xia FN , Avouris P , Lin YM (2012) State-of-the-art graphene high-frequency electronics. Nano Lett 12 3062 -
Kim S , Nah J , Jo I , Shahrjerdi D , Colombo L , Yao Z , Tutuc E , Banerjee SK (2009) Realization of a high mobility dual-gated graphenefield-effect transistor with Al2O3dielectric. Appl Phys Lett 94 062107 -    DOI : 10.1063/1.3077021
Shin WC , Kim TY , Sul O , Choa BJ (2012) Seeding atomic layer depositionof high-k dielectric on graphene with ultrathin poly (4-vinylphenol)layer for enhanced device performance and reliability. Appl Phys Lett 101 033507 -    DOI : 10.1063/1.4737645
Xia F , Perebeinos V , Lin Y , Wu Y , Avouris P (2011) The origins and limitsof metal-graphene junction resistance. Nat Nanotechnol 6 179 -    DOI : 10.1038/nnano.2011.6
Moon JS , Antcliffe M , Seo HC , Curtis D , Lin S , Schmitz A , Milosavljevic I , Kiselev AA , Ross RS , Gaskill DK , Campbell PM , Fitch RC , Lee KM , Asbeck P (2012) Ultra-low resistance ohmic contacts in graphene field effect transistors. Appl Phys Lett 100 203512 -    DOI : 10.1063/1.4719579
Farmer DB , Lin YM , Avouris P (2010) Graphene field-effect transistorswith self-aligned gates. Appl Phys Lett 97 013103 -    DOI : 10.1063/1.3459972
Liu Z , Bol AA , Haensch W (2010) Large-scale graphene transistors withenhanced performance and reliability based on interface engineering by phenylsilane self-assembled monolayers. Nano Lett 11 523 -    DOI : 10.1021/nl1033842
Nagashio K , Nishimura T , Kita K , Toriumi A (2009) Metal/graphene contactas a performance Killer of ultra-high mobility graphene analysisof intrinsic mobility and contact resistance. IEEE International Electron Devices Meeting Baltimore, MD 1 -    DOI : 10.1109/IEDM.2009.5424297
Blake P , Yang R , Morozov S , Schedin F , Ponomarenko L , Zhukov A , Nair R , Grigorieva I , Novoselov K , Geim A (2009) Influence of metal contacts and charge inhomogeneity on transport properties of graphene near the neutrality point. Solid State Commun 149 1068 -    DOI : 10.1016/j.ssc.2009.02.039
Murali R , Yang Y , Brenner K , Beck T , Meindl JD (2009) Breakdown currentdensity of graphene nanoribbons. Appl Phys Lett 94 243114 -    DOI : 10.1063/1.3147183
Xia F , Farmer DB , Lin Y , Avouris P (2010) Graphene field-effect transistorswith high on/off current ratio and large transport band gapat room temperature. Nano Lett 10 715 -    DOI : 10.1021/nl9039636
Russo S , Craciun M , Yamamoto M , Morpurgo A , Tarucha S (2010) Contact resistance in graphene-based devices. Physica E 42 677 -    DOI : 10.1016/j.physe.2009.11.080
Venugopal A , Colombo L , Vogel E (2010) Contact resistance in few andmultilayer graphene devices. Appl Phys Lett 96 013512 -    DOI : 10.1063/1.3290248
Nagashio K , Nishimura T , Kita K , Toriumi A (2010) Contact resistivityand current flow path at metal/graphene contact. Appl Phys Lett 97 143514 -    DOI : 10.1063/1.3491804
Schwierz F (2010) Graphene transistors. Nat Nanotechnol 5 487 -    DOI : 10.1038/nnano.2010.89
Leonard F , Talin AA (2011) Electrical contacts to one-and two-dimensionalnanomaterials. Nat Nanotechnol 6 773 -    DOI : 10.1038/nnano.2011.196
Giovannetti G , Khomyakov P , Brocks G , Karpan V , Van den Brink J , Kelly P (2008) Doping graphene with metal contacts. Phys Rev Lett 101 26803 -    DOI : 10.1103/PhysRevLett.101.026803
Khomyakov P , Starikov A , Brocks G , Kelly P (2010) Nonlinear screeningof charges induced in graphene by metal contacts. Phys Rev B 82 115437 -    DOI : 10.1103/PhysRevB.82.115437
Yu YJ , Zhao Y , Ryu S , Brus LE , Kim KS , Kim P (2009) Tuning the graphenework function by electric field effect. Nano Lett 9 3430 -    DOI : 10.1021/nl901572a
Yan L , Punckt C , Aksay IA , Mertin W , Bacher G (2011) Local voltagedrop in a single functionalized graphene sheet characterized byKelvin probe force microscopy. Nano Lett 11 3543 -    DOI : 10.1021/nl201070c
Lee EJH , Balasubramanian K , Weitz RT , Burghard M , Kern K (2008) Contact and edge effects in graphene devices. Nat Nanotechnol 3 486 -    DOI : 10.1038/nnano.2008.172
Xia F , Mueller T , Golizadeh-Mojarad R , Freitag M , Lin Y , Tsang J , Perebeinos V , Avouris P (2009) Photocurrent imaging and efficient photon detection in a graphene transistor. Nano Lett 9 1039 -    DOI : 10.1021/nl8033812
Mueller T , Xia F , Freitag M , Tsang J , Avouris P (2009) Role of contactsin graphene transistors: A scanning photocurrent study. Phys Rev B 79 245430 -    DOI : 10.1103/PhysRevB.79.245430
Knoch J , Chen Z , Appenzeller J (2011) Properties of metal-graphene contacts. IEEE Trans Nanotechnol 11 513 -    DOI : 10.1109/TNANO.2011.2178611
Low T , Hong S , Appenzeller J , Datta S , Lundstrom MS (2009) Conductanceasymmetry of graphene pn junction. IEEE Trans Electron Devices 56 1292 -    DOI : 10.1109/TED.2009.2017646
Nagashio K , Toriumi A (2011) Density-of-states limited contact resistancein graphene field-effect transistors. Jpn J Appl Phys 50 070108 -    DOI : 10.1143/jjap.50.070108
Nouchi R , Tanigaki K (2010) Charge-density depinning at metal contactsof graphene field-effect transistors. Appl Phys Lett 96 253503 -    DOI : 10.1063/1.3456383
Huard B , Stander N , Sulpizio J , Goldhaber-Gordon D (2008) Evidenceof the role of contacts on the observed electron-hole asymmetryin graphene. Phys Rev B 78 121402 -    DOI : 10.1103/PhysRevB.78.121402
Chen Z , Appenzeller J (2009) Gate modulation of graphene contacts-onthe scaling of graphene FETs. Symposium on VLSI Technology Honolulu, HI 128 -
Song SM , Park JK , Sul OJ , Cho BJ (2012) Determination of work functionof graphene under a metal electrode and its role in contactresistance. Nano Lett 12 3887 -    DOI : 10.1021/nl300266p
Wang QJ , Che JG (2009) Origins of distinctly different behaviors of Pdand Pt contacts on graphene. Phys Rev Lett 103 66802 -    DOI : 10.1103/PhysRevLett.103.066802
Ran Q , Gao M , Guan X , Wang Y , Yu Z (2009) First-principles investigationon bonding formation and electronic structure of metal-graphenecontacts. Appl Phys Lett 94 103511 -    DOI : 10.1063/1.3095438
Berdebes D , Low T , Sui Y , Appenzeller J , Lundstrom MS (2011) Substrategating of contact resistance in graphene transistors. IEEE Trans Electron Devices 58 3925 -    DOI : 10.1109/TED.2011.2163800
Farmer DB , Golizadeh-Mojarad R , Perebeinos V , Lin YM , Tulevski GS , Tsang JC , Avouris P (2008) Chemical doping and electron– hole conduction asymmetry in graphene devices. Nano Lett 9 388 -    DOI : 10.1021/nl803214a
Grosse KL , Bae MH , Lian F , Pop E , King WP (2011) Nanoscale Jouleheating, Peltier cooling and current crowding at graphenemetalcontacts. Nat Nanotechnol 6 287 -    DOI : 10.1038/nnano.2011.39
Xu HT , Wang S , Zhang ZY , Wang ZX , Xu HL , Peng LM (2012) Contactlength scaling in graphene field-effect transistors. Appl Phys Lett 100 103501 -    DOI : 10.1063/1.3691629
Murrmann H , Widmann D (1969) Current crowding on metal contactsto planar devices. IEEE Trans Electron Devices 16 1022 -    DOI : 10.1109/T-ED.1969.16904
Cheianov VV , Fal’ko VI (2006) Selective transmission of Dirac electronsand ballistic magnetoresistance of n-p junctions in graphene. Phys Rev B 74 041403 -    DOI : 10.1103/PhysRevB.74.041403
Katsnelson MI , Novoselov KS , Geim AK (2006) Chiral tunnelling andthe Klein paradox in graphene. Nat Phys 2 620 -    DOI : 10.1038/Nphys384
Matsuda Y , Deng WQ , Goddard WA (2010) Contact resistance for “endcontacted”metal− graphene and metal− nanotube interfaces fromquantum mechanics. J Phys Chem C 114 17845 -    DOI : 10.1021/jp806437y
Song SM , Cho BJ (2010) Investigation of interaction between grapheneand dielectrics. Nanotechnology 21 335706 -    DOI : 10.1088/0957-4484/21/33/335706.
Oh JG , Shin YS , Shin WC , Sul OJ , Cho BJ (2011) Dirac voltage tunabilityby Hf1-xLaxO gate dielectric composition modulation forgraphene field effect devices. Appl Phys Lett 99 193503 -    DOI : 10.1063/1.3610476
Martin J , Akerman N , Ulbricht G , Lohmann T , Smet J , Von Klitzing K , Yacoby A (2007) Observation of electron–hole puddles in grapheneusing a scanning single-electron transistor. Nat Phys 4 144 -    DOI : 10.1038/nphys781
Zhang Y , Brar VW , Girit C , Zettl A , Crommie MF (2009) Origin of spatialcharge inhomogeneity in graphene. Nat Phys 5 722 -    DOI : 10.1038/nphys1365
Liu H , Liu Y , Zhu D (2011) Chemical doping of graphene. J Mater Chem 21 3335 -    DOI : 10.1039/C0JM02922J
Levesque PL , Sabri SS , Aguirre CM , Guillemette J , Siaj M , Desjardins P , Szkopek T , Martel R (2010) Probing charge transfer at surfacesusing graphene transistors. Nano Lett 11 132 -    DOI : 10.1021/nl103015w
Casiraghi C , Pisana S , Novoselov K , Geim A , Ferrari A (2007) Ramanfingerprint of charged impurities in graphene. Appl Phys Lett 91 233108 -    DOI : 10.1063/1.2818692
Berciaud S , Ryu S , Brus LE , Heinz TF (2008) Probing the intrinsic propertiesof exfoliated graphene: Raman spectroscopy of free-standingmonolayers. Nano Lett 9 346 -    DOI : 10.1021/nl8031444
Lafkioti M , Krauss B , Lohmann T , Zschieschang U , Klauk H , Klitzing K , Smet JH (2010) Graphene on a hydrophobic substrate: dopingreduction and hysteresis suppression under ambient conditions. Nano Lett 10 1149 -    DOI : 10.1021/nl903162a
Pirkle A , Chan J , Venugopal A , Hinojos D , Magnuson C , Mc-Donnell S , Colombo L , Vogel E , Ruoff R , Wallace R (2011) The effectof chemical residues on the physical and electrical properties ofchemical vapor deposited graphene transferred to SiO2. Appl Phys Lett 99 122108 -    DOI : 10.1063/1.3643444
Lin YC , Lu CC , Yeh CH , Jin C , Suenaga K , Chiu PW (2011) Grapheneannealing: how clean can it be? Nano Lett 12 414 -    DOI : 10.1021/nl203733r
Cheng Z , Zhou Q , Wang C , Li Q , Wang C , Fang Y (2011) Toward intrinsicgraphene surfaces: a systematic study on thermal annealing andwet-chemical treatment of SiO2-supported graphene devices. Nano Lett 11 767 -    DOI : 10.1021/nl103977d
Dean CR , Young AF , Meric I , Lee C , Wang L , Sorgenfrei S , Watanabe K , Taniguchi T , Kim P , Shepard KL , Hone J (2010) Boron nitridesubstrates for high-quality graphene electronics. Nat Nanotechnol 5 722 -    DOI : 10.1038/nnano.2010.172
Xue J , Sanchez-Yamagishi J , Bulmash D , Jacquod P , Deshpande A , Watanabe K , Taniguchi T , Jarillo-Herrero P , LeRoy BJ. (2011) Scanningtunnelling microscopy and spectroscopy of ultra-flat graphene onhexagonal boron nitride. Nat Mater 10 282 -    DOI : 10.1038/nmat2968.
Decker R , Wang Y , Brar VW , Regan W , Tsai H-Z , Wu Q , Gannett W , Zettl A , Crommie MF. (2011) Local electronic properties of grapheneon a BN substrate via scanning tunneling microscopy. Nano Lett 11 2291 -    DOI : 10.1021/nl2005115
Kim K , Choi JY , Kim T , Cho SH , Chung HJ. (2011) A role for graphenein silicon-based semiconductor devices. Nature 479 338 -    DOI : 10.1038/nature10680
Robinson JA , LaBella M , Zhu M , Hollander M , Kasarda R , Hughes Z , Trumbull K , Cavalero R , Snyder D. (2011) Contacting graphene. Appl Phys Lett 98 053103 -    DOI : 10.1063/1.3549183
Choi MS , Lee SH , Yoo WJ. (2011) Plasma treatments to improve metal contacts in graphene field effect transistor. J Appl Phys 110 073305 -    DOI : 10.1063/1.3646506
Liu W , Li M , Xu S , Zhang Q , Zhu Y , Pey K , Hu H , Shen Z , Zou X , Wang J. (2010) Understanding the contact characteristics in singleor multi-layer graphene devices: the impact of defects (carbonvacancies) and the asymmetric transportation behavior. IEEE International Electron Devices Meeting San Francisco, CA 23.3.1    DOI : 10.1109/IEDM.2010.5703420
Matsubara K , Sugihara K , Tsuzuku T. (1990) Electrical-resistance in theC-direction of graphite. Phys Rev B 41 969 -    DOI : 10.1103/PhysRevB.41.969
Khatami Y , Li H , Xu C , Banerjee K. (2012) Metal-to-multilayer-graphenecontact—Part I: Contact resistance modeling. IEEE Trans Electron Devices 59 2444 -    DOI : 10.1109/TED.2012.2205256.
Franklin AD , Han SJ , Bol AA , Perebeinos V. (2012) Double Contactsfor Improved Performance of Graphene Transistors. IEEE Electron Device Lett 33 17 -    DOI : 10.1109/Led.2011.2173154
Smith JT , Franklin AD , Farmer DB , Dimitrakopoulos CD. (2013) Reducingcontact resistance in graphene devices through contact areapatterning. ACS Nano 7 3661 -    DOI : 10.1021/nn400671z.