# Band structure of graphene, massless Dirac fermions as low-energy quasiparticles, Berry phase, and all that

## 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 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 , which can be expressed as:

, , , etc.

Among six corners, only two are inequivalent. For example, notice that .

**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 and . The even states are derived from carbon and , orbitals (i.e., their hybridized orbitals discussed above), while the odd states are derived from carbon orbitals. These latter are cylindrically symmetric in the 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 electrons. The two sublattices are shown below in different colors, with the sublattice at vectors , and the sublattice at vectors with .

The tight-binding model is then defined as:

```
```

Here is a wave function of orbital on the sublattice atom at site , while is similar state on the sublattice atom, and 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 sublattice only to states on the sublattice, and vice versa. The hopping is chosen as 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 orbital per atom, as long as we cay about the band near the Fermi level. Since there are two atoms in the unit cell, we can form two candidate Bloch orbitals:

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.