第 22 课 · 量子 ODE/PDE 求解器——从哈密顿量模拟到 Schrödingerization · 模块四:应用与前沿专题

情况二 a:时变线性 ODE 与 Magnus 展开

  • 命题框:基本解 $\Phi(t)$:$d\Phi/dt = A(t)\cdot\Phi$、$\Phi(0) = I$;拟设 $\Phi = e^{\Omega(t)}$——$\exp(\int_0^t A)$ 一般不是解(不同时刻的 $A$ 不可交换)
  • 推导步骤(逐行):
    1. 指数求导恒等式:$d/dt\, e^{\Omega} = \int_0^1 e^{s\Omega}\cdot\Omega'\cdot e^{(1-s)\Omega}\, ds$(由级数定义逐项求导重排)
    2. 两边右乘 $e^{-\Omega}$:$A(t) = \sum_n \mathrm{ad}_\Omega^n(\Omega')/(n+1)!$,其中 $\mathrm{ad}_\Omega(X) := [\Omega, X]$
    3. 按 $t$ 的幂次迭代求解:零阶 $\Omega' = A$ → $\Omega_1 = \int_0^t A(\tau)\, d\tau$(量级 $O(t)$)
    4. 交叉项 $(1/2)[\Omega_1, \Omega_1'] = O(t^2)$ 必须被抵消 → $\Omega_2 = (1/2)\int_0^t ds\int_0^s [A(s), A(\tau)]\, d\tau$
  • 结果框:$u(t) = \exp(\int_0^t A + (1/2)\int\int[A(s),A(\tau)] + O(t^3))\cdot u_0$;$[A(s),A(\tau)] = 0$ 时对易子项消失、退化回常系数;截断到 $O(t^p)$ 项给出 $p$ 阶方法(各指数因子再由 Trotter–Suzuki 实现)
类型definition-derivation