SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series 11th Trends in NanoTechnology International Conference (TNT2010).

This article is part of the series 13th Trends in NanoTechnology International Conference.

Open Access Nano Express

Electronic and transport properties in circular graphene structures with a pentagonal disclination

Esther Jódar1*, Antonio Pérez–Garrido1 and Fernando Rojas12

Author affiliations

1 Departamento de Física Aplicada, Antiguo Hospital de Marina, Campus Muralla del Mar, UPCT, Cartagena, Murcia 30202, Spain

2 Departamento de Física Teórica, Centro de Nanociencias y Nanotecnologías, UNAM, Apdo. 356, Ensenada, Baja California 22830, Mexico

For all author emails, please log on.

Citation and License

Nanoscale Research Letters 2013, 8:258  doi:10.1186/1556-276X-8-258

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

Received:16 November 2012
Accepted:24 April 2013
Published:29 May 2013

© 2013 Jódar et al.; licensee Springer.

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.


We investigate the electronic and transport properties of circular graphene structures (quantum dots) that include a pentagonal defect. In our calculations, we employ a tight-binding model determining total and local density of states, transmission function and participation number. For the closed structure, we observe that the effect of the defect is concentrated mainly on energies near to zero, which is characteristic of edge states in graphene. The density of states and transmission functions for small energies show several peaks associated with the presence of quasi-bound states generated by the defect and localized edge states produced by both the circular boundaries of the finite lattice and induced by the presence of the pentagonal defect. These results have been checked by calculating the participation number, which is obtained from the eigenstates. We observe changes in the available quasi-bound states due to the defect and the creation of new peaks in the transmission function.

Graphene; Transport; Defects; 72.80.Vp; 72.15.Rn; 73.20.At


The theoretical and experimental study of properties of graphene has attracted the attention of many authors in the last few years since a method to isolate single graphene layers was developed (the authors Geim and Novoselov were awarded with the Nobel prize). These graphene sheets may be stable enough to be freely suspended [1], which allows us to use them in solid state experiments. Besides, the electronic properties of graphene are surprising: one finds new quasi-particles described by the Dirac equation at low energies that behave like massless particles. This opens the possibility to study quantum electrodynamics properties in solid-state devices and to carry out new developments, e. g., biosensors (see other studies [2-10]).

The influence of defects and edges in graphene properties has been widely studied [11-13]. Other authors made similar studies to ours but considered different geometries: Zhang et al. [14,15] worked on transport with narrow ballistic ribbon of graphene with zigzag edges including topological defects. Carpio et al. [12] studied the electronic properties in a similar geometry but with dislocations consisting of heptagon-pentagon pairs in an hexagon lattice. The aim of this article is to study the combined effects, on the electronic properties, of the geometric structure and the presence of topological defects. We then consider a model with a pentagonal defect (disclination), henceforth PD, at the centre of a graphene sheet with a circular shape (see Figure 1). We characterize the electronic and transport properties with the local and total density of states, participation number and transmission function. This work can be useful for the search of structures suitable for confinement of Dirac electrons, which are the basis for the construction of nanoelectronic devices with graphene.

thumbnailFigure 1. Graphene sheet with the topological defect. Schematic geometry of the graphene sheet studied in this work. Note the pentagonal defect placed at its centre (in red colour). This structure is connected to two semi-infinite graphene leads, which are partially shown in the figure (red colour).


Our geometry consists of a finite circular graphene quantum dot with 1,011 carbon atoms. For electronic transport, the quantum dot is connected to two semi-infinite leads. In Figure 1, we show the quantum dot and, partially, the semi-infinite leads. We employ a tight-binding model that only takes into account one π-orbital per atom. The overlap energy between nearest neighbours is taken as t=2.66 eV, where second-neighbour interactions are neglected. The advantage of using a single-band π-orbital model resides in its simplicity, being the general features of electronic transport in very good agreement with those obtained by more sophisticated approaches. The hamiltonian can then be written as

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


where <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M2">View MathML</a> are the creation/annihilation operators of an electron in site i. We expand the wave function in terms of the site base. <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M3">View MathML</a>, where <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M4">View MathML</a> is the amplitude probability that the electron is to be in site i for the eigenstate k. We need to solve <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M5">View MathML</a>. Four quantities are calculated to characterize the nature of the electronic and transport properties on two-circled structures, with PD and defect-free (ND) structures: the total density of states N(E), the local density of states ρ(i,E), the participation number P(E) and the transmission function T(E).

Electronic properties for the closed system

The density of states is determined from the energy spectrum as

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


Another useful property is the local density of states:

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


which measures how each site i contributes to the complete spectrum. For a fixed E, it characterizes the spatial nature of the state: it is localized when only few sites contribute to that energy, or extended when more sites participate. Finally, the participation number is defined as [16]

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


It assesses the wave function spreading so it can help to find out the localized or extended nature of an electronic state. For a completely localized wave function Ψk(i) is approximately δkiP≈1 while for a typical delocalized wave function on D atoms, Ψk(i) is approximately <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M9">View MathML</a>, and then PD.

Electronic transport properties: Green’s function method

We calculate transport properties for our graphene structures such as the density of states and the transmission function using Green’s function method. In order to obtain Green’s function, we use the following expression [17]:

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


where <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M11">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M12">View MathML</a> are the self-energy terms of left and right leads, respectively, and <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M13">View MathML</a> is the Hamiltonian of the conductor, i.e., in our case, the circular graphene sheet plus a few unit cells of the leads. In our approach, the contact leads at opposite sides of the circular graphene sheet is the graphene sheet itself extended to make the leads semi-infinite. This is equivalent to have reflectionless contacts in macroscopic conductors. Self-energy terms are calculated using the prescription <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M14">View MathML</a>, where <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M15">View MathML</a> is Green’s function of the semi-infinite lead (right or left) evaluated on sites k and l, which are in contact with sites i and j in the circular graphene sheet.

We only need to calculate <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M16">View MathML</a> in the sites in contact with the conductor. To do that, we use the formalism developed by López Sancho et al. [18]. This method has the advantage that the number of iterations close to singularities is very low compared to other transfer matrix methods, so it converges very fast and has been applied to graphene layers by other authors (see e.g. [19]). In this scheme, Green’s function is <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M17">View MathML</a>, where <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M18">View MathML</a> is the Hamiltonian of one isolated graphene cell in the lead, and <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M19">View MathML</a> is the matrix that takes into account the interaction between two consecutive cells. For the calculation of T, we use the iterative method described in [18].

From Green’s function of the graphene structure, we calculate the transmission function and the density of states as [17]

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


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


In Equation 6, GR/A are the retarded and advanced Green’s functions, respectively, and <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/258/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/258/mathml/M22">View MathML</a>. We denote the trace of the matrix considered by “Tr”, which is extended over the whole matrix.

Results and discussion

We have obtained different properties of graphene structures with and without pentagonal defects, in order to evaluate the influence of the defect and the geometry on their electronic properties. For the closed structure, we have calculated the total density of states, which is shown in Figure 2, for both the defect-free structure (dashed line) and with PD (continuous line). We see that the density for the structure with PD shows a shoulder near E=0, indicating the existence of additional edge states induced by the presence of the PD and the circular shape of the structure. The behaviour of the participation number confirmes these findings (see Figure 3a for the ND and Figure 3b for the PD structures). One can observe that PPD<PND, which suggests that the PD induces localized states beyond the known edge states in the structure with ND. In addition, near E=0, the structure of the participation number is slightly different in both structures. This shows that the nature and the number of the edge states are different, which can be seen in the local density function for this energy.

thumbnailFigure 2. Density of states for the closed structure. Total density of states for the structure with no defects (dashed line) and with a pentagonal defect (continuous line). Note the non-null density near zero as a manifestation of the edge defects.

thumbnailFigure 3. Participation number for the closed structure. Participation number P(E) of the available energy states for the structure with no defects (a) and with the pentagonal defect in the centre (b).

The edge states are localized so only few states contribute to a certain site; this is shown in Figure 4 for the local density of states at E=0 and ρ(i,E=0) for both the ND (Figure 4a) and PD (Figure 4b). Clearly, these are edge states, and the PD structure shows contribution from two zones, compared to the ND structure with one. The effect of the PD on the density of states near E=0 is of geometrical nature; the whole structure is affected by the presence of the pentagon since it changes the relative orientation of the edge sites and induces the creation of edge states. This has to do mainly with the atom rearrangement in the lower part of the structure, which creates new edge states and, clearly, the PD sites do not have an explicit contribution to such sites. For larger values of E, in the local density of ρ(i,E=2.6), more sites contribute to that energy (see Figure 5). Specifically, we see the contribution of sites around the PD as it can be seen in Figure 5b, where a star shape appears. The rest of the sites contribute more or less similarly to the structure with ND (Figure 5a).

thumbnailFigure 4. Local density of states for E = 0. Spatial distribution of the local density for ρ(i,E) for the energy E close to zero in (a) a structure with no defect and (b) one with the pentagonal defect in the centre. Due to single-bond atoms (see Figure 1), the quantum dot is not fully symmetric around a central vertical axis.

thumbnailFigure 5. Local density of states for E = 2.6 eV. Same as Figure 6 but for the energy E = 2.6 eV.

thumbnailFigure 6. Density of states for the open structure. Density of states for the graphene sheet with the pentagon at its centre (red line) and without it (black line). Note the displacement of the different peaks.

As the change in behaviour with the presence of PD is near zero energy (around the Fermi energy), we concentrate in the analysis of the transport properties around such energy. We have also checked our previous results in the open structure calculating the density of states (Figure 6) and the transmission function (Figure 7). The density of states shows several peaks associated with both the presence of quasi-bound states (due to the circular confinement and the defect) and localized edge states due to circular boundaries of the finite lattice. These results are clearly observed in the peak structure of the transmission function (Figure 7), where we observe changes in the quasi-bound states available to transport and the creation of new peaks in the transmission function.

thumbnailFigure 7. Transmission function for the open structure. Transmission function with the pentagon at its centre (red line) and without it (black line). Note the additional peaks for the case with the defect.


We have investigated the electronic and transport properties of circular graphene layers with a pentagonal disclination. In particular, using a tight-binding model, we have calculated the density of states, transmission function, participation number and local density of states of the structure with and without defects. The density of states for the structure with the PD shows several peaks that are associated with new localized states, which have been checked by calculating the local density of states and the participation number. We observe changes in the available quasi-bound states due to the defect and new peaks of the transmission function. Comparing these results, we conclude that there are more quasi-bound states in the structure with the defect, states associated with both the presence of quasi-bound states related to the atoms belonging to the defect and others due to the circular confinement and edge states due to circular boundaries of the finite lattice and the defect.


ND: Defect-free structure; PD: Pentagonal defect.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The work presented here was carried out in collaboration among all authors. FR defined the research theme. EJ carried out the calculations under APG’s supervision. All of them have discussed the results and wrote the manuscript. All authors read and approved the final manuscript.


FR would like to acknowledge the DGAPA project PAPPIT IN112012 for their financial support and sabbatical scholarship at the UPCT.


  1. Meyer JC, Geim AK, Katsnelson MI, Novoselov KS, Booth TJ, S R: The structure of suspended graphene sheets.

    Phys Rev Lett 1994, 72:1878. PubMed Abstract | Publisher Full Text OpenURL

  2. Castro Neto AH, Guinea F, Peres NMR, Novoselov KS, Geim AK: The electronic properties of graphene.

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

  3. Geim AK: Graphene: Status and prospects.

    Science 2009, 324:1530. PubMed Abstract | Publisher Full Text OpenURL

  4. Ihn T, Güttinger J, Molitor F, Schnez S, Schurtenberger E, Jacobsen A, Hellmüller S, Frey T, Dröscher S, Stampfer C, Ensslin K: Graphene single electron transistors.

    Mater Today 2010, 13:44. OpenURL

  5. Molitor F, Güttinger J, Stampfer C, Dröscher S, Jacobsen A, Ihn T, Ensslin K: Electronic properties of graphene nanostructures.

    J Phys: Condens Matter 2011, 23:243201. Publisher Full Text OpenURL

  6. Cooper DR, D’Anjou B, Ghattamaneni N, Harack B, Hilke M, Horth A, Majlis N, Massicotte M, Vandsburger L, Whiteway E, Yu V: Experimental review of Graphene.

    ISRN Condens Matter Phys 2012, 2012:501686. OpenURL

  7. Kim JH, Jung JM, Kwak JY, Jeong JH, Choi BC, Lim KT: Preparation of properties of SWNT/Graphene oxide type flexible transparent conductive film.

    J Nanosci Nanotechnol 2011, 11:7424. PubMed Abstract | Publisher Full Text OpenURL

  8. Yun JS, Yang KS, Kim DH: Multifunctional polydiacetylene-Graphene nanohybrids for biosensor application.

    J Nanosci Nanotechnol 2011, 11:5663. PubMed Abstract | Publisher Full Text OpenURL

  9. Zhang L, Xing Y, He N, Zhang Y, Lu Z, Zhang J, Zhang Z: Preparation of Graphene quantum dots for bioimaging application.

    J Nanosci Nanotechnol 2012, 12:2924. PubMed Abstract | Publisher Full Text OpenURL

  10. Islam MS, Kouzani AZ, Dai XJ, Michalski WP, Gholamhosseini H: Design and analysis of a multilayer localized surface plasmon resonance Graphene biosensor.

    J Nanosci Nanotechnol 2012, 8:380. OpenURL

  11. Meyer JC, Kisielowski C, Erni R, Rossell MD, Crommie MF, Zettl A: Direct imaging of lattice atoms and topological defects in Graphene membranes.

    Nano Lett 2008, 8:3582. PubMed Abstract | Publisher Full Text OpenURL

  12. Carpio A, Bonilla LL, de Juan F, Vozmediano MAH: Dislocations in graphene.

    New J Phys 2008, 10:053021. Publisher Full Text OpenURL

  13. Rycerz A: Electron transport and quantum-dot energy levels in Z-shaped graphene nanoconstriction with zigzag edges.

    Acta Phys Polon A 2010, 118:238. OpenURL

  14. Zhang Y, Hu JP, Bernevig BA, Wang XR, Xie XC, Liu WM: Quantum blockade and loop currents in graphene with topological defects.

    Phys Rev B 2008, 78:155413. OpenURL

  15. Zhang Y, Hu JP, Bernevig BA, Wang XR, Xie XC, Liu WM: Impurities in graphene.

    Phys Status Solidi A 2010, 207:2726. Publisher Full Text OpenURL

  16. Wegner FJ: Inverse participation ratio in 2+Epsilon dimensions.

    Z Phys B 1980, 36:209. Publisher Full Text OpenURL

  17. Datta S: Electronic Transport in Mesoscopic Systems. Cambridge: Cambridge University Press; 1995. OpenURL

  18. López Sancho MP, López Sancho JM, Rubio J: Quick iterative scheme for the calculation of transfer matrices: application to Mo (100).

    J Phys F: Met Phys 1984, 14:1205. Publisher Full Text OpenURL

  19. Li TC, Lu SP: Quantum conductance of graphene nanoribbons with edge defects.

    Phys Rev B 2008, 77:085408. OpenURL