This article is part of the series International Conference on Superlattices, Nanostructures and Nanodevices (ICSNN 2010).

Open Access Nano Express

Defect symmetry influence on electronic transport of zigzag nanoribbons

Hui Zeng12*, Jean-Pierre Leburton234, Yang Xu5 and Jianwei Wei6

Author Affiliations

1 College of Physical Science and Technology, Yangtze University, Jingzhou, Hubei 434023, China

2 Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA

3 Department of Electrical and Computer Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA

4 Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA

5 Department of Information Science and Electronic Engineering, Zhejiang University, Hangzhou, Zhejiang 310027, China

6 College of Mathematics and Physics, Chongqing University of Technology, Chongqing 400054, China

For all author emails, please log on.

Nanoscale Research Letters 2011, 6:254  doi:10.1186/1556-276X-6-254


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


Received:19 August 2010
Accepted:24 March 2011
Published:24 March 2011

© 2011 Zeng 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.

Abstract

The electronic transport of zigzag-edged graphene nanoribbon (ZGNR) with local Stone-Wales (SW) defects is systematically investigated by first principles calculations. While both symmetric and asymmetric SW defects give rise to complete electron backscattering region, the well-defined parity of the wave functions in symmetric SW defects configuration is preserved. Its signs are changed for the highest-occupied electronic states, leading to the absence of the first conducting plateau. The wave function of asymmetric SW configuration is very similar to that of the pristine GNR, except for the defective regions. Unexpectedly, calculations predict that the asymmetric SW defects are more favorable to electronic transport than the symmetric defects configuration. These distinct transport behaviors are caused by the different couplings between the conducting subbands influenced by wave function alterations around the charge neutrality point.

Introduction

As a truly two-dimensional nanostructure, graphene has attracted considerable interest, mainly because of its peculiar electronic and transport properties described by a massless Dirac equation [1,2]. As such, it is regarded as one of the most promising materials since its discovery [3-5] because charge carriers exhibit giant intrinsic mobility and long mean-free path at room temperature [6,7], suggesting broad range of applications in nanoelectronics [8-11]. Several experimental [4,8,12,13] and theoretical [2,14,15] studies are presently devoted to the electronic, transport, and optical properties [16] of graphene. By opening an energy gap between valence and conduction bands, narrow graphene nanoribbons (GNR) are predicted to have a major impact on transport properties [17,18]. Most importantly, GNR-based nano-devices are expected to behave as molecular devices with electronic properties similar to those of carbon nanotubes (CNTs) [19,20], as for instance, Biel et al. [21] reported a route to overcome current limitations of graphene-based devices through the fabrication of chemically doped GNR with boron impurities.

The investigation of transport properties of GNRs by various experimental methods such as vacancies generation [22], topological defects [23], adsorption [24], doping [25], chemical functionalization [26-28], and molecular junctions [29] have been reported. Meanwhile, defective GNR with chemically reconstructed edge profiles also have been experimentally evidenced [30] and have recently received much attention [31,32]. In particular, Stone-Wales (SW) defects, as one type of topological defects, are created by 90° rotation of any C-C bond in the hexagonal network [33], as shown by Hashimoto et al. [34]. More recently, Meyer et al. [35] have investigated the formation and annealing of SW defects in graphene membranes and found that the existence of SW defects is energetically more favorable than in CNTs or fullerenes. Therefore, the influences of SW defects on electronic transport of GNRs is crucial for the understanding of the physical properties of this novel material and for its potential applications in nanoelectronics.

In this brief communication, we investigate the influence of SW defects on the electronic transport of zigzag-edged graphene nanoribbons (ZGNRs). It is found that the electronic structures and transport properties of ZGNRs with SW defects can very distinctively depend on the symmetry of SW defects. The transformation energies obtained for symmetric SW defects and asymmetric SW defects are 5.95 and 3.34eV, respectively, and both kinds of defects give rise to quasi-bound impurity states. Our transport calculations predict different conductance behavior between symmetric and asymmetric SW defects; asymmetric SW defects are more favorable for electronic transport, while the conductance is substantially decreased in the symmetric defects configuration. These distinct transport behaviors result from the different coupling between the conducting subbands influenced by the wave function symmetry around the charge neutrality point (CNP).

Model and methods

The optimization calculations are done by using the density functional theory utilized in the framework of SIESTA code [36,37]. We adopt the standard norm-conserving Toullier-Martins [38] pseudopotentials orbital to calculate the ion-electron interaction. The numerical double-ζ polarized is used for basis set and the plane cutoff energy is chosen as 200 Ry. The generalized gradient approximation [39] proposed by Perdew and Burke and Ernzerhof was employed to calculate exchange correction term. All nanostructure geometries were converged until no forces acting on all atoms exceeded 0.01eV/Å.

The electronic transport properties of the nanoribbon device have been performed by using non-equilibrium Green's function (NEGF) methodology [40,41]. In order to self-consistently calculate the electrical properties of nanodevices, we construct the two-probe device geometry where the central region contains the SW defects and both leads consist each of the two supercell pristine ZGNR, as shown in Figure 1. The equilibrium conductance G is obtained from the Landauer formula such that G = G0T (E), where G0 is the quantum conductance with relationship . The transmission coefficient T as a function of the electron energy E is given by

(1)

where Σl (Σr) represents the self-energies of the left (right) electrode, GR (GA) is retard (advanced) Green's function. It is calculated from the relation:

(2)

where HS is the Hamiltonian of the system. More details about the NEGF formalism can be found in Ref. [42].

thumbnailFigure 1. (Color online) Schematics of the molecular device considered in the calculation of ZGNR. The whole device is composed of scattering region and two electrodes containing the corresponding pristine ZGNRs. The SW defects are highlighted by yellow atoms (a) ZGNR with symmetric SW defects, some C-C bond lengths and angles are shown by (c); (b) ZGNR with asymmetric SW defects, while some C-C bond lengths and angles are shown by (d).

In this study, we consider symmetric and asymmetric SW defects contained in 6-ZGNRs, where 6 denotes the number of zigzag chains (dimers) across the ribbon width [18]. Taking into account screening effects between electrodes and central molecules, we use 10-unit cell's length as scattering regions, and 2 units as electrodes to perform transport calculation. The electron temperature in the calculation is set to be 300 K.

Results and discussions

In Figure 1, we show the geometry of defective ZGNR after relaxation. After introducing symmetric SW defects, the GNR shrinks along the width axis, by 0.526 Å, and correspondingly, the nearest four H atoms move toward the central region by 0.21 Å. As a result, the bond angles of the edge near the SW defects are reduced from 120 to 116°, as shown in Figure 1c, d. In contrast to the shrinking along the width axis, the SW defects stretch from 4.88 to 5.38 Å along the length axis direction. No distinct change for the H-C bond length at the edge is observed. Thus, the effect of symmetric SW defects on the geometry modification is limited to the defective area, with mirror reflection around their axis. However, the presence of asymmetric SW defects to the geometry modification is far more complex. They twist the whole structure by shifting the left side upward, while the right side is downward shifted. Hence, the mirror symmetry is broken because of the asymmetric SW. The transformation energies for symmetric and asymmetric SW defects are 5.95 and 3.34eV, respectively. These results imply that the asymmetric SW defects are energetically more favorable than the symmetric SW defects.

Wave functions of electronic states at the Gamma point of the highest-occupied electronic states (HOES) and the lowest-unoccupied electronic states (LUES) are depicted in Figure 2. As expected, the wave functions of the pristine even-index ZGNR at the Gamma-point associated to the HOES and LUES exhibit the well-defined parity with respect to the mirror plane, and their eigenstates in the case of symmetric SW defects, the HOES and LUES, keep the same parity because of the potential induced by the symmetric defects [25,43]. Note that, although the wave functions of both the pristine and symmetric SW defects have well-defined parity, the sign of their wave functions, especially for the electronic states below the CNP, are precisely opposite. For the asymmetric SW defects, the well-defined parity of the wave functions is not preserved. Moreover, the wave function symmetry in this configuration is broken leading to substantial electron backscattering below and above the CNP.

thumbnailFigure 2. (Color online) Wave functions at the Gamma point of defective ZGNRs. Wave functions at the Gamma point associated with the LUES above the CNP (top Panel) and the HOES below the CNP (bottom Panel) for ZGNR with no defects (a, d), symmetric SW defects (b, e) asymmetric SW defects (c, f). Dark gray (blue online) and light gray (red online) colors correspond to the opposite signs of the wave function.

The central issue of this study is to investigate the influence of SW defects in the ZGNRs on their electronic and transport behavior. ZGNRs are known to present very peculiar electronic structure, that is, strong edge effects at low energies originated from the wave functions localized along the GNR edges [44]. Spin-unpolarized calculations reveal that all ZGNRs are metallic with the presence of sharply localized edge states at the CNP [25,43,44], while ab initio calculation with spin effect taken into consideration found that a small band gap opens up [18]. The electronic band structures of defective nanoribbons and the corresponding pristine GNRs are shown for comparison. In the case of pristine GNR, zone-folded effects give rise to nondegenerated bands for α- and β-spin states, and the corresponding spin bands shift upward and downward with respect to the CNP, respectively. It also leads to gapless electronic structure as well as 3G0 conductance in the vicinity of CNP (see Figure 3). Meanwhile, zone-folded effects create more subbands near the CNP, namely, four α-spin subbands around 0.4eV and four β-spin subbands around 0.4eV. The presence of symmetric SW defects substantially split the electronic bands, especially for the β-spin bands above the CNP, resulting from the bands anticrossing at Γ or π point. More importantly, the symmetric defects open a band gap of about 0.12eV for α-spin bands and 0.09eV for β-spin bands, which is attributed to the mismatch coupling between its LUES and HOES wave functions due to the presence of defects. It is interesting to note that a defect state deriving from the α-spin subband is located at about 1.15eV above the CNP producing a localized state, where complete backscattering is obtained (see red dashed line in Figure 3). Thus, these changes in the band structures arising from introducing symmetric SW defects are unfavorable to electronic transport. In contrast to the extensive split produced by the symmetric SW defects, the electronic structure modification due to the asymmetric SW defects is slight. Except for some bands splitting that could be unfavorable to electron transport, the band structure away from the CNP does not experience much change. Similar to the emergence of defect states induced by the symmetric SW configuration, two defect states are observed in the asymmetric SW configurations; one defect state arising from the α-spin subband locates at about 0.62eV above the CNP, and the other one from the β-spin is -1.20eV below the CNP. Both defects give rise to localized states that lead to conductance gaps (see, dotted line in Figure 3). Overall, the band structure results reveal that the SW defect states near the CNP lead to complete electron backscattering region, where the location depends on the spatial symmetry of the defects.

thumbnailFigure 3. (Color online) Electronic band structures of defective ZGNRs. (a) for the pristine, (b) for the symmetric SW defects and (c) for the asymmetric SW defects. The solid red (dotted blue) line denotes the α-spin (β-spin) bands. The dashed black line indicates the CNP, and solid circles indicate defect states.

The electronic transport results are displayed in Figure 3. The states induced by H atoms at the edge produce a conductance peak in the vicinity of CNP in the pristine ZGNR. In this study, our results show a good agreement with previous studies [43-45]. The first conductance plateau corresponding to the occupied and unoccupied states is G0. In the case of symmetric SW defects in the ZGNR, the conductance in the vicinity of CNP is decreased as a result of the four H atoms shrinking. The conductance with symmetric defects remarkably decreases below the CNP, manifesting monotonous reduction of conductance with increasing electron energy. We attribute this effect to the antisymmetry (opposite sign at every position) of wave functions, with respect to the pristine GNR in the wave functions (see Figure 4e) that block the electronic transport. On the other hand, the orientation of about 50% of all wave functions corresponding to LUES is reversed, which gives rise to a conducting plateau (about 0.5G0) that ranged from 0.04 to 0.8eV above the CNP. More importantly, strong electron backscattering induced by the coupling between all states are expanded to lead to full suppression of the conduction channel at particular resonance energies. Accordingly, a smooth conductance valley around 1.12eV corresponding to complete electron backscattering is observed. Concerning the transport properties of asymmetric SW configuration, we find that the absence of conductance peak at the CNP is due to the breaking of edge states. In addition, localized states in the vicinity of CNP lead to reduced conductance. The main feature of the first conducting plateau below the CNP is preserved except for the smooth conductance valley located at about -1.2eV. This illustrates the obvious different transport behaviors between the symmetric and asymmetric SW defects. We indeed found that such different transport behaviors result from different coupled electronic states supported by the wave function results. The HOES and LUES wave functions of asymmetric SW defect configuration are very similar to that of the pristine GNR except for the defective area. Therefore, the first conducting plateau near the CNP is preserved for the asymmetric configuration. Naturally, the asymmetric SW defects are responsible for the two conductance valleys, namely, a smooth valley at -1.2eV and a sharp valley at 1.48eV. The large reduction of conductance at these areas induced by the asymmetric SW defects corresponds to complete electron backscattering region. which is different from the situation in CNTs, where SW defects induce suppression of only half of the conductance channels [46]. However, the impact of the two conductance valleys on the ZGNRs is limited because they are far away from the CNP. The transport properties of asymmetric SW configuration are predicted to be comparable with that of the pristine GNR in spite of non-preservation of the geometry and wave function symmetry for the former. We note that similar results have been obtained under spin-dependent calculation by Ren et al. [47] very recently. Overall, the electronic transport calculations predict that it is more likely to be observed for asymmetric SW defects in the ZGNR, since these defects are more favorable for electronic transport in contrast to the substantially transport degradation in the symmetric defects configuration.

thumbnailFigure 4. (Color online) Conductance of defective ZGNRs as a function of electron energy. The black solid line, red dashed line, and blue dotted line denote the results of pristine, symmetric SW, and asymmetric SW defects.

Conclusion

In summary, we investigate the influence of local structural defects on the electronic transport of ZGNR using first principles calculations. The transformation energies reveal that the asymmetric SW defects is energetically more favorable than the symmetric SW defects. Both defects give rise to complete electron backscattering region that depends on the spatial symmetry of the defects. Our transport calculations predict that the asymmetric SW defects are more favorable for electronic transport in contrast to the substantially decreased in the symmetric defects configuration. We attribute these distinct transport behaviors to the different coupling between the conducting subbands influenced by the wave function modification around the CNP.

Abbreviations

CNP: charge neutrality point; GNR: graphene nanoribbons; HOES: highest-occupied electronic states; LUES: lowest-unoccupied electronic states; SW: Stone-Wales; ZGNR: zigzag-edged graphene nanoribbon.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

HZ carried out molecular dynamic studies, participated in the sequence alignment and drafted the manuscript. JL participated in the design of the study and the sequence alignment. YX participated in the sequence alignment. JW took part in the simulation and participated in the sequence alignment. All the authors discussed the results.

Acknowledgements

The authors gratefully thank Prof. K.-L. Yao and Dr. M. A. Kuroda for their technical assistance with performing ab initio transport properties and the relax calculation in the Mac OS X Turing cluster. This study is financially supported by the Scientific Research Foundation of Yangtze University (Grant No.801080010111) and the Chongqing University of Technology (Grant No.2008EDJ01), and the Natural Science Foundation of China under Grant No.11047176.

References

  1. Geim AK, Novoselov KS: The rise of graphene.

    Nat Mater 2007, 6:183-191. 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-162. Publisher Full Text OpenURL

  3. Novoselov KS, Geim AK, Morozov SV, Jiang D, Zhang Y, Dubonos SV, Grig-orieva IV, Firsov AA: Electric field effect in atomically thin carbon films.

    Science 2004, 306:666-669. PubMed Abstract | Publisher Full Text OpenURL

  4. Novoselov KS, Geim AK, Morozov SV, Jiang D, Katsnelson MI, Grigorieva IV, Dubonos SV, Firsov AA: Two-dimensional gas of massless Dirac fermions in graphene.

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

  5. Zhang YB, Tan Y-W, Stormer HL, Kim P: Experimental observation of the quantum Hall effect and Berry's phase in graphene.

    Nature 2005, 438:201-204. PubMed Abstract | Publisher Full Text OpenURL

  6. Berger C, Song Z, Li X. Wu X, Brown N, Naud C, Mayou D, Li T, Hass J, Marchenkov AN, Conrad EH, First PN, de Heer WA: Electronic confinement and coherence in patterned epitaxial graphene.

    Science 2006, 312:1191-1196. PubMed Abstract | Publisher Full Text OpenURL

  7. Orlita M, Faugeras C, Plochocka P, Neugebauer P, Martinez G, Maude DK, Barra A-L, Sprinkle M, Berger C, de Heer WA, Potemskil M: Approaching the dirac point in high-mobility multilayer epitaxial graphene.

    Phys Rev Lett 2008, 101:267601. PubMed Abstract | Publisher Full Text OpenURL

  8. Bunch JS, van der Zande AM, Verbridge SS, Frank IW, Tanenbaum DM, Parpia JM, Craighead HG, McEuen PL: Electromechanical resonators from graphene sheets.

    Science 2007, 315:490-493. PubMed Abstract | Publisher Full Text OpenURL

  9. Yan QM, Huang B, Yu J, Zheng FW, Zang J, Wu J, Gu BL, Liu F, Duan WH: Intrinsic current-voltage characteristics of graphene nanoribbon transistors and effect of edge doping.

    Nano Lett 2007, 7:1469-1473. PubMed Abstract | Publisher Full Text OpenURL

  10. Martins TB, Miwa RH, da Silva AJR, Fazzio A: Electronic and transport properties of Boron-doped graphene nanoribbons.

    Phys Rev Lett 2007, 98:196803. PubMed Abstract | Publisher Full Text OpenURL

  11. Fiori G, Iannaccone G: Simulation of graphene nanoribbon field-effect transistors.

    IEEE Electron Dev Lett 2007, 28:760-762. Publisher Full Text OpenURL

  12. Novoselov KS, Jiang D, Schedin F, Booth TJ, Khotkevich VV, Morozov SV, Geim AK: Two-dimensional atomic crystals.

    Proc Natl Acad Sci USA 2005, 102:10451-10453. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Lee C, Wei X, Kysar JW, Hone J: Measurement of the elastic properties and intrinsic strength of monolayer graphene.

    Science 2008, 321:385-388. PubMed Abstract | Publisher Full Text OpenURL

  14. Fujita M, Wakabayashi K, Nakada K, Kusakabe K: Peculiar localized state at zigzag graphite edge.

    J Phys Soc Jpn 1996, 65:1920-1923. Publisher Full Text OpenURL

  15. Miyamoto Y, Nakada K, Fujita M: First-principles study of edge states of H-terminated graphitic ribbons.

    Phys Rev B 1999, 59:9858-9861. Publisher Full Text OpenURL

  16. Dresselhaus MS, Jorio A, Hofmann M, Dresselhaus G, Saito R: Perspectives on carbon nanotubes and graphene Raman spectroscopy.

    Nano Lett 2010, 10:751-758. PubMed Abstract | Publisher Full Text OpenURL

  17. Nakada K, Fujita M, Dresselhaus G, Dresselhaus MS: Edge state in graphene ribbons: Nanometer size effect and edge shape dependence.

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

  18. Son Y-W, Cohen ML, Louie SG: Energy gaps in graphene nanoribbons.

    Phys Rev Lett 2006, 97:216803. PubMed Abstract | Publisher Full Text OpenURL

  19. Avouris Ph: Molecular electronics with carbon nanotubes.

    Acc Chem Res 2002, 35:1026-1034. PubMed Abstract | Publisher Full Text OpenURL

  20. Weitz RT, Zschieschang U, Eenberger F, Klauk H, Burghard M, Kern K: High-performance carbon nanotube field effect transistors with a self-assembled monolayer gate dielectric.

    Nano Lett 2007, 7:22-27. PubMed Abstract | Publisher Full Text OpenURL

  21. Biel B, Triozon F, Blase X, Roche S: Chemically-induced mobility gaps in graphene nanoribbons.

    Nano Lett 2009, 9:2725-2729. PubMed Abstract | Publisher Full Text OpenURL

  22. Carlsson JM, Scheffer M: Structural, electronic, and chemical properties of nanoporous carbon.

    Phys Rev Lett 2006, 96:046806. PubMed Abstract | Publisher Full Text OpenURL

  23. Cortijo A, Vozmediano MAH: Effects of topological defects and local curvature on the electronic properties of planar graphene.

    Nucl Phys B 2007, 763:293-308. Publisher Full Text OpenURL

  24. Choi SM, Jhi SH: Self-assembled metal atom chains on graphene nanoribbons.

    Phys Rev Lett 2008, 101:266105. PubMed Abstract | Publisher Full Text OpenURL

  25. Biel B, Blase X, Triozon F, Roche S: Anomalous doping effects on charge transport in graphene nanoribbons.

    Phys Rev Lett 2009, 102:096803. PubMed Abstract | Publisher Full Text OpenURL

  26. Cantele G, Lee Y-S, Ninno D, Marzari N: Spin channels in functionalized graphene nanoribbons.

    Nano Lett 2009, 9:3425-3429. PubMed Abstract | Publisher Full Text OpenURL

  27. López-Bezanilla A, Triozon F, Roche S: Chemical functionalization effects on armchair graphene nanoribbon transport.

    Nano Lett 2009, 9:2537-2541. PubMed Abstract | Publisher Full Text OpenURL

  28. Cervantes-Sodi F, Csányi G, Piscanec S, Ferrari AC: Edge-functionalized and substitutionally doped graphene nanoribbons: Electronic and spin properties.

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

  29. Wang B, Wang J, Guo H: Ab initio calculation of transverse spin current in graphene nanostructures.

    Phys Rev B 2009, 79:165417. Publisher Full Text OpenURL

  30. Koskinen P, Malola S, Häkkinen H: Evidence for graphene edges beyond zigzag and armchair.

    Phys Rev B 2009, 80:073401. Publisher Full Text OpenURL

  31. Wassmann T, Seitsonen AP, Saitta AM, Lazzeri M, Mauri F: Structure, stability, edge states, and aromaticity of graphene ribbons.

    Phys Rev Lett 2008, 1011:096402. Publisher Full Text OpenURL

  32. Dubois SM-M, Lopez-Bezanilla A, Cresti A, Triozon F, Biel B, Charlier J-C, Roche S: Quantum transport in graphene nanoribbons: Effects of edge reconstruction and chemical reactivity.

    ACS Nano 2010, 4:1971-1976. PubMed Abstract | Publisher Full Text OpenURL

  33. Stone AJ, Wales DJ: Theoretical studies of icosahedral C60 and some related structures.

    Chem Phys Lett 1986, 128:501-503. Publisher Full Text OpenURL

  34. Hashimoto A, Suenaga K, Gloter A, Urita K, Iijima S: Direct evidence for atomic defects in graphene layers.

    Nature 2004, 430:870-873. PubMed Abstract | Publisher Full Text OpenURL

  35. 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-3586. PubMed Abstract | Publisher Full Text OpenURL

  36. Ordejón P, Artacho E, Soler JM: Self-consistent order-N density-functional calculations for very large systems.

    Phys Rev B 1996, 53:R10441-R10444. OpenURL

  37. Soler JM, Artcho E, Gale JD, Garía A, Junquera J, Ordejón P, Sánchez-Portal D: The SIESTA method for ab initio order-N materials simulation.

    J Phys Condensed Matter 2002, 14:2745-2779. Publisher Full Text OpenURL

  38. Toullier N, Martins JL: Efficient pseudopotentials for plane-wave calculations.

    Phys Rev B 1993, 43:1993-2006. Publisher Full Text OpenURL

  39. Perdew JP, Burke K, Ernzerhof M: Generalized gradient approximation made simple.

    Phys Rev Lett 1996, 77:3865-3868. PubMed Abstract | Publisher Full Text OpenURL

  40. Taylor J, Guo H, Wang J: Ab initio modeling of quantum transport properties of molecular electronic devices.

    Phys Rev B 2001, 63:245407. Publisher Full Text OpenURL

  41. Brandbyge M, Mozos JL, Ordejón P, Taylor J, Stokbro K: Density-functional method for nonequilibrium electron transport.

    Phys Rev B 2002, 65:165401. Publisher Full Text OpenURL

  42. Datta S: Transmission function, S-matrix and Green's functions. In Electronic Transport in Mesoscopic Systems. New York: Cambridge University Press; 1995:117-163. OpenURL

  43. Li Z, Qian H, Wu J, Gu B-L, Duan W: Role of symmetry in the transport properties of graphene nanoribbons under bias.

    Phys Rev Lett 2008, 100:206802. PubMed Abstract | Publisher Full Text OpenURL

  44. Areshkin DA, Gunlycke D, White CT: Ballistic transport in graphene nanostrips in the presence of disorder: Importance of edge effects.

    Nano Lett 2007, 7:204-210. PubMed Abstract | Publisher Full Text OpenURL

  45. Oeiras P, Araújo-Moreira FM, da Silva FM: Defect-mediated half-metal behavior in zigzag graphene nanoribbons.

    Phys Rev B 2009, 80:073405. Publisher Full Text OpenURL

  46. Wei JW, Hu H-F, Zeng H, Wang Z-Y, Wang L, Peng P: Effects of nitrogen in Stone-Wales defect on the electronic transport of carbon nanotube.

    Appl Phys Lett 2007, 91:092121. Publisher Full Text OpenURL

  47. Ren Y, Chen K-Q: Effects of symmetry and Stone-Wales defect on spin-dependent electronic transport in zigzag graphene nanoribbons.

    J Appl Phys 2010, 107:044514. Publisher Full Text OpenURL