Next Article in Journal
Activated Carbon/ZnFe2O4 Nanocomposite Adsorbent for Efficient Removal of Crystal Violet Cationic Dye from Aqueous Solutions
Next Article in Special Issue
Vanadium-Doped FeBP Microsphere Croissant for Significantly Enhanced Bi-Functional HER and OER Electrocatalyst
Previous Article in Journal
Optically-Induced Symmetry Switching in a Reconfigurable Kagome Photonic Lattice: From Flatband to Type-III Dirac Cones
Previous Article in Special Issue
Simulation of Figures of Merit for Barristor Based on Graphene/Insulator Junction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Independent Dual-Channel Approach to Mesoscopic Graphene Transistors

1
Instituto de Investigaciones en Materiales, Universidad Nacional Autónoma de México, Mexico City 04510, Mexico
2
Departamento de Física, Facultad de Ciencias, Universidad Nacional Autónoma de México, Mexico City 04510, Mexico
*
Author to whom correspondence should be addressed.
Nanomaterials 2022, 12(18), 3223; https://doi.org/10.3390/nano12183223
Submission received: 22 August 2022 / Revised: 9 September 2022 / Accepted: 13 September 2022 / Published: 16 September 2022
(This article belongs to the Special Issue Nanotechnologies and Nanomaterials: Selected Papers from CCMR)

Abstract

:
Graphene field-effect transistors (GFETs) exhibit unique switch and sensing features. In this article, GFETs are investigated within the tight-binding formalism, including quantum capacitance correction, where the graphene ribbons with reconstructed armchair edges are mapped into a set of independent dual channels through a unitary transformation. A new transfer matrix method is further developed to analyze the electron transport in each dual channel under a back gate voltage, while the electronic density of states of graphene ribbons with transversal dislocations are calculated using the retarded Green’s function and a novel real-space renormalization method. The Landauer electrical conductance obtained from these transfer matrices was confirmed by the Kubo–Greenwood formula, and the numerical results for the limiting cases were verified on the basis of analytical results. Finally, the size- and gate-voltage-dependent source-drain currents in GFETs are calculated, whose results are compared with the experimental data.

1. Introduction

Graphene is a one-atom-thick planar sheet of carbon atoms in a honeycomb structure, whose free-standing monolayer was isolated and characterized for the first time in 2004 by Andre Geim and Konstantin Novoselov [1]. This material possesses outstanding properties such as high specific surface area [2], very high carrier mobility with null effective mass [3], huge thermal conductivity, and giant Young’s modulus [4], making it an ideal candidate for a wide range of applications [5].
When the semiconductor layer of a traditional metal–oxide–semiconductor field-effect transistor (MOSFET) is replaced by a graphene sheet, it is called as graphene field-effect transistor (GFET), whose source-drain current is extremely sensitive to the adsorbed molecules, and it is employable for use in gas and bio sensing [6,7,8]. Moreover, GFET constitutes a central component of far-infrared radiation and terahertz devices [9] since graphene exhibits a quick optical response based on high carrier mobility [10]. In addition, the graphene nanoribbon field-effect transistor (GNRFET) presents an elevated I o n / I o f f ratio around the threshold voltage compared to the conventional Si-based CMOS [11]. Among possible disadvantages of GFETs, there are difficulties in fabricating wafer-sized high-quality graphene, as well as the apparent degradation of carrier mobility when it is placed over an oxide [12].
Additionally, GFET has a singular capacitor behavior, since the graphene ribbon, being one of the parallel plates, possesses a very low charge density around the Dirac point, which may induce large variations in its chemical potential when a small amount of electric charge is introduced into the capacitor. This fact could significantly increase the potential difference between the back gate and the graphene, because it is equal to the difference between the chemical potentials at each plate. Such small quantities of charge generating a large potential difference lead to a lower capacitance compared to geometric or classic ones ( C C ) obtained using ideal metal plates. Hence, the capacitance of GFET ( C t o t ) can be calculated by including a so-called quantum capacitance ( C Q ) through C t o t 1 = C C 1 + C Q 1 [13].
On the theoretical side, GFETs can be modeled at the atomic scale by means of first-principle or semiempirical methods. For example, N/B/P co-doped seven-atom-wide armchair GNRFET was analyzed by using the density functional theory (DFT) [14], while ten-atom-width GNRFET has been studied by self-consistently solving Schrödinger–Poison equations within the tight-binding formalism [15]. Recently, mesoscopic GFETs have been extensively investigated for the detection of micron-scale biomolecules [16] and their graphene ribbons contain several billions of atoms, which impede atomic-scale first-principle modeling. In addition, such graphene ribbons usually possess numerous structural defects, which requires an innovative real-space approach [17].
In this article, the source-drain current ( I s d ) in mesoscopic GFETs is studied by means of an independent dual-channel method that transforms the graphene ribbon with edge reconstruction defects and transversal dislocations into a set of double chains with first and second neighbor interactions within the tight-binding formalism. The electrical conductance and electronic density of states in these double chains or dual channels connected to two semi-infinite periodic leads at their ends are further addressed by means of a new transfer matrix method within the Landauer formalism and a novel real-space renormalization procedure based on the retarded Green’s function, respectively. The theoretical I s d , as a function of the gate voltage, is compared with experimental data obtained from GFETs of different lengths.

2. The Model

In GFETs, an extra carrier density ( n 2 D ) can be induced in graphene ribbon by the application of a gate voltage ( V G ), which is related to the chemical potentials of graphene ( μ ) and of back gate metal ( μ M ) through e V G = μ μ M . Figure 1 shows a schematic representation of (Figure 1a) top and (Figure 1b,c) side views of a backgated GFET with energy band diagrams for (b) V G = 0 and (c) V G > 0 , where the diminution of μ M due to V G is denoted as e V C . Hence, we have:
e V G = μ μ M = μ + e V C ,
which is equivalent to C t o t 1 = C C 1 + C Q 1 , since C t o t = e   n / V G , C C = e   n / V C , C Q = e   n / V Q and e   V Q = μ μ 0 with μ 0 = 0 or at the Dirac point in our case [18].
In general, the capacitance per unit of area (C) of an ideal parallel-plate capacitor with a dielectric material of thickness d and a relative permittivity ε r is C = ε r ε 0 / d = Q / V C in units of MKS [19], where ε 0 is the vacuum permittivity, Q = e n 2 D is the electric charge per unit of area in the capacitor and V C is the potential difference between the parallel plates. Hence:
e V C = e 2 n 2 D d ε r ε 0 ,
where e is the magnitude of the electron charge.
Let us consider an armchair-edged graphene ribbon with an even number of atoms (2W) in each transversal zig-zag line, as illustrated in Figure 2a. This ribbon of width W   = 3 a S + 3 ( 2 W 3 )   a 0 / 2 and length L = ( 3 / 2 ) L a 0 a 0 is found in a backgated GFET, where W and L are integer numbers, a 0 = 1.41 Å and a S = 1.36 Å are bond lengths between interior atoms and between the edge and next-edge atoms, respectively denoted by blue and thick blue lines in Figure 2a. This variation in interatomic distance is originated from the edge reconstruction of graphene ribbons [20]. The single-electron tight-binding Hamiltonian (H) of this graphene ribbon with null on-site energies can be written as
H = l   =   1 L t s l , 1 l , 2 + t j   =   2 2 W 2 l , j l , j + 1 + t s l , 2 W 1 l , 2 W + l   =   1 ( L 1 ) / 2 j   =   1 W t 2 l 1 2 l 1 , 2 j 2 l , 2 j + t 2 l 2 l , 2 j 1 2 l + 1 , 2 j 1 ,
where l , j is the Wannier function at site ( l , j ) , t is the hopping integral of pristine graphene sheet and t s   = 1.2   t is that between the edge and next-edge atoms in a graphene ribbon with edge reconstructions [20] and t l with l = 1 , 2 , , L 1 is the hopping integral between transversal zig-zag lines l and l + 1 , representing possible transversal dislocations when t l   t .
There is a unitary transformation discussed in Appendix A that converts the Hamiltonian (3), including the edge reconstruction defects, into a set of W independent dual channels, whose Hamiltonian contains hopping integrals t   = t / 2 , t l   = t l / 2 with l = 1 , 2 , , L 1 and on-site energies ±   ε j with j = 1 ,   2 , ,   W , which are eigenvalues obtained from the Hamiltonian of an arbitrary transversal zig-zag line in Figure 2.
The extra carrier density ( n 2 D ) induced by V G in the graphene ribbon can be calculated by [21]
n 2 D =   0 μ D O S ( E )   d E = 0 e V G   e V C D O S ( E )   d E ,
where
D O S ( E ) = 1 π Im T r   [ G + ( E ) ] ,
is the density of states with G + ( E )   = lim η   0 + ( E + i η H ) 1 the retarded Green’s function and η the imaginary part of energy [22]. Combining Equations (2) and (4), V C can be determined as the solution of following self-consistent equation,
e V C   = e 2 d ε r ε 0 0 e V G   e V C D O S ( E )   d E .
On the other hand, the source-drain current ( I s d ) of GFETs at temperature T is given by [23]
I s d ( T ) = V s d g ( μ , T ) ,
where V s d is the source-drain voltage and
g ( μ , T ) = + d E f F D ( E ) E g ( E , 0 ) = g 0 + d E f F D ( E ) E   T ( E )
is the electrical conductance of graphene ribbon at chemical potential μ within the Landauer formalism. In Equation (8), f F D = { exp [ ( E μ ) / k B T ] + 1 } 1 is the Fermi–Dirac distribution, g 0 = 2 e 2 / h is the conductance quantum, T ( E ) is the carrier transmittance of graphene ribbon, and at zero temperature we have [24]
g ( μ , 0 ) = g 0 T ( μ ) .
By means of the independent dual-channel method presented in Appendix A, an armchair-edged graphene ribbon can be transformed into a set of double chains or dual channels (see Figure 2) and then
T ( E ) = j = 1 W T j ( E ) ,
where T j ( E ) is the carrier transmittance along the j-th dual channel. We further developed a new transfer matrix method to calculate T j ( E ) , which is thoroughly presented in Appendix B, and the result can be written as
T j ( E ) = u 0 ( j )   d 0 ( j ) d 0 ( j )   u 0 ( j ) m 11 ( j , E )   d 0 ( j )   u 0 ( j ) + m 12 ( j , E ) | d 0 ( j )   | 2   m 21 ( j , E )   | u 0 ( j )   | 2   m 22 ( j , E )   u 0 ( j )   d 0 ( j ) 2 ,
where m μ ν ( j , E ) are elements of the total transfer matrix, u 0 ( j ) and d 0 ( j ) are respectively initial wavefunction amplitudes at the up and down channels, both from the j-th dual channel (see Appendix B).
One of structural defects commonly observed in graphene ribbons is linear dislocation [25,26], because its presence slightly alters the free energy but may significantly modify the carrier transport. In this article, we consider a distribution of transversal dislocations following the Rudin–Shapiro (RS) sequence, which constitutes one of the closest aperiodic arrays to the random distribution [27]. A Rudin–Shapiro lattice with bond disorder can be built by using four kinds of dislocation hopping integrals t l = t A ,   t B ,   t C   or   t D , respectively denoted by letters A, B, C or D, whose positions along the ribbon are ordered by the substitution rules [28]
A A B = A B , B A C = A C , C D B = D B , D D C = D C
or by the addition rules expressed as
S A ( k + 1 ) = S A ( k ) S B ( k ) ,   S B ( k + 1 ) = S A ( k ) S C ( k ) ,   S C ( k + 1 ) =   S D ( k ) S B ( k ) ,   S D ( k + 1 ) =   S D ( k ) S C ( k ) ,
where the symbol represents the concatenation operation and S α ( k ) is the α-type RS sequence of generation k with α =   A , B , C   or   D . As a consequence, a graphene ribbon based on S α ( k ) has 2 k + 1 transversal zig-zag lines. To have a limited number of dislocations, let us introduce dilute RS lattices, which are built by using periodic segments to construct the RS lattice until generation m 1 and a single transversal dislocation is placed at the initial of next generation m. For example, when k = m = 3 , S A ( 3 ) = A P P P , S B ( 3 ) = B P P P , S C ( 3 ) = C P P P and S D ( 3 ) = D P P P , where P denotes a periodic RS sequence with hopping integral t connecting all transversal zig-zag lines. A detailed discussion of this dilute RS lattice and the real-space renormalization method can be found in Appendix C.
Another quantity that should be considered for determining the electrical conductance is the contact resistance between the graphene ribbon and periodic leads. In the diffusive regime, the total resistance ( R t o t a l ) can be written as [29,30]
R t o t a l = R + 2 R c / W ,
where R = R s h e e t L / W is the electrical resistance of a graphene ribbon with length L and width W , being R s h e e t that of a square graphene sheet, and R c is the specific contact resistance in units of Ω   · μm, which can be experimentally determined by the transfer length method [31]. Hence, the total conductance defined as G t o t a l = R t o t a l 1 can be calculated from Equation (14) through
G t o t a l = G 2 ( R c / W )   G + 1 ,
where G = R 1 is the conductance of graphene ribbon. For a uniform sample, the electrical resistance of a p-micron segment, R( L = p μm), should be p / q times of R( L = q μm) obtained from a q micrometer portion. Therefore, Equation (15) can be rewritten as
G t o t a l ( p   μ m ) = G ( p   μ m ) 2 ( R c / W )   G ( p   μ m ) + 1 = G ( q   μ m ) 2 ( R c / W )   G ( q   μ m ) + p / q ,
where G ( p μ m ) = q G ( q μ m ) / p . This equation will be used in the length variation analysis.

3. Results

The electronic density of states, D O S α ( k , E ) , of a graphene ribbon built by W dual channels based on an α-type dilute RS sequence of generation k can be written as [22]
D O S α ( k , E ) = j = 1 W D O S j α ( k , E ) = j = 1 W 1 π lim η 0 + l = 1 L ( k ) Im   [ G l , l α , j ( k , z ) ] ,
where L ( k ) = 2   ( 2 k + 1 ) is the number of sites in each dual channel and G l , l α , j ( k , z ) is the l-th diagonal element of the retarded Green’s function evaluated at z = E + i η . The D O S j α ( k , E ) of the j-th dual channel can be rewritten as
D O S j α ( k , E ) = 1 π lim η 0 + Im   M α ( k )   G L , L α , j ( k , z )   + P α ( k )   G L _ , L _ α , j ( k , z )   +   Q α ( k )   G R , R α , j ( k , z )   +   R α ( k )   G R _ , R _ α , j ( k , z )   + S α ( k )   G L , L _ α , j ( k , z ) + U α ( k )   G L , R α , j ( k , z ) + V α ( k )   G L , R _ α , j ( k , z ) + W α ( k )   G L _ , R α , j ( k , z ) +   X α ( k )   G L _ , R _ α , j ( k , z ) + Y α ( k )   G R , R _ α , j ( k , z ) + Z α ( k ) ,
where M α ( k )   , , Z α ( k ) are the renormalization coefficients and G μ , ν α , j ( k , z ) with μ , ν = L , L _ , R   or   R _ are matrix elements of the Green’s function corresponding to the left up (L), left down (L), right up (R) and right down (R) sites of the renormalized four-site j-th dual channel, which constitutes the basic element of the real-space renormalization method presented in Appendix C.
Figure 3 shows the electronic density of states (DOS) and the electrical conductance (g) at zero temperature in units of the conductance quantum g 0 = 2 e 2 / h obtained from Equation (9) as functions of the chemical potential (μ) for free-standing graphene ribbons without dislocations, whose widths are W = 0.122     μ m (red lines), 0.073     μ m (dark yellow lines) and 0.037     μ m (blue lines) as well as an arbitrary length of L transversal zig-zag lines. All the ribbons studied in this article are connected to two semi-infinite periodic leads at their ends, as illustrated in Figure 2a, while a surface hopping integral t s = 1.2   t due to the edge reconstruction and an imaginary part of energy η = 10 3 |   t   | are also included.
Observe the Van Hove singularities in the inset of Figure 3a, which are located at the band limits of each dual channel. The DOS spectra of Figure 3a can be analytically verified through [32]
D O S ( E ) = d S / d k | d E / d k   | = j =   1 W 2 | E | π 4 t 2 ε j 2 ( E 2 4 t 2 ε j 2 ) 2 ,
where d S / d k = L a / ( 2 π ) and d E / d k can be obtained from Equation (A9) in Appendix B. Note also the quantized conductance spectra in the inset of Figure 3b and the maximum values of g ( E , 0 ) around E = ±   t in Figure 3b is equal to the number of dual channels (W).
To analyze the electrical conductance (g) of GFETs, we first self-consistently calculate V C from Equation (6) for a given V G and then, the chemical potential μ and g can be respectively determined using Equations (1) and (8). Figure 4 shows the electrical conductance (g) at T   = 0 (blue lines) and at k B T = 0.02 |   t   | (red lines) versus the gate voltage ( V G ) for GFETs of W = 0.122   μ m , i.e., W = 500 dual channels, and L = 4097 zig-zag lines without dislocations, whose capacitances per unit area are (a) C = 10 nF / cm 2 and (b) C = 50 nF / cm 2 . Note the quantized conductance spectra for T   = 0 , in contrast to the smoothed ones when T increases, as well as the growth of conductance with the capacitance.
In Figure 5, the electrical conductance ( g ) at k B T   = 0.026 |   t   | is plotted as a function of both the back gate voltage ( V G ) and the dislocation hopping integral t A for the same graphene ribbon of Figure 4 except by the value of t A maintaining t B = t C = t D = t . This ribbon contains 72 A-type transversal dislocation lines obtained from a A-type dilute RS lattice with k = 12 and dislocations introduced from generation m = 5 , and it is placed on a GFET with a capacitance per unit area C = 10   nF / cm 2 . Observe the V-shaped dependence between g and V G , even for t A = 0.9   t , and this dependence can be noted from Figure 3b, since V G modifies the value of chemical potential μ, as shown in Figure 1c.
Figure 6 shows the electrical conductance (g) at (Figure 6a) k B T = 0.01 | t | and (Figure 6b) k B T = 0.02 | t | with t = 2.70   eV [20] as functions of the gate voltage ( V G ) for GFETs of W = 0.122   μ m (i.e., W = 500 dual channels) and L = 0.217   μ m (red lines), 0.433   μ m (green lines), and 0.866   μ m (blue lines), respectively corresponding to generations k = 10 , 11 and 12 of A-type dilute RS sequences with dislocations introduced from generation m = 3 , hopping integrals t A = 2.389   eV , t B = t C = 2.74   eV and t D = 2.63   eV , where the capacitance per unit area of GFET is C = 10 nF / cm 2 .
The results reveal a more emphasized V-shaped conductance behavior at lower temperatures, while the electrical conductance diminishes with increasing RS generation, being almost inversely proportional to the GFET length, in accordance with the Ohm’s law as well as the length variation measurements carried out in GFET [33].
Finally, let us consider a specific GFET [29] built by an armchair-edged graphene ribbon placed on a substrate of silicon dioxide (SiO2) of thickness d = 285 nm , i.e., a geometric capacitance per unit area of C = ε 0 ε r / d 12.11 nF / cm 2 . We first calculated the electrical conductance g of a graphene ribbon having a width of W   = 9009 dual channels ( W = 2.2   μ m ) and a length of L = 4097 zig-zag lines ( L = 0.866   μ m corresponding to generation 12 of the RS sequence), which contains 1024 transversal dislocation lines from a A-type dilute RS lattice with m = 3 and hopping integrals t = 2.7   eV , t s = 1.2   t , t A = 2.389   eV , t B = t C = 2.74     eV and t D = 2.63     eV as in Figure 6, while the analysis of system length ( L ) effects on the conductance was carried out by means of Equation (16). Additionally, the contact resistance described in ref. [29] is included in the total conductance ( G t o t a l ) using Equation (15). In Figure 7, we present the source-drain current ( I s d ) at T = 300     K as a function of the gate voltage ( V G ) and graphene ribbon length ( L ) for the considered GFET with a source-drain voltage of V s d   = 0.05     volts , where I s d ( T ) (red spheres) was obtained from I s d = V s d G t o t a l and Equation (16) with q = 0.866 μ m and p = L .
The accordance between the theoretical ( I s d T ) and experimental ( I s d E ) source-drain currents can be quantified by means of the standard deviation (σ), given by
σ = 1 N j = 1 N I s d T [ V G ( j ) ] I s d E [ V G ( j ) ] 2 ,
where V G ( j ) = j 51   volts is the j-th analyzed gate voltage and I s d E [ V G ( j ) ] are digitalized experimental data (cyan spheres in Figure 7) from reference [29] with N   = 101 . The resulting σ for graphene ribbons of L = 1 , 2 , 3 , 4 , 5 and 6 μ m are respectively 3.4423, 3.0119, 2.5801, 2.7210, 1.9327 and 1.4795 μA.
Note in Figure 7 the asymmetry of I s d values with respect to V G   = 0 , whose magnitude diminishes as the ribbon length increases. This asymmetry is derived from that of the contact-resistance ( R c ) spectrum in ref. [29], which is probably related to the slight doping of the used graphene ribbon [30]. In fact, the contribution of R c to I s d diminishes as the resistance of graphene ribbon (R) grows [see Equation (15)], while R in turn increases with the ribbon length. Therefore, the asymmetry of I s d respect to V G   = 0 should decrease as the ribbon length grows.

4. Conclusions

The correlation between source-drain current and gate voltage in graphene field-effect transistors (GFETs) is investigated by means of an atomic-scale tight-binding model, where mesoscopic graphene ribbons with edge reconstruction defects and transversal dislocations are addressed by an independent dual-channel transformation, which converts armchair-edged graphene ribbons into a set of independent dual channels.
The electronic transport was studied within the Landauer formalism, and a new transfer matrix method has been further developed for each dual channel with nearest- and next-nearest-neighbor interactions, including aperiodically placed dislocations following a dilute Rudin–Shapiro (RS) sequence, which is one of the closest aperiodic arrangements to the random distribution. The Landauer conductance was confirmed by that obtained from the Kubo–Greenwood formula (see Appendix B), while the density of states from the Green’s function has been verified by analytical solutions of Equation (19) for the case without dislocations.
Quantized electrical conductance spectra are observed even in presence of edge reconstruction defects along mesoscopic graphene ribbons without dislocations, whose maximum value corresponds to the number of dual channels, while the multi-step behavior is softened at finite temperatures. In addition, the back gate voltage ( V G ) creates a shift of the chemical potential in graphene ribbons, which produces a V-shape correlation between the source-drain current ( I s d ) and V G originated from the V-shape DOS spectrum of graphene around the Dirac point. Finally, the realistic calculation of I s d reveals an excellent agreement between the theoretical prediction from a single set of parameters and multiple experimental data obtained from GFETs of six different lengths [29].
The approach of independent dual-channel plus renormalization method presented in this article can be used to study biosensor devices based on GFETs by considering adsorbed molecules as Fano impurities [34]. This study is currently in process.

Author Contributions

Conceptualization, F.S., V.S. and C.W.; methodology, F.S., V.S. and C.W; software, F.S., V.S. and C.W; validation, F.S., V.S. and C.W.; formal analysis, F.S., V.S. and C.W; investigation, F.S., V.S. and C.W; resources, F.S., V.S. and C.W; data curation, F.S., V.S. and C.W; writing—original draft preparation, F.S., V.S. and C.W; writing—review and editing, F.S., V.S. and C.W; visualization, F.S., V.S. and C.W; supervision, F.S., V.S. and C.W; project administration, F.S., V.S. and C.W; funding acquisition, V.S. and C.W All authors have read and agreed to the published version of the manuscript.

Funding

This work has been partially supported by the Consejo Nacional de Ciencia y Tecnología of Mexico (CONACyT) through grant 252943 and by the National Autonomous University of Mexico through projects PAPIIT-IN110020 and PAPIIT-IN112522. The computations were performed at Miztli of DGTIC-UNAM.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Acknowledgments

The technical assistance of Alberto López (R.I.P.), Alejandro Pompa, Cain González, Silvia E. Frausto, and Yolanda Flores are fully appreciated. F.S. acknowledges the postdoctoral fellowship from CONACyT.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Independent Dual-Channel Method

Let us consider an armchair-edged graphene ribbon with a width of four atoms and an arbitrary length of L transversal zig-zag lines described by Hamiltonian (3), which in the matrix form following the atom numbering of Figure A1a can be written as
H = h t 1 0 t 1 h t 2 0 t 2 h 0 t L 1 0 t L 1 h ,
where
h = 0 t s 0 0 t s 0 t 0 0 t 0 t s 0 0 t s 0 ,   t 1 = 0 0 0 0 0 t 1 0 0 0 0 0 0 0 0 0 t 1 ,   t 2 = t 2 0 0 0 0 0 0 0 0 0 t 2 0 0 0 0 0 ,   and   t L 1 = t L 1 0 0 0 0 0 0 0 0 0 t L 1 0 0 0 0 0
being t and t s the interior- and surface-hopping integrals, whose difference is originated from the edge reconstruction of graphene ribbons [20]. In Equation (A2), t l with l = 1 , 2 , , L 1 are hopping integrals that connect zig-zag lines and through them transverse dislocations can be introduced. The eigenvalues of matrix h are ±   ε 1 = ± 1 2 t t 2 + 4 t s 2 and ±   ε 2 = 1 2 t + t 2 + 4 t s 2 .
Let us further introduce a unitary transformation given by
U = u 0 0 0 u 0 0 0 u   with   u = 1 2 υ 1   υ 2 υ 2   υ 1 ε 1 υ 1 / t s ε 2 υ 2 / t s ε 2 υ 2 / t s ε 1 υ 1 / t s ε 1 υ 1 / t s ε 2 υ 2 / t s ε 2 υ 2 / t s ε 1 υ 1 / t s υ 1 υ 2 υ 2 υ 1 ,
where υ j = ( ε j / t s ) 2 + 1 1 / 2 being j = 1   or   2 .
Figure A1. (a) Sketch of an armchair-edged graphene nanoribbon with a width of 4 atoms, arbitrary length L, null on-site energies, hopping integrals t (light blue lines), t 1 (violet lines), t 2 (red lines), t L 1 (green lines) and t s = 1.2   t (thick blue lines); (b) by applying a unitary transformation U, the nanoribbon (a) becomes into two independent dual-channels with on-site energies ( ± ε 1 and ± ε 2 ) and hopping integrals t 1 / 2 (light violet lines), ± t 2 / 2 (light red lines) and ± t L 1 / 2 (light green lines). Numbers inside spheres denote the positions in Hamiltonian matrices (A1) and (A4) respectively before and after the application of U.
Figure A1. (a) Sketch of an armchair-edged graphene nanoribbon with a width of 4 atoms, arbitrary length L, null on-site energies, hopping integrals t (light blue lines), t 1 (violet lines), t 2 (red lines), t L 1 (green lines) and t s = 1.2   t (thick blue lines); (b) by applying a unitary transformation U, the nanoribbon (a) becomes into two independent dual-channels with on-site energies ( ± ε 1 and ± ε 2 ) and hopping integrals t 1 / 2 (light violet lines), ± t 2 / 2 (light red lines) and ± t L 1 / 2 (light green lines). Numbers inside spheres denote the positions in Hamiltonian matrices (A1) and (A4) respectively before and after the application of U.
Nanomaterials 12 03223 g0a1
Applying this transformation (U) to H we obtain
H ˜ = U T H   U = α β 1 0 β 1 α β 2 0 β 2 α 0 β L 1 0 β L 1 α ,
where
α = ε 1 0 0 0 0 ε 2 0 0 0 0 ε 2 0 0 0 0 ε 1 ,   β 1 = t 1 0 0 t 1 0 t 1 t 1 0 0 t 1 t 1 0 t 1 0 0 t 1 ,   β 2 = t 2 0 0 t 2 0 t 2 t 2 0 0 t 2 t 2 0 t 2 0 0 t 2   and   β L 1 = t L 1 0 0 t L 1 0 t L 1 t L 1 0 0 t L 1 t L 1 0 t L 1 0 0 t L 1
being t l = t l / 2 with l = 1 , 2 , , L   1 .
Hence, for a graphene ribbon with a width of four atoms and arbitrary length L, Hamiltonian H can be visualized as two independent dual channels, as illustrated in Figure A1b with on-site energies ± ε j and hopping integrals t 1 / 2 , ±   t 2 / 2 and ± t L 1 / 2 . It is worth mentioning that the independent dual-channel method presented in this appendix includes the ribbon edge reconstruction, and a particular case of this method can be found in reference [35] for t s = t .

Appendix B. Transfer Matrix Method for Dual Channels

Let us consider a dual channel or double chain obtained from the unitary transformation of Appendix A, as illustrated in Figure A2, whose Hamiltonian can be written as
H = ε l u , l u , l d , l d , l + l t l d , l d , l + 1 + u , l u , l + 1 + c . c . +   l t l d , 2 l 1 u , 2 l + u , 2 l 1 d , 2 l d , 2 l u , 2 l + 1 u , 2 l d , 2 l + 1 + c . c . ,
where t l   = t l / 2 is the hopping integral between transversal lines l and l + 1 , u   , l and d , l are respectively the Wannier functions at the up and down sites of transversal line l. In contrast, there is no variation in hopping integrals t l   = t   = t / 2 in periodic leads. Moreover, the system possesses an odd number L of transversal lines, since a Rudin–Shapiro (RS) distribution of bonds is considered.
Figure A2. Dual-channel system of 2L sites with hopping integrals ±   t l connected to two semi-infinite periodic leads with ±   t , where the on-site energies are ε u = ε and ε d = ε . Green dashed lines indicate unit cells whose wavefunction amplitudes without phases are denoted on the spheres.
Figure A2. Dual-channel system of 2L sites with hopping integrals ±   t l connected to two semi-infinite periodic leads with ±   t , where the on-site energies are ε u = ε and ε d = ε . Green dashed lines indicate unit cells whose wavefunction amplitudes without phases are denoted on the spheres.
Nanomaterials 12 03223 g0a2
In these periodic dual-channel leads, the unit cell of length a is constituted by four sites and the electronic wavefunction with wavenumber k can be written as
Ψ ( k ) = m u 0 ( k )   u , 0 + u 1 ( k )   e i k a / 2 u , 1 + d 0 ( k )   d , 0 + d 1 ( k )   e i k a / 2 d , 1 e i k m a ,
where m is the index of unit cells. The stationary Schrödinger equation H Ψ ( k ) = E ( k ) Ψ ( k ) leads to
ε 0 2   t cos   ( k a / 2 ) 2   i   t sin   ( k a / 2 ) 0 ε 2   i   t sin   ( k a / 2 ) 2   t cos   ( k a / 2 ) 2   t cos   ( k a / 2 ) 2   i   t sin   ( k a / 2 ) ε 0 2   i   t sin   ( k a / 2 ) 2   t cos   ( k a / 2 ) 0 ε u 0 ( k ) d 0 ( k ) u 1 ( k ) d 1 ( k ) = E ( k ) u 0 ( k ) d 0 ( k ) u 1 ( k ) d 1 ( k ) ,
whose eigenvalues or dispersion relations are
E ( k ) = { E + ( k ) = + ε 2 + 4 ( t ) 2 4   t ε   cos   ( k a / 2 ) ,  if  ε 2 + 4   ( t ) 2 E ε 2   t E + + ( k ) = + ε 2 + 4 ( t ) 2 + 4   t ε   cos   ( k a / 2 ) ,  if  ε + 2   t E ε 2 + 4   ( t ) 2 E + ( k ) = ε 2 + 4 ( t ) 2 + 4   t ε   cos   ( k a / 2 ) ,  if  ε 2 + 4   ( t ) 2 E ε + 2   t E ( k ) = ε 2 + 4 ( t ) 2 4   t ε   cos   ( k a / 2 ) ,  if    ε 2   t E ε 2 + 4   ( t ) 2
where a is the lattice constant indicated in Figure A2. These four dispersion relations of (A9) are plotted in Figure A3 for t =   0.5 |   t   | , t s =   1.2 |   t   | and ε u = ε d = 1 2 |   t   | + t 2 + 4 t s 2 = 1.8   |   t   | in the first dual channel analyzed in Appendix A.
Figure A3. Dispersion relations E ± + ( k ) and E ± ( k ) versus ka for periodic dual-channel leads with ε u =   ε d   = 1.8 |   t   | and t   =   0.5 |   t   | .
Figure A3. Dispersion relations E ± + ( k ) and E ± ( k ) versus ka for periodic dual-channel leads with ε u =   ε d   = 1.8 |   t   | and t   =   0.5 |   t   | .
Nanomaterials 12 03223 g0a3
Wavefunction amplitudes u l ( k ) and d l ( k ) for l = 0   and   1 obtained from Equation (A8) are summarized in Table A1 for each energy band, shown for example in Figure A3.
Table A1. Amplitudes u l and d l in periodic dual-channel leads with ε u = ε d = ε > 0 and t < 0 .
Table A1. Amplitudes u l and d l in periodic dual-channel leads with ε u = ε d = ε > 0 and t < 0 .
Energy Bandu0 (k)d0 (k)u1 (k)d1 (k)
ε + 2 t E ε 2 + 4 ( t ) 2 i   [   C + ε + E + + ] / S 1 i   [   C + ε + E + + ] / S 1
ε 2 + 4 ( t ) 2 E ε + 2 t i   [ C +   ε + E + ] / S 1 i   [   C + ε + E + ] / S 1
ε 2 + 4 ( t ) 2 E 2 t ε i   [   C ε E + ] / S 1 i   [   C ε E + ] / S 1
  2 t ε E ε 2 + 4 ( t ) 2 i   [   C ε E ] / S 1 i   [ C   ε E ] / S 1
where C = 2 t cos   ( k a / 2 ) and S = 2 t sin   ( k a / 2 ) .
On the other hand, for a dual-channel system of L transversal lines connected to two periodic leads represented by lines zero and L + 1 , the wavefunction of both system and leads can be written as
Ψ = l = 0 L + 1 A l u , l + B l d , l
and then, the stationary Schrödinger equation from Hamiltonian (A6) is given by
{ [ ε u E ]   A 1 + t 0 A 0 + t 1 A 2 t 0 B 0 + t 1 B 2 = 0 [ ε d E ]   B 1 + t 0 B 0 + t 1 B 2 t 0 A 0 + t 1 A 2 = 0 [ ε u E ]   A 2 + t 1 A 1 + t 2 A 3 + t 1 B 1 t 2 B 3 = 0 [ ε d E ]   B 2 + t 1 B 1 + t 2 B 3 + t 1 A 1 t 2 A 3 = 0 [ ε u E ]   A L 1 + t L 2 A L 2 + t L 1 A L + t L 2 B L 2 t L 1 B L = 0 [ ε d E ]   B L 1 + t L 2 B L 2 + t L 1 B L + t L 2 A L 2 t L 1 A L = 0 [ ε u E ]   A L + t L 1 A L 1 + t L A L + 1 t L 1 B L 1 + t L B L + 1 = 0 [ ε d E ]   B L + t L 1 B L 1 + t L B L + 1 t L 1 A L 1 + t L A L + 1 = 0 ,
where ε u = ε , ε d = ε and t 0   = t L   = t . Equations (A11) can be rewritten in matrix form as
A L + 1 B L + 1 = M L ( E ) M 2 j ( E ) M 2 j 1 ( E ) M 0 ( E ) A 0 B 0 = M T ( E ) A 0 B 0 = m 11 ( E ) m 12 ( E ) m 21 ( E ) m 22 ( E ) A 0 B 0 ,
where j = 1 , 2 , , ( L 1 ) / 2 , A 0 ( B 0 ) and A L + 1 ( B L + 1 ) are respectively wavefunction coefficients at the up (down) site of transversal lines zero and L + 1 ,
M 2 l ( E ) = 1 4   ε   t 2 l E 2 ε 2 4   t 2 l 2 ( E + ε ) 2 + 4   t 2 l 2 ( E ε ) 2 4   t 2 l 2 E 2 + ε 2 + 4   t 2 l 2   for   l = 0 , 1 , , L 1 2
M 2 s 1 ( E ) = 1 4   ε   t 2 s 1 E 2 ε 2 4   t 2 s 1 2 ( E ε ) 2 + 4   t 2 s 1 2 ( E ε ) 2 + 4   t 2 s 1 2 E 2 ε 2 4   t 2 s 1 2   for   s = 1 , 2 , , L + 1 2
and the total transfer matrix M T including connections to periodic leads M L and M 0 .
For the case of an electronic plane wave incident from the left lead of Figure A1, the scattered ones by the system are a reflected wave moving to the left and another transmitted one going to the right. Hence, the wavefunction coefficients at transversal line zero and line L + 1 can be respectively written as
A 0 B 0 = u 0 + r   u 0 d 0 + r   d 0   and   A L + 1 B L + 1 = τ   u 0 e i k ( L + 1 ) a τ   d 0 e i k ( L + 1 ) a ,
where r and τ are respectively the reflection and transmission coefficients satisfying |   r   | 2 + |   τ   | 2 = 1 . From Equations (A12) and (A15), we have
u 0 τ e i k ( L + 1 ) a d 0 τ e i k ( L + 1 ) a = M T ( E ) u 0 + u 0 r d 0 + d 0 r ,
which can be rewritten as
τ e i k   L a 0 = u 0 e i k a 0 d 0 e i k a 1 1 m 11 ( E ) m 12 ( E ) m 21 ( E ) m 22 ( E ) u 0 + r   u 0 d 0 + r   d 0 .
Therefore, the transmittance (T) of electronic waves in a dual channel is given by
T ( E ) = | τ | 2 = u 0   d 0 d 0   u 0 m 11 ( E )   d 0 u 0 + m 12 ( E )   | d 0   | 2   m 21 ( E )   | u 0   | 2   m 22 ( E )   u 0   d 0 2
where u 0 and d 0 are given in Table A1.
The electrical conductance (g) at 0 K obtained from Equations (9) and (A18) as a function of the chemical potential (μ) is shown in Figure A4 for the graphene ribbon of Figure A1 containing L = 9 transversal lines (see Table A3 for generation k = 3 ) (a) without and (b) with a single transversal dislocation t A = 1.1   t , while t B =     t . These Landauer conductance spectra are further compared with those obtained from the Kubo formalism within the linear response theory. In general, the electric conductivity (σ) can be calculated by means of the Kubo–Greenwood formula [22],
Figure A4. Electrical conductance (g) at 0 K in units of g 0 = 2 e 2 / h obtained from the Landauer (red lines) and Kubo–Greenwood formulas using η = 10 3 |   t   | (blue lines) and η = 10 4 |   t   | (green lines) as functions of the chemical potential (μ) for graphene ribbons with W = 4 , L   = 9 , t s = 1.2   t ,   t B   = t , (a,a′) t A = t and (b,b’) t A = 1.1   t , whose magnifications are shown in (a’) and (b’).
Figure A4. Electrical conductance (g) at 0 K in units of g 0 = 2 e 2 / h obtained from the Landauer (red lines) and Kubo–Greenwood formulas using η = 10 3 |   t   | (blue lines) and η = 10 4 |   t   | (green lines) as functions of the chemical potential (μ) for graphene ribbons with W = 4 , L   = 9 , t s = 1.2   t ,   t B   = t , (a,a′) t A = t and (b,b’) t A = 1.1   t , whose magnifications are shown in (a’) and (b’).
Nanomaterials 12 03223 g0a4
σ ( μ , ω , T ) = 2   e 2 π   m 2 Ω d E f F D ( E ) f F D ( E + ω ) ω T r p ^ x Im G + ( z + ω ) p ^ x Im G + ( z ) ,
where e and m are respectively the electrical charge and the mass of electron, Ω is the area of graphene ribbon, and f F D ( E ) =   1 + exp ( E μ ) / k B T   1 is the Fermi–Dirac distribution with the chemical potential μ, the Boltzmann constant k B and temperature T. p ^ x   = i m   [ H ^ , x ] / is the projection of the momentum operator along the oscillating electrical field with frequency ω, and G + ( z ) with z = E + i η is the retarded one-particle Green’s function. For the case of ω = T = 0 , Equation (A19) can be rewritten as [36].
g ( μ , 0 , 0 ) = σ ( μ , 0 , 0 ) W L = 2   e 2 π   m 2 L 2 T r p ^ x Im G + ( μ ) p ^ x Im G + ( μ ) .
In Figure A4 the Kubo conductance spectra are obtained from Equation (A20) for two imaginary parts of energy η = 10 3 |   t   | (blue lines) and η = 10 4 |   t   | (green lines), where the retarded Green’s function was calculated through G ( z ) = ( z H ) 1 from the Dyson equation for a Hamiltonian of 36 atoms since W   = 4 and L = 9 for a dilute A-type RS lattice with k = 4 and m = 3 containing only dislocation hopping integrals t A and t B , whose boundary conditions were obtained from periodic leads of 200,000 transversal lines. Observe in Figure A4a’,b’ the coincidence between electrical conductance spectra from Landauer and Kubo formalisms when η 0 in the latter one.

Appendix C. Real-Space Renormalization Method for Dilute Rudin–Shapiro Dual Channels

The Rudin–Shapiro (RS) sequence of generation k, S ( k ) , can be built using the substitution rules (12) or addition rules (13). The latter can be summarized into a single equation given by
S α ( k ) = S β ( k 1 ) S γ ( k 1 ) ,
whose subindices α, β and γ denote A, B, C or D-type RS sequence according to Table A2.
Table A2. Addition rules for all type Rudin–Shapiro sequences.
Table A2. Addition rules for all type Rudin–Shapiro sequences.
α-Typeβ-Typeγ-Type
AAB
BAC
CDB
DDC
Let us consider an RS dual channel with bond disorder formed by alternating hopping integrals t A ,   t B ,   t C and t D , instead of letters A, B, C and D in the RS sequences, being t α = t α / 2 with α = A , B , C   or   D . Moreover, a dilute RS dual channel can be built using periodic RS segments of generation m 1 formed by hopping integrals t = t / 2 and placing a single hopping integral defect t α at the beginning of S α ( m ) , while subsequent generations of dilute RS dual channels are constructed following the standard addition rule (A21) and Table A2. For example, dilute RS dual channels with m = 2 are presented in Table A3.
Table A3. Dilute Rudin–Shapiro dual channels with m = 2.
Table A3. Dilute Rudin–Shapiro dual channels with m = 2.
GenerationA-TypeB-TypeC-TypeD-Type
1Nanomaterials 12 03223 i001Nanomaterials 12 03223 i001Nanomaterials 12 03223 i001Nanomaterials 12 03223 i001
2Nanomaterials 12 03223 i002Nanomaterials 12 03223 i003Nanomaterials 12 03223 i004Nanomaterials 12 03223 i005
3Nanomaterials 12 03223 i006Nanomaterials 12 03223 i007Nanomaterials 12 03223 i008Nanomaterials 12 03223 i009
Note: Symbol meaning Nanomaterials 12 03223 i010.
In this appendix, a new real-space renormalization method is developed for dilute RS dual channels, as an extension of that for single channels developed in reference [37], where the freedom degrees of interior sites are removed but their exact participation is maintained in the results. As an example, we present in Figure A5 this renormalization procedure for A-type RS dual channels by starting from generation one, which consists of six sites (blue spheres) with on-site energies ε u   = ε and ε d   =   ε , respectively, for up and down sites, as well as hopping integrals t = t / 2 (violet lines) and   t (red lines). The central sites located at the second transversal line can be renormalized by reducing the six-site dual channel into a four-site one with effective on-site energies ε L A ( 1 ) , ε L _ A ( 1 ) , ε R A ( 1 ) and ε R _ A ( 1 ) , respectively, for left-up, left-down, right-up and right-down sites, as well as effective hopping integrals t h ( 1 ) ,   t h _ ( 1 ) ,   t v ( 1 ) ,   t v _ ( 1 ) ,   t d ( 1 ) and t d _ ( 1 ) , correspondingly denoting up horizontal, down horizontal, left vertical, right vertical, decreasing diagonal and increasing diagonal bonds. Using the addition rule for A-type RS dual channels given by Equation (A21) and Table A2, this four sites dual channel is concatenated to another four-site renormalized B-type RS dual channel to build a six-site dual channel of generation 2, where the on-site energies of central sites are ε C A ( 2 ) = ε R A ( 1 ) + ε L B ( 1 ) ε and ε C _ A ( 2 ) = ε R _ A ( 1 ) + ε L _ B ( 1 ) + ε respectively for up and down ones connected by the hopping integral t C A ( 2 ) = t v _ A ( 1 ) + t v B ( 1 ) , represented by a double line.
Figure A5. Schematic representation of the renormalization procedure to build an A-type RS dual channel of generation K connected to two periodic leads of generation K’ (gray spheres), where the effective on-site energies ε l α ( k , z ) and hopping integrals t j α ( k , z ) are respectively denoted as ε l α ( k ) and t j α ( k ) for α = A   or   B , l = L , L _ , C , C _ , R   or   R _ and j = h , h _ , v , v _ , d , d _ , C   or   C _ .
Figure A5. Schematic representation of the renormalization procedure to build an A-type RS dual channel of generation K connected to two periodic leads of generation K’ (gray spheres), where the effective on-site energies ε l α ( k , z ) and hopping integrals t j α ( k , z ) are respectively denoted as ε l α ( k ) and t j α ( k ) for α = A   or   B , l = L , L _ , C , C _ , R   or   R _ and j = h , h _ , v , v _ , d , d _ , C   or   C _ .
Nanomaterials 12 03223 g0a5
In general, for an α-type RS dual channel of generation k, similar addition and renormalization procedures were carried out. The effective on-site energies ε l α ( k , E ) with l = L , L _ , R   or   R _ and hopping integrals t j α ( k , E ) , being l = h , h _ , v , v _ , d   or   d _ , of four-site dual channel can be obtained from the stationary Schrödinger equation:
( H E I ) Ψ = 0
where I is the identity matrix, H is the Hamiltonian of a concatenated six-site dual channel built from two dual channels of generations k 1 and Ψ = j b j j is the electronic wavefunction with j = L , L _ , C , C _ , R , or   R _ . Equation (A22) can be rewritten as
ε L β ( k 1 , E ) E t v β ( k 1 , E ) t h β ( k 1 , E ) t d β ( k 1 , E ) 0 0 t v β ( k 1 , E ) ε L ̲ β ( k 1 , E ) E t d ̲ β ( k 1 , E ) t h ̲ β ( k 1 , E ) 0 0 t h β ( k 1 , E ) t d ̲ β ( k 1 , E ) ε C α ( k , E ) E t c α ( k , E ) t h ̲ γ ( k 1 , E ) t d γ ( k 1 , E ) t d β ( k 1 , E ) t h ̲ β ( k 1 , E ) t c α ( k , E ) ε C ̲ α ( k , E ) E t d ̲ γ ( k 1 , E ) t h ̲ γ ( k 1 , E ) 0 0 t h γ ( k 1 , E ) t d ̲ γ ( k 1 , E ) ε R γ ( k 1 , E ) E t v ̲ γ ( k 1 , E ) 0 0 t d γ ( k 1 , E ) t h ̲ γ ( k 1 , E ) t v ̲ γ ( k 1 , E ) ε R ̲ γ ( k 1 , E ) E b L β b L ̲ β b C α b C ̲ α b R γ b R ̲ γ = ( 0 0 0 0 0 0 ) ,
whose two central equations can be considered as a set of equations for b C α and b C _ α , with solutions given by
b C α = Γ 1 α ( k , E )   b L β + Γ 2 α ( k , E )   b L _ β + Γ 3 α ( k , E )   b R γ + Γ 4 α ( k , E )   b R _ γ b C _ α = Θ 1 α ( k , E )   b L β + Θ 2 α ( k , E )   b L _ β + Θ 3 α ( k , E )   b R γ + Θ 4 α ( k , E )   b R _ γ ,
where
Γ 1 α ( k , E ) = θ 1 α ( k , E )   [ t h β ( k   1 , E )   +   ϕ 2 α ( k , E )   t d β ( k   1 , E ) ] Γ 2 α ( k , E ) = θ 1 α ( k , E )   [ t d _ β ( k   1 , E )   +   ϕ 2 α ( k , E )   t h _ β ( k   1 , E ) ] Γ 3 α ( k , E ) = θ 1 α ( k , E )   [ t h γ ( k 1 , E ) + ϕ 2 α ( k , E )   t d _ γ ( k 1 , E ) ] Γ 4 α ( k , E ) = θ 1 α ( k , E )   [ t d γ ( k 1 , E ) + ϕ 2 α ( k , E )   t h _ γ ( k 1 , E ) ]
and
Θ 1 α ( k , E ) = θ 2 α ( k , E )   [ t d β ( k 1 , E ) + ϕ 1 α ( k , E )   t h β ( k 1 , E ) ] Θ 2 α ( k , E ) = θ 2 α ( k , E )   [ t h _ β ( k 1 , E ) + ϕ 1 α ( k , E )   t d _ β ( k 1 , E ) ] Θ 3 α ( k , E ) = θ 2 α ( k , E )   [ t d _ γ ( k 1 , E ) + ϕ 1 α ( k , E )   t h γ ( k 1 , E ) ] Θ 4 α ( k , E ) = θ 2 α ( k , E )   [ t h _ γ ( k 1 , E ) + ϕ 1 α ( k , E )   t d γ ( k 1 , E ) ] ,
where
θ 1 α ( k , E ) = E ε c α ( k , E ) [ t c α ( k , E ) ] 2 E ε c _ α ( k , E ) 1 θ 2 α ( k , E ) = E ε c _ α ( k , E ) [ t c α ( k , E ) ] 2 E ε c α ( k , E ) 1   and   ϕ 1 α ( k , E ) = t c α ( k , E ) E ε c α ( k , E ) ϕ 2 α ( k , E ) = t c α ( k , E ) E ε c _ α ( k , E )
with
ε c α ( k , E ) = ε R β ( k 1 , E ) + ε L γ ( k 1 , E ) ε ε c _ α ( k , E ) = ε R _ β ( k 1 , E ) + ε L _ γ ( k 1 , E ) + ε t c α ( k , E ) = t v _ β ( k 1 , E ) + t v γ ( k 1 , E ) .
Substituting Equation (A24) into (A23), we obtain the stationary Schrödinger equation of a four-site dual channel given by
ε L α ( k , E ) E t v α ( k , E ) t h α ( k , E ) t d α ( k , E ) t v α ( k , E ) ε L _ α ( k , E ) E t d _ α ( k , E ) t h _ α ( k , E ) t h α ( k , E ) t d _ α ( k , E ) ε R α ( k , E ) E t v _ α ( k , E ) t d α ( k , E ) t h _ α ( k , E ) t v _ α ( k , E ) ε R _ α ( k , E ) E b L α b L _ α b R α b R _ α = 0 0 0 0 ,
where
{ ε L α ( k , E ) = ε L β ( k   1 , E ) + Γ 1 α ( k , E )   t h β ( k   1 , E ) + Θ 1 α ( k , E )   t d β ( k 1 , E ) ε L _ α ( k , E ) = ε L _ β ( k   1 , E ) + Θ 2 α ( k , E )   t h _ β ( k   1 , E ) + Γ 2 α ( k , E )   t d _ β ( k 1 , E ) ε R α ( k , E ) = ε R γ ( k   1 , E ) + Γ 3 α ( k , E )   t h γ ( k   1 , E ) + Θ 3 α ( k , E )   t d _ γ ( k 1 , E ) ε R _ α ( k , E ) = ε R _ γ ( k   1 , E ) + Θ 4 α ( k , E )   t h _ γ ( k   1 , E ) + Γ 4 α ( k , E )   t d γ ( k 1 , E )
and
{ t h α ( k , E ) = Γ 3 α ( k , E )   t h β ( k 1 , E ) + Θ 3 α ( k , E )   t d β ( k 1 , E ) t h _ α ( k , E ) = Γ 4 α ( k , E )   t d _ β ( k 1 , E ) + Θ 4 α ( k , E )   t h _ β ( k 1 , E ) t d α ( k , E ) = Γ 4 α ( k , E )   t h β ( k 1 , E ) + Θ 4 α ( k , E )   t d β ( k 1 , E ) t d _ α ( k , E ) = Γ 3 α ( k , E )   t d _ β ( k 1 , E ) + Θ 3 α ( k , E )   t h _ β ( k 1 , E ) t v α ( k , E ) = t v β ( k 1 , E ) + Γ 2 α ( k , E )   t h β ( k 1 , E ) + Θ 2 α ( k , E )   t d β ( k 1 , E ) t v _ α ( k , E ) = t v _ γ ( k 1 , E ) + Γ 4 α ( k , E )   t h γ ( k 1 , E ) + Θ 4 α ( k , E )   t d _ γ ( k 1 , E ) .
These effective on-site energies and hopping integrals are those indicated in Figure A5.
On the other hand, the electronic density of states (DOS) for an α-type RS dual channel of generation k evaluated at energy E can be written as [22]
D O S α ( k , E ) = 1 π lim η 0 + j = 1 L ( k ) Im   [ G j , j α ( z ) ] ,
where L ( k ) = 2   ( 1 + 2 k ) is the number of sites in the dual channel and G j , j α ( z ) is the j-th diagonal element of retarded Green’s function evaluated at z = E + i η . Extending the real-space renormalization method developed for Fibonacci chains [37,38], the D O S α ( k , E ) of (A32) can be rewritten as
D O S α ( k , E ) = 1 π lim η 0 + Im   [ M α ( k , z )   G L , L α ( k , z )   + P α ( k , z )   G L _ , L _ α ( k , z )   +   Q α ( k , z )   G R , R a ( k , z )       + R α ( k , z )   G R _ , R _ α ( k , z ) + S α ( k , z )   G L , L _ α ( k , z ) + U α ( k , z )   G L , R α ( k , z ) + V α ( k , z )   G L , R _ α ( k , z )   + W α ( k , z )   G L _ , R α ( k , z ) + X α ( k , z )   G L _ , R _ α ( k , z ) + Y α ( k , z )   G R , R _ α ( k , z ) + Z α ( k , z ) ] ,
where M α ( k , z ) ,   , Z α ( k , z ) are the renormalization coefficients and G l , j α ( k , z ) with l , j = L , L _ , R   or   R _ are elements of the Green’s function matrix, with L ( L _ ) and R ( R _ ) respectively being the up (down) left and right sites of the renormalized four-site dual channel, as illustrated in Figure A5. From Equation (A32) and the addition rule (A21), we have
D O S α ( k , E ) = D O S β ( k 1 , E ) + D O S γ ( k 1 , E ) G C , C α ( k , z ) G C _ , C _ α ( k , z ) ,
where subindex C denotes the central sites at the union of two renormalized dual channels. Equation (A34) can be rewritten using Equation (A33) as
D O S α ( k , E ) = 1 π lim η 0 + Im   [ M β ( k 1 , z ) G L , L β ( k , z ) + P β ( k 1 , z ) G L _ , L _ β ( k , z ) + Q β ( k 1 , z ) G C , C β ( k , z ) + R β ( k 1 , z ) G C _ , C _ β ( k , z ) + S β ( k 1 , z ) G L , L _ β ( k , z ) + U β ( k 1 , z ) G L , C β ( k , z ) + V β ( k 1 , z ) G L , C _ β ( k , z ) + W β ( k 1 , z ) G L _ , C β ( k , z ) + X β ( k 1 , z ) G L _ , C _ β ( k , z ) + Y β ( k 1 , z ) G C , C _ β ( k , z ) + Z β ( k 1 , z ) + M γ ( k 1 , z ) G C , C γ ( k , z ) + P γ ( k 1 , z ) G C _ , C _ γ ( k , z ) + Q γ ( k 1 , z ) G R , R γ ( k , z ) + R γ ( k 1 , z ) G R _ , R _ γ ( k , z ) + S γ ( k 1 , z ) G C , C _ γ ( k , z ) + U γ ( k 1 , z ) G C , R γ ( k , z ) + V γ ( k 1 , z ) G C , R _ γ ( k , z ) + W γ ( k 1 , z ) G C _ , R γ ( k , z ) + X γ ( k 1 , z ) G C _ , R _ γ ( k , z ) + Y γ ( k 1 , z ) G R , R _ γ ( k , z ) + Z γ ( k 1 , z ) G C , C γ ( k , z ) G C _ , C _ γ ( k , z ) ] .
The Dyson equation of a six-site α-type RS dual channel (see Figure A5 for α = A ) given by ( z I H ) G = I can be rewritten as
{ [ z ε L β ( k 1 ) ] G L , j β ( k 1 , z ) = δ L , j + t v β ( k 1 ) G L _ , j β ( k 1 , z ) + t h β ( k 1 ) G C , j α ( k , z ) + t d β ( k 1 ) G C _ , j α ( k , z ) [ z ε L _ β ( k 1 ) ] G L _ , j β ( k 1 , z ) = δ L _ , j + t v β ( k 1 ) G L , j β ( k 1 , z ) + t d _ β ( k 1 ) G C , j α ( k , z ) + t h _ β ( k 1 ) G C _ , j α ( k , z )   [ z ε C α ( k ) ] G C , j α ( k , z ) = δ C , j + t h β ( k 1 ) G L , j β ( k 1 , z ) + t d _ β ( k 1 ) G L _ , j β ( k 1 , z ) + t h _ γ ( k 1 ) G R , j γ ( k 1 , z ) + t d γ ( k 1 ) G R _ , j γ ( k 1 , z ) + t c α ( k ) G C _ , j α ( k , z ) [ z ε C _ α ( k ) ] G C _ , j α ( k , z ) = δ C _ , j + t d β ( k 1 ) G L , j β ( k 1 , z ) + t h _ β ( k 1 ) G L _ , j β ( k 1 , z ) + t d _ γ ( k 1 ) G R , j γ ( k 1 , z ) + t h _ γ ( k 1 ) G R _ , j γ ( k 1 , z ) + t c α ( k ) G C , j α ( k , z ) [ z ε R γ ( k 1 ) ] G R , j β ( k 1 , z ) = δ R , j + t v _ γ ( k 1 ) G R _ , j β ( k 1 , z ) + t h γ ( k 1 ) G C , j α ( k , z ) + t d _ γ ( k 1 ) G C _ , j α ( k , z ) [ z ε R _ γ ( k 1 ) ] G R _ , j β ( k 1 , z ) = δ R _ , j + t v _ γ ( k 1 ) G R , j β ( k 1 , z ) + t d γ ( k 1 ) G C , j α ( k , z ) + t h _ γ ( k 1 ) G C _ , j α ( k , z )
where j = L , L _ , C , C _ , R , or   R _ , and then Equation (A36) represents the 36 equations of six effective sites. The two central equations of (A36) can be visualized as a system of equations for G C , j α ( k , z ) and G C _ , j α ( k , z ) , whose solutions are
{ G C , j α ( k , z ) = Γ 1 α ( k , z )   G L , j β ( k 1 , z ) + Γ 2 α ( k , z )   G L _ , j β ( k 1 , z ) + Γ 3 α ( k , z )   G R , j γ ( k 1 , z ) + Γ 4 α ( k , z )   G R _ , j γ ( k 1 , z ) + δ C , j + θ 1 α ( k , z ) δ C _ , j G C _ , j α ( k , z ) = Θ 1 α ( k , z )   G L , j β ( k 1 , z ) + Θ 2 α ( k , z )   G L _ , j β ( k 1 , z ) + Θ 3 α ( k , z )   G R , j γ ( k 1 , z ) + Θ 4 α ( k , z )   G R _ , j γ ( k 1 , z ) + δ C _ , j + θ 2 α ( k , z ) δ C , j .
Substituting Equation (A37) into (A35) and comparing it to Equation (A33), we obtain the following iteration relations for renormalization coefficients given by
M α ( k , z )   = M β ( k 1 , z ) + Γ 1 α ( k , z )   [ Ξ 1 α ( k , z ) + U β ( k 1 , z ) ] + Θ 1 α ( k , z )   [ V β ( k 1 , z ) + Ω 1 α ( k , z ) ]
P α ( k , z )   =   P β ( k 1 , z ) + Γ 2 α ( k , z )   [ Ξ 2 α ( k , z ) + W β ( k 1 , z ) ] + Θ 2 α ( k , z )   [ X β ( k 1 , z ) + Ω 2 α ( k , z ) ] ,
Q α ( k , z )   = Q γ ( k 1 , z ) + Γ 3 α ( k , z )   [ Ξ 3 α ( k , z ) + U γ ( k 1 , z ) ] + Θ 3 α ( k , z )   [ W γ ( k 1 , z ) + Ω 3 α ( k , z ) ] ,
R α ( k , z ) =   R γ ( k 1 , z ) + Γ 4 α ( k , z )   [ Ξ 4 α ( k , z ) + V γ ( k 1 , z ) ] + Θ 4 α ( k , z )   [ X γ ( k 1 , z ) + Ω 4 α ( k , z ) ] ,
S α ( k , z ) = S β ( k 1 , z ) + Γ 2 α ( k , z )   [ Ξ 1 α ( k , z ) + U β ( k 1 , z ) ] + Γ 1 α ( k , z )   [ Ξ 2 α ( k , z ) + W β ( k 1 , z ) ]   + Θ 2 α ( k , z )   [ V β ( k 1 , z ) + Ω 1 α ( k , z ) ] + Θ 1 α ( k , z )   [ X β ( k 1 , z ) + Ω 2 α ( k , z ) ] ,
U α ( k , z ) = Γ 3 α ( k , z )   [ Ξ 1 α ( k , z ) + U β ( k 1 , z ) ] + Γ 1 α ( k , z )   [ Ξ 3 α ( k , z ) + U γ ( k 1 , z ) ]   + Θ 3 α ( k , z )   [ V β ( k 1 , z ) + Ω 1 α ( k , z ) ] + Θ 1 α ( k , z )   [ W γ ( k 1 , z ) + Ω 3 α ( k , z ) ] ,
V α ( k , z ) = Γ 4 α ( k , z )   [ Ξ 1 α ( k , z ) + U β ( k 1 , z ) ] + Γ 1 α ( k , z )   [ Ξ 4 α ( k , z ) + V γ ( k 1 , z ) ]   + Θ 4 α ( k , z )   [ V β ( k 1 , z ) + Ω 1 α ( k , z ) ] + Θ 1 α ( k , z )   [ X γ ( k 1 , z ) + Ω 4 α ( k , z ) ] ,
W α ( k , z ) = Γ 3 α ( k , z )   [ Ξ 2 α ( k , z ) + W β ( k 1 , z ) ] + Γ 2 α ( k , z )   [ Ξ 3 α ( k , z ) + U γ ( k 1 , z ) ]   + Θ 3 α ( k , z )   [ X β ( k 1 , z ) + Ω 2 α ( k , z ) ] + Θ 2 α ( k , z )   [ W γ ( k 1 , z ) + Ω 3 α ( k , z ) ] ,
X α ( k , z ) = Γ 4 α ( k , z )   [ Ξ 2 α ( k , z ) + W β ( k 1 , z ) ] + Γ 2 α ( k , z )   [ Ξ 4 α ( k , z ) + V γ ( k 1 , z ) ]   + Θ 4 α ( k , z )   [ X β ( k 1 , z ) + Ω 2 α ( k , z ) ] + Θ 2 α ( k , z )   [ X γ ( k 1 , z ) + Ω 4 α ( k , z ) ] ,
Y α ( k , z ) =   Y γ ( k 1 , z ) + Γ 4 α ( k , z )   [ Ξ 3 α ( k , z ) + U γ ( k 1 , z ) ] + Γ 3 α ( k , z )   [ Ξ 4 α ( k , z ) + V γ ( k 1 , z ) ]   + Θ 4 α ( k , z )   [ W γ ( k 1 , z ) + Ω 3 α ( k , z ) ] + Θ 3 α ( k , z )   [ X γ ( k 1 , z ) + Ω 4 α ( k , z ) ] ,
Z α ( k , z ) = Z β ( k 1 , z ) + Z γ ( k 1 , z ) + θ 1 α ( k , z )   [ Q β ( k 1 , z ) + M γ ( k 1 , z ) 1 ]   + θ 2 α ( k , z )   [ R β ( k 1 , z ) + P γ ( k 1 , z ) 1 ] + ϕ 2 α ( k , z )   [ Y β ( k 1 , z ) + S γ ( k 1 , z ) ] ,
with
Ξ ρ α ( k , z ) = Γ ρ α ( k , z )   [ Q β ( k 1 , z ) + M γ ( k 1 , z ) 1 ] + Θ ρ α ( k , z )   [ Y β ( k 1 , z ) + S γ ( k 1 , z ) ]
and
Ω ρ α ( k , z )   =   Θ ρ α ( k , z )   [ R β ( k 1 , z ) + P γ ( k 1 , z ) 1 ] ,
where ρ = 1 , 2 , 3   or   4 . The initial conditions for this renormalization method are
M α ( 1 , z ) = P α ( 1 , z ) = Q α ( 1 , z ) = R α ( 1 , z )   = 1 + ( ξ 2 + ξ + 2 )   t α 2 U α ( 1 , z )   = V α ( 1 , z ) = W α ( 1 , z ) = X α ( 1 , z ) = 2 ( ξ 2 ξ + 2 )   t α 2 S α ( 1 , z )   = Y α ( 1 , z ) = 2 ( ξ 2 + ξ + 2 )   t α 2 Z α ( 1 , z ) = ( ξ + ξ + )
and
t h α ( 1 , z ) = t h _ α ( 1 , z ) = t d α ( 1 , z ) = t d _ α ( 1 , z ) = ( ξ ξ + )   t α 2 t v α ( 1 , z ) = t v _ α ( 1 , z ) = ( ξ + ξ + )   t α 2 ε L α ( 1 , z ) = ε R α ( 1 , z ) = ( ξ + ξ + )   t α 2 ε ε L _ α ( 1 , z ) = ε R _ α ( 1 , z ) = ( ξ + ξ + )   t α 2 + ε ,
where ξ ±   = ( z ± ε ) 1 . The elements of Green’s function in Equation (A33) are numerically calculated from the Dyson equation, written as an 8 × 8 matrix for the dual channel of the K-th generation connected to two periodic leads at its ends, as illustrated at the end of Figure A5.

References

  1. 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. [Google Scholar] [CrossRef]
  2. Bolotin, K.I.; Sikes, K.J.; Jiang, Z.; Klima, M.; Fudenberg, G.; Hone, J.; Kim, P.; Stormer, H.L. Ultrahigh electron mobility in suspended graphene. Solid State Commun. 2008, 146, 351–355. [Google Scholar] [CrossRef]
  3. Balandin, A.A.; Ghosh, S.; Bao, W.; Calizo, I.; Teweldebrhan, D.; Miao, F.; Lau, C.N. Superior thermal conductivity of single-layer graphene. Nano Lett. 2008, 8, 902–907. [Google Scholar] [CrossRef]
  4. Lee, C.; Wei, X.; Kysar, J.W.; Hone, J. Measurement of the elastic properties and intrinsic strength of monolayer graphene. Science 2008, 321, 385–388. [Google Scholar] [CrossRef]
  5. Dhinakaran, V.; Lavanya, M.; Vigneswari, K.; Ravichandran, M.; Vijayakumar, M.D. Review on exploration of graphene in diverse applications and its future horizon. Mater. Today Proc. 2020, 27, 824–828. [Google Scholar] [CrossRef]
  6. Schedin, F.; Geim, A.K.; Morozov, S.V.; Hill, E.W.; Blake, P.; Katsnelson, M.I.; Novoselov, K.S. Detection of individual gas molecules adsorbed on graphene. Nat. Mater. 2007, 6, 652–655. [Google Scholar] [CrossRef] [PubMed]
  7. Béraud, A.; Sauvage, M.; Bazán, C.M.; Tie, M.; Bencherifa, A.; Bouilly, D. Graphene field-effect transistors as bioanalytical sensors: Design, operation and performance. Analyst 2021, 146, 403–428. [Google Scholar] [CrossRef] [PubMed]
  8. Wang, S.; Qi, X.; Hao, D.; Moro, R.; Ma, Y.; Ma, L. Recent advances in graphene-based field-effect-transistor biosensors: A review on biosensor designing strategy. J. Electrochem. Soc. 2022, 169, 027509. [Google Scholar] [CrossRef]
  9. Wang, M.; Yang, E.-H. THz applications of 2D materials: Graphene and beyond. Nano-Struct. Nano-Objects 2018, 15, 107–113. [Google Scholar] [CrossRef]
  10. Liu, J.; Li, X.; Jiang, R.; Yang, K.; Zhao, J.; Khan, S.A.; He, J.; Liu, P.; Zhu, J.; Zeng, B. Recent progress in the development of graphene detector for terahertz detection. Sensors 2021, 21, 4987. [Google Scholar] [CrossRef]
  11. Lone, S.; Bhardwaj, A.; Pandit, A.K.; Gupta, S.; Mahajan, S. A review of graphene nanoribbon field-effect transistor structures. J. Electron. Mater. 2021, 50, 3169–3186. [Google Scholar] [CrossRef]
  12. Freitag, M.; Steiner, M.; Martin, Y.; Perebeinos, V.; Chen, Z.; Tsang, J.C.; Avouris, P. Energy dissipation in graphene field-effect transistors. Nano Lett. 2009, 9, 1883–1888. [Google Scholar] [CrossRef]
  13. Yu, G.L.; Jalil, R.; Belle, B.; Mayorov, A.S.; Blake, P.; Schedin, F.; Morozov, S.V.; Ponomarenko, L.A.; Chiappini, F.; Wiedmann, S.; et al. Interaction phenomena in graphene seen through quantum capacitance. Proc. Natl. Acad. Sci. USA 2013, 110, 3282–3286. [Google Scholar] [CrossRef] [PubMed]
  14. Wen, R.; Jiang, Z.; Miao, R.; Wang, L.; Liang, Y.; Deng, J.; Shao, Q.; Zhang, J. Electronic transport properties of B/N/P co-doped armchair graphene nanoribbon field effect transistor. Diam. Relat. Mater. 2022, 124, 108893. [Google Scholar] [CrossRef]
  15. Radsar, T.; Khalesi, H.; Ghods, V.; Izadbakhsh, A. Effects of channel dimension and doping concentration of source and drain contacts on GNRFET performance. Silicon 2021, 13, 3337–3350. [Google Scholar] [CrossRef]
  16. Saltzgaber, G.; Wojcik, P.; Sharf, T.; Leyden, M.R.; Wardini, J.L.; Heist, C.A.; Adenuga, A.A.; Remcho, V.T.; Minot, E.D. Scalable graphene field-effect sensors for specific protein detection. Nanotechnology 2013, 24, 355502. [Google Scholar] [CrossRef]
  17. Sánchez, V.; Wang, C. Real space theory for electron and phonon transport in aperiodic lattices via renormalization. Symmetry 2020, 12, 430. [Google Scholar] [CrossRef]
  18. Shylau, A.A.; Kłos, J.W.; Zozoulenko, I.V. Capacitance of graphene nanoribbons. Phys. Rev. B 2009, 80, 205402. [Google Scholar] [CrossRef]
  19. Griffiths, D.J. Introduction to Electrodynamics, 4th ed.; Cambridge University Press: Cambridge, UK, 2017; pp. 106–190. [Google Scholar]
  20. Son, Y.-W.; Cohen, M.L.; Louie, S.G. Energy gaps in graphene nanoribbons. Phys. Rev. Lett. 2006, 97, 216803. [Google Scholar] [CrossRef]
  21. Fernández-Rossier, J.; Palacios, J.J.; Brey, L. Electronic structure of gated graphene and graphene ribbons. Phys. Rev. B 2007, 75, 205441. [Google Scholar] [CrossRef] [Green Version]
  22. Economou, E.N. Green’s Functions in Quantum Physics, 3rd ed.; Springer: Berlin/Heidelberg, Germany, 2006; pp. 14–16, 184. [Google Scholar]
  23. Datta, S. Lessons from Nanoelectronics. A New Perspective on Transport—Part A: Basic Concepts, 2nd ed.; World Scientific: Singapore, 2017; pp. 32, 123. [Google Scholar]
  24. Landauer, R. Electrical resistance of disordered one-dimensional lattices. Philos. Mag. 1970, 21, 863–867. [Google Scholar] [CrossRef]
  25. Cao, Q.; Geng, X.; Wang, H.; Wang, P.; Liu, A.; Lan, Y.; Peng, Q. A review of current development of graphene mechanics. Crystals 2018, 8, 357. [Google Scholar] [CrossRef]
  26. Yang, G.; Li, L.; Lee, W.B.; Ng, M.C. Structure of graphene and its disorders: A review. Sci. Technol. Adv. Mater. 2018, 19, 613–648. [Google Scholar] [CrossRef] [PubMed]
  27. Pires, M.A.; Duarte-Queirós, S.M. Quantum walks with sequential aperiodic jumps. Phys. Rev. E 2020, 102, 012104. [Google Scholar] [CrossRef]
  28. Maciá, E. Aperiodic Structures in Condensed Matter: Fundamentals and Applications; CRC Press: Boca Raton, FL, USA, 2009; p. 130. [Google Scholar]
  29. Zhong, H.; Zhang, Z.; Xu, H.; Qiu, C.; Peng, L.-M. Comparison of mobility extraction methods based on field-effect measurements for graphene. AIP Adv. 2015, 5, 057136. [Google Scholar] [CrossRef]
  30. Xia, F.; Perebeinos, V.; Lin, Y.; Wu, Y.; Avouris, P. The origins and limits of metal-graphene junction resistance. Nat. Nanotechnol. 2011, 6, 179–184, see also supplementary information. [Google Scholar] [CrossRef]
  31. Schroder, D.K. Semiconductor and Material and Device Characterization, 3rd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2005; p. 146. [Google Scholar]
  32. Sutton, A.P. Electronic Structure of Materials; Clarendon Press: Oxford, UK, 1993; p. 58. [Google Scholar]
  33. Abhilash, T.S.; De Alba, R.; Zhelev, N.; Craighead, H.G.; Parpia, J.M. Transfer printing of CVD graphene FETs on patterned substrates. Nanoscale 2015, 7, 14109–14113. [Google Scholar] [CrossRef]
  34. Sánchez, V.; Sánchez, F.; Wang, C. Independent channel method for nanoribbons with dislocation and Fano defects. Phys. Status Solidi B 2021, 258, 2100095. [Google Scholar] [CrossRef]
  35. Nadri, F.; Mardaani, M.; Rabani, H. Semi-analytic study on the conductance of a lengthy armchair honeycomb nanoribbon including vacancies, defects, or impurities. Chin. Phys. B 2019, 28, 017202. [Google Scholar] [CrossRef]
  36. Sánchez, V.; Wang, C. Kubo conductivity in two-dimensional Fibonacci lattices. J. Non-Cryst. Solids 2003, 329, 151–154. [Google Scholar] [CrossRef]
  37. Sánchez, V.; Wang, C. Application of renormalization and convolution methods to the Kubo-Greenwood formula in multidimensional Fibonacci systems. Phys. Rev. B 2004, 70, 144207. [Google Scholar] [CrossRef]
  38. Sánchez, F.; Sánchez, V.; Wang, C. Renormalization approach to the electronic localization and transport in macroscopic generalized Fibonacci lattices. J. Non-Cryst. Solids 2016, 450, 194–208. [Google Scholar] [CrossRef]
Figure 1. Schematic representations of (a) top and (b,c) side views of a backgated GFET with Dirac cones and energy band diagrams for (b) VG = 0 and (c) VG > 0, where eVG = eVC + µ being µ the chemical potential of graphene ribbon with respect to the Dirac point.
Figure 1. Schematic representations of (a) top and (b,c) side views of a backgated GFET with Dirac cones and energy band diagrams for (b) VG = 0 and (c) VG > 0, where eVG = eVC + µ being µ the chemical potential of graphene ribbon with respect to the Dirac point.
Nanomaterials 12 03223 g001
Figure 2. Sketch of (a) a graphene ribbon of width W and length L with hopping integrals t, t s and t l , being l = 1 , 2 , , L 1 , connected to two semi-infinite periodic leads (light grey spheres) and (b) double chains obtained from the independent dual-channel method (Appendix A) with hopping integrals ±   t , ±   t l and on-site energies ±   ε j being j = 1 ,   2 , ,   W . Red highlight and gray zones respectively indicate zig-zag lines and possible transversal dislocations.
Figure 2. Sketch of (a) a graphene ribbon of width W and length L with hopping integrals t, t s and t l , being l = 1 , 2 , , L 1 , connected to two semi-infinite periodic leads (light grey spheres) and (b) double chains obtained from the independent dual-channel method (Appendix A) with hopping integrals ±   t , ±   t l and on-site energies ±   ε j being j = 1 ,   2 , ,   W . Red highlight and gray zones respectively indicate zig-zag lines and possible transversal dislocations.
Nanomaterials 12 03223 g002
Figure 3. (a) Electronic density of states (DOS) and (b) electrical conductance (g) at T = 0 in units of g 0 versus the chemical potential (μ) for graphene ribbons of W = 0.122   μ m (red lines), 0.073   μ m (dark yellow lines), and 0.037   μ m (blue lines) with an arbitrary length of L transversal zig-zag lines and hopping integrals of t and t s = 1.2   t .
Figure 3. (a) Electronic density of states (DOS) and (b) electrical conductance (g) at T = 0 in units of g 0 versus the chemical potential (μ) for graphene ribbons of W = 0.122   μ m (red lines), 0.073   μ m (dark yellow lines), and 0.037   μ m (blue lines) with an arbitrary length of L transversal zig-zag lines and hopping integrals of t and t s = 1.2   t .
Nanomaterials 12 03223 g003
Figure 4. Normalized electrical conductance (g) by the conductance quantum g 0 at T = 0 and at k B T = 0.02 |   t   | as a function of the gate voltage ( V G ) for GFETs with (a) C = 10 nF / cm 2 and (b) C = 50 nF / cm 2 , whose dimensions are W = 0.122   μ m and L = 4097 zig-zag lines without dislocations.
Figure 4. Normalized electrical conductance (g) by the conductance quantum g 0 at T = 0 and at k B T = 0.02 |   t   | as a function of the gate voltage ( V G ) for GFETs with (a) C = 10 nF / cm 2 and (b) C = 50 nF / cm 2 , whose dimensions are W = 0.122   μ m and L = 4097 zig-zag lines without dislocations.
Nanomaterials 12 03223 g004
Figure 5. Electrical conductance (g) at k B T   = 0.026 |   t   | as a function of gate voltage ( V G ) and dislocation hopping integral t A for the same GFET of Figure 4a except by 72 A-type transversal dislocations lines placed following the dilute RS sequence with m = 5 and t B = t C = t D = t .
Figure 5. Electrical conductance (g) at k B T   = 0.026 |   t   | as a function of gate voltage ( V G ) and dislocation hopping integral t A for the same GFET of Figure 4a except by 72 A-type transversal dislocations lines placed following the dilute RS sequence with m = 5 and t B = t C = t D = t .
Nanomaterials 12 03223 g005
Figure 6. Conductance (g) at (a) k B T = 0.01 |   t   | and (b) k B T = 0.02 |   t   | as functions of the gate voltage ( V G ) for GFETs with dislocations and a capacitance C = 10 nF / cm 2 , W = 0.122   μ m (i.e., W = 500 dual channels) and L = 0.217   μ m (red lines), 0.433   μ m (green lines), and 0.866   μ m (blue lines), whose hopping integrals are t A = 2.389   eV , t B = t C = 2.74   eV and t D = 2.63   eV .
Figure 6. Conductance (g) at (a) k B T = 0.01 |   t   | and (b) k B T = 0.02 |   t   | as functions of the gate voltage ( V G ) for GFETs with dislocations and a capacitance C = 10 nF / cm 2 , W = 0.122   μ m (i.e., W = 500 dual channels) and L = 0.217   μ m (red lines), 0.433   μ m (green lines), and 0.866   μ m (blue lines), whose hopping integrals are t A = 2.389   eV , t B = t C = 2.74   eV and t D = 2.63   eV .
Nanomaterials 12 03223 g006
Figure 7. Theoretical source-drain current (red spheres) versus gate voltage ( V G ) and graphene ribbon length ( L ) for a GFET of a width W = 2.2   μ m and a capacitance per unit area C = 12.11 nF / cm 2 , in comparison with experimental data (cyan spheres) reported from Ref. [28] for L = 1 , 2 , , 6 μ m .
Figure 7. Theoretical source-drain current (red spheres) versus gate voltage ( V G ) and graphene ribbon length ( L ) for a GFET of a width W = 2.2   μ m and a capacitance per unit area C = 12.11 nF / cm 2 , in comparison with experimental data (cyan spheres) reported from Ref. [28] for L = 1 , 2 , , 6 μ m .
Nanomaterials 12 03223 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sánchez, F.; Sánchez, V.; Wang, C. Independent Dual-Channel Approach to Mesoscopic Graphene Transistors. Nanomaterials 2022, 12, 3223. https://doi.org/10.3390/nano12183223

AMA Style

Sánchez F, Sánchez V, Wang C. Independent Dual-Channel Approach to Mesoscopic Graphene Transistors. Nanomaterials. 2022; 12(18):3223. https://doi.org/10.3390/nano12183223

Chicago/Turabian Style

Sánchez, Fernando, Vicenta Sánchez, and Chumin Wang. 2022. "Independent Dual-Channel Approach to Mesoscopic Graphene Transistors" Nanomaterials 12, no. 18: 3223. https://doi.org/10.3390/nano12183223

APA Style

Sánchez, F., Sánchez, V., & Wang, C. (2022). Independent Dual-Channel Approach to Mesoscopic Graphene Transistors. Nanomaterials, 12(18), 3223. https://doi.org/10.3390/nano12183223

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop