View source for Conductance and shot noise of pseudodiffusive and diffusive electron transport through graphene nanoribbons

Using KWANT or your own Python code, setup a two-terminal graphene nanoribbon with zigzag edges (ZGNR) where semi-infinite ZGNR leads of width <math> W </math> are attached to central region of length <math> L </math> and same width  (in other words, if the central region is clean your two-terminal device is just an infinite homogeneous ZGNR). In the calculations below, you can fix the width at some value that would allow you to perform calculations on available computational resources while length will be varied. 

Compute the zero-temperature conductance <math> G = \frac{2e^2}{h} \sum_n T_n </math>, and the Fano factor, <math> F = \frac{\sum_n T_n(1-T_n)}{\sum_n T_n} </math>,    quantifying the shot noise power, <math> S(0)=2F\langle I \rangle</math>, as a function of the length of the central region for two different device setups discussed below.  Here <math> T_n </math> is the transmission eigenvalue which is obtained by diagonalizing Hermitian matrix <math> \mathbf{t} \mathbf{t}^\dagger </math>, where <math> \mathbf{t} </math> is the transmission submatrix of the scattering matrix of the device. 

==Pseudodiffusive transport in clean ZGNRs==
For this part of the Project, the central region in your two-terminal ZGNR device should have geometry <math> W > L </math> and be covered by a gate electrode whose voltage is modeled as an on-site potential. The parameters like the gate voltage and Fermi energy can be set similarly to those in Fig. 5 of [https://doi.org/10.1103/PhysRevB.76.205433 Phys. Rev. B '''76''', 205433 (2007)].

==Diffusive transport and localization in ZGNRs with vacancies==
For the second part of the Project, introduce vacancies by removing carbon atoms from the central region of your two-terminal ZGNR device (where you '''do not''' use anymore gate potential as in the first part of the Project). You can start by using some large concentration <math> n_\mathrm{ad} =10\%</math> of vacancies, and later also explore dependence on this number akin to the analysis  in [https://doi.org/10.1103/PhysRevLett.115.106601 Phys. Rev. Lett. '''115''', 106601 (2015)] (where conductivity of bulk graphene rather than conductance of GNRs was studied). Use the Fermi energy '''very close''' to the Dirac point at <math> E_F=0 </math>, as well as '''exactly''' at this energy, and compute <math> G </math> and <math> F </math> as a function of length <math> L </math> of the central region, as well as as a function of <math> n_\mathrm{ad} =10\%</math> at some fixed length. Due to random distribution of vacancies, results for <math> G </math> and <math> F </math> should be obtained by disorder averaging them over at least 10 different vacancy configurations. Due to non-Gaussian distribution of these quantities over an ensemble of impurity configurations, the best is to plot [https://wiki.physics.udel.edu/wiki_qttg/images/0/0b/Tmt_localization.pdf typical value] as geometric average (instead of the usual arithmetic average):

<math> \overline G = \sqrt[N]{G_1 \cdot G_2 \cdots G_N} </math>    

Note that wire of length <math> L </math> will enter the diffusive transport regime, defined by electronic mean free path <math> \ell \ll L </math> being smaller than <math> L </math>, once the Fano factor reaches <math> F \simeq 1/3 </math> [at least sufficiently away from the Dirac point, see [https://doi.org/10.1103/PhysRevB.77.081410 Phys. Rev. B '''77''', 081410(R) (2008)]]. Discuss whether electrons become [https://physicstoday.scitation.org/doi/10.1063/1.3206091 '''Anderson-localized'''] or not at any of the Fermi energies you examine.