The Central Idea
Landauโs Fermi liquid theory is one of the great conceptual triumphs of twentieth-century physics. It explains why the low-temperature properties of a strongly interacting Fermi system -- such as liquid He-3 or electrons in a metal -- can be described in terms of weakly interacting quasiparticles that stand in one-to-one correspondence with the bare particles of the non-interacting Fermi gas.
The key assumption is adiabatic continuity: as interactions are turned on from zero, the ground state and low-lying excited states of the free Fermi gas evolve smoothly into those of the interacting system. No phase transition occurs, and the quantum numbers labelling each state remain valid. The Fermi surface survives, though it may be deformed, and the quasiparticle lifetime diverges at the Fermi surface as $\tau_k \sim (E_k - E_F)^{-2}$, guaranteeing that quasiparticles are well-defined near $E_F$.
The power of the theory lies in the fact that all thermodynamic and transport properties at low temperature are determined by a small number of phenomenological parameters -- the Landau interaction function -- without requiring a microscopic calculation of the ground-state energy.
Landau Energy Functional
Landauโs starting point is to write the total energy of the system as a functional of the quasiparticle distribution function. Let $n_k$ denote the quasiparticle occupation number for state $\mathbf{k}\sigma$, and define the deviation from the ground-state distribution:
$$\delta n_k \;=\; n_k - n_k^{(0)}$$
where $n_k^{(0)} = \Theta(k_F - |\mathbf{k}|)$ is the ground-state Fermi-Dirac step function at $T=0$. The total energy is expanded to second order in $\delta n_k$:
$$E[\delta n] \;=\; E_0 + \sum_{\mathbf{k}\sigma} \epsilon_k^{(0)}\,\delta n_{\mathbf{k}\sigma} + \frac{1}{2V}\sum_{\mathbf{k}\sigma,\,\mathbf{k}'\sigma'} f_{\mathbf{k}\sigma,\,\mathbf{k}'\sigma'}\,\delta n_{\mathbf{k}\sigma}\,\delta n_{\mathbf{k}'\sigma'}$$
Here $E_0$ is the ground-state energy, $\epsilon_k^{(0)}$ is the bare quasiparticle energy measured from the chemical potential, and$f_{\mathbf{k}\sigma,\mathbf{k}'\sigma'}$ is the Landau interaction function -- a symmetric matrix describing how the energy of quasiparticle$\mathbf{k}\sigma$ is affected by the presence of a quasiparticle in state $\mathbf{k}'\sigma'$.
In a translationally invariant, isotropic system (such as liquid He-3), both$\mathbf{k}$ and $\mathbf{k}'$ lie on the Fermi surface, and $f$ depends only on the angle $\theta$ between them and the spin channel. We decompose:
$$f_{\mathbf{k}\sigma,\mathbf{k}'\sigma'} = f^s(\theta) + f^a(\theta)\,\boldsymbol{\sigma}\cdot\boldsymbol{\sigma}'$$
where $f^s$ and $f^a$ are the spin-symmetric and spin-antisymmetric parts. These are further expanded in Legendre polynomials:$f^{s,a}(\theta) = \sum_l f_l^{s,a}\,P_l(\cos\theta)$. The dimensionless Landau parameters are defined as$F_l^{s,a} = N^*(0)\,f_l^{s,a}$, where $N^*(0)$ is the density of states at the Fermi level (per unit volume, per spin).
Quasiparticle Energy
The energy of a single quasiparticle is obtained by the functional derivative of the total energy with respect to the occupation number. Since the interaction term is quadratic in $\delta n$, the quasiparticle energy itself depends on the distribution:
$$\tilde{\epsilon}_{\mathbf{k}\sigma} \;=\; \frac{\delta E}{\delta n_{\mathbf{k}\sigma}} \;=\; \epsilon_k^{(0)} + \frac{1}{V}\sum_{\mathbf{k}'\sigma'} f_{\mathbf{k}\sigma,\,\mathbf{k}'\sigma'}\,\delta n_{\mathbf{k}'\sigma'}$$
This is a crucial result: the quasiparticle energy is not a fixed single-particle spectrum but a functional of the distribution. When the system is driven out of equilibrium (e.g., by a temperature gradient or magnetic field), the quasiparticle energies shift self-consistently. This feedback mechanism is responsible for the renormalization of thermodynamic quantities.
At equilibrium and $T=0$, the quasiparticle energy reduces to$\tilde{\epsilon}_k = \epsilon_k^{(0)}$ since $\delta n = 0$. The Fermi surface is therefore determined by$\epsilon_{k_F}^{(0)} = \mu$, exactly as in the free gas.
Effective Mass
The quasiparticle effective mass $m^*$ is defined through the quasiparticle velocity at the Fermi surface:
$$v_F^* \;=\; \frac{1}{\hbar}\frac{\partial \tilde{\epsilon}_k}{\partial k}\bigg|_{k=k_F} \;=\; \frac{\hbar k_F}{m^*}$$
The relation between the effective mass and the Landau parameters is derived by considering a state with the Fermi sphere uniformly displaced by a small wavevector$\delta\mathbf{q}$. The current carried by this state can be computed in two ways: from the quasiparticle group velocities and from the physical particle current. Equating the two yields the fundamental relation:
$$\boxed{\frac{m^*}{m} \;=\; 1 + \frac{F_1^s}{3}}$$
This is a remarkable result. In liquid He-3, $F_1^s \approx 6$, giving$m^*/m \approx 3$. The effective mass enhancement directly controls the specific heat, since $C_V \propto m^*$ for a degenerate Fermi system.
Stability condition:
The effective mass must be positive for the Fermi liquid to be stable. This requires$F_1^s > -3$. When $F_1^s \to -3$, the effective mass diverges and the quasiparticle description breaks down, signalling a quantum phase transition (Pomeranchuk instability in the $l=1$ channel).
Galilean Invariance Constraint
For a system of particles with a well-defined mass $m$ (such as He-3 atoms), Galilean invariance imposes a powerful constraint relating the effective mass to the interaction function. Under a Galilean boost by velocity $\mathbf{v}$, the momentum of each particle changes by $m\mathbf{v}$. The current must transform as:
$$\mathbf{j} = \frac{1}{V}\sum_{\mathbf{k}\sigma} \frac{\hbar\mathbf{k}}{m}\,n_{\mathbf{k}\sigma} = \frac{1}{V}\sum_{\mathbf{k}\sigma} \mathbf{v}_{\mathbf{k}}^*\,n_{\mathbf{k}\sigma}$$
where the left-hand side uses the bare mass (since the physical current is carried by real particles of mass $m$) and the right-hand side uses quasiparticle velocities. This identity, combined with the self-consistent shift of quasiparticle energies, yields:
$$\frac{m}{m^*} = 1 - \frac{1}{V}\sum_{\mathbf{k}'\sigma'} f_{\mathbf{k}\sigma,\mathbf{k}'\sigma'}\,\frac{\partial n_{\mathbf{k}'\sigma'}^{(0)}}{\partial \epsilon_{\mathbf{k}'}}\,\frac{\hat{\mathbf{k}}\cdot\hat{\mathbf{k}}'}{d}\,\frac{\hbar^2 k_F}{m}$$
In three dimensions ($d=3$), evaluating the angular integral produces exactly the $m^*/m = 1 + F_1^s/3$ relation. This is not merely a definition of $m^*$ but a constraint: Galilean invariance fixes the relationship between the effective mass and the $l=1$Landau parameter. Any microscopic calculation must satisfy this sum rule.
Note that this constraint does not hold for electrons in a lattice, where Galilean invariance is broken by the periodic potential. In that case, the effective mass and the Landau parameters are independently determined.
Compressibility
The isothermal compressibility measures the response of the density to a change in pressure. In the Fermi liquid, we compute the change in chemical potential when the density is uniformly increased by $\delta n$. The quasiparticle energy shifts both because the Fermi energy changes and because of the interaction term:
$$\delta\mu = \frac{\delta n}{N^*(0)} + f_0^s\,\delta n = \frac{\delta n}{N^*(0)}\left(1 + F_0^s\right)$$
Since $\kappa^{-1} = n^2\,\partial\mu/\partial n$ and the free-gas compressibility is $\kappa_0^{-1} = n^2/N_0^*(0)$, we obtain:
$$\boxed{\frac{\kappa}{\kappa_0} = \frac{m^*/m}{1 + F_0^s}}$$
The numerator $m^*/m$ arises because the density of states is enhanced. The denominator $(1 + F_0^s)$ reflects the cost of adding particles against the mean-field repulsion.
Pomeranchuk instability ($l=0$, symmetric):
When $F_0^s \to -1$, the compressibility diverges, signalling a mechanical instability toward phase separation. The system spontaneously develops density inhomogeneities. Stability requires $F_0^s > -1$.
Spin Susceptibility
The spin susceptibility is computed by considering the response to an external magnetic field $\mathbf{B}$, which splits the quasiparticle energies by spin. A field along $\hat{z}$ produces$\delta n_{\uparrow} = -\delta n_{\downarrow} = \delta n$, with the magnetization $M = \mu_B(\delta n_{\uparrow} - \delta n_{\downarrow})$.
The self-consistent quasiparticle energy shift in the presence of the magnetization is:
$$\tilde{\epsilon}_{\uparrow} - \tilde{\epsilon}_{\downarrow} = 2\mu_B B + 2f_0^a\,\delta n$$
The factor of $f_0^a$ appears because the spin-antisymmetric interaction couples to the magnetization density. Solving self-consistently with$\delta n = N^*(0)\,\delta\epsilon/2$, we find:
$$\boxed{\frac{\chi}{\chi_0} = \frac{m^*/m}{1 + F_0^a}}$$
This has the same structure as the compressibility formula with $F_0^s$replaced by $F_0^a$. In liquid He-3, $F_0^a \approx -0.67$, leading to a large enhancement of the spin susceptibility:$\chi/\chi_0 \approx 3/(1 - 0.67) \approx 9$.
Stoner instability:
When $F_0^a \to -1$, the spin susceptibility diverges, indicating a ferromagnetic instability. The system spontaneously polarizes. This is the Landau analogue of the Stoner criterion. Stability requires $F_0^a > -1$. Liquid He-3 with $F_0^a \approx -0.67$ is fairly close to this instability.
Specific Heat Enhancement
The specific heat at low temperatures is determined by the density of states at the Fermi level, which in the interacting system involves the effective mass$m^*$ rather than the bare mass $m$:
$$N^*(0) = \frac{m^* k_F}{\pi^2 \hbar^2}$$
The Sommerfeld expansion for the specific heat remains valid in the Fermi liquid:
$$\boxed{C_V = \frac{\pi^2}{3}\,k_B^2\,N^*(0)\,T = \frac{m^*}{m}\,C_V^{(0)}}$$
The specific heat is enhanced by exactly the effective mass ratio $m^*/m$. Importantly, the specific heat depends only on $m^*$ (or equivalently $F_1^s$), not on $F_0^s$ or $F_0^a$. This is because the specific heat measures the number of thermally excited quasiparticles, which depends on the density of states but not on the interactions between them.
The Wilson ratio, defined as $R_W = (\pi^2 k_B^2/3\mu_B^2)(\chi/\gamma)$, provides a useful diagnostic. For the free Fermi gas $R_W = 1$, while for the Fermi liquid:
$$R_W = \frac{1}{1 + F_0^a}$$
General Stability Conditions
Landau showed that the Fermi liquid is stable against arbitrary deformations of the Fermi surface if and only if all the Pomeranchuk stability conditions are satisfied. For each angular momentum channel $l$ and spin channel $s/a$:
$$\boxed{F_l^{s} > -(2l+1), \qquad F_l^{a} > -(2l+1)}$$
These conditions ensure that the energy cost of deforming the Fermi surface in any channel is positive. The most important cases are:
| Channel | Condition | Instability |
|---|---|---|
| $l=0$, symmetric | $F_0^s > -1$ | Phase separation (collapse) |
| $l=0$, antisymmetric | $F_0^a > -1$ | Ferromagnetism (Stoner) |
| $l=1$, symmetric | $F_1^s > -3$ | Negative effective mass |
| $l=1$, antisymmetric | $F_1^a > -3$ | Spin-current instability |
| $l=2$, symmetric | $F_2^s > -5$ | Nematic (d-wave Pomeranchuk) |
The Pomeranchuk instabilities correspond to spontaneous deformations of the Fermi surface, breaking the rotational symmetry of the ground state. These have attracted renewed interest in the context of strongly correlated electron systems, where nematic order ($l=2$ instability) has been observed in several materials.
Summary of Fermi Liquid Relations
The following table collects the key results of Landau Fermi liquid theory, expressing the ratios of physical quantities to their free-gas values in terms of the dimensionless Landau parameters:
$$\frac{m^*}{m} = 1 + \frac{F_1^s}{3}$$
$$\frac{C_V}{C_V^0} = \frac{m^*}{m}$$
$$\frac{\kappa}{\kappa_0} = \frac{m^*/m}{1 + F_0^s}$$
$$\frac{\chi}{\chi_0} = \frac{m^*/m}{1 + F_0^a}$$
$$c_1^2 = \frac{v_F^2}{3}\left(1 + F_0^s\right) \cdot \frac{m}{m^*}$$
$$R_W = \frac{1}{1 + F_0^a}$$
Detailed Derivations
Derivation of the effective mass relation $m^*/m = 1 + F_1^s/3$
This fundamental relation follows from Galilean invariance. Consider a state where the entire Fermi sphere is displaced by a small wavevector $\delta\mathbf{q}$.
Step 1: The physical current density must equal $\mathbf{j} = n\hbar\delta\mathbf{q}/m$, since each of the $n$ particles carries momentum $\hbar\delta\mathbf{q}$ with bare mass $m$.
Step 2: In the quasiparticle picture, the current is carried by quasiparticles with velocity $\mathbf{v}^* = \hbar\mathbf{k}/m^*$. But the quasiparticle energy is shifted by the Landau interaction:
Step 3: For the displaced Fermi sphere, $\delta n_{\mathbf{k}\sigma} = n^0(|\mathbf{k} - \delta\mathbf{q}|) - n^0(k) \approx -\delta\mathbf{q}\cdot\hat{\mathbf{k}}\,\delta(k - k_F)$. The quasiparticle velocity acquires an interaction correction:
Step 4: The current carried by the quasiparticles is $\mathbf{j}_\text{qp} = \sum_{\mathbf{k}\sigma} \mathbf{v}^*_\mathbf{k}\,\delta n_{\mathbf{k}\sigma}/V$. Equating with the physical current $\mathbf{j} = n\hbar\delta\mathbf{q}/m$:
Step 5: In 3D, the angular average picks out the $l = 1$ Legendre component: $\langle \hat{\mathbf{k}}\cdot\hat{\mathbf{k}}'\rangle_{P_1} = \cos\theta/3$. Converting the sum to the density of states and using $F_1^s = N^*(0)\,f_1^s$:
Rearranging:
Derivation of the compressibility $\kappa/\kappa_0 = (m^*/m)/(1 + F_0^s)$
Step 1: Consider a uniform change in density $n \to n + \delta n$. This means every quasiparticle state acquires a shift $\delta n_{\mathbf{k}\sigma} = \delta n / (2N^*(0))$ (uniformly distributed, factor of 2 for both spins).
Step 2: The chemical potential shift has two contributions โ from the kinetic energy (density of states) and from the interaction:
Step 3: Using $F_0^s = N^*(0)\,f_0^s$ (where $N^*(0) = 2N^*(0)_\text{per spin}$ is the total DOS), and recalling the convention $\delta\mu = \delta n / N^*(0) + f_0^s \delta n$:
Step 4: The inverse compressibility is $\kappa^{-1} = n^2\,\partial\mu/\partial n = n^2(1 + F_0^s)/N^*(0)$. For the free gas, $\kappa_0^{-1} = n^2/N_0(0)$ where $N_0(0)$ uses the bare mass. Since $N^*(0)/N_0(0) = m^*/m$:
Derivation of the spin susceptibility $\chi/\chi_0 = (m^*/m)/(1 + F_0^a)$
Step 1: Apply a magnetic field $B$ along $\hat{z}$. The Zeeman energy shifts spin-up and spin-down quasiparticles in opposite directions: $\delta\epsilon_\uparrow = -\mu_B B$, $\delta\epsilon_\downarrow = +\mu_B B$.
Step 2: The occupation shifts are $\delta n_\uparrow = N^*(0)\,\delta\mu_\uparrow/2$ and $\delta n_\downarrow = -N^*(0)\,\delta\mu_\downarrow/2$, where $\delta\mu_\sigma$ includes the self-consistent Landau interaction. The magnetization is $\delta m = \delta n_\uparrow - \delta n_\downarrow$.
Step 3: The self-consistent equation for the spin splitting is:
The factor of $f_0^a$ (not $f_0^s$) appears because the spin-antisymmetric interaction couples opposite spins.
Step 4: Using $\delta m = N^*(0)(\tilde{\epsilon}_\uparrow - \tilde{\epsilon}_\downarrow)/2$ and substituting:
Step 5: Solving for $\delta m$:
Step 6: The susceptibility $\chi = \mu_0\mu_B\,\delta m / B$, and comparing with the free-gas Pauli susceptibility $\chi_0 = \mu_0\mu_B^2 N_0(0)$:
Derivation: Landau Energy Functional and Equilibrium Distribution
We begin from the total energy of the interacting Fermi system expressed as a functional of the quasiparticle distribution. Let $\delta n_{\mathbf{k}\sigma} = n_{\mathbf{k}\sigma} - n_{\mathbf{k}\sigma}^{(0)}$ denote the deviation of the quasiparticle occupation from the ground-state step function. The energy is expanded to second order:
$$E[\delta n] = E_0 + \sum_{\mathbf{k}\sigma} \epsilon_{\mathbf{k}}^0\,\delta n_{\mathbf{k}\sigma} + \frac{1}{2V}\sum_{\mathbf{k}\sigma,\,\mathbf{k}'\sigma'} f_{\mathbf{k}\sigma,\,\mathbf{k}'\sigma'}\,\delta n_{\mathbf{k}\sigma}\,\delta n_{\mathbf{k}'\sigma'}$$
Step 1 โ The Landau interaction function. The quantity $f_{\mathbf{k}\sigma,\,\mathbf{k}'\sigma'}$ is defined as the second functional derivative of the energy:
$$f_{\mathbf{k}\sigma,\,\mathbf{k}'\sigma'} = V\,\frac{\delta^2 E}{\delta n_{\mathbf{k}\sigma}\,\delta n_{\mathbf{k}'\sigma'}}$$
It is symmetric under exchange of $(\mathbf{k}\sigma) \leftrightarrow (\mathbf{k}'\sigma')$ and encodes all residual interactions between quasiparticles near the Fermi surface.
Step 2 โ Spin-channel decomposition. For a rotationally invariant system with spin-$\tfrac{1}{2}$ quasiparticles, the interaction function depends on the relative angle $\theta$ between $\mathbf{k}$ and $\mathbf{k}'$(both on the Fermi surface) and can be decomposed into spin-symmetric ($s$) and spin-antisymmetric ($a$) channels:
$$f_{\mathbf{k}\sigma,\,\mathbf{k}'\sigma'} = f^s(\theta) + f^a(\theta)\,\boldsymbol{\sigma}\cdot\boldsymbol{\sigma}'$$
Equivalently, for parallel spins $f_{\uparrow\uparrow} = f^s + f^a$ and for antiparallel spins $f_{\uparrow\downarrow} = f^s - f^a$. Each channel is expanded in Legendre polynomials:
$$f^{s,a}(\theta) = \sum_{\ell=0}^{\infty} f_\ell^{s,a}\,P_\ell(\cos\theta)$$
The dimensionless Landau parameters are $F_\ell^{s,a} = N^*(0)\,f_\ell^{s,a}$, where$N^*(0) = m^* k_F / (\pi^2\hbar^2)$ is the quasiparticle density of states at the Fermi surface (per spin).
Step 3 โ Equilibrium from free-energy minimization. At finite temperature, we minimize the Helmholtz free energy $F = E - TS$ with respect to $\delta n_{\mathbf{k}\sigma}$, subject to fixed total particle number. The entropy is the standard combinatorial expression:
$$S = -k_B \sum_{\mathbf{k}\sigma}\left[n_{\mathbf{k}\sigma}\ln n_{\mathbf{k}\sigma} + (1 - n_{\mathbf{k}\sigma})\ln(1 - n_{\mathbf{k}\sigma})\right]$$
Setting $\delta F / \delta n_{\mathbf{k}\sigma} = \mu$ (the chemical potential), we obtain:
$$\tilde{\epsilon}_{\mathbf{k}\sigma} + k_B T \ln\frac{n_{\mathbf{k}\sigma}}{1 - n_{\mathbf{k}\sigma}} = \mu$$
where $\tilde{\epsilon}_{\mathbf{k}\sigma} = \epsilon_{\mathbf{k}}^0 + V^{-1}\sum_{\mathbf{k}'\sigma'} f_{\mathbf{k}\sigma,\mathbf{k}'\sigma'}\,\delta n_{\mathbf{k}'\sigma'}$ is the self-consistent quasiparticle energy. Solving for the occupation:
$$\boxed{n_{\mathbf{k}\sigma} = \frac{1}{\exp\!\left[(\tilde{\epsilon}_{\mathbf{k}\sigma} - \mu)/(k_B T)\right] + 1}}$$
This is a Fermi-Dirac distribution, but with the self-consistent quasiparticle energy$\tilde{\epsilon}_{\mathbf{k}\sigma}$ replacing the bare energy. The distribution and the quasiparticle energies must be determined simultaneously, making this a nonlinear self-consistency problem. At $T = 0$, the solution reduces to a sharp Fermi step at $\tilde{\epsilon} = \mu$.
Derivation: Effective Mass from Galilean Invariance
We derive $m^*/m = 1 + F_1^s/3$ rigorously from Galilean invariance, connecting the quasiparticle velocity to the bare particle velocity through the $\ell = 1$ Landau parameter.
Step 1 โ Current from bare particles. For a system of $N$ particles each with bare mass $m$, the particle current density is exactly:
$$\mathbf{j} = \frac{1}{V}\sum_{\mathbf{k}\sigma} \frac{\hbar\mathbf{k}}{m}\,n_{\mathbf{k}\sigma}$$
This is an exact operator identity that holds regardless of interactions, because the current operator $\hat{\mathbf{j}} = \hat{\mathbf{P}}/(mV)$ involves only the total momentum and the bare mass.
Step 2 โ Current from quasiparticles. In the Landau framework, the same current must also be expressible in terms of quasiparticle group velocities:
$$\mathbf{j} = \frac{1}{V}\sum_{\mathbf{k}\sigma} \mathbf{v}_{\mathbf{k}\sigma}^*\,n_{\mathbf{k}\sigma}$$
where the quasiparticle velocity is $\mathbf{v}_{\mathbf{k}\sigma}^* = \nabla_{\mathbf{k}}\tilde{\epsilon}_{\mathbf{k}\sigma}/\hbar$. Since $\tilde{\epsilon}$ depends on $\delta n$, the velocity itself is distribution-dependent.
Step 3 โ Displaced Fermi sphere. Consider a state where the entire Fermi sea is shifted by a small wavevector $\delta\mathbf{q}$. Then:
$$\delta n_{\mathbf{k}\sigma} = n^{(0)}(|\mathbf{k} - \delta\mathbf{q}|) - n^{(0)}(k) \approx -(\delta\mathbf{q}\cdot\hat{\mathbf{k}})\,\delta(k - k_F)$$
This is a $P_1(\cos\theta)$ distortion of the Fermi surface. The bare-particle current for this state is simply $\mathbf{j} = n\hbar\delta\mathbf{q}/m$.
Step 4 โ Quasiparticle velocity correction. The gradient of the quasiparticle energy gives:
$$\hbar\mathbf{v}_{\mathbf{k}}^* = \frac{\hbar^2\mathbf{k}}{m^*} + \sum_{\mathbf{k}'\sigma'}\frac{f_{\mathbf{k}\sigma,\mathbf{k}'\sigma'}}{V}\,\nabla_{\mathbf{k}'}\!\left(\delta n_{\mathbf{k}'\sigma'}\right)\cdot\frac{\partial\mathbf{k}'}{\partial\mathbf{k}'}$$
For the displaced Fermi sphere, the interaction correction to the velocity along $\hat{\mathbf{k}}$at the Fermi surface becomes:
$$v_F^{*,\text{corr}} = \frac{1}{V}\sum_{\mathbf{k}'\sigma'} f_{\mathbf{k}\sigma,\mathbf{k}'\sigma'}\,(\hat{\mathbf{k}}\cdot\hat{\mathbf{k}}')\,\frac{\delta(k' - k_F)}{\hbar}$$
Step 5 โ Angular integration selects $\ell = 1$. Converting the sum over $\mathbf{k}'$ to an integral in 3D and using $\hat{\mathbf{k}}\cdot\hat{\mathbf{k}}' = \cos\theta = P_1(\cos\theta)$, the Legendre orthogonality relation picks out only the $\ell = 1$ component:
$$\int \frac{d\Omega'}{4\pi}\,f^s(\theta)\,\cos\theta = \frac{f_1^s}{3}$$
Step 6 โ Equating currents. The condition $\hbar\mathbf{k}/m = \mathbf{v}_{\mathbf{k}}^*$for every $\mathbf{k}$ on the Fermi surface gives:
$$\frac{\hbar k_F}{m} = \frac{\hbar k_F}{m^*} + \frac{\hbar k_F}{m^*}\cdot\frac{N^*(0)\,f_1^s}{3} = \frac{\hbar k_F}{m^*}\left(1 + \frac{F_1^s}{3}\right)$$
Rearranging immediately yields:
$$\boxed{\frac{m^*}{m} = 1 + \frac{F_1^s}{3}}$$
This is not merely a definition but a constraint imposed by Galilean symmetry: the $\ell = 1$ Landau parameter and the effective mass are not independent. Any microscopic theory of a Galilean-invariant Fermi liquid must satisfy this relation exactly. Note that for electrons in a crystal lattice, Galilean invariance is broken and this relation does not hold โ the band mass and Landau parameters become independent quantities.
Derivation: Compressibility and Spin Susceptibility Renormalization
Part A: Compressibility $\kappa/\kappa_0 = (m^*/m)/(1 + F_0^s)$
Step 1 โ Thermodynamic definition. The isothermal compressibility is defined as $\kappa = -(1/V)\partial V/\partial P|_T = (1/n^2)\partial n/\partial\mu|_T$. We compute $\partial\mu/\partial n$ for the Fermi liquid.
Step 2 โ Uniform density change. Increase the total density by $\delta n$. This corresponds to a spin-independent, angle-independent ($\ell = 0$) perturbation of the Fermi surface. Every quasiparticle state on the Fermi surface shifts outward uniformly, so$\delta n_{\mathbf{k}\sigma} = (\delta n / 2)\,\delta(\epsilon_k - \mu)/N^*(0)$(the factor of 2 accounts for both spins). The net change in density is:
$$\delta n = 2N^*(0)\,\delta\mu_\text{kin}$$
where $\delta\mu_\text{kin}$ is the shift in the kinetic part of the quasiparticle energy at the Fermi surface.
Step 3 โ Landau interaction correction. The total chemical potential shift includes the mean-field interaction with all other quasiparticles whose occupations have changed:
$$\delta\mu = \frac{\delta n}{2N^*(0)} + f_0^s\,\delta n$$
Only the $\ell = 0$ component $f_0^s$ survives because the perturbation is isotropic. Using $F_0^s = 2N^*(0)\,f_0^s$ (with the total density of states):
$$\delta\mu = \frac{\delta n}{2N^*(0)}\left(1 + F_0^s\right)$$
Step 4 โ Compressibility ratio. The inverse compressibility is$\kappa^{-1} = n^2\,\partial\mu/\partial n$. For the Fermi liquid:
$$\kappa_\text{FL}^{-1} = \frac{n^2}{2N^*(0)}(1 + F_0^s)$$
For the free Fermi gas, $\kappa_0^{-1} = n^2/(2N_0(0))$ where $N_0(0)$ uses the bare mass. Since $N^*(0)/N_0(0) = m^*/m$:
$$\boxed{\frac{\kappa}{\kappa_0} = \frac{m^*/m}{1 + F_0^s}}$$
The numerator $m^*/m$ enhances the compressibility (more states available at the Fermi surface), while $(1 + F_0^s)$ in the denominator suppresses it when $F_0^s > 0$ (repulsive interactions resist compression). When $F_0^s \to -1$, the compressibility diverges โ signalling a Pomeranchuk instability toward phase separation.
Part B: Spin Susceptibility $\chi/\chi_0 = (m^*/m)/(1 + F_0^a)$
Step 1 โ Applied magnetic field. Apply a uniform field$\mathbf{B} = B\hat{z}$. The Zeeman coupling shifts quasiparticle energies by $\mp\mu_B B$ for spin up/down. The resulting distribution change is purely in the spin-antisymmetric channel: $\delta n_\uparrow = -\delta n_\downarrow \equiv \delta m / 2$.
Step 2 โ Self-consistent energy shift. The quasiparticle spin splitting receives a Landau correction from the spin-antisymmetric interaction:
$$\Delta\tilde{\epsilon} \equiv \tilde{\epsilon}_\uparrow - \tilde{\epsilon}_\downarrow = -2\mu_B B + 2f_0^a\,\delta m$$
The key is that only $f_0^a$ (not $f_0^s$) enters, because the perturbation changes spin populations without changing total density.
Step 3 โ Self-consistency relation. The magnetization is determined by the density of states: $\delta m = N^*(0)\,|\Delta\tilde{\epsilon}|/2$. Substituting:
$$\delta m = N^*(0)\left(\mu_B B - f_0^a\,\delta m\right)$$
Step 4 โ Solving for $\delta m$. Collecting terms:
$$\delta m\left(1 + N^*(0)\,f_0^a\right) = N^*(0)\,\mu_B B$$
$$\delta m = \frac{N^*(0)\,\mu_B B}{1 + F_0^a}$$
Step 5 โ Susceptibility ratio. The Pauli susceptibility is $\chi = \mu_0\mu_B^2\,\delta m / (\mu_B B) = \mu_0\mu_B^2 N^*(0)/(1 + F_0^a)$. For the free gas, $\chi_0 = \mu_0\mu_B^2 N_0(0)$. Taking the ratio:
$$\boxed{\frac{\chi}{\chi_0} = \frac{N^*(0)/N_0(0)}{1 + F_0^a} = \frac{m^*/m}{1 + F_0^a}}$$
The structure is identical to the compressibility formula with $F_0^s \to F_0^a$. A negative$F_0^a$ enhances the susceptibility (the exchange interaction favors alignment). In liquid He-3,$F_0^a \approx -0.67$ and $m^*/m \approx 3$, giving$\chi/\chi_0 \approx 3/0.33 \approx 9$ โ a dramatic enhancement. When $F_0^a \to -1$, the susceptibility diverges, signalling a ferromagnetic (Stoner) instability.
Historical Context
Lev Davidovich Landau published his theory of the Fermi liquid in two landmark papers in 1956 and 1957 in the Soviet journal Zh. Eksp. Teor. Fiz. (translated in Soviet Physics JETP). Working at the Institute for Physical Problems in Moscow, Landau was motivated by the puzzle of liquid helium-3: unlike helium-4 (a boson which undergoes Bose-Einstein condensation), helium-3 atoms are fermions, and the liquid remained normal down to the lowest accessible temperatures of the time. The question was: how can a strongly interacting Fermi system be described theoretically?
Landau's genius was to realize that one did not need to solve the full many-body Schrรถdinger equation. Instead, the low-energy physics is captured by a phenomenological energy functional depending on a few Landau parameters $F_l^{s,a}$. These parameters encode all the effects of interactions and can be determined from experiment. The theory predicted the existence of zero sound, which was subsequently confirmed experimentally by Abel, Anderson, and Wheatley in 1966 in liquid He-3 at millikelvin temperatures.
The microscopic justification of Landau's phenomenological theory was provided by Abrikosov and Khalatnikov (1957) using diagrammatic perturbation theory, and by Noziรจres and Luttinger using the Green's function formalism. More recently, Shankar (1994) and Polchinski (1992) showed that Fermi liquid theory emerges as the stable low-energy fixed point of a Wilsonian renormalization group analysis, placing Landau's physical insights on the firmest possible theoretical footing. Landau received the Nobel Prize in Physics in 1962 for his contributions to the theory of condensed matter, particularly liquid helium.
Landau's 1956โ1958 Papers and Adiabatic Continuity
Landau's original 1956 paper (Zh. Eksp. Teor. Fiz. 30, 1058) introduced the concept of the quasiparticle energy functional and derived the thermodynamic relations. His 1957 follow-up paper extended the theory to transport phenomena and predicted zero sound โ a collisionless collective mode distinct from ordinary (hydrodynamic) sound. A third paper in 1958 treated the case of a charged Fermi liquid, incorporating the long-range Coulomb interaction and establishing the connection to the electron gas in metals.
The concept of adiabatic continuity is the cornerstone of the theory. Landau argued that if interactions are turned on infinitely slowly starting from the non-interacting Fermi gas, the eigenstates evolve one-to-one into the eigenstates of the interacting system, provided no phase transition (level crossing) occurs. Each non-interacting single-particle state $|\mathbf{k}\sigma\rangle$ maps onto a quasiparticle state with the same quantum numbers. The Fermi surface topology is preserved, and the quasiparticle residue $Z_k$ โ the overlap between the bare particle and the dressed quasiparticle โ remains nonzero ($0 < Z_{k_F} \leq 1$).
This adiabatic mapping is guaranteed to work in three dimensions (where the phase space for scattering vanishes as $(\epsilon - \epsilon_F)^2$), but breaks down in one dimension, where the Fermi liquid is replaced by a Luttinger liquid with qualitatively different properties (spin-charge separation, power-law correlations).
Migdal's Theorem
Arkady Migdal (1957) proved a key result underpinning the microscopic validity of Fermi liquid theory: the quasiparticle residue $Z_{k_F}$ is related to the discontinuity of the momentum distribution function $n(k)$ at the Fermi surface. Specifically, even in a strongly interacting system, the momentum distribution has a finite jump:
$$n(k_F^-) - n(k_F^+) = Z_{k_F} > 0$$
This discontinuity is the defining signature of a Fermi liquid. Migdal's theorem further establishes that vertex corrections to the electron-phonon coupling are suppressed by a factor of $m/M$ (electron-to-ion mass ratio), justifying the neglect of such corrections in the theory of conventional superconductivity. The theorem ensures that the quasiparticle picture remains valid even when the self-energy is large, as long as it is smooth near the Fermi surface.
Experimental Verification in Helium-3
The experimental confirmation of Landau's theory came primarily from measurements on liquid He-3 at millikelvin temperatures. Wheatley and collaborators at the University of California, San Diego, performed precision measurements of the specific heat, showing the linear-in-$T$ behavior with an enhanced Sommerfeld coefficient consistent with $m^*/m \approx 3$. The spin susceptibility was measured by NMR techniques, confirming the large enhancement predicted by the negative $F_0^a$ parameter.
The most dramatic confirmation was the observation of zero sound by Abel, Anderson, and Wheatley (1966). At temperatures below about 50 mK, where the quasiparticle mean free path exceeds the acoustic wavelength, ordinary (first) sound gives way to zero sound โ a collisionless mode propagating through coherent oscillations of the Fermi surface shape. The measured zero-sound velocity agreed quantitatively with the prediction from the known Landau parameters.
Later measurements by Greywall (1983) provided the most precise determination of the He-3 Landau parameters by combining specific heat, susceptibility, and sound velocity data. These parameters ($F_0^s \approx 10.07$, $F_1^s \approx 6.05$, $F_0^a \approx -0.67$) remain the benchmark values used today and provide a complete characterization of the normal state of liquid He-3 within the Landau framework.
Applications
Applications of Landau Fermi Liquid Theory
Advanced Derivation: Landau Energy Functional
Complete Structure of the Energy Functional
The Landau energy functional is the central object of Fermi liquid theory. We present here a detailed derivation that carefully distinguishes between bare quasiparticle energies, the interaction function, and the resulting renormalized quantities. The total energy of the interacting Fermi system, expanded to second order in deviations $\delta n_{\mathbf{k}\sigma}$ from the ground-state distribution, reads:
$$E[\delta n] = E_0 + \sum_{\mathbf{k}\sigma} \epsilon_{\mathbf{k}}^0\,\delta n_{\mathbf{k}\sigma} + \frac{1}{2V}\sum_{\mathbf{k}\sigma,\,\mathbf{k}'\sigma'} f_{\mathbf{k}\sigma,\,\mathbf{k}'\sigma'}\,\delta n_{\mathbf{k}\sigma}\,\delta n_{\mathbf{k}'\sigma'}$$
Here $E_0$ is the exact ground-state energy of the fully interacting system,$\epsilon_{\mathbf{k}}^0 = \delta E/\delta n_{\mathbf{k}\sigma}\big|_{n^{(0)}}$ is the bare quasiparticle energy (the first variational derivative evaluated at the ground state), and the Landau interaction function is the second variational derivative:
$$f_{\mathbf{k}\sigma,\,\mathbf{k}'\sigma'} = V\,\frac{\delta^2 E}{\delta n_{\mathbf{k}\sigma}\,\delta n_{\mathbf{k}'\sigma'}}\bigg|_{n^{(0)}}$$
The factor of $V$ (system volume) ensures that $f$ has dimensions of energy times volume and remains finite in the thermodynamic limit. The truncation at second order is exact for low-energy excitations because the density of excited quasiparticles scales as $T/\epsilon_F$ at low temperatures.
Definition of $f_{\mathbf{k}\mathbf{k}'}$ and Spin-Channel Decomposition
For a rotationally invariant system with spin-$\frac{1}{2}$ quasiparticles, the interaction function on the Fermi surface depends only on the angle $\theta$ between $\hat{\mathbf{k}}$ and$\hat{\mathbf{k}}'$. It can be decomposed into spin-symmetric and spin-antisymmetric channels:
$$f_{\mathbf{k}\sigma,\,\mathbf{k}'\sigma'} = f^s(\theta) + f^a(\theta)\,\boldsymbol{\sigma}\cdot\boldsymbol{\sigma}'$$
Equivalently, the parallel-spin and antiparallel-spin interaction functions are:
$$f^{\uparrow\uparrow}(\theta) = f^s(\theta) + f^a(\theta), \qquad f^{\uparrow\downarrow}(\theta) = f^s(\theta) - f^a(\theta)$$
Each channel is expanded in Legendre polynomials:
$$f^{s,a}(\theta) = \sum_{\ell=0}^{\infty} f_\ell^{s,a}\,P_\ell(\cos\theta)$$
The dimensionless Landau parameters are obtained by multiplying by the density of states at the Fermi surface:
$$F_\ell^{s,a} = N^*(0)\,f_\ell^{s,a}$$
where $N^*(0) = m^* k_F / (\pi^2 \hbar^2)$ is the renormalized density of states per unit volume (for both spin species). The $F_\ell^{s,a}$ are the physical quantities that enter all thermodynamic and transport relations. In practice, only $F_0^s$, $F_0^a$, and$F_1^s$ are needed for most observables, as higher harmonics are small in typical Fermi liquids.
Stability Conditions from the Energy Functional
Thermodynamic stability requires that the energy functional be a minimum at the ground state. This demands that all eigenvalues of the quadratic form remain positive, yielding the Pomeranchuk stability conditions for each angular momentum channel:
$$\boxed{F_\ell^s > -(2\ell + 1), \qquad F_\ell^a > -(2\ell + 1)}$$
Violation of any of these inequalities signals a phase transition: $F_0^s \to -1$ drives a density instability (phase separation), $F_0^a \to -1$ drives a ferromagnetic (Stoner) instability, and $F_1^s \to -3$ causes the effective mass to vanish, destabilizing the Fermi surface. These conditions impose rigorous bounds on the allowed values of the Landau parameters in any stable Fermi liquid.
Advanced Derivation: Effective Mass from Galilean Invariance
Rigorous Derivation of $m^*/m = 1 + F_1^s/3$
The relation between the effective mass and the $\ell = 1$ Landau parameter is not a model assumption but an exact consequence of Galilean invariance. We present here the complete derivation, emphasizing each logical step.
Step 1 โ Exact current operator. In a Galilean-invariant system of $N$ particles with bare mass $m$, the exact particle current density operator is:
$$\hat{\mathbf{j}} = \frac{\hat{\mathbf{P}}}{mV} = \frac{1}{mV}\sum_{i=1}^{N}\hat{\mathbf{p}}_i$$
This follows from $\hat{\mathbf{P}} = m\,d\hat{\mathbf{R}}_{\text{CM}}/dt$ and is an exact identity independent of interactions. In second quantization for a translationally invariant system:
$$\langle\hat{\mathbf{j}}\rangle = \frac{1}{V}\sum_{\mathbf{k}\sigma}\frac{\hbar\mathbf{k}}{m}\,\langle\hat{n}_{\mathbf{k}\sigma}\rangle$$
Step 2 โ Quasiparticle current. Within Landau theory, the same current must be computable from quasiparticle degrees of freedom. The quasiparticle group velocity is:
$$\mathbf{v}_{\mathbf{k}}^* = \frac{1}{\hbar}\nabla_{\mathbf{k}}\tilde{\epsilon}_{\mathbf{k}\sigma} = \frac{\hbar\mathbf{k}}{m^*} + \frac{1}{\hbar V}\sum_{\mathbf{k}'\sigma'} f_{\mathbf{k}\sigma,\mathbf{k}'\sigma'}\,\nabla_{\mathbf{k}'}\delta n_{\mathbf{k}'\sigma'}$$
The first term is the free quasiparticle velocity; the second is the backflow correction arising from the interaction-induced redistribution of quasiparticles. The quasiparticle current is:
$$\mathbf{j}_{\text{qp}} = \frac{1}{V}\sum_{\mathbf{k}\sigma}\mathbf{v}_{\mathbf{k}}^*\,\delta n_{\mathbf{k}\sigma}$$
Step 3 โ Galilean boost. Consider boosting the entire system by an infinitesimal velocity $\delta\mathbf{v}$, corresponding to shifting the Fermi sphere by$\delta\mathbf{q} = m\,\delta\mathbf{v}/\hbar$. The deviation from equilibrium is:
$$\delta n_{\mathbf{k}\sigma} = n^{(0)}(|\mathbf{k} - \delta\mathbf{q}|) - n^{(0)}(k) \approx -\hbar(\delta\mathbf{v}\cdot\hat{\mathbf{k}})\,\frac{m}{\hbar}\,\delta(k - k_F)$$
This is a purely $\ell = 1$ (dipolar) deformation. The bare-particle current for this state is exactly $\mathbf{j} = n\,\delta\mathbf{v}$, where $n = k_F^3/(3\pi^2)$ is the particle density.
Step 4 โ Interaction correction selects $F_1^s$. Substituting $\delta n$ into the quasiparticle velocity and performing the angular integration over $\hat{\mathbf{k}}'$, the orthogonality of Legendre polynomials ensures that only the $\ell = 1$ harmonic of$f^s(\theta)$ contributes (the antisymmetric channel $f^a$ drops out because the boost is spin-independent):
$$\int\frac{d\Omega'}{4\pi}\,f^s(\cos\theta)\,\cos\theta' = \frac{f_1^s}{3}\,\cos\theta$$
Step 5 โ Equating the two expressions for the current. Setting$\mathbf{j}_{\text{bare}} = \mathbf{j}_{\text{qp}}$ and using $N^*(0) = m^* k_F/(\pi^2\hbar^2)$:
$$\frac{\hbar k_F}{m} = \frac{\hbar k_F}{m^*} + \frac{\hbar k_F}{m^*}\cdot\frac{N^*(0) f_1^s}{3} = \frac{\hbar k_F}{m^*}\left(1 + \frac{F_1^s}{3}\right)$$
Canceling $\hbar k_F$ and inverting:
$$\boxed{\frac{m^*}{m} = 1 + \frac{F_1^s}{3}}$$
This result is a rigorous theorem for any Galilean-invariant Fermi liquid. It reflects the fact that when a quasiparticle moves through the medium, it drags along a backflow of surrounding quasiparticles, effectively increasing its inertia. The $F_1^s$ parameter quantifies this backflow. For liquid He-3,$F_1^s \approx 6.05$ gives $m^*/m \approx 3.02$, in excellent agreement with specific heat measurements. Crucially, this relation does not hold in crystalline solids where Galilean invariance is broken by the lattice โ there the band effective mass and Landau parameters are independent quantities.
Advanced Derivation: Compressibility and Spin Susceptibility
Compressibility: $\kappa/\kappa_0 = (m^*/m)/(1 + F_0^s)$
The isothermal compressibility is defined as $\kappa = -\frac{1}{V}\frac{\partial V}{\partial P}\big|_T = \frac{1}{n^2}\frac{\partial n}{\partial \mu}\big|_T$. To evaluate $\partial n/\partial \mu$ in the Landau framework, we consider a uniform change in the chemical potential $\mu \to \mu + \delta\mu$ and compute the resulting change in density.
Step 1 โ Self-consistency condition. When the chemical potential shifts by $\delta\mu$, the quasiparticle energy also shifts through the interaction term. A quasiparticle at the Fermi surface sees a renormalized energy:
$$\tilde{\epsilon}_{\mathbf{k}} = \epsilon_{\mathbf{k}}^0 + \frac{1}{V}\sum_{\mathbf{k}'\sigma'} f_{\mathbf{k}\sigma,\mathbf{k}'\sigma'}\,\delta n_{\mathbf{k}'\sigma'}$$
For a uniform compression (isotropic, spin-independent perturbation), $\delta n_{\mathbf{k}\sigma} = \delta(k - k_F)\,\delta k_F$ is independent of angle. The angular integration selects only the $\ell = 0$ component:
$$\delta\tilde{\epsilon} = \delta\epsilon^0 + f_0^s \cdot \frac{\delta n}{2} = \delta\epsilon^0 + \frac{F_0^s}{N^*(0)}\cdot\frac{\delta n}{2}$$
Step 2 โ Equilibrium condition. At equilibrium, the Fermi surface adjusts so that$\tilde{\epsilon}_{k_F} = \mu + \delta\mu$. The change in density is$\delta n = N^*(0)\,\delta\tilde{\epsilon}$, but the self-consistency requires:
$$\delta\mu = \delta\tilde{\epsilon} = \frac{\delta n}{N^*(0)} + F_0^s\,\frac{\delta n}{N^*(0)} = \frac{\delta n}{N^*(0)}\,(1 + F_0^s)$$
Therefore $\partial n/\partial\mu = N^*(0)/(1 + F_0^s)$. Using $N^*(0)/N_0(0) = m^*/m$ and$\kappa_0 = N_0(0)/n^2$:
$$\boxed{\frac{\kappa}{\kappa_0} = \frac{m^*/m}{1 + F_0^s}}$$
A repulsive interaction in the $\ell = 0$ symmetric channel ($F_0^s > 0$) stiffens the system against compression. Conversely, $F_0^s \to -1$ drives the compressibility to infinity โ the system becomes mechanically unstable and undergoes spinodal decomposition (phase separation).
Spin Susceptibility: $\chi/\chi_0 = (m^*/m)/(1 + F_0^a)$
The spin susceptibility derivation is the magnetic analogue of the compressibility calculation. We apply a uniform magnetic field $\mathbf{B} = B\hat{z}$ and compute the induced magnetization.
Step 1 โ Spin-dependent Fermi surface shift. The magnetic field splits the Fermi surfaces for spin-up and spin-down quasiparticles. The deviation from equilibrium is:
$$\delta n_{\mathbf{k}\uparrow} = +\frac{1}{2}N^*(0)\,\delta\epsilon_B, \qquad \delta n_{\mathbf{k}\downarrow} = -\frac{1}{2}N^*(0)\,\delta\epsilon_B$$
where $\delta\epsilon_B$ is the effective spin splitting to be determined self-consistently. The key point is that this perturbation is isotropic but spin-antisymmetric: $\delta n_\uparrow = -\delta n_\downarrow$. Consequently, only the spin-antisymmetric channel $f^a$ contributes.
Step 2 โ Self-consistency with exchange interaction. The quasiparticle energy shift for spin-up particles includes the bare Zeeman term plus the molecular field from the spin polarization:
$$\delta\tilde{\epsilon}_\uparrow = -\mu_B B + f_0^a \cdot N^*(0)\,\delta\epsilon_B = -\mu_B B + F_0^a\,\delta\epsilon_B$$
Self-consistency requires $\delta\epsilon_B = \delta\tilde{\epsilon}_\uparrow$, giving:
$$\delta\epsilon_B = \frac{-\mu_B B}{1 + F_0^a}$$
Step 3 โ Magnetization and susceptibility. The magnetization is $M = \mu_B(n_\uparrow - n_\downarrow) = \mu_B N^*(0)\,|\delta\epsilon_B|$. Taking the ratio with the free-gas result $\chi_0 = \mu_0\mu_B^2 N_0(0)$:
$$\boxed{\frac{\chi}{\chi_0} = \frac{N^*(0)/N_0(0)}{1 + F_0^a} = \frac{m^*/m}{1 + F_0^a}}$$
The physical interpretation is transparent: the numerator $m^*/m$ enhances the susceptibility through the increased density of states, while the denominator $(1 + F_0^a)$ represents the exchange molecular field. A negative $F_0^a$ (attractive spin-exchange interaction) enhances the susceptibility beyond the density-of-states effect. In liquid He-3, $F_0^a \approx -0.67$ and $m^*/m \approx 3$, giving $\chi/\chi_0 \approx 9$ โ a ninefold enhancement. When $F_0^a \to -1$, the susceptibility diverges, signaling a ferromagnetic Stoner instability where the system spontaneously spin-polarizes.
Extended Applications of Fermi Liquid Theory
Liquid Helium-3: The Canonical Fermi Liquid
Liquid He-3 remains the purest realization of Landau's theory. At temperatures below approximately 100 mK, the quasiparticle description becomes quantitatively accurate, with a mean free path scaling as $\ell \propto T^{-2}$ due to the phase-space restriction on quasiparticle scattering. The experimentally determined Landau parameters (Greywall 1983) are:
$$F_0^s \approx 10.07, \quad F_1^s \approx 6.05, \quad F_0^a \approx -0.67$$
These parameters yield: effective mass $m^*/m = 1 + F_1^s/3 \approx 3.02$, compressibility ratio$\kappa/\kappa_0 = 3.02/11.07 \approx 0.27$ (the liquid is much stiffer than the free gas due to strong repulsion), and spin susceptibility ratio $\chi/\chi_0 = 3.02/0.33 \approx 9.15$(dramatically enhanced by the nearly ferromagnetic exchange interaction).
The prediction and observation of zero sound โ a collisionless collective oscillation of the Fermi surface shape propagating at velocity $c_0 \approx v_F\sqrt{F_0^s/3}$ in the collisionless regime $\omega\tau \gg 1$ โ provided the most dramatic confirmation of the theory. The crossover from first sound (hydrodynamic, $\omega\tau \ll 1$) to zero sound occurs near 50 mK and was first observed by Abel, Anderson, and Wheatley (1966).
Heavy Fermion Systems: CeAl$_3$, UPt$_3$, and Beyond
Heavy fermion compounds represent Fermi liquid theory taken to its extreme. In these materials, localized $f$-electrons hybridize with conduction electrons through the Kondo effect, forming composite quasiparticles with effective masses $m^*/m \sim 100$โ$1000$.
CeAl$_3$ was one of the first identified heavy fermion metals (Andres, Graebner, Ott 1975), with a Sommerfeld coefficient$\gamma \approx 1600\;\text{mJ}\,\text{mol}^{-1}\text{K}^{-2}$ โ roughly 1000 times larger than in copper. The $T^2$ resistivity below about 0.3 K confirms Fermi liquid behavior, with the coefficient $A$ in $\rho = \rho_0 + AT^2$ satisfying the Kadowaki-Woods ratio $A/\gamma^2 \approx 1.0 \times 10^{-5}\;\mu\Omega\,\text{cm}\,(\text{mol}\,\text{K}/\text{mJ})^2$.
UPt$_3$ exhibits heavy Fermi liquid behavior with$\gamma \approx 450\;\text{mJ}\,\text{mol}^{-1}\text{K}^{-2}$ and undergoes unconventional superconductivity at $T_c \approx 0.53$ K with two distinct superconducting phases in the$H$โ$T$ phase diagram โ strong evidence for a multi-component order parameter arising from the heavy quasiparticles. De Haasโvan Alphen measurements directly map the heavy Fermi surface, confirming $m^*/m \approx 180$ on certain orbits.
The universal Kadowaki-Woods scaling $A/\gamma^2 \approx \text{const}$ across different heavy fermion compounds demonstrates that the scattering rate and specific heat both scale with the same quasiparticle density of states, as predicted by Fermi liquid theory. Deviations from this scaling (as in CeCu$_6$ near its quantum critical point) signal the breakdown of the Fermi liquid paradigm.
Normal State of Superconductors
The BCS theory of superconductivity is built upon the assumption that the normal state above$T_c$ is a well-defined Fermi liquid. The electron-phonon interaction renormalizes the quasiparticle effective mass through $m^*/m_b = 1 + \lambda$, where $\lambda$ is the electron-phonon coupling constant and $m_b$ is the band mass. Typical values are$\lambda \approx 0.43$ for aluminum, $\lambda \approx 1.55$ for lead, and$\lambda \approx 1.6$ for mercury.
The quasiparticle residue $Z = 1/(1 + \lambda)$ is directly measurable through tunneling spectroscopy: the tunneling density of states in the superconducting state reveals the energy-dependent self-energy $\Sigma(\omega)$, from which $\lambda$ and the spectral function$\alpha^2F(\omega)$ are extracted by McMillan-Rowell inversion. These experiments provide the most direct microscopic validation of the Fermi liquid picture in metals.
In the cuprate high-temperature superconductors, the situation is more complex. The optimally doped materials show a $T$-linear resistivity (the "strange metal" phase) that violates the Fermi liquid prediction of $T^2$ scattering, while the overdoped side recovers conventional Fermi liquid behavior. Understanding this crossover remains one of the central open problems in condensed matter physics.
Nuclear Matter
The interior of atomic nuclei and the outer core of neutron stars constitute dense Fermi liquids of nucleons. The nuclear Landau parameters are determined from nuclear structure data and nucleon-nucleon scattering phase shifts. A distinctive feature of nuclear matter is that the effective mass isreduced below the bare nucleon mass: $m^*/m \approx 0.7$โ$0.8$, implying a negative $F_1^s$.
The nuclear symmetry energy โ the energy cost per nucleon of converting protons to neutrons โ is directly related to the isospin-antisymmetric Landau parameter $F_0^a$ through:
$$E_{\text{sym}} = \frac{\epsilon_F}{3}\,\frac{m^*/m}{1 + F_0^a}$$
This quantity is a critical input for neutron star equation-of-state calculations and determines the proton fraction, cooling rate, and maximum mass of neutron stars. Giant dipole and quadrupole resonances in nuclei provide experimental access to the nuclear Landau parameters, connecting nuclear structure physics to the astrophysics of compact objects.
Two-Dimensional Electron Gas
Electrons confined at semiconductor interfaces (GaAs/AlGaAs heterostructures, Si-MOSFETs) form high-mobility 2D Fermi liquids. The two-dimensional analogue of Landau's theory replaces the Legendre expansion with a Fourier series in the azimuthal angle, and the effective mass relation becomes:
$$\frac{m^*}{m_b} = 1 + \frac{F_1^s}{2}$$
where the denominator $2$ (instead of $3$) reflects the two-dimensional geometry. At low electron densities (large $r_s = 1/(a_B^*\sqrt{\pi n})$), interactions become strong and the effective mass enhancement grows. Measurements by Tan et al. (2005) and Padmanabhan et al. (2008) show $m^*/m_b$ increasing by factors of 2โ3 as $r_s \to 8$โ$10$, approaching the Wigner crystallization threshold.
At high magnetic fields, the 2DEG enters the fractional quantum Hall regime where Fermi liquid theory breaks down entirely. However, Jain's composite fermion theory (1989) reintroduces a Fermi liquid description at even-denominator filling fractions (notably $\nu = 1/2$), where the composite fermions โ electrons bound to two flux quanta โ form their own Fermi sea with a well-defined Fermi surface observed in commensurability oscillations and surface acoustic wave experiments.
Expanded Historical Context
Landau's Program: 1956โ1958
Landau's theory of the Fermi liquid was developed in three landmark papers that progressively expanded the scope of the framework:
1956: The first paper (Zh. Eksp. Teor. Fiz. 30, 1058; JETP 3, 920) introduced the quasiparticle concept and the energy functional formalism. Landau derived the thermodynamic identities โ specific heat, compressibility, and susceptibility โ in terms of the interaction function. The physical motivation was liquid He-3, which had recently been liquefied and whose properties were anomalous compared to classical liquid theory.
1957: The second paper (Zh. Eksp. Teor. Fiz. 32, 59; JETP 5, 101) extended the theory to non-equilibrium phenomena and predicted the existence of zero sound โ a collisionless collective mode of the Fermi surface. This was a genuinely new prediction with no classical analogue: in contrast to ordinary sound, which requires frequent collisions to maintain local equilibrium, zero sound propagates through the coherent forward scattering of quasiparticles and exists only in the collisionless regime $\omega\tau \gg 1$.
1958: The third paper (Zh. Eksp. Teor. Fiz. 35, 97; JETP 8, 70) treated the charged Fermi liquid โ the electron gas in metals. Here Landau showed how the long-range Coulomb interaction modifies the formalism: the $\ell = 0$ Landau parameter acquires a Coulomb contribution that is screened at finite wavevector (Thomas-Fermi screening), while higher harmonics are essentially short-range. This paper established the connection between Fermi liquid theory and the established theory of metals by Sommerfeld and Bloch.
Adiabatic Continuity: The Foundation
The concept of adiabatic continuity is the logical cornerstone of Fermi liquid theory. The argument proceeds as follows: start with a non-interacting Fermi gas with a well-defined ground state (filled Fermi sea) and sharp quasiparticle excitations. Now imagine turning on the interactions infinitely slowly (adiabatically). If no phase transition (level crossing) occurs during this process, then:
(i) The ground state evolves continuously into the interacting ground state. (ii) Each single-particle excitation $|\mathbf{k}\sigma\rangle$ evolves into a quasiparticle excitation with the samequantum numbers ($\mathbf{k}$, $\sigma$). (iii) The one-to-one correspondence between free-particle states and quasiparticle states preserves the Fermi surface topology and Luttinger's theorem ($k_F$ is unchanged by interactions). (iv) The quasiparticle residue$Z_{k_F}$ โ the overlap between the bare electron and the dressed quasiparticle โ satisfies $0 < Z_{k_F} \leq 1$.
This mapping is guaranteed to work in three dimensions because the quasiparticle decay rate vanishes as$\text{Im}\,\Sigma \propto (\epsilon - \epsilon_F)^2$ near the Fermi surface โ the quasiparticles are long-lived enough to be well-defined. In one dimension, the phase space argument fails (the decay rate only vanishes linearly) and adiabatic continuity breaks down, leading instead to a Luttinger liquid. In two dimensions, the situation is marginal: logarithmic corrections appear, but the Fermi liquid generally survives (except at van Hove singularities or in the presence of nesting).
Migdal's Theorem and Its Consequences
Arkady Migdal (1957) proved that the momentum distribution function$n(\mathbf{k})$ of an interacting Fermi system retains a discontinuity at the Fermi surface:
$$n(k_F^-) - n(k_F^+) = Z_{k_F}$$
where $Z_{k_F}$ is the quasiparticle residue. This result has two profound implications. First, it provides a sharp, experimentally testable criterion for the existence of a Fermi liquid: if the momentum distribution (measurable via Compton scattering or angle-resolved photoemission) shows a discontinuity, quasiparticles exist. Second, it connects the phenomenological Landau theory to the microscopic Green's function formalism through:
$$Z_{k_F} = \left(1 - \frac{\partial\text{Re}\,\Sigma}{\partial\omega}\bigg|_{\omega=0}\right)^{-1}$$
Migdal's theorem also establishes that vertex corrections to the electron-phonon coupling are suppressed by $\mathcal{O}(\omega_D/\epsilon_F) \sim \mathcal{O}(\sqrt{m/M})$, where$\omega_D$ is the Debye frequency and $M$ the ion mass. This result is essential for the validity of Eliashberg theory of strong-coupling superconductivity and explains why the quasiparticle picture survives even when the electron-phonon coupling constant $\lambda$ is of order unity.
Key Experiments in Helium-3
The experimental verification of Landau's theory in liquid He-3 unfolded over three decades:
1960s โ Specific heat and susceptibility: Measurements by Wheatley and collaborators at UC San Diego confirmed the linear-in-$T$ specific heat with an enhanced Sommerfeld coefficient, and NMR experiments measured the spin susceptibility enhancement, both consistent with $m^*/m \approx 3$ and $F_0^a \approx -0.67$.
1966 โ Zero sound: Abel, Anderson, and Wheatley observed the transition from first sound to zero sound at frequencies above 15 MHz and temperatures below 50 mK. The zero-sound velocity exceeded the first-sound velocity by the predicted ratio, providing a uniquely quantum-mechanical collective phenomenon with no classical analogue.
1972 โ Superfluid He-3: The discovery of superfluidity in He-3 by Osheroff, Richardson, and Lee at Cornell (Nobel Prize 1996) demonstrated that the heavy Fermi liquid quasiparticles undergo $p$-wave (spin-triplet) Cooper pairing at 2.5 mK โ a direct consequence of the strong spin fluctuations encoded in the large negative $F_0^a$.
1983 โ Precision Landau parameters: Greywall's comprehensive measurements combined specific heat, susceptibility, and sound velocity data to extract the most precise values of the Landau parameters: $F_0^s = 10.07 \pm 0.15$,$F_1^s = 6.05 \pm 0.10$, $F_0^a = -0.67 \pm 0.01$. These values remain the definitive benchmark and demonstrate the internal consistency of Landau's framework at the sub-percent level.
Simulation: Landau Parameter Dependence
The following simulation plots the effective mass, compressibility, spin susceptibility, and specific heat enhancement as functions of the Landau parameters $F_0^s$,$F_0^a$, and $F_1^s$. Observe how the thermodynamic quantities diverge as the stability boundaries are approached.
Fermi Liquid: Thermodynamic Properties vs Landau Parameters
PythonPlots m*/m, compressibility, spin susceptibility, and specific heat as functions of Landau parameters F_0^s, F_0^a, and F_1^s. Stability boundaries are shown as red dashed lines.
Click Run to execute the Python code
Code will be executed with Python 3 on the server
Computation: Helium-3 Thermodynamics
Liquid He-3 is the paradigmatic Fermi liquid. Using the experimentally determined Landau parameters from Greywall and Wheatley, the Fortran program below computes the Fermi surface parameters, effective mass, specific heat, spin susceptibility, compressibility, and first-sound velocity of liquid He-3 at low temperatures.
Helium-3 Fermi Liquid Thermodynamics
FortranComputes thermodynamic properties of liquid He-3 using experimental Landau parameters: effective mass, Sommerfeld coefficient, spin susceptibility, compressibility, and sound velocity.
Click Run to execute the Fortran code
Code will be compiled with gfortran and executed on the server