Skip to content

Electronic Band Structure | Physics

In the simplest model, conduction electrons move freely in a box of volume VV (the “jellium” model). The allowed wave vectors are:

k=2πL(nx,ny,nz),niZ\mathbf{k} = \frac{2\pi}{L}(n_x, n_y, n_z), \quad n_i \in \mathbb{Z}

The energy spectrum:

ε(k)=2k22me\varepsilon(\mathbf{k}) = \frac{\hbar^2 k^2}{2m_e}

The Fermi wave vector is determined by the electron density n=N/Vn = N/V:

kF=(3π2n)1/3k_F = (3\pi^2 n)^{1/3}

The Fermi energy:

εF=2kF22me\varepsilon_F = \frac{\hbar^2 k_F^2}{2m_e}

For a 3D free electron gas:

g(ε)=V2π2(2me2)3/2εg(\varepsilon) = \frac{V}{2\pi^2}\left(\frac{2m_e}{\hbar^2}\right)^{3/2}\sqrt{\varepsilon}

Derivation. The number of states with kk\lvert\mathbf{k}\rvert \leq k is:

N(k)=2V(2π)34πk33N(k) = 2 \cdot \frac{V}{(2\pi)^3} \cdot \frac{4\pi k^3}{3}

Where the factor of 2 accounts for spin. Differentiating: g(k)dk=dN/dkdk=(Vk2/π2)dkg(k)\,dk = dN/dk\,dk = (Vk^2/\pi^2)\,dk. Converting to energy: g(ε)=g(k)dk/dε=(Vk2/π2)(me/2k)g(\varepsilon) = g(k)\lvert dk/d\varepsilon\rvert = (Vk^2/\pi^2)(m_e/\hbar^2 k). \blacksquare

At the Fermi energy: g(εF)=3N2εFg(\varepsilon_F) = \frac{3N}{2\varepsilon_F}.

The Fermi surface is the surface in k\mathbf{k}-space defined by ε(k)=εF\varepsilon(\mathbf{k}) = \varepsilon_F. For the free electron gas, this is a sphere of radius kFk_F. The shape of the Fermi surface Strongly influences transport properties (conductivity, Hall effect, cyclotron resonance).

In real metals, the periodic potential distorts the Fermi surface from a sphere. At the Brillouin Zone boundaries, band gaps open and the Fermi surface can develop “necks” (connecting to adjacent Zones) or become multiply connected. The topology of the Fermi surface determines whether a material Is a metal or insulator: a material is metallic if the Fermi surface crosses any Brillouin zone Boundary.

The number of electrons per atom determines the filling: 1 electron/atom (e.g., Na, Cu) gives a Nearly spherical Fermi surface well within the first BZ. 2 electrons/atom (e.g., Mg) nearly fills The first BZ and the Fermi surface contacts the zone boundary. 3—4 electrons/atom (e.g., Al, Pb) Produce complex multiply-connected Fermi surfaces.

Theorem 5.1 (Bloch, 1928). The eigenstates of the one-electron Hamiltonian in a periodic Potential V(r+R)=V(r)V(\mathbf{r} + \mathbf{R}) = V(\mathbf{r}) can be written as:

ψnk(r)=eikrunk(r)\psi_{n\mathbf{k}}(\mathbf{r}) = e^{i\mathbf{k}\cdot\mathbf{r}} u_{n\mathbf{k}}(\mathbf{r})

Where unk(r)u_{n\mathbf{k}}(\mathbf{r}) has the periodicity of the lattice: unk(r+R)=unk(r)u_{n\mathbf{k}}(\mathbf{r} + \mathbf{R}) = u_{n\mathbf{k}}(\mathbf{r}).

Proof. The translation operators T^R\hat{T}_{\mathbf{R}} commute with the Hamiltonian H^=22m2+V(r)\hat{H} = -\frac{\hbar^2}{2m}\nabla^2 + V(\mathbf{r}) since VV is periodic. Therefore, the Eigenstates of H^\hat{H} can be chosen as simultaneous eigenstates of all T^R\hat{T}_{\mathbf{R}}:

T^Rψ(r)=ψ(r+R)=cRψ(r)\hat{T}_{\mathbf{R}}\psi(\mathbf{r}) = \psi(\mathbf{r} + \mathbf{R}) = c_{\mathbf{R}}\psi(\mathbf{r})

From the composition rule T^R1T^R2=T^R1+R2\hat{T}_{\mathbf{R}_1}\hat{T}_{\mathbf{R}_2} = \hat{T}_{\mathbf{R}_1 + \mathbf{R}_2}:

cR1+R2=cR1cR2c_{\mathbf{R}_1 + \mathbf{R}_2} = c_{\mathbf{R}_1} c_{\mathbf{R}_2}

The only solution of this functional equation is cR=eikRc_{\mathbf{R}} = e^{i\mathbf{k}\cdot\mathbf{R}}. Therefore ψ(r+R)=eikRψ(r)\psi(\mathbf{r} + \mathbf{R}) = e^{i\mathbf{k}\cdot\mathbf{R}}\psi(\mathbf{r})Which is Satisfied by ψ(r)=eikruk(r)\psi(\mathbf{r}) = e^{i\mathbf{k}\cdot\mathbf{r}}u_{\mathbf{k}}(\mathbf{r}) with uku_{\mathbf{k}} periodic. \blacksquare

Consequences:

  • k\mathbf{k} is defined only up to a reciprocal lattice vector: k\mathbf{k} and k+G\mathbf{k} + \mathbf{G} are equivalent.
  • The energy spectrum consists of bands εn(k)\varepsilon_n(\mathbf{k})Each labelled by a band index nn.
  • Band gaps appear between allowed energy bands.

Starting from the free electron model, a weak periodic potential V(r)=GVGeiGrV(\mathbf{r}) = \sum_{\mathbf{G}} V_{\mathbf{G}} e^{i\mathbf{G}\cdot\mathbf{r}} Opens gaps at the Brillouin zone boundaries where k=k+G\lvert\mathbf{k}\rvert = \lvert\mathbf{k} + \mathbf{G}\rvert (Bragg Condition).

At the zone boundary k=G/2\mathbf{k} = \mathbf{G}/2The gap is:

Δε=2VG\Delta\varepsilon = 2\lvert V_{\mathbf{G}}\rvert

Derivation. Near the zone boundary, the free electron states at k\mathbf{k} and kG\mathbf{k} - \mathbf{G} Are degenerate: εk0=εkG0\varepsilon_{\mathbf{k}}^0 = \varepsilon_{\mathbf{k} - \mathbf{G}}^0. Degenerate Perturbation theory gives:

det(εk0EVGVGεkG0E)=0\det\begin{pmatrix} \varepsilon_{\mathbf{k}}^0 - E & V_{\mathbf{G}} \\ V_{\mathbf{G}}^* & \varepsilon_{\mathbf{k} - \mathbf{G}}^0 - E \end{pmatrix} = 0

At k=G/2\mathbf{k} = \mathbf{G}/2: E=εG/20±VGE = \varepsilon_{\mathbf{G}/2}^0 \pm \lvert V_{\mathbf{G}}\rvert So the gap is 2VG2\lvert V_{\mathbf{G}}\rvert. \blacksquare

The Drude model (1900) treats conduction electrons as a classical ideal gas scattering off Static ions with a mean free time τ\tau (relaxation time).

Equation of motion. Under an electric field E\mathbf{E}:

medvdt=eEmevτm_e\frac{d\mathbf{v}}{dt} = -e\mathbf{E} - \frac{m_e\mathbf{v}}{\tau}

The second term represents a frictional drag with characteristic time τ\tau.

DC conductivity. In steady state (dv/dt=0d\mathbf{v}/dt = 0): vd=eτmeE\mathbf{v}_d = -\frac{e\tau}{m_e}\mathbf{E}. The current density: J=nevd=ne2τmeE\mathbf{J} = -ne\mathbf{v}_d = \frac{ne^2\tau}{m_e}\mathbf{E}.

σ=ne2τme\sigma = \frac{ne^2\tau}{m_e}

AC conductivity. For E(t)=E0eiωt\mathbf{E}(t) = \mathbf{E}_0\,e^{-i\omega t}The Drude model gives:

σ(ω)=ne2τ/me1iωτ=σ01iωτ\sigma(\omega) = \frac{ne^2\tau/m_e}{1 - i\omega\tau} = \frac{\sigma_0}{1 - i\omega\tau}

The real part Re[σ(ω)]=σ01+ω2τ2\mathrm{Re}[\sigma(\omega)] = \frac{\sigma_0}{1 + \omega^2\tau^2} describes absorption, Peaking at ω=0\omega = 0 (the Drude peak). This explains the metallic reflectivity in the infrared.

Hall effect. With B=Bz^\mathbf{B} = B\hat{z} applied, the steady-state equation becomes:

eEmevτev×B=0-e\mathbf{E} - \frac{m_e\mathbf{v}}{\tau} - e\mathbf{v} \times \mathbf{B} = 0

For current J=Jxx^\mathbf{J} = J_x\hat{x}A transverse field EyE_y develops:

RH=EyJxB=1neR_H = \frac{E_y}{J_x B} = -\frac{1}{ne}

This provides a direct measurement of the carrier density nn.

Successes: Ohm”s law (J=σE\mathbf{J} = \sigma\mathbf{E}), Wiedemann—Franz law (κ/σT=π2kB23e2\kappa/\sigma T = \frac{\pi^2 k_B^2}{3e^2}), Hall effect.

Failures: Predicts χT1\chi \propto T^{-1} (Curie law) for magnetic susceptibility, but real Metals have nearly temperature-independent Pauli paramagnetism. Predicts CV=32nkBC_V = \frac{3}{2}nk_B But experiments give CV32nkBC_V \ll \frac{3}{2}nk_B at room temperature.

The Sommerfeld model (1928) corrects the Drude model by treating electrons as a Fermi gas Obeying Fermi—Dirac …/4-statistics-and-probability/2_statistics:

f(ε)=1e(εμ)/kBT+1f(\varepsilon) = \frac{1}{e^{(\varepsilon - \mu)/k_B T} + 1}

At T=0T = 0The chemical potential equals the Fermi energy: μ(0)=εF\mu(0) = \varepsilon_F. At finite TT:

μ(T)=εF[1π212(kBTεF)2+]\mu(T) = \varepsilon_F\left[1 - \frac{\pi^2}{12}\left(\frac{k_B T}{\varepsilon_F}\right)^2 + \cdots\right]

Since εF/kB104\varepsilon_F/k_B \sim 10^4 K for metals, the correction at room temperature is negligible: The chemical potential is essentially constant.

Electronic specific heat. By the Sommerfeld expansion:

Ce=π23kB2g(εF)T=γTC_e = \frac{\pi^2}{3}k_B^2\,g(\varepsilon_F)\,T = \gamma T

Where γ=π22NkB2εF\gamma = \frac{\pi^2}{2}\frac{Nk_B^2}{\varepsilon_F}. At room temperature, only electrons within kBT\sim k_B T of εF\varepsilon_F can be thermally excited, which is a tiny fraction T/TF1/100\sim T/T_F \sim 1/100 of the total. This explains why Ce32NkBC_e \ll \frac{3}{2}Nk_B.

Pauli paramagnetism. The spin susceptibility of a degenerate electron gas:

χP=μ0μB2g(εF)=3μ0μB2N2εF\chi_P = \mu_0\mu_B^2\,g(\varepsilon_F) = \frac{3\mu_0\mu_B^2 N}{2\varepsilon_F}

This is independent of TT (up to corrections of order (T/TF)2(T/T_F)^2), in contrast to the Curie law χ1/T\chi \propto 1/T of the Drude model.

Derivation: Sommerfeld Expansion

To compute thermal averages at low TTWe integrate h(ε)f(ε)h(\varepsilon) f(\varepsilon) where f(ε)=1/(eβ(εμ)+1)f(\varepsilon) = 1/(e^{\beta(\varepsilon - \mu)} + 1) is the Fermi—Dirac distribution and h(ε)h(\varepsilon) Is any smooth function (e.g., density of states times energy).

Define H(ε)=0εh(ε)dεH(\varepsilon) = \int_0^\varepsilon h(\varepsilon')\,d\varepsilon'. Then:

I=0h(ε)f(ε)dε=0dHdεfdε=[Hf]0+0H(ε)(fε)dεI = \int_0^\infty h(\varepsilon)f(\varepsilon)\,d\varepsilon = \int_0^\infty \frac{dH}{d\varepsilon}\,f\,d\varepsilon = [Hf]_0^\infty + \int_0^\infty H(\varepsilon)\left(-\frac{\partial f}{\partial \varepsilon}\right)d\varepsilon

Since f(0)1f(0) \approx 1 and f()=0f(\infty) = 0 And f/ε-\partial f/\partial \varepsilon is sharply peaked At ε=μ\varepsilon = \mu with width kBT\sim k_B TWe expand H(ε)H(\varepsilon) about μ\mu:

I=0μh(ε)dε+π26(kBT)2h(μ)+I = \int_0^\mu h(\varepsilon)\,d\varepsilon + \frac{\pi^2}{6}(k_B T)^2 h'(\mu) + \cdots

For the total energy with h(ε)=εg(ε)h(\varepsilon) = \varepsilon\,g(\varepsilon):

E=0μ0εg(ε)dε+π26(kBT)2ddε[εg(ε)]ε=μ0+E = \int_0^{\mu_0} \varepsilon\,g(\varepsilon)\,d\varepsilon + \frac{\pi^2}{6}(k_B T)^2 \frac{d}{d\varepsilon}[\varepsilon g(\varepsilon)]_{\varepsilon = \mu_0} + \cdots

Differentiating with respect to TT gives the specific heat CV=π23kB2g(εF)TC_V = \frac{\pi^2}{3}k_B^2\,g(\varepsilon_F)\,T. \blacksquare

Worked Example: Fermi Energy of Sodium

Sodium has n=2.65×1028 m3n = 2.65 \times 10^{28}\ \mathrm{m}^{-3} conduction electrons (one per atom, BCC structure).

kF=(3π2n)1/3=(3π2×2.65×1028)1/3=(7.85×1029)1/3=9.23×109 m1k_F = (3\pi^2 n)^{1/3} = (3\pi^2 \times 2.65 \times 10^{28})^{1/3} = (7.85 \times 10^{29})^{1/3} = 9.23 \times 10^9\ \mathrm{m}^{-1}

εF=2kF22me=(1.055×1034)2×(9.23×109)22×9.11×1031=9.48×10581.82×1030=5.21×1019 J=3.25 eV\varepsilon_F = \frac{\hbar^2 k_F^2}{2m_e} = \frac{(1.055 \times 10^{-34})^2 \times (9.23 \times 10^9)^2}{2 \times 9.11 \times 10^{-31}} = \frac{9.48 \times 10^{-58}}{1.82 \times 10^{-30}} = 5.21 \times 10^{-19}\ \mathrm{J} = 3.25\ \mathrm{eV}

TF=εFkB=5.21×10191.381×1023=3.77×104 KT_F = \frac{\varepsilon_F}{k_B} = \frac{5.21 \times 10^{-19}}{1.381 \times 10^{-23}} = 3.77 \times 10^4\ \mathrm{K}

vF=kFme=1.055×1034×9.23×1099.11×1031=1.07×106 m/sv_F = \frac{\hbar k_F}{m_e} = \frac{1.055 \times 10^{-34} \times 9.23 \times 10^9}{9.11 \times 10^{-31}} = 1.07 \times 10^6\ \mathrm{m}/s

The electronic specific heat coefficient: γ=π22nkB2εF=π2×2.65×1028×(1.381×1023)22×5.21×1019=1.38×103 J/(m3K2)\gamma = \frac{\pi^2}{2}\frac{nk_B^2}{\varepsilon_F} = \frac{\pi^2 \times 2.65 \times 10^{28} \times (1.381 \times 10^{-23})^2}{2 \times 5.21 \times 10^{-19}} = 1.38 \times 10^3\ \mathrm{J}/(m^3\cdot K^2)

The tight-binding model starts from isolated atomic orbitals and treats the overlap between Neighbours as a perturbation. For a 1D chain with lattice constant aa and a single ss-orbital Of energy ε0\varepsilon_0:

ψk(r)=1Nneiknaϕ(rna)\psi_k(r) = \frac{1}{\sqrt{N}}\sum_n e^{ikna}\,\phi(r - na)

Where ϕ(rna)\phi(r - na) is the atomic orbital centred at site nn.

Dispersion relation (nearest-neighbour approximation):

ε(k)=ε02tcos(ka)\varepsilon(k) = \varepsilon_0 - 2t\cos(ka)

Where t=ϕ(rna)H^ϕ(r(n+1)a)drt = -\int \phi^*(r - na)\,\hat{H}\,\phi(r - (n+1)a)\,dr is the hopping integral (t>0t > 0 for typical ss-orbitals).

Key features:

  • Band width: W=4tW = 4t.
  • Minimum at k=0k = 0: εmin=ε02t\varepsilon_{\min} = \varepsilon_0 - 2t.
  • Maximum at k=±π/ak = \pm\pi/a: εmax=ε0+2t\varepsilon_{\max} = \varepsilon_0 + 2t.
  • Effective mass at band bottom: m=2/(2ta2)m^* = \hbar^2/(2ta^2).

Extension to 3D: For a simple cubic lattice with nearest-neighbour hopping:

ε(k)=ε02t(coskxa+coskya+coskza)\varepsilon(\mathbf{k}) = \varepsilon_0 - 2t(\cos k_x a + \cos k_y a + \cos k_z a)

The band width is W=12tW = 12t and the density of states develops a van Hove singularity at ε=ε0\varepsilon = \varepsilon_0.

Worked Example: Tight-Binding Band Structure of Graphene

Graphene has a honeycomb lattice with two carbon atoms per unit cell. Using pzp_z orbitals with Nearest-neighbour hopping t2.8t \approx 2.8 eV, the tight-binding Hamiltonian gives:

ε±(k)=±t1+eika1+eika2\varepsilon_{\pm}(\mathbf{k}) = \pm t\left\lvert 1 + e^{i\mathbf{k}\cdot\mathbf{a}_1} + e^{i\mathbf{k}\cdot\mathbf{a}_2}\right\rvert

Where a1\mathbf{a}_1 and a2\mathbf{a}_2 are the primitive vectors of the hexagonal lattice.

The two bands touch at the Dirac points K\mathbf{K} and K\mathbf{K}' in the Brillouin zone. Near these points, expanding to linear order:

ε(q)=±vFq\varepsilon(\mathbf{q}) = \pm \hbar v_F \lvert\mathbf{q}\rvert

Where vF=32ta106v_F = \frac{\sqrt{3}}{2}\frac{ta}{\hbar} \approx 10^6 m/s and q=kK\mathbf{q} = \mathbf{k} - \mathbf{K}.

This linear (Dirac-like) dispersion means graphene has zero effective mass and a density of states g(ε)εg(\varepsilon) \propto \lvert\varepsilon\rvert (linear in energy), unlike the ε\sqrt{\varepsilon} Dependence of a parabolic band.

Near a band extremum at k0\mathbf{k}_0The energy can be expanded:

ε(k)=ε0+22ij(m1)ij(kik0,i)(kjk0,j)\varepsilon(\mathbf{k}) = \varepsilon_0 + \frac{\hbar^2}{2}\sum_{ij}(m^{-1})_{ij}(k_i - k_{0,i})(k_j - k_{0,j})

The effective mass tensor (m1)ij=122εkikj(m^{-1})_{ij} = \frac{1}{\hbar^2}\frac{\partial^2 \varepsilon}{\partial k_i \partial k_j} Determines the response to external fields. For isotropic bands, m=2/(d2ε/dk2)m^* = \hbar^2/(d^2\varepsilon/dk^2).

A large effective mass means a flat band (small group velocity). A small effective mass means a Steep band (high mobility).

The effective mass can be negative near a band maximum (holes). Cyclotron resonance experiments Measure mm^* directly: the resonance frequency is ωc=eB/m\omega_c = eB/m^*.

flowchart TD
A[5_Electronic Band Structure] --> B[Key Concepts]
A --> C[Core Principles]
A --> D[Practical Applications]
B --> E[Fundamental definitions]
C --> F[Design patterns]
D --> G[Real-world usage]

Band structure is like a highway system for electrons. In free space, electrons can have any energy, but a crystal lattice acts like a periodic toll booth that blocks certain energy ranges entirely. These forbidden zones are band gaps, and they are why some materials conduct while others insulate. The Brillouin zone is the set of unique momentum states an electron can occupy in the lattice, analogous to how a repeating wallpaper pattern has a single tile that encodes the whole design. Effective mass captures how strongly the lattice potential slows or accelerates an electron, allowing us to treat it as a free particle with modified inertia.

Mistake 1: Assuming that the free electron model applies to all metals. The free electron model works well for simple metals (e.g., alkali metals) but fails for transition metals and other materials with complex band structures. Do not assume that the free electron model is universally applicable.

Mistake 2: Confusing the Fermi energy with the Fermi level. The Fermi energy εF\varepsilon_F is the energy of the highest occupied state at T=0T = 0, while the Fermi level EFE_F is the chemical potential at finite temperature. In metals, they are approximately equal, but in semiconductors they can differ significantly. Do not confuse the two concepts.

Mistake 3: Forgetting that band gaps arise from the periodic potential. Band gaps open at Brillouin zone boundaries due to the periodic potential of the lattice. Without a periodic potential (free electron model), there are no band gaps. Do not assume that band gaps exist in all materials; they require a periodic potential.

Mistake 4: Assuming that DFT gives exact band gaps. DFT with LDA or GGA functionals underestimates band gaps by 30—50%. Hybrid functionals or GW calculations are needed for accurate band gaps. Do not assume that DFT band gaps are quantitatively accurate.

Mistake 5: Confusing the effective mass with the rest mass. The effective mass mm^* describes how an electron responds to external fields in a periodic potential. It can be larger or smaller than the rest mass mem_e, and can even be negative. Do not assume that the effective mass equals the rest mass.