Abstract
We investigate the electronic and transport properties of circular graphene structures (quantum dots) that include a pentagonal defect. In our calculations, we employ a tightbinding 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 quasibound 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 quasibound states due to the defect and the creation of new peaks in the transmission function.
Keywords:
Graphene; Transport; Defects; 72.80.Vp; 72.15.Rn; 73.20.AtBackground
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 quasiparticles described by the Dirac equation at low energies that behave like massless particles. This opens the possibility to study quantum electrodynamics properties in solidstate devices and to carry out new developments, e. g., biosensors (see other studies [210]).
The influence of defects and edges in graphene properties has been widely studied [1113]. 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 heptagonpentagon 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.
Figure 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 semiinfinite graphene leads, which are partially shown in the figure (red colour).
Methods
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 semiinfinite leads. In Figure 1, we show the quantum dot and, partially, the semiinfinite leads. We employ a tightbinding model that only takes into account one πorbital per atom. The overlap energy between nearest neighbours is taken as t=2.66 eV, where secondneighbour interactions are neglected. The advantage of using a singleband π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
where are the creation/annihilation operators of an electron in site i. We expand the wave function in terms of the site base. , where is the amplitude probability that the electron is to be in site i for the eigenstate k. We need to solve . Four quantities are calculated to characterize the nature of the electronic and transport properties on twocircled structures, with PD and defectfree (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
Another useful property is the local density of states:
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]
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 δ_{ki}→P≈1 while for a typical delocalized wave function on D atoms, Ψ^{k}(i) is approximately , and then P≈D.
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]:
where and are the selfenergy terms of left and right leads, respectively, and 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 semiinfinite. This is equivalent to have reflectionless contacts in macroscopic conductors. Selfenergy terms are calculated using the prescription , where is Green’s function of the semiinfinite 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 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 , where is the Hamiltonian of one isolated graphene cell in the lead, and 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]
In Equation 6, G^{R/A} are the retarded and advanced Green’s functions, respectively, and . 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 defectfree 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 P_{PD}<P_{ND}, 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.
Figure 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 nonnull density near zero as a manifestation of the edge defects.
Figure 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).
Figure 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 singlebond atoms (see Figure 1), the quantum dot is not fully symmetric around a central vertical axis.
Figure 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 quasibound 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 quasibound states available to transport and the creation of new peaks in the transmission function.
Figure 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.
Conclusion
We have investigated the electronic and transport properties of circular graphene layers with a pentagonal disclination. In particular, using a tightbinding 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 quasibound states due to the defect and new peaks of the transmission function. Comparing these results, we conclude that there are more quasibound states in the structure with the defect, states associated with both the presence of quasibound 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.
Abbreviations
ND: Defectfree 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.
Acknowledgements
FR would like to acknowledge the DGAPA project PAPPIT IN112012 for their financial support and sabbatical scholarship at the UPCT.
References

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

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

Geim AK: Graphene: Status and prospects.
Science 2009, 324:1530. PubMed Abstract  Publisher Full Text

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.

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

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.

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

Yun JS, Yang KS, Kim DH: Multifunctional polydiacetyleneGraphene nanohybrids for biosensor application.
J Nanosci Nanotechnol 2011, 11:5663. PubMed Abstract  Publisher Full Text

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

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

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

Carpio A, Bonilla LL, de Juan F, Vozmediano MAH: Dislocations in graphene.
New J Phys 2008, 10:053021. Publisher Full Text

Rycerz A: Electron transport and quantumdot energy levels in Zshaped graphene nanoconstriction with zigzag edges.

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

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

Wegner FJ: Inverse participation ratio in 2+Epsilon dimensions.
Z Phys B 1980, 36:209. Publisher Full Text

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

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

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