📚 课程目录第3章哈密顿量基础

3. 哈密顿量基础

在数字量子模拟中，我们真正要“模拟”的对象，几乎总是一个哈密顿量。它决定了系统能量如何随状态分布、时间如何演化、测量会得到什么统计结果。因此，理解哈密顿量不是数学准备，而是整个模拟任务的起点。本节从物理意义出发，逐步过渡到泡利算符、泡利字符串、张量积结构，最后落到 Ising、Heisenberg、XY 三类常见自旋模型的哈密顿量形式。

一、哈密顿量的物理意义

哈密顿量 H 是一个厄米算符，它的本征值对应系统允许的能量，本征态对应能量确定的状态。在量子力学中，系统的时间演化由薛定谔方程给出：

iℏ d|ψ(t)⟩/dt = H |ψ(t)⟩
|ψ(t)⟩ = e^{-iHt/ℏ} |ψ(0)⟩

这意味着，只要给定 H 和初始态，演化就完全确定。数字量子模拟的核心任务，就是把 e^-iHt 这类连续演化拆成量子门可以执行的离散操作。因此，哈密顿量的结构越清晰，越容易设计模拟线路。

从物理上看，H 通常由两部分组成：单体项和相互作用项。单体项描述每个自旋或粒子在外场中的能量，例如 Z 方向磁场；相互作用项描述粒子之间的耦合，例如相邻自旋的 ZZ 或 XX+YY+ZZ 耦合。不同模型的区别，本质上就是这两部分的形式不同。

二、泡利算符与泡利字符串

在自旋系统中，最常用的局部算符是泡利算符：

X = [[0,1],[1,0]]
Y = [[0,-i],[i,0]]
Z = [[1,0],[0,-1]]
I = [[1,0],[0,1]]

它们满足 X²=Y²=Z²=I，且不同算符反对易。泡利算符的物理意义很直接：Z 测量自旋在 z 方向的投影，X 对应 x 方向翻转，Y 对应带相位的翻转。在量子模拟中，我们通常把哈密顿量写成泡利字符串的线性组合：

H = Σ_k c_k P_k
P_k ∈ {I, X, Y, Z}^{⊗n}

其中 P_k 是作用在多个量子比特上的张量积，例如 X⊗Z⊗I 表示第 1 个比特作用 X，第 2 个比特作用 Z，第 3 个比特不作用。这种表示法的好处是：每一项都可以直接映射为量子线路中的旋转门，系数 c_k 决定旋转角度。

三、张量积结构与局部性

张量积结构决定了哈密顿量的“局部性”。如果一项只作用在少数相邻比特上，它就是局部的；如果作用范围很广，模拟代价会显著增加。以一维链为例，常见的局部项包括：

Z_i：第 i 个比特的纵向场；
X_i：第 i 个比特的横向场；
Z_iZ_i+1：相邻比特的 ZZ 耦合；
X_iX_i+1 + Y_iY_i+1：相邻比特的 XY 型交换耦合。

在数字量子模拟中，我们通常利用 Trotter 分解把 e^-iHt 近似为多个局部项的指数乘积。局部项越多、耦合范围越大，Trotter 步数和线路深度就越高。因此，张量积结构不仅是数学形式，也直接决定模拟资源。

四、三类常见自旋模型的哈密顿量形式

Ising 模型。一维横场 Ising 模型的哈密顿量为：

H_Ising = -J Σ_{i=1}^{n-1} Z_i Z_{i+1} - h Σ_{i=1}^{n} X_i

其中 J 是相邻自旋的耦合强度，h 是横向磁场强度。当 J>0 时，ZZ 项倾向于让相邻自旋同向排列；当 h 较大时，X 项倾向于让自旋沿 x 方向翻转。两者竞争会在临界点附近产生量子相变。在数字模拟中，通常把 e^-iHt 拆成 ZZ 项和 X 项交替执行，例如：

U(t) ≈ [e^{-i(-J)ZZΔt} e^{-i(-h)XΔt}]^m,  Δt = t/m

这里 m 是 Trotter 步数，Δt 是每步时间。步数越多，误差越小，但线路越长。

Heisenberg 模型。各向同性 Heisenberg 模型写作：

H_Heis = J Σ_{i=1}^{n-1} (X_i X_{i+1} + Y_i Y_{i+1} + Z_i Z_{i+1})

它描述相邻自旋之间的完整交换耦合。与 Ising 不同，Heisenberg 模型同时包含三个方向的耦合，因此 Trotter 分解通常需要把 XX、YY、ZZ 三组项分别处理。若 J>0，系统倾向于反铁磁排列；若 J<0，则倾向于铁磁排列。由于三项耦合强度相同，系统具有 SU(2) 对称性，这会影响可观测量的守恒性质。

XY 模型。XY 模型只保留 XY 平面内的耦合：

H_XY = J Σ_{i=1}^{n-1} (X_i X_{i+1} + Y_i Y_{i+1})

它没有 ZZ 项，因此对称性从 SU(2) 降为 U(1)。在量子模拟中，XY 模型常用于研究自旋输运、量子相变和拓扑激发。其 Trotter 分解只需处理 XX 和 YY 两组项，线路结构比 Heisenberg 更简单。

五、从哈密顿量到模拟线路的对应关系

把哈密顿量写成泡利字符串后，每一项 c_kP_k 都可以通过基变换和旋转门实现。例如，e^-iθZ_iZ_i+1 可以用 CNOT 加 R_z(2θ) 再加 CNOT 实现；e^-iθX_i 直接是 R_x(2θ)。因此，哈密顿量的项数、作用范围和系数，直接决定线路中的门类型、门数量和深度。

在动手写模拟线路前，建议先把哈密顿量整理成“泡利字符串 + 系数”的列表。例如 Ising 模型可写成：
[("ZZ", 1, 2, -J), ("ZZ", 2, 3, -J), ..., ("X", 1, -h), ("X", 2, -h), ...]
这样每一项都对应一个可执行的旋转操作，便于后续 Trotter 分解和线路生成。

注意符号约定：本文采用 H = -J Σ ZZ - h Σ X 的 Ising 形式，其中 J>0 表示铁磁耦合。若你参考的文献使用 H = J Σ ZZ + h Σ X，则 J 的符号含义相反。在模拟前务必统一约定，否则相变位置和基态会判断错误。

总结来看，哈密顿量是数字量子模拟的“输入规格书”。泡利字符串给出每一项的算符形式，张量积结构给出作用范围，系数给出强度，而 Ising、Heisenberg、XY 则是三类最常用的自旋模型模板。掌握这三类形式，就能把大多数自旋链模拟任务转化为可执行的量子线路。