in the last three decades transport through molecular junctions @xcite has attracted much attention because of fundamental aspects of the processes involved , as well as of potential practical prospects . in particular , studies in this direction
have intensified since the discovery of new forms of carbon allotropes , starting with the fullerenes @xcite and carbon nanotubes @xcite ( cnts ) in the 1980s and 1990s respectively , and the isolation of graphene @xcite in 2004 .
the above carbon allotropes differ in shape ( curvature ) , topology , and dimensionality , with the fullerenes being zero - dimensional ( 0d ) with spherical or prolate shape , the cnts being one - dimensional ( 1d ) cylinders , and graphene being a two - dimensional ( 2d ) plane ( or 1d planar ribbons @xcite ) . in the fullerenes
the carbon - atom network is made of ( non - adjacent ) hexagons and pentagons , whereas the cnts and graphene are entirely hexagonal lattices ( curved in cnts ) described in terms of a unit cell with a two - atom basis with the two carbon atoms occupying two sublattices ( @xmath0 and @xmath1 , also mapped into the up and down pseudospin states @xcite ) . in the absence of defects , in - plane ( @xmath2 )
bonding occurs through @xmath3 hybrid orbitals and out - of - plane bonding ( @xmath4 ) involves the @xmath5 orbital ; in the following only the physics of @xmath4-states is considered .
the hexagonal network topology of graphene gives rise to relativistic behavior of the low - energy excitations which is captured by the ultra - relativistic massless dirac - weyl ( dw ) equation , with the fermi velocity of graphene ( @xmath6 ) replacing the velocity of light .
@xcite among the many manifestations of the relativistic behavior in graphene is the klein paradox , that is `` ... unimpeded penetration of relativistic particles through high and wide potential barriers  is one of the most exotic and counterintuitive consequences of quantum electrodynamics.''@xcite the surprising relativistic behavior in graphene has indeed been recognized in the 2010 nobel award in physics to a. geim and k. novosolev .
another carbon - based system that was the subject of an earlier ( 2000 ) chemistry nobel award ( to a.j .
heeger , a.g .
macdiarmid , and h. shirakawa ) is conductive polymers , with polyacetelyne ( pac ) being a representative example .
@xcite pac is a 1d chain of carbon atoms forming a conjugated polymer with @xmath7 hybridization that leads to one unpaired electron per carbon atom ( half - filled @xmath4-band ) as in graphene . linearizing the spectrum of this 1d equally - spaced carbon chain at the fermi level ( that is at the dirac points , i.e. , the zero - energy points of the energy vs momentum dispersion relation ) results in a 1d massless dirac - equation description of the low energy excitations of the system
@xcite    however , the equally - spaced 1d system is unstable and consequently it distorts ( structurally ) spontaneously ( peierls distortion @xcite ) yielding a modulation ( alternation ) of the spacing between successive sites that results in dimerization of successive atoms along the chain and the opening of a gap in the electronic spectrum .
this dimerization can occur in two energetically degenerate , but spatially distinct , patterns , termed as equivalent domain structures . either of these domains is a realization of the 1d dirac equation with a constant mass term , @xmath8 .
disruption of the dimerization pattern ( e.g. , by transforming from one domain to another along the chain ) creates a domain wall which can be described in the 1d dirac formulation through the use of a position - dependent mass term of alternating sign ( with the spatial mass - profile connecting @xmath9 with @xmath10 ) .
the solution of this generalized 1d dirac equation is a soliton characterized by having zero energy and by being localized at the domain wall . in this paper , the dirac equation with position - dependent mass will be used in investigations of the electronic structure , transport , and magnetic - field - induced phenomena ( aharonov - bohm ) in graphene nanostructures such as nanoribbon junctions and rings .
@xcite    from the above we conclude that in the two extreme size - domains , that is a 2d infinite graphene sheet and a 1d carbon chain ( pac ) , the systems exhibit behavior that is relativistic in nature . in an attempt to bridge between these two size - domains ,
we briefly discuss in the following systems of successively larger width , starting from the polyacene ( a quasi-1d chain of fused benzene rings ) which may be regarded as the narrowest graphene nanoribbon ( gnr ) with zigzag edge terminations .
polyacene was investigated @xcite first in 1983 .
it was found that , as in graphene and pac ( without distortion ) , the valence and conduction bands of undistorted polyacene touch at the edge of the brillouin zone .
however , unlike 2d graphene and pac , the dispersion relation about the touching point is quadratic , conferring a non - relativistic ( schrdinger equation ) character .
we show in this paper that this surprising finding persists for sufficiently narrow gnrs with zigzag termination ( zgnrs ) . however , narrow armchair - terminated gnrs ( agnrs ) are found here to maintain relativistic behavior , with metallic ones being massless and semiconducting ones being massive ( both classes obeying the einstein energy relation ) .
it this paper , we discuss mainly manifestations of relativistic and/or nonrelativistic quantum behavior explored through theoretical considerations of transport measurements in segmented graphene nanoribbons of variable width , and spectral and topological effects in graphene nanorings in the presence of magnetic fields .
transport through narrow graphene channels @xmath11 particularly bottom - up fabricated and atomically - precise graphene nanoribbons @xcite @xmath11 is expected to offer ingress to unique behavior of dirac electrons in graphene nanostructures . in particular ,
the wave nature of elementary particles ( e.g. , electrons and photons ) is commonly manifested and demonstrated in transport processes .
because of an exceptionally high electron mobility and a long mean - free path , it has been anticipated that graphene @xcite devices hold the promise for the realization , measurement , and possible utilization of fundamental aspects of coherent and ballistic transport behavior , which to date have been observed , with varying degrees of success , mainly at semiconductor interfaces @xcite , quantum point contacts @xcite , metallic wires @xcite , and carbon nanotubes @xcite .
another manifestation of coherent ballistic transport are interference phenomena , reflecting the wave nature of the transporting physical object , and associated most often with optical ( electromagnetic waves , photons ) systems or with analogies to such systems ( that is , the behavior of massless particles , as in graphene sheets ) . measurements of interference patterns are commonly made with the use of interferometers , most familiar among them the multi - pass optical fabry - prot ( ofp ) interferometer @xcite .
the advent of 2d forms of carbon allotropes has motivated the study of optical - like interference phenomena associated with relativistic massless electrons , as in the case of metallic carbon nanotubes @xcite and graphene 2d @xmath12-@xmath13 junctions @xcite .
( we note that the hallmark of the ofp is that the energy separation between successive maxima of the interference pattern varies as the inverse of the cavity length @xmath14 . )    for gnrs with segments of different widths , our investigations reveal diverse fabry - prot transport modes beyond the ofp case , with conductance quantization steps ( @xmath15 , with @xmath16 ) found only for uniform gnrs .
in particular , three distinct categories of fabry - prot interference patterns are identified :    1 .   _
fp - a : _ an _ optical _ fp pattern corresponding to _ massless _ graphene electrons exhibiting _ equal spacing _ between neighboring peaks .
this pattern is associated with metallic armchair nanoribbon central segments .
this category is subdivided further to _ fp - a1 _ and _ fp - a2 _ depending on whether a valence - to - conduction gap is absent ( _ fp - a1 _ , associated with metallic armchair leads ) , or present ( _ fp - a2 _ , corresponding to semiconducting armchair leads ) .
2 .   _ fp - b : _
a _ massive relativistic _ fp pattern exhibiting a shift in the conduction onset due to the valence - to - conduction gap and _ unequal peak spacings_. this pattern is associated with semiconducting armchair nanoribbon central segments , irrespective of whether the leads are metallic armchair , semiconducting armchair , or zigzag .
_ fp - c : _ a _ massive non - relativistic _ fp pattern with @xmath17 _ peak spacings _ , but with a vanishing valence - to - conduction gap , @xmath14 being the length of the central segment .
this pattern is the one expected for usual semiconductors described by the ( nonrelativistic ) schrdinger equation , and it is associated with zigzag nanoribbon central segments , irrespective of whether zigzag or metallic armchair leads are used
.    we report in this paper on the unique apects of transport through segmented gnrs obtained from tight - binding non - equilibrium green s function @xcite ( tb - negf ) calculations in conjunction with an analysis based on a one - dimensional ( 1d ) relativistic dirac continuum model .
this continuum model goes beyond the physics of the massless dirac - weyl ( dw ) electron , familiar from two - dimensional ( 2d ) honeycomb carbon sheets @xcite , and it is referred to by us as the dirac - fabry - prot ( dfp ) theory ( see below for the choice of name ) . in particular , it is shown that the valence - to - conduction energy gap in armchair gnr ( agnr ) segments , as well as the barriers at the interfaces between nanoribbon segments , can be incorporated in an effective position - dependent mass term in the dirac hamiltonian ; the transport solutions associated with this hamiltonian exhibit conductance patterns comparable to those obtained from the microscopic negf calculations . for zigzag graphene nanoribbon ( zgnr ) segments , the valence - to - conduction energy gap vanishes , and the mass term is consonant with excitations corresponding to massive nonrelativistic schrdinger - type carriers .
the faithful reproduction of these unique tb - negf conductance patterns by the dfp theory , including mixed armchair - zigzag configurations ( where the carriers transit from a relativistic to a nonrelativistic regime ) , provides a unifying framework for analysis of coherent transport phenomena and for interpretation of experiments targeting fundamental understanding of transport in gnrs and the future development of graphene nanoelectronics .    to demonstrate the aforementioned soliton formation due to structural topological effects ( discussed by us above in the context of polyacetylene ) , we explore with numerical tight - binding calculations and a dirac - kronig - penney ( dkp ) approach , soliton formation and charge fractionization in graphene rhombic rings ; this dkp approach is based on a generalized dirac equation with alternating - sign position - dependent masses .    before leaving the introduction
, we mention that , due to their importance as fundamental phenomena , aharonov - bohm - type effects in graphene - nanoribbons systems have attracted ( in addition to refs .   ) considerable theoretical attention . @xcite these latter theoretical papers , however , based their analysis exclusively on tight - binding and/or dft calculations@xcite , or they used in addition a two - dimensional dirac equation with infinite - mass boundary conditions.@xcite transcending the level of current understanding which explores direct similaritiess with the aharonov - bohm physics in semiconductor and metallic mesoscopic rings , our work here analyzes the tb results in conjunction with a continuum 1d generalized dirac equation ( that incorporates a position - dependent mass term ) , and thus it enables investigations of until - now unexplored topological aspects and relativistic quantum - field analogies of the ab effect in graphene nanosystems .
furthermore we note that oscillations in the conductance of graphene nanoribbons in the presence of magnetic barriers were found in a theoretical study @xcite ( using exclusively a tb - negf approach ) , as well as in an experimental investigation of high quality bilayer nanoribbons , @xcite and they were attributed to fabry - prot - type interference . in the absence of a continuum
dirac analysis , the precise relation of such oscillatory patterns to our fabry - prot theory ( based on the incorporation of a position - dependent mass term in the dirac equation at zero magnetic field ) warrants further investigation .
finally , gap engineering in graphene ribbons under strong external fields was studied in ref .  .
we stress again that one of the main results in this paper is the appearance of `` forbidden '' solitonic states inside the energy gap in the context of the low - magnetic - field aharonov - bohm spectra of graphene nanorings ; see the part titled `` aharonov - bohm spectra of rhombic graphene rings '' in the results section .
the energy of a fermion ( with one - dimensional momentum @xmath18 ) is given by the einstein relativistic relation @xmath19 , where @xmath8 is the rest mass and @xmath20 is the fermi velocity of graphene . in a gapped uniform graphene nanoribbon
, the mass parameter is related to the particle - hole energy gap , @xmath21 , as @xmath22 . following the relativisitic quantum - field lagrangian formalism ,
the mass @xmath8 is replaced by a position - dependent scalar higgs field @xmath23 , to which the relativistic fermionic field @xmath24 couples through the yukawa lagrangian @xcite @xmath25 ( @xmath26 being a pauli matrix ) . for @xmath27 ( constant ) @xmath28 , and the massive fermion dirac theory
is recovered .
the dirac equation is generalized as ( here we do not consider applied electric or magnetic fields ) @xmath29 \psi + i \hbar v_f \alpha \frac{\partial \psi}{\partial x } - \beta \phi(x ) \psi=0 . \label{direq}\ ] ] in one dimension , the fermion field is a two - component spinor @xmath30 ; @xmath31 and @xmath32 stand , respectively , for the upper and lower component and @xmath33 and @xmath26 can be any two of the three pauli matrices .
note that the higgs field enters in the last term of eq .
( [ direq ] ) .
@xmath34 in the first term is the usual electrostatic potential , which is inoperative due to the klein phenomenon and thus is set to zero for the case of the armchair nanoribbons ( where the excitations are relativistic ) .
the generalized dirac eq .
( [ direq ] ) is used in the construction of the transfer matrices of the dirac - fabry - prot model described below .
the building block of the dfp model is a 2@xmath352 wave - function matrix @xmath36 formed by the components of two independent spinor solutions ( at a point @xmath37 ) of the onedimensional , first - order generalized dirac equation [ see eq .
( 3 ) in the main paper ] .
@xmath36 plays @xcite the role of the wronskian matrix @xmath38 used in the second - order nonrelativistic kronig - penney model .
following ref .   and
generalizing to @xmath39 regions , we use the simple form of @xmath36 in the dirac representation ( @xmath40 , @xmath41 ) , namely @xmath42 where @xmath43 the transfer matrix for a given region ( extending between two matching points @xmath44 and @xmath45 specifying the potential steps @xmath46 ) is the product @xmath47 ; this latter matrix depends only on the width @xmath48 of the region , and not separately on @xmath44 or @xmath45 .
the transfer matrix corresponding to a series of @xmath39 regions can be formed @xcite as the product @xmath49 where @xmath50 is the width of the @xmath51th region [ with @xmath52 .
the transfer matrix associated with the transmission of a free fermion of mass @xmath8 ( incoming from the right ) through the multiple mass barriers is the product @xmath53 with @xmath54 , @xmath55 ; for armchair leads @xmath56 , while for zigzag leads @xmath57 .
naturally , in the case of metallic armchair leads , @xmath58 , since @xmath59 .
then the transmission coefficient @xmath60 is @xmath61 while the reflection coefficient is given by @xmath62    at zero temperature , the conductance is given by @xmath63 ; @xmath60 is the transmission coefficient in eq .
( [ trans ] ) .
+      the transfer matrix corresponding to either half of the rhombus can be formed @xcite as the product @xmath64 with @xmath14 being the length of half of the perimeter of the rhombus ; @xmath65 , with @xmath66 and @xmath67 .
the transfer matrix associated with the complete unit cell ( encircling the rhombic ring ) is the product @xmath68    following refs .
, we consider the superlattice generated from the virtual periodic translation of the unit cell as a result of the application of a magnetic field @xmath1 perpendicular to the ring .
then the aharonov - bohm energy spectra are given as solutions of the dispersion relation @xmath69 = \operatorname{tr}[{\bf t}(e)]/2 , \label{disrel}\ ] ] where we have explicitly denoted the dependence of the r.h.s .
on the energy @xmath70 ; @xmath71 for a rhombus with type - i corners and @xmath72 for a rhombus with type - ii corners .    the energy spectra and single - particle densities do not depend on a specific representation .
however , the wave functions ( upper and lower spinor components of the fermionic field @xmath73 ) do depend on the representation used . to transform the initial dkp wave functions to the ( @xmath74 , @xmath75 ) representation , which corresponds to the natural separation of the tight - binding amplitudes into the @xmath0 and @xmath1 sublattices ,
we apply successively the unitary transformations @xmath76 and @xmath77 .
to describe the properties of graphene nanostructures in the tight - binding approximation , we use the hamiltonian @xmath78 with @xmath79 indicating summation over the nearest - neighbor sites @xmath80 .
@xmath81 ev is the hopping parameter of two - dimensional graphene .    to calculate the tb - negf transmission coefficients , the hamiltonian ( [ htb ] )
is employed in conjunction with the well known transport formalism which is based on the nonequilibrium green s functions @xcite .    according to the landauer theory , the linear conductance is @xmath82 , where the transmission coefficient is calculated as @xmath83 $ ] .
the green s function @xmath84 is given by @xmath85 with @xmath86 being the hamiltonian of the isolated device ( junction without the leads ) .
the self - energies @xmath87 are given by @xmath88^{-1}\tau_{l(r)}^\dagger$ ] , where the hopping matrices @xmath89 describe the left ( right ) device - to - lead coupling , and @xmath90 is the hamiltonian of the semi - infinite left ( right ) lead .
the broadening matrices are given by @xmath91 $ ] .
+   +   +   +   +
our results for a 3-segment ( @xmath92 , where @xmath93 is the lead width and @xmath94 is the width of the central segment ) _ all - semiconducting _ agnr , aaa ( sss ) , are portrayed in fig .
[ fig1 ] [ see schematic lattice diagrams in figs .
[ fig1](i ) and [ fig1](ii ) ] . a uniform semiconducting armchair gnr
[ see fig .
[ fig1](i ) ] exhibits ballistic quantized - conductance steps [ see fig .  [
fig1](a ) ] . in contrast
, conductance quantization is absent for a nonuniform 3-segment ( @xmath95 ) agnr ; see figs .
[ fig1](b ) @xmath11 [ fig1](d ) . here
oscillations appear instead of quantized steps .
the first oscillation appears at an energy @xmath96 [ fig .
[ fig1](b ) ] , which reflects the intrinsic gap @xmath97 of the semiconducting central segment belonging to the class ii of agnrs , specified @xcite by a width @xmath98 , @xmath99 .
we recall that as a function of their width , @xmath100 , the armchair graphene nanoribbons fall into three classes : ( i ) @xmath101 ( semiconducting , @xmath102 ) , ( ii ) @xmath103 ( semiconducting , @xmath102 ) , and ( iii ) @xmath104 ( metallic , @xmath105 ) , @xmath99 .    that the leads are semiconducting does not have any major effect .
this is due to the fact that @xmath106 , and as a result the energy gap @xmath107 of the central segment is larger than the energy gap @xmath108 of the semiconducting leads [ see schematic in fig .
[ fig1](iii ) ] . in the opposite case ( central segment wider than the leads )
, the energy gap of the semiconducting leads would have determined the onset of the conductance oscillations .
the armchair gnr case with interchanged widths ( i.e. , @xmath109 instead of @xmath95 ) is portrayed in figs .
[ fig1](e ) @xmath11 [ fig1](f ) . in this case
the energy gap of the semiconducting leads ( being the largest ) determines the onset of the conductance oscillations .
it is a testimonial of the consistency of our dfp method that it can reproduce [ see fig .  [
fig1](d ) and fig .
[ fig1](f ) ] both the @xmath95 and @xmath109 tb - negf conductances ; this is achieved with very similar sets of parameters taking into consideration the central - segment - leads interchange .
we note that the larger spacing between peaks ( and also the smaller number of peaks ) in the @xmath109 case is due to the smaller mass of the central segment ( @xmath110 instead of @xmath111 ) .
further insight can be gained by an analysis of the discrete energies associated with the humps of the conductance oscillations in fig .
[ fig1](c ) and the resonant spikes in fig .
[ fig1](e ) .
indeed a simplified approximation for the electron confinement in the continuum model consists in considering the graphene electrons as being trapped within a 1d infinite - mass square well ( imsw ) of length @xmath112 ( the mass terms are infinite outside the interval @xmath112 and the coupling to the leads vanishes ) .
the discrete spectrum of the electrons in this case is given @xcite by @xmath113 where the wave numbers @xmath114 are solutions of the transcendental equation @xmath115 when @xmath59 [ massless dirac - weyl ( dw ) electrons ] , one finds for the spectrum of the imsw model : @xmath116 with @xmath117 .
@xmath118 which is twice the energy @xmath119 of the lowest state .
as is well known , a constant energy separation of the intensity peaks , inversely proportional to the length of the resonating cavity [ here @xmath112 , see eqs .
( [ edis ] ) and ( [ eoff ] ) above ] is the hallmark of the optical fabry - prot , reflecting the linear energy dispersion of the photon in optics @xcite or a massless dw electron in graphene structures.@xcite most revealing is the energy offset away form zero of the first conductance peak , which equals exactly one - half of the constant energy separation between the peaks . in one - dimension
, this is the hallmark of a massless fermion subject to an infinite - mass - barrier confinement @xcite .
naturally , in the case of a semiconducting segment ( see below ) , this equidistant behavior and @xmath120-offset of the conductance peaks do not apply ; this case is accounted for by the dirac - fabry - prot model presented in methods , and it is more general than the optical fabry - prot theory associated with a photonic cavity .
@xcite    in the nonrelativistic limit , i.e. , when @xmath121 , one gets @xmath122 which yields the well known relations @xmath123 and @xmath124 for a massive relativistic electron , as is the case with the semiconducting agnrs in this paper , one has to numerically solve eq .
( [ eimsw ] ) and then substitute the corresponding value of @xmath114 in einstein s energy relation given by eq .
( [ erel ] ) .    from an inspection of fig .
[ fig1 ] , one can conclude that the physics associated with the all - semiconducting aaa junction is that of multiple reflections of a _ massive _ relativistic dirac fermion bouncing back and forth from the edges of a `` quantum box '' created by a double - mass barrier [ see the schematic of the double - mass barrier in figs .
[ fig1](iii ) and [ fig1](v ) ] .
in particular , to a good approximation the energies of the conductance oscillation peaks are given by the imsw eq.([eimsw ] ) with @xmath125 ( @xmath95 ) or @xmath126 ( @xmath109 ) .
in this respect , the separation energy between successive peaks in figs .  [ fig1](b ) , [ fig1](c ) , [ fig1](e ) , and [ fig1](f ) is not a constant , unlike the case of an all - metallic junction @xcite ( or a photonic cavity @xcite ) .    the conductance patterns in figs .
[ fig1](d ) and [ fig1](f ) correspond to the category _ fp - b _
( see the introductory section ) .
these patterns can not be accounted for by the optical fabry - prot theory , but they are well reproduced by the generalized dirac - fabry - prot model introduced by us in the methods section .          our results for a 3-segment ( @xmath127 ) _ semiconducting - metallic - semiconducting _ agnr , aaa ( sms ) , are portrayed in fig .
[ fig2 ] [ left , see schematic lattice diagram in fig .
[ fig2](i ) ] .
the first fp oscillation in the tb - negf conductance displayed in fig .
[ fig2](a ) appears at an energy @xmath96 , which reflects the intrinsic gap @xmath97 of the semiconducting leads ( with @xmath128 ) .
the energy spacing between the peaks in fig .
[ fig2](a ) is constant in agreement with the metallic ( massless dw electrons ) character of the central segment with @xmath129 .
the tb - negf pattern in fig .
[ fig2](a ) corresponds to the fabry - prot category _
fp - a2_. as seen from fig .
[ fig2](b ) , our generalized dirac - fabry - prot theory is again capable of faithfully reproducing this behavior .    a deeper understanding of the aaa ( sms ) case can be gained via an inspection of the density of states ( dos ) plotted in fig .
[ fig2](c ) for the total segmented agnr ( central segment plus leads ) and in fig .
[ fig2](d ) for the the isolated leads . in fig .
[ fig2](c ) , nine equidistant resonance lines are seen .
their energies are close to those resulting from the imsw eq .
( [ imsw0 ] ) ( with @xmath130 , see the caption of fig .
[ fig2 ] ) for a massless dw electron .
out of these nine resonances , the first five do not conduct [ compare figs .
[ fig2](a ) and [ fig2](c ) ] because their energies are lower than the minimum energy ( i.e. , @xmath131 ) of the incoming electrons in the leads [ see the onset of the first band ( marked by an arrow ) in the dos curve displayed in fig .  [
fig2](d ) ] .      as shown in refs .
, a detailed description of hygrogen passivation requires that the hopping parameters @xmath132 for the nearest - neighbor c - c bonds at the armchair edges be given by @xmath133 . taking this modification into account , our results for a 3-segment _ semiconducting - metallic - semiconducting _ agnr
are portrayed in fig .
[ fig2 ] [ right , see schematic lattice diagram in fig .
[ fig2](iii ) ] ; this lattice configuration is denoted as `` aaa ( sms ) h - passivation . ''
the first fp oscillation in the tb - negf conductance displayed in fig .
[ fig2](e ) appears at an energy @xmath134 , which reflects the intrinsic gap @xmath97 of the properly passivated semiconducting leads ( with @xmath128 ) .
the energy spacing between the peaks in fig .
[ fig2](e ) is slightly away from being constant in agreement with the small mass @xmath135 acquired by the central segment with @xmath129 , due to taking @xmath133 . as seen from fig .
[ fig2](f ) , our generalized dirac - fabry - prot theory is again capable of faithfully reproducing this behavior .    a deeper understanding of the aaa ( sms)-h - passivation case can be gained via an inspection of the dos plotted in fig .
[ fig2](g ) for the total segmented agnr ( central segment plus leads ) and in fig .
[ fig2](h ) for the isolated leads . in fig .
[ fig2](g ) , eight ( almost , but not exactly , equidistant ) resonance lines are seen .
their energies are close to those resulting from the imsw eq .
( [ eimsw ] ) ( with @xmath136 and @xmath135 ; see the caption of fig .
[ fig2 ] ) for a dirac electron with a small mass . out of these eight resonances , the first six
do not conduct [ compare figs .
[ fig2](e ) and [ fig2](g ) ] because their energies are lower than the minimum energy ( i.e. , @xmath137 ) of the incoming electrons in the leads [ see the onset of the first band ( marked by an arrow ) in the dos curve displayed in fig .
[ fig2](h ) ] . from the above we conclude that hydrogen passivation of the agnr resulted in a small shift of the location of the states , and opening of a small gap for the central metallic narrower ( with a width of @xmath138 ) segment , but did not modify the conductance record in any qualitative way .
moreover , the passivation effect can be faithfully captured by the dirac fp model by a small readjustment of the model parameters .    , with @xmath139 ( @xmath140 ) being the number of carbon atoms specifying the width of the ribbon segments .
the armchair leads in the aza junction are metallic ( @xmath141 , class iii agnr ) .
( a)-(b ) tb - negf conductance for the zzz and aza junction , respectively .
( c ) dfp conductance reproducing the tb - negf result in ( a ) for the zzz junction .
( d ) dfp conductance reproducing the tb - negf result in ( b ) for the aza junction . in spite of the different edge morphology , the fabry - prot patterns in ( a ) and ( b ) are very similar .
the central zigzag segment controls the fabry - prot patterns . according to the continuum dfp analysis , the physics underlying such patterns is that of a massive nonrelativistic schrdinger fermionic carrier performing multiple reflections within a cavity defined by a double - mass barrier [ see diagram in ( iii ) ] , but with the additional feature that @xmath142 and @xmath143 are also considered for segments or leads with zigzag edge terminations ( see text for details ) .
the mass and @xmath144 parameters used in the dfp calculations were @xmath145 , @xmath146 , with @xmath147 , @xmath142 , @xmath148 , @xmath149 , @xmath150 , @xmath151 , @xmath143 in ( c ) and @xmath152 , @xmath153 , with @xmath154 , @xmath142 , @xmath155 , @xmath156 , @xmath157 , @xmath158 , @xmath159 in ( d ) .
@xmath160 and @xmath161 denote parameters of the leads .
@xmath70 is the energy in units of @xmath162 .
@xmath163 nm is the graphene lattice constant ; @xmath81 ev is the hopping parameter.,width=434 ]      it is interesting to investigate the sensitivity of the interference features on the edge morphology .
we show in this section that the relativistic transport treatment applied to segmented armchaie gnrs does not maintain for the case of a nanoribbon segment with zigzag edge terminations .
in fact zigzag gnr ( zgnr ) segments exhibit properties akin to the well - known transport in usual semiconductors , i.e. , their excitations are governed by the nonrelativistic schrdinger equation .    before discussing segmented gnrs with
zigzag edge terminations , we remark that such gnrs with uniform width exhibit stepwise quantization of the conductance , similar to the case of a uniform armchair - edge - terminated gnr [ see fig .  [ fig1](a ) ] .    in fig .
[ fig3](a ) , we display the conductance in a three - segment junction [ see lattice schematic in fig .
[ fig3](i ) ] when all three segments have zigzag edge terminations ( denoted as zzz ) , but the central one is narrower than the lead segments .
the main finding is that the central segment behaves again as a resonant cavity that yields an oscillatory conductance pattern where the peak spacings are unequal [ fig .
[ fig3](a ) ] .
this feature , which deviates from the optical fabry - prot behavior , appeared also in the dfp patterns for a three - segment armchair junction whose central segment was semiconducting , albeit with a different dependence on @xmath14 [ see figs .
[ fig1 ] and eq .
( [ eimsw ] ) ] .
moreover , from a set of systematic calculations ( not shown ) employing different lengths and widths , we found that the energy of the resonant levels in zgnr segments varies on the average as @xmath164 , where the integer @xmath13 counts the resonances and @xmath14 indicates the length of the central segment .
however , a determining difference with the armchair gnr case in fig .
[ fig1 ] is the vanishing of the valence - to - conductance gap in the zigzag case of fig .
[ fig3](a ) .
it is well known that the above features are associated with resonant transport of electronic excitations that obey the nonrelativistic second - order schrdinger equation .
naturally , one could formulate a continuum transport theory based on transfer matrices ( see methods ) that use the 1d schrdinger equation instead of the generalized dirac eq .
( [ direq ] ) . such a schrdinger - equation continuum approach , however , is unable to describe mixed armchair - zigzag interfaces ( see below ) , where the electron transits between two extreme regimes , i.e. , an ultrarelativistic ( i.e. , including the limit of vanishing carrier mass ) dirac regime ( armchair segment ) and a nonrelativistic schrdinger regime ( zigzag segment ) .
we have thus been led to adopt the same dirac - type transfer - matrix approach as with the armchair gnrs , but with nonvanishing potentials @xmath57 .
this amounts to shifts ( in opposite senses ) of the energy scales for particle and hole excitations , respectively , and it yields the desired vanishing value for the valence - to - conduction gap of zigzag gnrs .
the calculated dfp conductance that reproduces well the tb - negf result for the zzz junction [ fig.[fig3](a ) ] is displayed in fig .
[ fig3](c ) ; the parameters used in the dfp calculation are given in the caption of fig .
we note that the carrier mass ( @xmath165 ) in the central zigzag segment exhibits an energy dependence .
this is similar to a well known effect ( due to nonparabolicity in the @xmath166 dispersion ) in the transport theory of usual semiconductors @xcite .
we further note that the average mass associated with a zigzag segment is an order of magnitude larger than that found for semiconducting armchair segments of similar width ( see caption in fig .
[ fig1 ] ) , and this yields energy levels @xmath164 close to the nonrelativistic limit [ see eq .
( [ ennr ] ) ] .
we note that the fp pattern of the zzz junction belongs to the category _ fp - c_.      fig .
[ fig3 ] ( right column ) presents an example of a mixed armchair - zigzag - armchair ( aza ) junction , where the central segment has again zigzag edge terminations [ see lattice schematic in fig .
[ fig3](ii ) ] .
the corresponding tb - negf conductance is displayed in fig .
[ fig3](b ) . in spite of the different morphology of the edges between the leads ( armchair ) and the central segment ( zigzag ) , the conductance profile of the aza junction [ fig .
[ fig3](b ) ] is very similar to that of the zzz junction [ fig .
[ fig3](a ) ] .
this means that the characteristics of the transport are determined mainly by the central segment , with the left and right leads , whether zigzag or armchair , acting as reservoirs supplying the impinging electrons .
the dfp result reproducing the tb - negf conductance is displayed in fig .
[ fig3](d ) , and the parameters used are given in the caption .
we stress that the mixed aza junction represents a rather unusual physical regime , where an ultrarelativistic dirac - weyl massless charge carrier ( due to the metallic armchair gnrs in the leads ) transits to a nonrelativistic massive schrdinger electron in the central segment .
we note that the fp pattern of the aza junction belongs to the _ fp - c _ caregory .    .
( b ) tb spectrum for a nanoring with type - ii corners and the same width @xmath167 .
these armchair graphene rings are semiconducting ( type - i ) and metallic ( type - ii ) .
the three ( four ) lowest - in - energy two - membered bands are shown . the hole states ( with @xmath168 , not shown ) are symmetric to the particle states ( with @xmath169 ) .
( c , d ) dkp spectra reproducing the tb ones in ( a ) and ( b ) , respectively .
insets in ( c ) and ( d ) : schematics of the higgs fields ( position - dependent mass ) @xmath170 employed in the dkp modeling .
@xmath170 is approximated by steplike functions @xmath46 ; @xmath51 counts the three regions of each half of the rhombus ( @xmath171 and @xmath172 ) , and @xmath13 ( @xmath173 ) .
the non - zero ( constant ) variable - mass values of @xmath170 are indicated by yellow ( red ) color when positive ( negative ) .
the parameters used in the dkp modeling are : ( c ) @xmath174 , @xmath175 , @xmath176 , @xmath177 [ see corresponding schematic inset in ( c ) ] and ( d ) @xmath178 , @xmath179 , @xmath176 , @xmath180 with @xmath181 [ see schematic inset in ( d ) ] .
note the two - membered braided bands and the `` forbidden '' band [ within the gap , in ( b ) and ( d ) ] .
the twofold forbidden band with @xmath182 appears as a straight line due to the very small amplitude oscillations of its two members .
@xmath163 nm is the graphene lattice constant and @xmath81 ev is the hopping parameter .
the edge terminations of both the inside and outside sides of the ring are armchair.,width=510 ]      the energy of a particle ( with onedimensional momentum @xmath18 ) is given by the einstein relativistic relation @xmath19 , where @xmath8 is the rest mass .
as aforementioned , in armchair graphene ribbons , the mass parameter is related to the particle - hole energy gap , @xmath21 , as @xmath22 . in relativistic quantum field theory ,
the mass of elementary particles is imparted through interaction with a scalar field known as the higgs field .
accordingly , the mass @xmath8 is replaced by a position - dependent higgs field @xmath183 , to which the relativistic fermionic field @xmath24 couples through the yukawa lagrangian @xcite @xmath25 ( @xmath26 being a pauli matrix ) . in the elementary - particles standard model ,
@xcite such coupling is responsible for the masses of quarks and leptons . for @xmath27 ( constant ) @xmath28 , and the massive fermion dirac theory
is recovered .
we exploit the generalized dirac physics governed by a total lagrangian density @xmath184 , where the fermionic part is given by @xmath185 and the scalar - field part has the form @xmath186 with the potential @xmath187 ( second term ) assumed to have a double - well @xmath188 form ; @xmath189 and @xmath190 are constants .
henceforth , the dirac equation ( see methods ) is generalized as @xmath191 in one dimension , the fermion field is a two - component spinor @xmath30 ; @xmath31 and @xmath32 stand , respectively , for the upper and lower component and @xmath33 and @xmath26 can be any two of the three pauli matrices .
-sublattice ( red ) and @xmath1-sublattice ( blue ) components of the tb state with energy @xmath192 at @xmath193 , belonging to the forbidden solitonic band of the type - ii nanoring with @xmath167 [ see fig .
[ fig44](b ) ] .
( b ) upper ( red ) and lower ( blue ) spinor components for the corresponding state ( forbidden band ) according to the dkp spectrum [ see fig .
[ fig44](d ) ] , reproducing the tb behavior of the type - ii nanoring with @xmath167 ( @xmath181 ) .
the tb and dkp wave functions for all states of the solitonic band are similar to those displayed here .
the wave functions here represent a pair of solitons .
for contrast , see fig .  10 in ref .   which portrays schematically the spinor @xmath194 for a _ single _ fermionic soliton attached to a higgs field with a smooth kink - soliton analytic shape @xmath195 .
@xmath196 is a solution @xcite of the lagrangian in eq .
( [ lagrphi ] ) .
dkp densities in units of @xmath197 . , width=472 ]    a graphene polygonal ring can be viewed as made of connected graphene - nanoribbon fragments ( here we consider agnrs ) .
the excitations of an infinite agnr are described by the 1d massive dirac equation , see eq .
( [ direq1 ] ) with @xmath74 , @xmath75 , and @xmath198 .
the two ( in general ) unequal hopping parameters @xmath199 and @xmath200 are associated with an effective 1d tight - binding problem ) and are given @xcite by @xmath201 $ ] , @xmath202 and @xmath203 ; @xmath100 is the number of carbon atoms specifying the width of the nanoribbon and @xmath81 ev is the hopping parameter for 2d graphene .
the effective @xcite tb hamiltonian of an agnr has a form similar to that used in _
trans_-polyacetylene ( a single chain of carbon atoms ) . in _
trans_-polyacetylene , the inequality of @xmath199 and @xmath200 ( referred to as dimerization ) is a consequence of the aforementioned peierls distortion induced by the electron - phonon coupling . for an armchair graphene ring ,
this inequality is a topological effect associated with the geometry of the edge and the width of the ribbon .
we recall that as a function of their width , @xmath100 , the armchair graphene nanoribbons fall into three classes : ( i ) @xmath204 ( semiconducting , @xmath205 ) , ( ii ) @xmath98 ( semiconducting , @xmath205 ) , and ( iii ) @xmath206 ( metallic @xmath207 ) , @xmath99 .
we adapt the `` crystal '' approach @xcite to the aharonov - bohm ( ab ) effect , and introduce a virtual dirac - kronig - penney @xcite ( dkp ) relativistic superlattice ( see methods ) .
charged fermions in a perpendicular magnetic field circulating around the ring behave like electrons in a spatially periodic structure ( period @xmath208 ) with the magnetic flux @xmath209 ( @xmath210 ) playing the role of the bloch wave vector @xmath211 , i.e. , @xmath212 [ see the cosine term in eq .
( [ disrel ] ) ] .
naturally , nanorings with arms made of nanoribbon segments belonging to the semiconducting classes may be expected to exhibit a particle - hole gap ( particle - antiparticle gap in rqf theory ) .
indeed this is found for a rhombic armchair graphene ring ( agr ) [ see gap @xmath21 in fig .
[ fig44](a ) ] with a width of @xmath167 carbon atoms having type - i corners .
suprisingly , a rhombic armchair graphene nanoring of the same width @xmath167 , but having corners of type - ii , demonstrates a different behavior , showing a `` forbidden '' band ( with @xmath182 ) in the middle of the gap region [ see fig .
[ fig44](b ) ] .
this behavior of rhombic armchair graphene rings with type - ii corners can be explained through analogies with rqf theoretical models , describing single zero - energy fermionic solitons with fractional charge @xcite or their modifications when forming soliton / anti - soliton systems .
@xcite ( a solution of the equation of motion corresponding to eq .
( [ lagrphi ] ) , is a @xmath213 kink soliton , @xmath196 .
the solution of eq .
( [ direq1 ] ) with @xmath214 is the fermionic soliton . )
we model the rhombic ring with the use of a continuous 1d kronig - penney @xcite model ( see methods ) based on the generalized dirac equation ( [ direq1 ] ) , allowing variation of the scalar field @xmath170 along the ring s arms .
we find that the dkp model reproduces [ see fig .
[ fig44](d ) ] the spectrum of the type - ii rhombic ring ( including the forbidden band ) when considering alternating masses @xmath215 associated with each half of the ring [ see inset in fig .  [ fig44](d ) ] .    in analogy with the physics of _ trans_-polyacetylene ( see remarks in the introductory section ) , the positive and negative masses correspond to two degenerate domains associated with the two possible dimerization patterns @xcite @xmath216 and @xmath217 , which are possible in a single - atom chain .
the transition zones between the two domains ( here two of the four corners of the rhombic ring ) are referred to as the domain walls .    for a single soliton ,
a ( precise ) zero - energy fermionic excitation emerges , localized at the domain wall . in the case of soliton - antisoliton pairs ,
paired energy levels with small positive and negative values appear within the gap . the tb spectrum in fig .
[ fig44](b ) exhibits a forbidden band of two paired @xmath218 levels , a property fully reproduced by the dkp model that employs two alternating mass domains [ fig .
[ fig44](d ) .
the twofold forbidden band with @xmath182 appears as a straight line due to the very small amplitude oscillations of its two members .
@xmath163 nm is the graphene lattice constant and @xmath81 ev is the hopping parameter .    the strong localization of a fraction of a fermion at the domain walls ( two of the rhombus corners ) , characteristic of fermionic solitons @xcite and of soliton / anti - soliton pairs ,
@xcite is clearly seen in the tb density distributions ( modulus of single - particle wave functions ) displayed in fig .
[ fig55](a ) .
the tb @xmath0 ( @xmath1 ) sublattice component of the tight - binding wave functions localizes at the odd numbered corners .
these alternating localization patterns are faithfully reproduced [ see fig .
[ fig55](b ) ] by the upper , @xmath219 , and lower , @xmath220 , spinor components of the continuum dkp model .
the soliton - antisoliton pair in fig .
[ fig55](b ) generates an @xmath221 charge fractionization at each of the odd - numbered corners , which is similar to the @xmath221 fractionization familiar from polyacetylene .
the absence of a forbidden band ( i.e. , solitonic excitations within the gap ) in the spectrum of the type - i rhombic nanorings [ see fig .
[ fig44](a ) ] indicates that the corners in this case do not induce an alternation between the two equivalent dimerized domains ( represented by @xmath215 in the dkp model ) . here
the corners do not act as topological domain walls .
nevertheless , direct correspondence between the tb and dkp spectra is achieved here too by using a variable higgs field defined as @xmath222 with @xmath223 and @xmath177 [ see the schematic inset in fig .  [ fig44](c ) and the dkp spectrum plotted in the same figure ] .
in this paper we focused on manifestations of relativistic and/or nonrelativistic quantum behavior explored through theoretical considerations of transport in graphene nanostructures and spectral and topological effects in graphene nanorings in the presence of magnetic fields .
in particular , we investigated the emergence of new behavior of electrons in atomically precise segmented graphene nanoribbons ( gnrs ) with different edge terminations ( armchair , zigzag and mixed ones ) , and in graphene rings . to these aims
we have employed tight - binding calculations of electronic states with and without applied magnetic fields , the non - equilibrium green s function transport theory , and a newly developed dirac continuum model that absorbs the valence - to - conductance energy gaps as position - dependent masses , including topological - in - origin mass - barriers at the contacts between segments .
the electronic conductance has been found to exhibit fabry - prot oscillations , or resonant tunneling , associated with partial confinement and formation of a quantum box ( resonant cavity ) in the junction . along with the familiar optical fp oscillations , exhibiting equal spacing between neighboring peaks , that we find for massless electrons in gnrs with metallic armchair central segments , we find other fp categories that differ from the optical one . in particular ,
our calculations reveal : ( a ) a massive relativistic fp pattern exhibiting a valence - to - conduction gap and unequal peak spacings .
this pattern is associated with semiconducting armchair nanoribbon central segments , irrespective of whether the armchair leads are matallic or semiconducting .
( b ) a massive non - relativistic fp pattern with @xmath17 peak spacings , but with a vanishing valence - to - conduction gap .
this pattern is the one expected for the carriers in usual semiconductors described by the ( nonrelativistic ) schrdinger equation , and it is associated with zigzag nanoribbon central segments , regardless of whether zigzag or metallic armchair leads are used .
perfect quantized - conductance flat steps were found only for uniform gnrs . in the absence of extraneous factors , like disorder ,
in our theoretical model , the deviations from the perfect quantized - conductance steps were unexpected .
however , this aforementioned behavior obtained through tb - negf calculations is well accounted for by a 1d contimuum fermionic dirac - fabry - prot interference theory ( see methods ) .
this approach employs an effective position - dependent mass term in the dirac hamiltonian to absorb the finite - width ( valence - to - conduction ) gap in armchair nanoribbon segments , as well as the barriers at the interfaces between nanoribbon segments forming a junction . for zigzag nanoribbon segments the mass term in the dirac equation
reflects the nonrelativistic schrodinger - type behavior of the excitations .
the carrier mass in zigzag - terminated gnr segments is much larger than the particle mass in semiconducting armchair - terminated gnr segments .
furthermore in the zigzag gnr segments ( which are always characterized by a vanishing valence - to - conduction energy gap ) , the mass corresponds simply to the carrier mass . in the armchair gnr segments , the carrier mass endows ( in addition ) the segment with a valence - to - conduction energy gap , according to einstein s relativistic energy relation [ see eq .
( 1 ) ] .
we concluded with a brief discussion of the physics of electrons in segmented polygonal rings , which may be regarded as constructed by connecting gnr segments .
evaluation of the electronic states in a rhombic graphene nanoring under the influence of an applied magnetic field in the aharonov - bohm regime , and their analysis with the use of a relativistic quantum - field theoretical model , unveils development of a topological - in - origin zero - energy soliton state and charge fractionization .
the above findings point to a most fundamental underlying physics , namely that the topology of disruptions of the regular honeycomb lattice ( e.g. , variable width segments , corners , edges ) generates a scalar - potential field ( position - dependent mass , identified @xcite also as a higgs - type field ) , which when integrated into a generalized dirac equation for the electrons provides a unifying framework for the analysis of transport processes through graphene segmented junctions and the nature of electronic states in graphene nanorings .    with growing activities and further improvements in the areas of bottom - up fabrication and manipulation of atomically precise @xcite graphene nanostructures and the anticipated measurement of conductance through them
, the above findings could serve as impetus and implements aiding the design and interpretation of future experiments .
999 nitzan , a. ; ratner , m. a. electron transport in molecular wire junctions .
_ science _ * 2003 * , _ 300 _ , 1384 - 1389 .
joachim , c. ; ratner , m. a. molecular electronics : some views on transport junctions and beyond . _ proc .
( usa ) _ * 2005 * , _ 102 _ , 8801 - 8808 .
lindsay , s. m. ; ratner , m. a. molecular transport junctions : clearing mists .
. maters . _ * 2007 * , _ 19 _ , 23 - 31 .
hod , o. ; rabani , e. ; baer , r. magnetoresistance of nanoscale molecular devices .
_ * 2006 * , _ 39 _ , 109 - 117 .
bergfield , j. p. ; ratner , m. a. forty years of molecular electronics : non - equilibrium heat and charge transport at the nanoscale . _ phys .
solidi b _ * 2013 * , _ 250 _ , 2249 - 2266 .
rai , d. ; hod , o. ; nitzan , a. magnetic field control of the current through molecular ring junctions .
_ j. phys .
lett . _ * 2011 * , _ 2 _ , 2118 - 2124 .
kroto , h. w. ; heath , j. r. ; obrien , s. c. ; curl r. f. ; smalley , r. e. c60 : buckminsterfullerine .
_ nature _ * 1985 * , _ 318 _ , 162 - 163 .
iijima , s. helical microtubules of graphitic carbon .
_ nature _ * 1991 * , _ 354 _ , 56 - 58 .
novoselov , k. s. ; geim , a. k. ; morozov , s. v. ; jiang , d. ; zhang , y. ; dubonos , s. v. ; grigorieva , i. v. ; firsov , a. a. electric field effect in atomically thin carbon films .
_ science _ * 2004 * , _ 306 _ , 666 - 669 .
nakada , k. ; fujita , m. ; dresselhaus , g. ; dresselhaus , m. s. edge state in graphene ribbons : nanometer size effect and edge shape dependence .
b : condens . matter mater .
phys . _ * 1996 * , _ 54 _ , 17954 - 17961 .
son , y. w. ; cohen , m. l. ; louie , s. g. energy gaps in graphene nanoribbons .
_ * 2006 * , _ 97 _ , 216803 .
wang , z. f. ; li , q. ; zheng , h. ; ren , h. ; su , h. ; shi , q. w. ; chen , j. tuning the electronic structure of graphene nanoribbons through chemical edge modification : a theoretical study . _
phys . rev .
b : condens . matter mater .
_ * 2007 * , _ 75 _ , 113406 .
wakabayashi , k. ; sasaki , k. ; nakanishi , t ; enoki , t. electronic states of graphene nanoribbons and analytical solutions .
technol . adv .
mater . _ * 2010 * , _ 11 _ , 054504 , and references therein .
castro neto , a. h. ; guinea , f. ; peres , n. m. r. ; novoselov , k. s. ; geim , a. k. the electronic properties of graphene .
phys . _ * 2009 * , _ 81 _ , 109 - 162 .
katsnelson , m. i. ; novoselov , k. s. ; geim , a. k. chiral tunnelling and the klein paradox in graphene .
_ nature phys .
_ * 2006 * , _ 2 _ , 620 - 625 .
heeger , a. j. ; kivelson , s. ; schrieffer , j. r. ; su , w .- p .
solitons in conducting polymers .
phys . _ * 1988 * , _ 60 _ , 781 - 847 .
jackiw , r. ; schrieffer , j. r. solitons with fermion number @xmath120 in condensed matter and relativistic field theories .
b _ * 1981 * , _ 190 [ fs3 ] _ , 253 - 265 .
jackiw , r. fractional and majorana fermions : the physics of zero - energy modes . _ phys .
scr . _ * 2012 * , _
t146 _ , 014005 .
peierls , r. f. _ quantum theory of solids _ ; clarendon press : oxford , * 1955*. romanovsky , i. ; yannouleas , c. ; landman , u. patterns of the aharonov - bohm oscillations in graphene nanorings . _
b : condens
. matter mater .
* 2012 * , _ 85 _ , 165434 .
romanovsky , i. ; yannouleas , c. ; landman , u. topological effects and particle physics analogies beyond the massless dirac - weyl fermion in graphene nanorings .
b : condens . matter mater .
phys . _ * 2013 * , _ 87 _ , 165431 .
yannouleas , c. ; romanovsky , i. ; landman , u. beyond the constant - mass dirac physics : solitons , charge fractionization , and the emergence of topological insulators in graphene rings .
b : condens . matter mater .
* 2014 * , _ 89 _ , 035432 .
kivelson , s. ; chapman , o. i. polyacene and a new class of quasi - one - dimensional conductors .
b : condens . matter mater .
phys . _ * 1983 * , _ 28 _ , 7236 - 7243 .
cai , j. m. ; ruffieux , p. ; jaafar , r. ; bieri , m. ; braun , t. ; blankenburg , s. ; muoth , m. ; seitsonen , a. p. ; saleh , m. ; feng , x .-
l . ; et al .
atomically precise bottom - up fabrication of graphene nanoribbons .
_ nature _ * 2010 * , _ 466 _ , 470 - 473 .
m. s. graphene : ribbons piece - by - piece .
mat . _ * 2010 * , _ 9 _ , 611 - 612 .
huang , h. ; wei , d. ; sun , j. ; wong , s. l. ; feng , y. p. ; castro neto , a. h. ; wee , a. t. s. spatially resolved electronic structures of atomically precise armchair graphene nanoribbons . _
rep . _ * 2012 * , _ 2 _ , 983 .
van der lit , j. ; boneschanscher , m. p. ; vanmaekelbergh , d. ; ijs , m. ; uppstu , a. ; ervasti , m. ; harju , a. ; liljeroth , p. ; swart , i. suppression of electron - vibron coupling in graphene nanoribbons contacted via a single atom .
commun . _ * 2013 * , _ 4 _ , 2023 .
narita , a. ; feng , x .-
; hernandez , y. ; jensen , s. a. ; bonn , m. ; yang , h. f. ; verzhbitskiy , i. a. ; casiraghi , c. ; hansen , m. r. ; koch , a. h. r. ; et al . synthesis of structurally well - defined and liquid - phase - processable graphene nanoribbons .
chem . _ * 2014 * , _ 6 _ , 126 - 132 .
hartley , c. s. graphene synthesis : nanoribbons from the bottom - up .
chem . _ * 2014 * , _ 6 _ , 91 - 92 .
vo , t. h. ; shekhirev , m. ; kunkel , d. a. ; morton , m. d. ; berglund , e. ; kong , l. ; wilson , p. m. ; dowben , p. a. ; enders , a. ; sinitskii , a. large - scale solution synthesis of narrow graphene nanoribbons .
commun . _ * 2014 * , _ 5 _ , 3189 .
blankenburg , s. ; cai j. m. ; ruffieux , p. ; jaafar , r. ; passerone , d. ; feng , x .-
l . ; mllen , k. ; fasel , r. ; pignedoli , c. a. intraribbon heterojunction formation in ultranarrow graphene nanoribbons .
_ acs nano _ _ 2012 _ , _ 6 _ , 2020 - 2025 .
van houten , h. ; beenakker , c. w. j. ; williamson , j. g. ; broekaart , m. e. i. ; van loosdrecht , p. h. m. coherent electron focusing with quantum point contacts in a two - dimensional electron gas .
b : condens . matter mater .
_ * 1989 * , _ 39 _ , 8556 - 8575 .
ji , y. ; chung , y. c. ; sprinzak , d. ; heiblum , m. ; mahalu , d. ; shtrikman , h. an electronic mach - zehnder interferometer .
_ nature _ * 2003 * , _ 422 _ , 415 - 418 .
vanwees , b. j. ; van houten , h. ; beenakker , c. w. j. ; williamson , j. g. ; kouwenhoven , l. p. ; van der marel , d. ; foxon , c. t. quantized conductance of point contacts in a two - dimensional electron gas .
lett . _ * 1988 * , _ 60 _ , 848 - 850 . pascual , j. i. ; mndez , j. ; gmez - herrero , j. ; bar , a. m. ; garcia , n. ; landman , u. ; luedtke , w. d. ; bogachek , e. n. ; cheng , h .- p .
properties of metallic nanowires : from conductance quantization to localization .
_ science _ * 1995 * , _ 267 _ , 1793 - 1795 .
liang , w. ; bockrath , m. ; bozovic , d. ; hafner , j. h. ; tinkham , m. ; park , h. fabry - prot interference in a nanotube electron waveguide .
_ nature _ * 2001 * , _ 411 _ , 665 - 669 .
vaughan , j.m .
; _ the fabry - prot interferometer _ ; hilger : bristol , * 1989*. lipson , s. g. ; lipson , h. ; tannhauser , d. s. _ optical physics _ ; cambridge univ . press : london , * 1995 * , 3rd ed . ;
rickhaus , p. ; maurand , r. ; liu , m .- h . ; weiss , m. ; richter , k. ; schnenberger , ch .
ballistic interferences in suspended graphene .
commun . _ * 2013 * , _ 4 _ , 2342 .
datta , s. ; _ quantum transport : atom to transistor _ ; cambridge university press : cambridge , * 2005*. lake , r. ; datta , s. nonequilibrium greens - function method applied to double - barrier resonant - tunneling diodes .
b : condens . matter mater .
phys . _ * 1992 * , _ 45 _ , 6670 - 6685 .
ezawa , m. peculiar band gap structure of graphene nanoribbons .
( c ) _ * 2007 * , _ 4 _ , 489 - 492 .
recher , p. ; trauzettel , b. ; rycerz , a. ; blanter , ya .
m. ; beenakker , c. w. j. ; morpurgo , a. f. aharonov - bohm effect and broken valley degeneracy in graphene rings .
b : condens . matter mater .
phys . _ * 2007 * , _ 76 _ , 235404 .
bahamon , d. a. ; pereira , a. l. c. ; schulz , p. a. inner and outer edge states in graphene rings : a numerical investigation .
b : condens . matter mater .
phys . _ * 2009 * , _ 79 _ , 125414 .
wurm , j. ; wimmer , m. ; baranger , h. u. ; richter , k. graphene rings in magnetic fields : aharonov - bohm effect and valley splitting .
_ semicond .
technol . _ * 2010 * , _ 25 _ , 034003 . hung nguyen , v. ; niquet , y. m. ; dollfus , p. aharonov - bohm effect and giant magnetoresistance in graphene nanoribbon rings .
b : condens . matter mater .
phys . _ * 2013 * , _ 88 _ , 035408 . hung nguyen , v. ; niquet , y. m. ; dollfus , p. the interplay between the aharonov - bohm interference and parity selective tunneling in graphene nanoribbon rings .
_ j. phys . : condens
. matter _ * 2014 * , _ 26 _ , 205301 .
xu , h. ; heinzel , t. ; evaldsson , m. ; zozoulenko , i. v. magnetic barriers in graphene nanoribbons : theoretical study of transport properties .
_ j. phys . : condens
. matter _ * 2008 * , _ 77 _ , 245401 .
jiao , l. ; wang , x. ; diankov , g. ; wang , h. ; dai , h. facile synthesis of high - quality graphene nanoribbons .
nanotech . _ * 2010 * , _ 5 _ , 321 - 325 .
ritter , c. ; makler , s. s. ; latg , a. energy - gap modulations of graphene ribbons under external fields : a theoretical study .
_ j. phys . : condens
. matter _ * 2008 * , _ 77 _ , 195443 .
mckellar , b. h. j. ; stephenson , jr . ,
g. j. relativistic quarks in one - dimensional periodic structures .
c _ * 1987 * , _ 35 _ , 2262 - 2271 .
bttiker , m. ; imry , y. ; landauer , r. josephson behavior in small normal one - dimensional rings .
* 1983 * , _ 96 _ , 365 - 367 .
alberto , p. ; fiolhais , c. ; gil , v. m. s. relativistic particle in a box .
j. phys . _ * 1996 * , _ 17 _ , 19 - 24 .
lpez - villanueva , j. a. ; melchor , i. ; cartujo , p. ; carceller , j. e. modified schrdinger - equation including nonparabolicity for the study of a two - dimensional electron gas .
b : condens . matter mater .
phys . _ * 1993 * , _ 48 _ , 1626 - 1631 .
mackenzie , r. ; palmer , w. f. bags via bosonization .
d _ * 1990 * , _ 42 _ , 701 - 709 .
bednyakov , v. a. ; giokaris , n. d. ; bednyakov , a. v. on higgs mass generation mechanism in the standard model .
nucl . _ * 2008 * , _ 39 _ , 13 - 36 ; arxiv : hep - ph/0703280 .
griffith , d. ; _ introduction to elementary particles _ ; wiley - vch : weinheim , germany , * 2008 * ; 2nd edition .
zheng , h. x. ; wang , z. f. ; luo , t. ; shi , q. w. ; chen , j. analytical study of electronic structure in armchair graphene nanoribbons . _ phys . rev .
b : condens . matter mater .
phys . _ * 2007 * , _ 75 _ , 165414 .
jackiw , r. ; rebbi , c. solitons with fermion number @xmath120 .
d _ * 1976 * , _ 13 _ , 3398 - 3409 .
jackiw , r. ; kerman , a. k. ; klebanov , i. ; semenoff , g. fluctuations of fractional charge in soliton anti - soliton systems .
b _ * 1983 * , _ 225 _ , 233 - 246 .
+ * acknowledgements : * this work was supported by a grant from the office of basic energy sciences of the us department of energy under contract no .
fg05 - 86er45234 .
computations were made at the gatech center for computational materials science .
+   + * author contributions : * i.r . & c.y .
performed the computations .
c.y . , i.r . & u.l .
analyzed the results .
& u.l . wrote the manuscript .
+   + * competing financial interests : * the authors declare no competing financial interests .
+   + * correspondence * should be addressed to u.l .
: ( uzi.landman@physics.gatech.edu ) .