上一讲我们讲了自由费米气体：无相互作用时，基态是填满的费米海，激发是一目了然的粒子、空穴。但真实体系有相互作用。这一讲介绍朗道（Landau）1956 年提出的费米液体理论——一种把强相互作用的费米子系统用费米气体的语言来理解的唯象理论。它描述了 $^3\text{He}$ 和几乎所有已知金属中的电子系统，是量子多体理论的基石之一。

相互作用的费米子与费米液体思想

典型的相互作用费米子哈密顿量是：

$H = H_0 + V_{ee} = \sum_{k\sigma}\varepsilon_{k\sigma}c^\dagger_{k\sigma}c_{k\sigma} + \sum_{k_1\sigma_1,k_2\sigma_2,q} V_q\, c^\dagger_{k_1-q\sigma_1}c^\dagger_{k_2+q\sigma_2}c_{k_2\sigma_2}c_{k_1\sigma_1}$

其中 $V_q$ 是相互作用势的傅里叶变换。对库仑相互作用，在三维空间中 $V_q = \dfrac{4\pi e^2}{q^2}$，对嵌入三维空间的二维系统 $V_q = \dfrac{2\pi e^2}{q}$。下面我们用”电子”泛指费米子（凝聚态中最常见）。

没有相互作用时，基态是填满的费米球。对于有相互作用、且电子数和总动量守恒的体系，基态 $|\psi_0\rangle$ 必须是电子气体能量本征态的叠加（在总电子数 $N$、总动量 $p=0$ 的子空间内）：

$|\psi_0\rangle = C_1\,|\text{费米海}\rangle + C_2\,|\text{费米海 + 粒子-空穴对}\rangle + C_3\,|\text{费米海 + 两个粒子-空穴对}\rangle + \cdots$

【配图】图 3.0.1：基态的叠加示意。费米海中挖去一个空位称为”空穴”（hole），图示为基态中的各种涨落成分。

为什么没有单粒子-空穴激发？动量守恒//为什么动量=0是基态？能量最低？

$|\psi_0\rangle$ 中有无穷多项，激发态更难描述。我们怎么对相互作用电子做预言？答案就是费米液体理论：它表明，我们可以用费米气体的同一种语言来理解强相互作用的费米子体系。

费米海中粒子的寿命

先做一些背景知识：费米面附近一个电子的寿命。用 $|G\rangle$ 表示费米气体基态。在费米海之上、费米面附近加一个动量 $k$、能量 $\varepsilon$ 的电子：$c^\dagger_k|G\rangle$。在 $t=0$ 时刻突然打开相互作用，之后初始态会演化成别的东西：

$|\psi(t)\rangle = c^\dagger_k|G\rangle + \underbrace{|\text{一阶散射项}\rangle}_{t \text{不太长时是主体}} + |\text{二阶项}\rangle + \cdots$

一阶散射项主要是激发单个电子-空穴对：初粒子散射到 $k_1$，同时另一个粒子从 $k_3$（费米面以下）被激发到 $k_2$（费米面以上）。后者是动量为 $k_2$、能量为 $\xi_{k_2}$ 的粒子激发和动量为 $-k_3$、能量为 $-\xi_{k_3}$ 的空穴激发。定义 $\xi_k \equiv \varepsilon_k - E_F$ 为相对费米能的单粒子能量。

用费米黄金规则估算衰变率：

$\begin{aligned} \gamma_\xi &= \frac{2\pi}{\hbar}\sum_{\substack{k_1+k_2-k_3=k\\ k_1,k_2>k_F;\ k_3<k_F}} |a(0;1,2,3)|^2\, \delta(\xi - (\xi_{k_1}+\xi_{k_2}-\xi_{k_3}))\\ &\sim \langle|a|^2\rangle\,\nu^3 \int_0^\infty d\xi_1 d\xi_2 d\xi_3\, \delta(\xi-(\xi_1+\xi_2-\xi_3))\\ &\sim \langle|a|^2\rangle\,\nu^3 \int_0^\xi d\xi_1 \int_0^{\xi_1} d\xi_2 \sim c_a\, \xi^2 \end{aligned}$

其中 $a$ 是散射振幅（相互作用矩阵元），$\nu$ 是单电子态密度。第二步把动量守恒替换为对散射振幅的适当平均，并把动量积分换成三个（能量都为正的）末态电子/空穴的能量积分，假定态密度为常数——只要初电子靠近费米面，这相当合理。第三步的积分范围来自能量守恒：末态电子与空穴的能量不能超过初粒子的能量。虽然近似看起来草率，但足以得到正确的标度 $\gamma_\xi \sim \xi^2$，这正是费米液体理论的动机。

结论：靠近费米面的电子衰变得非常慢。随着 $\xi \to 0$，其寿命发散！

这个估算太”快而糙”了：它在 3D 给出正确的 $\gamma_\xi \sim \xi^2$ 标度，但在 2D 会漏掉一个 $\ln(E_F/\xi)$ 因子。更精细的计算见 Giuliani & Vignale 的 §8.4 或 Coleman 的 §6.7.2。

对于长程相互作用的费米子（如金属中的库仑电子），散射振幅 $a$ 应当是被所有电子屏蔽后的相互作用矩阵元 $V_q/\epsilon(q) \sim 1/\nu$，几乎与原始相互作用强度无关（详见第六讲介电函数）。考虑到 $\gamma_\xi$ 具有能量量纲，且只应涉及初电子能量 $\xi$ 和费米面性质（$k_F, v_F$），由上面的估算可得 3D 下
$\gamma_\xi = c_\gamma \frac{v_F \xi^2}{k_F} \sim \frac{\xi^2}{E_F}$
奇妙的是，屏蔽相互作用算出的衰变率对相互作用强度几乎是非微扰的、与其无关。可以这样记忆：费米面越大，屏蔽使有效相互作用越弱，超过了散射相空间增大的效应。

准粒子

绝热对应

绝热定理（Born-Fock）：如果哈密顿量随时间足够缓慢地变化，且体系初始处于其本征态，那么演化后的状态仍是哈密顿量对应的本征态——只要本征值与谱的其余部分有能隙。例如无限深势阱缓慢变为谐振子阱时，单粒子本征态光滑地演化到新本征态，节点数（标签）不变。

【配图】图 3.2.1：绝热对应示意。(a) 单粒子势阱缓慢变形，本征态平滑演化；(b) 费米气体与费米液体谱之间的绝热对应：在时间 $T$ 内缓慢打开相互作用 $V_{ee}$，费米气体的基态演化成费米液体的基态（费米面变得模糊），单粒子激发态对应含一个额外费米子的激发态。

朗道受此启发：如果我们缓慢地打开相互作用 $V_{ee}$，非相互作用体系（费米气体）的能量本征态也会演化成相互作用体系（费米液体）的本征态，从而在两者之间建立一一对应。形式上，定义打开相互作用过程中的时间演化算符 $\hat{U}(T)$：

$|\psi_0\rangle = \hat{U}(T)|G\rangle$

对多一个电子的情形：

$|\psi_1\rangle = \hat{U}c^\dagger_k|G\rangle = \hat{U}c^\dagger_k\hat{U}^\dagger \hat{U}|G\rangle = \Psi^\dagger_k|\psi_0\rangle$

这里我们定义了新费米子的产生算符 $\Psi^\dagger_k$——这就是准粒子（quasi-particle）。可以验证幺正变换不改变对易/反对易关系。准粒子描述费米液体中”加一个粒子”的激发态（准空穴类似），是费米气体中普通粒子的绝热对应，具有相同的电荷、自旋 $\sigma$ 和动量 $k$（这些量在绝热打开过程中守恒）。

准粒子图像何时成立

【配图】图 3.2.2：准粒子成立的条件。纵轴为能量 $\xi$，画两条线：$\gamma_\xi \sim \xi^2/E_F$（粒子衰变率）和 $\Gamma = 1/T$（相互作用打开速率）。蓝色区域满足 $\gamma_\xi \ll \Gamma \ll \xi$。

假设相互作用在时间 $T$ 内逐渐打开，打开速率为 $\Gamma = 1/T$。打开过程给末态引入约 $\Gamma$ 量级的能量不确定度。要成为好的能量本征态，需要 $\Gamma \ll \xi_k$。但如 3.1 节所述，$c^\dagger_k|G\rangle$ 在此过程中以 $\gamma_\xi \sim \xi^2$ 的速率衰变；若演化时间无穷长，我们找不到具有确定动量 $k$ 的粒子，因此还需 $\Gamma \gg \gamma_\xi$。综上：

$\gamma_\xi \ll \Gamma \ll \xi_k$

这对靠近费米面的粒子没有问题，因为 $\gamma_\xi \sim \xi^2/E_F \ll \xi$。固定一个打开协议 $\Gamma$，只有能量在蓝色区域的粒子能变成良定义的准粒子；若对不同粒子用不同 $\Gamma$，则每个满足 $\gamma_\xi \ll \xi$ 的粒子都有良定义的准粒子对应（但这时它们没有良定义的费米反对易关系）。

朗道能量泛函与朗道参数

零温情形

既然费米液体的本征谱由基态 $|G\rangle$ 之上的准粒子激发标记，能量必是准粒子占据数的泛函。准粒子密度低时，把能量按偏离基态的占据数 $\delta n_{k\sigma} = n_{k\sigma} - n^{(0)}_{k\sigma}$ 展开（$n^{(0)}_{k\sigma} = \Theta(k_F - k)$）：

$E = E_G + \sum_{k\sigma}\xi_{k\sigma}\delta n_{k\sigma} + \frac{1}{2}\sum_{\substack{k_1k_2\\ \sigma_1\sigma_2}} f_{k_1\sigma_1,k_2\sigma_2}\delta n_{k_1\sigma_1}\delta n_{k_2\sigma_2} + \cdots$

其中 $k>k_F$ 时 $\delta n_{k\sigma}>0$，$k<k_F$ 时 $\delta n_{k\sigma}<0$。这就是朗道能量泛函。

第二项是单粒子贡献，系数就是单准粒子能量。由此定义准粒子的费米速度、有效质量和态密度：

$v_F^* \equiv \left.\partial_k\xi_k\right|_{k=k_F}, \qquad m^* \equiv \frac{k_F}{v_F^*}, \qquad \nu^* \equiv g_s\Omega_{D-1}k_F^{D-1}\frac{dk}{d\xi} = g_s\Omega_{D-1}k_F^{D-1}\frac{1}{v_F^*} = \frac{m^*}{m}\,\nu$

费米液体效应（相互作用效应）对单粒子激发的作用，仅仅是把费米速度 $v_F$ 重正化为 $v_F^*$。

第三项是两个准粒子之间的”相互作用”（注意它不同于式 (3.0.1) 中的相互作用项）。由自旋旋转对称性，相互作用强度一般可写成：

$f_{k\sigma,k'\sigma'} = f_{kk'}^s + f_{kk'}^a\,\sigma\sigma'$

第一项是自旋对称部分，第二项是自旋反对称部分。对费米面附近的准粒子，可忽略 $f$ 对 $|k|$ 的依赖、只保留角度依赖。由旋转不变性，$f$ 只依赖 $k$ 与 $k’$ 的夹角 $\theta$。定义无量纲相互作用强度 $F^{s,a}(\cos\theta) = \nu^* f^{s,a}(\cos\theta)$，用勒让德多项式展开：

$F^{s,a}(\cos\theta) = \sum_l (2l+1)\,F_l^{s,a}\,P_l(\cos\theta)$

其中 $P_l(x)$ 满足 $\int_{-1}^1 dx\,P_l(x)P_m(x) = \frac{2}{2l+1}\delta_{lm}$。系数 $F_l^{s,a}$ 就是朗道参数：

$F_l^{s,a} = \langle F^{s,a}(\Omega)P_l(\Omega)\rangle_\Omega \equiv \frac{1}{4\pi}\int \sin\theta\,d\theta d\varphi\, F^{s,a}(\cos\theta)P_l(\cos\theta) = \frac{1}{2}\int_{-1}^1 d\cos\theta\, F^{s,a}(\cos\theta)P_l(\cos\theta)$

$F_l^s$（自旋对称）参数化相互作用的非磁性部分，$F_l^a$（自旋反对称）定义磁性部分。利用 $\delta(\xi_k)$ 相当于向费米面的投影、$\sum_k g_s\delta(\xi_k) = \nu^*$，可写为：

$F_l^{s,a} = \frac{g_s}{\nu^*}\sum_{k'} F_{kk'}^{s,a}P_l(\cos\theta_{kk'})\delta(\xi_{k'}) = g_s\sum_{k'} f_{kk'}^{s,a}P_l(\cos\theta_{kk'})\delta(\xi_{k'})$

例子：对式 (3.0.1) 定义的体系，计算 $V_{ee}$ 对朗道参数的领头阶贡献（忽略费米速度重正化）。两体相互作用系数为

$f_{k_1\sigma_1,k_2\sigma_2} = \langle\Phi|\hat{V}_{ee}|\Phi\rangle = \sum_q V_q\langle\Phi|c^\dagger_{k_1+q\sigma_1}c^\dagger_{k_2+q\sigma_2}c_{k_2\sigma_2}c_{k_1\sigma_1}|\Phi\rangle = V(0) - V(k_1-k_2)\delta_{\sigma_1\sigma_2}$

其中 $|\Phi\rangle = c^\dagger_{k_1\sigma_1}c^\dagger_{k_2\sigma_2}|G\rangle$。两项分别称为”直接”（Hartree）和”交换”（Fock）项。对自旋 $1/2$，利用 $\delta = \frac12\begin{pmatrix}1&1\\1&1\end{pmatrix} + \frac12\begin{pmatrix}1&-1\-1&1\end{pmatrix}$（自旋对称/反对称分解），得

$f_{k_1k_2}^s = V(0) - \frac12 V(k_1-k_2), \qquad f_{k_1k_2}^a = -\frac12 V(k_1-k_2)$

于是朗道参数（费米面上 $\cos\theta = \dfrac{k_1\cdot k_2}{k_F^2}$，$|k_1-k_2| = 2k_F\sin\frac{\theta}{2}$）：

$\begin{cases} F_0^s = \nu\int_{-1}^{1}d\cos\theta\left[V(0) - \frac12 V(2k_F\sin\frac{\theta}{2})\right]\\[0.4em] F_1^s = \nu\int_{-1}^{1}d\cos\theta\left[V(0) - \frac12 V(2k_F\sin\frac{\theta}{2})\right]\cos\theta\\[0.4em] F_0^a = \nu\int_{-1}^{1}d\cos\theta\left[-\frac12 V(2k_F\sin\frac{\theta}{2})\right]\\[0.4em] \cdots \end{cases}$

对屏蔽库仑势 $V(q) = \dfrac{4\pi e^2}{q^2+q_{TF}^2}$（$q_{TF}^2 \sim 4\pi e^2 2\nu \sim \alpha_M k_F^2$）：

$F_0^s = 1 - \frac{\alpha_M}{8}\ln\left(1+\frac{4}{\alpha_M}\right), \qquad F_0^a = -\frac{\alpha_M}{8}\ln\left(1+\frac{4}{\alpha_M}\right)$

$F_0^a$ 的负号意味着什么？请思考（提示：它使自旋方向上的响应被增强，有利于磁性）。

非零温度

有了式 (3.3.1) 的本征谱，非零温度下能量 $E$ 的状态概率是玻尔兹曼因子 $\frac1Z e^{-E}$。但这是个不平凡的统计问题：不同 $k$ 之间存在关联，密度矩阵不再是各子系统 $\hat{\rho}_{k\sigma}$ 的直积。做平均场近似（假设密度矩阵仍是乘积态），自由能可写成准粒子占据数期望值的简单泛函：

$\begin{aligned} F &= \langle E\rangle_{\hat{\rho}} - TS = E_G + \sum_{k\sigma}\xi_{k\sigma}\delta n_{k\sigma} + \frac12\sum_{\substack{k_1k_2\\\sigma_1\sigma_2}} f_{k_1\sigma_1,k_2\sigma_2}\delta n_{k_1\sigma_1}\delta n_{k_2\sigma_2} + \cdots\\ &\qquad + T\sum_{k\sigma}\left(n_{k\sigma}\ln n_{k\sigma} + (1-n_{k\sigma})\ln(1-n_{k\sigma})\right) \end{aligned}$

对 $\delta n_{k\sigma}$ 最小化（$\partial_{\delta n_{k\sigma}}F = 0$），得平衡条件：

$\xi_{k\sigma} + k_B T\ln\frac{n_{k\sigma}}{1-n_{k\sigma}} = 0 \quad\Rightarrow\quad n_{k\sigma} = \frac{1}{e^{\xi_{k\sigma}/k_BT}+1} = f_F(\xi_{k\sigma})$

其中

$\xi_{k\sigma} = \frac{\partial\langle E\rangle_{\hat{\rho}}}{\partial\delta n_{k\sigma}} = \xi^{(0)}_{k\sigma} + \sum_{k'\sigma'} f_{k\sigma,k'\sigma'}\delta n_{k'\sigma'} + \cdots$

是被其他激发准粒子通过相互作用项重正化后的有效准粒子能量。注意这个 $\xi_{k\sigma}$ 与第二章费米气体中的记号相同但含义不同。结论：准粒子的热分布仍是费米-狄拉克分布，但其中的有效能量包含了相互作用的反馈效应。

相互作用对热力学量的修正

比热

与费米气体相比，相互作用通过态密度的同一个因子改变比热：

$C_V = \sum_{k\sigma}\xi_{k\sigma}\partial_T n_{k\sigma} \xrightarrow{T\ll\varepsilon_F} \nu^*\int d\xi\,\xi\,\partial_T f(\xi,T) = \nu^*\int d\xi\,(-\partial_\xi f) \sim \nu^* T$

电荷磁化率

电荷磁化率指化学势增加 $\delta\mu$ 引起的粒子数增加 $\delta N$。更方便的做法是已知 $\delta N$ 反求 $\delta\mu$。有效准粒子能量由于反馈势被修正：

$\delta\mu = \xi = \delta\mu^{(0)} + \sum_{k'}f_{kk'}^s g_s\delta n_{k'} \xrightarrow{T\to 0} (1+F_0^s)\delta\mu^{(0)} = (1+F_0^s)\frac{\delta N}{\nu^*}$

其中 $\delta\mu^{(0)}$ 由 $\delta N = \nu^*\delta\mu^{(0)}$ 定义，$F_0^s$ 项反映费米球整体膨胀对单准粒子能量的修正（反馈势）。因此电荷磁化率为：

$\chi_c = \frac{\delta N}{\delta\mu} = \frac{\nu^*}{1+F_0^s}$

【配图】图 3.4.1：(a) 电荷磁化率与 $F_0^s$：粒子数增加 $\delta N$（红色区域）使费米能级从 0 升到 $\xi$；(b) 自旋磁化率与 $F_0^a$：磁场使自旋上、下的能级相对移动。

$F_0^s>0$（排斥相互作用）使液体变”硬”，更难压缩。例如 $^3\text{He}$ 在某个区域 $F_0^s = 10.8$。

自旋磁化率与 Stoner 失稳

有磁场 $B$ 时，自旋上与下的能量移动等大反号。对自旋向下电子：

$\delta\xi_\downarrow = \mu_b B + \sum_{k'}f_{k\downarrow,k'\sigma'}\delta n_{k'\sigma'} = \mu_b B - \delta\xi_\downarrow F_0^a$

因此自旋磁化率：

$\chi_s = \frac{\delta M}{B} = \frac{\mu_b\nu^*\delta\xi_\downarrow}{B} = \frac{\nu^*\mu_b^2}{1+F_0^a}$

$F_0^a$ 反映自旋极化对准粒子的反馈场。当 $F_0^a \to -1$，$\chi \to \infty$，表明铁磁失稳，即 Stoner 失稳。

一般多极失稳：Pomeranchuk 失稳

一般地，考虑对具有某种多极对称性的准粒子施加外势，使裸准粒子能量改变 $\delta\xi^{(0)}_{p\sigma} = v_l Y_{lm}(\hat{p})$。准粒子分布响应并达到自洽稳定态后，总势 $\delta\xi_{p\sigma} = t_l Y_{lm}(\hat{p})$ 与裸势通过朗道参数联系：反馈势为 $-F_l t_l$，于是

$t_l = v_l - F_l t_l \quad\Leftrightarrow\quad t_l = \frac{v_l}{1+F_l}$

$F_l>0$ 时总势被抑制（排斥相互作用），$F_l<0$ 时被增强（吸引）。临界情形 $F_l = -1$ 时费米球对某种形变失稳，称为 Pomeranchuk 失稳。

有效质量与 F₁ᵈ：伽利略不变性

费米液体参数并非全部独立。例如在有伽利略不变性的体系中，$F_1^s$ 与费米面上的有效质量有密切关系（朗道最先讨论）。伽利略不变性指式 (3.0.1) 中的裸粒子具有二次色散 $\varepsilon_k = k^2/(2m)$，在移动参考系（$t’ = t,\ x’ = x + vt$）中服从相同的物理规律。

考虑处于基态的费米液体，在以速度 $v = p/m$ 运动的参考系中观察，它表现为费米面发生形变的激发态。用两种方式计算其激发能：

质心动能（唯一可能的宏观答案）：

$\Delta E_G = \frac12 nmv^2 = n\frac{p^2}{2m}$

朗道能量泛函：把准粒子激发数展开到 $p$ 的二阶：

$\Delta E_L = \frac{p^2}{2m^*}n\,(1+F_1^s)$

（推导中用到 $\nu^v_F^{2}/D = n/m^*$。）

两者应相等，于是得到有效质量重正化与 $F_1^s$ 的关系：

$\frac{m^*}{m} = \frac{\nu^*}{\nu} = \frac{v_F}{v_F^*} = 1 + F_1^s$

【配图】图 3.5.1：(a) 伽利略不变性示意：实验室系（左）与移动参考系（右）中的费米液体，上下分别是实空间与动量空间（准粒子占据）图；(b) 有空间涨落的费米液体实空间示意：每个空间位置的小圆代表该处的费米面形变。

理解：$F_1^s$ 是费米面 $l=1$ 形变对准粒子的反馈势。当 $F_1^s>0$ 使 $m^>m$ 时，朗道泛函的单准粒子能量项给出 $n\frac{p^2}{2m^}$，低估了实际能量 $n\frac{p^2}{2m}$；激发准粒子之间的正 $F_1^s$ 相互作用能恰好补上差额。

另一个理解角度是电荷流。伽利略不变性要求总电流 $j = np/m$；但单准粒子能量 $\xi^{(0)}_k$ 给出群速度 $k/m$ 和电流 $j = np/m^*$，低估了电流。解释是相互作用项 $F_1^s p\cos\theta_k$ 的反馈势给整个液体额外速度，恰好补上差额。

问题 3.5.1：单个准电子态 $\Psi^\dagger_k|G\rangle$ 的电荷流是多少？是 $k/m$ 还是 $k/m^*$？

解答：是 $e\mathbf k/m$（裸质量 $m$），不是 $e\mathbf k/m^$。$1/m^$ 只出现在”准粒子相对随流背景的速度”里，而总电流永远由裸质量决定。下面三层论证。

(1) 算符层面：伽利略不变性 $\Rightarrow$ 电流算符是纯抗磁的。 记 $\mathbf K$ 为伽利略推动算符。对 $\hat H_0=\sum_i\hat p_i^2/2m$ 与只依赖 $\mathbf r_i-\mathbf r_j$ 的 $\hat V_{ee}$（平移与推动均不变），
$U(\mathbf u)HU^\dagger(\mathbf u)=\hat H-\mathbf u\cdot\hat{\mathbf P}+\tfrac12Nm u^2$
（$\mathbf K$ 中 $\propto\sum_i\mathbf r_i$ 的部分对固定 $N$ 只是常数，不影响电流）。因 $H(\mathbf u)$ 与 $H$ 幺正等价、谱相同，而它只比 $H$ 多出 $-\mathbf u\cdot\hat{\mathbf P}$，故由定理其电流与 $H$ 相同，即
$\mathbf J=e\frac{\partial\mathbf K}{\partial\mathbf u}=\frac{e}{m}\hat{\mathbf P}_{tot}\quad\Longrightarrow\quad \boxed{\ \mathbf J\big|_{\Psi^\dagger_k|G\rangle}=\frac{e\mathbf k}{m}\ }$
这里 $m$ 只能是裸质量：相互作用的平移不变性使 $\hat{\mathbf P}$ 在电流算符中的”顺磁”贡献与逆场（diamagnetic）部分恰好抵消（Nielsen–Luttinger 定理）。而绝热打开相互作用守恒动量（§3.2），$|\psi_0\rangle$ 有 $\mathbf P=0$，故 $\Psi^\dagger_k|\psi_0\rangle$ 的总动量精确等于 $\mathbf k$。这一步严格，不含任何费米液体参数——它与本节用能量 $\frac12Nmv^2$ 定出 $m^*/m=1+F_1^s$ 用的是同一个伽利略不变性输入，只是电荷流这个可观测量取的是 $m$ 而已。

(2) 为什么”准粒子速度 $\Rightarrow$ 电流”的naïve 做法给出错误的 $k/m^*$。 试算 $\mathbf j=e\sum_k\mathbf v_k^n_k$：取整体以 $\mathbf u$ 平移的液体，Landau 第二定律给运动液体中准粒子能量 $\xi_k-\mathbf k\cdot\mathbf u$，故实验室系速度为 $\mathbf v_k^-\mathbf u$；又平衡占据 $n_k=f_F(\xi_k-\mathbf k\cdot\mathbf u)$ 使实验室系下无净流动。代入并做变量平移 $\mathbf k\to\mathbf k-m^*\mathbf u$（只保留 $O(u)$）：
$\mathbf j=e\sum_kn_k\left(\frac{\mathbf k}{m^*}-\mathbf u\right)=e\left(\frac{\mathbf P}{m^*}+n\mathbf u-n\mathbf u\right)=\frac{e\mathbf P}{m^*}$
正是错误答案。矛盾的出现说明 $\mathbf j\neq e\sum_k\mathbf v_kn_k$：单靠”准粒子速度乘占据数”凑不出电流算符，缺的正是 Landau 碰撞积分中的反流项（backflow）——被推动的准粒子带动周围液体反向流动，把差额补齐：
$\mathbf j=\frac{e\mathbf P}{m^*}+e\mathbf P_{\rm bf}=\frac{e\mathbf P}{m},\qquad \underbrace{\frac{e\mathbf k}{m^*}}_{\text{准粒子自身}}+\underbrace{\left(\frac{e\mathbf k}{m}-\frac{e\mathbf k}{m^*}\right)}_{=-\frac{eF_1^s\mathbf k}{m^*}\,\text{（反流，与 }\mathbf k\text{ 反向）}}$
这就是本节”反馈势给整个液体额外速度”的定量版本。反流项同时保证准粒子云输运总动量的速度为 $\mathbf v^*/(1+F_1^s)=\mathbf k/m$，与上面一致。

(3) $k/m^*$ 到底是什么。 它是准粒子相对随流液体的速度，出现在 Landau 动力学方程中：$\dot r=\partial_k\xi_k-\mathbf u$（本节式中的 $\dot r=\partial_k\xi^{(0)}_k=\mathbf v_k^$ 只是 $\mathbf u=0$ 的特例）。于是两种”质量”分工明确：$m^$ 刻画能量—动量关系（比热 $C_V\sim\nu^T$、态密度、de Haas–van Alphen 测 $\nu^/m$），$m$ 刻画电荷—动量关系（电流、Drude 权重 $n/m$、超流刚度）。二者被伽利略不变性以 $m^/m=1+F_1^s$ 绑在一起——**热力学量看到 $m^$，电荷输运看到 $m$*（$^3$He 中 $m^/m\approx2.6$、$F_1^s\approx1.6$，差别显著）。

例外与易错点：① 若伽利略不变性被破坏（晶格、能带非抛物、自旋轨道耦合、非平庸拓扑能带），则 $\mathbf J\neq e\mathbf P/m$，电流由能带速度 $\mathbf v_{n\mathbf k}=\hbar^{-1}\nabla_k\varepsilon_{n\mathbf k}$ 决定，抛物能带时关系改为 $m^/m=1+\dfrac{F_1^s}{1+F_1^s/D}$；② 若问的是”相对共动背景的电流”，那是 $e(\mathbf v_k^-\mathbf u)$ 那部分，与 $e\mathbf k/m$ 无关——但总电流恒为 $e\mathbf P/m$；③ $F_1^s\to0$（弱相互作用）时 $k/m\to k/m^*$，两者才无从分辨。

中性费米液体中的集体模式

前面我们处理的是空间均匀的费米液体。现在允许它有时空依赖，讨论大尺度动力学。若空间变化相对粒子间距足够缓慢（特征波矢 $q\ll k_F$），每个空间位置的局域液体都可以用前面均匀费米液体的准粒子占据函数描述（见上图 3.5.1(b)）。

在大尺度（$q\ll k_F$）、慢（$\omega\ll\varepsilon_F$）动力学中，每个准粒子可看作位置和动量都相对确定的小波包，在平均场中按经典运动方程运动。体系动力学由玻尔兹曼动力学方程描述：

$\frac{d}{dt}n_{k\sigma} = I[n_{k\sigma}]$

其中 $n_{k\sigma}(r,t)$ 是空间-时间依赖的分布函数，$I[n]$ 是碰撞积分。全导数 $d_t = \partial_t + \dot{r}\cdot\partial_r + \dot{k}\cdot\partial_k$，其中群速度 $\dot{r} = \partial_k\xi_k(r)$，受力 $\dot{k} = -\partial_r\varepsilon_k(r)$。

无碰撞区：零声

关注碰撞项可忽略的情形（频率 $\omega$ 远大于典型准粒子散射率 $\gamma_{ee}$）。考虑均匀平衡费米液体上的小涨落：

$n_k(r) = f_F(\xi_k) + \delta n_k(r,t)$

（为简单起见只考虑自旋对称通道 $\delta n_{k,\sigma} = \delta n_k$。）涨落描述在自身反馈场中集体运动的准粒子。代入并保留 $\delta n_k$ 的一阶项，得到线性化输运方程：

$(\partial_t + \dot{r}\cdot\partial_r)\delta n_k + (\partial_k f_F)\cdot\dot{k} = 0$

速度与力为：

$\dot{r} = \partial_k\xi^{(0)}_k = v_k^*, \qquad \dot{k} = -\partial_r\varepsilon_{k,\sigma} = -\partial_r\sum_{k',\sigma'}f_{k\sigma,k'\sigma'}\delta n_{k'\sigma'}(r)$

力来自小涨落的反馈势，对 $\delta n_k(r)$ 线性。由平移对称性，取平面波解 $\delta n_k(r) = \alpha_k e^{i(q\cdot r-\omega t)}$。由于准粒子占据涨落只在费米面附近，可写 $\alpha_k = \eta_{\hat{k}}(-\partial_{\xi_k}f_F)$。定义无量纲波速 $s = \omega/(v_F^*q)$ 和 $\cos\theta_k = \hat{k}\cdot\hat{q}$，得本征模方程：

$\eta_{\hat{k}} = \frac{\cos\theta_k}{s-\cos\theta_k}\int d\hat{k}'\,F_{kk'}\,\eta_{\hat{k}'}$

若朗道参数中只有 $F_0^s$ 非零（$F_{kk’} = F_0^s$），则积分 $\int d\hat{k}’F_{kk’}\eta_{k’}$ 必为常数，本征模的角度分布为 $\eta_{\hat{k}} = A\frac{\cos\theta_k}{s-\cos\theta_k}$。代入得：

$1 = F_0^s\int_{-1}^{1}\frac{d\cos\theta}{2}\frac{\cos\theta}{s-\cos\theta} \quad\Rightarrow\quad \frac{1}{F_0^s} = -1 + \frac{s}{2}\ln\frac{s+1}{s-1} \quad\Rightarrow\quad \begin{cases} \frac{1}{3s^2}, & s\gg 1\\[0.3em] \frac12\ln\frac{2}{s-1}, & s\to 1^+ \end{cases}$

这个波叫零声（zero sound），是中性费米液体在无碰撞区的集体模式：准粒子之间的相互作用使费米面整体集体振动。当 $F_0^s\gg 1$ 使 $s\gg 1$ 时，零声速度

$u = sv_F^* = \sqrt{\frac{F_0^s}{3}}v_F^*$

相互作用很弱时 $F_0^s\to 0^+$，零声速度趋近费米速度 $u\to v_F^$。当 $s<1$，集体模被 *朗道阻尼（Landau damping）强烈衰减（集体模共振激发粒子-空穴对，把能量耗散掉，寿命极短，不再是良定义的集体模）。

补充：零声推导——先认变量，再走步骤，最后说人话

一、先把变量认清楚

真正需要解的未知量只有一类：涨落的占据 $\delta n_k$；其余符号都是”谁在哪、跑多快”的标签。

符号含义备注

$k,\ \mathbf k$ 准粒子动量，$\lvert\mathbf k\rvert\simeq k_F$ 只关心费米面附近

$\hat k$，$\cos\theta_k$ 动量方向；与波方向 $\hat q$ 的夹角余弦全篇唯一的几何变量

$\xi_k$ 相对费米能的单准粒子能量 $\xi_{k_F}=0$

$f_F(\xi_k)$ 平衡占据（费米-狄拉克分布） $T\to0$ 时是阶跃函数

$n_k(r,t)$ 时空依赖的准粒子占据朗道理论的动力学变量

$\delta n_k$ 小涨落，$n_k=f_F(\xi_k)+\delta n_k$ 要解的未知量；取自旋对称通道 $\delta n_{k\sigma}=\delta n_k$

$\varepsilon_{k\sigma}(r,t)=\xi^{(0)}_{k\sigma}+\sum_{k’\sigma’}f_{k\sigma,k’\sigma’}\delta n_{k’\sigma’}$ 涨落自己产生的反馈势 扮演”外场”，自洽机制的关键

$f^s_{kk’},\ f^a_{kk’}$ 朗道泛函中的准粒子间相互作用与 $V_{ee}$ 无关，只是被重参数化

$F_l^{s,a}$，$F_{kk’}=\nu^*f^s_{kk’}$ 无量纲朗道参数 $l$ 是勒让德展开的阶

$\nu^*$ 费米面态密度 $\propto1/v_F^*$，第 4 步靠它换算

$\mathbf v_k^=\partial_k\xi^{(0)}_k=\frac{k}{m^}$ 准粒子群速度 $v_F^=k_F/m^$ 是费米面上的值

$q,\ \omega$ 涨落的波矢、频率平面波解

$\alpha_k\ \to\ \eta_{\hat k}$ 平面波振幅；费米面上的角分布二者关系见第 3 步

$-\partial_\xi f_F=\frac{1}{4T}\mathrm{sech}^2(\xi/2T)$ 态密度权重线性理论里涨落的”有效大小”

$s=\frac{\omega}{qv_F^*}$ 无量纲波速，即 $u/v_F^*$ 主角：$s>1$ 还是 $<1$ 决定生死

$I[n_{k\sigma}]$，$\gamma_{ee}\sim T^2/\varepsilon_F$ 碰撞积分、散射率无碰撞区（$\omega\gg\gamma_{ee}$）置 $I=0$

二、推导过程

第 0 步 · 适用条件。零声只存在于中性液体：带电体系的反馈势被等离子体屏蔽（$F_l^{s,a}\to0$），集体模退化为声子（plasmon）。两个尺度分离支撑”小波包”图像：$q\ll k_F$ 使位置、动量都确定，经典运动方程可用；$\omega\ll\varepsilon_F$ 使涨落只在费米面附近。

第 1 步 · 反馈势就是外场。有效单准粒子能量 $\xi_{k\sigma}(r,t)=\xi^{(0)}_{k\sigma}+\sum_{k’\sigma’}f_{k\sigma,k’\sigma’}\delta n_{k’\sigma’}$ 依赖 $r$，于是它扮演外势，波包按经典方程运动 $\dot r=\partial_k\xi_{k\sigma}$、$\dot k=-\partial_r\xi_{k\sigma}$；输运方程的左端正是沿这条轨迹的占据变化率。

第 2 步 · 线性化，得到自洽方程。设 $n_k=f_F(\xi_k)+\delta n_k$，$\lvert\delta n\rvert\ll1$。平衡部分不贡献（$\frac{d}{dt}f_F=0$），又因无外场时 $\partial_r\xi^{(0)}_k=0$，力完全来自涨落：
$\dot k=-\sum_{k'\sigma'}f_{k\sigma,k'\sigma'}\nabla\delta n_{k'\sigma'}$
代回后得到整个问题的骨架——左边自由流，右边被自己推：
$\underbrace{(\partial_t+\mathbf v_k^*\cdot\nabla)\delta n_k}_{\text{涨落随准粒子自由流过}}\;=\;\underbrace{(\nabla_kf_F)\cdot\sum_{k'\sigma'}f_{k\sigma,k'\sigma'}\nabla\delta n_{k'\sigma'}}_{\text{涨落产生的场反过来推它}}$
只取自旋对称通道时 $f^a$ 项消失（它乘 $\sum_{\sigma’}\delta n_{k’\sigma’}=0$），场中只剩 $f^s_{kk’}$。

第 3 步 · 平面波解与投影。平移对称性允许取 $\delta n_k(r,t)=\alpha_k e^{i(\mathbf q\cdot\mathbf r-\omega t)}$，再写 $\alpha_k=\eta_{\hat k}(-\partial_{\xi_k}f_F)$，理由有三条：① 涨落只能非零于 $\lvert\xi_k\rvert\lesssim\omega$ 的薄壳层内（低频激发离不开费米面）；② 方程中 $\delta n$ 只线性出现，其作用权重恰是态密度，壳层内变化平缓可提到壳外；③ 壳层极薄，剩下的自由度只有方向 $\hat k$。于是未知量从”每个 $k$ 各一个数”压缩为”费米面上的角分布 $\eta_{\hat k}$”。

第 4 步 · 代入（关键代数）。代入 $\partial_t\delta n_k=-i\omega\delta n_k$、$\nabla\delta n_k=i\mathbf q\,\delta n_k$ 并移项：
$-(\omega-\mathbf q\cdot\mathbf v_k^*)\delta n_k=-\mathbf q\cdot(\nabla_kf_F)\sum_{k'}f_{kk'}\delta n_{k'}$
再用 $\mathbf q\cdot\nabla_kf_F=-q\cos\theta_k(-\partial_{\xi_k}f_F)$ 与 $\mathbf q\cdot\mathbf v_k^=qv_F^\cos\theta_k$，代入投影式后 $(-\partial_\xi f_F)$ 两侧约掉：
$(\omega-qv_F^*\cos\theta_k)\eta_{\hat k}=q\cos\theta_k\int d\hat k'\,f_{kk'}\eta_{\hat k'}$
两边除以 $qv_F^$，$1/v_F^$ 由投影定义 $\sum_k g_s\delta(\xi_k)=\nu^$ 吸收成 $F_{kk’}=\nu^f_{kk’}$，即本节的本征模方程：
$\eta_{\hat k}=\frac{\cos\theta_k}{s-\cos\theta_k}\int d\hat k'\,F_{kk'}\eta_{\hat k'}$
这里藏着两处物理信息：反馈场只沿 $\mathbf q$ 方向（$\cos\theta_k$，即 $l=1$ 分量被激发），而角分布被 $1/(s-\cos\theta_k)$ 放大，其共振点正是 $\cos\theta_k=s$。

第 5 步 · 只留 $F_0^s$（各向同性）。$F_{kk’}=F_0^s$ 与方向无关时，右边积分与 $\hat k$ 无关（记作常数 $A$），角分布被唯一确定：$\eta_{\hat k}=A\frac{\cos\theta_k}{s-\cos\theta_k}$。可见零声的密度形变是”偶极式”的：顺波方向与背波方向不对称。

第 6 步 · 求根。把上一步的角分布代回积分，令 $\mu=\cos\theta$ 并作代数分解 $\frac{\mu}{s-\mu}=\frac{s}{s-\mu}-1$：
$\int_{-1}^{1}\frac{d\mu}{s-\mu}=\Big[-\ln\lvert s-\mu\rvert\Big]_{-1}^{1}=\ln\frac{s+1}{s-1},\qquad \int_{-1}^{1}d\mu\,\frac{\mu}{s-\mu}=s\ln\frac{s+1}{s-1}-2$
于是得零声色散关系 $\frac{1}{F_0^s}=-1+\frac{s}{2}\ln\frac{s+1}{s-1}$。

第 7 步 · 两个极限。$s\gg1$ 时展开 $\ln\frac{s+1}{s-1}=\frac2s+\frac{2}{3s^3}+\cdots$，得 $\frac{1}{F_0^s}\simeq\frac{1}{3s^2}$，即 $u_0=sv_F^=\sqrt{\frac{F_0^s}{3}}v_F^$；$s\to1^+$ 时 $\ln\frac{s+1}{s-1}\approx\ln\frac{2}{s-1}\to+\infty$，故 $F_0^s\to0^+$ 对应 $s\to1^+$、$u_0\to v_F^$。*$s<1$ 没有实根：$s-1<0$ 时对数取复值，方程右边成为复数。

第 8 步 · $s<1$ 的阻尼，以及”声速 vs 费米速度”。极点 $\cos\theta=s$ 落进积分区间 $(-1,1)$，对应费米面上有一组准粒子满足 $\mathbf q\cdot\mathbf v_k^=\omega$：它们群速度与波同相，被波卷入集体运动并被波带走能量，即集体模共振激发粒子-空穴对，这些能量在 $\gamma_{ee}\sim T^2/\varepsilon_F$ 的碰撞中被耗散，寿命 $\sim1/\gamma_{ee}$，不再是良定义的集体模。严格地说无碰撞方程在 $s<1$ 无实解，定量阻尼率必须把碰撞项 $I$ 加回来、用 Sokhotski–Plemelj 公式处理极点（解为复频率 $\omega=\omega_1-i\omega_2$，$\omega_2$ 即阻尼率）；由于极点紧邻积分区间，$s$ 略小于 1 时阻尼极强。另外，由笔记式 3.5.14 的 $v_F/v_F^=m^/m=1+F_1^s$（注意方向：$v_F^/v_F=1/(1+F_1^s)$）可把”声速超过费米速度”写成定量判据：
$\frac{u_0}{v_F}=\sqrt{\frac{F_0^s}{3}}\cdot\frac{v_F^*}{v_F}=\sqrt{\frac{F_0^s}{3}}\ \frac{1}{1+F_1^s}>1\quad\iff\quad F_0^s>3(1+F_1^s)^2$
液氦-3（$F_0^s\approx10.8$、$m^/m\approx2.6$）给出 $s\simeq\sqrt{F_0^s/3}\approx1.9$（展开勉强成立）、$u_0\approx1.9v_F^\approx1.2\times10^2$ m/s，实验值约 190 m/s，理想气体 $v_F=166$ m/s——$u_0$ 与 $v_F$ 同一量级但略小，而 $u_0>v_F^*$ 是显然的。

三、用大白话重讲一遍

想象一锅静止的液氦。有人往里扔了一个波：某处密度多了一点。

多出来的这些准粒子自己会造出一个”势”，像水面波峰的落差那样，把旁边的准粒子推一把。波因此不是自由跑的——波和它自己造成的推力必须对得上账（这就是第 2 步那个等式，也是”自洽”两个字的全部含义）。

推力只沿波前进的方向有分量，所以形变不是均匀鼓包，而是”前面挤、后面松”（第 4 步的 $\cos\theta_k$），像一个波前薄、一个波后厚的单向波。

致命之处在于：总有一部分准粒子的速度恰好等于波速，它们会被波”搭上便车”、被卷进集体运动、能量被波带走——这就是阻尼的来源。波跑得越慢，能搭上便车的准粒子越多，衰减越狠。

所以只有当波跑得比最快的准粒子还快（$s>1$，即 $u>v_F$）时，才没有任何准粒子追得上它，模式才能无阻尼地传下去——这就是零声。相互作用越强（$F_0^s$ 越大），它跑得越快：$u_0=\sqrt{F_0^s/3}\,v_F^$；相互作用趋于零时它退化到 $v_F$ 而*不趋零，这也正是朗道把它与普通声波（一声）在命名上区分开的关键。

补充：$F_0^s\to0$ 的无相互作用气体为什么仍然有声速——简并压

本节有个容易卡住的反常：零声方程给出 $F_0^s\to0^+\Rightarrow s\to1^+$、$u_0\to v_F^$，不是 $u_0\to0$。既然相互作用都没了，声靠什么传？而且理想费米气体（$F_l\equiv0$）确实有声（一声，$u_1=v_F/\sqrt3$）。两件事的答案是同一件：*简并压。

(1) 没有相互作用 $\neq$ 没有压力。 理想费米气体的压强完全由泡利原理给出——每态最多填一个粒子，密度一涨就得把外层电子推到更高能量：
$P=\frac{2}{D+2}n\varepsilon_F\xrightarrow{D=3}\frac25n\varepsilon_F,\qquad \varepsilon_F=\frac{\hbar^2k_F^2}{2m}$
这是 $T=0$ 的”量子压”（与中子星的简并压同源），与 $V_{ee}$ 无关。它给出的体弹模量
$\kappa=-V\frac{\partial P}{\partial V}=n\frac{\partial P}{\partial n}=\frac{2}{D}n\varepsilon_F$
于是自由费米气体的一声
$u_1^2=\frac{\kappa}{n}=\frac{2\varepsilon_F}{D}=\frac{v_F^2}{D}\quad\xrightarrow{D=3}\quad u_1=\frac{v_F}{\sqrt3}$
与本节弱相互作用极限 $u_1=v_F^/\sqrt D$ 完全吻合（$F_l=0\Rightarrow v_F^=v_F$）。

(2) 在朗道框架里这是显然的。 朗道能量泛函（式 3.3.1）分成两半：
$E=E_G+\underbrace{\sum_{k\sigma}\xi^{(0)}_{k\sigma}\delta n_{k\sigma}}_{\text{自由气体的动能项，永远存在}}+\frac12\underbrace{\sum_{k\sigma,k'\sigma'}f_{k\sigma,k'\sigma'}\delta n_{k\sigma}\delta n_{k'\sigma'}}_{\text{相互作用：只提供修正}}$
$\xi^{(0)}_k=\varepsilon_k-E_F$ 这一项与相互作用无关，它单独就已经给出 $\nu^\propto1/v_F^$、$f_F$ 的结构、以及后面 $1+F_l$ 里那个 $1$。换句话说：朗道参数 $F_l$ 度量的是”相互作用在简并压之上额外加了多少刚度”；$F_l=0$ 不是”零压强”极限，而是”只有简并压”极限。这正是 $s\to1$ 而不是 $s\to0$ 的根源：零声的色散关系里 $F_0^s$ 出现在分母上（$\frac1{F_0^s}=g(s)$），$F_0^s\to0$ 意味着解跑到方程的另一端（$s\to1^+$），而不是跑到 $s\to0$。

(3) 但零声确实会消失——两个极限说的是不同的波。 要把两件事分开：

一声在碰撞弛豫区（$\omega\ll\gamma_{ee}$），分布始终保持局域热平衡，回复力是体弹模量 $\kappa$。它对一切 $F_l$ 都存在，$F_l=0$ 时就是 $v_F/\sqrt3$——纯简并压的声波。

零声在无碰撞区，除了简并压还额外需要排斥反馈 $F_0^s>0$ 把集体模推离粒子-空穴对连续谱的边缘（$s=1$）。$F_0^s\to0^+$ 时 $s\to1^+$，模正好贴上连续谱边缘：极点 $\cos\theta_k=s$ 落进积分区间，阻尼发散、寿命 $\to0$。所以严格说 $F_0^s=0$ 时没有良定义的零声，$s\to1^+$ 只是一个奇异极限。 $\underbrace{u_1=\frac{v_F^*}{\sqrt D}}_{\text{一声：一切 }F_l\text{，纯简并压}}\qquad \underbrace{u_0=\sqrt{\frac{F_0^s}{3}}\,v_F^*}_{\text{零声：需 }F_0^s>0\text{（强相互作用展开）}}$ 一句话：简并压是费米液体的”底色声速”，相互作用是它之上的增量。零声是”底色 + 排斥反馈”的产物，$F_0^s\to0$ 时它退回底色（$s\to1$）并因贴上粒子-空穴对连续谱而失去良定性；而底色本身始终以一声（$v_F/\sqrt3$）的形式存在。

易错点：① 把”$F_l=0$”读成”零压强、零声速”——错，$P=\frac25n\varepsilon_F\neq0$；② 以为”零声与一声是同一支波的连续过渡”——错，一支在碰撞区、一支在无碰撞区，中间隔着 Landau 阻尼区（见上图 3.6.2）；③ 不要把 $u_0=\sqrt{F_0^s/3}\,v_F^*$ 外推到 $F_0^s\to0$——该式来自 $s\gg1$（即 $F_0^s\gg3$）的展开，外推会给出 $u_0\to0$，与色散关系矛盾；$F_0^s\to0^+$ 那一支的正确行为由 $\frac{1}{F_0^s}\simeq\frac12\ln\frac{2}{s-1}$ 给出，即 $s\to1^+$、$u_0\to v_F^*$。

补充：空气、固体、液体的声速与零声是同一个结构——“声速本底 + 相互作用刚度”

把三个看起来毫无关系的公式摆在一起，会发现它们是同一件事：
$\text{理想单原子气体（空气）}\quad c=\sqrt{\frac{\gamma p}{\rho}}=\sqrt{\frac{D+2}{D-2}\frac{k_BT}{m}}\quad\xrightarrow{D=3}\ \sqrt{\frac{5k_BT}{3m}}$ $\text{自由费米气体}\quad u_1=\sqrt{\frac{\kappa}{n}}=\frac{v_F}{\sqrt D}$ $\text{固体}\quad c_L=\sqrt{\frac{B}{\rho}},\qquad c_T=\sqrt\frac{\mu}{\rho}$
共同结构：$c=\sqrt{\text{刚度}/\text{质量密度}}$。而”刚度”都可以拆成本底（无相互作用时粒子的统计动能提供的回复力）$+$ 相互作用增量（把粒子拉在一起的力使体系更难压缩）。

(1) 三种”本底”分别是谁。

空气：本底是热运动 $k_BT$。理想气体毫无相互作用，照样有 $c=\sqrt{5k_BT/3m}$——声波是分子热运动的宏观表现。

自由费米气体：本底是简并压 $\frac{2}{D}n\varepsilon_F$（见上一段补充），量级只由 $v_F$ 决定、与温度无关：$u_1=v_F/\sqrt D$。

固体：本底仍只有费米统计那一项（金属里 $\frac{2}{5}n\varepsilon_F\sim10^{5}$ Pa 量级），而实际固体 $B\sim10^{10}$ Pa——相互作用贡献压倒性地大（差五个数量级）。而且固体的剪切模量 $\mu$ 完全来自键合：把相互作用关掉，”固体”就散成气体，$c_T$ 直接消失（$c_L$ 只剩 $\sqrt{B_{\rm 本底}/\rho}$）。

(2) 相互作用越强（越排斥），声速越大。 这是三者共同的规律，朗道参数把它写成了乘法：
$u_1^2=\frac{\kappa}{n}=\frac{n^2/\chi_c}{n}=\frac{n(1+F_0^s)}{\nu^*}=\frac{v_F^{*2}}{D}\,(1+F_0^s)$
即 $F_0^s$ 就是”在简并压之上再加的相对刚度”（$F_0^s>0$ 排斥使流体更难压缩；液氦-3 的 $F_0^s=10.8$ ⟹ 刚度增强约 11 倍）。零声则把它直接写进速度：
$u_0=sv_F^*,\qquad s\ \text{由}\ \frac{1}{F_0^s}=-1+\frac{s}{2}\ln\frac{s+1}{s-1}\ \text{单调决定}\ \Longrightarrow\ F_0^s\uparrow\Rightarrow s\uparrow\Rightarrow u_0\uparrow$
空气里同一个规律由范德瓦耳斯方程给出：$p=nk_BT(1+b\bar n-a\bar n^2)$，故
$c^2=\frac{\gamma p}{\rho}\ \propto\ \frac{1+2b\bar n-a\bar n^2}{1+b\bar n-a\bar n^2}\cdot\frac{k_BT}{m}$
排斥（$b$，短程硬核）抬高分子有效浓度 ⟹ $c$ 变大；吸引（$a$，范德华）压低它。固体中那个”$b$ 项”就是键合。

(3) 反过来，吸引/软化使声速下降，并在相变点归零。 这是三种体系共有的”软化”信号：
$c^2=\frac{n}{\chi}\quad\Longrightarrow\quad c\to0\ \text{当且仅当}\ \chi\to\infty\ (\text{负压缩系数})$
气体到临界点、液体到临界点、固体到结构相变、费米液体到 $F_0^s\to-1$（零声则对应 $F_0^s\to0^+$ 时模贴上粒子-空穴对连续谱边缘），都是同一个式子在起作用。声速是”软硬度”的量度计。

(4) 定量对照表（3D，除固体外 $\rho=n$）

体系声速本底（无相互作用）相互作用的效果备注

空气（理想气体） $\sqrt{\frac{5k_BT}{3m}}\approx340$ m/s $b\uparrow\Rightarrow c\uparrow$ 本底是热运动

自由费米气体 $\frac{v_F}{\sqrt D}=\sqrt{\frac23\varepsilon_F}$ $F_0^s>0\Rightarrow u_1\uparrow$、$u_0\uparrow$ 本底是简并压

弱相互作用液体（$^3$He） $v_F/\sqrt3\approx96$ m/s $F_0^s=10.8$ ⟹ $u_1\approx2.4\times10^2$ m/s $u_0\approx1.9\times10^2$ m/s，与 $u_1$ 同量级

固态金属/固体 $B_{\rm 本底}/\rho$ 可忽略键合主导：$c_L\sim10^4$ m/s，$c_T\ne0$ 相互作用贡献 ~99.999%

一句话：理想气体靠热运动、费米气体靠简并压，两者都是”无相互作用时的声速本底”；一切真实物质都在这个本底之上叠加”相互作用刚度”，排斥越强刚度越大、声速越快（空气的硬核、液体的 $F_0^s$、固体的键合是同一件事的三个名字），而零声是费米液体里这个结构在无碰撞区的化身——它的速度 $u_0=sv_F^*$ 直接由 $F_0^s$ 决定。

易错点：① 把”无相互作用”当成”无声速”——空气和理想费米气体都是反例，本底来自统计而非相互作用；② 忽略 $c_T$ 的存在意义：它只能由相互作用产生，是固体键合刚度的直接探针；③ 以为 $u_0=\sqrt{F_0^s/3}\,v_F^*$ 可以外推到弱相互作用——它来自 $s\gg1$（即 $F_0^s\gg3$）的展开，$F_0^s\to0^+$ 那一支的正确行为是 $s\to1^+$（见上一段补充的易错点 ③）。

【配图】图 3.6.1：零声在实空间某一时刻的快照。每个空间位置的圆代表该处的费米面形变 $\delta n_k(r,t) \propto -\partial_\xi f\cdot\frac{\cos\theta_k}{s-\cos\theta_k}\cos(qx-\omega t)$；局域流动造成压缩-稀疏，密度涨落产生的反馈场又加速流体，维持流动。

碰撞与一声：准粒子以 $\gamma_{ee} \sim T^2/\varepsilon_F$ 的速率彼此散射。在超大尺度、低频区（$\omega\ll\gamma_{ee}$，$q\ll\gamma_{ee}/v_F$），碰撞使体系始终保持局域平衡，称为流体力学区。此时局域密度矩阵是热的：
$\hat{\rho}(r) = \frac1Z e^{-\frac1{T(r)}(\hat{H}-\mu(r)\hat{N}-\mathbf{u}(r)\cdot\hat{\mathbf{p}})}$
自由度缩减为几个”碰撞零模”（碰撞下守恒的量）：局域温度 $T(r)$、化学势 $\mu(r)$ 和流速 $\mathbf{u}(r)$。忽略温度涨落，动力学由封闭的流体力学方程描述：
$\begin{cases} \partial_t\rho + \nabla\cdot\mathbf{j} = 0\\ (\partial_t + \mathbf{u}\cdot\nabla)\mathbf{u} = -\frac1\rho\nabla P + \frac13(\xi+\nu)\nabla(\nabla\cdot\mathbf{u}) + \nu\nabla^2\mathbf{u} \end{cases}$
第二式是 Navier-Stokes 方程。流体力学方程有密度涨落的声波解，即一声（first sound，类比水中的声波），速度：
$u_1^2 = \frac{\kappa}{\rho} = \frac{v_F^{*2}}{D}(1+F_0^s)(1+F_1^s)$
其中 $\kappa = -V\partial_V P = n^2/\chi_c$ 是体弹模量。$1+F_0^s$ 因子纯粹来自费米液体效应使体弹模量（电荷磁化率的倒数）增强——排斥相互作用使流体更难压缩。弱相互作用极限下 $u_1 = v_F/\sqrt D$ 小于费米速度。因此从低频到高频，声模从一声过渡到零声。

待核：上面 $u_1^2$ 中那个 $(1+F_1^s)$ 因子按理应是分母形式（流的惯性/反流修正），即
$u_1^2=\frac{\kappa}{n}\cdot\frac{1}{1+F_1^s/D}=\frac{v_F^{*2}}{D}\cdot\frac{1+F_0^s}{1+F_1^s/D}$
理由：$u_1^2=\kappa/n=n(1+F_0^s)/\nu^$ 只用压缩性；$F_1^s$ 描述的是准粒子流的惯性（$m^/m=1+F_1^s$，$v_F^/v_F=1/(1+F_1^s)$），它改变的是”流”的惯性而非”压”，故应进分母。用液氦-3 检验：$F_0^s=10.8$、$m^/m\approx2.6$ 时，乘法给 $u_1\approx530$ m/s、分母式给 $u_1\approx270$ m/s，而实验值约 238 m/s（零声 190 m/s）——分母式才对。$F_l\to0$ 时两式都回到 $v_F/\sqrt D$，故弱耦合极限的结论不受影响。请对照讲义确认 $F_1^s$ 的位置（也可能讲义把流动惯性与压缩性合并进了 $(1+F_1^s)$ 的另一种写法）。

【配图】图 3.6.2：中性费米液体中的集体模式。横轴 $q$，纵轴 $\omega$：$q\ll\gamma_{ee}/v_F$ 处为一声（声速 $u_1$），高频处为无阻尼的零声（$u_0 = \sqrt{F_0^s/3}v_F^*$），两者之间是朗道阻尼区；在 $q\sim 2k_F$ 附近是粒子-空穴对激发区。

思考题

单个准电子态 $\Psi^\dagger_k|G\rangle$ 的电荷流是 $k/m$ 还是 $k/m^*$？

如何从动力学方程出发推导流体力学方程，从而理解 $\chi_c = \nu^/(1+F_0^s)$ 和 $m^/m = 1+F_1^s$？（即讲义式 3.6.14）

$F_0^a<0$ 的物理含义是什么？当 $F_0^a\to-1$ 时体系会发生什么？

零声与一声有什么区别？什么时候”声速”大于费米速度？

符号	含义	备注
$k,\ \mathbf k$	准粒子动量，$\lvert\mathbf k\rvert\simeq k_F$	只关心费米面附近
$\hat k$，$\cos\theta_k$	动量方向；与波方向 $\hat q$ 的夹角余弦	全篇唯一的几何变量
$\xi_k$	相对费米能的单准粒子能量	$\xi_{k_F}=0$
$f_F(\xi_k)$	平衡占据（费米-狄拉克分布）	$T\to0$ 时是阶跃函数
$n_k(r,t)$	时空依赖的准粒子占据	朗道理论的动力学变量
$\delta n_k$	小涨落，$n_k=f_F(\xi_k)+\delta n_k$	要解的未知量；取自旋对称通道 $\delta n_{k\sigma}=\delta n_k$
$\varepsilon_{k\sigma}(r,t)=\xi^{(0)}_{k\sigma}+\sum_{k’\sigma’}f_{k\sigma,k’\sigma’}\delta n_{k’\sigma’}$	涨落自己产生的反馈势	扮演”外场”，自洽机制的关键
$f^s_{kk’},\ f^a_{kk’}$	朗道泛函中的准粒子间相互作用	与 $V_{ee}$ 无关，只是被重参数化
$F_l^{s,a}$，$F_{kk’}=\nu^*f^s_{kk’}$	无量纲朗道参数	$l$ 是勒让德展开的阶
$\nu^*$	费米面态密度	$\propto1/v_F^*$，第 4 步靠它换算
$\mathbf v_k^=\partial_k\xi^{(0)}_k=\frac{k}{m^}$	准粒子群速度	$v_F^=k_F/m^$ 是费米面上的值
$q,\ \omega$	涨落的波矢、频率	平面波解
$\alpha_k\ \to\ \eta_{\hat k}$	平面波振幅；费米面上的角分布	二者关系见第 3 步
$-\partial_\xi f_F=\frac{1}{4T}\mathrm{sech}^2(\xi/2T)$	态密度权重	线性理论里涨落的”有效大小”
$s=\frac{\omega}{qv_F^*}$	无量纲波速，即 $u/v_F^*$	主角：$s>1$ 还是 $<1$ 决定生死
$I[n_{k\sigma}]$，$\gamma_{ee}\sim T^2/\varepsilon_F$	碰撞积分、散射率	无碰撞区（$\omega\gg\gamma_{ee}$）置 $I=0$

体系	声速本底（无相互作用）	相互作用的效果	备注
空气（理想气体）	$\sqrt{\frac{5k_BT}{3m}}\approx340$ m/s	$b\uparrow\Rightarrow c\uparrow$	本底是热运动
自由费米气体	$\frac{v_F}{\sqrt D}=\sqrt{\frac23\varepsilon_F}$	$F_0^s>0\Rightarrow u_1\uparrow$、$u_0\uparrow$	本底是简并压
弱相互作用液体（$^3$He）	$v_F/\sqrt3\approx96$ m/s	$F_0^s=10.8$ ⟹ $u_1\approx2.4\times10^2$ m/s	$u_0\approx1.9\times10^2$ m/s，与 $u_1$ 同量级
固态金属/固体	$B_{\rm 本底}/\rho$ 可忽略	键合主导：$c_L\sim10^4$ m/s，$c_T\ne0$	相互作用贡献 ~99.999%

附：第一次作业解答（题 2：费米气中库仑相互作用的一阶微扰）

本题是对相互作用的第一次接触——在费米海上计算库仑相互作用的一阶能量修正，得到 Hartree 与 Fock（交换）两项，为理解朗道费米液体中相互作用的作用做准备。

(a) 库仑相互作用的第二量子化形式

实空间：第一量子化中 $\hat H_{ee}=\frac{\kappa}{2}\sum_{i\neq j}\frac{1}{|\mathbf r_i-\mathbf r_j|}$（$\kappa=e^2$ 为原子/高斯单位，SI 中 $\kappa=e^2/4\pi\varepsilon_0$；因子 $\frac12$ 避免同一对算两遍）。一般二体势 $v(\mathbf r-\mathbf r’)$ 的第二量子化形式为

$\hat H_{ee} = \frac{\kappa}{2}\sum_{\sigma,\sigma'}\int d^3r\,d^3r'\, \frac{\hat\psi_\sigma^\dagger(\mathbf r)\hat\psi_{\sigma'}^\dagger(\mathbf r')\hat\psi_{\sigma'}(\mathbf r')\hat\psi_\sigma(\mathbf r)}{|\mathbf r-\mathbf r'|}$

算符顺序”先湮灭后产生”（从右往左作用：先在 $\mathbf r,\mathbf r’$ 移除两个电子，再放回），因此测量的是两体概率密度。

动量空间：场算符 $\hat\psi_\sigma(\mathbf r)=\frac1{\sqrt\Omega}\sum_ke^{i\mathbf k\cdot\mathbf r}\hat c_{\mathbf k\sigma}$，三维库仑势的傅里叶变换为 $v(\mathbf q)=4\pi\kappa/|\mathbf q|^2$（用收敛因子 $e^{-\eta r}$ 正则化后取 $\eta\to0^+$）。代入并对 $\mathbf r,\mathbf r’$ 积分（$\int_\Omega d^3r\,e^{i(\mathbf k-\mathbf p)\cdot\mathbf r}=\Omega\delta_{\mathbf k,\mathbf p}$），得

$\boxed{\ \hat H_{ee} = \frac{1}{2\Omega}\sum_{\mathbf k,\mathbf k',\mathbf q}\sum_{\sigma,\sigma'} v(\mathbf q)\,\hat c_{\mathbf k+\mathbf q,\sigma}^\dagger\hat c_{\mathbf k'-\mathbf q,\sigma'}^\dagger\hat c_{\mathbf k',\sigma'}\hat c_{\mathbf k,\sigma},\qquad v(\mathbf q)=\frac{4\pi\kappa}{|\mathbf q|^2}\ }$

物理图像：动量 $\mathbf k,\mathbf k’$ 的两个电子散射成 $\mathbf k+\mathbf q,\mathbf k’-\mathbf q$，总动量守恒。

补充：用泊松方程 + 高斯定理算库仑势的傅里叶变换。上面用收敛因子正则化后直接做径向积分，还有一个更优雅的”课上版”做法：把库仑势当拉普拉斯算符的格林函数，用高斯定理绕开积分。

三维：库仑势满足泊松方程
$\nabla^2\frac{1}{|\mathbf r|} = -4\pi\delta(\mathbf r)$
两边做傅里叶变换。右边 $\int d^3r\,e^{-i\mathbf q\cdot\mathbf r}(-4\pi\delta(\mathbf r)) = -4\pi$。左边用高斯定理（散度定理 $\int_V\nabla\cdot\mathbf F\,d^3r=\oint_{\partial V}\mathbf F\cdot d\mathbf S$）把拉普拉斯算符一步步”搬到”平面波上（边界项在 $R\to\infty$ 时消失，可用收敛因子 $e^{-\eta r}$ 严格化）：
$\int d^3r\,e^{-i\mathbf q\cdot\mathbf r}\nabla^2\frac{1}{|\mathbf r|} = \int d^3r\,\frac{1}{|\mathbf r|}\nabla^2e^{-i\mathbf q\cdot\mathbf r} = -q^2\int d^3r\,\frac{e^{-i\mathbf q\cdot\mathbf r}}{|\mathbf r|} = -q^2 v(\mathbf q)$
其中用了 $\nabla^2 e^{-i\mathbf q\cdot\mathbf r} = -q^2e^{-i\mathbf q\cdot\mathbf r}$。两式联立即得
$v(\mathbf q) = \frac{4\pi\kappa}{|\mathbf q|^2}$
高斯定理在这里的用武之地，是证明 $\nabla^2(1/r) = -4\pi\delta$：对半径为 $\varepsilon$ 的小球积分 $\int_{V_\varepsilon}\nabla^2\frac1r\,d^3r = \oint_{S_\varepsilon}\nabla\frac1r\cdot d\mathbf S = \oint\left(-\frac{\hat r}{r^2}\right)\cdot\hat r\,dS = -4\pi$，这正是 $\delta$ 函数的归一化。

二维：两种情况要分清。

二维体系嵌入三维空间（本课程第 3 讲式 3.0.1 的约定）：相互作用仍是 $\kappa/r$（$r$ 为面内距离），其二维傅里叶变换为 $v(\mathbf q) = \int d^2r\,e^{-i\mathbf q\cdot\mathbf r}\frac{\kappa}{r} = 2\pi\kappa\int_0^\infty dr\,J_0(qr) = \frac{2\pi\kappa}{q}$ （用了 $\int_0^\infty J_0(x)dx=1$）。注意这里的 $1/r$ 不是二维拉普拉斯算符的格林函数，所以不能直接套泊松方程 + 高斯定理。一个偷懒的等价做法是”维数约化”：对三维结果沿 $q_z$ 积分 $\int\frac{dq_z}{2\pi}\frac{4\pi\kappa}{q^2+q_z^2} = \frac{2\pi\kappa}{q}$。

真正的二维库仑势（对数势）：二维拉普拉斯算符的格林函数是 $\nabla^2\ln\frac{a}{r} = -2\pi\delta(\mathbf r)$（用二维高斯定理验证：$\oint\nabla\ln(a/r)\cdot d\mathbf l = \oint\left(-\frac{\hat r}{r}\right)\cdot\hat r\,d\ell = -2\pi$），对应的二维泊松方程给 $v(\mathbf q) = \frac{2\pi\kappa}{|\mathbf q|^2}$ 注意它与”嵌入三维的 2D”结果 $2\pi\kappa/q$ 差一个 $1/q$ 的幂——二者物理不同，不要混淆。

(b) 一阶能量密度修正

写 $\hat H(\lambda)=\hat H_0+\lambda\hat H_{ee}$、$|\Psi(\lambda)\rangle=|\Psi_0\rangle+\lambda|\Psi_1\rangle+\cdots$、$E(\lambda)=E^{(0)}+\lambda E^{(1)}+\cdots$，代入薛定谔方程取 $\lambda$ 一阶并左乘 $\langle\Psi_0|$，含 $|\Psi_1\rangle$ 的两项相消，得一阶修正即取期望值：

$E^{(1)} = \langle\Psi_0|\hat H_{ee}|\Psi_0\rangle$

对零温无自旋极化电子气，基态是费米海 $|FS\rangle$（$n_{\mathbf k\sigma}=\Theta(k_F-k)$，$n=k_F^3/3\pi^2$）。关键是四费米算符的期望值：为让末态与费米海有非零重叠，两个湮灭算符制造的空穴必须被两个产生算符补回，只有两种配对：

$\langle\hat c_a^\dagger\hat c_b^\dagger\hat c_c\hat c_d\rangle_0 = n_an_b\left(\delta_{ad}\delta_{bc}-\delta_{ac}\delta_{bd}\right)$

直接配对（$d=a,c=b$）给 $+n_an_b$；交换配对（$c=a,d=b$）因 $\hat c_a\hat c_b=-\hat c_b\hat c_a$ 给 $-n_an_b$——负号来自反对易关系。代回动量指标：

$\langle\hat c_{\mathbf k+\mathbf q,\sigma}^\dagger\hat c_{\mathbf k'-\mathbf q,\sigma'}^\dagger\hat c_{\mathbf k',\sigma'}\hat c_{\mathbf k,\sigma}\rangle_0 = n_{\mathbf k\sigma}n_{\mathbf k'\sigma'}\left(\delta_{\mathbf q,0}-\delta_{\sigma\sigma'}\delta_{\mathbf q,\mathbf k'-\mathbf k}\right)$

于是 $E^{(1)}=E_H+E_x$：

Hartree（直接）项：要求 $\mathbf q=0$，自旋任意，

$E_H = \frac{v(0)}{2\Omega}\sum_{\mathbf k,\mathbf k',\sigma,\sigma'}n_{\mathbf k\sigma}n_{\mathbf k'\sigma'} = \frac{v(0)N^2}{2\Omega}$

$v(0)=4\pi\kappa/q^2|_{q=0}$ 发散——这是全体电子与自身的长程静电自能。均匀正电背景（jellium）使体系电中性：电子-电子、电子-背景、背景-背景三项能量密度系数分别为 $+\frac12n^2,\ -n^2,\ +\frac12n^2$，之和为零，恰好抵消 Hartree 项。物理上 $q=0$ 模必须结合电中性条件处理，不能无解释地删除。

Fock（交换）项：要求 $\mathbf q=\mathbf k’-\mathbf k$ 且 $\sigma=\sigma’$（交换项只发生在同自旋之间），

$E_x = -\frac{1}{2\Omega}\sum_{\mathbf k,\mathbf k',\sigma}v(\mathbf k-\mathbf k')n_{\mathbf k\sigma}n_{\mathbf k'\sigma}$

无自旋极化时 $\sum_\sigma$ 给出因子 2、恰好抵消 $\frac12$。取热力学极限并做积分（交换积分中的 $1/q^2$ 奇点在 3D 可积：$\int_{q<\delta}d^3q/q^2=4\pi\delta<\infty$；先角积分 $8\pi^2$，再令 $p=kx$ 做径向二重积分），得

$\frac{E_x}{\Omega} = -\int_{k<k_F}\frac{d^3k}{(2\pi)^3}\int_{k'<k_F}\frac{d^3k'}{(2\pi)^3}\frac{4\pi\kappa}{|\mathbf k-\mathbf k'|^2} = -\frac{\kappa}{2\pi^3}\int_0^{k_F}\!dk\int_0^{k_F}\!dp\,kp\ln\frac{k+p}{|k-p|} = -\frac{\kappa k_F^4}{4\pi^3}$

最终答案：电中性均匀电子气的一阶能量密度修正就是交换能：

$\boxed{\ \frac{E^{(1)}}{\Omega}=\frac{E_x}{\Omega}=-\frac{\kappa k_F^4}{4\pi^3} = -\frac34\left(\frac3\pi\right)^{1/3}\kappa\,n^{4/3}\ }$

每粒子的交换能：

$\frac{E^{(1)}}{N} = -\frac{3\kappa k_F}{4\pi}$

若再计入无相互作用动能 $E_0/\Omega=\hbar^2k_F^5/(10\pi^2m)$：

$\frac{E}{\Omega}=\frac{\hbar^2k_F^5}{10\pi^2m}-\frac{\kappa k_F^4}{4\pi^3}+O(\kappa^2), \qquad \frac{E}{N}=\frac35E_F-\frac{3\kappa k_F}{4\pi},\quad E_F=\frac{\hbar^2k_F^2}{2m}$

为什么结果是负的？ 库仑势本身始终是排斥的（$\kappa/|\mathbf r-\mathbf r’|>0$）。负的交换能不是电子间相互作用变吸引，而是费米统计的量子效应：同自旋电子因波函数反对称形成”交换洞”（exchange hole），彼此靠近的概率被压低，平均库仑排斥能低于经典平均场估计。Fock 项的负号来自反对易关系，不是库仑力变吸引。

易错清单：① 忘了 $\frac12$（二体项会把电子对算两遍）；② 忘了自旋指标（交换项只在同自旋）；③ 把 $\hat\psi(\mathbf r)$ 当普通波函数（它是场算符）；④ 无解释删掉 $q=0$ 项（正确理由是电中性背景抵消）；⑤ 把交换能负号解释为吸引相互作用（错）；⑥ 混淆能量密度 $E/\Omega$ 与每粒子能量 $E/N$（$E=\Omega n\varepsilon$）。