Band structure of graphene, massless Dirac fermions as low-energy quasiparticles, Berry phase, and all that: Difference between revisions
No edit summary |
|||
| (15 intermediate revisions by the same user not shown) | |||
| Line 1: | Line 1: | ||
== Graphene as the first truly two-dimensional crystal == | == Graphene as the first truly two-dimensional crystal == | ||
The surprising [ | The surprising [https://www.scientificamerican.com/article/carbon-wonderland/ experimental discovery] of a two-dimensional (2D) allotrope of carbon, '''termed graphene''', has ushered unforeseen avenues to explore transport and interactions of low-dimensional electron system, build quantum-coherent carbon-based nanoelectronic devices, and probe high-energy physics of "charged neutrinos" in table-top experiments. Graphene represents one-atom–thick layer of carbon atoms tightly packed into a honeycomb crystal lattice. Its discovery is unexpected in the view of [https://www.pnas.org/content/114/8/1856 Mermin-Wagner theorem] which prevents breaking of continuous symmetry in less than three-dimensions, which means that at first sight 2D crystals (of astronomical size though) would be unstable at any finite temperature. | ||
== Basic chemistry of graphene == | == Basic chemistry of graphene == | ||
| Line 123: | Line 123: | ||
</math> | </math> | ||
Note that the phase factors <math> 1 + e^{i \mathbf{k} \cdot \mathbf{a}} + e^{i \mathbf{k} \cdot \mathbf{a}} </math> become the 1/3 roots of unity, <math> 1 + e^{i2pi/3} + e^{-i2pi/3} = 0 </math>, when <math> \mathbf{k} </math> lies on | Note that the phase factors <math> 1 + e^{i \mathbf{k} \cdot \mathbf{a}} + e^{i \mathbf{k} \cdot \mathbf{a}} </math> become the 1/3 roots of unity, <math> 1 + e^{i2pi/3} + e^{-i2pi/3} = 0 </math>, when <math> \mathbf{k} </math> lies on corner points of the Brillouin zone, K and K' (the other four are equivalent to those as they can be reached from them by adding the reciprocal lattice vector). Therefore, the energy <math> e(\mathbf{k}) = 0 </math> at the zone corners. Everywhere else in k-space, <math> e(\mathbf{k}) \neq 0 </math> and the splitting of two graphene bands [<math> \pi </math> and <math> \pi^* </math> or conduction and valence, respectively, in Fig. 4(c)] lying symmetrically above and below the Fermi energy <math> \varepsilon_F=0 </math> is <math> 2 |e(\mathbf{k})| </math>. | ||
A perfect undoped graphene sheet has one electron per carbon in the <math> \pi </math> level. Therefore, the Fermi level is between the two symmetrical bands, with zero excitation energy needed to excite an electron from just below the Fermi energy to just above at the <math> K </math> points in the corners of the Brillouin zone. Note that the <math> (\pi \pi^*) </math> degeneracy at isolated K points forming the Fermi "surface" is a general to the single particle description of graphene. It follows from symmetry, and it is not just accidental result of the simply nearest-neighbor tight-binding model. For example, in density functional theory calculations (see [http://www.physics.udel.edu/~bnikolic/teaching/phys824/lectures/intro_dft.pdf Lecture | A perfect undoped graphene sheet has one electron per carbon in the <math> \pi </math> level. Therefore, the Fermi level is between the two symmetrical bands, with zero excitation energy needed to excite an electron from just below the Fermi energy to just above at the <math> K </math> points in the corners of the Brillouin zone. Note that the <math> (\pi \pi^*) </math> degeneracy at isolated K points forming the Fermi "surface" is a general to the single particle description of graphene. It follows from symmetry, and it is not just accidental result of the simply nearest-neighbor tight-binding model. For example, in density functional theory calculations (see [http://www.physics.udel.edu/~bnikolic/teaching/phys824/lectures/intro_dft.pdf Lecture 5]), there is no particle-hole symmetry, but the degeneracy at the K and K' points persists. | ||
== | == Massless Dirac fermions as low-energy quasiparticles and their Berry phase == | ||
The wavefunctios of low-energy quasiparticles (i.e., for energy close to the Fermi energy or charge neutrality point <math> | The wavefunctios of low-energy quasiparticles (i.e., for energy close to the Fermi energy or charge neutrality point <math> \varepsilon_F=0 </math>) have attracted a lot of interest. Let us consider a circular path in k-space around <math> K </math> point. The energy is a linear function of <math> \delta \mathbf{k} = \mathbf{k}-\mathbf{K} </math>. If we move adiabatically in k-space around this point, the wavefunctions acquire a "Berry phase" <math> e^{i\pi}=-1 </math> when completing a circle. This can be seen by expanding the wavefunctions to first order in <math> \delta \mathbf{k} </math>: | ||
<math> | <math> | ||
| Line 137: | Line 137: | ||
where <math> \theta </math> is the angle from the x-axis in the <math> k_x, k_y </math> plane. Here the formulas <math> \mathbf{a} + \mathbf{b} = \sqrt{3} a \mathbf{i} </math> and <math> \mathbf{a} - \mathbf{b} = - a \mathbf{j} </math> were employed. Thus, <math> |e(\mathbf{k})| = \sqrt{3} |\delta\mathbf{k}| a t </math>, and <math> e(\mathbf{k}) </math> has phase <math> \phi = - \theta - \pi/2 </math>. Since the wavefunctions have phases <math> \pm \phi/2 </math>, they change phase by <math> \pi </math> when <math> \theta </math> increase by <math> 2 \pi </math>, that is, when the <math> \mathbf{k}</math>-vector goes once around a loop surrounding a <math> K </math> point. | where <math> \theta </math> is the angle from the x-axis in the <math> k_x, k_y </math> plane. Here the formulas <math> \mathbf{a} + \mathbf{b} = \sqrt{3} a \mathbf{i} </math> and <math> \mathbf{a} - \mathbf{b} = - a \mathbf{j} </math> were employed. Thus, <math> |e(\mathbf{k})| = \sqrt{3} |\delta\mathbf{k}| a t </math>, and <math> e(\mathbf{k}) </math> has phase <math> \phi = - \theta - \pi/2 </math>. Since the wavefunctions have phases <math> \pm \phi/2 </math>, they change phase by <math> \pi </math> when <math> \theta </math> increase by <math> 2 \pi </math>, that is, when the <math> \mathbf{k}</math>-vector goes once around a loop surrounding a <math> K </math> point. | ||
The code word | The code word '''massless Dirac fermion''' is used to refer not only to linear <math> \varepsilon(\delta\mathbf{k}) = \pm v_F \hbar \sqrt {\delta \mathbf{k}_x^2 + \delta \mathbf{k}_y^2} </math> dispersion, but it also takes into account that the quasiparticle Hamiltonian around the '''Dirac points''' K and K' has a form of the Weyl Hamiltonian, i.e., the Dirac Hamiltonian of a relativistic particle with zero mass: | ||
<math> | <math> | ||
\hat{H}(\delta \mathbf{k}) = -\frac{\sqrt{3}}{2} a t \frac{1}{\hbar} | \hat{H}(\delta \mathbf{k}) = -\frac{\sqrt{3}}{2} a t \frac{1}{\hbar} | ||
| Line 143: | Line 143: | ||
0 & \hbar (\delta k_x - i \delta k_y)\\ | 0 & \hbar (\delta k_x - i \delta k_y)\\ | ||
\hbar (\delta k_x + i \delta k_y) & 0 \\ | \hbar (\delta k_x + i \delta k_y) & 0 \\ | ||
\end{pmatrix} = -v_F \hat{\boldsymbol \sigma} \cdot \hat{\mathbf{p}} | \end{pmatrix} = -v_F \hat{\boldsymbol \sigma} \cdot \hat{\mathbf{p}} | ||
</math> | </math> | ||
where the Fermi velocity | where the Fermi velocity <math> v_F = \frac{\sqrt{3}}{2} a t \approx c/300 </math> replaces the velocity of light <math> c </math> when the same Hamiltonian is used in quantum field theory of elementary particles (such as neutrinos). | ||
== Pseudospin, isospin and chirality == | == Pseudospin, isospin and chirality of massless Dirac fermions == | ||
Although graphene linear spectrum is important, it is not the only essential feature. Above zero energy, the current carrying | Although graphene linear spectrum is important, it is not the only essential feature. Above zero energy, the current carrying | ||
states in graphene are, as usual, electron-like and negatively charged. At negative energies, if the valence band is not full, unoccupied electronic states behave as positively charged quasiparticles (holes), which are often viewed as a condensed matter equivalent of positrons. Note, however, that electrons and holes in condensed matter physics are normally described by separate Schrödinger equations, which are not in any way connected (as a consequence of the so-called Seitz sum rule, the equations should also involve different effective masses). | states in graphene are, as usual, electron-like and negatively charged. At negative energies, if the valence band is not full, unoccupied electronic states behave as positively charged quasiparticles (holes), which are often viewed as a condensed matter equivalent of positrons. Note, however, that electrons and holes in condensed matter physics are normally described by separate Schrödinger equations, which are not in any way connected (as a consequence of the so-called Seitz sum rule, the equations should also involve different effective masses). | ||
In contrast, electron and hole states in graphene should be interconnected, exhibiting properties analogous to the charge-conjugation symmetry in quantum electrodynamics (QED). For the case of graphene, the latter symmetry is a consequence of the crystal symmetry, because graphene low-energy quasiparticles have to be described by two-component wave functions discussed above, which are needed to define the relative contributions of the <math> \mathrm{A} </math> and <math> \mathrm{B} </math> sublattices in the quasiparticles make-up. The two-component description for graphene is very similar to the spinor wave functions in QED, but the ‘spin’ index for graphene indicates the sublattice rather than the real spin of the electrons and is usually referred to as | In contrast, electron and hole states in graphene should be interconnected, exhibiting properties analogous to the charge-conjugation symmetry in quantum electrodynamics (QED). For the case of graphene, the latter symmetry is a consequence of the crystal symmetry, because graphene low-energy quasiparticles have to be described by two-component wave functions discussed above, which are needed to define the relative contributions of the <math> \mathrm{A} </math> and <math> \mathrm{B} </math> sublattices in the quasiparticles make-up. The two-component description for graphene is very similar to the spinor wave functions in QED, but the ‘spin’ index for graphene indicates the sublattice rather than the real spin of the electrons and is usually referred to as isospin. | ||
[[Image:isospin_pseudospin.png||thumb|center|700px|FIG 6: Effective continuous (i.e., long wave length) Hamiltonian of graphene with both isospin and pseudospin degrees of freedom taken into account. '''Due to the isospin conservation, A-B symmetric potential cannot backward scatter chiral fermions.''']] | [[Image:isospin_pseudospin.png||thumb|center|700px|FIG 6: Effective continuous (i.e., long wave length) Hamiltonian of graphene with both isospin and pseudospin degrees of freedom taken into account. '''Due to the isospin conservation, A-B symmetric potential cannot backward scatter chiral fermions.''']] | ||
This allows one to introduce ''chirality'' – formally a projection of isospin on the direction of motion – which is positive and negative for electrons and holes, respectively. The description of the electron spectrum of graphene in terms of Dirac massless fermions is a kind of continuum-medium description applicable for electron wavelengths much larger than interatomic distances. However, even at these length scales, there is still some retention of the structure of the elementary cell, that is, the existence of two sublattices. In terms of continuum field theory, this can be described only as an internal degree of freedom of the charge carriers, which is just the chirality. | This allows one to introduce ''chirality'' – formally a projection of isospin on the direction of motion – which is positive and negative for electrons and holes, respectively. ''The description of the electron spectrum of graphene in terms of Dirac massless fermions is a kind of continuum-medium description applicable for electron wavelengths much larger than interatomic distances. However, even at these length scales, there is still some retention of the structure of the elementary cell, that is, the existence of two sublattices.'' In terms of continuum field theory, this can be described only as an internal degree of freedom of the charge carriers, which is just the chirality. | ||
The concepts of | The concepts of isospin and chirality, i.e., the fact that isospin direction is linked to the axis determined by the electron momentum, allows us to understand some transport experiments. For example, QED predicts the so-called Klein paradox - '''unimpeded penetration of relativistic particles through high and wide potential barriers''' - as one of the most exotic and counterintuitive phenomena. Although discussed in many contexts in particle, nuclear and astro-physics, direct observation of the Klein paradox using elementary particles have so far proved impossible. On the other hand, it has been observed in recent [http://www.nature.com/nphys/journal/v5/n3/abs/nphys1198.html experiments] on graphene. | ||
The matching between directions of | The matching between directions of chirality illustrated in the Fig. 6 for quasiparticles | ||
inside and outside the barrier results in perfect | inside and outside the barrier results in perfect tunneling. In the strictly one-dimensional case, such perfect transmission of Dirac fermions has been discussed in the context of electron transport in carbon nanotubes as well. | ||
[[Image:klein.png||thumb|center|700px|FIG 6: Tunnelling through a potential barrier in graphene: (a) Schematic diagrams of the spectrum of quasiparticles in single-layer graphene. The spectrum is linear at low Fermi energies (< 1 eV). The red and green curves emphasize the origin of the linear spectrum, which is the crossing between the energy bands associated with crystal sublattices <math> \mathrm{A} </math> and <math> \mathrm{B} </math>. (b) Potential barrier of height <math> V_0 </math> and width <math> D </math>. The three diagrams in a schematically show the positions of the Fermi energy E across such a barrier. The Fermi level (dotted lines) lies in the conduction band outside the barrier and the valence band inside it. The blue filled areas indicate occupied states. The | [[Image:klein.png||thumb|center|700px|FIG 6: Tunnelling through a potential barrier in graphene: (a) Schematic diagrams of the spectrum of quasiparticles in single-layer graphene. The spectrum is linear at low Fermi energies (< 1 eV). The red and green curves emphasize the origin of the linear spectrum, which is the crossing between the energy bands associated with crystal sublattices <math> \mathrm{A} </math> and <math> \mathrm{B} </math>. (b) Potential barrier of height <math> V_0 </math> and width <math> D </math>. The three diagrams in a schematically show the positions of the Fermi energy E across such a barrier. The Fermi level (dotted lines) lies in the conduction band outside the barrier and the valence band inside it. The blue filled areas indicate occupied states. The isospin denoted by vector <math> \sigma </math> is parallel (antiparallel) to the direction of motion of electrons (holes), which also means that <math> \sigma </math> keeps a fixed direction along the red and green branches of the electronic spectrum. Adapted from [http://www.nature.com/nphys/journal/v2/n9/full/nphys384.html Nature Physics '''2''', 620 (2006).]]] | ||
Latest revision as of 15:34, 8 October 2020
Graphene as the first truly two-dimensional crystal
The surprising experimental discovery of a two-dimensional (2D) allotrope of carbon, termed graphene, has ushered unforeseen avenues to explore transport and interactions of low-dimensional electron system, build quantum-coherent carbon-based nanoelectronic devices, and probe high-energy physics of "charged neutrinos" in table-top experiments. Graphene represents one-atom–thick layer of carbon atoms tightly packed into a honeycomb crystal lattice. Its discovery is unexpected in the view of Mermin-Wagner theorem which prevents breaking of continuous symmetry in less than three-dimensions, which means that at first sight 2D crystals (of astronomical size though) would be unstable at any finite temperature.
Basic chemistry of graphene
Carbon has four valence electrons, of which three are used for the bonds in graphene, graphite, and carbon nanotubes. The origin of covalent bonding between two carbon atoms in the Fig. 1 is tunneling of electrons through the potential barrier between the atoms which removes the degeneracy of the original two-level system, making spatially symmetric wave function (the bonding one) to have lower energy eigenvalue of the corresponding Hamiltonian.
When many carbon atoms are brought into a crystal, each bonding or antibonding orbital will acquire dispersion to become an energy band. The analysis below is concerned with the band structure which originates from orbital hosting the fourth valence electron. The bands which correspond to the dispersion of bonding and antibonding molecular orbital (constructed from orbitals on two carbon atoms) are called and Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \pi^* } bands.
Lattice structure of graphene
Real space honeycomb lattice: The honeycomb lattice of graphene show in Fig. 2(a) is not a Bravais lattice. Instead, it can be viewed as bipartite lattice composed of two interpenetrating triangular sublattices. Thus, such lattice with a basis contains two identical carbon atoms per primitive unit cell.
Reciprocal space honeycomb lattice: The corners of the Brillouin zone are at the points Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{K}_i } , which can be expressed as:
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{K}_1 = (\mathbf{G}_2 - \mathbf{G}_1)/3 }
,
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{K}_2 = (2\mathbf{G}_2 + \mathbf{G}_1)/3 }
,
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{K}_3 = (2\mathbf{G}_1 + \mathbf{G}_2)/3 }
,
etc.
Among six corners, only two are inequivalent. For example, notice that Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{K}_3 = \mathbf{K}_1 + \mathbf{G}_1 } .
Symmetries: Because the lattice is two-dimensional, all translations commute with reflection in the plane of the lattice, so all electron (or vibrational) eigenstates can be chosen to be either even or odd under this reflection. For this reason, the single-particle electron states are rigorously separated into two classes, called Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \sigma } and Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \pi } . The even Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \sigma } states are derived from carbon Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle s } and Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle p_x } , Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle p_y } orbitals (i.e., their hybridized Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle sp^2 } orbitals discussed above), while the odd Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \pi } states are derived from carbon Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle p_z } orbitals. These latter are cylindrically symmetric in the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle xy } plane, lie near the Fermi level (half-filled) and are the electrically active states of interest in low energy experimental probing of graphene.
Tight-binding Hamiltonian of graphene
A useful picture of electron behavior can be derived by using the phenomenological nearest-neighbor tight-binding model to look at the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \pi } electrons. The two sublattices are shown below in different colors, with the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathrm{A} } sublattice at vectors Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{R} = n_1 \mathbf{a} + n_2 \mathbf{b} } , and the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathrm{B} } sublattice at vectors Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{R} + \mathbf{\tau} } with Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{\tau} = -(\mathbf{a} + \mathbf{b})/3 } .
The tight-binding model is then defined as:
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \hat{H} = -t \sum_\mathbf{R} \left( |\mathbf{R}\rangle \langle \mathbf{R} + \mathbf{\tau} | + |\mathbf{R}\rangle \langle \mathbf{R} + \mathbf{a} + \mathbf{\tau} | + |\mathbf{R}\rangle \langle \mathbf{R} + \mathbf{b} + \mathbf{\tau} | + \mathrm{H.c.} \right). }
Here Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \langle \mathbf{r}|\mathbf{R} \rangle = \Psi_{p_z}(\mathbf{R}-\mathbf{r}) } is a wave function of Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle p_z } orbital on the sublattice Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathrm{A} } atom at site Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{R} } , while Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \langle \mathbf{r}|\mathbf{R}+\mathbf{\tau} \rangle } is similar state on the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathrm{B} } sublattice atom, and Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle t } is the "hopping integral" from a state to an adjacent similar state.
In the simplest model discussed here, the hopping matrix element couples states on the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathrm{A} } sublattice only to states on the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathrm{B} } sublattice, and vice versa. The hopping is chosen as Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle t \approx 2.7 } eV in order to match the first principles computed band structure around the corners of the Brillouin zone (to match the density functional theory computed band structure over the whole Brillouin zone, one has to use the third-nearest neighbor tight-binding model - see Lecture 5 for more details).
Diagonalization of the tight-binding model of graphene: LCAO method
The symmetry arguments above show why we were able to truncate the basis to a single Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle p_z } orbital per atom, as long as we cay about the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \pi } band near the Fermi level. Since there are two atoms in the unit cell, we can form two candidate Bloch orbitals:
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle |\mathbf{k}A\rangle = \frac{1}{\sqrt{N}} \sum_{\mathbf{R}} e^{i\mathbf{k}\cdot\mathbf{R}} |\mathbf{R}\rangle, }
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle |\mathbf{k}B\rangle = \frac{1}{\sqrt{N}} \sum_{\mathbf{R}} e^{i\mathbf{k}\cdot\mathbf{R}} |\mathbf{R} + \mathbf{\tau} \rangle. }
The terms "Bloch" indicates that they are eingestates of the translation operator. Thus, for the basis functions from which we can construct approximate (in the sense that we assume that no other orbitals can mix into the energy eigenstates near the Fermi energy) Bloch eigenstates of the Hamiltonian:
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle |\mathbf{k} \rangle = \alpha |\mathbf{k} A\rangle + \beta |\mathbf{k} B\rangle. }
where the expansion coefficients (normalized by Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle |\alpha|^2 + |\beta|^2 =1 } are eigenvectors of the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle 2 \times 2 } Hamiltonian matrix:
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \hat{H}(\mathbf{k}) = \begin{pmatrix} \langle \mathbf{k} A |\hat{H}|\mathbf{k}A\rangle & \langle \mathbf{k} A |\hat{H}|\mathbf{k}B\rangle \\ \langle \mathbf{k} B |\hat{H}|\mathbf{k}A\rangle & \langle \mathbf{k} B |\hat{H}|\mathbf{k}B\rangle \\ \end{pmatrix} }
The diagonal matrix elements are zero:
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \langle \mathbf{k} A |\hat{H}|\mathbf{k}A\rangle = \langle \mathbf{k}B |\hat{H}|\mathbf{k}B\rangle = 0, }
while the off-diagonal are given by:
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \langle \mathbf{k} A |\hat{H}|\mathbf{k} B \rangle = \frac{1}{N} \sum_\mathbf{RR^{\prime}} e^{i\mathbf{k} \cdot (\mathbf{R^\prime} - \mathbf{R})} \langle \mathbf{R} | \hat{H} | \mathbf{R} + \mathbf{\tau} \rangle = - t \left( 1 + e^{i\mathbf{k} \cdot \mathbf{a}} + e^{i\mathbf{k}\cdot \mathbf{b}} \right) = e(\mathbf{k}). }
The single particle Bloch eigenenergies are thus Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \varepsilon(k_x,k_y) = \pm |e(\mathbf{k})| }
, where
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle |e(\mathbf{k})|/t = \sqrt{1 + 4 \cos \left(\frac{\sqrt{3} k_x a}{2} \right)\cos\left(\frac{k_y a}{2} \right) + 4 \cos^2 \left( \frac{k_y a}{2} \right)}. }
This band structure can be plotted in a usual way (the so-called "spaghetti diagram") shown in Fig. 4(b), along lines defined by the high symmetry points on the edges of the "irreducible wedge" in Fig. 4(a). Energy at any other Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{k} } -vector outside of the wedge is obtained from Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \varepsilon(U\mathbf{k})=\varepsilon(\mathbf{k}) } , where Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle U } is the point group symmetry operation of the crystal. Similarly, energy at any Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{k} } -vector outside of the First Brillouin zone is obtained from Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \varepsilon(T\mathbf{k})=\varepsilon(\mathbf{k}) } , where Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle T } is the translation group symmetry operation of the crystal. In fact, since band dispersion Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \varepsilon(k_x,k_y) } for 2D crystal is just a function of two variables, we can plot its complete surface as shown in Fig. 4(c).
Let us write Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle e(\mathbf{k}) = |e(\mathbf{k})| e^{i\phi(\mathbf{k})} } . The Schrödinger equation is then:
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \hat{H}(\mathbf{k}) |\Psi_\mathbf{k}\rangle = \varepsilon(\mathbf{k}) |\Psi_\mathbf{k}\rangle }
with the Hamiltonian matrix being:
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \hat{H}(\mathbf{k}) = \begin{pmatrix} 0 & |e(\mathbf{k})|e^{i\phi(\mathbf{k})} \\ |e(\mathbf{k})|e^{-i\phi(\mathbf{k})} & 0 \\ \end{pmatrix}. }
The eigenvectors Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle |\Psi_\mathbf{k} \rangle } are given by:
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle | \mathbf{k} - \rangle = \frac{1}{\sqrt{2}} \begin{pmatrix} e^{i\phi(\mathbf{k})/2} \\ -e^{-i\phi(\mathbf{k})/2} \\ \end{pmatrix}, | \mathbf{k}+ \rangle = \frac{1}{\sqrt{2}} \begin{pmatrix} e^{i\phi(\mathbf{k})/2} \\ e^{-i\phi(\mathbf{k})/2} \\ \end{pmatrix}, }
Note that the phase factors Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle 1 + e^{i \mathbf{k} \cdot \mathbf{a}} + e^{i \mathbf{k} \cdot \mathbf{a}} } become the 1/3 roots of unity, Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle 1 + e^{i2pi/3} + e^{-i2pi/3} = 0 } , when Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{k} } lies on corner points of the Brillouin zone, K and K' (the other four are equivalent to those as they can be reached from them by adding the reciprocal lattice vector). Therefore, the energy Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle e(\mathbf{k}) = 0 } at the zone corners. Everywhere else in k-space, Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle e(\mathbf{k}) \neq 0 } and the splitting of two graphene bands [Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \pi } and Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \pi^* } or conduction and valence, respectively, in Fig. 4(c)] lying symmetrically above and below the Fermi energy Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \varepsilon_F=0 } is Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle 2 |e(\mathbf{k})| } .
A perfect undoped graphene sheet has one electron per carbon in the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \pi } level. Therefore, the Fermi level is between the two symmetrical bands, with zero excitation energy needed to excite an electron from just below the Fermi energy to just above at the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle K } points in the corners of the Brillouin zone. Note that the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle (\pi \pi^*) } degeneracy at isolated K points forming the Fermi "surface" is a general to the single particle description of graphene. It follows from symmetry, and it is not just accidental result of the simply nearest-neighbor tight-binding model. For example, in density functional theory calculations (see Lecture 5), there is no particle-hole symmetry, but the degeneracy at the K and K' points persists.
Massless Dirac fermions as low-energy quasiparticles and their Berry phase
The wavefunctios of low-energy quasiparticles (i.e., for energy close to the Fermi energy or charge neutrality point Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \varepsilon_F=0 } ) have attracted a lot of interest. Let us consider a circular path in k-space around Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle K } point. The energy is a linear function of Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \delta \mathbf{k} = \mathbf{k}-\mathbf{K} } . If we move adiabatically in k-space around this point, the wavefunctions acquire a "Berry phase" Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle e^{i\pi}=-1 } when completing a circle. This can be seen by expanding the wavefunctions to first order in Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \delta \mathbf{k} } :
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle -\frac{e(\mathbf{k})}{t} = 1 + e^{(\mathbf{K}+\delta \mathbf{k})\cdot \mathbf{a}} + e^{(\mathbf{K}+\delta \mathbf{k}) \cdot \mathbf{b}} \approx 1 + \left( -\frac{1}{2} -i\frac{\sqrt{3}}{2} \right) (1+i\delta \mathbf{k} \cdot \mathbf{a}) + \left( - \frac{1}{2} + i \frac{\sqrt{3}}{2} \right) ( 1+ i \delta\mathbf{k} \cdot \mathbf{b}) \approx -\frac{1}{2} i \delta \mathbf{k} \cdot (\mathbf{a} + \mathbf{b}) + \frac{\sqrt{3}}{2} i \delta \mathbf{k} \cdot (\mathbf{a} - \mathbf{b}) = \frac{\sqrt{3}}{2}|\delta \mathbf{k}| a e^{-i\pi/2} e^{-i\theta}, }
where Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \theta } is the angle from the x-axis in the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle k_x, k_y } plane. Here the formulas Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{a} + \mathbf{b} = \sqrt{3} a \mathbf{i} } and Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{a} - \mathbf{b} = - a \mathbf{j} } were employed. Thus, Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle |e(\mathbf{k})| = \sqrt{3} |\delta\mathbf{k}| a t } , and Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle e(\mathbf{k}) } has phase Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \phi = - \theta - \pi/2 } . Since the wavefunctions have phases Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \pm \phi/2 } , they change phase by Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \pi } when Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \theta } increase by Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle 2 \pi } , that is, when the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathbf{k}} -vector goes once around a loop surrounding a Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle K } point.
The code word massless Dirac fermion is used to refer not only to linear Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \varepsilon(\delta\mathbf{k}) = \pm v_F \hbar \sqrt {\delta \mathbf{k}_x^2 + \delta \mathbf{k}_y^2} } dispersion, but it also takes into account that the quasiparticle Hamiltonian around the Dirac points K and K' has a form of the Weyl Hamiltonian, i.e., the Dirac Hamiltonian of a relativistic particle with zero mass:
Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \hat{H}(\delta \mathbf{k}) = -\frac{\sqrt{3}}{2} a t \frac{1}{\hbar} \begin{pmatrix} 0 & \hbar (\delta k_x - i \delta k_y)\\ \hbar (\delta k_x + i \delta k_y) & 0 \\ \end{pmatrix} = -v_F \hat{\boldsymbol \sigma} \cdot \hat{\mathbf{p}} }
where the Fermi velocity Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle v_F = \frac{\sqrt{3}}{2} a t \approx c/300 } replaces the velocity of light Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle c } when the same Hamiltonian is used in quantum field theory of elementary particles (such as neutrinos).
Pseudospin, isospin and chirality of massless Dirac fermions
Although graphene linear spectrum is important, it is not the only essential feature. Above zero energy, the current carrying states in graphene are, as usual, electron-like and negatively charged. At negative energies, if the valence band is not full, unoccupied electronic states behave as positively charged quasiparticles (holes), which are often viewed as a condensed matter equivalent of positrons. Note, however, that electrons and holes in condensed matter physics are normally described by separate Schrödinger equations, which are not in any way connected (as a consequence of the so-called Seitz sum rule, the equations should also involve different effective masses).
In contrast, electron and hole states in graphene should be interconnected, exhibiting properties analogous to the charge-conjugation symmetry in quantum electrodynamics (QED). For the case of graphene, the latter symmetry is a consequence of the crystal symmetry, because graphene low-energy quasiparticles have to be described by two-component wave functions discussed above, which are needed to define the relative contributions of the Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathrm{A} } and Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \mathrm{B} } sublattices in the quasiparticles make-up. The two-component description for graphene is very similar to the spinor wave functions in QED, but the ‘spin’ index for graphene indicates the sublattice rather than the real spin of the electrons and is usually referred to as isospin.
This allows one to introduce chirality – formally a projection of isospin on the direction of motion – which is positive and negative for electrons and holes, respectively. The description of the electron spectrum of graphene in terms of Dirac massless fermions is a kind of continuum-medium description applicable for electron wavelengths much larger than interatomic distances. However, even at these length scales, there is still some retention of the structure of the elementary cell, that is, the existence of two sublattices. In terms of continuum field theory, this can be described only as an internal degree of freedom of the charge carriers, which is just the chirality.
The concepts of isospin and chirality, i.e., the fact that isospin direction is linked to the axis determined by the electron momentum, allows us to understand some transport experiments. For example, QED predicts the so-called Klein paradox - unimpeded penetration of relativistic particles through high and wide potential barriers - as one of the most exotic and counterintuitive phenomena. Although discussed in many contexts in particle, nuclear and astro-physics, direct observation of the Klein paradox using elementary particles have so far proved impossible. On the other hand, it has been observed in recent experiments on graphene.
The matching between directions of chirality illustrated in the Fig. 6 for quasiparticles inside and outside the barrier results in perfect tunneling. In the strictly one-dimensional case, such perfect transmission of Dirac fermions has been discussed in the context of electron transport in carbon nanotubes as well.
