# Gibbs 态制备详解：量子热力学、虚时间演化与采样

*前置阅读：[Trotterization 与哈密顿量模拟](../ch05-hamiltonian-qsp/trotterization-tutorial.md)、[量子奇异值变换（QSVT）详解](../ch05-hamiltonian-qsp/qsvt-tutorial.md)；与 [广义 QSP 与 QETU](../ch05-hamiltonian-qsp/gqsp-qetu.md) 的滤波路线直接衔接。*

**课程目标：**

1. 理解 Gibbs 态 $\rho_\beta \propto e^{-\beta H}$ 为什么是"被低估的基础件"：热力学模拟、SDP 求解器、量子 Metropolis 的共同入口。
2. 掌握三条主路线的原理与代价：**虚时间滤波**（非酉 → 后选择/LCU/QSVT）、**热场双态/ purification**、**量子 Metropolis 采样**。
3. 理解 $\Omega(\beta)$ 查询下界的物理来源（制备热平衡至少要等"系统弛豫那么久"），以及 $\beta$ 与精度 $\varepsilon$ 的依赖为什么不同路线差别巨大。
4. 会按"谱隙已知否、需要态还是需要采样、容错时代"选择路线。

:::{admonition} 本课知识点
:class: tip

1. **[Gibbs 态与制备难点](#gibbs-state-motivation)**——能写出 Gibbs 态 $\rho_\beta = e^{-\beta H}/Z(\beta)$ 的定义，列出它作为公共入口的至少三类应用，并解释难点在于 $e^{-\beta H}$ 是非酉算符。
2. **[虚时间滤波](#imaginary-time-filtering)**——能计算虚时间投影后激发权重的衰减因子 $e^{-2\beta\Delta}$，并解释纯滤波为何只给基态、给不出有限 $\beta$ 的 Gibbs 态。
3. **[QSVT 实现与复杂度](#qsvt-gibbs-complexity)**——能写出切比雪夫截断次数 $d = O(\beta\|H\| + \log(1/\varepsilon))$ 的来源，并比较后选择与 LCU 两种组装方式的代价与失败模式。
4. **[热场双态 purification](#thermofield-double)**——能写出 $|\mathrm{TFD}_\beta\rangle$ 的显式形式并通过对环境求部分迹验证它给出 $\rho_\beta$，同时说明双倍寄存器与幅度结构的代价。
5. **[量子 Metropolis 采样](#quantum-metropolis)**——能解释提议-接受-拒绝随机游走为何以 $\rho_\beta$ 为不动分布、为何不需要配分函数，以及拒绝步骤为何不能直接测量能量。
6. **[三条路线比较与选型](#three-routes-comparison)**——能按"要态还是要采样、精度类型、资源时代"比较三条路线对 $\beta$ 与 $\varepsilon$ 的依赖，并针对给定任务选出合适路线。
7. **[查询下界与精度分层](#beta-lower-bound)**——能解释 $\Omega(\beta)$ 查询下界的 Lieb–Robinson 物理来源，并区分"输出态距离"与"输出严格样本"两种精度依赖。
8. **[应用接口与验证](#gibbs-applications)**——能用振幅估计说明 $\log Z$ 的 $O(1/\varepsilon)$ 估计，并对制备出的 $\rho_\beta$ 写出能量涨落 $\langle H^2\rangle - \langle H\rangle^2 = \partial_\beta^2 \ln Z$ 的自洽检查。

:::

(gibbs-state-motivation)=
## 1. 问题：为什么制备一个"热平衡态"这么重要

温度 $\beta = 1/T$ 下的热平衡态（Gibbs 态）是

$$
\rho_\beta \;=\; \frac{e^{-\beta H}}{Z(\beta)}, \qquad Z(\beta) = \mathrm{tr}\, e^{-\beta H}.
$$

它在至少四类问题里出现：

- **材料与化学**：有限温性质（自由能、热容、输运系数）都从 $\rho_\beta$ 的期望值来；经典方法（量子蒙特卡洛）在符号问题面前经常失效，量子计算天然免疫。
- **热力学量与配分函数**：$Z(\beta)$ 给出自由能 $F = -\beta^{-1}\ln Z$；有了 Gibbs 态制备 + 量子幅度/相位估计（第 3 章振幅估计），可以把自由能估到 $\varepsilon$ 精度而只花 $O(1/\varepsilon)$（对比经典蒙特卡洛的 $1/\varepsilon^2$）。
- **半正定规划 (SDP)**：量子机器学习中的 SDP 求解器以 Gibbs 态为中间解的原像（Gibbs 态是相对熵投影的解），制备 Gibbs 态的复杂度直接决定求解器复杂度。
- **量子 Metropolis 采样**：一般系综采样的量子加速也需要热态作为不变分布。

**难点一眼可见**：$e^{-\beta H}$ 是**非酉**的衰减算符——量子线路只能做酉演化，"按能量压低权重"这件事必须绕道。所有 Gibbs 制备算法，本质都是在回答"非酉的 $e^{-\beta H}$ 怎么被酉线路逼出来"。

## 2. 路线一：虚时间滤波——把 $e^{-\beta H}$ 当多项式做

(imaginary-time-filtering)=
### 2.1 从虚时间演化到滤波器

**与虚时间演化的关系**：实时间 $e^{-iHt}$ 在 $t = -i\beta$ 处解析延拓就是 $e^{-\beta H}$。所以"制备 Gibbs 态"的第一直觉是**虚时间演化** (imaginary-time evolution, ITE)：对任意初态 $\rho_0$，

$$
\frac{e^{-\beta H}\,\rho_0\, e^{-\beta H}}{\mathrm{tr}(e^{-2\beta H}\rho_0)} \;\xrightarrow{\ \beta \text{ 增大}\ }\; |\mathrm{GS}\rangle\langle\mathrm{GS}|
$$

低能态按 $e^{-2\beta\Delta}$ 被指数压制——ITE 是**滤波器**，与 [QETU 滤波](../ch05-hamiltonian-qsp/gqsp-qetu.md)是同一件事的两种语言。但要注意：**纯滤波只给基态，不给 Gibbs 态**；有限 $\beta$ 的 Gibbs 态是滤波后的态**再除以配分函数权重的混合**，需要更精细的构造——这正是"用多项式实现 $e^{-\beta x}/Z$ 型函数"的信号处理问题。

(qsvt-gibbs-complexity)=
### 2.2 多项式实现：切比雪夫截断与两条组装路线

**QSVT 式实现**：把 $f(x) = e^{-\beta x}$（在 $H$ 的谱区间 $[-\|H\|, \|H\|]$ 上重标到 $[-1,1]$）做切比雪夫截断（[QSP 详解](../ch05-hamiltonian-qsp/qsp-tutorial.md)"逼近精度"节的 Bernstein 椭圆分析直接适用：截断次数 $d = O(\beta\|H\| + \log(1/\varepsilon))$），随后有两条组装方式：

1. **后选择**：信号处理线路以振幅 $\sim e^{-\beta E}/\sqrt{Z'}$ 保留能量本征态分支，对辅助比特测量成功即得（近似）Gibbs 态——成功率 $1/\mathrm{poly}$，用振幅放大（第 3 章）拉回；
2. **LCU/线性组合**：把 $e^{-\beta H}$ 与"剩余部分"组成一个更大的酉（块编码思想），确定性输出混合态——无失败分支，常数更大。

**关键复杂度事实**：查询复杂度对 $\beta$ 是**线性** $\Theta(\beta)$——这是物理必然（下界），见第 4 节；对精度，截断只贡献 $\log(1/\varepsilon)$，但"混合态输出的精度"（迹距离）通常再引入 $1/\varepsilon$ 或 $\log(1/\varepsilon)$ 依赖，取决于是否允许后选择与重复。

### 2.3 NISQ 变分捷径

**变分捷径 (variational QITE / QMETTS)**：NISQ 时代不追求严格的 QSVT，而用参数化线路做 McLachlan 方程的虚时间流、或"METTS：最大纠缠态 + 虚时间投影交替"的经典-量子混合采样。精度无严格保证，但小系统上好用，是近中期实验的主力。

## 3. 路线二与三：purification 与 Metropolis

(thermofield-double)=
### 3.1 路线二：热场双态 (TFD)

**热场双态 (thermofield double, TFD / purification)**。Gibbs 态 = 一个**纯态的半边约化**：

$$
|\mathrm{TFD}_\beta\rangle \;=\; \frac{1}{\sqrt{Z}}\sum_E e^{-\beta E/2}\, |E\rangle_S |E\rangle_E, \qquad \mathrm{tr}_E\,|\mathrm{TFD}\rangle\langle\mathrm{TFD}| = \rho_\beta .
$$

这把"制备混合态"转化为"制备纯态"——而纯态制备是谱滤波的老本行（QETU/GQSP 或 QPE + 幅度放大，第 5、6 章）。代价：(1) 需要系统两倍大的寄存器（$E$ 侧是"热浴影子"）；(2) 逐能级的幅度 $e^{-\beta E/2}$ 是无奇偶性的光滑目标——正是 [GQSP](../ch05-hamiltonian-qsp/gqsp-qetu.md) 拆掉奇偶约束后大显身手的场景。TFD 路线在现代资源估计里常给出最优常数。

(quantum-metropolis)=
### 3.2 路线三：量子 Metropolis 采样

**量子 Metropolis (quantum Metropolis sampling)**。Temme–Osborne–Vollbrecht–Lidar (2009) 把经典 Metropolis–Hastings 算法搬进量子：在能量本征值间做提议-接受-拒绝的随机游走，且用"相位回踢式"的量子化拒绝步骤避免测量坍缩破坏记忆。不动分布恰为 $\rho_\beta$，因此**不需要知道配分函数**，天然适合"只要采样、不要波函数"的应用（关联函数、热观测量的系综平均）。Chowdhury–Somma (2017) 的容错版本给出多项式复杂度；缺点是混合时间 (mixing time) 依赖问题的物理（无一般保证），且线路结构比滤波路线复杂。

(three-routes-comparison)=
### 3.3 三条路线对比与选型

**三条路线画像**：

| 路线 | 输出 | 对 $\beta$ | 对 $\varepsilon$ | 需要 | 适合 |
|---|---|---|---|---|---|
| 虚时间滤波 / QSVT | Gibbs 态（迹距离 $\varepsilon$） | $O(\beta)$ | $\mathrm{poly}\log$ 到 $1/\varepsilon$（视变体） | 块编码或演化黑盒、能隙或谱信息 | 已知要"态"本身 |
| TFD / purification | 纯态双倍寄存器 | $O(\beta)$ | $\log(1/\varepsilon)$ 友好 | 谱滤波 (QETU/GQSP) | 追最优常数、容错时代 |
| 量子 Metropolis | 系综采样 | $\Omega(\beta)$ + 混合时间 | $\log(1/\varepsilon)$/样本 | 受控演化 | 只要观测值、不需知道 $Z$ |

(beta-lower-bound)=
## 4. 下界与极限：$\beta$ 是硬通货

**查询下界 $\Omega(\beta)$**：无论何种算法，制备 $\rho_\beta$（或从它采样）必须调用黑盒至少 $\Omega(\beta)$ 次。物理来源很直观：热平衡是信息从系统扩散到"热浴"的过程，而量子传播（Lieb–Robinson 型速度极限）以 $O(t)$ 的光锥展开——要建立温度 $\beta$ 对应的关联结构，至少需要"演化 $\beta$ 那么久"。滤波路线恰好做到 $\Theta(\beta)$：**理论已经顶格**。

**精度依赖的分层**：若允许输出"任意 $\varepsilon$-接近的态"，滤波型算法对 $\varepsilon$ 只需 $\mathrm{poly}\log(1/\varepsilon)$（截断多项式足够光滑）；若要求输出严格 Gibbs 态的**样本**（如 Metropolis 的马尔可夫输出），每样本精度一般带 $1/\varepsilon$ 或混合时间因子。资源估计时务必分清自己买的是哪一种"精度"。

**低温极限**：$\beta \to \infty$ 时 Gibbs 态退化为基态，问题变成第 6 章的基态制备（[QETU 滤波](../ch05-hamiltonian-qsp/gqsp-qetu.md)路线），复杂度由谱隙 $\Delta$ 与初态重叠决定而非 $\beta$——两章在低温处无缝衔接；反过来，基态制备的"模糊二分"也可以看作 Gibbs 制备在 $\beta\Delta \to \infty$ 的极限情形。

(gibbs-applications)=
## 5. 应用接口

- **配分函数与自由能**：$\rho_\beta$ 制备 + 振幅估计给出 $\log Z$ 的 $O(1/\varepsilon)$ 估计；相变点定位（自由能交叉）由此而来。
- **量子 SDP 求解器**：以 Gibbs 态为对偶变量的原像，制备复杂度 $\times$ 迭代次数 = 求解器复杂度；Gibbs 的 $\Theta(\beta)$ 下界直接翻译为 SDP 的规模依赖。
- **有限温关联与输运**：Metropolis/TFD 输出的系综样本 + 影子读出（[经典影子应用](classical-shadow-applications.md)）= 完整的有限温物性流水线。
- **验证**：制备出的 $\rho_\beta$ 应通过影子保真度见证（对角熵分布、能量涨落 $\langle H^2\rangle - \langle H\rangle^2 = \partial_\beta^2 \ln Z$ 的自洽检查）。

## 本课总结

- Gibbs 态是热力学模拟、配分函数、SDP、Metropolis 的公共入口；核心难点是 $e^{-\beta H}$ 非酉。
- 三条主路线：虚时间滤波（QSVT/GQSP 实现，$\Theta(\beta)$ 顶格）、TFD purification（纯态化 + 谱滤波，常数友好）、量子 Metropolis（免配分函数，付混合时间）。
- $\Omega(\beta)$ 查询下界由量子传播速度限制背书；精度依赖分层（态距离 $\mathrm{poly}\log$ vs 采样 $1/\varepsilon$）。
- 低温极限连接基态制备，读出端连接经典影子——Gibbs 制备是第 5、6 章技术在中端温度的汇合点。

## 练习题

**练习 1【Gibbs 态与制备难点】**（→ [第 1 节](#gibbs-state-motivation)）

1. 基础：写出温度 $\beta$ 下 Gibbs 态 $\rho_\beta$ 与配分函数 $Z(\beta)$ 的定义，列出正文给出的至少三类应用，并写出 $Z(\beta)$ 与自由能 $F$ 的关系。
2. 进阶：设 $H = \mathrm{diag}(0, \Delta)$，计算 $Z(\beta)$ 与 $\rho_\beta$ 的两个本征权重，并说明 $\beta\Delta \to \infty$ 时 $\rho_\beta$ 趋于什么。

> 提示：低温极限下全部权重应集中到最低能级上，与第 4 节的"低温极限"相衔接。

**练习 2【虚时间滤波】**（→ [2.1 节](#imaginary-time-filtering)）

1. 基础：写出实时间演化 $e^{-iHt}$ 经解析延拓得到 $e^{-\beta H}$ 的对应关系；设谱隙为 $\Delta$、初态与基态重叠为 $\eta$，计算投影后激发权重相对基态权重的衰减因子。
2. 进阶：证明 ITE 的滤波极限：对 $\rho_0$ 与基态重叠 $\eta>0$ 的初态，虚时间投影后的激发权重以 $e^{-2\beta\Delta}$ 衰减；由此估算"用 ITE 把重叠从 $\eta$ 提到 $1-\varepsilon$"需要的 $\beta$（设谱隙 $\Delta$）。

> 提示：归一化前基态权重与激发权重之比为 $\eta : (1-\eta)e^{-2\beta\Delta}$，令归一化后的基态权重超过 $1-\varepsilon$ 解出 $\beta$。

**练习 3【QSVT 实现与复杂度】**（→ [2.2 节](#qsvt-gibbs-complexity)）

1. 基础：写出 $f(x) = e^{-\beta x}$ 重标到 $[-1,1]$ 后的切比雪夫截断次数 $d = O(\beta\|H\| + \log(1/\varepsilon))$，并各用一句话比较"后选择"与"LCU"两种组装方式的成功模式与代价。
2. 进阶：用 Bernstein 椭圆分析（[QSP 详解](../ch05-hamiltonian-qsp/qsp-tutorial.md)）证明 $f(x)=e^{-\beta x}$ 的切比雪夫截断次数 $d = O(\beta\|H\| + \log(1/\varepsilon))$：写出 $M = \max_{E_\rho}|f|$ 与 $\rho$ 的选取。

> 提示：$e^{-\beta x}$ 是整函数，可取固定大小的 Bernstein 椭圆（如 $\rho = e$），并用 $\max_{E_\rho}|\mathrm{Re}\,z| = \frac{\rho+\rho^{-1}}{2}$ 估计 $M$。

**练习 4【热场双态 purification】**（→ [3.1 节](#thermofield-double)）

1. 基础：写出 $|\mathrm{TFD}_\beta\rangle$ 的一般形式，解释它如何把"制备混合态"转化为"制备纯态"，并指出这条路线的两个代价。
2. 进阶：写出二维系统 $H = \mathrm{diag}(0, \Delta)$ 的 $|\mathrm{TFD}_\beta\rangle$ 显式形式；求 $\mathrm{tr}_E$ 后验证为 Gibbs 态；计算制备它需要的幅度比 $\frac{\cos\text{-权重}}{\sin\text{-权重}}$ 并与 $\beta\Delta$ 联系。

> 提示：把 $|\mathrm{TFD}_\beta\rangle$ 参数化为 $\cos\theta\,|00\rangle + \sin\theta\,|11\rangle$，则 $\tan\theta = e^{-\beta\Delta/2}$。

**练习 5【量子 Metropolis 采样】**（→ [3.2 节](#quantum-metropolis)）

1. 基础：说明量子 Metropolis 在能量本征值之间如何做"提议-接受-拒绝"的随机游走，写出它的不动分布，并解释为什么不需要知道配分函数。
2. 进阶：思考：为什么量子 Metropolis 的"拒绝"步骤不能直接测量能量再回退？查 Temme 等人的相位记录解决方案并用自己的话复述。

> 提示：直接测量会让波函数坍缩，破坏提议分布的记忆。

**练习 6【三条路线比较与选型】**（→ [3.3 节](#three-routes-comparison)）

1. 基础：按三条路线画像表回答：哪条路线输出纯态双倍寄存器？哪条天然不需要知道 $Z$？哪条对 $\beta$ 做到 $\Theta(\beta)$ 顶格？并各举一类适合的应用。
2. 进阶：能隙已知时，比较"QETU 模糊二分到基态"（大 $\beta$ 极限）与"QSVT 直接做 $e^{-\beta x}/Z$"两条路线的调用次数公式，找出交叉点（$\beta\Delta \approx ?$）并解释物理含义。

> 提示：把两条路线的调用次数都写成 $\beta\Delta$ 与 $\log(1/\varepsilon)$ 的函数，再联立求交点。

**练习 7【查询下界与精度分层】**（→ [第 4 节](#beta-lower-bound)）

1. 基础：复述 $\Omega(\beta)$ 查询下界的物理来源（Lieb–Robinson 型速度极限），并解释滤波路线的 $\Theta(\beta)$ 为什么说明"理论已经顶格"。
2. 进阶：说明 $\beta \to \infty$ 时 Gibbs 制备如何退化为基态制备、复杂度改由哪些量决定；并解释"输出 $\varepsilon$-接近的态"与"输出严格 Gibbs 态的样本"两种精度目标在 $\varepsilon$ 依赖上的差别。

> 提示：对照第 4 节"低温极限"与"精度依赖的分层"两段。

**练习 8【应用接口与验证】**（→ [第 5 节](#gibbs-applications)）

1. 基础：写出"制备 $\rho_\beta$ + 振幅估计"把自由能估到 $\varepsilon$ 精度的代价量级，并与经典蒙特卡洛对比。
2. 进阶：证明 $\langle H^2\rangle - \langle H\rangle^2 = \partial_\beta^2 \ln Z$，并说明如何用它对制备出的 $\rho_\beta$ 做能量涨落的自洽检查。

> 提示：先证 $\partial_\beta \ln Z = -\langle H\rangle$，再对 $\langle H\rangle = \mathrm{tr}(H e^{-\beta H})/Z(\beta)$ 关于 $\beta$ 求导。