3. 量子模拟基础

在分子能量面与反应路径的量子模拟中，我们最终要把“电子在原子核骨架中的能量”翻译成量子硬件能执行的语言。这一翻译过程依赖五个基本构件：量子比特、量子门、量子线路、测量，以及态制备与哈密顿量模拟的总体范式。理解它们之间的分工，比记住单个门的矩阵更重要。

量子比特：信息的载体

经典比特只能处于 0 或 1，而量子比特（qubit）的状态是二维复向量空间中的单位向量，写作 |ψ⟩ = α|0⟩ + β|1⟩，其中 α、β 为复数且 |α|² + |β|² = 1。|α|² 和 |β|² 分别是测量得到 0 和 1 的概率。对分子模拟而言，一个 qubit 通常不直接对应一个电子，而是对应某个自旋轨道（spin orbital）的占据数：|0⟩ 表示未占据，|1⟩ 表示占据。若有 N 个自旋轨道，就需要 N 个 qubit 来编码一个电子组态。

这种编码方式称为 Jordan-Wigner（JW）变换。它把费米子的产生/湮灭算符映射为 Pauli 算符串，例如 a†_j → (X_j − iY_j)/2 ⊗ Z_0 ⊗ … ⊗ Z_{j−1}。JW 变换的代价是引入 Z 串来维持费米子反对易关系，因此一个分子哈密顿量会展开成大量 Pauli 项。以 H₂ 在 STO-3G 基组下为例，4 个自旋轨道对应 4 个 qubit，其电子哈密顿量经 JW 变换后约含 15 个 Pauli 项（系数随核间距变化），这是后续线路设计的基本输入。

量子门：对量子比特的受控操作

量子门是作用在 qubit 上的酉变换。单比特门包括 X（翻转）、Y、Z（相位翻转）、H（叠加）、S、T（相位门）以及绕任意轴的旋转门 R_x(θ)、R_y(θ)、R_z(θ)。其中旋转门在变分量子算法中尤为关键，因为 θ 是可优化的连续参数。双比特门主要是 CNOT（受控非）和 CZ，它们负责制造纠缠，是表达电子关联的必要资源。

在 JW 编码下，一个费米子激发项 a†_p a_q + a†_q a_p 会映射为若干 Pauli 项，其线路实现通常由 CNOT 链、单比特旋转和逆 CNOT 链组成。例如模拟 exp(−iθ Z_0 Z_1) 时，标准做法是：CNOT(q0, q1) → R_z(2θ) 作用在 q1 → CNOT(q0, q1)。这个“CNOT—旋转—CNOT”结构是哈密顿量模拟中最常见的积木，理解它就能读懂大多数量子化学线路。

量子线路：门的有序排列

量子线路是把量子门按时间顺序排列的图形化表示，从左到右表示时间推进。每条水平线代表一个 qubit，门用方块或符号标出，CNOT 用控制点和目标点加竖线表示。线路中通常包含三个区域：态制备区（把初态 |00…0⟩ 变成目标参考态）、演化/变分区（施加参数化门或 Trotter 步）、测量区（把 qubit 投影到计算基）。

以 4-qubit H₂ 的变分线路为例，一种常见结构是：先对每个 qubit 施加 R_y(θ_i) 制备 Hartree-Fock 参考态附近的叠加态，再用 CNOT 链和 R_z 门引入纠缠，最后测量。线路深度直接决定在当前含噪硬件上的可行性，因此实际设计时要在表达能力和噪声之间权衡。

测量：从量子态到经典数据

测量把量子态坍缩到计算基，得到 0 或 1 的经典结果。要估计某个 Pauli 项 ⟨P⟩ 的期望值，需要多次重复制备同一线路并测量，统计平均值。例如估计 ⟨Z_0 Z_1⟩，可分别测量 q0 和 q1，计算 (n_{00} + n_{11} − n_{01} − n_{10})/N。对于 X 或 Y 项，则需在测量前施加基变换门（如 H 门把 X 基转到 Z 基）。

测量次数（shots）决定统计误差，误差约按 1/√N 下降。若哈密顿量有 M 个 Pauli 项，逐项测量需要 M 组线路；实际中会把可对易的项分组，用同一组测量设置同时估计，以减少总 shots。这是量子化学模拟中一个重要的工程优化点。

态制备与哈密顿量模拟的基本范式

整个量子模拟流程可以概括为：编码 → 态制备 → 演化 → 测量 → 能量估计。编码把分子哈密顿量映射为 qubit 哈密顿量 H = Σ c_k P_k；态制备给出一个参考态 |φ₀⟩（常用 Hartree-Fock 态）；演化通过 e^{−iHt} 或变分线路实现；测量得到各 Pauli 项的期望值；最后按 E = Σ c_k ⟨P_k⟩ 合成能量。

演化有两种主流范式。其一是 Trotter 分解：把 e^{−iHt} 近似为 ∏ e^{−i c_k P_k t/n} 重复 n 次，适合时间演化类算法。其二是变分量子本征求解器（VQE）：用参数化线路 |ψ(θ)⟩ 逼近基态，通过经典优化器最小化 E(θ) = ⟨ψ(θ)|H|ψ(θ)⟩。VQE 对线路深度要求较低，是当前含噪设备上分子能量计算的主要方案。

把上述五个环节串起来，就得到分子能量面量子模拟的最小工作流。实际执行时，先确定基组和活性空间，写出二次量子化哈密顿量，用 JW 或 Bravyi-Kitaev 变换映射到 qubit；再选择参考态和线路模板，设定参数初值；然后通过测量估计各 Pauli 项期望值，合成能量；最后用经典优化器更新参数，直到能量收敛。每一步的精度和资源开销都会影响最终能量面的形状，尤其是键长附近曲率较大的区域，对测量误差更敏感。

需要强调的是，量子比特数并不等于电子数，而是等于自旋轨道数；线路深度受限于当前硬件的相干时间和门保真度；测量次数决定了统计误差下限。理解这三条约束，才能在具体分子体系上做出合理的编码与线路选择。