Abstract
In this work, we present a theoretical study of the transport properties of two finite and parallel armchair graphene nanoribbons connected to two semiinfinite leads of the same material. Using a single Πband tight binding Hamiltonian and based on Green’s function formalisms within a real space renormalization techniques, we have calculated the density of states and the conductance of these systems considering the effects of the geometric confinement and the presence of a uniform magnetic field applied perpendicularly to the heterostructure. Our results exhibit a resonant tunneling behaviour and periodic modulations of the transport properties as a function of the geometry of the considered conductors and as a function of the magnetic flux that crosses the heterostructure. We have observed AharonovBohm type of interference representing by periodic metalsemiconductor transitions in the DOS and conductance curves of the nanostructures.
Keywords:
Graphene nanostructures; Transport properties; Magnetic field effectsBackground
Graphene is a single layer of carbon atoms ordered in a twodimensional hexagonal lattice. In the literature, it is possible to find different experimental techniques in order to obtain graphene such as mechanical peeling, epitaxial growth or assembled by atomic manipulation of carbon monoxide molecules over a conventional twodimensional electron system at a copper surface [14]. The physical properties of this crystal have been studied over the last 70 years; however, the recent experimental breakthroughs have revealed that there are still a lot of open questions, such as timedependent transport properties of graphenebased heterostructures, the thermoelectric and thermal transport properties of graphenebased systems in the presence of external perturbations, the thermal transport properties of graphene under timedependent gradients of temperatures, etc.
On the other hand, graphene nanoribbons (GNRs) are quasi onedimensional systems based on graphene which can be obtained by different experimental techniques [58]. The electronic behaviour of these nanostructures is determined by their geometric confinement which allows the observation of quantum effects. The controlled manipulation of these effects, by applying external perturbations to the nanostructures or by modifying the geometrical confinement [913], could be used to develop new technological applications, such as graphenebased composite materials [14], molecular sensor devices [1517] and nanotransistors [18].
One important aspect of the transport properties of these quasi onedimensional systems is the resonant tunneling behaviour which, for certain configurations of conductors or external perturbations, appears into the system. It is has been reported that in S and Ushaped ribbons, and due to quasibound states present in the heterostructure, it is possible to obtain a rich structure of resonant tunneling peaks by tuning through the modification of the geometrical confinement of the heterostructure [19]. Another way to obtain resonant tunneling in graphene is considering a nanoring structure in the presence of external magnetic field. It has been reported that these annular structures present resonance in the conductance at defined energies, which can be tuned by gate potentials, the intensity of the magnetic field or by modifying their geometry [20]. From the experimental side, the literature shows the possibility of modulating the transport response as a function of the intensity of the external magnetic field. In some configuration of gate potential applied to the rings, it has been observed that the AharonovBohm oscillations have good resolution [2123].
In this context, in this work, we present a theoretical study of the transport properties of GNRbased conductors composed of two finite and parallel armchair nanoribbons (AGNRs) of widths N_{d} and N_{u}, and length L (measured in unit cell units), connected to two semiinfinite contacts of width N made of the same material. We have thought this system as two parallel ‘wires’ connected to the same reservoirs, whether the the leads are made of graphene or another material. This consideration allows us to study the transport of a hypothetical circuit made of graphene ‘wires’ in different scenarios. A schematic view of a considered system is shown in Figure 1. We have focused our analysis on the electronic transport modulations due to the geometric confinement and the presence of an external magnetic field. In this sense, we have studied the transport response due to variations of the length and widths of the central ribbons, considering symmetric and asymmetric configurations. We have obtained interference effects at low energies due to the extra spatial confinement, which is manifested by the apparition of resonant states at this energy range, and consequently, a resonant tunneling behaviour in the conductance curves. On the other hand, we have considered the interaction of electrons with a uniform external magnetic field applied perpendicular to the heterostructure. We have observed periodic modulations of the transport properties as function of the external field, obtaining metalsemiconductor transitions as function of the magnetic flux.
Figure 1. Schematic view of the conductor. Two finite armchair graphene ribbons (red lines). The length L of the conductor is measured in unitary cell units.
Methods
All considered systems have been described using a single Πband tight binding Hamiltonian, taking into account only the nearest neighbour interactions with a hopping γ_{0} = 2.75eV[24]. We have described the heterostructures using surface Green’s function formalism within a renormalization scheme [16,17,25]. In the linear response approach, the conductance is calculated using the Landauer formula. In terms of the conductor Green’s function, it can be written as [26]:
where , is the transmission function of an electron crossing the conductor region, is the coupling between the conductor and the respective lead, given in terms of the selfenergy of each lead: Σ_{L/R} = V_{C,L/R}g_{L/R }V_{L/R,C}. Here, V_{C,L/R} are the coupling matrix elements and g_{L/R }is the surface Green’s function of the corresponding lead [16]. The retarded (advanced) conductor Green’s function is determined by [26]: , where H_{C} is the hamiltonian of the conductor. Finally, the magnetic field is included by the Peierls phase approximation [2731]. In this scheme, the magnetic field changes the unperturbed hopping integral to , where the phase factor is determined by a line integral of the vector potential A by:
Using the vectors exhibited in Figure 1, R_{1} = (1, 0)a, and , where a = R_{n,m} = 1.42 Å, the phase factors for the armchair configuration in the Landau gauge A = (0, Bx) are given by:
where y_{n} is the carbon atom position in the transverse direction of the ribbons. In what follows, the Fermi energy is taken as the zero energy level, and all energies are written in units of γ_{0}.
Results and discussion
Unperturbed systems
Let us begin the analysis by considering the effects of the geometrical confinement. In Figure 2, we present results of (a) Local density of sates (LDOS) and (b) conductance for a conductor composed of two AGNRs of widths N_{d} = N_{u} = 5 connected to two leads of width N = 17 for different conductor lengths (L = 5, 10 and 20 unit cells). The most evident result is reflected in the LDOS curves at energies near the Fermi level. There are several sharp states at defined energies, which increase in number and intensity as the conductor length L is increased. These states that appear in the energy range corresponding to the gap of a pristine N = 5 AGNRs [24,32] correspond to a constructive interference of the electron wavefunctions inside the heterostructure, which can travel forth and back generating stationary (welllike) states. In this sense, the finite length of the central ribbons imposes an extra spatial confinement to electrons, as analogy of what happens in open quantum dot systems [16,17,19,33,34]. Independently of their sharp line shape, these discrete levels behave as resonances in the system allowing the conduction of electrons at these energies, as it is shown in the corresponding conductance curves of Figure 2b. It is clear that as the conductor length is increased, the number of conductance peaks around the Fermi level is also increased, tending to form a plateau of one quantum of conductance (G_{0} = 2e^{2}/h) at this energy range. These conductance peaks could be modulated by the external perturbations, as we will show further in this work.
Figure 2. LDOS and conductance for different geometries. (a) LDOS (black line) and (b) conductance of two AGRNs (red line) of widths N_{d} = N_{u} = 5, connected to two leads of widths N = 17 for different conductor lengths: L = 5, 10, 20 u.c. (c) Conductance of a system composed of two parallel N_{d} = 5 and N_{u} = 7 AGNRs of lengths L = 15. As a comparison, we have included the pristine cases (black and blue curves, respectively).
At higher energies, the conductance plateaus appear each as 2G_{0}, which is explained by the definition of the transmission probability T(E) of an electron passing through the conductor. In these types of heterostructures, if the conductor is symmetric (N_{u} = N_{d}), the number of allowed transverse channels are duplicated; therefore, electrons can be conduced with the same probability through both finite ribbons. On the other hand, in Figure 2c, we present results of conductance for a conductor of length L = 15 and composed of two AGNRs of widths N_{d} = 5 and N_{u} = 7, connected to two leads of widths N = 17. As a comparison, we have included the corresponding pristine cases. As it is expected, the conductance for an asymmetric configuration (red curve) reflects the exact addition of the transverse channels of the constituent ribbons, with the consequent enhancement of the conductance of the systems. Nevertheless, there is still only one quantum of conductance near the Fermi energy due to the resonant states of the finite system, whether the constituent ribbons are semiconductor or semimetal. We have obtained these behaviours for different configurations of conductor, considering variations in length and widths of the finite ribbons and leads.
Magnetic field effects
In what follows, we will include the interaction of a uniform external magnetic field applied perpendicularly to the conductor region. We have considered in our calculations that the magnetic field could affect the ends of the leads, forming an effective ring of conductor. The results of LDOS and conductance as a function of the Fermi energy and the normalized magnetic flux (ϕ/ϕ_{0}) for three different conductor configurations are displayed in the contour plots of Figure 3. The left panels correspond to a symmetric system composed of two metallic AGNRs of widths N_{u} = N_{d} = 5. The central panels correspond to an asymmetric conductor composed of two AGNRs of widths N_{d} = 5 (metallic) and N_{u} = 7 (semiconductor). The right panels correspond to a symmetric system composed of two semiconductor AGNRs of widths N_{u} = N_{d} = 7. All configurations have been considered of the same length L = 10 and connected to the same leads of widths N = 17. Finally, we have included as a reference, the plots of LDOS versus Fermi energy for the three configurations.
Figure 3. Magnetic field effects on LDOS and conductance. Contour plots of LDOS (lower panels) and conductance (upper panels) as a function of the Fermi energy and the magnetic flux crossing the hexagonal lattice for three different configurations of conductor. As a comparison, we have included the LDOS curves of the corresponding system without the magnetic field (bottom plots).
From the observation of these plots, it is clear that the magnetic field strongly affects the electronic and transport properties of the considered heterostructures, defining and modelling the electrical response of the conductor. In this sense, we have observed that in all considered systems, periodic metalsemiconductor electronic transitions for different values of magnetic flux ratio ϕ/ϕ_{0}, which are qualitatively in agreement with the experimental reports of similar heterosructures [2123]. Although the periodic electronic transitions are more evident in symmetric heterostructures (left and right panels), it is possible to obtain a similar effect in the asymmetric configurations. These behaviours are direct consequences of the quantum interference of the electronic wave function inside this kind of annular conductors, which in general present an AharonovBohm period as a function of the magnetic flux.
The evolution of the electronic levels of the system, depending of their energy, exhibits a rich variety of behaviours as a function of the external field. In all considered cases, the LDOS curves exhibit electronic states pinned at the Fermi Level, at certain magnetic flux values. This state corresponds to a nondispersive band, equivalent with the supersymmetric Landau level of the infinite twodimensional graphene crystal [30,35]. At low energy region and for low magnetic field, it is possible to observe the typical squareroot evolution of the relativistic Landau levels [36]. The electronic levels at highest energies of the system evolve linearly with the magnetic flux, like regular Landau levels. This kind of evolution is originated by the massive bands in graphene, which is expected for these kinds of states in graphenebased systems [37,38].
By comparing the LDOS curves and the corresponding conductance curves, it is possible to understand and define which states contribute to the transport of the systems (resonant tunneling peaks), and which ones only evolve with the magnetic flux but remain as localized states (quasibond states) of the conductor. These kind of behaviour has been reported before in similar systems [19,20]. This fact is more evident in the symmetric cases, where there are several states in the ranges ϕ/ϕ_{0} ∈ [0.1, 0.9] and E(γ_{0}) ∈ [1.0, 1.0] of the LDOS curves which evolve linearly with the magnetic flux, but are not reflected in the conductance curves. In fact, at these ranges, the conductance curves exhibit marked gaps with linear evolution as a function of the magnetic flux. For the asymmetric case, it is more difficult to define which states behave similarly; however, there are still some regions at which the conductance exhibits gaps with linear evolution as a function of the magnetic flux. All these electronic modulations could be useful to generate on/off switches in electronic devices, by changing in a controlled way the magnetic field intensity applied to the heterostructures. We have obtained these behaviours for different configurations of conductor, considering variations in length and widths of the finite ribbons and leads.
Conclusions
In this work, we have analysed the electronic and transport properties of a conductor composed of two parallel and finite AGNRs, connected to two semiinfinite lead, in the presence of an external perturbation. We have thought these systems as two parallel wires of an hypothetical circuit made of graphene, and we have studied the transport properties as a function of the separation and the geometry of these ‘wires’, considering the isolated case and the presence of an external magnetic field applied to the system. We have observed resonant tunneling behaviour as a function of the geometrical confinement and a complete AharonovBohm type of modulation as a function of the magnetic flux. These two behaviours are observed even when the two AGNRs have different widths, and consequently, different transverse electronic states. Besides, the magnetic field generates a periodic metalsemiconductor transition of the conductor, which can be used in electronics applications. We want to note that our results are valid only in low temperature limits and in the absence of strong disorder into the systems. In the case of nonzero temperature, it is expected that the resonances in the conductance curves will become broad and will gradually vanish at room temperature [20].
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
LR and JWG have worked equally in all results presented in this paper. Both authors read and approved the final manuscript.
Authors’ information
LR is a professor at the Physics Department, Technical University Federico Santa Maria, Valparaiso, Chile. JWG is a postdoctoral researcher at the International Iberian Nanotechnology Laboratory, Braga, Portugal.
Acknowledgements
Authors acknowledge the financial support of FONDECYT (grant no.: 11090212) and USMDGIP (grant no.: 11.12.17).
References

Novoselov KS, Geim AK, Morozov SV, Jiang D, Zhang Y, Dubonos SV, Grigorieva IV, Firsov AA: Electric field effect in atomically thin carbon films.
Science 2004, 306:666. PubMed Abstract  Publisher Full Text

Berger C, Song Z, Li T, Li X, Ogbazghi AY, Feng R, Dai Z, Marchenkov AN, Conrad EH, First PN, de Heer WA: Ultrathin epitaxial graphite: 2D electron gas properties and a route toward graphenebased nanoelectronics.
J Phys Chem B 2004, 108:19912. Publisher Full Text

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

Gomes KK, Mar W, Ko W, Guinea F, Manoharan HC: Designer Dirac fermions and topological phases in molecular graphene.
Nature 2012, 483:306. PubMed Abstract  Publisher Full Text

Li X, Wang X, Zhang L, Lee S, Dai H: Chemically derived, ultrasmooth graphene nanoribbon semiconductors.
Science 2008, 319:1229. PubMed Abstract  Publisher Full Text

Ci L, Xu Z, Wang L, Gao W, Ding F, Kelly KF, Yakobson BI, Ajayan PM: Controlled nanocutting of graphene.
Nano Res 2008, 1:116. Publisher Full Text

Kosynkin D, Higginbotham AL, Sinitskii A, Lomeda JR, Dimiev A, Price BK, Tour JM: Longitudinal unzipping of carbon nanotubes to form graphene nanoribbons.
Nature 2009, 458:872. PubMed Abstract  Publisher Full Text

Terrones M: Materials science: nanotubes unzipped.
Nature 2009, 458:845. PubMed Abstract  Publisher Full Text

Oezyilmaz B, JarilloHerrero P, Efetov D, Abanin D, Levitov LS, Kim P: Electronic transport and quantum Hall effect in bipolar graphene pnp junctions.
Phys Rev Lett 2007, 99:166804. PubMed Abstract  Publisher Full Text

Ponomarenko LA, Schedin F, Katsnelson MI, Yang R, Hill EW, Novoselov KS, Geim A: Chaotic Dirac billiard in graphene quantum dots.
Science 2008, 320:356. PubMed Abstract  Publisher Full Text

González JW, Santos H, Pacheco M, Chico L, Brey L: Electronic transport through bilayer graphene flakes.

Pedersen TG, Flindt C, Pedersen J, Mortensen N, Jauho A, Pedersen K: Graphene antidot lattices: designed defects and spin qubits.
Phys Rev Lett 2008, 100:136804. PubMed Abstract  Publisher Full Text

Oezyilmaz B, JarilloHerrero P, Efetov D, Kim P: Electronic transport in locally gated graphene nanoconstrictions.

Stankovich S, Dikin DA, Dommett GHB, Kohlhaas KM, Zimney EJ, Stach EA, Piner RD, Nguyen ST, Ruoff RS: Graphenebased composite materials.
Nature 2006, 442:282. PubMed Abstract  Publisher Full Text

Schedin F, Geim A, Morozov S, Hill E, Blake P, Katsnelson M, Novoselov K: Detection of individual gas molecules adsorbed on graphene.
Nat Mat 2007, 6:652. Publisher Full Text

Rosales L, Pacheco M, Barticevic Z, Latgé A, Orellana P: Transport properties of graphene nanoribbons with sideattached organic molecules.
Nanotechnology 2008, 19:065402. PubMed Abstract  Publisher Full Text

Rosales L, Pacheco M, Barticevic Z, Latgé A, Orellana P: Conductance gaps in graphene ribbons designed by molecular aggregations.
Nanotechnology 2009, 20:095705. PubMed Abstract  Publisher Full Text

Schurtenberger E, Molitor F, Gttinger J, Ihn T, Ensslin K: Tunable graphene single electron transistor.

Zhang ZZ, Wu ZH, Chang K, Peteers F M: Resonant tunneling through S and Ushaped graphene nanoribbons.
Nanotechnology 2009, 20:415203. PubMed Abstract  Publisher Full Text

Wu ZH, Zhang ZZ, Chang K, Peteers FM: Quantum tunneling through graphene nanorings.
Nanotechnology 2010, 21:185201. PubMed Abstract  Publisher Full Text

Smirnov D, Schmidt H, Haug RJ: AharonovBohm effect in an electronhole graphene ring system.
Appl Phys Lett 2012, 100:203114. Publisher Full Text

Russo S, Oostinga JB, Wehenkel D, Heersche HB, Sobhani SS, Vandersypen LMK, Morpurgo AF: Observation of AharonovBohm conductance oscillations in a graphene ring.

Huefner M, Molitor F, Jacobsen A, Pioda A, Stampfer C, Ensslin K, Ihn T: The AharonovBohm effect in a sidegated graphene ring.
New J Phys 2010, 12:043054. Publisher Full Text

Son YW, Cohen ML, Louie SG: Energy gaps in graphene nanoribbons.
Phys Rev Lett 2006, 97:216803. PubMed Abstract  Publisher Full Text

Nardelli M: Electronic transport in extended systems: application to carbon nanotubes.
Phys Rev B 1999, 60:7828. Publisher Full Text

Datta S: Electronic Transport Properties of Mersoscopic Systems. Cambridge: Cambridge University Press; 1995.

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. Publisher Full Text

Ritter C, Makler SS, Latgé A: Energygap modulations of graphene ribbons under external fields: a theoretical study.
Phys Rev B 2008, 77:195443.
A published erratum appears in Phys Rev B 2010, 82:089903(E)

Wakabayashi K, Fujita M, Ajiki H, Sigrist M: Electronic and magnetic properties of nanographite ribbons.
Phys Rev B 1999, 59:8271. Publisher Full Text

Nemec N, Cuniberti G: Hofstadter butterflies of carbon nanotubes: pseudofractality of the magnetoelectronic spectrum.

Rocha CG, Latgé A, Chico L: Metallic carbon nanotube quantum dots under magnetic fields.

Wakabayashi K: Electronic transport properties of nanographite ribbon junctions.

González JW, Rosales L, Pacheco M: Resonant states in heterostructures of graphene nanoribbons.

González JW, Pacheco M, Rosales L, Orellana PA: Transport properties of graphene quantum dots.

Nemec N, Cuniberti G: Surface physics, nanoscale physics, lowdimensional systemsHofstadter butterflies of bilayer graphene.

Zhang ZZ, Chang K, Peteers FM: Tuning of energy levels and optical properties of graphene quantum dots.

Nemec N: Quantum Transport in Carbonbased Nanostructures: Theory and Computational Methods. New York: Simon & Schuster; 2008.

Katsnelson M: Graphene: Carbon in Two Dimensions. Cambridge: Cambridge University Press; 2012.