Skip to content

Fermi-Liquid Response and Zero Sound

In a neutral or screened three-dimensional Fermi liquid, the same forward-scattering parameters that set the equilibrium compressibility and spin susceptibility also reshape the dynamical particle–hole response. A sufficiently repulsive density interaction pulls a collective pole above the particle–hole continuum: zero sound. The result is simple only after the route toward (q,ω)=(0,0)(q,\omega)=(0,0), the response convention, and the collision hierarchy have all been specified.

Required background. Use Landau theory for the quasiparticle interaction and total-density-of-states convention, and the Kubo-response framework for retarded correlators. Helpful background. The Ward-consistent response and diffusion, conductivity, and susceptibility pages explain why conservation laws distinguish the limits below.

Four paths to the origin give four observables

Section titled “Four paths to the origin give four observables”

Assume an isotropic, spin-rotation-invariant normal Fermi liquid at low temperature. Let N(0)N(0) be the total two-spin quasiparticle density of states at the Fermi surface, vF∗v_F^* the quasiparticle Fermi velocity, and τcoll\tau_{\mathrm{coll}} the relaxation time of the nonconserved angular distortions relevant to the longitudinal mode. The small-qq, small-ω\omega plane contains several physically different experiments:

ProbeLimiting path or hierarchyLeading density result
Equilibrium compressionω→0\omega\to0 before q→0q\to0Rs=N(0)/(1+F0s)\mathcal R_s=N(0)/(1+F_0^s)
Spatially uniform driveq→0q\to0 at nonzero ω\omegaRs=0\mathcal R_s=0 by number conservation
Collisionless rayq,ω→0q,\omega\to0 at fixed z=(ω+i0+)/(qvF∗)z=(\omega+i0^+)/(qv_F^*) and ωτcoll≫1\omega\tau_{\mathrm{coll}}\gg1Particle–hole continuum and possible zero-sound pole
Hydrodynamic rayq,ω→0q,\omega\to0 with ωτcoll≪1\omega\tau_{\mathrm{coll}}\ll1First sound in a clean neutral fluid

Thus zero sound is neither of the two iterated limits. At any fixed nonzero temperature, taking q→0q\to0 along ω∝q\omega\propto q eventually makes ωτcoll≪1\omega\tau_{\mathrm{coll}}\ll1; the strict infrared mode is then hydrodynamic. The collisionless limit instead keeps

q≪kF,τcoll−1≪ω=sqvF∗≪EF,qℓcoll≫1,q\ll k_F, \qquad \tau_{\mathrm{coll}}^{-1}\ll\omega=sqv_F^*\ll E_F, \qquad q\ell_{\mathrm{coll}}\gg1,

with ℓcoll=vF∗τcoll\ell_{\mathrm{coll}}=v_F^*\tau_{\mathrm{coll}}. One may reach this window by lowering the temperature before taking the long-wavelength limit.

The two panels below separate collisionless spectral kinematics from the paths used to approach the origin. In panel (a), inspect the narrow gap between the one-pair edge and the F0s=1F_0^s=1 zero-sound ray. In panel (b), notice that at fixed finite τcoll\tau_{\mathrm{coll}} the collisionless route stops before the crossover, whereas the hydrodynamic route reaches the origin.

At positive momentum and frequency, the particle–hole continuum fills the wedge below its linear edge and a collisionless zero-sound branch lies above it; a second panel distinguishes the static, uniform, fixed-z collisionless, and hydrodynamic first-sound paths toward the same origin.

Small-qq response of a neutral or statically screened, isotropic three-dimensional Fermi liquid. The hatched region in panel (a) is the leading collisionless one-particle–one-hole continuum, 0<ω<qvF∗0<\omega<qv_F^*; the displayed F0s=1F_0^s=1 zero-sound root is s0=1.044382s_0=1.044382 and lies just above its edge. More generally, the F0sF_0^s-only kernel has a separated real root s0>1s_0>1 for F0s>0F_0^s>0 in the collisionless hierarchy ωτcoll≫1\omega\tau_{\mathrm{coll}}\gg1. Panel (b) shows why static, uniform, collisionless, and hydrodynamic approaches to (q,ω)=(0,0)(q,\omega)=(0,0) are inequivalent. First sound belongs to ωτcoll≪1\omega\tau_{\mathrm{coll}}\ll1 and is not a pole of the collisionless kernel; its possible location inside the collisionless one-pair wedge does not by itself imply collisionless Landau damping. The figure is schematic apart from the equation-checked slope ratio in panel (a), and it omits recoil, thermal and multipair broadening, lattice anisotropy, and the unscreened Coulomb plasmon.

Hatching denotes continuum support, the heavy solid ray denotes zero sound, dash–dot denotes first sound, and the pale dotted band denotes a collision crossover rather than a phase boundary. The downloadable SVG and machine-readable semantic record preserve the visual and scientific distinctions.

For orientation, the chapter’s validity and failure map places this response calculation among its microscopic inputs and instabilities, while the claim table states the assumptions behind the principal Fermi-liquid formulas.

Define density and spin-number channels by

ns=n↑+n↓,na=n↑−n↓,n_s=n_\uparrow+n_\downarrow, \qquad n_a=n_\uparrow-n_\downarrow,

and couple sources through

HU=H0−∫d3x [Usns+Uana].H_U=H_0-\int d^3x\,[U_s n_s+U_a n_a].

With the site’s convention

GcR(t,x)=−iθ(t)⟨[nc(t,x),nc(0)]⟩,c=s,a,G_c^R(t,\mathbf x) =-i\theta(t)\langle[n_c(t,\mathbf x),n_c(0)]\rangle, \qquad c=s,a,

the physical response is Rc=δnc/δUc=−GcR\mathcal R_c=\delta n_c/\delta U_c=-G_c^R. This minus sign is worth making explicit: the static physical response is positive. In the total-DOS convention,

lim⁡q→0lim⁡ω→0Rs(q,ω)=N(0)1+F0s=n2κT,\lim_{q\to0}\lim_{\omega\to0}\mathcal R_s(q,\omega) =\frac{N(0)}{1+F_0^s} =n^2\kappa_T, lim⁡q→0lim⁡ω→0Ra(q,ω)=N(0)1+F0a.\lim_{q\to0}\lim_{\omega\to0}\mathcal R_a(q,\omega) =\frac{N(0)}{1+F_0^a}.

If Ua=μmagBU_a=\mu_{\mathrm{mag}}B and M=μmagnaM=\mu_{\mathrm{mag}}n_a, then χM=μmag2Ra\chi_M=\mu_{\mathrm{mag}}^2\mathcal R_a. A convention using Sz=na/2S_z=n_a/2 instead has the corresponding factor of 1/41/4. The divergences at F0s=−1F_0^s=-1 and F0a=−1F_0^a=-1 mark density and ferromagnetic Pomeranchuk instabilities, respectively. These thermodynamic relations follow from varying the Landau energy functional; see Dupuis 2025, § 4.2.2–4.2.3, pp. 296–297.

The kinetic equation transports a surface deformation

Section titled “The kinetic equation transports a surface deformation”

Choose q\mathbf q as the polar axis, write μ=p^⋅q^\mu=\hat{\mathbf p}\cdot\hat{\mathbf q}, and parametrize a Fermi-surface displacement by

δnpσ=−∂n0∂ϵp[νs(μ)+σνa(μ)]ei(q⋅r−ωt),σ=±1.\delta n_{\mathbf p\sigma} =-\frac{\partial n_0}{\partial\epsilon_{\mathbf p}} \bigl[\nu^s(\mu)+\sigma\nu^a(\mu)\bigr] e^{i(\mathbf q\cdot\mathbf r-\omega t)}, \qquad \sigma=\pm1.

The dimensionless Landau functions are

Fc(p^⋅p^′)=N(0)fc(p^⋅p^′)=∑ℓ=0∞FℓcPℓ(p^⋅p^′),c=s,a.F^c(\hat{\mathbf p}\cdot\hat{\mathbf p}') =N(0)f^c(\hat{\mathbf p}\cdot\hat{\mathbf p}') =\sum_{\ell=0}^{\infty}F_\ell^c P_\ell(\hat{\mathbf p}\cdot\hat{\mathbf p}'), \qquad c=s,a.

Keeping only F0cF_0^c is a solvable truncation, not a claim that higher harmonics vanish in a real liquid. In this model the quasiparticle-energy shift is

δϵpσ=F0sνˉs+σF0aνˉa−Us−σUa,νˉc=12∫−11νc(μ) dμ.\delta\epsilon_{\mathbf p\sigma} =F_0^s\bar\nu^s+\sigma F_0^a\bar\nu^a -U_s-\sigma U_a, \qquad \bar\nu^c=\frac12\int_{-1}^{1}\nu^c(\mu)\,d\mu.

The linearized Landau–Boltzmann equation is

(∂t+vp∗⋅∇r)δnpσ−∂n0∂ϵp vp∗⋅∇rδϵpσ=Icoll[δn].\left( \partial_t+\mathbf v_{\mathbf p}^*\cdot\nabla_{\mathbf r} \right)\delta n_{\mathbf p\sigma} -\frac{\partial n_0}{\partial\epsilon_{\mathbf p}}\, \mathbf v_{\mathbf p}^*\cdot\nabla_{\mathbf r} \delta\epsilon_{\mathbf p\sigma} =I_{\mathrm{coll}}[\delta n].

Dropping its collision integral in the hierarchy above reduces it to one equation per channel,

(z−μ)νc(μ)=μ[F0cνˉc−Uc],z=ω+i0+qvF∗.(z-\mu)\nu^c(\mu) =\mu\bigl[F_0^c\bar\nu^c-U_c\bigr], \qquad z=\frac{\omega+i0^+}{qv_F^*}.

This equation says something geometrically simple: quasiparticles with velocity component μvF∗\mu v_F^* stream the angular distortion, while the mean field F0cνˉcF_0^c\bar\nu^c tries to make the whole Fermi surface move coherently. Landau’s original derivation proceeds from the full kinetic equation to this angular eigenproblem in Landau 1957, printed pp. 102–103, Eqs. (4)–(14).

The angular integral exposes the continuum and the pole

Section titled “The angular integral exposes the continuum and the pole”

Solving for νc(μ)\nu^c(\mu) and averaging once gives

νˉc=Φ(z)[F0cνˉc−Uc],\bar\nu^c =\Phi(z)\bigl[F_0^c\bar\nu^c-U_c\bigr],

where

Φ(z)=12∫−11μ dμz−μ=z2log⁡(z+1z−1)−1.\Phi(z) =\frac12\int_{-1}^{1}\frac{\mu\,d\mu}{z-\mu} =\frac{z}{2}\log\left(\frac{z+1}{z-1}\right)-1.

The retarded logarithm is analytic for Im⁡z>0\operatorname{Im}z>0 and is real for z>1z>1. Since δnc=N(0)νˉc\delta n_c=N(0)\bar\nu^c, the correlator and physical response are

GcR(q,ω)=N(0)Φ(z)1−F0cΦ(z),Rc(q,ω)=−N(0)Φ(z)1−F0cΦ(z).G_c^R(q,\omega) =N(0)\frac{\Phi(z)}{1-F_0^c\Phi(z)}, \qquad \mathcal R_c(q,\omega) =-N(0)\frac{\Phi(z)}{1-F_0^c\Phi(z)}.

Two immediate checks recover the first table: Φ(0)=−1\Phi(0)=-1 gives the static susceptibility, whereas Φ(z→∞)→0\Phi(z\to\infty)\to0 makes the uniform finite-frequency density response vanish.

For 0<s<10<s<1, the upper rim of the retarded branch cut is

Φ(s+i0+)=s2log⁡(1+s1−s)−1−iπs2.\Phi(s+i0^+) =\frac{s}{2}\log\left(\frac{1+s}{1-s}\right) -1-\frac{i\pi s}{2}.

This cut is the linearized particle–hole continuum, 0<ω<qvF∗0<\omega<qv_F^* at positive frequency. The leading edge qvF∗qv_F^* is universal at small qq, but its O(q2)O(q^2) correction depends on the curvature of the quasiparticle dispersion; for a strictly quadratic dispersion it is q2/(2m∗)q^2/(2m^*). The Lindhard-response page displays this recoil explicitly. At T=0T=0 and positive frequency, use the same spectral normalization as that page:

Sc(q,ω)=−2Im⁡GcR(q,ω)=2Im⁡Rc(q,ω).S_c(q,\omega) =-2\operatorname{Im}G_c^R(q,\omega) =2\operatorname{Im}\mathcal R_c(q,\omega).

The continuum weight is manifestly positive:

Sccont(q,ω)=N(0)πs[1−F0cRe⁡Φ(s)]2+[πF0cs/2]2,0<s<1.S_c^{\mathrm{cont}}(q,\omega) =\frac{N(0)\pi s} { [1-F_0^c\operatorname{Re}\Phi(s)]^2 +[\pi F_0^c s/2]^2 }, \qquad 0<s<1.

At fixed qq, a density probe therefore sees a broad particle–hole continuum plus, when it is separated, the sharp collective peak derived next. At finite temperature, −2Im⁡GcR-2\operatorname{Im}G_c^R remains the spectral density, but the unsymmetrized structure factor contains detailed balance:

Sc>(q,ω)=−2Im⁡GcR(q,ω)1−e−βω.\mathcal S_c^>(q,\omega) =\frac{-2\operatorname{Im}G_c^R(q,\omega)} {1-e^{-\beta\omega}}.

The response function, continuum, and collective contribution are derived together in Dupuis 2025, § 4.3.4.1, printed pp. 305–307, Eqs. (4.100)–(4.108).

The zero-sound pole has a computable residue

Section titled “The zero-sound pole has a computable residue”

A source-free density oscillation requires the denominator to vanish:

1=F0sΦ(s0),s0>1,ω0(q)=s0qvF∗.1=F_0^s\Phi(s_0), \qquad s_0>1, \qquad \omega_0(q)=s_0qv_F^*.

For every F0s>0F_0^s>0, Φ(s)\Phi(s) decreases monotonically from +∞+\infty at s=1+s=1^+ to 00 at infinity, so there is one real root above the one-pair continuum. It is undamped only within the collisionless, linearized, single-pair theory. Collisions, thermally excited quasiparticles, multipair states, and finite-qq curvature can all broaden or obscure it.

The pole’s weight is not determined by its speed. Since

Φ′(s)=12log⁡(s+1s−1)−ss2−1<0(s>1),\Phi'(s) =\frac12\log\left(\frac{s+1}{s-1}\right) -\frac{s}{s^2-1}<0 \qquad (s>1),

the pole part of the correlator is

Gs,poleR(q,ω)=Z0(q)[1ω−ω0+i0+−1ω+ω0+i0+],G_{s,\mathrm{pole}}^R(q,\omega) =\mathcal Z_0(q) \left[ \frac1{\omega-\omega_0+i0^+} -\frac1{\omega+\omega_0+i0^+} \right],

with positive collective-response residue

Z0(q)=−N(0)qvF∗(F0s)2Φ′(s0)>0.\mathcal Z_0(q) =-\frac{N(0)qv_F^*} {(F_0^s)^2\Phi'(s_0)}>0.

This derivative formula is the long-wavelength form of Dupuis 2025, printed pp. 306–307, Eq. (4.108).

This Z0\mathcal Z_0 is not the one-particle quasiparticle residue conventionally called ZZ. At T=0T=0 and positive frequency the pole contributes

Sspole(q,ω)=2πZ0(q)δ(ω−ω0).S_s^{\mathrm{pole}}(q,\omega) =2\pi\mathcal Z_0(q)\delta(\omega-\omega_0).

For example, F0s=1F_0^s=1 gives

s0=1.044382…,Z0(q)N(0)qvF∗=0.104217….s_0=1.044382\ldots, \qquad \frac{\mathcal Z_0(q)}{N(0)qv_F^*} =0.104217\ldots.

The limiting forms provide strong diagnostic checks:

s0−1=2exp⁡[−2(1+1F0s)][1+o(1)],0<F0s≪1,s_0-1 =2\exp\left[-2\left(1+\frac1{F_0^s}\right)\right] [1+o(1)], \qquad 0<F_0^s\ll1, s02=F0s3+35+O((F0s)−1),F0s≫1.s_0^2 =\frac{F_0^s}{3}+\frac35 +O\bigl((F_0^s)^{-1}\bigr), \qquad F_0^s\gg1.

At weak coupling the pole exists but is exponentially close to the continuum and has exponentially small weight. Resolving it requires the separation (s0−1)qvF∗(s_0-1)qv_F^* to exceed collision broadening, thermal smearing, and the recoil correction. Landau 1957, printed pp. 103–104 and Dupuis 2025, p. 302, Eqs. (4.83)–(4.84) give the pole equation and leading weak- and strong-coupling limits; the subleading 3/53/5 above follows by expanding the same equation one order further.

For a stable attractive channel, −1<F0s≤0-1<F_0^s\le0, there is no separated real pole above the continuum in this truncation. Any pole obtained by continuing through the branch cut lies on an unphysical sheet and is damped; it must not be described as a real undamped mode “inside” the continuum. At F0s<−1F_0^s<-1, an upper-half-plane mode instead diagnoses instability of the assumed normal state. See Dupuis 2025, § 4.3.3.1, printed p. 302.

Spin modes and higher harmonics change the eigenproblem

Section titled “Spin modes and higher harmonics change the eigenproblem”

In an exactly spin-rotation-invariant liquid, the same derivation applies to the spin-number response after F0s→F0aF_0^s\to F_0^a. An undamped F0F_0-only spin-zero-sound pole therefore requires F0a>0F_0^a>0. The common range −1<F0a<0-1<F_0^a<0 enhances the static spin response but leaves the collective feature in the particle–hole continuum; F0a→−1+F_0^a\to-1^+ approaches a ferromagnetic instability. Spin–orbit, dipolar, and other spin-relaxation processes add decay, and the hydrodynamic spin mode is generally diffusive rather than ordinary first sound. Landau gives the corresponding spin kinetic equation in Landau 1957, printed p. 107, Eqs. (25)–(29).

When Fℓ>0cF_{\ell>0}^c is retained, νc(p^)\nu^c(\hat{\mathbf p}) must be expanded in angular harmonics and the kinetic equation becomes a matrix eigenvalue problem. Higher harmonics can shift the mode speed, redistribute spectral weight, and create additional angular collective modes. The one-line F0F_0 pole equation is then no longer the complete dispersion relation.

Conservation requires the matching current vertex

Section titled “Conservation requires the matching current vertex”

At large zz, the F0F_0 response has

Φ(z)=13z2+15z4+O(z−6),\Phi(z)=\frac1{3z^2}+\frac1{5z^4}+O(z^{-6}),

so

GsR(q,ω)=N(0)(vF∗)2q23ω2+O(q4ω4).G_s^R(q,\omega) =\frac{N(0)(v_F^*)^2q^2}{3\omega^2} +O\left(\frac{q^4}{\omega^4}\right).

The explicit q2q^2 verifies that a spatially uniform scalar drive cannot change a conserved total particle number. For a spherical Fermi surface, however, the displayed coefficient is n/m∗n/m^*. A Galilean-invariant system must instead have n/mn/m and obey

∫0∞dω2π ωSs(q,ω)=nq22m.\int_0^\infty \frac{d\omega}{2\pi}\, \omega S_s(q,\omega) =\frac{nq^2}{2m}.

The missing factor is supplied by the F1sF_1^s current vertex and quasiparticle backflow. The Galilean Ward identity

1+F1s3=m∗m1+\frac{F_1^s}{3}=\frac{m^*}{m}

restores the bare-mass coefficient. Thus the F0F_0 model is a controlled illustration of the density pole, but it is not a complete Galilean-invariant current response when m∗≠mm^*\ne m. On a lattice, the corresponding optical or stress sum rule replaces n/mn/m, and interband or incoherent weight may be essential. The Ward-response page develops this vertex logic in general.

Collisions turn zero sound into first sound

Section titled “Collisions turn zero sound into first sound”

The collisionless equation discarded Icoll[δn]I_{\mathrm{coll}}[\delta n] from the Landau–Boltzmann equation. In a clean neutral fluid, a usable collision operator must conserve particle number, momentum, and energy:

∑pσIcoll=0,∑pσp Icoll=0,∑pσϵpIcoll=0.\sum_{\mathbf p\sigma}I_{\mathrm{coll}}=0, \qquad \sum_{\mathbf p\sigma}\mathbf p\,I_{\mathrm{coll}}=0, \qquad \sum_{\mathbf p\sigma}\epsilon_{\mathbf p}I_{\mathrm{coll}}=0.

A naive replacement ω→ω+i/τcoll\omega\to\omega+i/\tau_{\mathrm{coll}} does not enforce these constraints and cannot be trusted to produce the hydrodynamic limit. A relaxation-time model must project out the local-equilibrium density, momentum, and energy modes. This issue is explicit in Abrikosov and Khalatnikov 1958, § 9, printed p. 79.

When ωτcoll≪1\omega\tau_{\mathrm{coll}}\ll1, collisions establish local equilibrium and the longitudinal density mode becomes first sound. Its speed is c12=(mnκS)−1c_1^2=(mn\kappa_S)^{-1}, where κS\kappa_S is the adiabatic compressibility. Using the T→0T\to0 thermodynamic coefficients of a clean, neutral, Galilean-invariant liquid—while retaining a small nonzero temperature or another momentum-conserving equilibration mechanism so that the hydrodynamic hierarchy still holds—gives

c12=(vF∗)23(1+F0s)(1+F1s3),c_1^2 =\frac{(v_F^*)^2}{3} (1+F_0^s)\left(1+\frac{F_1^s}{3}\right),

and, before dissipative corrections,

Rshyd(q,ω)=nq2/mc12q2−(ω+i0+)2.\mathcal R_s^{\mathrm{hyd}}(q,\omega) =\frac{nq^2/m} {c_1^2q^2-(\omega+i0^+)^2}.

At strictly T=0T=0 in a clean normal Fermi liquid, τcoll−1→0\tau_{\mathrm{coll}}^{-1}\to0 and this hydrodynamic window closes; c1c_1 is the T→0T\to0 hydrodynamic extrapolation.

First and zero sound solve different eigenproblems, so their speeds need not agree. For example, F0s=1F_0^s=1 and F1s=0F_1^s=0 give

c1vF∗=23=0.816497…,c0vF∗=s0=1.044382….\frac{c_1}{v_F^*}=\sqrt{\frac23}=0.816497\ldots, \qquad \frac{c_0}{v_F^*}=s_0=1.044382\ldots.

There is no contradiction in the first-sound speed lying inside the collisionless one-pair wedge: frequent collisions have already reorganized the dynamics into local equilibrium, and hydrodynamic attenuation is controlled by transport coefficients rather than collisionless Landau damping.

Viscosity and thermal conduction add sound attenuation. The intermediate regime ωτcoll∼1\omega\tau_{\mathrm{coll}}\sim1 is often strongly damped and requires the full collision integral. The relevant angular-harmonic relaxation rates scale as T2/EFT^2/E_F at low temperature, with probe-frequency corrections of order ω2/EF\omega^2/E_F and channel-dependent coefficients; they need not equal the on-shell one-particle lifetime. Cooling at fixed probe frequency can therefore drive first sound into zero sound, although pairing may intervene first. Landau 1957, printed p. 105 separates the hydrodynamic, crossover, and collisionless regimes, while Dupuis 2025, § 4.3.3.2, printed pp. 303–304 gives the first-sound speed in the present convention.

The derivation above assumes short-range or statically screened interactions, a smooth three-dimensional Fermi surface, a normal state, and a probe inside the Landau window. Several changes alter the conclusion:

  • In a three-dimensional charged fluid with unscreened Coulomb interactions, the longitudinal density oscillation is lifted to a gapped plasmon; see Silin 1959, printed p. 872, Eqs. (18)–(20) and the plasmon page. In two dimensions with three-dimensional Coulomb fields, the long-wavelength plasmon instead scales as q\sqrt q.
  • On a lattice, momentum relaxation, Umklapp, and anisotropy can make the hydrodynamic density mode diffusive and turn the angular equation into a Fermi-surface integral rather than a Legendre problem.
  • At finite temperature or disorder, a collisionless pole must be separated from both the continuum and its finite width. “Above the continuum” is not by itself an experimental visibility criterion.
  • Near pairing, density-wave, or Pomeranchuk instabilities, the normal-state Landau parameters and collision kernel acquire strong scale dependence; the assumed normal Fermi liquid may fail before the formal pole is reached.

Taking limits without naming the path. Compressibility, a uniform dynamic response, zero sound, and first sound approach the same origin in the (q,ω)(q,\omega) plane but are not interchangeable. State the order of limits and the size of ωτcoll\omega\tau_{\mathrm{coll}} before interpreting a formula.

Calling every denominator zero an undamped mode. A real pole on the physical sheet above the one-pair continuum is different from a resonance reached through a branch cut. Also compare its separation and residue with all sources of broadening.

Using vFv_F when the theory needs vF∗v_F^*. The kinetic equation transports quasiparticles, so its velocity is the interacting Fermi velocity. The bare mass reappears in a Galilean sum rule only after the matching current vertex and backflow are included.

Adding i/τcolli/\tau_{\mathrm{coll}} without conserving slow densities. This shortcut can erase the hydrodynamic zero modes of the collision operator. Projecting out the conserved local-equilibrium distortions is part of the approximation, not an optional refinement.

  1. Starting from (z−μ)ν=μ(F0νˉ−U)(z-\mu)\nu=\mu(F_0\bar\nu-U), derive the response R=−N(0)Φ/(1−F0Φ)\mathcal R=-N(0)\Phi/(1-F_0\Phi). Check both iterated limits at the origin.
Solution

Solve for the angular distortion and average:

νˉ=(F0νˉ−U)12∫−11μ dμz−μ=Φ(z)(F0νˉ−U).\bar\nu =(F_0\bar\nu-U) \frac12\int_{-1}^{1}\frac{\mu\,d\mu}{z-\mu} =\Phi(z)(F_0\bar\nu-U).

Therefore νˉ/U=−Φ/(1−F0Φ)\bar\nu/U=-\Phi/(1-F_0\Phi) and R=N(0)νˉ/U\mathcal R=N(0)\bar\nu/U. For ω→0\omega\to0 first, z→0z\to0 and Φ→−1\Phi\to-1, giving R=N(0)/(1+F0)\mathcal R=N(0)/(1+F_0). For q→0q\to0 at fixed nonzero ω\omega, z→∞z\to\infty and Φ→0\Phi\to0, giving R=0\mathcal R=0.

  1. Derive the weak- and strong-coupling limits of the zero-sound speed through the first nontrivial term.
Solution

For s=1+δs=1+\delta with δ≪1\delta\ll1,

Φ(1+δ)=12log⁡2δ−1+o(1).\Phi(1+\delta) =\frac12\log\frac{2}{\delta}-1+o(1).

Setting F0Φ=1F_0\Phi=1 gives δ=2e−2(1+1/F0)[1+o(1)]\delta=2e^{-2(1+1/F_0)}[1+o(1)]. For s≫1s\gg1, use

Φ(s)=13s2+15s4+O(s−6).\Phi(s)=\frac1{3s^2}+\frac1{5s^4}+O(s^{-6}).

Writing s2=F0/3+c+O(F0−1)s^2=F_0/3+c+O(F_0^{-1}) in F0Φ=1F_0\Phi=1 yields c=3/5c=3/5.

  1. For F0s=1F_0^s=1, verify the quoted speed and residue numerically, and explain why the residue is positive.
Solution

Solving Φ(s0)=1\Phi(s_0)=1 for s0>1s_0>1 gives s0=1.044382…s_0=1.044382\ldots. At this point Φ′(s0)<0\Phi'(s_0)<0, so

Z0N(0)qvF∗=−1Φ′(s0)=0.104217…>0.\frac{\mathcal Z_0}{N(0)qv_F^*} =-\frac1{\Phi'(s_0)} =0.104217\ldots>0.

The sign is also required by the positive-frequency spectral weight

Sspole=2πZ0δ(ω−ω0).S_s^{\mathrm{pole}} =2\pi\mathcal Z_0\delta(\omega-\omega_0).
  1. Show why an F0F_0-only model with m∗≠mm^*\ne m misses the Galilean density ff-sum coefficient, and identify the missing ingredient.
Solution

At large zz, the model gives

GsR≃N(0)(vF∗)2q23ω2=nq2m∗ω2.G_s^R\simeq \frac{N(0)(v_F^*)^2q^2}{3\omega^2} =\frac{nq^2}{m^*\omega^2}.

The Galilean commutator sum rule requires nq2/(mω2)nq^2/(m\omega^2). The ℓ=1\ell=1 interaction dresses the current through quasiparticle backflow; using 1+F1s/3=m∗/m1+F_1^s/3=m^*/m restores the required bare-mass coefficient and hence the ff-sum rule.

  1. Take F0s=1F_0^s=1, F1s=0F_1^s=0, and qvF∗=10−2EFqv_F^*=10^{-2}E_F. Classify the response when (a) τcollEF=10\tau_{\mathrm{coll}}E_F=10 and (b) τcollEF=103\tau_{\mathrm{coll}}E_F=10^3. In case (b), use τcoll−1\tau_{\mathrm{coll}}^{-1} as a broadening estimate and decide whether the zero-sound peak is cleanly separated from the continuum.
Solution

In case (a), the first-sound estimate gives

ωτcoll≃23 (10−2)(10)=0.0816≪1.\omega\tau_{\mathrm{coll}} \simeq \sqrt{\frac23}\,(10^{-2})(10) =0.0816\ll1.

The probe is hydrodynamic, so first sound is the appropriate mode. In case (b), the collisionless zero-sound root gives

ω0τcoll=s0(10−2)(103)=10.44≫1,\omega_0\tau_{\mathrm{coll}} =s_0(10^{-2})(10^3) =10.44\gg1,

so the kinetic regime is collisionless. Nevertheless, its separation from the leading continuum edge is only

Δω=(s0−1)qvF∗=4.44×10−4EF,\Delta\omega =(s_0-1)qv_F^* =4.44\times10^{-4}E_F,

whereas the assumed broadening scale is τcoll−1=10−3EF\tau_{\mathrm{coll}}^{-1}=10^{-3}E_F. The response can satisfy ω0τcoll≫1\omega_0\tau_{\mathrm{coll}}\gg1 without showing a cleanly resolved peak; regime classification and experimental visibility are separate tests.

  • A. A. Abrikosov and I. M. Khalatnikov, “Theory of the Fermi Fluid (The Properties of Liquid He³ at Low Temperatures),” Soviet Physics Uspekhi 1, 68–90 (1958), doi:10.1070/PU1958v001n01ABEH003086, Open PDF.
  • Nicolas Dupuis, Field Theory of Condensed Matter and Ultracold Gases, Vol. 1 (World Scientific, 2023), ch. 4; updated chapter PDF, January 16, 2025, doi:10.1142/q0409, author PDF.
  • L. D. Landau, “Oscillations in a Fermi Liquid,” Soviet Physics JETP 5, 101–108 (1957), official article, official PDF.
  • V. P. Silin, “The Oscillations of a Degenerate Electron Fluid,” Soviet Physics JETP 8, 870–875 (1959), official PDF.
  • Gordon Baym and Christopher Pethick, Landau Fermi-Liquid Theory: Concepts and Applications, Wiley, New York (1991), ch. 1, doi:10.1002/9783527617159.

Original QFT.org content:CC BY 4.0, unless an item supplies different terms. Third-party material retains its own terms.