SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Online First articles in Volume 6 (2011). These articles were published as Online First on SpringerLink in 2010. They should be cited with a 2011 publication year.

Open Access Highly Accessed Open Badges Nano Express

Conductance of Graphene Nanoribbon Junctions and the Tight Binding Model

Y Wu and PA Childs*

Author Affiliations

School of Electronic, Electrical and Computer Engineering, University of Birmingham, B15 2TT, Birmingham, UK

For all author emails, please log on.

Nanoscale Res Lett 2011, 6:62  doi:10.1007/s11671-010-9791-y

The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/6/1/62

Received:9 July 2010
Accepted:9 September 2010
Published:7 October 2010

© 2010 Wu and Childs.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Planar carbon-based electronic devices, including metal/semiconductor junctions, transistors and interconnects, can now be formed from patterned sheets of graphene. Most simulations of charge transport within graphene-based electronic devices assume an energy band structure based on a nearest-neighbour tight binding analysis. In this paper, the energy band structure and conductance of graphene nanoribbons and metal/semiconductor junctions are obtained using a third nearest-neighbour tight binding analysis in conjunction with an efficient nonequilibrium Green's function formalism. We find significant differences in both the energy band structure and conductance obtained with the two approximations.

Graphene nanoribbon junction; Tight binding; Conductance; Band structure


Since the report of the preparation of graphene by Novoselov et al. [1] in 2004, there has been an enormous and rapid growth in interest in the material. Of all the allotropes of carbon, graphene is of particular interest to the semiconductor industry as it is compatible with planar technology. Although graphene is metallic, it can be tailored to form semiconducting nanoribbons, junctions and circuits by lithographic techniques. Simulations of charge transport within devices based on this new technology exploit established techniques for low dimensional structures [2,3]. The current flowing through a semiconducting nanoribbon formed between two metallic contacts has been established using a nonequilibrium Green's Function (NEGF) formalism based coupled with an energy band structure derived using a tight binding Hamiltonian [4-7]. To minimise computation time, the nearest-neighbour tight binding approximation is commonly used to determine the energy states and overlap is ignored. This assumption has also been used for calculating the energy states of other carbon-based materials such as carbon nanotubes [8] and carbon nanocones [9]. Recently, Reich et al. [10] have demonstrated that this approximation is only valid close to the K points, and a tight binding approach including up to third nearest-neighbours gives a better approximation to the energy dispersion over the entire Brillouin zone.

In this paper, we simulate charge transport in a graphene nanoribbon and a nanoribbon junction using a NEGF based on a third nearest-neighbour tight binding energy dispersion. For transport studies in nanoribbons and junctions, the formulation of the problem differs from that required for bulk graphene. Third nearest-neighbour interactions introduce additional exchange and overlap integrals significantly modifying the Green's function. Calculation of device characteristics is facilitated by the inclusion of a Sancho-Rubio [11] iterative scheme, modified by the inclusion of third nearest-neighbour interactions, for the calculation of the self-energies. We find that the conductance is significantly altered compared with that obtained based on the nearest-neighbour tight binding dispersion even in an isolated nanoribbon. Hong et al. [12] observed that the conductance is modified (increased as well as decreased) by the presence of defects within the lattice. Our results show that details of the band structure can significantly modify the observed conductivities when defects are included in the structure.


The basis for our analysis is the hexagonal graphene lattice shown in Figure 1. a1 and b1 are the principal vectors of the unit cell containing two carbon atoms belonging to the two sub-lattices. Atoms on the concentric circles of increasing radius correspond to the nearest-neighbours, second nearest-neighbours and third nearest-neighbours, respectively.

thumbnailFigure 1. Armchair-edge graphene metal (index N = 23)/semiconductor (index N = 13) junction. The rectangle shows the semiconductor unit cell, and the concentric circles of increasing radius show first, second and third nearest-neighbours, respectively.

Saito et al. [8] derived the dispersion relation below using a nearest-neighbour tight binding analysis including the overlap integral s0.

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M1">View MathML</a>


Here, f(k) = 3 + 2 cos k · a1 +2 cos k · b1 + 2 cos k · (a1 - b1) and the parameters, ε2p, γ0 and s0 are obtained by fitting to experimental results or ab initio calculations.

Most analyses of charge transport in graphene-based structures simplify the result further by ignoring s0. Reich et al. [10] derived the dispersion relation for graphene based on third nearest-neighbours. In this work, the energy band structure of a graphene nanoribbon including third nearest-neighbour interactions is obtained from the block Hamiltonian and overlap matrices given below for the unit cell defined by the rectangle in Figure 1.

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M2">View MathML</a>


For the nth row of the above equation, we have

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M3">View MathML</a>


Considering the energy dispersion in the direction of charge transport, the Bloch form of the wavefunction ensures that φn~eikn. Substitution of φn into the above equation yields the secular equation

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M4">View MathML</a>


In the case of first nearest approximation without orbital overlap, Sn,n-1 and Sn,n+1 are empty matrices. To facilitate comparison with published results, we use an armchair-edge with index [13]N = 13 as our model nanoribbon. In the paper by Reich, tight binding parameters were obtained by fitting the band structure to that obtained by ab initio calculations. Recently, Kundo [14] has reported a set of tight biding parameters based on fitting to a first principle calculation but more directly related to the physical quantities of interest. These parameters have been utilised in our calculation and are presented below for third nearest-neighbour interactions (Table 1).

Table 1. Tight binding parameters [14]

Figure 2 compares the energy band structure of the modelled armchair-edge graphene nanoribbon obtained from the first nearest-neighbour tight binding method with that obtained by including up to third nearest-neighbours. Agreement is reasonable close to the K point but significant discrepancies occur at higher energies.

thumbnailFigure 2. Energy band structure of an N = 13 armchair graphene nanoribbon, (a) obtained from the first nearest-neighbour tight binding method and (b) including third nearest-neighbours.

Conductance of Graphene Nanoribbons and Junctions

Conductance in graphene nanoribbons and metal/semiconductor junctions is determined using an efficient nonequilibrium Green's function formalism described by Li and Lu [15]. The retarded Green's function is given by

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M5">View MathML</a>


Here, E+ = E + iη and η is a small positive energy value (10-5 eV in this simulation) which circumvents the singular point of the matrix inversion [16]. H is a tight binding Hamiltonian matrix including up to third nearest-neighbours, and S is the overlap matrix. Open boundary conditions are included through the left and right self-energy matrix elements, σL.R. The self-energies are independently evaluated through an iterative scheme described by Sancho et al. [11], modified to include third nearest-neighbour interactions. Determination of the retarded Green's function through equation 5 is facilitated by the inclusion of the body of the device in the right-hand contact through the recursive scheme described in ref. [15]. We will now outline the numerical procedure for deriving the conductance with third nearest-neighbour interactions included. Figure 3 shows a schematic of the unit cell labelling used to formulate the Green's function.

thumbnailFigure 3. Schematic showing the unit cell labelling used to formulate the Green's function.

We calculate the surface retarded Green's functions of the left and right leads by

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M6">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M7">View MathML</a>


where θ and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M8">View MathML</a> are the appropriate transfer matrices calculated from the following iterative procedure.

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M9">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M10">View MathML</a>


where ti and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M11">View MathML</a> are defined by

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M12">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M13">View MathML</a>



<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M14">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M15">View MathML</a>


The process is repeated until <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M16">View MathML</a> with δ arbitrarily small. The nonzero elements of the self-energies <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M17">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M18">View MathML</a> can be then obtained by

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M19">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M20">View MathML</a>


The conductance is obtained from the relation

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M21">View MathML</a>


where the transmission coefficient is obtained from

<a onClick="popup('http://www.nanoscalereslett.com/content/6/1/62/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/6/1/62/mathml/M22">View MathML</a>


with ΓL,R = iL,R - (ΣL,R).

Figure 4a, b compares the conductance of a graphene armchair-edge nanoribbon of index N = 13 and metal/semiconductor junction formed with the nanoribbon assuming first and third nearest-neighbour interactions, respectively. For graphene nanoribbons, differences are observed in the step-like structure, reflecting differences in the calculated band structure. When only first nearest-neighbour interactions are considered, the conductance of the conduction and valence bands is always symmetrical as determined by the formulation of the energy dispersion relation, equation 1. In the case of graphene nanoribbons, the conductance within a few electron volts of the Fermi energy is symmetrical for both first and third nearest-neighbour interactions. However, it is notable that at higher energies, overlap integrals introduced by third nearest-neighbour interactions result in asymmetry between the conductance in the conduction and valence bands. For metal/semiconductor junctions, significant differences in conductivity occur even at low energies due to mismatches of the sub-bands. Asymmetry in the conduction and valence band conductance (absent for first nearest-neighbour interactions) is also apparent when third nearest-neighbour interactions are included in the Green's function. Differences are also seen when defects are incorporated within a metal/semiconductor junction, an interesting system explored by Hong et.al. [12]. In this work, vacancies are introduced in the lattice at the positions marked by the solid rectangle and triangle in Figure 1 and the conductance obtained in each case. Hong et al. derive a coupling term associated with differences in band structure. For third nearest-neighbour, the solution to the coupling strength must be derived numerically.

thumbnailFigure 4. Conductance vs Energy for the junction shown in Figure 1, a using first nearest-neighbour parameters and b using third nearest-neighbours parameters. Dotted lines are for N = 13 armchair nanoribbon, solid lines are for ideal metal/semiconductor junctions, dotdash lines and dash lines are for junctions with a single defect type A (triangle in Figure 1) and type B (rectangle in Figure 1) respectively.


In this paper, we have determined the energy band structure of graphene nanoribbons and conductance of nanoribbons and graphene metal/semiconductor junctions using a NEGF formalism based on the tight binding method approximated to first nearest-neighbour and third nearest-neighbour. Significant differences are observed, suggesting the commonly used first nearest-neighbour approximation may not be sufficiently accurate in some circumstances. The most notable differences are observed when defects are introduced in the metal/semiconductor junctions.


  1. Novoselov KS, et al.: Two-dimensional gas of massless Dirac fermions in graphene.

    Nature 2005, 438:197-200. PubMed Abstract | Publisher Full Text OpenURL

  2. Cresti A, et al.: Charge transport in disordered graphene-based low dimensional materials.

    Nano Res 2008, 1:361-394. Publisher Full Text OpenURL

  3. Neto A, et al.: The electronic properties of graphene.

    Rev Mod Phys 2009, 81:109-162. Publisher Full Text OpenURL

  4. Golizadeh-Mojarad R, et al.: Atomistic non-equilibrium Green's function simulations of Graphene nano-ribbons in the quantum hall regime.

    J Comput Electron 2008, 7:407-410. Publisher Full Text OpenURL

  5. Liang G, et al.: Ballistic graphene nanoribbon metal-oxide-semiconductor field-effect transistors: a full real-space quantum transport simulation.

    J Appl Phys 2007, 102:054307. Publisher Full Text OpenURL

  6. Zhao P, Guo J: Modeling edge effects in graphene nanoribbon field-effect transistors with real and mode space methods.

    J Appl Phys 2009, 105:4503. OpenURL

  7. Odili D, et al.: Modeling charge transport in graphene nanoribbons and carbon nanotubes using a Schrödinger-Poisson solver.

    J Appl Phys 2009, 106:024509. Publisher Full Text OpenURL

  8. Saito R, et al.: Physical Properties of Carbon Nanotubes. Imperial College Press, London; 1998. Publisher Full Text

  9. Chen J, et al.: Low-energy electronic states of carbon nanocones in an electric field.

    Nano Micro Lett 2010, 2:121-125. OpenURL

  10. Reich S, et al.: Tight-binding description of graphene.

    Phys Rev B 2002, 66:035412. Publisher Full Text OpenURL

  11. Sancho M, et al.: Quick iterative scheme for the calculation of transfer matrices: application to Mo (100).

    J Phys F: Metal Phys 1984, 14:1205-1215. Publisher Full Text OpenURL

  12. Hong S, et al.: Metal-semiconductor junction of graphene nanoribbons.

    Appl Phys Lett 2008, 92:083107. Publisher Full Text OpenURL

  13. Nakada K, et al.: Edge state in graphene ribbons: nanometer size effect and edge shape dependence.

    Phys Rev B 1996, 54:17954-17961. Publisher Full Text OpenURL

  14. Kundu R: Tight binding parameters for graphene.


    Arxiv preprint arXiv:0907.4264

  15. Li TC, Lu S-P: Quantum conductance of graphene nanoribbons with edge defects.

    Phys Rev B 2008, 77:085408. Publisher Full Text OpenURL

  16. Datta S: Quantum Transport: Atom to Transistor. Cambridge University Press, Cambridge; 2005. OpenURL