Physics / Condensed Matter Solid State Physics I 100% Free Open Access
Chapter 6 • Theory & Derivations

Electronic Band Structure & Bloch Electron Dynamics

Comprehensive energy band theory in crystals: periodic potential, Bloch theorem, Kronig-Penney solvable model, origin of energy band gaps, reduced and extended zone schemes, nearly free electron approximation, tight-binding LCAO method, semiclassical electron dynamics, group velocity, effective mass tensor, concept of positive holes, and Fermi surfaces.

§6.1 Periodic Crystal Potential & Bloch's Theorem

1. Electrons in a Periodic Crystal Potential

An electron in a perfect crystalline solid experiences a potential energy $V(\vec{r})$ that possesses the full discrete translational periodicity of the Bravais lattice:

$$V(\vec{r} + \vec{R}) = V(\vec{r}) \quad \forall \vec{R} = n_1 \vec{a}_1 + n_2 \vec{a}_2 + n_3 \vec{a}_3$$

The single-electron Schrödinger equation is:

$$\hat{H} \psi(\vec{r}) = \left[ - \frac{\hbar^2}{2m} \nabla^2 + V(\vec{r}) \right] \psi(\vec{r}) = E \psi(\vec{r})$$

2. Formal Statement & Proof of Bloch's Theorem

Because the Hamiltonian commutes with all lattice translation operators $\hat{T}_{\vec{R}}$ ($[\hat{H}, \hat{T}_{\vec{R}}] = 0$), simultaneous eigenstates of $\hat{H}$ and $\hat{T}_{\vec{R}}$ can be constructed. Bloch's Theorem states that the eigenstates of an electron in a periodic potential can be chosen in the form of a plane wave modulated by a periodic function $u_{n,\vec{k}}(\vec{r})$ having the periodicity of the lattice:

$$\psi_{n,\vec{k}}(\vec{r}) = e^{i \vec{k} \cdot \vec{r}} u_{n,\vec{k}}(\vec{r})$$

where $u_{n,\vec{k}}(\vec{r} + \vec{R}) = u_{n,\vec{k}}(\vec{r})$ for all lattice vectors $\vec{R}$. Equivalently, translating by any lattice vector $\vec{R}$ merely shifts the phase of the wavefunction:

$$\psi_{n,\vec{k}}(\vec{r} + \vec{R}) = e^{i \vec{k} \cdot \vec{R}} \psi_{n,\vec{k}}(\vec{r})$$

The vector $\vec{k}$ is the crystal wavevector, and the integer $n = 1, 2, 3, \dots$ is the band index.

3. Crystal Momentum vs. Physical Momentum

The quantity $\hbar \vec{k}$ is known as crystal momentum. It is not the eigenvalue of the physical momentum operator $-i\hbar\nabla$ (since momentum is not conserved in the presence of the lattice potential). Instead, crystal momentum is conserved in electron-photon, electron-phonon, and electron-electron scattering events modulo an arbitrary reciprocal lattice vector $\vec{G}$:

$$\vec{k}' = \vec{k} + \vec{q} + \vec{G}$$

§6.2 The Kronig-Penney Model & Origin of Energy Band Gaps

1. The Kronig-Penney 1D Solvable Model

R. de L. Kronig and W. G. Penney (1931) introduced an idealized one-dimensional periodic array of rectangular potential barriers of height $V_0$, barrier width $b$, and well width $a$ (lattice constant $d = a + b$):

$$V(x) = \begin{cases} 0, & 0 < x < a \\ V_0, & -b < x < 0 \end{cases}$$

In Region I ($0 < x < a$), $\psi_1(x) = A e^{i K x} + B e^{-i K x}$, where $K = \sqrt{2mE/\hbar^2}$.

In Region II ($-b < x < 0$), for $E < V_0$, $\psi_2(x) = C e^{Q x} + D e^{-Q x}$, where $Q = \sqrt{2m(V_0 - E)/\hbar^2}$.

2. The Dirac Delta-Function Barrier Limit

Taking the limit where the barriers become infinitely thin and tall ($b \to 0, V_0 \to \infty$) while their barrier area $P = \lim \frac{m V_0 b a}{\hbar^2}$ remains constant, matching wavefunctions and their derivatives across the boundary using Bloch's theorem yields the famous Kronig-Penney dispersion relation:

$$P \frac{\sin(K a)}{K a} + \cos(K a) = \cos(k a)$$

where $P$ is the dimensionless barrier strength (representing the strength of the periodic crystal potential).

3. Origin of Energy Bands and Band Gaps

Because the right-hand side is $\cos(k a)$, real solutions for crystal wavevector $k$ exist if and only if the left-hand side satisfies:

$$-1 \le P \frac{\sin(K a)}{K a} + \cos(K a) \le +1$$
  • Allowed Energy Bands: Values of $E = \frac{\hbar^2 K^2}{2m}$ where the condition is satisfied represent allowed energy bands.
  • Forbidden Band Gaps ($E_g$): Energy ranges where the magnitude of the left-hand side exceeds unity ($> +1$ or $< -1$). No traveling Bloch waves can propagate through the crystal at these energies; electron waves undergo destructive interference and total Bragg reflection.
  • Limits:
    • $P \to 0$ (Free Electrons): $\cos(Ka) = \cos(ka) \implies K = k \implies E = \frac{\hbar^2 k^2}{2m}$ (continuous parabolic spectrum).
    • $P \to \infty$ (Isolated Atoms): $\sin(Ka) = 0 \implies Ka = n\pi \implies E_n = \frac{n^2 \pi^2 \hbar^2}{2m a^2}$ (discrete bound atomic levels).

§6.3 Zone Schemes & Nearly Free Electron (NFE) Approximation

1. Representation of Energy Bands: Zone Schemes

Because energy eigenvalues satisfy $E_n(k + G) = E_n(k)$, the band structure can be represented in three equivalent representations:

  1. Extended Zone Scheme: Band $n$ is plotted in the $n$-th Brillouin zone: band 1 in $[-\pi/a, +\pi/a]$, band 2 in $[-2\pi/a, -\pi/a] \cup [+\pi/a, +2\pi/a]$, etc. Most closely resembles the free-electron parabola $E = \hbar^2 k^2/2m$.
  2. Reduced Zone Scheme: All energy bands are mapped back into the First Brillouin Zone ($-\pi/a \le k \le +\pi/a$) by translating by appropriate reciprocal lattice vectors $G = 2\pi n / a$. Standard representation used in modern solid-state physics.
  3. Periodic (Repeated) Zone Scheme: The First Brillouin Zone dispersion is repeated periodically across all $k$-space with period $2\pi/a$. Convenient for visualizing semiclassical electron trajectories.

2. The Nearly Free Electron (NFE) Approximation

When the periodic crystal potential $V(x)$ is weak compared to electron kinetic energies, it can be treated as a perturbation. Expanding the potential in a Fourier series over reciprocal lattice vectors $G$:

$$V(x) = \sum_{G} V_G e^{i G x} \quad (V_0 = 0)$$

Away from the Brillouin zone boundaries, standard non-degenerate perturbation theory shifts energies insignificantly. However, near the zone boundary $k = \pm G/2$, the two unperturbed plane-wave states $|k\rangle$ and $|k - G\rangle$ are degenerate ($E_0(k) \approx E_0(k - G)$).

Applying degenerate perturbation theory with trial state $\psi = c_1 |k\rangle + c_2 |k - G\rangle$ yields the $2 \times 2$ secular equation:

$$\begin{pmatrix} E_0(k) - E & V_G \\ V_G^* & E_0(k - G) - E \end{pmatrix} \begin{pmatrix} c_1 \\ c_2 \end{pmatrix} = 0$$

At the exact zone boundary $k = G/2$, $E_0(k) = E_0(k - G) = \frac{\hbar^2 (G/2)^2}{2m} = E_0$:

$$(E_0 - E)^2 - |V_G|^2 = 0 \implies E_{\pm} = E_0 \pm |V_G|$$

A band gap of magnitude $\Delta E_g = 2 |V_G|$ opens up at every Brillouin zone boundary. The standing wave solutions $\psi_+ \sim \cos(Gx/2)$ and $\psi_- \sim \sin(Gx/2)$ pile up electronic charge density either at the ion cores (lowering energy by $-|V_G|$) or midway between the ion cores (raising energy by $+|V_G|$).

§6.4 The Tight-Binding Approximation (LCAO for Crystals)

1. Physical Foundation of Tight-Binding

In contrast to the nearly free electron model, the Tight-Binding Approximation assumes that electrons are tightly bound to individual atomic cores, spending most of their time in localized atomic orbitals $\phi(\vec{r} - \vec{R}_n)$. When atoms are brought together into a crystal lattice, the overlap between adjacent atomic wavefunctions broadens discrete atomic energy levels into continuous energy bands.

2. Formulation of the Bloch Sum

To satisfy Bloch's theorem, the trial crystal wavefunction is formed as a coherent linear combination of atomic orbitals (LCAO) summed over all $N$ lattice sites $\vec{R}_n$:

$$\psi_k(\vec{r}) = \frac{1}{\sqrt{N}} \sum_{n} e^{i \vec{k}\cdot\vec{R}_n} \phi(\vec{r} - \vec{R}_n)$$

3. Derivation of 1D and 3D Tight-Binding Dispersion

Evaluating the expectation value of the crystal Hamiltonian $\hat{H} = \hat{H}_{\text{atom}} + \Delta U(\vec{r})$:

$$E(k) = \frac{\langle \psi_k | \hat{H} | \psi_k \rangle}{\langle \psi_k | \psi_k \rangle} \approx \frac{\sum_{n, m} e^{i \vec{k}\cdot(\vec{R}_m - \vec{R}_n)} \int \phi^*(\vec{r} - \vec{R}_n) \hat{H} \phi(\vec{r} - \vec{R}_m) d^3r}{N}$$

Retaining only on-site and nearest-neighbor matrix elements:

  • On-site atomic energy: $\varepsilon_0 = \int \phi^*(\vec{r}) \hat{H} \phi(\vec{r}) d^3r \approx - E_0 - \alpha$.
  • Nearest-neighbor transfer (hopping) integral: $t = - \int \phi^*(\vec{r}) \hat{H} \phi(\vec{r} - \vec{a}) d^3r > 0$.

For a one-dimensional chain with nearest-neighbor distance $a$:

$$E(k) = \varepsilon_0 - 2 t \cos(k a)$$

The total bandwidth of the tight-binding band is $W = E_{\text{max}} - E_{\text{min}} = (\varepsilon_0 + 2t) - (\varepsilon_0 - 2t) = 4t$.

For a 3D Simple Cubic lattice with nearest-neighbor spacing $a$:

$$E(\vec{k}) = \varepsilon_0 - 2 t [\cos(k_x a) + \cos(k_y a) + \cos(k_z a)] \quad (\text{Bandwidth } W = 12 t)$$

§6.5 Semiclassical Electron Dynamics, Effective Mass & Positive Holes

1. Semiclassical Equations of Motion for Bloch Electrons

An electron in an energy band $E_n(\vec{k})$ is described by a localized wave packet centered at position $\vec{r}$ and mean wavevector $\vec{k}$. Its dynamical evolution under external electromagnetic fields $\vec{E}$ and $\vec{B}$ is governed by the semiclassical equations of motion:

  1. Group Velocity:
    $$\vec{v}(\vec{k}) = \frac{1}{\hbar} \nabla_{\vec{k}} E(\vec{k})$$
  2. Rate of Change of Crystal Momentum:
    $$\hbar \frac{d\vec{k}}{dt} = \vec{F}_{\text{ext}} = - e [\vec{E} + \vec{v}(\vec{k}) \times \vec{B}]$$

2. Derivation of the Effective Mass Tensor $m^*$

Differentiating the group velocity with respect to time:

$$\frac{d\vec{v}}{dt} = \frac{1}{\hbar} \frac{d}{dt} \nabla_{\vec{k}} E(\vec{k}) = \frac{1}{\hbar} \sum_{j} \left( \frac{\partial^2 E}{\partial k_i \partial k_j} \right) \frac{dk_j}{dt}$$

Substituting $\hbar \frac{dk_j}{dt} = F_j$:

$$\frac{dv_i}{dt} = \sum_{j} \left( \frac{1}{\hbar^2} \frac{\partial^2 E}{\partial k_i \partial k_j} \right) F_j = \sum_{j} \left(\frac{1}{m^*}\right)_{ij} F_j$$

Comparing with Newton's second law $\vec{a} = (m^*)^{-1} \vec{F}$ defines the Effective Mass Tensor $(m^*)_{ij}$:

$$\left(\frac{1}{m^*}\right)_{ij} = \frac{1}{\hbar^2} \frac{\partial^2 E(\vec{k})}{\partial k_i \partial k_j}$$

In an isotropic 1D band: $m^* = \hbar^2 / \left(\frac{d^2E}{dk^2}\right)$.

  • Near the bottom of a band ($\frac{d^2E}{dk^2} > 0$): $m^* > 0$. The electron accelerates in the direction of the applied electric force like a free electron.
  • Near the top of a band ($\frac{d^2E}{dk^2} < 0$): $m^* < 0$. The electron accelerates in the direction opposite to the external force because of Bragg reflections from the lattice potential.

3. The Concept of Positive Holes

In a nearly filled band, the absence of an electron from a state with wavevector $\vec{k}_e$, energy $E_e$, and charge $-e$ is described as a quasiparticle called a hole:

  • Wavevector: $\vec{k}_h = - \vec{k}_e$
  • Energy: $E_h(\vec{k}_h) = - E_e(\vec{k}_e)$
  • Velocity: $\vec{v}_h = \vec{v}_e$
  • Charge: $q_h = + e$ (positive charge)
  • Effective Mass: $m_h^* = - m_e^* > 0$ (positive effective mass at the band top)

Treating valence band transport in terms of positive holes simplifies the quantum statistical mechanics of semiconductors and metals.

Solved Problem Example 6.1: Energy Band Gap Opening in the Nearly Free Electron Model

An electron moves in a 1D crystal with lattice spacing $a = 3.0 \text{ Å}$ under a weak periodic potential $V(x) = 2 V_1 \cos(2\pi x / a)$ with Fourier amplitude $V_1 = 0.75 \text{ eV}$. (a) Calculate the unperturbed free-electron kinetic energy $E_0$ at the First Brillouin Zone boundary $k = \pi / a$. (b) Determine the energy values $E_+$ and $E_-$ at the zone boundary and the magnitude of the band gap $\Delta E_g$.

Step 1: Calculate the Reciprocal Lattice Vector G
$$G = \frac{2\pi}{a} = \frac{2\pi}{3.0 \times 10^{-10}\text{ m}} \approx 2.094 \times 10^{10}\text{ m}^{-1} \implies k_{\text{edge}} = \frac{G}{2} = \frac{\pi}{a} \approx 1.047 \times 10^{10}\text{ m}^{-1}$$

Compute the First Brillouin Zone boundary wavevector.

Step 2: Calculate the Unperturbed Free-Electron Energy E_0
$$E_0 = \frac{\hbar^2 k_{\text{edge}}^2}{2m} = \frac{(1.0546 \times 10^{-34}\text{ J}\cdot\text{s})^2 (1.047 \times 10^{10}\text{ m}^{-1})^2}{2(9.109 \times 10^{-31}\text{ kg})} = \frac{1.2185 \times 10^{-47}}{1.8218 \times 10^{-30}}\text{ J} \approx 6.688 \times 10^{-19}\text{ J} \approx 4.174\text{ eV}$$

Compute unperturbed kinetic energy in electron-volts.

Step 3: Calculate the Band Gap Delta E_g
$$\Delta E_g = 2 |V_G| = 2 V_1 = 2 \times 0.75\text{ eV} = 1.50\text{ eV}$$

In degenerate perturbation theory, the band gap equals twice the Fourier component of the potential.

Step 4: Determine the Band Edge Energies E_+ and E_-
$$E_- = E_0 - V_1 = 4.174\text{ eV} - 0.75\text{ eV} = 3.424\text{ eV}, \quad E_+ = E_0 + V_1 = 4.174\text{ eV} + 0.75\text{ eV} = 4.924\text{ eV}$$

The lower band terminates at E_- and the upper band begins at E_+.

Final Answer & Physical Insight

E_0 = 4.17 \text{ eV}, \quad \Delta E_g = 1.50 \text{ eV}, \quad E_- = 3.42 \text{ eV}, \quad E_+ = 4.92 \text{ eV}

Solved Problem Example 6.2: Effective Mass and Group Velocity in a 1D Tight-Binding Band

The energy dispersion of a 1D tight-binding conduction band is given by $E(k) = - E_0 - 2t \cos(ka)$, where $E_0 = 1.20 \text{ eV}$, hopping integral $t = 0.85 \text{ eV}$, and lattice spacing $a = 2.50 \text{ Å}$. (a) Calculate the total bandwidth $W$. (b) Derive the expression for effective mass $m^*(k)$ and evaluate its numerical value in units of free electron mass $m_0$ at the band bottom ($k = 0$) and band top ($k = \pi/a$). (c) Find the wavevector $k$ at which the electron group velocity $v_g$ reaches its maximum value.

Step 1: Calculate the Total Bandwidth W
$$W = E_{\text{max}} - E_{\text{min}} = 4t = 4 \times 0.85\text{ eV} = 3.40\text{ eV}$$

In a 1D tight-binding cosine band, the bandwidth equals 4 times the transfer integral t.

Step 2: Derive the Effective Mass Formula
$$\frac{dE}{dk} = 2 t a \sin(k a) \implies \frac{d^2E}{dk^2} = 2 t a^2 \cos(k a) \implies m^*(k) = \frac{\hbar^2}{2 t a^2 \cos(k a)}$$

Differentiate E(k) twice with respect to wavevector k.

Step 3: Evaluate m* at k = 0 (Band Bottom)
$$m^*(0) = \frac{\hbar^2}{2 t a^2} = \frac{(1.0546 \times 10^{-34}\text{ J}\cdot\text{s})^2}{2 \times (0.85 \times 1.6022 \times 10^{-19}\text{ J}) \times (2.50 \times 10^{-10}\text{ m})^2} = \frac{1.1122 \times 10^{-68}}{2.7237 \times 10^{-19} \times 6.25 \times 10^{-20}} = \frac{1.1122 \times 10^{-68}}{1.7023 \times 10^{-38}}\text{ kg} \approx 6.533 \times 10^{-31}\text{ kg} \implies \frac{m^*(0)}{m_0} = \frac{6.533 \times 10^{-31}}{9.109 \times 10^{-31}} \approx 0.717$$

Substitute numerical values to find the effective mass at k = 0.

Step 4: Evaluate m* at k = pi / a (Band Top)
$$m^*(\pi/a) = \frac{\hbar^2}{2 t a^2 \cos(\pi)} = - m^*(0) = - 0.717 m_0$$

At the zone boundary, cos(pi) = -1, yielding a negative effective mass corresponding to a positive hole mass m_h* = +0.717 m_0.

Step 5: Determine Maximum Group Velocity
$$v_g(k) = \frac{1}{\hbar}\frac{dE}{dk} = \frac{2ta}{\hbar}\sin(ka) \implies \left. v_g \right|_{\text{max}} \text{ occurs when } \sin(ka) = 1 \implies ka = \frac{\pi}{2} \implies k = \frac{\pi}{2a}$$

The maximum group velocity occurs at the inflection point k = pi / (2a) where effective mass diverges to infinity.

Final Answer & Physical Insight

W = 3.40 \text{ eV}, \quad m^(0) = +0.717 m_0, \quad m^(\pi/a) = -0.717 m_0, \quad k_{\text{max}} = \pi / (2a)

Solved Problem Example 6.3: Period of Bloch Oscillations in an External Electric Field

An electron in a crystal with lattice constant $a = 3.50 \text{ Å}$ is subjected to a uniform external electric field $E_x = 1.50 \times 10^5 \text{ V/m}$. Assuming negligible scattering (mean free path exceeds the oscillation length): (a) Calculate the rate of change of crystal wavevector $dk/dt$. (b) Determine the time period $T_B$ and cyclical frequency $\nu_B$ of the resulting Bloch oscillations.

Step 1: Calculate the Rate of Change of Wavevector dk / dt
$$\hbar \frac{dk}{dt} = - e E_x \implies \left|\frac{dk}{dt}\right| = \frac{e E_x}{\hbar} = \frac{(1.6022 \times 10^{-19}\text{ C})(1.50 \times 10^5\text{ V/m})}{1.0546 \times 10^{-34}\text{ J}\cdot\text{s}} \approx 2.279 \times 10^{20}\text{ m}^{-1}\cdot\text{s}^{-1}$$

Apply the semiclassical equation of motion in an electric field.

Step 2: Determine the Wavevector Traversed in One Full Oscillation
$$\Delta k = \frac{2\pi}{a} = \frac{2\pi}{3.50 \times 10^{-10}\text{ m}} \approx 1.795 \times 10^{10}\text{ m}^{-1}$$

One full cycle corresponds to traversing the entire First Brillouin Zone of width 2 pi / a.

Step 3: Calculate the Bloch Oscillation Period T_B
$$T_B = \frac{\Delta k}{|dk/dt|} = \frac{2\pi \hbar}{e E_x a} = \frac{h}{e E_x a} = \frac{6.626 \times 10^{-34}\text{ J}\cdot\text{s}}{(1.6022 \times 10^{-19}\text{ C})(1.50 \times 10^5\text{ V/m})(3.50 \times 10^{-10}\text{ m})} = \frac{6.626 \times 10^{-34}}{8.4116 \times 10^{-24}}\text{ s} \approx 7.877 \times 10^{-11}\text{ s} = 78.8\text{ ps}$$

Evaluate the Bloch period T_B = h / (e E a).

Step 4: Calculate the Bloch Oscillation Frequency nu_B
$$\nu_B = \frac{1}{T_B} = \frac{1}{7.877 \times 10^{-11}\text{ s}} \approx 1.269 \times 10^{10}\text{ Hz} = 12.7\text{ GHz}$$

Inverted period gives the Bloch oscillation frequency in the microwave / terahertz regime.

Final Answer & Physical Insight

T_B = 78.8 \text{ ps}, \quad \nu_B = 12.7 \text{ GHz}

EXAM SUCCESS WORKSHOP

Solved University Examination Problems

Step-by-step mathematical solutions to classic university honors examination questions.