Efficient Classical Simulation of Weakly Interacting Fermion Dynamics

技术摘要：弱相互作用费米子动力学的有效经典模拟

问题陈述
本文解决了模拟弱相互作用费米子系统实时动力学的计算挑战。具体而言，目标是计算在给定时间 $t$ 下局部可观测量 $A$ 的期望值，即 $\text{Tr}(\rho_0 e^{iHt} A e^{-iHt})$，其中 $\rho_0$ 是初始态（例如，费米子高斯态），且哈密顿量 $H = H_0 + \lambda V$ 由占主导地位的二次项（自由项）$H_0$ 和弱四次相互作用 $\lambda V$ 组成。虽然一般的相互作用费米子系统被认为在经典计算机上是不可解的，但作者研究了可以进行高效模拟的结构化机制，重点关注在 $D$ 维晶格上具有几何局域性或具有有界度局域性的系统。

方法论
作者提出了一种基于海森堡绘景中可观测量在相互作用强度 $\lambda$ 下的微扰展开的经典算法。其核心方法包括：

1. 相互作用绘景公式化： 通过定义 $\tilde{A}(t) = e^{-iH_0 t} e^{iHt} A e^{-iHt} e^{iH_0 t}$，将演化的可观测量从自由动力学中分离出来。这使得自由演化 $H_0$ 可以被精确处理（对于二次哈密顿量，这在经典上是高效的），同时将相互作用 $\lambda V$ 作为微扰处理。
2. 戴森级数展开： 将可观测量 $\tilde{A}(t)$ 展开为戴森级数（式 3）。第 $k$ 阶项涉及相互作用 $V$ 与在自由哈密顿量 $H_0$ 下演化的可观测量之间的嵌套对易子。
3. 收敛性的范数分析： 作者使用两种不同的范数来分析该级数的收敛性，以建立不同的运行时间保证：
   • 归一化 Frobenius 范数： 用于证明对于有界度局域哈密顿量的收敛性。该分析依赖于自由演化的酉性，这防止了 Frobenius 范数因 $H_0$ 而增长，从而仅由相互作用结构来限制增长。
   • Majorana 1-范数： 定义为算符在 Majorana 单项式基底下的系数的 $\ell_1$-范数。该范数上界化了算符范数，并且对于设计随机采样算法至关重要。作者通过分析在几何局域性下（利用 Lieb-Robinson 界）的算符扩散来限制该范数的增长。
4. 随机采样算法： 为了避免计算完整的（呈指数级增长的）多项式展开，该算法对 Dyson 级数中的单个 Majorana 单项式进行采样。
   • 算法从时间有序单纯形中均匀地采样时间变量 $(t_1, \dots, t_k)$。
   • 通过递归地采样算符的分量，从嵌套对易子结构中采样单个 Majorana 单项式。
   • 每个样本被分配一个带符号的权重。期望值通过对这些加权样本求平均值来估计。
   • 通过在收敛性分析中推导出的 Majorana 1-范数界限，严格控制采样方差。

主要贡献与结果

1. 有界度局域性的扩展时间机制：
   对于具有有界度局域性（即每个模式参与常数个项）的哈密顿量，作者证明了当 $t = O(1/\lambda)$ 时，Dyson 级数在归一化 Frobenius 范数下呈指数收敛。这指数级地扩展了 Facelli 等人 [17] 建立的已知机制 $t = O(\log(1/\lambda))$。

   • 结果： 对于在归一化 Frobenius 范数下逼近可观测量，具有运行时复杂度为 $\tilde{O}(N^{O(\log(1/\epsilon))})$ 的拟多项式时间算法。

2. 几何局域性的多项式时间算法：
   对于在 $D$ 维晶格上具有几何局域性的哈密顿量，作者建立了在 $\lambda |t| (1+|t|)^{2D} = O(1)$ 条件下，在 Majorana 1-范数下的指数收敛性。

   • 结果： 一种随机采样算法，能在 $\tilde{O}(\epsilon^{-2} M^3 e^{2d_0} N^3)$ 的时间内计算出具有加性误差 $\epsilon$ 的期望值，其中 $M$ 是可观测量中的 Majorana 单项式数量，$d_0$ 是可观测量的阶数，$\tilde{O}(\cdot)$ 隐藏了关于 $M, 1/\epsilon,$ 和 $1/\delta$ 的对数因子。只要相互作用足够弱，使得采样方差保持独立于系统规模，这便代表了相对于现有技术中拟多项式运行时间的超多项式改进。

3. 安德森局域化增强：
   作者将结果扩展到自由部分 $H_0$ 表现出安德森局域化（例如，具有足够无序度的安德森模型）的系统。在这种机制下，算符扩散受到抑制，导致 Majorana 1-范数的增长界限更紧（随时间呈多项式增长变为多对数增长）。

   • 结果： 收敛机制扩展到 $\lambda |t| \cdot \text{polylog}(1+|t|) = O(1)$。该算法对于 $N$ 仍保持多项式复杂度，有效地模拟了无序系统中比洁净系统更长时间的动力学。

意义与主张
本文声称，通过利用弱相互作用、局域性和局域化特性，识别出了广泛的经典可处理实时间费米子动力学的机制。其主要意义在于：

• 严格的方差控制： 不同于许多遭受动力学符号问题（即方差随时间和系统规模指数增长）的量子蒙特卡洛（QMC）方法，该算法提供了严格证明，表明在指定的参数机制内，采样方差与系统规模 $N$ 保持独立且有界。
• 算法改进： 该工作通过将可达时间尺度（从对数级扩展到与 $1/\lambda$ 成反比）以及将计算复杂度（从拟多项式降低到多项式）从 Majorana 传播算法 [17] 进行了改进。
• 理论框架： 通过结合连续时间 QMC、图解 QMC 以及在海森堡绘景下算符增长分析内的 Majorana 传播思想，作者提供了一个针对弱相互作用系统的可证明高效框架。

作者指出，其结果并不排除在更强相互作用或更长时间下存在量子优势的可能性，这表明高效模拟的时间 $t$ 从根本上受限于相互作用强度 $\lambda$ 和精度 $\epsilon$。他们还建议未来的研究方向，例如将其方法与“蠕虫”（inchworm）方法或量子嵌入技术相结合，但并未声称目前已实现这些技术。