Fast convergence of Majorana Propagation for weakly interacting fermions

技术摘要：弱相互作用费米子马约拉纳传播（Majorana Propagation）的快速收敛性

问题陈述
模拟哈密顿量演化下的局部可观测量的时间动力学是一个 BQP 完全问题，这意味着除了在受限设定下，通常不期望存在高效的经典算法。本研究关注的是由作用于 $N$ 个马约拉纳模 $\gamma_1, \dots, \gamma_N$ 的四次哈密顿量 $H$ 描述的相互作用费米子系统。目标是近似可观测量 $A$ 的海森堡演化，即 $A(t) = e^{iHt}Ae^{-iHt}$。虽然精确演化在 $t>0$ 时会立即产生具有最大度 $N$ 的可观测量（其中“度”指非零系数的马约拉纳长字符串的长度），但本文研究了在存在低度近似的情况下如何高效找到该近似。

方法论：马约拉纳传播 (MP)
本文分析了结合了基于度数截断（degree-based truncation）与 Trotter 化方法的马约拉纳传播 (MP) 算法。该算法流程如下：

1. 哈密顿量分解： 将四次哈密顿量 $H$ 分解为 $G \le 4\Delta$ 个组 $H = H_1 + \dots + H_G$，其中 $\Delta$ 是稀疏参数（每个马约拉纳模最多出现在 $\Delta$ 项中）。在每个组 $H_g$ 内，马约拉纳字符串具有不相交的支持集，因此它们两两对易。
2. Trotter 化： 将时间演化离散化为步长为 $\delta t$ 的步长。演化算符通过一阶 Trotter 乘积进行近似：$\tau^{\text{trot}}_{\delta t}(A) = e^{iH_G \delta t} \dots e^{iH_1 \delta t} A e^{-iH_1 \delta t} \dots e^{-iH_G \delta t}$。
3. 传播与截断： 可观测量通过每个 Trotter 步进行传播。至关重要的是，在每个步骤之后，应用截断映射 $\text{Trunc}_\ell$，该映射丢弃可观测量展开式中所有度数大于阈值 $\ell$ 的马约拉纳字符串。
4. 迭代： 重复此过程 $t/\delta t$ 次以近似 $A(t)$。

核心技术贡献
本文建立了支持 MP 算法的两个主要理论支柱：

1. 依赖于度的范数不等式： 作者推导出了新的关于对易子和 Trotter 误差的界限，这些界限取决于可观测量度数而非系统规模 $N$。

   • 对易子界限： 对于 $\Delta$-稀疏的四次哈密顿量和度数为 $d$ 的可观测量 $A$，$\|[H, A]\|_F \le 2\Delta\sqrt{d(d+2)}\|A\|_F$。这与传统的随 $N$ 缩放的界限形成对比。
   • Trotter 误差界限： 对于低度可观测量，证明了精确演化与 Trotter 近似之间的误差与 $N$ 无关，而是随度数 $d$、稀疏度 $\Delta$ 和组数 $G$ 缩放。具体而言，$\|\tau_t(A) - \tau^{\text{trot}}_t(A)\|_F \le C \cdot (d\Delta G)^2 t^2$。

2. 收敛性分析： 主定理（定理 1.1/3.3）证明，MP 算法的误差受限于 Trotter 离散化误差与最佳可能的度数为 $\ell$ 的真实时间演化可观测量近似值之间的误差之和。
   $$\|A^{\ell, \delta t}_{\text{MP}}(t) - A(t)\|_F = O\left( \ell^2 t \delta t + (t/\delta t)\eta^\star_F(H, A, t, \ell) \right)$$
   其中 $\eta^\star_F$ 量化了与最佳度数为 $\ell$ 的近似项之间的距离。这一结果提供了第一个可证明的保证，即只要存在低度近似，MP 就能高效地找到它。

弱相互作用哈密顿量的结果
本文将这些结果应用于形式为 $H = H_0 + uV$ 的哈密顿量，其中 $H_0$ 是二次项，$V$ 是四次项，$u \le 1$ 是小的相互作用强度。

• 低度近似的存在性： 利用微扰理论，作者证明了对于时间 $t < t_{\max}(u) \propto \log(e/u)/\Delta$，时间演化的可观测量 $A(t)$ 可以由一个低度多项式近似，且误差随度数 $\ell$ 指数级衰减。
• 高效模拟： 通过结合这些近似的存在性与 MP 的收敛保证，作者证明（推论 1.3）MP 可以模拟任何 $\Delta$-稀疏弱相互作用哈密顿量在时间 $t_{\max}(u)$ 内的动力学，其运行时间为 $N^{O(\log((1+t)/\varepsilon))}$，以实现归一化 Frobenius 范数下的误差 $\varepsilon$。
• 极限行为： 当相互作用强度 $u \to 0$ 时，$t_{\max}(u) \to +\infty$，恢复了已知事实，即二次（自由）费米子系统可以用常数度数精确模拟。

数值验证
作者使用一维和二维晶格上的 Fermi-Hubbard 模型验证了其理论发现。数值实验表明：

• 随着截断度 $\ell$ 的增加，MP 输出向理想 Trotter 演化的指数级收敛。
• 特征时间 $t_{\max}(u)$ 的行为，即当相互作用强度 $U$ 增大时，固定 $\ell$ 的模拟开始偏离精确解。
• 在各种相互作用强度下，模拟二维晶格中空穴密度动力学的高精度。

意义与主张
本文声称提供了关于马约拉纳传播在哈密顿量演化中高效性的第一个可证明的保证。它使“基于截断高阶项的算法对于弱相互作用系统在长时间跨度内是准确的”这一直觉形式化。这项工作将经典模拟技术的适用范围扩展到了超越短时动力学或特定几何局部设置的领域，其依据是相互作用图的稀疏性和相互作用的微弱性。作者指出，虽然这些结果是针对费米子系统推导的，但可以很容易地推广到自旋系统（泡利传播）。其意义在于识别出了一个特定的机制（弱相互作用、稀疏四次哈密顿量），在该机制下经典模拟仍然是高效的，这为量子优势研究提供了基准。