notes / periodic systems AUG 9, 2022 · 2 MIN READ

FCC and Graphene Bands

From a tight-binding problem set, computed with PythTB.

Tight binding is the cheapest honest model of a band structure. Put one orbital on each lattice site, let electrons hop only to nearest neighbors with amplitude tt, and the Hamiltonian in k-space is already diagonal in the band index. For a Bravais lattice with one site per cell the single band is written as

E(k)=ε0+tδeikδE(\mathbf{k}) = \varepsilon_0 + t \sum_{\boldsymbol{\delta}} e^{i\mathbf{k}\cdot\boldsymbol{\delta}}

where the sum runs over the nearest-neighbor vectors δ\boldsymbol{\delta} and ε0\varepsilon_0 is the onsite energy. Everything that follows comes from evaluating this sum on different lattices.

The FCC band

The face-centered cubic lattice has primitive vectors (12,12,0)(\tfrac{1}{2},\tfrac{1}{2},0), (0,12,12)(0,\tfrac{1}{2},\tfrac{1}{2}), (12,0,12)(\tfrac{1}{2},0,\tfrac{1}{2}) and twelve nearest neighbors. With ε0=0\varepsilon_0 = 0 and t=1t = -1, sweeping k\mathbf{k} along the standard path through the high-symmetry points gives the single band below.

Single tight-binding band of the FCC lattice along the Gamma-X-W-L-Gamma-K path

The bandwidth is set by the coordination. More neighbors means more terms in the sum, so close-packed FCC produces a wide band, and the shape between the labeled points is a fingerprint of the lattice geometry that shows up in real FCC metals well before any serious theory is applied.

Graphene and its density of states

Graphene needs the two-site honeycomb cell, and its two bands touch at the K points of the Brillouin zone. Histogramming the eigenvalues over a uniform k-mesh gives the density of states.

Density of states of graphene tight binding, with van Hove peaks flanking a vanishing density at zero energy

Two features carry the physics. The density of states drops to zero at zero energy without a gap, which is the signature of the Dirac points, where the bands touch only at isolated k-points. And the two sharp peaks at E=±tE = \pm t are van Hove singularities, coming from saddle points of the dispersion at the M points. Adding a next-nearest-neighbor hopping (a tenth of tt in the problem set) breaks the symmetry between the two bands, tilting the spectrum so the electron and hole sides are no longer mirror images.

Zigzag nanotubes

Rolling graphene into a cylinder quantizes the momentum around the circumference, so a nanotube’s bands are one-dimensional slices of the graphene dispersion. Whether a slice passes through a Dirac point decides everything. For a (n,0)(n, 0) zigzag tube the allowed transverse momenta hit the Dirac point only when nn is divisible by 3, so those tubes are metallic and the rest are semiconductors. Plotting the lowest subband for n=1,2,4,10,16n = 1, 2, 4, 10, 16, none of which is divisible by 3, the dispersion stays pinched off from zero energy with a small gap, and zooming in near the minimum confirms the gap never closes. Two thirds of all zigzag tubes are semiconductors by this counting argument alone, which is a remarkable amount of device physics to fall out of one hopping parameter.

Feel free to contact me if there is an error or confusion!

tight bindingelectronic structureband theory