Efficient time-evolution of matrix product states using average Hamiltonians

技术摘要：利用平均哈密顿量实现矩阵乘积态的高效时间演化

问题陈述
模拟由随时间变化的哈密顿量驱动的量子多体系统（QMBS）动力学仍然是一个重大挑战，因为其希尔伯特空间呈指数级增长。虽然张量网络（特别是矩阵乘积态，MPS）为一维系统提供了高效的框架，但标准的演化算法（如时间演化块解离法 TEBD 和 MPO $W_{II}$）主要针对时不变哈密顿量进行了优化。当应用于随时间变化的情景时，标准实现通常依赖于“黎曼步进器”（Riemann stepper），该步进器将时间间隔开始时的瞬时哈密顿量（$H(0)$）代入演化算符中。这种方法会降低高阶积分器的精度，使演化算符的收敛速率从随时间变化情况下的三阶降至二阶，从而限制了给定步长下的模拟精度。

方法论
作者提出了一种对现有 MPS 时间步进器的改进方案，旨在不改变张量网络库基本架构的前提下，恢复对于随时间变化的哈密顿量的高阶精度。其核心方法论涉及使用通过数值积分理论精心构建的“平均哈密顿量”，来替换单步传播器中使用的标准瞬时哈密顿量。

具体而言，作者使用基于辛普森法则（Simpson's rule）计算的哈密顿量取代了标准的黎曼近似：
$$H_{\text{Simp}}(\Delta t) \equiv \frac{1}{6} \left( H(0) + 4H(\Delta t/2) + H(\Delta t) \right)$$
随后将此平均哈密顿量输入到标准的二阶 TEBD 或 MPO $W_{II}$ 算法中。理论分析表明，虽然黎曼步进器（$H(0)$）产生的演化算符误差界为 $O(\Delta t^2)$，但辛普森步进器保留了底层积分器固有的三阶局部精度（$O(\Delta t^3)$）。该证明依赖于使用 Hölder-$\infty$ 范数来界定精确时间演化算符与近似算符之间的差异，表明辛普森近似相比黎曼近似显著降低了积分误差项。

主要贡献与结果
所提方法通过模拟受微波脉冲驱动的钻石中氮-空位（NV）色心的自旋-1 链进行了验证。该系统使用包含漂移项、相互作用项和控制项的哈密顿量进行建模，其参数源自实验数据。

1. 精度提升： 模拟结果对比了“辛普森步进器”与标准“黎曼步进器”在不同系统规模（3 到 8 个 NV 色心）和积分步数下的表现。

   • 对于小规模系统（例如 3 个 NV 色心），在适中步长下，与黎曼步进器相比，辛普森步进器将平均误差降低了约 1000 倍。
   • 对于更大规模的系统（高达 8 个 NV 色心），该方法始终产生更低的传播误差，其平均误差的改进幅度在 50 到 4000 倍之间，具体取决于步数和系统规模。
   • 误差比率（$E_{\text{Riemann}} / E_{\text{Simpson}}$}$）在所有测试配置中均保持在高位，证实了该方法的鲁棒性。

2. 计算成本： 精度的提高伴随着运行时间的轻微增加。辛普森步进器在每个时间步需要三次哈密顿量评估（分别在 $t=0$, $t=\Delta t/2$, 以及 $t=\Delta t$），而黎曼步进器仅需一次。作者指出，这种运行时间的增加非常小，且完全归因于额外的函数评估，使得这种精度与成本的权衡非常有利。

3. 可扩展性： 该方法在系统规模增加时表现出稳定的性能提升，在保证键维数（bond dimension）足够大以使截断误差可以忽略不计的前提下，能够维持其相对于标准方法的优势。

意义与局限性
本文声称，该方法为现有的 MPS 解算器提供了一种简单且高效的增强手段。通过利用经典的数值积分思想（辛普森法则），该方法恢复了那些在模拟随时间变化的动力学时原本会丢失的高阶收敛特性。这对于量子控制、光谱学以及 Floer 工程等随时间变化的驱动场普遍存在的应用领域尤为重要。

作者明确指出了几点局限性和实际考虑因素：

• 光滑性要求： 该方法依赖于哈密顿量 $H(t)$ 在每个区间内存在且有界的至少前四阶导数。
• 哈密顿量构建成本： 如果随时间变化的哈密顿量构建过程计算成本很高，那么每个步长需要三次评估的要求可能会显著增加计算成本。
• 纠缠极限： 虽然该方法改善了积分误差，但它并不能消除大规模多体系统中由纠缠增长带来的限制，这些限制是由 MPS 表示中的截断误差所决定的。

作者建议，辛森平均的概念可以被改编，并可能与更高阶的求积规则（如高斯-勒让德求积）或多阶段组合方法相结合，以构建更高阶的有效哈密顿量，供未来研究使用。