线性响应理论基础笔记#

0. 这篇笔记要做什么#

线性响应理论（Linear Response Theory, LRT） 是量子统计物理中最优雅也最实用的框架之一。它的核心问题极其简单：

给一个处于平衡态的系统施加一个非常弱的外场，系统会如何响应？

这个问题驱动了 Kubo 公式、涨落耗散定理、格林函数方法等一系列工具的发展。本笔记从大四物理学生的视角出发，假设你学过量子力学和统计物理，但不假设你熟悉多体理论、场论或高级凝聚态。

本文与前人综述的关键区别：不罗列结论，而是展示"每个概念为什么被需要"——从问题出发，逐步构建理论工具。

1. 线性响应理论的出发点：小扰动下系统如何响应#

1.1 核心问题#

想象一个物理系统——可以是一块金属、一个超导量子比特、或者一个分子——初始时处于热平衡。现在施加一个非常弱的外场（电场、磁场、激光……）。你测量某个可观测量的变化时，希望只保留线性于外场强度的项。

这之所以可行，是因为一个深刻的物理事实：

1.2 为什么要关注"线性"响应？#

三个原因：

1. 普适性：线性响应系数（如电导率 $\sigma$、磁化率 $\chi$）是材料的内禀性质，不依赖于外场的具体细节
2. 可测性：实验中测量的几乎都是线性响应——交流电导率、中子散射截面、核磁共振谱线
3. 理论基础：线性响应是通往非平衡统计力学的门户——理解"系统离平衡不远时如何行为"，是理解更复杂非平衡现象的前提

1.3 经典类比：阻尼谐振子#

在进入量子力学之前，先看一个经典例子。考虑受迫阻尼谐振子：

当 $F(t)$ 很小时，稳态解为 $x(t) = \int_{-\infty}^t \chi(t - t') F(t') dt'$，其中 $\chi(t) = \frac{1}{m\omega_1} e^{-\gamma t/2} \sin(\omega_1 t)$ （$\omega_1 = \sqrt{\omega_0^2 - \gamma^2/4}$）是响应函数。

关键特征：

• 因果性：$x(t)$ 只依赖于过去的力（$t' \le t$）→ $\chi(t) = 0$ for $t < 0$
• 卷积形式：$\chi(t - t')$ 只依赖于时间差 → 时间平移不变性（平衡态系统的特征）
• 频域简洁性：Fourier 变换后 $\tilde{x}(\omega) = \tilde{\chi}(\omega) \tilde{F}(\omega)$

量子线性响应的形式完全平行，只是 $\chi(t)$ 现在由系统的量子力学算符决定。

1.4 适用条件#

• 外场强度足够小（$|\Phi| \ll 1$，系统的特征能量尺度）
• 系统初始处于稳态（通常是热平衡态 $\rho_0 \sim e^{-\beta H_0}$）
• 外场是"经典"的（不是量子化的）——即用一个含时函数 $F(t)$ 而非量子算符来描述

2. 从经典响应到量子响应：为什么需要响应函数#

2.1 量子系统中的"响应"是什么？#

在经典物理中，响应函数联系力 $F(t)$ 和位移 $x(t)$。在量子力学中，我们关心的不是轨道，而是期望值的变化：

其中：

• $A$ 是我们要测量的可观测量（如电流 $\hat{J}$、磁化强度 $\hat{M}$）
• $\rho_0$ 是初始平衡态密度矩阵
• $\rho(t)$ 是外场作用下 $t$ 时刻的密度矩阵

2.2 量子响应函数的定义#

与经典类比，量子线性响应函数 $\phi_{AV}(t)$ 由以下关系定义：

这里 $F(t)$ 是外场强度（标量函数），$V$ 是与外场耦合的系统算符——即微扰哈密顿量为 $H_{\mathrm{ext}}(t) = F(t) V$（记号约定见下方 callout）。

2.3 常见误区#

3. 含时微扰论回顾：外场如何进入哈密顿量#

3.1 设定#

总哈密顿量：

其中 $H_0$ 是未微扰哈密顿量（已知其本征态和本征值），$H_{\mathrm{ext}}(t)$ 描述外场与系统的耦合。

3.2 相互作用表象：消除"自由演化"的干扰#

在薛定谔绘景中，密度矩阵按 Liouville-von Neumann 方程演化：

我们关心的是外场 $H_{\mathrm{ext}}$ 引起的偏离，而 $H_0$ 本身的自由演化在计算响应时需要和 $H_{\mathrm{ext}}$ 的效应分离开。为此转到相互作用表象（Dirac 表象）。对任意算符 $O$：

变换的目的是把 $H_0$ 吸收到算符的"自由转动"中，让密度矩阵只感受 $H_{\mathrm{ext}}$。直接对 $\rho_I(t) = e^{iH_0 t/\hbar} \rho(t) e^{-iH_0 t/\hbar}$ 求导：

其中 $H_{\mathrm{ext},I}(t) = e^{iH_0 t/\hbar} H_{\mathrm{ext}}(t) e^{-iH_0 t/\hbar} = F(t) V_I(t)$。注意：$H_0$ 项消掉了——这正是相互作用表象的设计目的。

3.3 一阶微扰展开#

将 Liouville-von Neumann 方程从初始时刻 $t_0$ 积分到 $t$：

这仍然是精确的（和原微分方程等价）——右边积分号内仍含有 $\rho_I(t')$。一阶微扰的关键近似：因为外场很弱（$F$ 是小量），$\rho_I(t')$ 与 $\rho_I(t_0)$ 的偏差是 $O(F)$ 量级。把 $\rho_I(t')$ 在积分号内替换为 $\rho_I(t_0)$ 只会引入 $O(F^2)$ 的误差——而线性响应本来就要丢弃 $O(F^2)$ 及更高阶项。因此：

这个近似等价于只保留 Dyson 级数的一阶项。它成立的条件是：$|F| \times (\text{微扰作用时间}) \ll \hbar$（具体地说，$\int_{t_0}^t |F(t')| dt' \ll \hbar / \|V\|$）。

3.4 期望值变化#

可观测量 $A$ 的期望值变化：

代入 $H_{\mathrm{ext},I}(t') = F(t') V_I(t')$：

3.5 关键假设#

1. 初始稳态：$[\rho_0, H_0] = 0$（平衡态与 $H_0$ 对易）→ 可推到第 3 行
2. 外场弱：仅保留 $F$ 的一阶项
3. $t_0 \to -\infty$：外场从无穷过去绝热开启（Kubo 的标准设定）

取 $t_0 \to -\infty$，利用关联函数的时移不变性 $\langle [A_I(t), V_I(t')] \rangle_0 = \langle [A_I(t - t'), V] \rangle_0$：

这是量子响应函数的最基本形式。$\theta(t)$ 是 Heaviside 阶跃函数（保证因果性——响应不能先于微扰）。

4. Kubo 公式的详细推导#

4.1 什么是 Kubo 公式？#

Kubo 公式是上节结果的频域版本。频域更方便，因为：

• 实验通常用单色（单频）外场
• 响应函数在频域中的解析性质与因果性直接相关
• 频域关联函数与谱函数的联系更直接

4.2 频域转换#

第一步：将 $F(t')$ 的 Fourier 表示代入卷积公式，并交换积分顺序（假定积分绝对收敛或理解为分布意义下的操作）：

第二步：令 $\tau = t - t'$（$dt' = -d\tau$，当 $t' = -\infty$ 时 $\tau = +\infty$，当 $t' = t$ 时 $\tau = 0$）：

代入上式：

第三步：第二个积分只依赖于 $\omega$ 而不再依赖于 $t$——这正是我们需要的。定义广义磁化率（susceptibility）：

则 $\Delta \langle A \rangle(t) = \frac{1}{2\pi} \int d\omega \, \chi_{AV}(\omega) \tilde{F}(\omega) e^{-i\omega t}$。这就是频域 Kubo 公式的基本形式：时域的因果卷积 = 频域的乘积——和经典阻尼谐振子完全平行（§1.3）。

4.3 从实时间响应函数到谱分解形式#

第一步：将 §3.5 的 $\phi_{AV}(t)$ 代入 $\chi_{AV}(\omega)$ 的定义（$\theta(t)$ 已被积分限 $[0, \infty)$ 吸收）：

第二步：在 $H_0$ 的本征基 $\{|m\rangle\}$（$H_0|m\rangle = E_m|m\rangle$）下展开迹。$\rho_0 = e^{-\beta H_0}/Z$ 在此基下是对角的：

展开对易于的迹（利用迹的循环性将 $\rho_0$ 移到最右边）：

利用 $\rho_0|m\rangle = (e^{-\beta E_m}/Z) |m\rangle$，以及 $A_I(t) = e^{iH_0 t/\hbar} A e^{-iH_0 t/\hbar}$ 给出 $\langle m|A_I(t)|n\rangle = \langle m|A|n\rangle e^{i(E_m - E_n)t/\hbar}$：

第三步：代入 (4.1)，交换 $m \leftrightarrow n$ 于第二项中，合并：

第四步：时间积分。被积函数是 $e^{i\Omega t}$ 形式（$\Omega = \omega + (E_m - E_n)/\hbar$），在 $t \to \infty$ 时振荡而不衰减——积分的上限不收敛。物理上，这是因为我们假设外场从 $t = -\infty$ 起一直存在；实际上 Kubo 考虑的是外场从遥远过去缓慢绝热开启，这等价于给频率加一个无穷小的虚部：$\omega \to \omega + i\eta$（$\eta \to 0^+$）。这个 $i\eta$ 在时域中对应衰减因子 $e^{-\eta t}$，保证积分收敛：

代入 $\Omega = \omega + (E_m - E_n)/\hbar$：

最后一步分子分母同乘 $\hbar$。把这个结果代回 $\chi_{AV}(\omega)$ 的表达式：

$-\frac{i}{\hbar}$ 与 $i\hbar$ 恰好消去（$-\frac{i}{\hbar} \cdot i\hbar = 1$）。剩下的 $\hbar\eta$ 量纲是能量且 $\hbar\eta \to 0^+$，将它重命名为 $\eta$（就是一个趋于 $0^+$ 的能量量纲无穷小量），得到 Kubo 公式的谱分解形式：

其中分母中的 $+i\eta$（$\eta \to 0^+$）来自时间积分的收敛因子——它的物理作用是选择推迟边界条件（响应只能出现在微扰之后），数学作用是告诉我们在遇到极点 $E_n = E_m + \hbar\omega$ 时应该从复平面的上半平面绕过。

4.4 推导中的关键步骤#

1. $\theta(t)$ 的 Fourier 变换：严格地说，为了让积分收敛，需要将外场看成从过去缓慢打开，引入收敛因子：

   $$ \int_0^{\infty} dt \, e^{i\omega t - 0^+ t} = \frac{i}{\omega + i0^+} $$

   其中 $0^+$ 是正无穷小量。它决定了极点从复平面的上半平面绕过——正是这个 $i0^+$ 选择了"推迟响应"（因果响应）。下面把所有 $i0^+$ 简写为 $i\eta$（$\eta \to 0^+$，不是物理耗散率，只是数学工具）。

2. 迹的展开：$\operatorname{Tr}[[A_I(t), V] \rho_0] = \sum_{m,n} (e^{-\beta E_m} - e^{-\beta E_n}) \langle m|A|n\rangle \langle n|V|m\rangle e^{i(E_n - E_m)t/\hbar}/Z$

3. 时间积分：$\int_0^{\infty} dt e^{i(\omega + (E_n - E_m)/\hbar)t - \eta t} = i\hbar/(E_m - E_n + \hbar\omega + i\hbar\eta)$

4.5 物理意义#

• 分母 $E_m - E_n + \hbar\omega$：外场频率 $\omega$ 必须匹配系统的能级差 $\Delta E_{nm} = E_n - E_m$ 才能产生显著响应 → 共振条件
• 权重 $(e^{-\beta E_m} - e^{-\beta E_n})$：只有不同布居数的态之间才有净响应（平衡态下布居数相同 $\Rightarrow$ 无净跃迁）
• 矩阵元 $\langle m|A|n\rangle \langle n|V|m\rangle$：跃迁的选择定则——$A$ 和 $V$ 都必须能连接 $|m\rangle$ 和 $|n\rangle$

4.6 完整计算示例：两能级系统的动态磁化率#

理论最怕空洞。让我们用 Kubo 公式实际算一个最简单的量子系统——两能级系统（TLS）在外场下的动态磁化率。

模型：$H_0 = \frac{\hbar\omega_0}{2} \sigma_z$，基态 $|g\rangle$（$E_g = -\hbar\omega_0/2$），激发态 $|e\rangle$（$E_e = +\hbar\omega_0/2$）。

外场耦合：外磁场 $B(t)$ 沿 $x$ 方向，$H_{\mathrm{ext}}(t) = -B(t) \mu \sigma_x$，即 $F(t) = B(t)$，$V = \mu \sigma_x$。我们测量 $A = \mu \sigma_x$（磁矩的 $x$ 分量）。

第一步：计算矩阵元。在 $\sigma_z$ 的本征基下 $\sigma_x = |e\rangle\langle g| + |g\rangle\langle e|$，非零矩阵元只有 $\langle e|\sigma_x|g\rangle = \langle g|\sigma_x|e\rangle = 1$。

第二步：写 Kubo 公式的谱分解。代入 $A = V = \mu\sigma_x$，$Z = 1 + e^{-\beta\hbar\omega_0}$：

$(m,n) = (g,e)$ 项：$\Delta E = E_g - E_e = -\hbar\omega_0$，权重 $= (p_g - p_e)$ $(m,n) = (e,g)$ 项：$\Delta E = E_e - E_g = +\hbar\omega_0$，权重 $= (p_e - p_g)$

其中 $p_g = 1/(1+e^{-\beta\hbar\omega_0})$，$p_e = e^{-\beta\hbar\omega_0}/(1+e^{-\beta\hbar\omega_0})$，且 $p_g - p_e = \tanh(\beta\hbar\omega_0/2)$。

第三步：整理：

第四步：取虚部（吸收谱）。需要用到 Sokhotski-Plemelj 公式（也叫 Plemelj 公式）：

其中 $\mathcal{P}$ 表示 Cauchy 主值。取虚部即得 $\Im[(x + i\eta)^{-1}] = -\pi\delta(x)$。这个公式的来源是：把 $1/(x + i\eta)$ 的实部和虚部分离：

实部的 $\frac{x}{x^2 + \eta^2}$ 在 $\eta \to 0^+$ 时趋近于 $\mathcal{P}(1/x)$（主值——在 $x = 0$ 处对称挖掉奇点）；虚部的 $\frac{\eta}{x^2 + \eta^2}$ 是一个宽度 $\sim \eta$ 的 Lorentzian，面积恒为 $\pi$，在 $\eta \to 0^+$ 时趋近于 $\pi\delta(x)$。合起来就是上式。

现在把 $\chi_{xx}(\omega)$ 的两个分式分别应用这个公式：

取虚部（丢掉实部的主值项），代入 $\chi_{xx}$：

物理诠释：

• 吸收峰在 $\omega = \pm\omega_0$：外场频率必须匹配 TLS 的能级差才能被吸收——共振条件是 Kubo 公式分母 $\hbar\omega + E_m - E_n$ 的必然结果
• 负 $\omega$ 的峰对应受激发射：系统向外场释放能量（$\omega < 0$ 意味着系统从高能级跃迁到低能级，把能量交给外场）——这是量子系统与经典谐振子的关键区别
• 温度因子 $\tanh(\beta\hbar\omega_0/2)$：高温下 $p_g \approx p_e$ → 吸收和发射几乎抵消 → 净响应趋近于零。物理上：两能级在高温下"饱和"了，对外场不敏感
• 零温极限：$\tanh \to 1$ → 只有基态布居 → 只有吸收峰（$\omega = +\omega_0$），无发射——因为 $T=0$ 时系统没有能量可以释放

验证 FDT：计算对称关联函数 $C_{xx}(t) = \frac{\mu^2}{2} \langle \{\sigma_x(t), \sigma_x(0)\} \rangle_0$。利用 $\sigma_x(t) = \sigma_+ e^{i\omega_0 t} + \sigma_- e^{-i\omega_0 t}$ 和 $\operatorname{Tr}[\sigma_+ \sigma_- \rho_0] = p_g$，得：

代入 FDT $\Im[\chi_{xx}] = -\frac{1}{\hbar}\tanh(\beta\hbar\omega/2) C_{xx}$，利用 $\delta$ 函数的性质 $\tanh(\beta\hbar\omega/2)|_{\omega = \pm\omega_0} = \pm \tanh(\beta\hbar\omega_0/2)$，等式恰好成立。

实部（色散）：$\Re[\chi_{xx}(\omega)] = \frac{\mu^2}{\hbar} \tanh(\beta\hbar\omega_0/2) \left[\mathcal{P}\frac{1}{\omega - \omega_0} - \mathcal{P}\frac{1}{\omega + \omega_0}\right]$——在 $\omega = \pm\omega_0$ 处有奇点，远离共振时 $|\Re[\chi]| \sim 1/|\omega - \omega_0|$。实验上测量的是折射率 $\propto \Re[\chi]$ 和吸收系数 $\propto \Im[\chi]$。

5. 因果性、推迟响应函数与对易子#

5.1 为什么响应函数必须是对易子？#

回到实时间响应函数：

对易子 $[A_I(t), V]$ 的来源：它是从 $\operatorname{Tr}[[A_I(t), V] \rho_0] = \operatorname{Tr}[A_I(t) [V, \rho_0]]$ 这条路径来的。物理上，对易子衡量的是 $t$ 时刻的 $A$ 与 $0$ 时刻的 $V$ 之间的量子非对易性——这正是"外场 $V$ 施加后，系统通过量子动力学将扰动传播到 $A$"的数学表达。

5.2 推迟响应函数#

上节定义的 $\phi_{AV}(t)$ 就是推迟响应函数——"推迟"意味着 $A$ 在 $V$ 之后才响应（$t > 0$ 时才有非零值）：

许多多体教材会定义不含 $1/\hbar$ 的推迟格林函数：

两者只差一个约定因子：$G^R_{AV} = \hbar \cdot \phi_{AV}$。本文统一使用 $\phi_{AV}(t)$ 表示响应核，$\chi_{AV}(\omega)$ 表示它的频域形式（见 §4.2），不再单独引入 $G^R$ 以免混淆。

5.3 因果性与 Kramers-Kronig 关系#

为什么因果性 $\Rightarrow$ 上半平面解析？#

回忆 $\chi_{AV}(\omega)$ 的定义（§4.2）：$\chi_{AV}(\omega) = \int_0^{\infty} dt \, \phi_{AV}(t) e^{i\omega t}$。把 $\omega$ 推广到复数 $\omega = \omega_R + i\omega_I$：

当 $\omega_I > 0$（上半平面）时，$e^{-\omega_I t}$ 是指数衰减的——积分的收敛性比实轴上更好。而且 $\chi_{AV}(\omega)$ 在上半平面的任意闭区域内可以逐点求导（因为 $e^{i\omega t} = e^{i\omega_R t - \omega_I t}$ 对 $\omega$ 解析，积分收敛一致）。因此 $\chi_{AV}(\omega)$ 在 $\Im[\omega] > 0$ 内解析。

如果 $\omega_I < 0$（下半平面），$e^{-\omega_I t}$ 是指数增长的，积分可能发散——$\chi_{AV}(\omega)$ 在下半平面可能有奇点。解析性仅在上半平面有保证，这是 $\phi_{AV}(t)$ 只在 $t \ge 0$ 非零（因果性）的直接后果。

从 Cauchy 积分公式到 K-K 关系#

由于 $\chi_{AV}(z)$ 在上半平面解析，对上半平面内任意一点 $\omega + i\epsilon$（$\epsilon > 0$），Cauchy 积分公式给出：

取闭合围道 $C$：沿实轴从 $-R$ 到 $+R$，再绕上半平面的大半圆回到 $-R$。当 $R \to \infty$ 时，大半圆上的积分趋于零（因为 $|\chi_{AV}(z)| \to 0$ 当 $|z| \to \infty$，$\Im[z] > 0$——这是 $\phi_{AV}(t)$ 在 $t=0$ 连续的推论）。剩下实轴积分：

令 $\epsilon \to 0^+$，利用 §4.6 的 Sokhotski-Plemelj 公式 $\lim_{\epsilon \to 0^+} \frac{1}{\omega' - \omega - i\epsilon} = \mathcal{P}\frac{1}{\omega' - \omega} + i\pi\delta(\omega' - \omega)$：

移项得 $\chi_{AV}(\omega) = \frac{1}{\pi i} \mathcal{P} \int_{-\infty}^{\infty} \frac{\chi_{AV}(\omega')}{\omega' - \omega} d\omega'$。分离实部和虚部：

物理含义：响应函数的实部（色散/折射）和虚部（耗散/吸收）不是独立的。知道所有频率的吸收谱 $\Im[\chi(\omega)]$，就能通过第一个公式算出任意频率的色散 $\Re[\chi(\omega)]$——无需再做实验。实验上通常测量 $\Im[\chi]$（吸收谱），再通过 K-K 关系积分得到 $\Re[\chi]$（折射率）。

记忆法：两个公式差一个负号。第一个公式 $\Re$ 由 $\Im$ 的积分给出（正号）；第二个公式 $\Im$ 由 $\Re$ 的积分给出（负号）。

5.4 常见误区#

• "响应函数就是对易子" → 不完全。响应函数包含 $\theta(t)$（因果性）和 $(-i/\hbar)$（来自 Schrödinger 方程）
• "格林函数就是响应函数" → 有多种格林函数（推迟、超前、因果、热……），推迟格林函数才对应物理响应
• "K-K 关系对所有函数成立" → 仅当函数在上半平面解析（等价于因果性），对无因果性的量（如关联函数）不直接成立

6. 响应函数、关联函数和格林函数的关系#

6.1 什么是关联函数？#

关联函数衡量的是两个可观测量在不同时刻的量子涨落相关性：

其中 $\{X, Y\} = XY + YX$ 是反对易子（anticommutator）。

在经典极限（$\hbar \to 0$）下：

• 对易子 $\frac{1}{i\hbar}[A, B]$ 过渡到经典泊松括号 $\{A, B\}_{\mathrm{PB}}$。因此响应并非"只有量子系统才有"——经典系统同样有线性响应，只是量子理论用对易子表达，经典理论用泊松括号表达
• 反对易子 $\frac{1}{2}\{A, B\}$ 在 $\hbar \to 0$ 时直接过渡为经典关联 $A(t)B(0)$——这是纯经典可测的量

6.2 三种函数的关系图#

核心关系（热平衡下）：

这是涨落耗散定理的最基本形式。

6.3 常见误区#

6.4 Onsager 倒易关系：响应函数的对称性#

Lars Onsager 在 1931 年获得诺贝尔奖的工作证明了一个深刻的结论：线性输运系数矩阵是对称的。在 Kubo 公式的语言中，这等价于：

其中 $B$ 代表外加磁场，$\varepsilon_A = \pm 1$ 是 $A$ 在时间反演下的宇称（如位置是偶、速度是奇）。

物理含义：如果你测量"用力 $F_V$ 推系统，观测 $A$ 如何响应"，和"用力 $F_A$ 推系统，观测 $V$ 如何响应"，两者在磁场反转下给出相同的结果。例如：电场产生温度梯度（Peltier 效应）和温度梯度产生电流（Seebeck 效应）的系数互为倒易。

从 Kubo 公式验证 Onsager 关系：在时间反演下 $\Theta H_0(B) \Theta^{-1} = H_0(-B)$，$\Theta \rho_0 \Theta^{-1} = \rho_0$（平衡态时间反演不变），$\Theta A \Theta^{-1} = \varepsilon_A A$。将这些代入 $\chi_{AV}$ 的定义，利用 $\langle [A(t), V] \rangle_0 = \varepsilon_A \varepsilon_V \langle [V(t), A] \rangle_0 |_{B \to -B}$，即可得证。

对 Berry 相研究的直接意义：Onsager 关系保证了 $\sigma_{xy}(B) = \sigma_{yx}(-B)$——霍尔电导率在磁场反转下是反对称的（$\sigma_{xy} = -\sigma_{yx}$）。你的 Berry 相/几何量子门论文中处理的参数空间，本质上具有和磁场类似的"时间反演破缺"结构——拓扑响应的反对称性来自同一个数学根源。

7. 谱表示与谱函数#

7.1 为什么需要谱函数？#

Kubo 公式涉及对能级求和。对于大系统（多体、连续谱），直接求和不可行。谱函数将这些求和打包成一个频率函数。先讨论自关联情形（$A = V$）：

$S_{AA}(\omega)$ 的物理含义：系统在频率 $\omega$ 处吸收或发射能量 $\hbar\omega$ 的跃迁强度权重。对于两个不同算符 $A$ 和 $V$ 的交叉谱函数，矩阵元一般是 $\langle m|A|n\rangle \langle n|V|m\rangle$，不能简单写成模平方。

7.2 与响应函数的关系#

这个关系不是额外的假设——直接从 Kubo 公式的谱分解形式和谱函数的定义推出来。

从 §4.3 的 Kubo 谱分解出发（令 $A = V$，用 §7.1 的记号 $S_{AA}$ 表示自关联谱函数）：

取虚部。对分母用 Sokhotski-Plemelj 公式（§4.6）：$\Im[(E_m - E_n + \hbar\omega + i\eta)^{-1}] = -\pi\delta(E_m - E_n + \hbar\omega)$。注意 $\delta$ 函数在 $\hbar\omega = E_n - E_m$ 时非零：

将分子拆成两项。第一项（含 $e^{-\beta E_m}$）的 $\delta$ 条件为 $\hbar\omega = E_n - E_m$；第二项（含 $e^{-\beta E_n}$）交换哑指标 $m \leftrightarrow n$，$\delta$ 变为 $\delta(\hbar\omega - (E_m - E_n)) = \delta(\hbar\omega - (E_n - E_m))$（$\delta$ 是偶函数），得：

第二个因子正好是 §7.1 定义的谱函数 $S_{AA}(\omega)$。因此：

物理含义：

• $\Im[\chi_{AA}] < 0$（$\omega > 0$，吸收）：因子 $(1 - e^{-\beta\hbar\omega}) > 0$，$S_{AA} > 0$ → 负号来自 $\chi$ 的约定
• $\omega < 0$（发射）：$1 - e^{-\beta\hbar\omega} < 0$ → $\Im[\chi_{AA}]$ 变号 → 受激发射对应 $\Im[\chi] > 0$
• 实验上：$\Im[\chi_{AA}]$ 是吸收系数（光学吸收谱 $\alpha(\omega) \propto \omega \Im[\chi]$），$S_{AA}$ 是非弹性散射截面（中子散射直接测 $S(\mathbf{q}, \omega)$）。两者通过温度因子 $(1 - e^{-\beta\hbar\omega})$ 联系——同一个物理信息的不同表示

7.3 谱函数与态密度的区别#

7.4 动态结构因子#

在凝聚态和中子散射中，常用动态结构因子（dynamical structure factor）：

其中 $\rho_{\mathbf{q}} = \sum_j e^{i\mathbf{q} \cdot \mathbf{r}_j}$ 是密度算符的 Fourier 分量。

它与响应函数的关系（涨落耗散定理的另一种形式）：

实验意义：$S(\mathbf{q}, \omega)$ 直接由非弹性中子散射或 X 射线散射测量 → 从中提取 $\Im[\chi]$ → 通过 K-K 关系得 $\Re[\chi]$ → 完整响应函数。

7.5 Lehmann 表示：从 Kubo 到谱函数的形式桥梁#

将谱分解形式的 Kubo 公式进一步写开，可以得到Lehmann 表示——这是连接 Kubo 公式和谱函数的标准形式。将 $\chi_{AV}(\omega)$ 的虚部显式写出：

当 $A = V$ 且 $A$ 为厄米算符时，在本文约定下，正频部分（$\omega > 0$）的 $\Im[\chi_{AA}(\omega)]$ 为非正。这对应外场向系统输入能量（吸收）。负频部分由 $\Im[\chi_{AA}(-\omega)] = -\Im[\chi_{AA}(\omega)]$（奇对称性）给出，符号相反——对应受激发射。整体上，$(e^{-\beta E_m} - e^{-\beta E_n})$ 与 $\delta(\hbar\omega - (E_n - E_m))$ 中隐含的因子共同保证了净耗散功率为非负（热力学第二定律在 LRT 层面的体现）。

Lehmann 表示的物理图像：想象把所有能级画在一条轴上，每个能级对的跃迁对应一条从 $E_m$ 到 $E_n$ 的"箭头"。$\Im[\chi_{AA}(\omega)]$ 在 $\omega$ 处的值，就是所有"能量差恰好等于 $\hbar\omega$ 的箭头"的加权和。低温下，从低能级出发的箭头（$E_m$ 小 → $e^{-\beta E_m}$ 大）权重更大。

与实验的直接对应：

• 光学吸收谱 $\alpha(\omega) \propto \omega \Im[\chi_{xx}(\omega)]$（电偶极跃迁，$A = V = e\hat{x}$）
• 非弹性中子散射截面 $\frac{d^2\sigma}{d\Omega d\omega} \propto S(\mathbf{q}, \omega) \propto \frac{1}{1-e^{-\beta\hbar\omega}} \Im[\chi_{\rho\rho}(\mathbf{q}, \omega)]$
• 核磁共振（NMR）线形 $\propto \Im[\chi_{S_+ S_-}(\omega)]$（自旋翻转跃迁，见 §13.5）

常见误区：$\Im[\chi]$ 和 $S(\omega)$ 只差一个温度因子 $(1 - e^{-\beta\hbar\omega})$——实验上测量的是 $S(\omega)$（散射截面），但从 $S(\omega)$ 提取 $\Im[\chi]$ 需要除以这个因子。在 $\hbar\omega \ll k_B T$ 时这个因子趋近于零 → 低频 $S(\omega)$ 信号很强，但对应的 $\Im[\chi]$ 可能被放大若干量级 → 数据处理中的常见陷阱。

7.6 常见实验技术及其对应的响应函数#

注意：ARPES 测量的是单粒子谱函数 $A(\mathbf{k}, \omega) = -\frac{1}{\pi} \Im[G^R(\mathbf{k}, \omega)]$（推迟格林函数的虚部），而中子散射和光学测量的是两粒子谱函数（密度-密度或电流-电流关联）。两者通过不同的格林函数与 LRT 联系——单粒子响应涉及 $G^R \sim \langle \{c(t), c^\dagger(0)\} \rangle$，两粒子响应涉及 $\chi \sim \langle [\rho(t), \rho(0)] \rangle$。

7.7 线性响应的边界：何时失效？#

线性响应理论有三个暗含假设，当它们被破坏时，Kubo 公式不再适用：

1. 微扰强度：当 $|F|$ 不满足 $|F| \ll$ 系统的特征能量尺度时，$O(F^2)$ 及以上项不可忽略。例如，强激光场下 TLS 的 Rabi 频率 $\Omega_R = \mu E/\hbar$ 与 $\omega_0$ 可比 → 非线性响应起主导
2. 响应时间尺度：LRT 假设系统始终接近平衡。如果外场驱动系统远离平衡（如脉冲激光在 fs 尺度内将电子加热到 $T_e \gg T_{\mathrm{lattice}}$），线性关系 $\Delta \langle A \rangle \propto F$ 在脉冲期间不成立
3. 量子相干性：在强关联系统（如 Luttinger 液体、分数量子霍尔态）中，低能激发是集体模式而非独立准粒子。需要注意：Lehmann 谱表示本身不是近似——它对任意量子系统都严格成立。真正可能失效的是"把谱函数理解成许多独立准粒子跃迁"的简单图像。在强关联系统中，谱函数仍然存在，但峰可能变宽、分裂，甚至不再对应清楚的单粒子准粒子 → 需要非微扰方法（Bethe ansatz、共形场论、张量网络……）来计算谱函数本身

与你的 Rabi 模型研究的关系：USC 区间 $(g/\omega \gtrsim 0.1)$ 的 Rabi 模型恰好触碰了边界 #2——系统-环境耦合不再"弱"，$\rho_{\mathrm{ss}}$（非平衡稳态）与 $\rho_0$（热平衡态）的差异不能忽略 → 需要 Ban 等 (2017) 的推广 LRT 或 Liouvillian 微扰论。

关键判断准则：LRT 的有效性等价于以下条件：

第一个条件（响应函数的 $L^1$ 可积性）在无耗散系统中可能不成立（如理想晶体在 $\omega = 0$ 时 $\sigma(\omega) \propto \delta(\omega)$ 不绝对可积——但这仅是分布意义下的 Dirac delta，物理上仍可处理）；第二个条件在相变临界点处失效（$\lim_{F \to 0} \Delta M/F = \infty$，磁化率发散）。

8. 涨落耗散定理#

8.1 核心问题#

平衡态下的热涨落（无外场时的自发涨落）与系统对外场的耗散响应（外场做功的耗散率）之间，有什么定量关系？

答案是涨落耗散定理（Fluctuation-Dissipation Theorem, FDT）。

8.2 物理图像#

考虑一个简单的例子：Brown 粒子在液体中的运动。

• 无外力时：粒子在热涨落下做随机运动，其均方位移 $\langle x^2(t) \rangle$ 由扩散系数 $D$ 决定
• 有外力时：粒子以迁移率 $\mu$ 移动，$v = \mu F$

Einstein 关系 $D = \mu k_B T$ 是最简单的 FDT：描述耗散的系数（$\mu$）和描述涨落的系数（$D$）被温度联系在一起。

量子 FDT 的完整形式：

8.3 经典极限与量子极限#

• 经典极限 $\hbar\omega \ll k_B T$（高温或低频）：$\tanh(\beta\hbar\omega/2) \approx \beta\hbar\omega/2$ → $\Im[\chi_{AA}] \approx -\frac{\beta\omega}{2} C_{AA}$，等价地 $C_{AA} \approx -\frac{2k_B T}{\omega} \Im[\chi_{AA}]$ → 恢复 Nyquist 定理和 Einstein 关系（负号来自本文约定，不同教材可能整体差一个符号）
• 量子极限 $T \to 0$（$\beta \to \infty$）：$\tanh(\beta\hbar\omega/2) \to \operatorname{sgn}(\omega)$ → $\Im[\chi_{AA}(\omega)] = -\frac{\operatorname{sgn}(\omega)}{\hbar} C_{AA}(\omega)$ → 零温下仍有量子涨落（零点运动）

8.4 FDT 的推导：从 Kubo 公式出发#

FDT 不是额外假设——它可以直接从 Kubo 公式导出。推导过程揭示了 FDT 的根源在于 $\rho_0$ 是 Boltzmann 分布。

从 Kubo 公式的谱分解出发，计算 $\Im[\chi_{AA}(\omega)]$：

另一方面，计算关联函数 $C_{AA}(\omega) = \int_{-\infty}^{\infty} dt \, e^{i\omega t} \frac{1}{2} \langle \{A(t), A(0)\} \rangle_0$ 的谱表示：

注意 $C_{AA}(\omega)$ 是 $\omega$ 的偶函数（来自反对易子的对称性），而 $\Im[\chi_{AA}(\omega)]$ 是奇函数（来自对易子的反对称性）。

现在比较两式。利用 $\delta$ 函数，将 $\Im[\chi_{AA}]$ 中对 $E_n, E_m$ 的求和拆分：

第二项交换 $m \leftrightarrow n$，利用 $|\langle n|A|m\rangle|^2 = |\langle m|A|n\rangle|^2$（$A$ 自伴），得：

与 $C_{AA}(\omega)$ 比较（注意 $C_{AA}$ 含 $\delta(\hbar\omega + \cdots)$ 项仅在 $\omega < 0$ 时贡献）：

推导中唯一用到的假设：$\rho_0 = e^{-\beta H_0}/Z$（Gibbs 态）。如果初始态不是热平衡态（如 Rabi 模型的 NESS），FDT 需要用更复杂的"非平衡涨落耗散关系"取代（见 §13.2）。

8.5 为什么 FDT 如此重要？#

1. 实验：通过测量平衡涨落（无外场！）就能预测系统对外场的响应——中子散射测 $S(\mathbf{q}, \omega)$ 就是利用这个原理
2. 理论：统一了输运理论（Kubo 公式）和统计力学（Einstein 关系、Nyquist 定理、Onsager 关系）
3. 计算：在很多情况下，计算关联函数比直接计算响应函数更容易（或反过来）

8.6 非微扰推广（选读）#

这一小节是选读。大四阶段只需要知道：FDT 的核心并不只是"一阶微扰"的结果——它更深地来自热平衡态的时间反演结构。下面的非微扰版本只是佐证这一点，不要求掌握推导。

Jacob & Goold (2025) 证明：即使不用 Kubo 的线性近似，散射理论中仍存在一个非微扰版本的 FDT：

其中 $\chi_{\Delta}$ 和 $C_{\Delta}$ 用精确散射振幅（而非 Born 近似）定义。这暗示 FDT 比 Kubo 公式本身更基本——它来自时间反演对称性和平衡态的性质，而不依赖弱耦合。

9. 静态响应、动态响应与零频极限#

9.1 静态响应#

静态响应是 $\omega \to 0$ 的极限：外场随时间变化极慢（绝热）。此时：

静态磁化率 $\chi = \partial M/\partial B$、静态电极化率 $\alpha = \partial P/\partial E$ 都是 $\chi(0)$ 的例子。

9.2 等温响应 vs 绝热响应#

需区分两种静态极限：

• 等温响应（isothermal）：$\omega \to 0$ 但系统始终与热库接触 → 温度恒定
• 绝热响应（adiabatic）：$\omega \to 0$ 且系统孤立 → 熵恒定

两者一般不相等。它们的差别来自系统是否能和热库交换热量——具体差多少依赖系统的热容、热膨胀系数以及被测响应量的性质。对气体声速等简单例子，会出现类似 $C_P/C_V$ 的因子；但在一般线性响应理论中，不能简单认为所有静态响应都只差一个 $C_P/C_V$。

9.3 Kubo-Strěda 公式（只需知道物理意义）#

这一节只需要理解物理结论，不要求掌握严格推导。Kubo-Strěda 公式是 Kubo 公式在零温、静态、Fermi 投影条件下对霍尔电导率的特殊形式。它最重要的结论是：横向电导率不仅和 Fermi 面附近的态有关，还和所有占据态的几何结构（Berry 曲率）有关。

当 $\rho_0 = \chi_{(-\infty, E_F]}(H_0)$（零温 Fermi 投影）时，静态霍尔电导率简化为：

这是双重对易子公式（DCF）——Kubo 公式在 $\omega = 0, T = 0$ 的静态极限。其物理意义是：霍尔电导等于 Fermi 海在实空间中的 Chern 数。

9.4 常见误区#

• "静态响应 = 取 $\omega = 0$ 代入任何公式" → 需要先算动态响应再取极限，否则可能丢失拓扑项（$\delta(\omega)$ 的贡献）
• "Kubo-Strěda = Kubo 的直接化简" → DCF 需要 $X$ 和 $Y$ 的非对角部分（OD），标准 Kubo 中 $X$ 无界带来的发散被抵消——这是非平凡的技术步骤

10. 有限温度下的线性响应#

10.1 温度进入的两种方式#

1. 通过 $\rho_0 = e^{-\beta H_0}/Z$（初始态）：每个本征态按 Boltzmann 权重贡献
2. 通过 $\tanh(\beta\hbar\omega/2)$（FDT）：量子涨落与热涨落的竞争

10.2 零温极限#

$T \to 0$ 时，$\rho_0 \to |\Omega\rangle\langle\Omega|$（基态投影）。Kubo 谱分解中只有 $m = \Omega$ 的项存活（$e^{-\beta E_\Omega}/Z \to 1$，其余 $e^{-\beta E_m}/Z \to 0$）：

第一项：外场将基态 $|\Omega\rangle$ 激发到 $|n\rangle$，然后 $A$ 测量该激发 → 共振吸收（$\omega \approx (E_n - E_\Omega)/\hbar$ 时分母趋近 $i\eta$，响应最强）。第二项：$|n\rangle$ 先被 $V$ 从基态"反激发"（$V$ 作用于 $|\Omega\rangle$ 产生 $|n\rangle$，但这是初态为 $|\Omega\rangle$ 时不可能的逆过程）→ 反共振项，在旋转波近似下通常被丢弃，但对 Berry 曲率和拓扑响应有不可忽略的贡献。

10.3 有限温度的核心效应：热布居与 detailed balance#

有限温度下，$\rho_0 = e^{-\beta H_0}/Z$ 不再是纯态。这意味着相比于零温，有三个本质变化：

1. 激发态也可以作为初态。Kubo 公式中 $m$ 不限于基态——任何 $E_m$ 的态都以权重 $e^{-\beta E_m}/Z$ 参与。这使得：

• 吸收峰（$\omega > 0$）旁出现受激发射峰（$\omega < 0$）——激发态布居可被外场"拉下来"，向外场释放能量
• 两者的净权重由布居差 $(e^{-\beta E_m} - e^{-\beta E_n})$ 控制——这正是 detailed balance 的体现

2. 热展宽。在 $T > 0$ 时，谱函数中的 $\delta$ 峰被替换为热平均——实际上是许多不同初态 $m$ 贡献的叠加，每个初态 $m$ 对应不同的跃迁能量 $E_n - E_m$。对连续谱系统，这导致谱峰的热展宽（thermal broadening）。

3. 静态响应发散被截断。对金属等无能隙系统，零温下某些静态响应（如 Pauli 磁化率 $\chi_{\mathrm{Pauli}} \propto \rho(E_F)$）是有限值，但 $T > 0$ 时 Fermi-Dirac 分布的尾部使得更多态参与响应 → 温度修正 $\sim (k_B T/E_F)^2$。

Detailed balance 条件：从 FDT $\Im[\chi_{AA}(-\omega)] = -\Im[\chi_{AA}(\omega)]$ 可以推出，正频（吸收）和负频（发射）的强度之比为：

这正是 Einstein 1917 年推导的受激发射与吸收系数之比——平衡态下低能级的布居多于高能级，所以吸收强于发射。

10.4 高温极限：恢复经典物理#

$k_B T \gg \hbar\omega$ 时（对室温 $T \approx 300\,\mathrm{K}$，$k_B T \approx 25\,\mathrm{meV}$，对应 $\omega \ll 6\,\mathrm{THz}$——远红外及更低频），$\tanh(\beta\hbar\omega/2) \approx \beta\hbar\omega/2$：

负号来自本文约定，物理上 $\Im[\chi_{AA}] < 0$（正频吸收）→ $C_{AA} > 0$。此时量子 FDT 退化为经典 Nyquist 形式：热涨落 $\propto k_B T$，响应/耗散由 $\Im[\chi]/\omega$ 决定。

经典极限的经验判据：$k_B T \gg \hbar \times$ （系统的特征频率）。对 NMR（$\omega \sim 100\,\mathrm{MHz}$，$\hbar\omega \sim 4 \times 10^{-7}\,\mathrm{eV}$），室温下经典 FDT 极好；对光学吸收（$\omega \sim 10^{15}\,\mathrm{Hz}$，$\hbar\omega \sim 4\,\mathrm{eV}$），室温下量子效应主导。

10.5 Matsubara 形式#

有限温度下，虚时（Matsubara）格林函数比实时格林函数更方便计算：

其中 $\omega_n = 2\pi n/\beta$（玻色型）是 Matsubara 频率。实时响应函数由解析延拓 $i\omega_n \to \omega + i\eta$ 得到。

10.6 Matsubara 计算实例：TLS 的虚时关联函数#

以 §4.6 的两能级系统为例，展示 Matsubara 形式的具体计算。定义虚时 Heisenberg 算符 $A(\tau) = e^{\tau H_0} A e^{-\tau H_0}$（注意没有 $i$——这是 Matsubara 和实时形式的关键区别）。

对于 $H_0 = \frac{\hbar\omega_0}{2}\sigma_z$，$\sigma_x(\tau) = e^{\tau\hbar\omega_0/2 \cdot \sigma_z} \sigma_x e^{-\tau\hbar\omega_0/2 \cdot \sigma_z} = \sigma_+ e^{\tau\hbar\omega_0} + \sigma_- e^{-\tau\hbar\omega_0}$。

Matsubara 格林函数（$A = V = \mu\sigma_x$）：

计算虚时排序期望值（对 $0 < \tau < \beta$）：

积分 $\int_0^{\beta} d\tau e^{i\omega_n\tau \pm \tau\hbar\omega_0} = \frac{e^{\pm\beta\hbar\omega_0} - 1}{i\omega_n \pm \hbar\omega_0}$（利用 $e^{i\omega_n\beta} = 1$，因为 $\omega_n = 2\pi n/\beta$）。得：

解析延拓 $i\omega_n \to \omega + i\eta$：

而实时推迟格林函数 $G^R(\omega) = -\chi_{xx}(\omega)$。分解因式 $(\omega + i\eta)^2 - \omega_0^2 = (\omega - \omega_0 + i\eta)(\omega + \omega_0 + i\eta)$，部分分式后恰好恢复 §4.6 的 $\chi_{xx}(\omega)$。

要点：Matsubara 的优势在于——虚时排序 $T_\tau$ 和 Boltzmann 因子 $e^{-\beta H_0}$ 的结合使得计算完全在有限温度下进行，无需先算零温再"加热"。代价是最后一步的解析延拓 $i\omega_n \to \omega + i\eta$ 在数值上可能是不适定问题（小噪声可能被放大）。对简单模型（如 TLS）解析延拓可精确完成；对多体系统，常需借助 Padé 近似或最大熵方法。

10.7 Matsubara 频率求和：从虚时格林函数到响应函数#

前两节展示了如何在虚时框架内计算 $\mathcal{G}(i\omega_n)$，但最终目标是得到实频率的 $\chi_{AV}(\omega)$。这需要两步：频率求和和解析延拓。

步骤 1：频率求和。Matsubara 格林函数通常以 $i\omega_n$ 的有理函数形式出现（如 $\mathcal{G}(i\omega_n) \propto 1/((i\omega_n)^2 - \omega_0^2)$）。在微扰展开中，更常见的任务是计算形如 $\frac{1}{\beta} \sum_n \mathcal{G}(i\omega_n) e^{i\omega_n 0^+}$ 的 Matsubara 频率和。标准技巧是将求和转化为围道积分：

其中 $n_B(z) = 1/(e^{\beta z} - 1)$ 是 Bose 分布函数，$C$ 是包围虚轴上所有 Matsubara 极点的围道。$n_B(z)$ 的极点恰好在 $z = i\omega_n$（$\omega_n = 2\pi n/\beta$）处，留数为 $1/\beta$——这正是 Matsubara 求和的来源。将围道变形到实轴上下，可把虚时求和转化为实频率积分，过程中自动产生 $\tanh(\beta\hbar\omega/2)$ 或 $\coth(\beta\hbar\omega/2)$ 型的温度因子。

步骤 2：解析延拓。得到 $\mathcal{G}(i\omega_n)$ 后，将 $i\omega_n$ 替换为 $\omega + i\eta$（$\eta \to 0^+$）即得推迟格林函数，进而得 $\chi_{AV}(\omega) = -\mathcal{G}(\omega + i\eta)$（差一个约定因子）。这一步对精确可解模型是直接的代数操作，但对数值数据（如量子 Monte Carlo 产生的 $\mathcal{G}(i\omega_n)$）是著名的数值解析延拓问题——一个不适定反问题。

为什么大费周章走虚时路径？ 对多体系统（如 Hubbard 模型、强关联电子系统），直接在实时框架下做微扰展开需要对连续谱做积分，通常发散或极难收敛。Matsubara 框架把时间轴旋转到虚轴，将被积函数的振荡行为（$e^{i\omega t}$）转换为指数衰减（$e^{-\omega_n \tau}$），使得积分绝对收敛。这就是 §10.5 的 TLS 例子里，虚时积分 $\int_0^{\beta} d\tau e^{i\omega_n\tau \pm \tau\hbar\omega_0}$ 能轻松求出的原因。

10.8 常见误区#

• "有限温度 = 加一个 $\coth(\beta\hbar\omega/2)$ 因子" → 仅对平衡态且线性响应成立。对非平衡稳态（如 Rabi 模型的耗散态），FDT 可能不成立或需推广
• "Matsubara 频率是物理频率" → 它们是辅助量，只在解析延拓后与物理频率联系
• "高温极限下量子效应消失" → 热展宽（$\sim k_B T$）可以大于量子能级间距，但零点涨落（$\sim \hbar\omega/2$）始终存在——高温掩盖了量子效应，但未消除它们。对 $\omega$ 很大的探测（如 X 射线吸收），即使 $T = 300\,\mathrm{K}$，$\hbar\omega \gg k_B T$ 时量子效应仍主导
• "Matsubara 求和只需算少数几项" → 实际上 Matsubara 频率求和通常收敛很慢（$\sim 1/\omega_n$），需要算 $O(10^3)$ 项或用围道积分技巧（§10.7）加速。直接截断会引入系统误差

11. 线性响应与 Berry 曲率 / 几何响应#

11.1 Berry 曲率如何从响应中出现？#

考虑一个系统受绝热慢变参数 $\mathbf{\lambda}(t) = (\lambda_1(t), \ldots, \lambda_d(t))$ 的调制。这些参数可以是磁场方向、外部磁通、晶格动量，或 Rabi 模型中的某个慢变控制参数。每个参数都有一个共轭的广义力：

线性响应关心的是：当参数 $\lambda_l$ 缓慢变化时，另一个广义力 $F_k$ 会不会产生正比于速度 $\dot{\lambda}_l$ 的响应——这就是"几何响应"的物理来源。

根据绝热微扰论，参数变化诱导的广义力期望值为：

其中 $E_0$ 是瞬时基态能量，$\Omega_{kl} = i(\langle \partial_k \psi_0 | \partial_l \psi_0 \rangle - \langle \partial_l \psi_0 | \partial_k \psi_0 \rangle)$ 是 Berry 曲率（Berry curvature）。

第一项是保守力（Hellmann-Feynman），第二项是几何响应力——正比于参数变化速度，垂直（横向）于变化方向。

11.2 从 Kubo 到 Berry 曲率#

将 Kubo 公式的参数取为 $H(t) = H_0 - \mathbf{F}(t) \cdot \mathbf{r}$（外电场），在直流极限 $\omega \to 0$ 和零温极限 $T \to 0$ 下计算霍尔电导率：

这就是 TKNN (Thouless-Kohmoto-Nightingale-den Nijs, 1982) 的核心结论——量子霍尔电导的量子化来自 Brillouin 区上 Berry 曲率的积分 = 陈数（Chern number）。

11.3 Henheik 与 Teufel 的贡献#

Henheik & Teufel (2021) 论证了一个微妙之处：当外场关闭谱 gap 时，Kubo 公式仍可能成立——但需要引入 NEASS（非平衡几乎稳态）的概念。传统的绝热定理假设瞬时本征态跟踪成立（要求 gap 保持开放），但在输运情境中外场会关闭全局 gap。NEASS 方法通过"局部 gap 结构"绕过了这个限制。

物理直觉（图 1 of Henheik & Teufel）：施加电压 $\Delta U$ 在距离 $L$ 上，虽然全局势 $\sim xE$ 关闭了 gap，但电子要从占据带跃迁到非占据带，仍需克服 gap $g$ 或隧穿距离 $g/E$——只要 $g/E$ 很大，原基态仍是"几乎稳态"。

11.4 和我的研究的关系#

• Berry 相论文中的 STA 测量（Zhang 2017, Leek 2007）：绝热捷径（shortcut to adiabaticity, STA）添加反绝热项以加速 Berry 相的积累。线性响应在"STA 轨迹附近"的微扰稳定性分析中，Kubo 公式给出了驱动噪声对 Berry 相位的影响
• 几何量子门（Zhao 2021, Deng 2026）：几何门的保真度受控于几何相位的环境响应 → 线性响应框架可用于计算退相干率（$1/T_2^\ast \propto$ 噪声谱密度 $\times$ Berry 曲率在参数空间中的涨落）

12. 线性响应在 Rabi / Spin-Boson 模型中的用法#

12.1 De Filippis 等 (2023) 的完整推导：耗散 Rabi 模型的 BKT 量子相变#

这篇 PRL 是 LRT 在强关联开放系统中应用的典范——它展示了线性响应如何直接探测量子相变。下面按作者的推导链逐步展开。

12.1.1 模型与关键映射#

哈密顿量 $H = H_{Q\text{-}O} + H_I$ 中，两能级-腔部分为：

环境由 $N$ 个谐振子组成，与腔模的坐标 $x = (a + a^\dagger)/\sqrt{2m\omega_0}$ 线性耦合：

其中 $M_i$ 是第 $i$ 个浴振子的质量，$k_i$ 是将该振子连接到腔模坐标 $x$ 的弹簧劲度系数。$\frac{k_i}{2}(x - x_i)^2$ 即"腔模偏离 $x$ 越远，拉回浴振子 $x_i$ 的力越大"——这是典型的 Caldeira-Leggett 模型 [Rev. Mod. Phys. 59, 1 (1987)]，也是 Spin-Boson 模型的核心。$\{k_i, M_i\}$ 这 $2N$ 个微观参数的物理效应全部封装在谱密度 $J(\omega)$ 中：

浴的微观参数 $\{k_i, M_i\}$ 封装为谱密度 $J(\omega) = \sum_i \frac{k_i \omega_i}{2m\omega_0} \delta(\omega - \omega_i) = \alpha_{\mathrm{cav}}\,\omega\,\Theta(\omega_c - \omega)$（Ohmic，$\alpha_{\mathrm{cav}}$ 无量纲，$\omega_c$ 截断）。

12.1.2 推导 $J_{\mathrm{eff}}(\omega)$：腔+浴对角的四个步骤#

这是全文最关键的技术步骤——把一个二次型玻色子环境精确消去，得到两能级系统"看到的"有效谱密度。

为什么可以单独处理腔+浴而不涉及两能级？ 完整哈密顿量为：

关键观察：$H_I$ 和腔模的 $\omega_0 a^\dagger a$ 合在一起，是纯二次型玻色子哈密顿量——只含 $a, a^\dagger, \{x_i\}, \{p_i\}$，不含任何 Pauli 矩阵。两能级系统通过 $g\sigma_z(a + a^\dagger)$ 只与一个玻色子自由度（腔模坐标 $x \propto a + a^\dagger$）耦合，$\sigma_z$ 本身不是玻色子算符，不能被对角化。

概念上，标准的 Caldeira-Leggett 操作是显式对角化：

1. 对角化玻色子部分——把 $\{a, \{x_i\}\}$ 这 $N+1$ 个耦合谐振子对角化为 $N+1$ 个独立正常模 $\{Q_\lambda\}$
2. 对角化后，$x = \sum_\lambda c_\lambda Q_\lambda$（腔模坐标是正常模的线性组合）
3. 两能级-腔耦合变成 $g\sigma_z(a + a^\dagger) \propto \sigma_z \sum_\lambda g_\lambda Q_\lambda$——原本只与腔模耦合的 $\sigma_z$，现在与所有正常模耦合
4. 正常模集合 $\{Q_\lambda\}$ 构成了有效环境，谱密度 $J_{\mathrm{eff}}(\omega) \equiv \frac{\pi}{2}\sum_\lambda g_\lambda^2 \delta(\omega - \omega_\lambda) \propto \sum_\lambda |c_\lambda|^2\delta$

但在本问题中，De Filippis 等并没有显式做这个对角化——对 $N \sim \infty$ 的连续浴，$N+1$ 维对角化不可行。他们利用了一个关键事实：两能级系统只与一个玻色子组合（$x$）耦合，不需要所有 $c_\lambda$，只需要 $\sum_\lambda |c_\lambda|^2\delta$——而这就是 $x$ 的推迟响应函数 $D^R(\omega)$ 的虚部（§7.2 的 $\Im[\chi] \leftrightarrow S$）。于是对角化被运动方程 + 谱表示替代：

实际执行：只处理玻色子部分 → 运动方程消去浴 → 得 $D^R(\omega)$ → $\Im[D^R]$ 直接给谱权重 → $J_{\mathrm{eff}} \propto g^2 \Im[D^R]$。以下四步即为这个"绕过对角化"的实际推导。

步骤 1：从哈密顿量到运动方程。取玻色子部分 $H_{CB} = \omega_0 a^\dagger a + H_I$，用坐标-动量改写（$x = (a+a^\dagger)/\sqrt{2m\omega_0}$，$p = i\sqrt{m\omega_0/2}(a^\dagger - a)$）：

对 $x$ 和 $p$ 用 Hamilton 方程 $\dot{x} = \partial H/\partial p$，$\dot{p} = -\partial H/\partial x$：

消去 $p$：$m\ddot{x} = \dot{p} = -m\omega_0^2 x - \sum_i k_i(x - x_i)$，即

对浴振子 $i$ 同理：$\dot{x}_i = p_i/M_i$，$\dot{p}_i = -k_i(x_i - x)$，消去 $p_i$ 得：

这两条方程构成一个封闭的线性常微分方程组——腔模和浴振子通过 $k_i(x - x_i)$ 力耦合。

步骤 2：消去浴振子。对第 $i$ 个浴振子做 Fourier 变换（$\partial_t \to -i\omega$）：

代入 $x$ 的运动方程（用 $k_i = M_i\omega_i^2$）：

最后一项化简：

第一个 $\sum_i k_i$ 项恰好与括号内的 $+\sum_i k_i$ 相消。第二个 $\sum_i k_i \frac{\omega^2}{\omega_i^2 - \omega^2}$ 项需要从离散和转换为连续积分。

从离散和到连续积分的转换：谱密度 $J(\omega) = \sum_i \frac{k_i \omega_i}{2m\omega_0} \delta(\omega - \omega_i)$ 的定义正是为了打包微观参数。对任意测试函数：

我们要转换 $\sum_i k_i \frac{\omega^2}{\omega_i^2 - \omega^2}$。先把 $k_i$ 写成 $k_i = 2m\omega_0 \cdot \frac{k_i \omega_i}{2m\omega_0} \cdot \frac{1}{\omega_i}$：

被积函数中 $J(\omega')/\omega'$ 自然出现。代入 $J(\omega') = \alpha_{\mathrm{cav}} \omega'$（Ohmic 谱密度，定义在 $\omega' \le \omega_c$）：

步骤 3：推迟边界条件与自能。积分 $\int_0^{\omega_c} d\omega'/(\omega'^2 - \omega^2)$ 在 $\omega' = \omega$ 处有极点。加上推迟边界条件（$\omega \to \omega + i0^+$——和本笔记 §4.4 的 $+i\eta$ 完全相同的操作），用 Plemelj 公式（§4.6）分离实部和虚部：

实部的 Cauchy 主值积分计算如下。第一步：部分分式分解

第二步：分别做主值积分。对于第一项（奇点在 $\omega' = \omega$）：

（左右对称挖去奇点，奇点两侧的无穷大相互抵消，剩下边界值的对数差）

第二项（奇点在 $\omega' = -\omega$，不在积分区间 $[0, \omega_c]$ 内——是普通积分）：

第三步：两项相减：

取绝对值并调换分子分母（绝对值内取负号不改变值）得常用形式：

负号被吸收进自能的定义 $h(\omega)$ 中（本文和 De Filippis 等的符号约定可能差一个负号，不影响 $J_{\mathrm{eff}}$ 的最终形式——分母中 $h(\omega)$ 会平方）。

虚部的计算更直接：Plemelj 公式 $\Im[(\omega'^2 - (\omega + i0^+)^2)^{-1}] = \frac{\pi}{2\omega}\delta(\omega' - \omega)$，积分后得 $i\frac{\pi}{2\omega}\Theta(\omega_c - \omega)$（$\Theta$ 来自积分上限 $\omega_c$ 的截断）。

因此腔模的自能：

因此腔模坐标的运动方程变为：

其中自能 $\Sigma(\omega) = h(\omega) - i\pi\alpha_{\mathrm{cav}}\omega_0\omega$，实部 $h(\omega) = \alpha_{\mathrm{cav}}\omega_0\omega \ln[(\omega_c + \omega)/(\omega_c - \omega)]$，虚部来自 Plemelj 公式。自能的虚部是耗散的标志——它来自浴的无穷多自由度产生的不可逆弛豫（完全平行于 §3.3 中外场 $F(t)$ 的 $i\eta$ 如何编码耗散）。

步骤 4：从腔的推迟格林函数到 $J_{\mathrm{eff}}$——谱表示的直接应用。

4a. 问题的等价形式。按上节概念步骤，对角化后 $x = \sum_\lambda c_\lambda Q_\lambda$，耦合变为 $\sigma_z \sum_\lambda g_\lambda Q_\lambda$（$g_\lambda = g\sqrt{2m\omega_0}\,c_\lambda$），谱密度 $J_{\mathrm{eff}}(\omega) \equiv \frac{\pi}{2}\sum_\lambda g_\lambda^2 \delta(\omega - \omega_\lambda)$（记为 (SD)）。目标是不做对角化求出 $\sum_\lambda |c_\lambda|^2\delta$。

4b. $D^R(\omega)$ 的谱分解——为何 $\Im[D^R]$ 包含 $\{c_\lambda\}$ 的信息。腔模坐标 $x$ 的推迟响应函数（用 §3-4 的 Kubo 公式框架，$\hbar = 1$）：

对频域做 Fourier-Laplace（§4.2，注意本文 §3-4 的 $\phi_{AV}$ 在此为 $D^R$）：

其谱分解（§4.3 的式 (4.1) 取 $A = V = x$）为：

4c. 零温极限——正常模的权重浮现。取 $T=0$（$\beta\to\infty$），基态 $|0\rangle$ 为真空，激发态 $|\lambda\rangle$ 为单个正常模的单量子态（$H_{\mathrm{normal}}|\lambda\rangle = \omega_\lambda|\lambda\rangle$）。对谐振子系统，只有 $x$ 能连接相差一个量子的态，即 $\langle 0|x|\lambda\rangle \neq 0$ 仅当 $|\lambda\rangle$ 是单量子激发。此时 $E_0 = 0$，$E_\lambda = \omega_\lambda$，权重 $e^{-\beta E_0}/Z \to 1$，$e^{-\beta E_\lambda}/Z \to 0$：

取虚部（§4.6 Plemelj：$\Im[(-\omega_\lambda + \omega + i0^+)^{-1}] = -\pi\delta(\omega - \omega_\lambda)$）：

4d. 矩阵元 $|\langle 0|x|\lambda\rangle|^2$ 与 $\{c_\lambda\}$ 的关系。正常模坐标 $Q_\lambda$ 满足标准谐振子对易关系：$Q_\lambda = (b_\lambda + b_\lambda^\dagger)/\sqrt{2\omega_\lambda}$。单量子态 $|\lambda\rangle = b_\lambda^\dagger|0\rangle$，故 $\langle 0|Q_\lambda|\lambda\rangle = 1/\sqrt{2\omega_\lambda}$。将 $x = \sum_\lambda c_\lambda Q_\lambda$ 代入：

代回式 (★)：

中间的 $1/\omega_\lambda$ 被 $\delta$ 函数替换为 $1/\omega$。因此：

4e. 从 $\{c_\lambda\}$ 到 $J_{\mathrm{eff}}$。将 $g_\lambda = g\sqrt{2m\omega_0}\,c_\lambda$ 代入谱密度定义 (SD)：

4f. 代入 $D^R$ 的显式。步骤 1-3 推导的是无外力的齐次方程 $[-\omega^2 + \omega_0^2 + \Sigma(\omega)]\hat{x}(\omega) = 0$。现在考虑施加一个微小的外力 $F_{\mathrm{ext}}(t)$ 到腔模坐标 $x$ 上——这正是 §1.3 经典阻尼谐振子受迫运动的量子版本。运动方程变为：

Fourier 变换后，加上浴的消去结果（步骤 2-3），得到非齐次方程：

响应函数（Green 函数）的定义：$\hat{x}$ 对外力的线性响应为 $\hat{x}(\omega) \equiv D^R(\omega) \cdot \frac{1}{m}\hat{F}_{\mathrm{ext}}(\omega)$。代入上式：

这就是 Green 函数的"逆算子"定义：运动方程写成 $\hat{L}\hat{x} = \hat{F}/m$，则响应函数 $D^R = \hat{L}^{-1}$。完全平行于 §1.3 中经典阻尼谐振子的 $\chi(\omega) = 1/[m(\omega_0^2 - \omega^2 - i\gamma\omega)]$——这里 $\Sigma(\omega) = h(\omega) - i\pi\alpha_{\mathrm{cav}}\omega_0\omega$ 扮演了经典阻尼 $\gamma$ 的角色（实部 $h(\omega)$ 为色散/Lamb 位移，虚部 $\pi\alpha_{\mathrm{cav}}\omega_0\omega$ 为耗散）。

代入 $\Sigma(\omega)$ 的显式：

取其虚部（分子分母乘共轭，和 §4.6 的 Plemelj 操作相同）得：

（注意分子的 $\omega_0\omega$ 和分母含 $\omega^2$ 项的来源——这是 $x$ 为位置坐标时 $D^R$ 的特殊形式，与 §1.3 阻尼谐振子的响应函数 $\chi(\omega)$ 结构完全平行。）

4g. 最终结果。代入 $-2g^2 m\omega_0\omega \Im[D^R]$，其中 $m\omega_0$ 和 $\Im[D^R]$ 分母的因子组合恰好整理为：

（系数 $2g^2\omega_0^2\alpha_{\mathrm{cav}}$ 来自 $2g^2 m\omega_0\omega \cdot \pi\alpha_{\mathrm{cav}}\omega_0\omega / (m\omega_0\cdots)$——正常模质量 $m$ 在公分母中消去，这是谐振子系统谱分解的必然结果。具体数值因子依赖 $D^R$ 和 $J_{\mathrm{eff}}$ 的归一化约定——De Filippis 等 [2023] 的式 (2) 采用了上述形式。）

线性响应理论的连接点——从 $D^R(t)$ 的定义到 $J_{\mathrm{eff}}$ 的最终表达式，每一步都是本笔记前面章节的直接应用：

12.1.3 $J_{\mathrm{eff}}$ 的物理结构#

在 $\omega \approx \omega_0$ 处：令 $\omega = \omega_0 + \delta$（$|\delta| \ll \omega_0$），$h(\omega_0) \approx 0$，分母化为 $(-2\omega_0\delta)^2 + (\pi\alpha_{\mathrm{cav}}\omega_0^2)^2$——这是一个中心在 $\omega_0$、半宽 $\Gamma = \pi\alpha_{\mathrm{cav}}\omega_0/2$ 的 Lorentzian。腔模被耗散展宽：不再是单一的 $\delta$ 峰，而是有一定宽度的共振。

在 $\omega \ll \omega_0$ 处：$\omega^2 \ll \omega_0^2$，$h(\omega) \approx 0$，分母 $\approx \omega_0^4$，分子 $\propto \omega$。因此：

这是低频 Ohmic 尾巴——有效耦合常数 $\alpha_{\mathrm{eff}}$ 同时依赖于两能级-腔耦合 $g$ 和腔-浴耦合 $\alpha_{\mathrm{cav}}$。即使 $\alpha_{\mathrm{cav}}$ 极小，只要 $g/\omega_0$ 足够大，$\alpha_{\mathrm{eff}}$ 可超过 BKT 临界值 1。

与 LRT 框架的联系：整个推导中，$i0^+$ 的推迟边界条件（§4.4）、Plemelj 公式（§4.6）、谱表示（§7.2）——都是本笔记前面构建的工具在真实研究中的精确应用。$J_{\mathrm{eff}}(\omega)$ 扮演的角色正是 §7.1 定义的谱函数：它决定了 $\Im[\chi(\omega)]$ 在 $\omega$ 处的权重，进而通过 FDT 完全确定了平衡涨落和线性响应。

12.1.4 路径积分与 BKT 相变的根源#

有了 $J_{\mathrm{eff}}$，消去浴自由度后，系统的平衡态由欧氏作用量描述（$\beta = 1/T$）：

核 $K_{\mathrm{eff}}(\tau) = \int_0^{\infty} d\omega J_{\mathrm{eff}}(\omega) \frac{\cosh[\omega(\beta/2 - \tau)]}{\sinh(\beta\omega/2)}$ 将浴的全部信息压缩为两能级自旋在虚时上的非局域相互作用。代入 §12.1.3 的低频 Ohmic 形式 $J_{\mathrm{eff}}(\omega) \simeq (\alpha_{\mathrm{eff}}/2)\omega$：

这是一个 $1/\tau^2$ 的长程（在虚时轴上！）铁磁相互作用——将两能级系统映射为一维经典自旋链，自旋分布在长度为 $\beta$ 的虚时链上。$1/\tau^2$ 的衰减指数恰好是 BKT 相变的边界：衰减比 $1/\tau^2$ 更快则无序相（量子隧穿主导），更慢则有序相（自旋冻结）。这就是BKT 量子相变的数学根源——它产生于 $J_{\mathrm{eff}}$ 的低频 Ohmic 尾巴，而 Ohmic 尾巴本身源于腔模被耗散浴展宽（§12.1.3 的 Lorentzian 峰的低频翼）。

12.1.5 相变判据：手征性方法#

定义自旋-自旋关联函数的零频分量（平方磁化强度）：

$M^2$ 是 BKT 相变的序参量——在 $T = 0$ 且 $g < g_c$ 时为零（无序相，量子隧穿主导），在 $g \ge g_c$ 时跃变为非零（有序相，局域化/冻结）。为精确确定临界耦合，引入标度量 $\Psi(\alpha_{\mathrm{eff}}, \beta) = \alpha_{\mathrm{eff}} M^2$ 和函数：

BKT 理论的普适跳变 $\Psi_c = 1$ 要求：在 $\alpha_{\mathrm{eff}} = \alpha_c$ 处，$G$ 渐近地与 $\beta$ 无关——这唯一确定了 $\alpha_c \approx 1.05$。代入 $\alpha_{\mathrm{eff}} = 4g^2\alpha_{\mathrm{cav}}/\omega_0^2$，得 $g_c/\Delta = (\omega_0/\Delta)\sqrt{\alpha_c/(4\alpha_{\mathrm{cav}})}$。

12.1.6 线性响应：弛豫函数与磁化率的精确关系#

这是论文中与本笔记主题最紧密的部分。沿着 $z$ 轴施加弱磁场 $h$，量子比特沿 $z$ 方向的磁化弛豫函数为：

$\bar{\sigma}_z$ 是微扰后的非平衡期望值。在 Mori-Zwanzig 投影算子形式下，$\Sigma_z(t) = \{[\bar{\sigma}_z(0), \sigma_z(t)]\} / \{[\sigma_z(0), \sigma_z(0)]\}$。其 Laplace 变换 $\Sigma_z(z)$ 与磁化率 $\chi(z) = -i\int_0^{\infty} dt e^{izt} \langle [\sigma_z(t), \sigma_z(0)] \rangle$ 严格相关：

这是固体物理中光学电导率与电流-电流关联函数关系的类比：弛豫函数 = 磁化率的"规整化差分"除以频率。

利用对易关系 $[\sigma_z, H] = -i\Delta\sigma_y$ 和对易子的 Kubo 内积 $(\sigma_y, \sigma_y)$，可推出两条精确恒等式：

恒等式 1（静态磁化率与序参量的关系）：

这意味着：磁化率的低频权重直接联系到相变序参量。$\beta \to \infty$ 时，$M^2\beta$ 在 $g < g_c$ 趋于有限值（无序相），在 $g \ge g_c$ 发散——这正是 BKT 相变在 LRT 中的直接信号。

恒等式 2（弛豫函数的 y-z 耦合）：

由此可定义有效能隙 $\Delta_{\mathrm{eff}}^2 = [(\sigma_y, \sigma_y)/(\sigma_z, \sigma_z)]\Delta^2$——$g = 0$ 时恢复 $\Delta$，随 $g \to g_c$ 趋于零。

12.1.7 实验可观测的 QPT 信号#

作者用 Feynman-Mori 变分法 + MPS 数值模拟计算了 $\Sigma_z(t)$ 和 $\Re[\Sigma_z(\omega)]$，得到三个区域的清晰信号：

$\omega = 0$ 处的 $\delta$ 峰是 QPT 的最锐利判据——它的出现意味着 $\Im[\chi(\omega)] \propto \delta(\omega)$，由 FDT 得 $\chi(0) = \infty$（静态磁化率发散）。在实验中，这表现为零频吸收峰的无限变窄——受限于有限温度和有限测量时间，实际观测到的是一个宽度 $\sim \max(k_B T, 1/t_{\mathrm{meas}})$ 的窄峰。

12.1.8 与本笔记 LRT 框架的对照#

De Filippis 等的推导每一步都有 LRT 对应物：

核心信息：De Filippis 等没有使用开放系统的主方程（Lindblad 形式），而是直接在路径积分框架中精确消去浴自由度，用 LRT 将响应函数与热力学序参量联系起来。这种方法不依赖 Born-Markov 近似，因此在 USC 区间——此时系统-环境的耦合强度与系统本身能量尺度相当——仍然准确。

12.2 De Filippis 方法 vs 标准 LRT 的关键区别#

上节推导展示了开放系统 LRT 的两种进路的对比：

核心教训：在 USC 区间，不能先用 Born-Markov 近似再算 LRT——Ban 等 (2017) 的四部分响应函数告诉我们，$\phi_{BA}^{(c)}$（热库关联修正）在非 Markov 区间不可忽略。De Filippis 等通过在路径积分层面精确消去浴，本质上把 $\phi_{BA}^{(c)}$ 完整保留了下来。

12.3 Spin-Boson 模型#

Spin-Boson 模型 ($H_{SB} = \frac{\Delta}{2}\sigma_x + \frac{\varepsilon}{2}\sigma_z + \sum_k \hbar\omega_k b_k^\dagger b_k + \sigma_z \sum_k g_k(b_k + b_k^\dagger)$) 是耗散两能级系统的"标准模型"。

线性响应在此模型中的应用：

1. 动态磁化率 $\chi_{zz}(\omega)$：探测自旋对外部偏置 $\varepsilon(t)$ 的响应 → 实验上对应介电损耗或量子比特的退相干谱
2. Ohmic 耗散的线性响应：当谱函数 $J(\omega) = \pi \sum_k g_k^2 \delta(\omega - \omega_k) = 2\pi\alpha\omega$（Ohmic）时，$\chi_{zz}(\omega)$ 在低温下的行为由 $\alpha$（无量纲耗散强度）决定：
   • $\alpha < 1/2$：相干振荡（coherent）
   • $\alpha > 1/2$：非相干衰减（incoherent/overdamped）
   • $\alpha = 1$：Kosterlitz-Thouless 量子相变（localization transition）

12.4 Ban 等的开放系统线性响应公式#

Ban 等 (2017) 给出了开放系统中线性响应函数的最一般形式（四个部分）：

各部分含义：

• $\phi_{BA}^{(0)}$：忽略所有关联（系统-热库、热库不同时刻之间）→ 这是 Markov 极限
• $\phi_{BA}^{(c)}$：热库关联修正（热库在不同时刻的态之间的关联）→ 非 Markov 性的直接体现
• $\phi_{BA}^{(i,0)}$：系统-热库初始关联的贡献
• $\phi_{BA}^{(i,c)}$：两种关联的合成效应

关键结论：量子回归定理（quantum regression theorem）仅在 Markov 极限下成立——如果热库有记忆，两时间关联函数 $\langle B(t) A(\tau) \rangle$ 不能仅由单时间传播子计算。这对 Rabi 模型在 USC 区间意味着：当 $g/\omega \gtrsim 0.1$ 时，标准的量子光学主方程（Born-Markov）失效——需要 Ban 等的推广公式或等效的非微扰方法。

12.5 具体应用：NMR 线形与动态磁化率#

核磁共振（NMR）是线性响应理论最经典的应用之一。下面展示如何从 Kubo 公式出发，推导 NMR 线形。

物理设定：核自旋 $I = 1/2$ 在静磁场 $B_0 \hat{z}$ 中（Larmor 频率 $\omega_L = \gamma B_0$），同时受到局部环境（其他核自旋、电子云屏蔽……）产生的随机涨落场 $\delta B(t)$ 的作用。加弱的横向射频场 $B_1(t)\cos(\omega t)\hat{x}$。

哈密顿量：$H_0 = -\hbar\omega_L I_z$，$H_{\mathrm{ext}}(t) = -\hbar\gamma B_1(t) I_x$（射频场），加上与环境的耦合 $H_{SB} = I_z \sum_k g_k(b_k + b_k^\dagger)$（简化 Spin-Boson 型）。

NMR 信号：测量的是横向磁化率 $\chi_{+-}(\omega)$，可观测量为 $I_+ = I_x + iI_y$，外场耦合到 $I_x = (I_+ + I_-)/2$。由 Kubo 公式：

Bloch 唯象结果的微观推导：在 Markov 极限（环境关联时间 $\tau_c$ 远短于自旋弛豫时间 $T_1, T_2$）下，上式给出：

即 Lorentzian 线形，半高宽 $= 1/T_2$——这就是 Bloch 方程的"横向弛豫"的微观起源。$T_2$ 由环境噪声谱密度 $J(\omega)$ 在 $\omega = 0$ 和 $\omega = \omega_L$ 处的值决定：

$T_2^\ast$（纯退相）来自环境在 $\omega \to 0$ 的低频噪声——这正是量子比特退相干的核心机制。

非 Markov 修正：当环境关联时间不可忽略时（如 USC 区间的 Rabi 模型或低温 Spin-Boson），Lorentzian 线形失效。此时需用 Ban 等 (2017) 的四部分响应函数 $\phi_{+-} = \phi_{+-}^{(0)} + \phi_{+-}^{(c)}$（假设无初始关联），其中 $\phi_{+-}^{(c)}$ 给出频率依赖的退相干率 $\Gamma_2(\omega)$：

这里 $\delta\omega(\omega)$ 是频率依赖的 Lamb 位移，$\Gamma_2(\omega)$ 正比于环境的量子噪声谱 $S_{BB}(\omega) = \frac{1}{2\pi} \int dt e^{i\omega t} \langle \{\delta B(t), \delta B(0)\} \rangle$——这就是线性响应理论的终极实验意义：你测量 NMR 线形 → 你得到 $\Im[\chi(\omega)]$ → 通过 FDT 你反推出环境的量子噪声谱 → 你知道了系统与环境的耦合强度和机制。

12.6 从 LRT 视角看你的论文收藏#

你的图书馆中有几篇论文可以从 LRT 视角重新审视：

• De Filippis 等 (2023)：Rabi 模型的耗散临界行为。临界点附近 Liouvillian gap 变小 → 系统弛豫时间变长 → 低频响应通常会增强。是否真的表现为 $\chi(\omega) \sim 1/\omega$ 或某个观测量发散，要看具体模型、极限顺序、噪声谱和稳态定义，不能只凭 LRT 形式直接断言。但 LRT 提供了一个有用的诊断：临界慢化必然伴随低频响应增强——这可以作为非平衡相变的一个判据
• Leek 等 (2007)：Berry 相观测。施加微波脉冲沿参数空间闭合路径驱动量子比特 → Berry 相 = 路径围成的立体角 = Kubo 公式在绝热极限下给出的 $\sigma_{xy}$。你的 Berry 相论文和 LRT 论文本质上在研究同一个物理量——Berry 曲率——只是前者从几何相位的时域积累切入，后者从频率依赖的响应函数切入
• Zhang 等 (2017)：STA 测量 Berry 相。反绝热项的作用是在不破坏Berry 相的前提下加速演化 → 从 LRT 视角：STA 改变了有效 Liouvillian 的谱，使得 $\chi(\omega)$ 在绝热频率 $\omega \sim \eta$ 处仍有良好定义（而不是被 $1/\eta$ 的发散所污染）

13. 输运系数：从 Kubo 公式到电导率#

13.1 Kubo-Greenwood 公式#

将 §4 的一般 Kubo 公式特化到电导率——这是 LRT 最经典也最重要的应用。

设定：电子系统 $H_0 = \sum_i \frac{\mathbf{p}_i^2}{2m} + V(\mathbf{r}_i)$，外加均匀电场 $\mathbf{E}(t)$。在长度规范下 $H_{\mathrm{ext}}(t) = -e \sum_i \mathbf{r}_i \cdot \mathbf{E}(t)$。

电流算符：$\mathbf{J} = \frac{e}{m} \sum_i \mathbf{p}_i$。电导率张量由 $J_\alpha(\omega) = \sum_\beta \sigma_{\alpha\beta}(\omega) E_\beta(\omega)$ 定义。由 Kubo 公式：

其中 $v_\alpha = p_\alpha/m$ 是速度算符，$n$ 是电子密度，第二项是抗磁项（diamagnetic term）——保证 $\sigma(\omega)$ 满足 gauge 不变性和求和规则 $\int_0^{\infty} d\omega \Re[\sigma_{xx}(\omega)] = \pi n e^2/2m$。

13.2 Kubo-Bastin 公式与拓扑#

在无序系统中使用位置算符不方便。改为用速度-速度关联：

这是 Kubo-Bastin 公式——无序系统中计算电导率的标准起点。

当 $\omega \to 0, T \to 0$ 时：

$P = \chi_{(-\infty, E_F]}(H_0)$，$\mathcal{T}$ 是体平均迹。这正是量子（反常）霍尔效应中电导率 = Chern 数 $\times \, e^2/h$ 的微观公式。

13.3 纵向 vs 横向响应#

• 纵向电导率 $\sigma_{xx}(\omega)$：电流 $\parallel$ 电场 → $\Re[\sigma_{xx}]$ = Joule 热耗散
• 横向电导率 $\sigma_{xy}$：电流 $\perp$ 电场 → 可非耗散 → 拓扑电流，时间反演破缺系统中量子化

13.4 常见误区#

• "电导率只依赖 Fermi 面附近" → 对纵向 $\sigma_{xx}$ 成立（Mott 公式），但霍尔电导率 $\sigma_{xy}$ 依赖于所有占据态（Fermi 海的积分）——它是热力学量而非仅输运量
• "速度-速度关联直接给出电导率" → 需要先算 $r$-$r$ 关联再取 $\omega \to 0$，或用 $\sigma_{\alpha\beta} = \lim_{\omega \to 0} \frac{i}{\omega} \chi_{J_\alpha J_\beta}(\omega)$（$1/\omega$ 因子来自 $r = v/(i\omega)$）

14. 初学者技术栈路线图#

步骤 1 (1-2天)
├── 经典阻尼谐振子的受迫响应 → 理解 χ(t), 因果性, K-K关系
├── 量子力学含时微扰论回顾 (Sakurai Ch. 5)
└── 正则系综基础 (Pathria Ch. 2-3)

步骤 2 (2-3天)
├── 从含时微扰论推导 Kubo 公式 (本笔记 §3-4)
├── 计算最简单的例子：两能级系统的动态磁化率
├── 理解 χ(ω) 的实部/虚部各自代表什么
└── 阅读 Jacob & Goold (2025) 附录 A

步骤 3 (3-5天)
├── 涨落耗散定理 (本笔记 §8)
├── 推迟格林函数入门 (Fetter & Walecka Ch. 3-4)
├── 谱函数与动态结构因子 (本笔记 §7)
└── 练习：从 χ(ω) 反推 S(ω)，从 S(ω) 理解中子散射截面

步骤 4 (1-2周)
├── 有限温度 Matsubara 形式 (Mahan Ch. 2)
├── Berry 曲率与量子霍尔效应的拓扑解释
├── Spin-Boson 模型的动态磁化率计算
└── 阅读 Henheik & Teufel (2021) — 理解非平衡几乎稳态

14.1 需要掌握的核心概念清单#

• 含时微扰论 → 一阶修正密度矩阵
• 响应函数 $\phi_{AV}(t)$ 的定义与物理解释
• 频域易感率 $\chi_{AV}(\omega)$ 与 $\phi_{AV}(t)$ 的关系
• Kramers-Kronig 关系的因果性起源
• 关联函数 $C_{AV}(t)$ 与涨落耗散定理
• 谱函数 $S_{AV}(\omega)$ 与可测量（中子散射、吸收谱）
• 推迟、超前、因果、热（Matsubara）四种格林函数的区别
• 零温极限 vs 经典高温极限
• Kubo-Strěda 公式 → Chern 数 → 量子霍尔效应
• 开放系统 LRT 中 Markov 性 vs 非 Markov 性的判据

14.2 最容易踩的坑#

1. 混淆 $\chi(\omega)$ 和 $C(\omega)$：前者是 $\int_0^{\infty}$ 的 Laplace-Fourier，后者是 $\int_{-\infty}^{\infty}$ 的 Fourier → 前者在上半平面解析，后者不一定
2. 忽略 $i\eta$ 收敛因子：$\chi(\omega) = \lim_{\eta \to 0^+} \cdots$ —— $\eta$ 的符号编码了推迟（$+i\eta$）vs 超前（$-i\eta$）边界条件
3. 以为线性响应只能用于平衡态：Ban 等 (2017) 和 Konopik & Lutz (2019) 已将 LRT 推广到非平衡稳态
4. 把 Matsubara 频率当物理频率：Matsubara $\omega_n = 2\pi n/\beta$ 是虚轴上离散点，物理 $\omega$ 是实轴上的连续变量
5. 以为"弱外场"等价于"弱系统-环境耦合"：前者是微扰参数小（$F \ll 1$），后者是耗散耦合小（$g \ll \omega$）——两个独立的"弱"条件

15. 参考文献与推荐阅读顺序#

15.1 入门教材（先读）#

1. L. E. Reichl, A Modern Course in Statistical Physics (4th ed., 2016)

   • 第 13-14 章：线性响应与涨落耗散定理的物理介绍
   • 中文读者友好：推导详细，不假设多体理论前置知识
   • 建议：读完本笔记 §1-8 后，读此书第 13 章巩固

2. J. J. Sakurai, Modern Quantum Mechanics (2nd ed.)

   • 第 5 章：含时微扰论（相互作用表象、Dyson 级数、跃迁几率）
   • Kubo 公式推导的前置技能

15.2 标准教材（同步读）#

3. R. Kubo, M. Toda, N. Hashitsume, Statistical Physics II: Nonequilibrium Statistical Mechanics (Springer, 1985)

   • 第 2 章：Kubo 自己的陈述——从物理动机到公式推导
   • 经典名著，细节丰富

4. G. D. Mahan, Many-Particle Physics (3rd ed., 2000)

   • 第 3 章：线性响应与 Kubo 公式（格林函数进路）
   • 第 7 章：有限温度 Matsubara 形式
   • 适合学完本笔记后进一步深入多体

5. A. Altland & B. Simons, Condensed Matter Field Theory (2nd ed., 2010)

   • 第 6 章：线性响应与涨落耗散定理（路径积分进路）
   • 从场论视角统一处理，适合后续提升

15.3 高级教材（后读）#

6. H. Bruus & K. Flensberg, Many-Body Quantum Theory in Condensed Matter Physics (2004)

   • 第 6 章：Kubo 公式与电导率计算
   • 强调实际计算（无序平均、Feynman 图）

7. G. De Nittis & M. Lein, Linear Response Theory: An Analytic-Algebraic Approach (Springer Briefs, 2017)

   • 本图书馆已有完整翻译笔记 content/papers/de-nittis2016-linear-response/
   • von Neumann 代数 + 非交换 $L^p$ 空间的统一数学框架
   • 对想深入理解 LRT 严格数学基础的读者不可或缺
   • 前提：需要泛函分析（Banach/Hilbert 空间、谱定理、$C^*$-代数基础）

15.4 经典论文#

8. R. Kubo, J. Phys. Soc. Jpn. 12, 570 (1957)

   • 原始 Kubo 公式论文。"不可逆过程的统计力学理论 I"
   • 建议：读完本笔记 §1-4 后去读，体会 Kubo 的原初动机

9. J. Bellissard, A. van Elst, H. Schulz-Baldes, J. Math. Phys. 35, 5373 (1994)

   • 非交换几何进路处理量子霍尔效应，证明了 Kubo 导出的 DCF 的量子化
   • $C^*$-代数框架在凝聚态中的经典应用

10. S. L. Jacob & J. Goold, arXiv:2504.12385 (2025)

    • 本图书馆已有完整翻译笔记 content/papers/jacob2025-collision-kubo/
    • Kubo 公式从散射理论中的自主涌现 + 附录 A 是优秀 LRT 教程
    • 建议：读完本笔记后去读附录 A——两者互补

15.5 与 Rabi / Spin-Boson / Berry 相相关的论文#

11. M. Ban, S. Kitajima, T. Arimitsu, F. Shibata, Phys. Rev. A 95, 022126 (2017)

    • 本图书馆已有完整翻译笔记 content/papers/ban2017-linear-response-open-systems/
    • 开放系统 LRT：初始态任意 + 含初始关联 + 四部分响应函数
    • 对理解 Rabi 模型 USC 区间非 Markov 响应至关重要

12. J. Henheik & S. Teufel, Comm. Math. Phys. 373, 621 (2020) & arXiv:2101.00717 (2021)

    • 本图书馆已有完整翻译笔记
    • NEASS 方法论证 Kubo 公式：处理外场关闭 gap 的情况
    • 对理解 Berry 相在输运实验中的绝热前提有直接价值

13. G. De Filippis et al., Phys. Rev. Lett. 130, 210401 (2023)

    • 本图书馆已有完整翻译笔记 content/papers/dissipation-driven-rabi-qpt/
    • Rabi 模型耗散驱动相变——Ness 处 Liouvillian gap 闭合
    • LRT 在非平衡相变临界点处的失效（发散）是理解 USC 物理的关键

14. P. J. Leek et al., Science 318, 1889 (2007)

    • 本图书馆已有完整翻译笔记
    • 固态量子比特中首次观测 Berry 相
    • Berry 相的 LRT 诠释：$\sigma_{xy}$ 作为 Berry 曲率在参数空间的积分

15.6 后续需要补充阅读的材料#

15. A. J. Leggett et al., Rev. Mod. Phys. 59, 1 (1987)

    • Spin-Boson 模型的综述——耗散两能级系统的"圣经"
    • 需要补充：LRT 在 Ohmic/Sub-Ohmic/Super-Ohmic 耗散下的具体计算

16. H.-P. Breuer & F. Petruccione, The Theory of Open Quantum Systems (Oxford, 2006)

    • 第 3-4 章：量子主方程的微观推导（投影算子方法）
    • 第 9 章：开放系统线性响应
    • 建议在读完 Ban 等 (2017) 笔记后补充阅读

17. D. A. Abanin et al., Rev. Mod. Phys. 91, 021001 (2019)

    • 多体局域化（MBL）综述
    • MBL 系统的 LRT（$\sigma(\omega)$ 在低频的行为）是当前活跃研究前沿

16. 补充：公式索引#

最后更新：2026-07-01。基于本图书馆书籍（De Nittis & Lein, Milonni）、论文笔记（Henheik & Teufel, Ban et al., Jacob & Goold, De Filippis et al., Leek et al.）及标准教材消化重组。