10 Condensed matter physics
10.1 Crystals & Reciprocal Space
Crystalline solids are periodic in space. Periodicity gives you enormous leverage: lattices turn sums into selection rules, Fourier transforms quantize , and diffraction fingerprints structure. This section lays down the lattice–basis language, reciprocal vectors, Brillouin zones, Miller indices, and the diffraction toolkit (Bragg, Laue, Ewald, structure factors). We’ll tie together fcc/bcc duality, selection rules for common crystals, and quick ways to not get lost in indices.
10.1.1 Lattice, basis, crystal structure
A Bravais lattice is the set of points
with three linearly independent primitive vectors . A crystal structure = lattice plus a basis of atoms at positions inside a primitive cell; the atomic positions are .
- Wigner–Seitz cell: the set of points closer to a given lattice point than any other.
- Conventional vs primitive: conventional cubic cells are nicer to look at; primitive cells are minimal-volume and great for theory.
Pro tip: never say “the diamond lattice.” Diamond has an fcc lattice with a two-point basis .
10.1.2 Symmetry: point groups and space groups (one-screenshot tour)
- Point group: rotations, reflections, inversion that leave at least one point fixed.
- Space group: point-group ops + translations + screw rotations and glide reflections.
- In 3D there are 14 Bravais lattices, 32 crystallographic point groups, and 230 space groups. For band theory you mostly need the lattice and the little groups of high-symmetry points.
10.1.3 Reciprocal lattice and Brillouin zone
Define the reciprocal primitive vectors by . In vector form
The reciprocal lattice consists of vectors with integers .
The first Brillouin zone (BZ) is the Wigner–Seitz cell of the reciprocal lattice. It is the primitive cell in -space and the natural domain for Bloch wavevectors.
- Duality: fcc bcc. The reciprocal of fcc is bcc and vice versa. Simple cubic is self-dual.
10.1.4 Bloch theorem in one line
In a periodic potential , single-particle eigenstates are Bloch waves:
Wavevectors can be chosen in the first BZ. Translating leaves physics invariant.
10.1.5 Planes, directions, and indices (Miller crash course)
- Directions use square brackets for a lattice vector proportional to . Families of equivalent directions: .
- Planes use parentheses , the Miller indices, defined as reciprocals of intercepts along . Families: .
- Spacing to plane family is tied to the reciprocal vector :
For cubic etc., .
Hexagonal often uses Miller–Bravais with to make basal-plane symmetry manifest.
10.1.6 Diffraction: Bragg, Laue, Ewald are the same vibe
Elastic scattering conserves and transfers crystal momentum equal to a reciprocal vector:
Bragg condition in real space:
Laue condition in -space is exactly . Ewald sphere construction: draw a sphere of radius anchored at the tail of ; reciprocal-lattice points on the sphere give allowed reflections.
Structure factor encodes the basis:
with atomic form factors . The diffracted intensity is
The lattice sum enforces belongs to the reciprocal lattice; decides if it actually lights up.
10.1.7 Selection rules for common structures
- Simple cubic (sc) with one atom at :
- Body-centered cubic (bcc), basis :
- Face-centered cubic (fcc), basis :
- NaCl (rock salt): fcc lattice with two-species basis . Then
Intensity alternates with depending on -contrast.
- Diamond: fcc lattice with basis . Extra phase leads to additional systematic zeros beyond fcc.
10.1.8 fcc–bcc duality and high-symmetry points
If are fcc vectors, the reciprocal generate a bcc lattice in -space (and vice versa). Learn the standard labels:
- fcc BZ: center, at face centers, at BZ corners, mid-edges.
- bcc BZ: , , , points.
These are where band extrema and degeneracies love to appear and where experimentalists measure ARPES and quantum oscillations.
10.1.9 Real ↔ reciprocal geometry cheats
- Volume relation:
- Plane normal: the reciprocal vector is normal to the planes.
- Angle between planes and :
- Spacing again: .
10.1.10 Worked mini-examples
(a) Reciprocal of bcc is fcc
Start with bcc conventional cubic cell of side and primitive vectors etc. Compute and show the reciprocal points sit at face centers of a cubic net with spacing .
(b) for cubic
Using , show . Bragg gives ; solve for the first-order reflection angle for given .
(c) fcc systematic absences
Evaluate for , , families and confirm only “all even” or “all odd” survive.
(d) Diamond zeros
Diamond has . Show that is extinct but is allowed, consistent with X-ray tables.
(e) Ewald construction sketch
For (Cu K), draw a circle of radius through the origin of the reciprocal fcc lattice and mark which vectors on the circle predict reflections for a given crystal orientation.
10.1.11 Common pitfalls
- Confusing lattice with structure: lattice is the periodic scaffold; the basis puts atoms on it.
- Wrong units in reciprocal space: crystallographers often use -free convention; here . Pick one and stay loyal.
- Miller vs direction indices: is a plane, is a direction. For cubic, ; not true in general.
- Forgetting the form factor: even if , decays with , reshaping intensities.
- BZ confusion: the “reduced zone” scheme folds into the first BZ; if bands look doubled, you probably mixed extended vs reduced plots.
10.1.12 Minimal problem kit
- Derive the reciprocal vectors from and prove
- Show that the Laue condition and Bragg’s law are equivalent via geometry
- Compute the structure factor and selection rules for NaCl and CsCl; compare which vanish in each
- For a hexagonal lattice with and , write and obtain in terms of
- Locate the high-symmetry points , , , for fcc in Cartesian -coordinates given , and sketch the first BZ
10.2 Lattice Vibrations & Phonons
Atoms in a crystal aren’t statues; they wiggle. Those collective wiggles are phonons: quantized normal modes of the lattice. In this section we go from 1D chains to full 3D dynamical matrices, split acoustic vs optical branches, compute dispersions, group velocities, and densities of states, then quantize. We also hit thermal transport, anharmonicity, thermal expansion via the Grüneisen parameter, electron–phonon coupling, and quick experimental fingerprints (neutron, Raman, IR).
10.2.1 Harmonic lattice: setting the stage
Let label Bravais sites and the basis atom index inside a unit cell. The displacement of atom is . In the harmonic approximation the potential energy is quadratic in displacements:
where are force constants and . Newton:
Try plane-wave solutions . You get an eigenproblem
with the dynamical matrix
For atoms per cell you get branches: 3 acoustic as , and optical.
10.2.2 1D monatomic chain: the cleanest demo
Masses connected by springs with lattice spacing . Let be the displacement of site . Equation of motion:
Ansatz gives the dispersion
Key takeaways:
- Acoustic branch only, linear near with sound speed
- Brillouin zone is ; is periodic with
- Group velocity controls energy flow
Density of states in 1D for a single branch:
So diverges at BZ edges where (van Hove singularity).
10.2.3 1D diatomic chain: acoustic + optical
Two alternating masses , same spring , unit cell length containing both atoms. With nearest-neighbor coupling, the two-branch dispersion is
- The minus sign is the acoustic branch with as
- The plus sign is the optical branch with
At the acoustic mode moves the two atoms in phase; the optical mode moves them out of phase. A gap opens at the zone boundary.
10.2.4 3D crystals: polarizations and sound speeds
For (monatomic Bravais) you have 3 acoustic branches: one longitudinal (LA) and two transverse (TA). Near they are sound waves with linear dispersion
In cubic crystals the speeds relate to elastic constants and mass density (continuum limit):
More generally, the Christoffel equation gives direction-dependent velocities.
For you add optical branches that typically have weak -dispersion near the zone center. If the basis carries opposite charges, some optical modes are IR active (see §10.14).
10.2.5 Quantization: phonons as bosons
Promote each normal mode to a harmonic oscillator:
with . The displacement operator is
Phonon occupation follows Bose statistics with chemical potential in equilibrium.
10.2.6 Heat capacity and Debye revisit
Total phonon energy per volume with . Low- acoustic modes Debye law in 3D; high Dulong–Petit (see §9.6 for the full Debye integral). Optical modes add Schottky-like bumps when their gaps become thermally accessible.
10.2.7 Thermal transport: from kinetic theory
In a crystal, heat is mostly carried by acoustic phonons. Kinetic theory gives
where is the mode heat capacity, the group velocity, and the lifetime. Scattering channels:
- Phonon–phonon (anharmonic). Split into Normal (N) processes conserving crystal momentum and Umklapp (U) with that degrade heat current
- Boundary scattering (dominant in nanoscale or at very low )
- Isotope/disorder scattering
- Phonon–electron and phonon–defect scattering in metals
Empirically, often rises as at very low (ballistic/Boundary), peaks, then falls roughly as at high (Umklapp-limited).
10.2.8 Anharmonicity, frequency shifts, and lifetimes
Real crystals aren’t perfectly harmonic. Cubic and quartic terms cause:
- Thermal expansion (quasi-harmonic effect)
- Frequency renormalization
- Finite lifetimes via three-phonon and four-phonon scattering
A handy bulk measure is the Grüneisen parameter
The volumetric thermal expansion coefficient relates to and the bulk modulus :
(Quasi-harmonic approximation; works surprisingly well for many simple solids.)
10.2.9 Electron–phonon coupling: the vibes behind resistivity and pairing
Displacements modulate the electronic Hamiltonian , producing an interaction . Two archetypes:
- Deformation potential (short-range): couples to LA phonons
- Fröhlich (long-range polar): in ionic crystals, LO phonons create macroscopic electric fields that couple strongly to electrons
Consequences:
- Metal resistivity: at low , Bloch–Grüneisen law ; at high roughly
- Superconductivity: virtual phonon exchange provides an effective attraction that can overcome Coulomb repulsion near the Fermi surface (BCS)
10.2.10 Instabilities and soft modes: Peierls in 1D
In 1D metals the electronic susceptibility peaks at . Electron–phonon coupling softens the phonon at that wavevector; below a critical the soft mode drives a lattice distortion with doubled period (dimerization), opening a gap: a Peierls transition. In higher dimensions, charge-density waves and structural transitions follow similar logic.
10.2.11 Experimental fingerprints of phonons
- Inelastic neutron scattering (INS): measures the dynamic structure factor , mapping full and eigenvectors
- Raman spectroscopy: probes zone-center optical phonons; symmetry selection rules apply
- Infrared (IR) absorption/reflectivity: IR-active optical phonons in polar crystals split LO vs TO at (Lyddane–Sachs–Teller relation)
- Brillouin light scattering: accesses long-wavelength acoustic phonons (GHz)
- Specific heat and thermal conductivity vs : bulk thermodynamic signatures of DOS and scattering
10.2.12 Worked mini-examples
(a) Group velocity in the monatomic chain
From , compute and show as .
(b) DOS edge singularity
Use to show the 1D DOS diverges at .
(c) Diatomic gap at
Evaluate and for the diatomic chain and interpret the eigenvectors.
(d) Debye
Assuming three acoustic branches with linear dispersion up to a Debye cutoff , integrate to get at low .
(e) Umklapp threshold
Show that an Umklapp process requires to leave the first BZ; estimate the minimum temperature where such become thermally populated in a simple cubic with Debye wavevector .
(f) Thermal expansion from
Starting from , show that in the quasi-harmonic approximation.
10.2.13 Common pitfalls
- “Optical requires ions.” No: you only need atoms per cell. IR activity requires polarity; optical phonons exist regardless
- Phase vs group velocity: energy transport uses , not
- Forgetting selection rules: Raman sees only Raman-active symmetries; IR sees IR-active modes; neutrons couple to nuclear positions and magnetic moments
- Mixing BZ schemes: reduced vs extended zone pictures change where branches appear but not the physics
- Overusing harmonic approximation: heat transport and expansion need anharmonicity; include N vs U processes when modeling
10.2.14 Minimal problem kit
- Derive the monatomic and diatomic 1D dispersions and sketch both within the first BZ
- Build a dynamical matrix for a 2-atom cubic basis with nearest-neighbor central forces and extract along a high-symmetry line
- Compute for a Debye model in and show the corresponding low- specific heat scalings
- Use Matthiessen’s rule with boundary, isotope, and Umklapp scattering rates to fit a synthetic curve; identify the peak temperature scaling with sample size
- Estimate the Bloch–Grüneisen temperature and the law crossover for a metal with given and
10.3 Drude–Sommerfeld Metals
Metals look messy at the atomic scale, but their electrons act surprisingly lawful. Drude gave the first kinetic picture (1900): classical electrons as a gas of charged billiard balls with occasional collisions. Sommerfeld (1928) fixed the big misses using Fermi–Dirac statistics. Together they still power a shocking amount of metal physics: DC/AC conductivity, Hall effect, magnetotransport, specific heat, Pauli paramagnetism, and the Wiedemann–Franz law.
10.3.1 Drude model: one relaxation time, huge mileage
Assume conduction electrons of (effective) mass and charge experience randomizing collisions with mean time .
Equation of motion (steady drive) for average velocity under fields :
- DC conductivity at with current density :
Define mobility so .
- AC (Drude) conductivity for field :
Real part is a Lorentzian of width ; imaginary part switches sign at .
- Plasma frequency (free-electron pole in optics):
For , metals reflect; for , they transmit like a dielectric.
- Skin depth in the good-conductor limit ():
10.3.2 Hall effect and the conductivity tensor
With , the Drude solution gives the conductivity tensor
Inverting gives the resistivity tensor. Two key outputs:
- Hall coefficient (weak field):
- Magnetoresistance in the single-band, isotropic Drude model is zero: despite shrinking, because the tensor inversion cancels it. Real materials show finite magnetoresistance due to multiple bands, anisotropy, or variations.
10.3.3 Drude thermal transport and Wiedemann–Franz (preview)
A simple kinetic estimate gives electronic thermal conductivity (once we introduce ). Drude naively predicts the Wiedemann–Franz ratio
But only with Sommerfeld statistics do we get the correct Lorenz number
We’ll derive this properly after quantizing the electron gas.
10.3.4 The fix: Sommerfeld free-electron gas
Electrons are fermions. Put of them in a box of volume and fill momentum states up to the Fermi sphere.
- Fermi wavevector and energy (spin degeneracy ):
- Density of states per volume for a parabolic band:
At low , only a shell of width around is thermally active. This explains why most classical Drude thermodynamics failed.
10.3.5 Electronic specific heat: linear in
Using Sommerfeld expansion,
Linear matches experiments and is tiny compared to at room temperature because only of electrons participate.
10.3.6 Wiedemann–Franz, done right
Kinetic theory with Fermi statistics yields
Combine with and to obtain, for elastic scattering with energy-independent ,
Deviations of from diagnose inelastic scattering, strong energy dependence of , or additional heat carriers (phonons, magnons).
10.3.7 Pauli paramagnetism and Landau diamagnetism (sketch)
- Pauli paramagnetism from spin polarization in a field :
Small, -independent to leading order.
- Landau diamagnetism from orbital motion in quantizing fields gives
for a free-electron gas. In real metals the net also includes core diamagnetism and possible Van Vleck terms.
10.3.8 Thermopower and the Mott formula (one-liner)
The Seebeck coefficient (voltage from a temperature gradient) vanishes in Drude’s classical picture but, with Fermi statistics,
This Mott formula makes a sensitive probe of particle–hole asymmetry and scattering near .
10.3.9 What Drude gets wrong (and why Sommerfeld wins)
- Heat capacity: classical Drude predicts ; data show . Sommerfeld fixes it.
- Hall sign: single-band Drude gives (electrons). Many materials have positive due to hole-like Fermi surfaces or multiband compensation.
- Magnetoresistance: zero in single-band Drude, but finite in reality from anisotropy, multiple pockets, or .
- High-frequency optics: Drude alone misses interband absorption and bound-charge contributions; add Lorentz oscillators for realistic .
10.3.10 Connecting to band structure
Replace by the band effective mass from the dispersion near the Fermi surface:
Anisotropic makes a tensor even at . Real transport uses Boltzmann theory linearized around the Fermi surface with (see §10.7).
10.3.11 Optical sum rule and spectral weight
The f-sum rule ties the integrated Drude weight to carrier density:
In correlated systems, Drude weight can be suppressed and shifted to mid-infrared interband-like features—often read as a fingerprint of reduced quasiparticle weight.
10.3.12 Worked mini-examples
(a) Numbers you can trust
Take , , . Compute , , , and at .
(b) Hall bar
For a sample with thickness , width , current along , and field , show the Hall voltage in the Drude single-band limit. Estimate for copper with the numbers from (a).
(c) Sommerfeld
Derive using the Sommerfeld expansion and express in terms of and .
(d) Wiedemann–Franz from Boltzmann
Linearize the Boltzmann equation with elastic impurity scattering and show when is constant near .
(e) Mott thermopower sign
Assuming near , show and interpret vs .
10.3.13 Common pitfalls
- Forgetting effective mass: always use from the band structure, not bare , in transport formulae.
- Reading magnetoresistance from : the measured can be -independent in single-band Drude even though shrinks; invert the tensor.
- Assuming universally: multiband compensation and warped Fermi surfaces can flip the sign or reduce magnitude.
- Overusing constant : energy- or angle-dependent scattering changes , , and MR; be cautious outside low- impurity-limited regimes.
- Mixing optical and transport masses: cyclotron, thermodynamic, and band masses can differ in anisotropic or interacting systems.
10.4 Bloch Theorem & Band Structure
Periodic potentials take messy atomic chaos and compress it into clean momentum-space rules. The headline: Bloch’s theorem says electrons in a crystal move as plane waves modulated by a lattice-periodic envelope. From that, gaps pop open at Brillouin-zone boundaries, dispersions define velocities and effective masses, and both “nearly free” and “tight-binding” lenses become the same story, zoomed in from opposite sides.
10.4.1 Bloch states and the central equation
For a single electron in a periodic potential ,
Expand the lattice-periodic part in reciprocal vectors :
and the potential as . Insert into Schrödinger’s equation to get the central equation
It’s an infinite linear system in the plane-wave coefficients . Truncate sensibly and you get bands.
10.4.2 Zone folding and reduced vs extended zone
Free electrons have parabolas that repeat every reciprocal vector: and are equivalent. You can:
- Extended-zone: keep as endlessly repeated parabolas in -space
- Reduced-zone: fold all momenta back into the first Brillouin zone; copies become separate bands indexed by
When , states that were degenerate under folding mix and gaps open at the crossings.
10.4.3 Gaps from Bragg reflection: a 2×2 demo
At a boundary where and are degenerate for the free electron, keep only those two plane waves in the central equation. The eigenproblem is
with . At the exact zone boundary where , the two energies split as
so the band gap is
Moral: stronger Fourier components of the lattice potential at reciprocal vectors carve bigger gaps at the corresponding Brillouin-zone planes.
10.4.4 Kronig–Penney in one line (1D intuition)
For a 1D periodic array of barriers (lattice constant ), the dispersion satisfies a transcendental relation
where in the wells and encodes barrier strength. Regions where the right-hand side has magnitude are allowed bands; where it exceeds 1 are forbidden gaps. Tuning moves you smoothly from nearly free () to tight-binding ().
10.4.5 Group velocity, effective mass, and DOS
Band dispersion is not just eye candy; its derivatives are physics:
- Group velocity
- Effective mass tensor near an extremum
Curvature ⇒ light mass; curvature ⇒ heavy mass. Negative curvature near a valence-band maximum means holes with positive charge and positive mass run the transport.
- Density of states for a single band in 3D
The DOS diverges where the gradient vanishes along the constant-energy surface (van Hove singularities).
10.4.6 Nearly free vs tight-binding: two zoom levels, same film
- Nearly free electron (NFE): start from parabolas and open gaps perturbatively where they cross. Good when .
- Tight-binding (TB): start from localized orbitals with hopping between sites; bands emerge from bonding/antibonding combinations. Best when orbitals are compact and is small compared to on-site energies.
Both are limits of the same Bloch machinery. NFE explains gap openings and light masses in simple metals; TB nails band shapes in covalent/ionic solids and will be our §10.5 focus (Wannier functions live here).
10.4.7 method: local expansion near a band edge
Near a high-symmetry , expand the Hamiltonian as . Second-order perturbation yields
where are Bloch states at . This connects effective mass to interband matrix elements and symmetry selection rules (e.g., heavy vs light holes in semiconductors).
10.4.8 Symmetry, degeneracies, and spin–orbit spice
- Crystal symmetry protects degeneracies at high-symmetry points and along lines. Breaking symmetry (strain, fields) splits them.
- Inversion and time reversal together enforce Kramers doublets with spin.
- Spin–orbit coupling (SOC) mixes spin and orbital character; with broken inversion you can get Rashba-like splittings. We stay light here and return to topology in §10.11.
10.4.9 Fermi surfaces and transport pointers
The Fermi surface is the constant-energy surface . Its curvature and connectivity control:
- Cyclotron masses and quantum oscillations (de Haas–van Alphen, Shubnikov–de Haas)
- Anisotropic conductivity and Hall effect (via velocities and lifetimes on the FS)
- Nesting vectors that can amplify density-wave instabilities
Bloch dispersion + scattering time is the whole transport game in §10.7.
10.4.10 Wannier functions: real-space avatars of bands
A unitary transform of Bloch states over the BZ gives Wannier functions localized near cell . Smooth gauge choices yield maximally localized Wanniers that act like tight-binding orbitals and are perfect for model building, interpolation, and computing real-space responses. We build on this in §10.5.
10.4.11 ARPES and band mapping (how we see )
Angle-resolved photoemission (ARPES) measures the occupied part of directly. In 2D materials it’s almost literal; in 3D there’s a wrinkle. Combine with quantum oscillations and optical conductivity to reconstruct the low-energy band structure that transport actually feels.
10.4.12 Worked mini-examples
(a) Gap size from a single Fourier component
Take a 1D cosine potential with . Do the 2×2 central-equation truncation at and show .
(b) Effective mass near a band minimum
For in 1D, extract and discuss how the term skews DOS but not at the minimum.
(c) DOS singularity in 2D
On a square lattice TB band , locate the saddle points and show the logarithmic van Hove singularity in .
(d) Zone-folding check
Start from a free-electron parabola and fold by into the first BZ. Draw where degeneracies occur and mark which get lifted when .
(e) for a direct-gap semiconductor
At , use symmetry to argue that conduction and valence states have opposite parity so the linear-in- coupling vanishes; show why the conduction-band mass is then dominated by second-order coupling to remote bands.
10.4.13 Common pitfalls
- Mixing reduced and extended zones: keep track of whether a “crossing” is real or just a folded copy.
- Equating curvature with speed: velocity is first derivative, not curvature; flat band heavy mass and slow carriers.
- Forgetting the basis: selection rules and gap sizes depend on the Fourier components , which encode the atomic basis and bonding.
- Assuming parabolic everywhere: effective-mass fits work only near extrema; far away, nonparabolicity and anisotropy matter.
- Ignoring SOC where it’s large: in heavy elements or broken-inversion crystals, SOC reshapes band order and splitting—critical for topology.
10.4.14 Minimal problem kit
- Derive the central equation from the Bloch expansion and write it as a matrix eigenproblem for a finite set of
- Perform the 2×2 gap-opening calculation at a generic zone plane with normal ; obtain
- Compute tensor for an ellipsoidal constant-energy surface and relate it to cyclotron mass in a field along an arbitrary axis
- For a 1D Kronig–Penney with delta barriers, plot allowed bands for several and identify the NFE and TB limits
- Using a simple TB model on honeycomb, show how Dirac cones arise and how a staggered sublattice potential opens a mass gap
10.5 Tight-Binding & Wannier: From Local Orbitals to Bands
When electrons are more “homebody” than “frequent flyer,” tight-binding (TB) is the lens. Start with localized atomiclike orbitals, let them hop between sites, and you get bands with shapes and symmetries carved by real-space chemistry. Wannier functions are the clean, orthonormal real-space avatars of Bloch bands, perfect for building models, computing responses, and adding fields or disorder. This section builds TB from scratch, adds multi-orbital and Slater–Koster rules, handles magnetic fields via Peierls phases, and shows how Wanniers tie to polarization and edge physics (SSH preview).
10.5.1 From orbitals to Bloch sums
Pick orbitals centered on basis sites in cell , with orbital label (e.g., ).
Form Bloch sums
Expand eigenstates as . Matrix elements
with hopping , overlap . Solve the generalized eigenproblem
In orthogonal TB we approximate .
10.5.2 1D monatomic chain (orthogonal TB)
One orbital per site, lattice constant , on-site energy , nearest-neighbor hopping .
Bandwidth . Effective mass near an extremum at is
Sign of sets band curvature and which end is bonding vs antibonding.
10.5.3 Diatomic chain and the SSH teaser
Two sites per cell, at , at , hoppings and alternating. In the basis
For the dispersion is
Gap at the zone boundary is . The SSH limit () supports topological edge states when with open boundaries; we revisit the topology in §10.11.
10.5.4 Square and honeycomb lattices (classics)
Square (1 orbital, nn)
Honeycomb (graphene, 1 per site, nn). Two-site basis and three nn vectors :
Eigenvalues
Touch at Dirac cones with linear dispersion and Fermi velocity .
10.5.5 Multi-orbital TB and Slater–Koster rules
Hoppings depend on orbital orientation relative to the bond. Slater–Koster parameterizes two-center integrals using direction cosines between sites.
Examples:
- –:
- –: projected on the axis
- –: and channels
For a bond along :
Cross-hoppings like times the appropriate direction cosine. This symmetry bookkeeping lets you write for real crystals quickly and keeps you from illegal hoppings.
10.5.6 Beyond nearest neighbors and van Hove moves
Adding next-nearest neighbors etc. reshapes curvature and DOS singularities. On the square lattice
Tuning slides the van Hove saddle through , a common route to instabilities (density waves, superconductivity) in layered materials.
10.5.7 Magnetic field: Peierls substitution
Minimal coupling in TB enters via Peierls phases on hoppings
So the product of phases around a plaquette equals the magnetic flux . On a square lattice with rational flux , translation symmetry enlarges to a magnetic unit cell and the spectrum fractures into Hofstadter’s butterfly.
10.5.8 Disorder and real-space Green’s functions
Real materials have impurities. Add site randomness or bond randomness . The TB Hamiltonian in real space is perfect for:
- Anderson localization studies (see §10.13)
- Surface/edge spectra via finite ribbons and recursive Green’s functions
- Transport with Landauer–Büttiker, where and comes from real-space TB + leads
10.5.9 Wannier functions: localized avatars of bands
Given Bloch states in a band (or group) subspace, define Wannier functions
where is a gauge unitary within the subspace. Properties:
- Orthonormal and cell-localized. The choice of controls spread
- Maximally localized Wanniers (MLWFs) choose to minimize total spread
- For isolated bands in trivial topology, exponentially localized Wanniers exist and reproduce the band exactly
Tight-binding from Wannier. Hoppings are matrix elements between Wanniers
The interpolation gives bands indistinguishable from the ab initio subspace but at TB cost.
10.5.10 Polarization and Wannier centers (Berry twist)
Electronic polarization per cell of an insulator relates to Wannier centers :
Equivalent Berry-phase formula uses the Bloch gauge. Practically, track Wannier center flow under adiabatic changes to compute polarization differences and piezoelectric responses (topology returns in §10.11).
10.5.11 Edge physics with TB: the SSH demo in 15 seconds
Open a finite chain with alternating . In the topological phase you get midgap edge modes exponentially localized with decay . Chiral symmetry pins them at until perturbations that break sublattice symmetry lift them.
10.5.12 Worked mini-examples
(a) Wannier from a 1D cosine band
For , show the Wannier is proportional to a modified Bessel asymptotically, with localization length set by the nearest gap.
(b) Graphene Fermi velocity
From the honeycomb nn Hamiltonian, linearize near to get and extract .
(c) Slater–Koster on a cubic bond
For a cubic crystal with one and three orbitals per site and nn bonds along axes, write using .
(d) Peierls phase check
On a square lattice in Landau gauge , show the hopping along picks a phase with and the -index, while hoppings stay real.
(e) Overlap matters
Include a nearest-neighbor overlap in the 1D chain and solve the generalized eigenproblem to show the dispersion becomes .
(f) SSH edge states
Diagonalize a 40-site open SSH chain numerically (TB) and plot the wavefunctions of the two midgap states; verify exponential localization and parity on opposite ends.
10.5.13 Common pitfalls
- Forgetting overlap: many intro TBs set ; when orbitals are not orthogonal, use the generalized problem or orthogonalize first
- Illegal hoppings: if symmetry forbids a matrix element, do not sneak it in. Use Slater–Koster or group theory to zero the right terms
- Gauge sloppiness in Wanniers: random phases across give delocalized Wanniers. MLWFs fix the gauge—do that before trusting real-space hoppings
- Peierls without flux accounting: ensure the product of link phases equals per plaquette, independent of gauge
- Overfitting TB: too many parameters hide physics and blow up outside the fit window. Start minimal; add neighbors only as data demand
10.5.14 Minimal problem kit
- Build the 1D and diatomic TB dispersions and extract effective masses around band edges
- Write the honeycomb nn Hamiltonian, find , linearize to obtain the Dirac Hamiltonian and
- Use Slater–Koster to construct a TB for graphene including and sectors; discuss why decouples in the flat sheet
- Implement Peierls phases for square-lattice TB at rational flux and plot the Hofstadter spectrum for
- Generate MLWFs for a two-band toy model, compute , and show Wannier-interpolated matches the exact band to machine precision
- Compute the polarization change from a sliding dimerization in the SSH model via Wannier center shift and compare with the Berry-phase result
10.6 Semiconductors & Doping
Semiconductors are metals on “hard mode”: carriers are scarce until you tune the band gap, sprinkle dopants, or shine photons. This section builds the effective-mass picture, intrinsic vs doped statistics, drift–diffusion with the Einstein relation, recombination channels, optical fingerprints, and the minimalist pn-junction toolkit. We keep it band-structure real while staying device-friendly.
10.6.1 Band edges, effective masses, and densities of states
Near a conduction minimum and a valence maximum,
Define the band gap . Effective-mass tensors may be anisotropic; often one uses DOS masses to define 3D densities of states
These feed the carrier statistics. e carrier statistics.
10.6.2 Intrinsic carriers and the mass-action law
With Fermi level , nondegenerate (Boltzmann) populations are
Multiply to get the mass-action law
In an intrinsic crystal and sits near midgap skewed by DOS masses.
10.6.3 Doping: donors, acceptors, and charge neutrality
Substitutional dopants add shallow hydrogenic levels. Effective-mass theory gives
Similarly for acceptors with . At room , shallow dopants are usually ionized.
For an n-type crystal with donor density and acceptors ,
Assuming full ionization in the nondegenerate regime, and . Then
In the majority-carrier limit , and .
At low , partial ionization leads to a freeze-out region where ; use Fermi–Dirac occupancy of dopant levels.
10.6.4 Mobilities, scattering, and conductivity
Each carrier has mobility with drift velocity . Conductivity is
Scattering sources:
- Lattice (phonons): typically in 3D
- Ionized impurities: at low
- Matthiessen’s rule for rates:
Two-carrier Hall coefficient (weak field) shows compensation:
Sign reveals the dominant carrier type.
10.6.5 Drift–diffusion and Einstein relation
Carrier continuity with drift and diffusion reads
In nondegenerate statistics,
This Einstein relation is the semiconductor face of fluctuation–dissipation.
10.6.6 Recombination, lifetime, and quasi-Fermi levels
Out of equilibrium (illumination, injection), electrons and holes can be described by quasi-Fermi levels with
Net recombination–generation channels:
- Radiative (band-to-band):
- Shockley–Read–Hall (trap-mediated): roughly with trap parameters
- Auger:
Define minority-carrier lifetime via small-signal . Lifetimes set diffusion lengths that control photodiode and solar-cell collection.
10.6.7 Direct vs indirect gaps, absorption, and excitons
Optical absorption coefficient near the band edge:
- Direct gap:
- Indirect gap: phonon-assisted, rises slowly and shows temperature dependence
Binding of electron–hole pairs gives excitons with hydrogenic series. Effective-mass theory estimates
with reduced mass . Excitonic peaks decorate the absorption edge in clean crystals and 2D materials.
10.6.8 Minimal pn junction: built-in potential, depletion, and diode law
Bring p and n regions together; carriers diffuse, leaving behind ionized dopants and forming a depletion region with built-in field. At equilibrium,
Ignoring free charge in depletion (abrupt junction), the total depletion width under applied bias is
Forward bias reduces the barrier; the ideal Shockley diode equation:
where scales inversely with diffusion lengths and doping. Reverse bias nearly clamps at until breakdown mechanisms kick in (Zener/tunneling or avalanche).
10.6.9 Metal–semiconductor contacts: ohmic vs Schottky in one breath
A metal on a lightly doped semiconductor forms either:
- Ohmic contact if the barrier is thin/low enough for tunneling or if heavy doping bends bands into easy injection
- Schottky barrier with rectification otherwise. Barrier height is set by metal work function, semiconductor electron affinity, and interface pinning
Transport across a Schottky contact at forward bias often follows thermionic emission with .
10.6.10 Temperature regimes (the three-act play)
For a doped semiconductor as rises:
- Freeze-out: dopants partially ionized,
- Extrinsic: plateau where (or ) and
- Intrinsic: once surpasses dopant densities exponentially
Resistivity and Hall data vs let you read off activation energies and dominant scattering.
10.6.11 Worked mini-examples
(a) Intrinsic carrier density
Given , , at , estimate for Si using the form.
(b) Fermi level in doped Si
For and at , find and locate relative to via .
(c) Einstein relation check
With measured in Si at , compute and a diffusion length for .
(d) Depletion width
For a Si pn junction with , , , compute at and at .
(e) Radiative vs SRH
Given , , estimate and compare to an SRH rate with .
10.6.12 Common pitfalls
- Forgetting mass-action: in equilibrium even in doped material; violations signal injection or illumination
- Using Einstein beyond its lane: standard needs nondegenerate statistics; degenerate semiconductors require Fermi–Dirac corrections
- Assuming full ionization: at low or for deeper dopants, solve dopant occupancy self-consistently
- Mixing DOS and conductivity masses: want DOS masses; transport needs conductivity (tensor) masses and scattering
- Ideal diode dogma: real diodes have series resistance, recombination in depletion (ideality factor ), and interface states
10.7 Fermi Surfaces & Transport
The Fermi surface (FS) is the VIP section of a metal: at low temperature almost everything transporty—conductivity, Hall, magnetoresistance, quantum oscillations—comes from quasiparticles living right on it. Geometry in -space becomes measurable physics in the lab. This section builds the semiclassical equations of motion, the Boltzmann/Chambers machinery, magnetotransport including multiband and compensation, and the quantum‐oscillation toolkit (Onsager + Lifshitz–Kosevich). We’ll also flag how FS topology (closed vs open orbits) controls magnetoresponse.
10.7.1 Semiclassical equations of motion (no Berry yet)
For a band in fields the semiclassical dynamics are
With you get cyclotron motion in -space along constant-energy contours perpendicular to ; real-space velocity is . We postpone Berry curvature and anomalous velocity to §10.11.
10.7.2 Cyclotron frequency and mass
Closed orbits have frequency
where the cyclotron mass is set by the extremal -space orbit area normal to :
For an isotropic parabola , ; for anisotropic ellipsoids, depends on field orientation.
10.7.3 Boltzmann transport: relaxation time and Chambers
Linearizing the Boltzmann equation in a steady field with relaxation time gives the current
In Chambers’ formula the vector mean free path is a history average along the orbit
leading to the conductivity tensor (schematically)
where is the orbit-averaged, decay-weighted velocity. For an isotropic, single- metal this collapses to Drude.
10.7.4 Hall angle, mobility tensor, mean free path
Define mobility and Hall angle (single band). The mean free path
Close to the Ioffe–Regel limit metals “bad” out, and semiclassical transport starts to crumble.
10.7.5 Magnetoresistance: why single-band Drude lies flat
In a single, isotropic band with constant , the resistivity is field-independent even though shrinks; tensor inversion cancels it. Real materials show MR because of (i) multiple carrier types, (ii) FS anisotropy/warping, (iii) angle-dependent , or (iv) open orbits. A practical summary is Kohler’s rule
If a single scattering time controls everything, curves at different collapse when plotted vs .
10.7.6 Two-band model (electrons + holes)
Let and be densities and mobilities; , charge signs. The conductivities sum:
Invert to get . Key regimes:
- Compensation with high mobilities ⇒ large, nearly quadratic unsaturating MR
- Dominant carrier ⇒ Hall sign fixed, MR smaller and tends to saturate
Multiband fits to and are the day-one way to read carrier types and mobilities.
10.7.7 Open vs closed orbits; angle dependence
- Closed orbits on closed FS pockets give conventional and oscillations.
- Open orbits (e.g., warped quasi-1D sheets) allow carriers to stream along a direction; MR becomes strongly anisotropic and can grow linearly or quadratically without saturation depending on geometry.
- AMRO (angle-dependent MR oscillations) in layered metals arise when tilting threads Landau orbits across a warped quasi-2D FS; minima occur at Yamaji angles set by interlayer warping.
10.7.8 Quantum oscillations: Onsager + Lifshitz–Kosevich
Landau quantization turns orbit areas into discrete levels. The Onsager relation ties the oscillation frequency (in ) to the extremal FS cross-section area perpendicular to :
The Lifshitz–Kosevich (LK) amplitude has three classic factors:
with
gives the mass plot to read , the Dingle factor for scattering, and the spin-splitting modulation. The phase can encode Berry phase for Dirac/Weyl systems (see §10.11).
10.7.9 Cyclotron resonance and optical probes
In clean samples, an AC field picks up resonances at giving directly. Optical conductivity decomposes into a Drude piece and interband contributions; the f-sum rule sets the total spectral weight (see §10.3). In layered conductors, interlayer conductivity vs tilt reveals FS warping.
10.7.10 Practical FS metrology (ARPES vs oscillations)
- ARPES maps occupied bands and FS contours directly in ; great for 2D and surfaces, trickier for bulk .
- Quantum oscillations see only extremal areas but with exquisite precision and bulk sensitivity.
- Quantum oscillation angle sweeps reconstruct full FS shapes from ; combine with ARPES to de-ambiguate pockets and 3D warping.
10.7.11 Thermoelectric and Nernst hints
Near the FS, the Mott formula relates the Seebeck coefficient to the energy derivative of conductivity at :
The Nernst effect (transverse voltage from a thermal gradient in ) is extremely sensitive to FS curvature and quasiparticle lifetime asymmetries—often a canary for small pockets.
10.7.12 Strong-field and quantum-limit quirks
At very large , carriers collapse into a few lowest Landau levels. In the quantum limit (only occupied), MR and thermoelectrics can behave non-Drude, sometimes linear in or with sign flips. Dirac/Weyl materials can show a chiral-anomaly contribution to longitudinal MR when (details in §10.11).
10.7.13 Worked mini-examples
(a) Onsager area from frequency
You measure for . Compute the extremal area
and an equivalent circular Fermi wavevector .
(b) LK mass plot
Given oscillation amplitudes vs at fixed , fit to extract from .
(c) Two-band MR
Take , , . Evaluate and show quadratic growth up to several tesla without saturation.
(d) AMRO/Yamaji
For a quasi-2D FS with interlayer dispersion , show that MR minima occur when with integer , giving from the angular period.
(e) Ioffe–Regel check
With , , and in a single band, estimate , compute , and comment on .
10.7.14 Common pitfalls
- Mixing masses: (band curvature), (cyclotron), and optical/thermodynamic masses differ in anisotropic or multiband systems
- Forgetting tensor inversion: trends in don’t translate 1:1 to ; always invert the full tensor
- Assuming single : hot/cold spots and anisotropic scattering wreck Kohler collapse and mimic multiband MR
- Misreading oscillation phase: Berry phase extraction needs careful background subtraction and index conventions; spin zeros can shift apparent phase
- Ignoring geometry: sample misalignment and current jetting can fake huge linear MR in high-mobility, compensated crystals
- Open-orbit blinders: seeing “saturating MR” and blaming impurities when the real culprit is FS topology and field orientation
10.8 Magnetism
Magnetism comes in two big flavors: localized moments that behave like tiny quantum spins coupled by exchange, and itinerant magnetism where the Fermi sea itself polarizes. From Curie–Weiss laws to spin waves and ferromagnetic resonance, this section builds the standard playbook: Heisenberg models and exchange mechanisms (direct, superexchange, double exchange, RKKY), Stoner ferromagnets, magnons in ferro/antiferromagnets, anisotropy and domains, and the Mermin–Wagner no-go in 1D/2D without symmetry breaking fields. We keep the math tight and the takeaways actionable.
10.8.1 Magnetic moments, factors, and units
An electron carries spin and possibly orbital contributions. The magnetic moment operator is
with Bohr magneton and . In solids, quenching of orbital angular momentum by the crystal field often leaves effective spin moments, sometimes with spin–orbit dressing that yields anisotropy.
The magnetization is moment per volume. Linear response defines susceptibility via (isotropic case), with .
10.8.2 Local-moment models: Ising/Heisenberg/XXZ
For localized spins on a lattice the minimal Hamiltonians are
- Ising
- Heisenberg (SU(2)-symmetric)
- XXZ anisotropy
Here favors ferromagnetism (FM) and antiferromagnetism (AFM) in our sign convention.
Curie–Weiss phenomenology for paramagnets: above the ordering temperature,
The sign of hints at dominant FM () or AFM () exchange.
10.8.3 Exchange mechanisms: where comes from
- Direct exchange: overlap of neighboring local orbitals; Pauli + Coulomb yield an energy difference between parallel and antiparallel spin alignments.
- Superexchange: virtual hopping via a nonmagnetic ligand (e.g., TM–O–TM). For a half-filled Hubbard-like situation,
often antiferromagnetic.
- Double exchange (mixed valence manganites): itinerant electrons align local core spins to maximize hopping, favoring FM and metallicity.
- RKKY (in metals with localized moments): conduction electrons mediate an oscillatory coupling
driving complex order or spin glass when disorder scrambles phases.
10.8.4 Itinerant ferromagnetism: Stoner criterion
In a simple band, exchange lowers energy if spin polarization splits the bands by . The Stoner mean-field yields a criterion
with an interaction parameter and the DOS per spin at the Fermi level. Consequences:
- Finite Pauli susceptibility enhanced to
- Spin-wave stiffness emerges in the broken-symmetry phase (see §10.8.6)
Itinerant magnets mix band structure with interactions; ARPES, quantum oscillations, and optical sum rules help disentangle.
10.8.5 Mean-field ordering temperatures
On a -coordination lattice with exchange between near neighbors, Weiss mean-field predicts
- FM Heisenberg
- AFM (bipartite)
Real are reduced by fluctuations, especially in low dimensions or small .
10.8.6 Spin waves (magnons) in ferromagnets
In a FM Heisenberg model, small transverse oscillations are magnons. Holstein–Primakoff mapping gives, to leading order,
with dispersion near
Here is the spin-wave stiffness set by and lattice geometry, and comes from anisotropy (otherwise Goldstone gapless). Thermal excitation of magnons reduces the magnetization:
This is the Bloch law in 3D FMs.
Specific heat of magnons in 3D scales as at low (quadratic dispersion → DOS ).
10.8.7 Antiferromagnets: two-sublattice spin waves
For a collinear AFM on a bipartite lattice, linear spin-wave theory yields two branches with linear dispersion at small :
Velocity depends on , , and lattice spacing. Anisotropy produces a small gap. The uniform susceptibility often shows a cusp at ; below, transverse and longitudinal components differ.
10.8.8 Magnetic anisotropy, demagnetizing fields, and domains
- Magnetocrystalline anisotropy from spin–orbit coupling selects easy axes/planes. Minimal model
affects coercivity and domain structure.
- Demagnetizing fields arise because at surfaces acts like “magnetic charge.” Shape defines a demag tensor with field .
- Domains form to balance wall energy and magnetostatic energy. Domain walls have width and energy where is exchange stiffness and anisotropy.
10.8.9 Ferromagnetic resonance (FMR) and dynamics
Magnetization dynamics obey Landau–Lifshitz–Gilbert (LLG)
with gyromagnetic ratio and damping . For a uniform mode in an in-plane thin film (Kittel formula, simple case),
FMR reads out anisotropy fields and damping.
10.8.10 Dimensionality: Mermin–Wagner and Kosterlitz–Thouless
Continuous-symmetry magnets in strictly 1D/2D with short-range exchange cannot order at finite (no gap, Goldstone modes kill long-range order): the Mermin–Wagner theorem. Workarounds:
- Add anisotropy (Ising-like discretization), long-range dipolar terms, or interlayer coupling
- In 2D XY, expect a Kosterlitz–Thouless transition with vortex unbinding and algebraic order rather than conventional
10.8.11 Magnetotransport cliffs notes
- Anomalous Hall effect (AHE) in FMs: , with intrinsic Berry-curvature and extrinsic skew/side-jump mechanisms.
- Spin waves and resistivity: additional or contributions at low in FMs from magnon scattering.
- Spin caloritronics: spin Seebeck and magnon-drag signals link magnon and charge/heat transport.
10.8.12 Worked mini-examples
(a) Curie–Weiss slope
Given , , number density , show that the high- slope of is with . Fit data to extract or .
(b) Spin-wave stiffness in a cubic FM
Nearest-neighbor Heisenberg on cubic lattice (spacing ): show
with and . Expand for small to extract .
(c) AFM two-sublattice dispersion
For a bipartite AFM with nn exchange , derive the linear small- magnon velocity up to geometry factors.
(d) Stoner enhancement
Given , estimate and comment on proximity to instability .
(e) Domain wall width
With and , estimate and .
(f) RKKY sign changes
For , compute the first few zeros/sign flips of to see at what separations the coupling toggles FM/AFM.
10.8.13 Common pitfalls
- Sign conventions for : check whether your Hamiltonian uses or ; FM/AFM flip with the choice.
- Curie–Weiss overreach: it is a high- law; near and below critical fluctuations break mean-field straightness.
- Ignoring anisotropy: in 2D, no finite- order without anisotropy or long-range terms; tiny matters a lot.
- Confusing Pauli and Curie: itinerant paramagnets (Pauli) have -independent ; local moments (Curie) go like . Real compounds can have both.
- Mixing demag and intrinsic: raw magnetization loops are geometry-dependent; correct for demagnetization to get intrinsic .
- Magnon gap blind spot: any observed finite-frequency intercept in FM/AFM spin-wave spectra signals anisotropy or Zeeman gaps—don’t fit as purely gapless.
10.8.14 Minimal problem kit
- Derive Curie–Weiss susceptibility from a Weiss field for spins and connect to and microscopic
- Compute and within mean-field for a cubic Heisenberg FM; compare to Bloch at low and discuss regimes of validity
- Do linear spin-wave theory (Holstein–Primakoff) for a 1D FM chain, obtain and show why 1D FM lacks finite- long-range order (infrared divergence)
- Starting from LLG, derive the Kittel formula for a thin film with uniaxial anisotropy and identify
- For a simple parabolic band with Stoner parameter , compute the magnetization vs at mean-field and sketch the phase diagram in space
10.9 Superconductivity
Superconductors are the drama queens of condensed matter: below a critical temperature they drop their DC resistance to (effectively) zero and expel magnetic fields (Meissner effect). Microscopically, electrons form Cooper pairs that condense into a phase-coherent quantum state. This section builds the phenomenology (London, Meissner), the mesoscopic field theory (Ginzburg–Landau), the microscopic BCS results (gap, , DOS), and the real-world signatures: Type I vs II, vortices, critical fields/currents, and Josephson effects.
10.9.1 London equations and the Meissner effect
Phenomenology first. The London equations for the supercurrent are
Combine with Maxwell to eliminate and get magnetic field screening
so a static field decays inside the superconductor as . This Meissner effect distinguishes superconductors from ideal (non-superconducting) perfect conductors, which would merely freeze whatever they had when cooled.
10.9.2 Ginzburg–Landau (GL) theory
GL is a Landau free energy for a complex order parameter (Cooper-pair wavefunction) minimally coupled to the vector potential:
Minimizing yields the GL equations
Key GL lengths:
The dimensionless GL parameter classifies superconductors:
- Type I if (single critical field , complete Meissner until abrupt breakdown)
- Type II if (mixed state with vortices between and )
10.9.3 BCS essentials: pairing, gap, and
BCS shows that any weak, retarded attraction between electrons near (e.g., phonon-mediated) causes a Cooper instability and a paired ground state with an energy gap .
Zero-temperature gap and critical temperature (isotropic -wave, weak coupling):
Universal BCS ratio:
The quasiparticle density of states (per spin) becomes
with coherence peaks at . The specific heat jumps by ; the low- behavior is activated for fully gapped -wave.
BCS coherence length (clean limit):
10.9.4 Critical fields, currents, and the mixed state
Type I has a thermodynamic critical field that satisfies .
Type II features three scales:
For , magnetic field penetrates via Abrikosov vortices carrying one flux quantum . Vortex cores have radius where is suppressed; circulating supercurrents decay over . Vortices form a triangular lattice in clean materials absent pinning.
Critical current density (GL depairing estimate near ):
In practice, vortex motion (flux flow) limits current; strong pinning raises the critical current used in applications.
10.9.5 Flux quantization and phase stiffness
Single-valuedness of the condensate phase around a loop gives
which implies flux quantization in a multiply connected superconductor. The phase rigidity is the origin of persistent currents and of the Josephson effects.
10.9.6 Josephson effects
Two superconductors separated by a thin barrier (SIS) support a supercurrent without voltage:
- DC Josephson:
where is the phase difference and the critical current. With magnetic field, interference across a junction of width gives a Fraunhofer pattern .
- AC Josephson: a DC voltage makes the phase run , yielding an AC current at
Microwave irradiation produces Shapiro steps at voltages . Two junctions in a loop form a SQUID with period vs applied flux, enabling ultrasensitive magnetometry.
In SNS junctions, transport below the gap proceeds via Andreev reflection, forming discrete Andreev bound states that carry the supercurrent.
10.9.7 Proximity and mesoscopic effects
A normal metal in good contact with a superconductor inherits pair correlations over a proximity length (dirty limit ). Subgap conductance of an NS interface is enhanced by Andreev reflection, approaching a factor of 2 in the ideal BTK limit at zero bias for a perfectly transparent contact.
10.9.8 Unconventional and high- notes
Not all superconductors are -wave. -wave pairing (cuprates) features line nodes, giving power-law thermodynamics ( in 3D) and a V-shaped tunneling DOS. -wave and other exotic states can host Majorana bound states in topological contexts. In many unconventional cases, the “glue” is not phonons but electron–electron interactions (spin fluctuations), and BCS ratios deviate from 3.53.
10.9.9 Experimental fingerprints
- Transport: zero resistance below ; finite voltage only when current exceeds or vortices move
- Magnetization: Meissner expulsion, irreversible hysteresis from vortex pinning in Type II
- Thermodynamics: specific-heat jump at ; condensation energy from integrating
- Spectroscopy: tunneling/STM coherence peaks at ; ARPES gap symmetry; optical conductivity with missing area transferred to the zero-frequency delta (Ferrell–Glover–Tinkham sum rule)
- Microwave/THz: AC conductivity shows a superfluid delta plus a finite-frequency quasiparticle response
- Muon/NMR: penetration depth, vortex lattice, and gap symmetry from relaxation rates
10.9.10 Worked mini-examples
(a) Penetration depth from carrier density
Using , estimate for and .
(b) estimate
A Type II superconductor has . Compute .
(c) Coherence length from gap and
With and , estimate and compare to the GL near .
(d) Josephson frequency
What microwave frequency results from across a junction Using , show .
(e) Flux quantization in a ring
A superconducting loop of area shows SQUID oscillations with period . Verify and compute the period.
(f) Depairing current (order of magnitude)
Given and , estimate and compare with a measured limited by pinning.
10.9.11 Common pitfalls
- Meissner vs perfect conductor: superconductors expel newly applied fields; perfect conductors merely lock in whatever field is present on cooling
- Forgetting the : Cooper pairs carry charge ; flux quantum is , not
- Type I/II mix-ups: classification uses ; it is not about “how strong the field is” but interfacial energy sign
- Ignoring vortices: in Type II, finite resistance under current + field often comes from vortex motion, not a destroyed gap
- Using clean-limit formulas in dirty samples: and renormalize with impurity scattering; check mean free path vs
- Assuming -wave: gap symmetry matters; power laws in , , and nodal quasiparticles betray unconventional pairing
10.9.12 Minimal problem kit
- Starting from the London equations, derive and show exponential decay of in a half-space
- From GL, compute the surface energy between normal and superconducting regions and reproduce the Type I/II boundary at
- Show that flux quantization follows from single-valuedness of and derive
- Use the weak-coupling BCS gap equation to obtain and the ratio
- Derive by linearizing the GL equation in a field and mapping to Landau levels for the order parameter
- Compute the Fraunhofer envelope for a rectangular Josephson junction and extract the junction area from the period
10.10 Strong Correlation & Mott Physics
If weakly interacting electrons are chill, strongly correlated ones are messy drama: localization without band gaps, gigantic effective masses, bad-metal transport, and spectral weight teleporting across eV scales when you tweak doping by a few percent. The mascot is the Mott insulator—an odd creature that should conduct by band theory (half-filled band) but refuses because electrons won’t share a site when on-site repulsion is large. This section builds the Hubbard family, superexchange and the – model, the Mott transition (Brinkman–Rice vs. Slater), spectral weight transfer, multi-orbital/Hund physics, and a DMFT crash course, with experimental fingerprints you can spot at 10 paces.
10.10.1 Hubbard model: the minimalist drama
Put electrons on a lattice, let them hop with amplitude , and penalize double occupancy by :
- At half filling and (bandwidth ), you get a metal.
- For , charge motion freezes: each site wants one electron to dodge , producing a Mott insulator with local moments.
- Doping away from half filling introduces mobile carriers in a strongly interacting background → strange transport and, in some materials, high- superconductivity.
10.10.2 From Hubbard to Heisenberg: superexchange
When at half filling, virtual hops create an effective antiferromagnetic exchange between neighboring spins:
This is superexchange. Intuition: parallel spins block exchange by Pauli, antiparallel can virtually double-occupy and return, lowering energy.
10.10.3 – model: doped Mott playbook
Doping a Mott insulator while keeping large gives the constrained model
where forbids double occupancy. The no-double-occupancy constraint is the core difficulty; mean-field, slave boson/fermion, or numerical methods tackle it. In 2D, competition between kinetic energy and spin singlet formation seeds -wave pairing tendencies.
10.10.4 Slater vs Mott: why insulators aren’t all the same
Two routes to a gap at half filling:
- Slater insulator: a weak-coupling antiferromagnet. Long-range order doubles the unit cell, folding bands and opening a mean-field gap , which closes above the Néel temperature .
- Mott insulator: gap from local physics; the material stays insulating without long-range order, even above .
How to tell in the lab: if the high- phase is still insulating with a large activation energy and ARPES shows lower/upper Hubbard bands far from the chemical potential, you’re in Mott-land (e.g., NiO). If the gap collapses right at with small spectral rearrangement, it’s Slater-ish.
10.10.5 Spectral function anatomy: Hubbard bands + quasiparticle
Define the one-electron spectral function . In a correlated metal near a Mott transition:
- A narrow quasiparticle (QP) peak of weight straddles the Fermi level.
- Incoherent Hubbard bands sit at , carrying weight .
The quasiparticle weight is
with the self-energy. As correlations strengthen, , effective mass , and at the Mott point —the QP dissolves, leaving only Hubbard bands.
10.10.6 Brinkman–Rice picture and the MIT
In a Fermi-liquid mean-field for correlations, the kinetic energy is renormalized by and the Mott metal–insulator transition (MIT) occurs when
The divergence of and collapse of Drude weight are hallmarks. In real 3D lattices, the transition at finite is often first order with hysteresis (two solutions coexist: metal and insulator), meeting at a critical endpoint—beautifully captured by DMFT.
10.10.7 Optical and ARPES fingerprints (how to spot Mottness)
- Optics: Drude weight shrinks as , while a mid-infrared feature and Hubbard-band absorption grow. Spectral weight transfer violates simple band sum rules—weight moves over eV scales when you dope by a few percent.
- ARPES: coherent QP peak near plus broad incoherent background and high-energy waterfalls. Doping a Mott insulator pulls new low-energy states from the upper Hubbard band and reshapes the dispersion (strong kink structures if electron–phonon also matters).
- Transport: resistivity above the Mott–Ioffe–Regel (MIR) limit (“bad metal”), huge thermopower asymmetries, Hall coefficient that changes sign with and doping.
10.10.8 Charge-transfer vs Mott insulators (Zaanen–Sawatzky–Allen)
In transition-metal oxides, the lowest unoccupied state might live on ligand orbitals rather than the metal . Compare (on ) to the charge-transfer energy (from O to TM ):
- Mott–Hubbard insulator: ; gap mostly between lower and upper -Hubbard bands.
- Charge-transfer insulator: ; gap is from O to upper band. Doping holes lands you on oxygen-like states (hello, cuprates).
This classification explains wildly different doping responses and optical edges across oxides with similar periodic tables.
10.10.9 Hund metals and multi-orbital spice
With several active orbitals, Hund’s coupling aligns spins on the same atom. Effects:
- Suppresses interorbital charge fluctuations, reducing even when is moderate → “Hund metal” with heavy quasiparticles and large scattering.
- Fosters orbital selectivity: some orbitals Mott-localize while others stay itinerant (orbital-selective Mott phase).
- Competes/cooperates with crystal fields and SOC, key to ruthenates, iron-based superconductors, and 4d/5d oxides.
10.10.10 DMFT in one page: local self-energy, global consequences
Dynamical Mean-Field Theory (DMFT) maps a lattice model to a quantum impurity embedded in a self-consistent bath. The self-energy is taken local (momentum-independent): —exact in infinite coordination, remarkably good in 3D.
Self-consistency (Bethe lattice with half bandwidth ) reads
where is the impurity hybridization, the non-interacting DOS. Iterate: solve impurity → get → recompute → update .
What DMFT nails
- Three-peak spectrum (QP + Hubbard bands), collapse at , first-order MIT with coexistence.
- Temperature evolution of resistivity including bad-metal behavior.
- Doping asymmetries and spectral-weight transfer.
- With LDA band inputs (LDA+DMFT): realistic materials, e.g., VO, SrVO, iron pnictides.
10.10.11 1D twist: Luttinger liquids and spin–charge separation
In strictly 1D, Fermi liquids are unstable: correlations produce a Luttinger liquid with separate spin and charge modes. Signatures: power-law tunneling DOS, different velocities for spin vs charge excitations seen in ARPES or cold-atom analogues. The 1D Hubbard model is exactly solvable (Bethe Ansatz) and shows a Mott gap at half filling for any .
10.10.12 Kondo, heavy fermions, and the lattice cousin (one-paragraph cameo)
A local moment in a metal is screened by conduction electrons below the Kondo temperature , forming a singlet with a narrow resonance at . A periodic array of such moments (Kondo lattice) hybridizes at low into heavy quasiparticles with ; depending on filling and symmetry you can get heavy Fermi liquids, unconventional superconductivity, or Kondo insulators. (Deep dive would live in its own section; headline here: “correlations can also light up bands rather than only split them.”)
10.10.13 Worked mini-examples
(a) Superexchange sign and scale
For a two-site Hubbard dimer at half filling, do second-order perturbation in and show and why it’s AFM.
(b) Brinkman–Rice mass blow-up
Assuming near the MIT, show diverges as and the Drude weight scales with .
(c) DMFT Bethe self-consistency
With semicircular DOS and , derive .
(d) Optical spectral weight transfer
Use the f-sum rule to argue why integrated low-energy spectral weight decreases with at fixed filling, and sketch how doping restores Drude weight while also creating mid-IR absorption.
(e) ZSA classification
Given and charge-transfer energy , identify the insulator type and predict whether hole doping creates O- or TM--like carriers.
(f) Hund metal toy
In a two-orbital Hubbard–Hund model with moderate but finite , explain qualitatively why interorbital charge fluctuations are suppressed and drops compared to the case.
10.10.14 Common pitfalls
- Calling every AFM insulator “Mott”: if the gap collapses at , it’s likely Slater. Check the paramagnetic phase.
- Equating large with insulator at any filling: away from half filling you can have a correlated metal with small , not necessarily an insulator.
- Forgetting spectral-weight transfer: simple band pictures conserve low-energy weight under doping; Mottness does not.
- Using band effective masses in correlated fits: transport and thermodynamics see ; use the renormalized values.
- Momentum-blind DMFT overreach: single-site DMFT can’t capture -dependent pseudogaps or short-range order well; cluster DMFT/DCA or diagrammatics are needed.
- Mixing Mott and Anderson: interaction-driven localization (Mott) differs from disorder-driven (Anderson); both can coexist (Mott–Anderson), but diagnostics differ (e.g., interaction tuning vs. disorder tuning).
10.10.15 Minimal problem kit
- Starting from the Hubbard model at half filling and , derive the Heisenberg via Schrieffer–Wolff transformation
- Implement a Hubbard-I self-energy and plot the emergence of upper/lower Hubbard bands with growing ; compare to DMFT’s three-peak structure
- For a semicircular DOS, solve the DMFT equations at within iterated perturbation theory and chart ; locate
- Compute the doping evolution of the optical conductivity in a toy three-peak model (QP + two Hubbard bands); verify spectral-weight transfer
- Classify a given oxide as Mott–Hubbard vs charge-transfer using tabulated and and predict the sign/magnitude trend of the Seebeck coefficient upon light hole doping
- Build a two-orbital Hubbard–Hund mean-field and map where an orbital-selective Mott phase appears as a function of , , and crystal-field splitting
10.11 Topological Band Theory
Some band insulators are “just” boring gaps. Others hide protected boundary states, quantized responses, and Berry-curvature superpowers. That’s topological band theory: use geometry in -space (Berry phase, curvature) plus symmetries (time reversal, inversion, chiral) to classify phases. Bulk invariants predict edge modes. Change the invariant? You must close the gap somewhere—cue phase transitions with Dirac or Weyl nodes. Let’s build the toolkit.
10.11.1 Berry phase, connection, and curvature
For a Bloch band smooth in :
- Berry connection
- Berry curvature
- Berry phase along a closed loop in the Brillouin zone (BZ)
Gauge note: depends on phase choice of , and loop phases modulo do not. In 1D, the Berry phase of occupied bands controls polarization (modern theory): is a Berry phase over the BZ (see §10.5.10).
10.11.2 Chern number and quantum Hall
In 2D, the Chern number of band is
For a fully filled set of bands,
This is the TKNN integer quantum Hall result: a bulk topological invariant predicts a boundary chiral edge count. Change → the gap must close somewhere.
10.11.3 Anomalous Hall: Berry curvature without
Even at , broken time-reversal symmetry (TRS) yields an intrinsic anomalous Hall effect
Here is the Fermi function. Think of as a “magnetic field in -space” deflecting wavepackets (§10.7 deferred Berry terms).
10.11.4 1D winding and the SSH warm-up
For a two-band, chiral-symmetric 1D model with , define . The winding number
counts how many times wraps the origin. Bulk gives edge zero modes in the SSH chain (§10.5.3).
10.11.5 2D Dirac mass flip → Chern insulator (Haldane vibe)
Consider a single massive Dirac cone
The valence-band Chern number is
Two cones (as on honeycomb) contribute; if their masses have opposite signs, total → a Chern insulator with one chiral edge channel.
10.11.6 Time-reversal symmetry and topological insulators
With TRS ( for electrons), bands come in Kramers pairs. In 2D and 3D, insulators are classified by a index (2D) or (3D). If inversion symmetry is present, use the Fu–Kane parity shortcut:
- At each TRIM , multiply inversion eigenvalues of occupied Kramers pairs to get
- Then
3D strong index is ; we omit the weak indices detail here. Strong TI () hosts an odd number of Dirac cones on any surface, protected by TRS.
Without inversion, compute a Wilson loop of the Berry connection (a.k.a. Wannier charge-center flow). If the hybrid Wannier centers wind across the BZ, the phase is topological.
10.11.7 Bulk–boundary correspondence (in one graphic sentence)
Nontrivial bulk invariant ⇒ robust boundary states that cannot be gapped without either closing the bulk gap or breaking the protecting symmetry.
- Chern insulator (): chiral edge modes (unidirectional).
- 2D TI: helical edge (Kramers pair counterpropagating).
- 3D strong TI: odd number of surface Dirac cones.
- Crystalline topological phases (protected by mirror, rotation) host symmetry-selective boundary states.
10.11.8 Dirac and Weyl semimetals: topology without a gap
Gapless cousins:
- Weyl nodes are monopoles of Berry curvature in 3D with integer chiral charge . They require breaking either TRS or inversion to separate nodes in -space. Consequences: Fermi arcs on surfaces and the chiral anomaly (longitudinal negative MR when ). The anomaly equation reads
- Dirac semimetals superpose opposite-chirality Weyl nodes pinned together by symmetry; break the symmetry and they split into Weyl pairs.
10.11.9 Rashba, spin–momentum locking, and surface Dirac cones
Spin–orbit coupling (SOC) plus broken inversion gives Rashba splitting
Topological surface states in TIs show spin–momentum locking (helical textures), suppressing backscattering from nonmagnetic disorder. Magnetic perturbations gap the surface Dirac cone, enabling quantized magneto-electric responses (axion physics below).
10.11.10 Axion electrodynamics and magneto-electric response
In 3D TIs, the effective EM action adds
with (trivial) or (strong TI) modulo if TRS holds. At a surface where jumps by , you get a half-quantized surface Hall conductance when TRS is explicitly broken at the boundary.
10.11.11 Wilson loops and Wannier flow (how we compute invariants)
Discretize along a closed path and build the Wilson loop
where is the non-Abelian Berry connection of the occupied subspace. Eigenphases of give Wannier centers; their flow vs transverse reveals Chern ( winds by per BZ) or (odd vs even partner-switching).
10.11.12 Orbital magnetization and modern polarization (Berry redux)
Two Berry-phase byproducts:
- Polarization (insulators) is a Berry phase of occupied states across the BZ; only changes are physical and quantized under adiabatic cycles (Thouless pump).
- Orbital magnetization in crystals can be expressed as BZ integrals involving and band energies; this links to the intrinsic AHE and to modern magneto-electric coefficients.
10.11.13 Minimal symmetry-class crib sheet (tiny slice of the tenfold way)
| Class | Symmetry (TRS, PHS, chiral) | 2D invariant | 3D invariant |
|---|---|---|---|
| A (no symm) | –, –, – | (Chern) | 0 |
| AII (TRS with ) | ✓, –, – | (strong ) | |
| AIII (chiral only) | –, –, ✓ | (winding) | 0 |
The full tenfold way covers superconductors (BdG) and more invariants; we focus on the workhorses for electrons.
10.11.14 Experimental fingerprints
- Transport: integer QHE plateaus; quantized edge conductance for Chern insulators; weak anti-localization from Berry phase; negative longitudinal MR in Weyls (careful: geometry artifacts exist).
- ARPES: surface Dirac cones, Fermi arcs, Rashba splitting; spin-resolved ARPES confirms textures.
- STM/STS: absence of backscattering in QPI for helical surfaces; Landau-level ladders revealing Berry phase .
- Optics: giant Kerr/Faraday angles with magnetized TI surfaces; cyclotron resonance with nontrivial phase offsets.
- Quantum oscillations: phase shifts (Landau index plot intercept) tied to Berry phase, with spin-splitting caveats.
10.11.15 Worked mini-examples
(a) Chern of a two-band Dirac
For , compute and integrate to show per cone; two inequivalent cones with opposite give .
(b) Parity test
Given inversion eigenvalues at four 2D TRIM as for the occupied Kramers pairs, evaluate ⇒ topological.
(c) Wilson loop wind
On a Chern band with , show the Wannier center angle advances by as goes across the BZ. Contrast with a trivial band.
(d) SSH end states
Compute the winding number for . Show for and plot edge spectra for an open chain.
(e) Weyl monopole charge
Evaluate the flux of through a small sphere around a Weyl node and show it equals (i.e., charge ).
(f) Surface Hall from axion
Assuming a TI with a magnetized surface so , integrate across the boundary to obtain .
10.11.16 Common pitfalls
- Gauge traps: is gauge-dependent; only and loop phases mod are physical. Use smooth gauges or Wilson loops for numerics.
- Using Fu–Kane without inversion: the parity shortcut demands inversion symmetry; otherwise rely on Wilson loops or flow of Wannier centers.
- Confusing Rashba with topology: Rashba splitting alone is not a TI; you need a bulk gap plus a nontrivial index.
- Misreading oscillation phases: extracting Berry phase from Landau fan diagrams requires accounting for Zeeman, warping, and 3D geometry.
- Edge ≠ surface every time: symmetry protects; break it at the boundary and you can gap the state (e.g., magnetic disorder gapping a TI surface).
- Chern in TRS systems: with unbroken TRS and no magnetic order, the total Chern number of occupied bands must be zero; the index can still be nontrivial.
10.11.17 Minimal problem kit
- Derive the TKNN formula starting from Kubo and show for a clean 2D insulator
- Implement a discrete-BZ Chern number computation via link variables and validate on the two-band Dirac model
- For a centrosymmetric SOC insulator, compute Fu–Kane from tabulated parities at TRIM; then confirm by Wilson-loop flow
- Build a two-node Weyl model, compute surface Green’s function for a slab, and plot Fermi arcs
- Add a surface Zeeman term to a 3D TI Dirac cone and show the opening of a gap and half-quantized Hall conductance
- Simulate an adiabatic Thouless pump in a 1D lattice and show the quantized charge transport equals the Chern number of the torus
10.12 Mesoscopic Transport & Quantum Coherence
Shrink a conductor until its size rivals coherence and mean free paths, and transport flips from “traffic flow” to wave interference. This is the mesoscopic regime: electrons remain phase coherent across the device, conductance becomes sample-specific yet reproducible, and quantum noise reveals transmission statistics. The workhorses here are Landauer–Büttiker transport, Aharonov–Bohm interference, weak localization/antilocalization, and universal conductance fluctuations. We’ll keep one eye on experiments (QPCs, rings, nanowires) and one on formulas you can deploy.
10.12.1 Landauer: conductance as transmission
For a phase-coherent two-terminal conductor at low temperature and small bias, current is carried by independent quantum channels with transmissions :
Factor 2 is spin degeneracy (lifted by Zeeman/SOC). For a quantum point contact (QPC) in the ballistic limit, transverse modes open one by one so and
giving quantized conductance steps vs gate voltage.
Multi-terminal (Büttiker) view. With contacts and transmission probabilities , linear response reads
This enforces current conservation and automatically includes contact geometry.
10.12.2 Aharonov–Bohm (AB) and phase memory
In a ring, two paths enclose magnetic flux . Wavefunctions pick up relative phase with . Result:
- AB oscillations in conductance periodic in
- Altshuler–Aronov–Spivak (AAS) oscillations at from time-reversed loop pairs
Visibility decays with phase-coherence length set by dephasing (electron–electron, electron–phonon, external noise). Roughly, signals scale like along an interfering loop of length .
10.12.3 Shot noise: counting statistics under the hood
Even at , current fluctuates because charge is discrete. For a phase-coherent conductor,
where is the zero-frequency noise spectral density and the bias. Define the Fano factor
Key cases: ballistic plateau ⇒ (noise suppressed); diffusive wire ⇒ ; tunnel junction ⇒ Poissonian .
10.12.4 Weak localization (WL) and antilocalization (WAL)
Quantum interference of time-reversed diffusive paths enhances backscattering:
- With preserved time-reversal and weak spin–orbit, you get WL: a negative magnetoconductance dip at (conductance reduced by quantum correction).
- With strong SOC (or Berry phase), you get WAL: a positive cusp (antilocalization).
In 2D, the magnetoconductance near follows the Hikami–Larkin–Nagaoka style correction (schematically)
where and encodes SOC. In 2D without SOC, the correction scales like
with mean free path (sign flips for WAL).
10.12.5 Universal conductance fluctuations (UCF)
Mesoscopic conductance or exhibits aperiodic but reproducible wiggles from sample-specific interference. Their variance is universal (independent of material specifics), of order
where labels orthogonal (TRS), unitary (broken TRS), and symplectic (strong SOC) symmetry classes. Precise coefficients depend on dimensionality and aspect ratio, but the order of magnitude is fixed—hence “universal.”
10.12.6 Dephasing and energy relaxation
Phase coherence dies due to inelastic processes:
- Electron–electron scattering (dominant at low ): often with depending on dimension (sublinear in 1D, linear-ish in 2D diffusive)
- Electron–phonon at higher : typically with
- Magnetic impurities: spin flips are coherence killers; Kondo physics may show up in trends
Extract by fitting WL/WAL magnetoconductance or AB amplitude decay vs circumference.
10.12.7 Contact resistance and the four-probe cure
Landauer reminds us that even a perfect 1D wire has a contact-limited resistance
per occupied spin-degenerate channel. Multi-terminal geometries (Hall bars, four-probe) separate contact/interface effects from intrinsic device transmission.
10.12.8 Random-matrix universality
Fully chaotic/strongly disordered but phase-coherent cavities (“quantum dots”) have transmission eigenvalues whose statistics follow Wigner–Dyson ensembles:
- GOE (β=1) with TRS, no SOC
- GUE (β=2) broken TRS
- GSE (β=4) strong SOC, TRS preserved
This governs distributions of , shot noise (e.g., in chaotic cavities), and level spacings (Wigner surmise). Symmetry crossovers are tuned by magnetic field and SOC strength.
10.12.9 Mesoscopic superconductivity: Andreev in nanoscale
At a clean normal–superconductor (NS) interface, subgap transport occurs via Andreev reflection: an electron retroreflects as a hole, adding a Cooper pair to the condensate. Consequences:
- Doubling of low-bias conductance for a perfectly transparent NS point contact (BTK limit)
- Multiple Andreev reflections (MAR) in short SNS junctions, producing subharmonic gap features at
- In proximitized nanowires with SOC and Zeeman field, Andreev bound states can morph into Majorana-like zero modes at tuned parameters (topological superconductivity cameo)
10.12.10 From mesoscopic to Anderson: the scaling bridge
As disorder grows or dimension sinks, interference corrections (WL/WAL, UCF) strengthen. In low dimensions:
- 1D, 2D: arbitrarily weak disorder localizes noninteracting electrons at in the orthogonal class (no SOC, no magnetic field)
- 3D: a true metal–insulator transition can occur at critical disorder
The dimensionless conductance obeys a scaling flow ; the sign of decides metal vs insulator. We detail full Anderson localization and scaling theory in §10.13.
10.12.11 Worked mini-examples
(a) Landauer 101
A QPC shows plateaus at as the gate opens. Estimate the number of occupied modes and the average on a sloped transition between plateaus.
(b) Shot-noise triad
Compute for (i) a tunnel junction with , (ii) a ballistic QPC on a plateau, (iii) a diffusive wire. Explain why distinguishes them even if is similar.
(c) AB vs AAS
A ring of circumference shows oscillations with amplitude and a weaker harmonic . Argue why decays as .
(d) WL fit
Given a 2D magnetoconductance trace with a sharp cusp at and curvature reversal after adding a top gate (stronger SOC), fit qualitatively for and trends.
(e) UCF scaling
Show how thermal averaging over an energy window reduces variance when exceeds the Thouless energy .
(f) Contact resistance
A ballistic single-mode nanowire of length shows . Explain why adding a fourth probe pair reports nearly zero internal resistance while two-probe still reads .
10.12.12 Common pitfalls
- Forgetting contacts: Landauer conductance includes contact resistance; don’t blame the “perfect wire” for a finite two-probe .
- Mixing ensemble with sample noise: UCF are reproducible fingerprints vs or , not random flicker; average over windows to expose WL on top.
- WL sign confusion: WAL needs strong SOC or Berry phase (e.g., topological surfaces); otherwise expect WL.
- Assuming Fano always: Poissonian shot noise appears only in tunneling; ballistic and diffusive devices suppress noise predictably.
- Phase vs energy relaxation: dephasing kills interference without necessarily equilibrating energy; don’t equate with energy-relaxation time.
- Overfitting HLN: the 2D WL/WAL formula has parameter correlations; report trends (, ) with care, not over-precise numbers.
10.13 Anderson Localization & Scaling Theory
Turn disorder up and interference stops being cute: wavefunctions localize, DC transport dies, and conductance becomes exponentially small. That’s Anderson localization—a single-particle quantum phase where disorder + interference trap eigenstates. This section builds the Anderson model, scaling theory and the -function, mobility edges and criticality in 3D, symmetry-class twists (orthogonal/unitary/symplectic), finite-size and multifractals, experimental fingerprints (VRH, Coulomb gap, WL/WAL magnetoresponse), and quick routes for calculations (transfer matrix, Thouless number, nonlinear sigma model cameo). We end with a postcard from many-body localization.
10.13.1 The minimal model
The tight-binding Anderson Hamiltonian on a lattice:
with random on-site energies drawn from a distribution of width (e.g., box ) and short-range hopping . Compete vs bandwidth .
Heuristics. When elastic mean free path drops to the Ioffe–Regel brink , waves can’t complete a phase-coherent mean-free-path without scrambling—interference runs the show.
10.13.2 From diffusion to localization: the scaling hypothesis
Define the dimensionless conductance of a hypercubic sample of size :
with density of states , diffusion constant , Thouless energy , and level spacing . The one-parameter scaling postulate (Abrahams–Anderson–Licciardello–Ramakrishnan) says
depends only on (and symmetry class). Consequences:
- and 2 (orthogonal class): for all ⇒ no true metallic fixed point at ; arbitrarily weak disorder localizes (2D logarithmically slowly).
- : crosses zero at ⇒ metal–insulator transition (MIT). Correlation/localization length
diverges at critical disorder ; on the metallic side .
Symmetry tweaks (unitary/symplectic) modify and can change the 2D verdict (see §10.13.6).
10.13.3 Mobility edge and 3D criticality
In 3D, eigenstates can be extended near band center and localized in the tails; the energy separating them is the mobility edge . Tuning disorder moves ; at , hits and DC transport vanishes.
Critical behavior (orthogonal class, noninteracting):
- Localization length with universal
- Conductivity scaling on the metallic side
At and , with an exponent tied to by scaling.
10.13.4 Quantum corrections: WL/WAL and drifts
In the diffusive regime (), interference of time-reversed paths yields weak localization (WL) corrections. In 2D (orthogonal class), the zero-field correction behaves as
A perpendicular magnetic field breaks time-reversal and suppresses WL, producing a negative magnetoresistance cusp around (Hikami–Larkin–Nagaoka line shape). With strong spin–orbit coupling (symplectic class), the sign flips to weak antilocalization (WAL).
Interactions add Altshuler–Aronov corrections with trends and zero-bias anomalies in tunneling DOS. Together, WL+AA often drive the observed in 2D metals.
10.13.5 Insulating side: hopping and Coulomb gap
When states are localized with length , DC transport proceeds by phonon-assisted variable-range hopping (VRH).
- Mott VRH (no long-range Coulomb gap):
- Efros–Shklovskii VRH (Coulomb interactions carve a soft gap in DOS):
Magnetoresponse can be large and positive in VRH (wavefunction shrinkage and interference of hopping paths).
10.13.6 Symmetry classes and 2D twists
The fate of 2D is not one-size-fits-all; it depends on global symmetries:
- Orthogonal (TRS, no SOC): all states localized at ; slides down logarithmically.
- Unitary (TRS broken by ): cooperon suppressed; localization still wins in the end in 2D, but is less negative.
- Symplectic (strong SOC, TRS intact): WAL makes at large and allows a genuine 2D MIT tuned by disorder or carrier density.
Other classes (chiral, Bogoliubov–de Gennes) add special cases: e.g., integer quantum Hall plateaus (class A) with topological transitions between localized Hall plateaus.
10.13.7 Finite-size scaling and the transfer-matrix route
Numerics extract and by computing a quasi-1D localization length for wide bars of width using transfer matrices, then collapsing
Curves for different cross at ; the slope yields . This is the workhorse for noninteracting universality classes.
10.13.8 Multifractality at the mobility edge
Critical wavefunctions are neither plane waves nor exponentially localized—they are multifractal. Moments of amplitudes scale as
with a nontrivial spectrum of generalized dimensions (not just the geometric dimension ). Multifractality shows up in participation ratios, LDOS fluctuations, and critical conductance distributions.
10.13.9 Field theory cameo: nonlinear sigma model
Diagrammatics and disorder averaging lead to a nonlinear sigma model for diffusive modes (diffusons/cooperons). Coupling constants track ; renormalization reproduces , WL/WAL signs, and interaction corrections (Finkel’stein extension). This is the analytic backbone behind scaling and universality.
10.13.10 Optical and AC transport
On the metallic side, optical conductivity has a Drude peak with disorder-broadened width . Approaching criticality, Drude weight collapses. In an Anderson insulator (non-interacting), low-frequency AC conductivity rises as a power law (Mott ac hopping); interactions and Coulomb gap modify exponents. Finite-frequency probes thus access and hopping scales even when DC is immeasurably small.
10.13.11 Many-body localization (MBL): the sequel in interacting systems
Add interactions without a bath and localization can persist at finite energy density:
- Emergent local integrals of motion (LIOMs), area-law entanglement in highly excited eigenstates
- Logarithmic entanglement growth after quenches
- Absence of DC transport and failure of eigenstate thermalization
MBL requires isolation; phonons or leads eventually delocalize. Still, it’s a striking “Anderson-but-for-many-body” phase for cold atoms and mesoscopic circuits.
10.13.12 Worked mini-examples
(a) Thouless number
Show by comparing energy sensitivity of levels under twisted boundary conditions (Thouless) to the mean level spacing.
(b) 2D WL slope
Starting from the Cooperon loop integral, derive in 2D and explain why a perpendicular cuts off the log via .
(c) VRH exponents
From a constant DOS in dimensions, optimize hopping distance and energy to obtain the Mott VRH exponent . Then include a Coulomb gap to get the Efros–Shklovskii law.
(d) Transfer-matrix crossing
Sketch vs for widths and mark the single crossing at . Explain how the slope ratio gives .
(e) Multifractal moment
Given eigenstate amplitudes on a lattice, compute and extract an effective fractal dimension from a log–log fit across sizes.
(f) Symmetry flip
Argue qualitatively how adding strong SOC flips the sign of the Cooperon contribution and yields WAL, and why a small restores WL-like behavior.
10.13.13 Common pitfalls
- “All 2D localizes, period.” Only in the orthogonal class without strong SOC. Symplectic 2D can show a genuine MIT; QHE adds topological criticality.
- Confusing VRH with activated gaps. VRH exponents are non-Arrhenius; check the vs slope for dimension and Coulomb-gap effects.
- Equating with the transition. It’s a heuristic boundary of diffusive physics, not the precise .
- Ignoring interactions in 2D. Altshuler–Aronov terms are often as large as WL. If your slope flips with gate-tuned SOC, that’s WAL, not an interaction sign change.
- Over-reading finite samples. A rising does not prove insulating if ; do finite-size scaling.
- MBL ≠ Anderson. MBL is an interacting, highly excited-state phenomenon; Anderson localization is single-particle (or noninteracting) at .
10.13.14 Minimal problem kit
- Starting from Drude diffusion, build and the -function to first quantum-interference order in 2D; show why it’s negative in orthogonal and positive in symplectic at large
- Compute the magnetoconductance near using the Hikami–Larkin–Nagaoka formula; extract from a mock dataset
- Fit synthetic to Mott vs Efros–Shklovskii VRH and decide the presence of a Coulomb gap
- Implement a 3D Anderson model transfer-matrix code, produce crossings, and estimate by data collapse
- Evaluate for at and verify multifractal scaling across sizes
- Show that reduces to by inserting and
10.14 Moiré Quantum Matter & 2D Materials
Stack atomically thin crystals with a tiny twist and the band structure goes full kaleidoscope. Van der Waals (vdW) heterostructures let us dial geometry (twist), electrostatics (dual gates), and symmetry (alignment, strain) with near-video-game control. The result: flat bands, correlated insulators, unconventional superconductivity, quantum anomalous Hall, and moiré excitons—all living in a designer superlattice you can reconfigure on a chip.
10.14.1 Why 2D, why moiré
- 2D crystals (graphene, hBN, TMDs like MoS, WSe) can be peeled, stacked, rotated, and wired.
- Weak vdW bonding between layers allows arbitrary heterostructures without nasty interfacial chemistry.
- A small twist angle between layers creates a long-wavelength moiré superlattice with period
where is the monolayer lattice constant. For , is tens of nanometers—huge compared to —so momentum space folds into mini Brillouin zones and bands flatten.
10.14.2 Continuum picture and “magic” flat bands
For twisted bilayer graphene (tBLG), the Bistritzer–MacDonald continuum model couples Dirac cones from the two layers with interlayer tunneling and relative rotation . Define
and a dimensionless coupling
At magic angles (first one near ), pushes the Dirac velocity to nearly zero; bandwidth collapses into a few meV. With interactions of similar size, correlation runs the show.
A cardboard back-of-the-envelope for interaction scale:
which easily lands in the 10–50 meV ballpark at small , i.e., .
10.14.3 Moiré Hubbard models and Wigner physics
The moiré potential defines effective lattices for Wannier functions:
- tBLG: fragile topology complicates a strictly local tight-binding, but low-energy physics often projects to valley-resolved Chern bands near charge neutrality.
- TMD heterobilayers (e.g., WSe/WS with small twist): well-localized moiré triangular or honeycomb lattices; a bona fide Hubbard model with tunable by twist/field.
At low fillings and large , kinetic energy loses and Wigner crystals or generalized charge order can appear. Competing orders (Mott, charge density waves, spin/valley polarization) are tuned by gates and displacement field.
10.14.4 Topology in moiré: valleys, Berry curvature, Chern
Monolayer TMDs carry valley-contrasting Berry curvature; stacking and moiré folding can engineer Chern bands even without net magnetic field. Valley-resolved Hall responses:
Time-reversal demands total Chern if both valleys are equally filled, but spontaneous valley polarization breaks that balance, yielding quantum anomalous Hall in zero .
Berry curvature dipole in non-centrosymmetric stacks drives nonlinear Hall signals at second order in .
10.14.5 Broken symmetries: correlated insulators, SC, QAHE, nematics
Dial carrier density per moiré cell (via dual gates) and the phase diagram lights up:
- Correlated insulators at integer (Mott-like or flavor-polarized band insulators).
- Unconventional superconductivity upon slight doping of correlated states; sensitivity to displacement field and Coulomb screening suggests pairing from electronic mechanisms in some platforms.
- Quantum anomalous Hall effect (QAHE) when a single-valley Chern band is exchange-polarized.
- Nematicity: rotational symmetry breaking evident in transport and STM, often tied to strain or interaction-driven Pomeranchuk-like instabilities.
A minimalist effective interaction for pairing in a flat band with DOS :
Even modest gives sizable if is tiny; details hinge on valley/spin structure and form factors.
10.14.6 TMD excitons, interlayer excitons, and moiré polaritons
Monolayer TMDs host tightly bound excitons (hundreds of meV) thanks to reduced screening and 2D kinematics. Stacking two different TMDs forms type-II band alignment and interlayer excitons with electron and hole on different layers:
- Long lifetimes (ns–s), large electric dipoles, drift under in-plane fields.
- With a twist, the moiré potential creates exciton minibands and trapped arrays (an “exciton lattice”).
A hydrogenic estimate for a 2D exciton (with effective masses and dielectric constant) still helps:
Real numbers require the Keldysh screening model, but the scaling with and is the vibe.
Couple the excitons to a cavity: you get moiré exciton–polaritons with tunable dispersion and interactions.
10.14.7 Spin–valley locking and valleytronics
In TMD monolayers, SOC locks spin to valley. Circularly polarized light addresses vs selectively, enabling valley polarization. With broken inversion in bilayers (or applied displacement field), you can steer valley Hall currents transverse to without . The nonlinear Hall effect emerges from a finite Berry curvature dipole :
measurable as a second-harmonic Hall voltage at zero magnetic field.
10.14.8 Experimental knobs and probes
- Knobs: twist , displacement field (dual gates), carrier density per moiré cell, dielectric environment (hBN thickness, nearby metal gates), strain.
- Transport: longitudinal/Hall resistivity vs , Landau fans under , quantum oscillations reading off moiré FS pockets.
- STM/STS: direct DOS, gap maps, visualization of moiré patterns and local order.
- Optics: PL/reflectance for excitons, circular dichroism for valley, pump–probe for dynamics.
- Capacitance/compressibility: thermodynamic DOS and symmetry breaking (spin/valley polarization plateaus).
- Magnetometry: QAHE plateaus, orbital magnetization from Chern bands, SQUID-on-tip imaging of edge modes.
10.14.9 Minimal modeling routes
- Continuum Dirac + periodic tunneling for graphene family; captures flat bands and topology with a few parameters .
- Moiré Hubbard for TMD bilayers on triangular/honeycomb lattices with on-site , nearest-neighbor , and possible Hund/Ising anisotropies from SOC.
- Hartree–Fock for flavor symmetry breaking; DMFT/ED for correlation and finite- spectra; Berry + Wilson loops (see §10.11) to track Chern.
10.14.10 Worked mini-examples
(a) Moiré length scale
Graphene with and . Compute and estimate for .
(b) Magic-angle parameter
With and , compute at and comment on flatness.
(c) Triangular-lattice Hubbard
Assume , for a TMD moiré. Estimate with and discuss Mott vs metal at half filling.
(d) Valley Hall sign
If a single-valley Chern band with is fully polarized, what is in units of at zero Compare to the case where both valleys with opposite are filled.
(e) Interlayer exciton dipole
For layer spacing , estimate the dipole moment and the Stark shift for .
(f) Nonlinear Hall estimate
Given a Berry curvature dipole and relaxation time , estimate the order of magnitude of the second-harmonic Hall current under (scaling answer is fine).
10.14.11 Common pitfalls
- “Flat band = magic only.” Many moiré platforms (TMDs, graphene on hBN) host narrow bands without the exact magic-angle tuning; displacement field and dielectric screening reshape bandwidths dramatically.
- Ignoring flavor: spin, valley, and layer give multiple “flavors”; broken-symmetry physics often picks one (or a combo). Models that compress flavors too soon miss real orders.
- Topology blinders: total Chern can cancel across valleys; you still get valley Hall and orbital magnetization if flavors are imbalanced.
- Wannierization traps: fragile topology in tBLG obstructs strictly local Wannier bases unless you enlarge the subspace—don’t force a too-minimal TB.
- Overfitting mean-field: Hartree–Fock loves order; cross-check with thermodynamics, compressibility, and disorder robustness.
- Exciton hydrogenic overuse: 2D screening is nonlocal (Keldysh); hydrogenic estimates are scale guides, not final numbers.