Chapter 1

Berry Phase in Solids

Geometric Phase in Quantum Mechanics

When a quantum system is transported adiabatically around a closed loop in parameter space, its state acquires a phase factor beyond the familiar dynamical phase. This additional contribution, discovered by Michael Berry in 1984, is purely geometric in origin -- it depends only on the path traversed in parameter space, not on the rate of traversal. In the context of Bloch bands, the Berry phase underlies a remarkable range of phenomena from electric polarization to the quantum Hall effect.

For a Hamiltonian $H(\mathbf{R})$ with slowly varying parameters $\mathbf{R}(t)$, the total phase splits into a dynamical part and a geometric part:

$$\Phi_n = \underbrace{-\frac{1}{\hbar}\int_0^T E_n(\mathbf{R}(t))\,dt}_{\text{dynamical phase}} + \underbrace{\gamma_n}_{\text{Berry phase}}$$

The Berry phase $\gamma_n$ is gauge-invariant and connects band topology, transport, and electromagnetic response in a unified framework.

The Berry Phase

The Berry phase acquired by the $n$-th eigenstate upon adiabatic transport around a closed loop $\mathcal{C}$ in parameter space is:

Berry Phase

$$\gamma_n = \oint_\mathcal{C} \langle n(\mathbf{R})|\nabla_\mathbf{R}|n(\mathbf{R})\rangle \cdot d\mathbf{R}$$

This line integral is gauge-invariant modulo $2\pi$: under a gauge transformation $|n\rangle \to e^{i\chi(\mathbf{R})}|n\rangle$, the Berry phase changes by $\oint \nabla_\mathbf{R}\chi \cdot d\mathbf{R}$, which is an integer multiple of $2\pi$ for a well-defined gauge around the loop.

For a two-level system $H = \mathbf{R}\cdot\boldsymbol{\sigma}$, the Berry phase of the ground state equals half the solid angle subtended by the loop:

$$\gamma_- = -\frac{1}{2}\,\Omega(\mathcal{C})$$

This exact result for any spin-$\frac{1}{2}$ system gives geometric intuition: the Berry phase measures the solid angle enclosed by the parameter trajectory on the Bloch sphere.

Berry Connection and Curvature

The Berry phase formalism naturally introduces a gauge structure analogous to electromagnetism. The Berry connection (or Berry potential) plays the role of a vector potential in parameter space:

Berry Connection

$$\mathcal{A}_n(\mathbf{k}) = i\langle n(\mathbf{k})|\nabla_\mathbf{k}|n(\mathbf{k})\rangle$$

In the Bloch band context, $\mathbf{k}$ is the crystal momentum. The Berry connection is gauge-dependent: under $|n\rangle \to e^{i\chi(\mathbf{k})}|n\rangle$, we have $\mathcal{A}_n \to \mathcal{A}_n - \nabla_\mathbf{k}\chi$, exactly like the electromagnetic vector potential.

Berry Curvature

$$\boldsymbol{\Omega}_n(\mathbf{k}) = \nabla_\mathbf{k} \times \mathcal{A}_n(\mathbf{k})$$

In two dimensions, the Berry curvature has a single component:

$$\Omega_n(k_x, k_y) = \frac{\partial \mathcal{A}_{n,y}}{\partial k_x} - \frac{\partial \mathcal{A}_{n,x}}{\partial k_y}$$

Unlike the Berry connection, the Berry curvature is gauge-invariant and therefore physically observable. It acts as a “magnetic field” in momentum space, deflecting wave packets through the anomalous velocity.

A practical Kubo-like formula expresses Berry curvature using all bands:

$$\Omega_n(\mathbf{k}) = -2\,\text{Im}\sum_{m \neq n} \frac{\langle n|\partial_{k_x} H|m\rangle\langle m|\partial_{k_y} H|n\rangle}{(E_m - E_n)^2}$$

Berry curvature peaks near band degeneracies and acts as a momentum-space “magnetic monopole” at points of exact degeneracy. Symmetry constraints: $\Omega_n(-\mathbf{k}) = -\Omega_n(\mathbf{k})$ under TRS,$\Omega_n(-\mathbf{k}) = \Omega_n(\mathbf{k})$ under inversion; when both hold,$\Omega_n = 0$ everywhere.

Chern Number

The integral of the Berry curvature over the entire Brillouin zone yields a topological invariant -- the first Chern number:

Chern Number (TKNN Invariant)

$$C_n = \frac{1}{2\pi}\int_\text{BZ} \Omega_n(\mathbf{k})\, d^2k$$

The Chern number is always an integer. This quantization follows from the fact that the Brillouin zone is a closed 2-torus (due to the periodicity $\mathbf{k} \sim \mathbf{k} + \mathbf{G}$), and the integral of the curvature of a U(1) bundle over a closed manifold is quantized: $C_n \in \mathbb{Z}$.

The Chern number determines the Hall conductivity via the TKNN formula (Thouless, Kohmoto, Nightingale, den Nijs, 1982): $\sigma_{xy} = (e^2/h)\sum_n C_n$, summed over filled bands. A non-zero Chern number implies:

Broken time-reversal symmetry (since $\int \Omega\,d^2k = 0$ under TRS)
Gapless chiral edge states at the boundary (bulk-boundary correspondence)
Quantized Hall conductance robust to disorder and interactions
Impossibility of constructing smooth, periodic Bloch functions across the BZ

Berry Phase in Two-Level Systems

The simplest non-trivial example of Berry phase arises in a generic two-level system, which also serves as the building block for understanding topological bands. Consider the Hamiltonian:

$$H(\mathbf{k}) = \mathbf{d}(\mathbf{k})\cdot\boldsymbol{\sigma} = d_x\sigma_x + d_y\sigma_y + d_z\sigma_z$$

where $\boldsymbol{\sigma} = (\sigma_x, \sigma_y, \sigma_z)$ are the Pauli matrices and $\mathbf{d}(\mathbf{k})$ is a three-component vector field over the BZ. The eigenvalues are $E_\pm = \pm|\mathbf{d}|$ and the Berry curvature of the lower band is:

$$\Omega_-(\mathbf{k}) = -\frac{1}{2}\,\hat{\mathbf{d}}\cdot\left(\frac{\partial\hat{\mathbf{d}}}{\partial k_x}\times\frac{\partial\hat{\mathbf{d}}}{\partial k_y}\right)$$

where $\hat{\mathbf{d}} = \mathbf{d}/|\mathbf{d}|$. The Chern number counts the number of times the map $\hat{\mathbf{d}}: \text{BZ} \to S^2$ wraps the unit sphere -- a winding number in $\pi_2(S^2) = \mathbb{Z}$.

Example: Haldane Model

The Haldane model on the honeycomb lattice has $\mathbf{d}(\mathbf{k})$ with components determined by nearest-neighbor hopping $t_1$, next-nearest-neighbor hopping $t_2 e^{\pm i\phi}$ (which breaks TRS), and a sublattice mass $M$ (which breaks inversion). The phase diagram features:

$C = 0$: trivial insulator when $|M/t_2| > 3\sqrt{3}\sin\phi$
$C = \pm 1$: Chern insulator with quantized Hall effect
Phase boundaries occur when the gap closes at K or K' points

Anomalous Velocity

The Berry curvature modifies the semiclassical equations of motion for Bloch electrons. The velocity of a wave packet in band $n$ acquires an anomalous contribution perpendicular to the applied force:

Semiclassical Equations of Motion

$$\dot{\mathbf{r}} = \frac{1}{\hbar}\nabla_\mathbf{k}\epsilon_n(\mathbf{k}) - \dot{\mathbf{k}}\times\boldsymbol{\Omega}_n(\mathbf{k})$$

$$\hbar\dot{\mathbf{k}} = -e\mathbf{E} - e\dot{\mathbf{r}}\times\mathbf{B}$$

The second term in the velocity equation is the anomalous velocity. Even without an external magnetic field ($\mathbf{B} = 0$), an electric field drives a transverse current via $\dot{\mathbf{r}}_\perp = (e\mathbf{E}/\hbar)\times\boldsymbol{\Omega}_n$. Summing over all filled states gives the intrinsic anomalous Hall conductivity.

For a fully filled band, summing the anomalous velocity over all occupied states yields the TKNN formula $\sigma_{xy} = C_n\,e^2/h$. For partially filled bands in ferromagnetic metals, the Berry curvature near the Fermi surface produces a non-quantized but intrinsic anomalous Hall conductivity.

Physical Applications

Electric Polarization (Modern Theory)

The bulk electric polarization of a crystalline insulator is a Berry phase of the occupied Bloch states. The polarization per unit cell along direction $\alpha$ is:

$$P_\alpha = \frac{e}{(2\pi)^d}\sum_n \int_\text{BZ} \mathcal{A}_{n,\alpha}(\mathbf{k})\,d^dk$$

This King-Smith--Vanderbilt formula resolves a long-standing puzzle: the polarization is only defined modulo a “quantum of polarization” $eR/V$ (where $R$ is a lattice vector), reflecting the gauge freedom of the Berry phase. Only changes in polarization (e.g., across a ferroelectric transition) are uniquely defined, which is precisely what experiments measure.

Orbital Magnetization

The orbital magnetization of Bloch electrons has a Berry-phase contribution involving $\text{Im}\langle \nabla_\mathbf{k} u_n| \times (H_\mathbf{k} - E_n)|\nabla_\mathbf{k} u_n\rangle$, which accounts for wave-packet self-rotation and is essential for magnetoelectric effects in topological insulators.

Connection to Gauge Theory

The mathematical structure of Berry phase is identical to that of a U(1) gauge theory. The correspondence is:

Berry Phase	Electromagnetism
$\mathcal{A}_n(\mathbf{k})$ (Berry connection)	$\mathbf{A}(\mathbf{r})$ (vector potential)
$\Omega_n(\mathbf{k})$ (Berry curvature)	$\mathbf{B}(\mathbf{r})$ (magnetic field)
Chern number $C_n$	Magnetic monopole charge
Gauge transformation $\|n\rangle \to e^{i\chi}\|n\rangle$	$\mathbf{A} \to \mathbf{A} + \nabla\chi$
Band degeneracy	Dirac monopole singularity

For multi-band systems, the Berry connection generalizes to a non-Abelian gauge field: $\mathcal{A}_{nm}(\mathbf{k}) = i\langle n|\nabla_\mathbf{k}|m\rangle$, and the relevant topological invariants involve Wilson loops and non-Abelian Berry phases. This non-Abelian structure is essential for understanding $\mathbb{Z}_2$ topological insulators and fragile topology.

Detailed Derivations

Derivation 1: Berry Phase from the Adiabatic Theorem

Consider a Hamiltonian $H(\mathbf{R}(t))$ depending on slowly varying parameters $\mathbf{R}(t)$. The instantaneous eigenstates satisfy $H(\mathbf{R})|n(\mathbf{R})\rangle = E_n(\mathbf{R})|n(\mathbf{R})\rangle$. The adiabatic theorem states that if the system starts in $|n(\mathbf{R}(0))\rangle$, it remains in the instantaneous eigenstate up to a phase factor.

Step 1: Write the general time-dependent state as an expansion in instantaneous eigenstates:

$$|\Psi(t)\rangle = \sum_n c_n(t)\,e^{-\frac{i}{\hbar}\int_0^t E_n(t')\,dt'}\,|n(\mathbf{R}(t))\rangle$$

Step 2: Substitute into the Schrodinger equation $i\hbar\partial_t|\Psi\rangle = H|\Psi\rangle$. After using $H|n\rangle = E_n|n\rangle$, the dynamical phase exponentials cancel, leaving:

$$\dot{c}_n(t) = -c_n(t)\,\langle n(\mathbf{R})|\dot{n}(\mathbf{R})\rangle - \sum_{m \neq n} c_m(t)\,e^{-\frac{i}{\hbar}\int_0^t (E_m - E_n)\,dt'}\,\langle n|\dot{m}\rangle$$

where $|\dot{n}\rangle = \frac{d}{dt}|n(\mathbf{R}(t))\rangle$.

Step 3: In the adiabatic limit, the oscillating exponentials suppress the $m \neq n$ terms (the non-adiabatic transitions). Keeping only the diagonal term:

$$\dot{c}_n = -c_n\,\langle n|\dot{n}\rangle$$

Step 4: Integrate directly to obtain:

$$c_n(T) = c_n(0)\,\exp\!\left(-\int_0^T \langle n(\mathbf{R}(t))|\frac{d}{dt}|n(\mathbf{R}(t))\rangle\,dt\right)$$

Step 5: Use the chain rule $\frac{d}{dt}|n(\mathbf{R})\rangle = \nabla_\mathbf{R}|n\rangle \cdot \dot{\mathbf{R}}$ to rewrite the integral as a line integral in parameter space:

$$\gamma_n = i\oint_\mathcal{C} \langle n(\mathbf{R})|\nabla_\mathbf{R}|n(\mathbf{R})\rangle \cdot d\mathbf{R}$$

The factor of $i$ ensures $\gamma_n$ is real: since $\langle n|n\rangle = 1$, differentiating gives $\langle n|\nabla n\rangle + \langle \nabla n|n\rangle = 0$, so $\langle n|\nabla n\rangle$ is purely imaginary, and $\gamma_n = i \times (\text{imaginary}) = \text{real}$.

Derivation 2: Berry Curvature for a Two-Level System

Consider the generic two-level Hamiltonian $H = \mathbf{d}\cdot\boldsymbol{\sigma}$ with $\mathbf{d} = (d_x, d_y, d_z)$. We derive the Berry curvature of the lower band in closed form.

Step 1: Parametrize $\hat{\mathbf{d}} = \mathbf{d}/|\mathbf{d}|$ using spherical angles $(\theta, \phi)$:

$$\hat{\mathbf{d}} = (\sin\theta\cos\phi,\;\sin\theta\sin\phi,\;\cos\theta)$$

Step 2: The ground state (eigenvalue $-|\mathbf{d}|$) of $\hat{\mathbf{d}}\cdot\boldsymbol{\sigma}$ is:

$$|-\rangle = \begin{pmatrix} \sin(\theta/2) \\ -\cos(\theta/2)\,e^{i\phi} \end{pmatrix}$$

Step 3: Compute the Berry connection components. First, $\mathcal{A}_\theta$:

$$\mathcal{A}_\theta = i\langle -|\partial_\theta|-\rangle = i\left[\cos(\theta/2)\cdot\frac{1}{2} \cdot (-\cos(\theta/2)\,e^{i\phi}) \cdot (-e^{-i\phi}) + \sin(\theta/2)\cdot\frac{1}{2}\cdot\sin(\theta/2)\right] \cdot \frac{1}{i}$$

Working through the algebra (the two terms cancel), we find $\mathcal{A}_\theta = 0$. For $\mathcal{A}_\phi$:

$$\mathcal{A}_\phi = i\langle -|\partial_\phi|-\rangle = i\left[\sin(\theta/2) \cdot 0 + (-\cos(\theta/2)\,e^{-i\phi})\cdot(-\cos(\theta/2)\,i\,e^{i\phi})\right]$$

$$= i \cdot i\cos^2(\theta/2) = -\cos^2(\theta/2)$$

Step 4: The Berry curvature on the unit sphere is:

$$\Omega_{\theta\phi} = \partial_\theta \mathcal{A}_\phi - \partial_\phi \mathcal{A}_\theta = \partial_\theta\left[-\cos^2(\theta/2)\right] = \frac{1}{2}\sin\theta$$

Step 5: To express this in $(k_x, k_y)$ coordinates, use the chain rule and the Jacobian of the map $\hat{\mathbf{d}}(\mathbf{k}): \text{BZ} \to S^2$:

$$\Omega_{k_x k_y} = \Omega_{\theta\phi}\left(\frac{\partial\theta}{\partial k_x}\frac{\partial\phi}{\partial k_y} - \frac{\partial\theta}{\partial k_y}\frac{\partial\phi}{\partial k_x}\right)$$

The term in parentheses is the Jacobian of the $(\theta,\phi)$ parametrization, which combines with $\frac{1}{2}\sin\theta$ to give the compact vector identity:

$$\Omega_-(\mathbf{k}) = -\frac{1}{2}\,\hat{\mathbf{d}}\cdot\left(\frac{\partial\hat{\mathbf{d}}}{\partial k_x}\times\frac{\partial\hat{\mathbf{d}}}{\partial k_y}\right)$$

Integrating over the BZ gives $C = \frac{1}{2\pi}\int \Omega_-\,d^2k$, which counts the number of times $\hat{\mathbf{d}}$ wraps the unit sphere.

Derivation 3: Kubo Formula for Berry Curvature

We derive the sum-over-states expression for Berry curvature that avoids taking derivatives of eigenstates.

Step 1: Start from $H|n\rangle = E_n|n\rangle$. Differentiate both sides with respect to $k_\mu$:

$$(\partial_{k_\mu} H)|n\rangle + H|\partial_{k_\mu} n\rangle = (\partial_{k_\mu} E_n)|n\rangle + E_n|\partial_{k_\mu} n\rangle$$

Step 2: Project onto $\langle m|$ with $m \neq n$. Using $\langle m|n\rangle = 0$ and $\langle m|H = E_m\langle m|$:

$$\langle m|\partial_{k_\mu} H|n\rangle + E_m\langle m|\partial_{k_\mu} n\rangle = E_n\langle m|\partial_{k_\mu} n\rangle$$

Step 3: Solve for the off-diagonal matrix element of the derivative:

$$\langle m|\partial_{k_\mu}|n\rangle = \frac{\langle m|\partial_{k_\mu} H|n\rangle}{E_n - E_m}, \quad m \neq n$$

Step 4: The Berry curvature involves $\langle \partial_{k_x} n|\partial_{k_y} n\rangle - (x \leftrightarrow y)$. Insert a complete set $\sum_m |m\rangle\langle m|$:

$$\Omega_n = -2\,\text{Im}\sum_{m\neq n}\langle \partial_{k_x} n|m\rangle\langle m|\partial_{k_y} n\rangle$$

Step 5: Substitute the result from Step 3, noting $\langle \partial_{k_x} n|m\rangle = \langle m|\partial_{k_x} n\rangle^*$:

$$\Omega_n(\mathbf{k}) = -2\,\text{Im}\sum_{m \neq n} \frac{\langle n|\partial_{k_x} H|m\rangle\langle m|\partial_{k_y} H|n\rangle}{(E_m - E_n)^2}$$

This Kubo-like formula shows that the Berry curvature is large near degeneracies (where $E_m - E_n$ is small) and involves all bands, not just the band of interest.

Real-World Applications of Berry Phase

1. Precision Metrology and Resistance Standards

The quantized Hall conductance $\sigma_{xy} = C e^2/h$, directly determined by the Chern number (an integral of Berry curvature), underpins the modern definition of the ohm. The von Klitzing constant $R_K = h/e^2 = 25\,812.807\,\Omega$ provides an exact resistance standard since the 2019 SI redefinition.

2. Anomalous Hall Effect in Magnetic Materials

The intrinsic anomalous Hall effect in ferromagnets such as SrRuO$_3$ and Fe arises from the integral of Berry curvature over occupied states at the Fermi level. This is exploited in spintronic devices for magnetic field sensing and spin-orbit torque memories.

3. Modern Theory of Ferroelectricity

The King-Smith-Vanderbilt formula expresses the bulk electric polarization as a Berry phase of occupied Bloch states. This is implemented in all modern density-functional theory codes (VASP, Quantum ESPRESSO, ABINIT) and is essential for computational prediction of ferroelectric and piezoelectric properties of materials used in sensors, actuators, and memory devices.

4. Valley Hall Effect and Valleytronics

In transition metal dichalcogenides (MoS$_2$, WSe$_2$), the Berry curvature has opposite sign at the K and K' valleys. Circularly polarized light selectively excites carriers in one valley, producing a transverse valley current. This is the basis of valleytronics, a proposed paradigm for information processing.

5. Topological Photonics and Phononics

Berry curvature concepts have been extended to photonic crystals and phononic metamaterials. Topological photonic waveguides guide light around sharp corners without backscattering, with applications in integrated photonic circuits and robust optical delay lines. Acoustic topological insulators similarly enable unidirectional sound propagation.

Historical Context

The geometric phase was first identified by S. Pancharatnam in 1956 in the context of polarization optics, where he showed that light acquires a phase when its polarization state is cycled around a closed path on the Poincare sphere. This work was largely overlooked until Michael V. Berry independently discovered the quantum mechanical geometric phase in 1984 in his landmark paper "Quantal Phase Factors Accompanying Adiabatic Changes." Berry demonstrated that the phase is a general consequence of adiabatic transport and connected it to the geometry of parameter space.

Shortly after Berry's work, Barry Simon (1983) recognized the Berry phase as the holonomy of a U(1) fiber bundle, placing it within the mathematical framework of differential geometry and gauge theory. Frank Wilczek and Alfred Shapere extended the concept to non-Abelian gauge structures in 1984. Meanwhile, Thouless, Kohmoto, Nightingale, and den Nijs (TKNN, 1982) had already -- before Berry's paper -- shown that the quantized Hall conductance is a topological invariant (Chern number), which was later understood as the integral of Berry curvature. The TKNN paper thus implicitly contained the Berry phase in the context of Bloch bands.

The modern theory of electric polarization was developed by R.D. King-Smith and David Vanderbilt in 1993, expressing the bulk polarization as a Berry phase. Raffaele Resta and Nicola Marzari further developed the formalism, connecting it to Wannier functions and orbital magnetization. These ideas transformed computational materials science by providing a rigorous framework for calculating properties that depend on the phase of the wavefunction, not just its modulus.

Python Simulation: Berry Curvature of the Haldane Model

Berry curvature over the Brillouin zone for the Haldane model on the honeycomb lattice. The Chern number is computed by numerical integration using the lattice gauge method of Fukui, Hatsugai, and Suzuki.

Berry Curvature & Chern Number (Haldane Model)

Python

Computes Berry curvature over the BZ using the lattice gauge method and band structure along high-symmetry path

script.py177 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

# ──────────────────────────────────────────────────────────
# Berry Curvature of a Two-Band (Haldane-like) Model
# Compute Berry curvature over the Brillouin zone and
# integrate to obtain the Chern number
# ──────────────────────────────────────────────────────────

# Lattice vectors for honeycomb lattice
a1 = np.array([1.0, 0.0])
a2 = np.array([0.5, np.sqrt(3)/2])

# Reciprocal lattice vectors
b1 = 2*np.pi * np.array([1.0, -1.0/np.sqrt(3)])
b2 = 2*np.pi * np.array([0.0, 2.0/np.sqrt(3)])

# Nearest-neighbor vectors (honeycomb)
delta1 = np.array([0.0, 1.0/np.sqrt(3)])
delta2 = np.array([0.5, -0.5/np.sqrt(3)])
delta3 = np.array([-0.5, -0.5/np.sqrt(3)])

# Next-nearest-neighbor vectors
nnn = [a1, a2, a2 - a1, -a1, -a2, a1 - a2]

# Haldane model parameters
t1 = 1.0       # NN hopping
t2 = 0.15      # NNN hopping
phi = np.pi/2  # Haldane phase (breaks TRS)
M = 0.0        # Sublattice mass (breaks inversion)

def haldane_hamiltonian(kx, ky):
    """Build the 2x2 Haldane Hamiltonian at momentum (kx, ky)."""
    k = np.array([kx, ky])

# Off-diagonal: nearest-neighbor hopping
    f_k = sum(np.exp(1j * np.dot(k, d)) for d in [delta1, delta2, delta3])
    h_x = t1 * np.real(f_k)
    h_y = t1 * np.imag(f_k)

# Diagonal: NNN hopping with complex phase + sublattice mass
    signs = [1, 1, -1, -1, -1, 1]  # chirality signs for NNN
    g_k = sum(s * np.cos(np.dot(k, d)) for s, d in zip(signs, nnn))
    h_z = M + 2*t2*np.sin(phi) * sum(s * np.sin(np.dot(k, d))
                                       for s, d in zip(signs, nnn))

# Constant energy shift (does not affect topology)
    h_0 = 2*t2*np.cos(phi) * sum(np.cos(np.dot(k, d)) for d in nnn)

H = np.array([
        [h_0 + h_z,   h_x - 1j*h_y],
        [h_x + 1j*h_y, h_0 - h_z]
    ])
    return H

# Discretize the Brillouin zone
Nk = 100
s_vals = np.linspace(0, 1, Nk, endpoint=False)
dk = 1.0 / Nk

# Compute Berry curvature via the lattice gauge method (Fukui-Hatsugai-Suzuki)
berry_curv = np.zeros((Nk, Nk))
kx_grid = np.zeros((Nk, Nk))
ky_grid = np.zeros((Nk, Nk))

def get_ground_state(kx, ky):
    H = haldane_hamiltonian(kx, ky)
    vals, vecs = np.linalg.eigh(H)
    return vecs[:, 0]  # lower band

for i in range(Nk):
    for j in range(Nk):
        k = s_vals[i] * b1 + s_vals[j] * b2
        kx_grid[i, j] = k[0]
        ky_grid[i, j] = k[1]

# Four corner states for the plaquette
        k00 = s_vals[i] * b1 + s_vals[j] * b2
        k10 = (s_vals[i] + dk) * b1 + s_vals[j] * b2
        k01 = s_vals[i] * b1 + (s_vals[j] + dk) * b2
        k11 = (s_vals[i] + dk) * b1 + (s_vals[j] + dk) * b2

u00 = get_ground_state(*k00)
        u10 = get_ground_state(*k10)
        u01 = get_ground_state(*k01)
        u11 = get_ground_state(*k11)

# Link variables (U(1) Berry connection on the lattice)
        U1 = np.vdot(u00, u10)
        U2 = np.vdot(u10, u11)
        U3 = np.vdot(u11, u01)
        U4 = np.vdot(u01, u00)

# Berry flux through plaquette
        F12 = np.log(U1 * U2 * U3 * U4)
        berry_curv[i, j] = -F12.imag

# Chern number = sum of Berry curvature / (2*pi)
chern = np.sum(berry_curv) / (2 * np.pi)

print(f"Haldane model parameters:")
print(f"  t1 = {t1}, t2 = {t2}, phi = {phi:.4f} (pi/2)")
print(f"  Sublattice mass M = {M}")
print(f"  Chern number C = {chern:.6f}  (expected: 1)")
print(f"  Quantized Hall conductance: sigma_xy = {chern:.4f} * e^2/h")
print()

# ─── Plot Berry curvature ───
fig, axes = plt.subplots(1, 2, figsize=(14, 6))
fig.patch.set_facecolor('#0a0a1a')

for ax in axes:
    ax.set_facecolor('#0a0a1a')
    for spine in ax.spines.values():
        spine.set_edgecolor('#334155')

# Left panel: Berry curvature heatmap
im = axes[0].pcolormesh(kx_grid, ky_grid, berry_curv,
                         cmap='RdBu_r', shading='auto')
cbar = plt.colorbar(im, ax=axes[0])
cbar.ax.yaxis.set_tick_params(color='#94a3b8')
cbar.ax.yaxis.label.set_color('white')
cbar.outline.set_edgecolor('#334155')
plt.setp(plt.getp(cbar.ax.axes, 'yticklabels'), color='#94a3b8')
cbar.set_label(r'$\Omega_n(\mathbf{k})$', color='white', fontsize=12)

axes[0].set_xlabel(r'$k_x$', color='#94a3b8', fontsize=12)
axes[0].set_ylabel(r'$k_y$', color='#94a3b8', fontsize=12)
axes[0].set_title('Berry Curvature (Lower Band)', color='white', fontsize=14)
axes[0].tick_params(colors='#94a3b8')
axes[0].set_aspect('equal')

# Right panel: Band structure along high-symmetry path
# Path: Gamma -> K -> M -> Gamma
G = np.array([0.0, 0.0])
K = (2*b1 + b2) / 3
M_pt = b1 / 2

path_points = [G, K, M_pt, G]
path_labels = [r'$\Gamma$', 'K', 'M', r'$\Gamma$']
N_path = 200

kpath = []
tick_positions = [0]
for seg in range(len(path_points)-1):
    for t in np.linspace(0, 1, N_path, endpoint=(seg == len(path_points)-2)):
        kpath.append((1-t)*path_points[seg] + t*path_points[seg+1])
    tick_positions.append(len(kpath)-1)

bands = np.zeros((len(kpath), 2))
for idx, k in enumerate(kpath):
    H = haldane_hamiltonian(k[0], k[1])
    vals = np.linalg.eigvalsh(H)
    bands[idx] = vals

x_path = np.arange(len(kpath))
axes[1].plot(x_path, bands[:, 0], color='#a855f7', linewidth=2, label='Lower band')
axes[1].plot(x_path, bands[:, 1], color='#22c55e', linewidth=2, label='Upper band')
axes[1].fill_between(x_path, bands[:, 0], alpha=0.15, color='#a855f7')

for tp in tick_positions:
    axes[1].axvline(x=tp, color='#334155', linewidth=0.8, linestyle='--')

axes[1].set_xticks(tick_positions)
axes[1].set_xticklabels(path_labels, color='#94a3b8', fontsize=12)
axes[1].set_ylabel('Energy / t', color='#94a3b8', fontsize=12)
axes[1].set_title(f'Haldane Model Bands (C = {chern:.0f})', color='white', fontsize=14)
axes[1].tick_params(colors='#94a3b8')
axes[1].legend(fontsize=10, facecolor='#1a1a2e', edgecolor='#334155', labelcolor='white')

plt.tight_layout()
plt.savefig('berry_curvature.png', dpi=130, bbox_inches='tight',
            facecolor=fig.get_facecolor())
print("Berry curvature plot saved.")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Fortran Implementation: Berry Phase of a Two-Level System

Numerical calculation of the Berry phase for a spin-$\frac{1}{2}$ system in a magnetic field tracing a cone on the unit sphere. The Berry phase equals half the solid angle subtended by the loop: $\gamma = -\pi(1 - \cos\theta_0)$.

Berry Phase: Two-Level System

Fortran

Compute Berry phase around closed loops in parameter space and verify against analytic solid-angle formula

program.f90129 lines

program berry_phase_two_level
  implicit none
  integer, parameter :: dp = selected_real_kind(15)
  integer, parameter :: N_points = 500
  real(dp), parameter :: pi = 4.0_dp * atan(1.0_dp)

! ──────────────────────────────────────────────────────────
  ! Berry Phase for a Two-Level System
  !
  ! H(R) = R_x*sigma_x + R_y*sigma_y + R_z*sigma_z
  !
  ! Parameter space: R = (sin(theta)cos(phi), sin(theta)sin(phi), cos(theta))
  ! Loop: theta = theta_0 fixed, phi: 0 -> 2*pi
  !
  ! Analytic result: gamma = -pi(1 - cos(theta_0))
  !   = -1/2 * solid angle subtended by the loop
  ! ──────────────────────────────────────────────────────────

integer :: i, n_loops
  real(dp) :: theta_0, phi, dphi
  real(dp) :: Rx, Ry, Rz, E_plus, E_minus
  complex(dp) :: psi_prev(2), psi_curr(2), overlap_product
  complex(dp) :: overlap
  real(dp) :: berry_phase_num, berry_phase_exact
  real(dp) :: norm_val
  real(dp) :: theta_vals(9), error

write(*,'(A)') '========================================================'
  write(*,'(A)') '  Berry Phase in a Two-Level System'
  write(*,'(A)') '  H = R . sigma,  |R| = 1 (unit sphere)'
  write(*,'(A)') '  Loop: fixed theta, phi = 0 -> 2*pi'
  write(*,'(A)') '========================================================'
  write(*,'(A)') ''

dphi = 2.0_dp * pi / real(N_points, dp)

! Test multiple values of theta_0
  theta_vals = (/ 0.1_dp, 0.3_dp, 0.5_dp, 0.8_dp, 1.0_dp, &
                  1.2_dp, 1.5_dp, 2.0_dp, 2.5_dp /)
  n_loops = 9

write(*,'(A)') '  theta_0 (rad)  |  Numerical      |  Exact           |  Error'
  write(*,'(A)') '  ──────────────────────────────────────────────────────────────'

do i = 1, n_loops
    theta_0 = theta_vals(i)

call compute_berry_phase(theta_0, N_points, berry_phase_num)

! Exact: gamma = -pi * (1 - cos(theta_0))
    berry_phase_exact = -pi * (1.0_dp - cos(theta_0))

! Normalize both to [-pi, pi]
    berry_phase_num = modulo(berry_phase_num + pi, 2.0_dp*pi) - pi
    berry_phase_exact = modulo(berry_phase_exact + pi, 2.0_dp*pi) - pi

error = abs(berry_phase_num - berry_phase_exact)

write(*,'(2X,F10.4,A,F14.8,A,F14.8,A,ES12.4)') &
      theta_0, '    | ', berry_phase_num, '  | ', berry_phase_exact, '  | ', error
  end do

write(*,'(A)') ''
  write(*,'(A)') 'Solid angle interpretation:'
  write(*,'(A)') '  Berry phase = -(1/2) * solid angle subtended by loop'
  write(*,'(A)') '  For theta_0 = pi/2: Omega = 2*pi => gamma = -pi'
  write(*,'(A)') '  For theta_0 = pi:   Omega = 4*pi => gamma =  0 (mod 2*pi)'
  write(*,'(A)') '========================================================'

contains

subroutine compute_berry_phase(theta_0, N, gamma)
    real(dp), intent(in) :: theta_0
    integer, intent(in) :: N
    real(dp), intent(out) :: gamma
    integer :: k
    real(dp) :: phi_k, Rx_k, Ry_k, Rz_k, dp_local
    complex(dp) :: u_prev(2), u_curr(2), prod, ov
    real(dp) :: nrm

dp_local = 2.0_dp * pi / real(N, dp)
    prod = cmplx(1.0_dp, 0.0_dp, dp)

! Ground state at phi = 0
    phi_k = 0.0_dp
    Rx_k = sin(theta_0) * cos(phi_k)
    Ry_k = sin(theta_0) * sin(phi_k)
    Rz_k = cos(theta_0)
    call ground_state(Rx_k, Ry_k, Rz_k, u_prev)

do k = 1, N
      phi_k = real(k, dp) * dp_local
      ! For the last point, wrap back to phi = 0
      if (k == N) phi_k = 0.0_dp

Rx_k = sin(theta_0) * cos(phi_k)
      Ry_k = sin(theta_0) * sin(phi_k)
      Rz_k = cos(theta_0)
      call ground_state(Rx_k, Ry_k, Rz_k, u_curr)

! Overlap <u_prev | u_curr>
      ov = conjg(u_prev(1))*u_curr(1) + conjg(u_prev(2))*u_curr(2)
      prod = prod * ov / abs(ov)  ! keep only the phase

u_prev = u_curr
    end do

gamma = aimag(log(prod))
  end subroutine

subroutine ground_state(Rx, Ry, Rz, u)
    real(dp), intent(in) :: Rx, Ry, Rz
    complex(dp), intent(out) :: u(2)
    real(dp) :: R_mag, theta_loc, phi_loc, nrm

R_mag = sqrt(Rx**2 + Ry**2 + Rz**2)
    theta_loc = acos(Rz / R_mag)
    phi_loc = atan2(Ry, Rx)

! Ground state (lower eigenvalue -|R|) of H = R.sigma:
    ! |-> = (sin(theta/2), -cos(theta/2)*exp(i*phi))
    u(1) = cmplx(sin(theta_loc/2.0_dp), 0.0_dp, dp)
    u(2) = -cos(theta_loc/2.0_dp) * cmplx(cos(phi_loc), sin(phi_loc), dp)

nrm = sqrt(abs(u(1))**2 + abs(u(2))**2)
    u = u / nrm
  end subroutine

end program berry_phase_two_level

Click Run to execute the Fortran code

Code will be compiled with gfortran and executed on the server

Key Takeaways

The Berry phase $\gamma_n = \oint \mathcal{A}_n \cdot d\mathbf{k}$ is a geometric phase depending only on the path in parameter space.
The Berry curvature $\Omega_n = \nabla_\mathbf{k}\times\mathcal{A}_n$ is gauge-invariant, peaks near degeneracies, and vanishes when both TRS and inversion hold.
The Chern number $C_n = \frac{1}{2\pi}\int_\text{BZ}\Omega_n\,d^2k \in \mathbb{Z}$ classifies bands and gives quantized Hall conductance via TKNN.
The anomalous velocity $\dot{\mathbf{r}}_\text{anom} = -\dot{\mathbf{k}}\times\boldsymbol{\Omega}_n$ produces intrinsic anomalous Hall effect.
Berry phase underlies the modern theory of polarization and orbital magnetization.

← Part IV Overview Chapter 2: Quantum Hall Effect →