4. 时间演化算符

在数字量子模拟中，我们真正关心的往往不是某个时刻的静态量子态，而是量子态如何随时间变化。时间演化算符正是描述这一变化的数学核心。它把“初始态”映射为“经过时间 t 后的态”，并且由于量子力学的基本约束，这个映射必须保持概率守恒，也就是幺正的。本节从薛定谔方程出发，逐步建立时间演化算符的定义、指数形式、幺正性，以及它与哈密顿量的关系。

一、薛定谔方程：时间演化的起点

量子力学中，态矢量 |ψ(t)⟩ 随时间的演化由薛定谔方程描述。在坐标表象下，其一般形式为：

iħ ∂/∂t |ψ(t)⟩ = H |ψ(t)⟩

其中 ħ 是约化普朗克常数，H 是系统的哈密顿量，通常是一个厄米算符。这个方程告诉我们：态随时间的变化率由哈密顿量决定。如果 H 不显含时间，即系统是保守系统，那么方程可以形式上直接积分。我们先把 ħ 吸收到 H 的定义中，或采用自然单位制 ħ = 1，这在量子模拟文献中非常常见。此时方程简化为：

i ∂/∂t |ψ(t)⟩ = H |ψ(t)⟩

这个形式更简洁，也便于后续讨论。需要强调的是，H 的厄米性（H† = H）是保证演化幺正的关键前提。

二、时间演化算符的定义

我们希望找到一个算符 U(t, t₀)，使得任意初始态 |ψ(t₀)⟩ 都能被映射到 |ψ(t)⟩：

|ψ(t)⟩ = U(t, t₀) |ψ(t₀)⟩

这个 U(t, t₀) 就是时间演化算符。它依赖于时间区间，而不只是某个绝对时刻。若取 t₀ = 0，则简写为 U(t)。将定义代入薛定谔方程，可以得到 U 本身满足的微分方程：

i ∂/∂t U(t, t₀) = H U(t, t₀)，  U(t₀, t₀) = I

初始条件 U(t₀, t₀) = I 表示“零时间演化不改变态”。这个方程与薛定谔方程等价，但把问题从态矢量转移到了算符上。对于不显含时间的 H，解可以直接写出。

三、指数形式与与哈密顿量的关系

当 H 不显含时间时，上述算符微分方程的解为：

U(t, t₀) = exp[ -i H (t - t₀) ]

这就是时间演化算符的指数形式。若令 Δt = t - t₀，则 U(Δt) = e^{-iHΔt}。这个表达式是量子模拟中最核心的公式之一：它把“时间演化”转化为“哈密顿量的指数”。在数字量子模拟中，我们通常无法直接实现 e^{-iHΔt}，而是通过 Trotter 分解、量子线路或变分方法去逼近它。

指数形式也直接揭示了 U 与 H 的关系：H 是 U 的生成元。换句话说，无穷小时间演化 U(dt) ≈ I - iH dt 中，H 就是生成演化的“方向”和“速率”。哈密顿量的本征值决定了演化的相位频率，本征态则是不随时间改变（仅积累相位）的特殊态。

四、幺正性：概率守恒的保证

时间演化算符必须是幺正的，即 U†U = UU† = I。这一点可以从 H 的厄米性直接推出：

U† = exp[ +i H† Δt ] = exp[ +i H Δt ]
U†U = exp[ +i H Δt ] exp[ -i H Δt ] = I

幺正性意味着内积保持不变：⟨ψ(t)|φ(t)⟩ = ⟨ψ(0)|φ(0)⟩，因此概率幅的模方总和恒为 1。对于量子模拟而言，这既是物理要求，也是数值稳定性的重要保障。如果我们在经典计算机上模拟时间演化，任何非幺正的近似都可能导致概率不守恒，从而产生非物理结果。

五、一个可计算的教学示例

考虑一个最简单的单量子比特系统，哈密顿量取为 H = (ω/2) Z，其中 Z 是泡利 Z 矩阵，ω 是实数频率。这个 H 是厄米的。时间演化算符为：

U(t) = exp[ -i (ω/2) Z t ]

利用泡利矩阵的性质 e^{-iθZ} = cosθ I - i sinθ Z，令 θ = ωt/2，得到：

U(t) = cos(ωt/2) I - i sin(ωt/2) Z

若初始态为 |0⟩，则 |ψ(t)⟩ = U(t)|0⟩ = [cos(ωt/2) - i sin(ωt/2)] |0⟩ = e^{-iωt/2}|0⟩。此时测量 Z 的期望值始终为 +1，但相位在旋转。若初始态为 |+⟩ = (|0⟩+|1⟩)/√2，则：

|ψ(t)⟩ = (e^{-iωt/2}|0⟩ + e^{+iωt/2}|1⟩)/√2

此时测量 X 的期望值为 ⟨X⟩ = cos(ωt)，测量 Y 的期望值为 ⟨Y⟩ = sin(ωt)。这说明在 H = (ω/2)Z 下，态在 Bloch 球的赤道面上绕 Z 轴旋转，角速度为 ω。这个例子虽然简单，却完整展示了“哈密顿量决定旋转轴和角速度”这一核心图像。

六、时变哈密顿量的情形

如果 H 显含时间，即 H(t)，那么指数形式不再直接成立，因为不同时刻的 H(t) 可能不对易。此时时间演化算符需要写成时间排序的指数：

U(t, t₀) = T exp[ -i ∫_{t₀}^{t} H(τ) dτ ]

其中 T 是时间排序算符。这个表达式在量子控制、绝热演化和含时量子模拟中非常重要。对于数字量子模拟，我们通常把总时间切成很多小步，在每一步内近似认为 H 不变，然后用 Trotter 公式拼接起来。这正是从连续时间演化走向离散量子线路的桥梁。

七、从连续演化到数字量子模拟

在数字量子模拟中，我们通常面对的是由多个项组成的哈密顿量，例如 H = H₁ + H₂ + … + Hₖ。直接实现 e^{-iHΔt} 往往很困难，因为各项可能不对易。Trotter 公式给出了一种近似方案：

e^{-i(H₁+H₂)Δt} ≈ e^{-iH₁Δt} e^{-iH₂Δt} + O(Δt²)

当 Δt 足够小时，误差可以控制。更精确的对称 Trotter 公式为：

e^{-i(H₁+H₂)Δt} ≈ e^{-iH₁Δt/2} e^{-iH₂Δt} e^{-iH₁Δt/2} + O(Δt³)

这些公式把连续的时间演化算符拆解为一系列可实现的量子门操作。理解 U(t) = e^{-iHt} 的指数形式和幺正性，是判断这些近似是否合理、误差如何积累的基础。

八、常见误区与检查要点

• 符号约定：本文采用 e^{-iHt}，部分教材写作 e^{+iHt/ħ}，差异来自方程中 i 的位置和单位制选择。使用前必须统一约定。
• 厄米性检查：若 H 不是厄米算符，U 不再幺正，概率不守恒。在构造有效哈密顿量时，必须验证 H† = H。
• 对易性检查：只有当 [H(t₁), H(t₂)] = 0 时，时变情形的指数才能写成普通积分形式；否则必须使用时间排序。
• Trotter 误差：误差阶数与分解方式有关。一阶 Trotter 误差为 O(Δt²)，二阶对称 Trotter 误差为 O(Δt³)。步长 Δt 越小，精度越高，但线路深度越大。