动力学平均场理论的实时迭代方案：面向近期量子模拟的框架

论文信息

标题: Real-Time Iteration Scheme for Dynamical Mean-Field Theory: A Framework for Near-Term Quantum Simulation

作者: Chakradhar Rangi, Aadi Singh, Ka-Ming Tam

发布日期: 2026-01-27

arXiv ID: 2601.19896v1

PDF 链接: 下载 PDF

3 分钟速览

研究问题：如何设计一种天然适配近期量子硬件的动力学平均场理论（DMFT）求解框架，既能描述金属-绝缘体转变，又不受传统虚时间/频率域求解器在量子计算机上难以实现实时间演化的问题限制。
核心方法：直接在实时域中构建 DMFT 自洽循环，把杂质问题映射到一个包含 5 个浴位点的有限一维链上，通过实时拟合杂化函数、用精确对角化求解并迭代更新，直至自洽收敛。
关键结果：仅使用 5 个浴位点和粗时间网格（100 个点），就成功重现了半满 Hubbard 模型在 Bethe 晶格上的金属-绝缘体转变，谱函数中清晰出现 Hubbard 带和费米能级处谱权重的抑制（见图 4）。
主要局限：在极弱耦合区（U ≲ 1），有限浴点的离散效应会引入人为能隙，无法准确还原无相互作用时的连续态密度；该方法在弱耦合极限下定量结果不可靠。
适合读者：研究强关联电子体系、DMFT 方法以及希望将量子计算用于凝聚态物理模拟的研究人员、研究生。

论文背景和研究动机

动力学平均场理论（DMFT）是处理强关联电子系统中的重要工具，能够自洽地将晶格问题映射为一个杂质的量子杂质模型。传统 DMFT 求解多采用虚时间或频率空间的数值方法，如量子蒙特卡罗、精确对角化等，但这些方法在近期含噪声中等规模量子（NISQ）设备上移植困难。量子硬件天然适合模拟实时间演化，而常规 DMFT 框架却很少直接使用实时格林函数。近年来虽然出现了少量基于变分量子本征求解器（VQE）的 DMFT 演示，但大多受限于 “两站点” 近似，只能抽取准粒子重量，无法呈现如 Hubbard 子带等细致谱学特征。

本文正是针对这一瓶颈，提出了一种全新的实时域迭代方案，直接在实时域中用推迟格林函数进行 DMFT 自洽计算。该方法不仅可以从原理上兼容量子计算机的实时动力学模拟，而且即便在经典计算机上用小浴位点数的精确对角化作概念验证，也能给出正确的金属-绝缘体转变和谱特征，从而突破了传统两站点 DMFT 的限制。

核心方法和技术细节

实时域自洽方程与系统映射

作者以半满 Hubbard 模型在无穷配位 Bethe 晶格为例，从虚时间路径积分出发，推导出 DMFT 自洽条件在实时间中的形式（见公式 12 到 14）。具体的自洽关系为：

\Delta^R(t) = (t^*)^2 G_{loc}^R(t)

其中 $\Delta^R(t)$ 是推迟杂化函数， $G_{loc}^R(t)$ 是局域推迟格林函数。与传统方法不同，整个自洽循环完全在实时域中执行，无需进行时频变换。

杂质问题被映射到一个链长 N_bath=5 的有限非相互作用浴上，加上一个相互作用的杂质位点（见图 1）。链上参数包括跃迁 $\{t_0, t_1, t_2, t_3, t_4\}$ 和位能 $\{\epsilon_1,\dots,\epsilon_5\}$ ，在半满情况下粒子-空穴对称性将独立参数简化为 5 个： $t_0$ 、 $t_1$ 、 $t_2$ 、 $\epsilon_1$ 、 $\epsilon_2$ ，其余由对称性确定。这是一个精心选择的离散化方案，兼顾了物理精度和计算效率。

迭代算法步骤

完整的实时域迭代流程包括：

初始化：利用无能隙情况下的解析解（包含贝塞尔函数 $J_1$ ）构造初始杂化函数 $\Delta^{(0)}(t)$ 。
浴参数拟合：在第 k 次迭代时，通过最小化误差函数 $\chi^2(\mathcal{P}) = \sum_n |\Delta^{(k)}(\bar t_n) - \Delta_{chain}^R(\bar t_n;\mathcal{P})|^2$ ，在当前时间网格上拟合出最优的链参数 $\mathcal{P}$ 。
杂质求解：用精确对角化作为概念演示求解器，计算相互作用的杂质格林函数 $G_{imp}^R(t)$ 。为抑制长时误差并模拟有限寿命，引入阻尼因子 $\eta=0.2$ ，即 $G_{imp}^R(t) \leftarrow G_{imp}^R(t)e^{-\eta t}$ 。
更新杂化函数：根据 DMFT 自洽条件 $\Delta_{new}^R(t) = (t^*)^2 G_{imp}^R(t)$ 得到新目标。
线性混合：为确保稳定性，用 $\Delta^{(k+1)}(t) = \alpha \Delta^{(k)}(t) + (1-\alpha)\Delta_{new}(t)$ ，其中 $\alpha=0.8$ ，返回到第 2 步直至收敛（误差 $\epsilon_{RMS}$ 小于某阈值）。

模拟在时间域 $t \in [0,20)$ 、均匀时间网格 $N_t=100$ 上进行，时间分辨力故意设置得较为粗糙，以检验算法的鲁棒性。

谱函数的提取

实时格林函数经过三次样条插值后，进行数值傅里叶变换得到谱函数 $A(\omega)= -\frac{1}{\pi} \operatorname{Im} G_{loc}^R(\omega)$ 。由于有限时间长度 $t_{max}$ 和阻尼的存在，低频分辨率受限于 $\delta \omega \sim 1/t_{max}$ ，但插值有效抑制了时域稀疏导致的频谱混叠。

创新点和贡献

实时间 DMFT 框架：首次将 DMFT 自洽循环完全建立在推迟格林函数的实时域上，与未来在量子计算机上进行实时演化（如 Trotter 分解或变分量子算法）天然兼容。
有限浴点下的谱特征解析：通过 5 个浴点的离散化，明显超越了文献中普遍采用的两站点近似，成功获得三峰结构（中心准粒子峰与两侧 Hubbard 带）的演化，充分展示了方法的物理保真度。
收敛稳定性：在金属到绝缘体的全参数区间（U=2~8）均稳定收敛，约 20-30 次迭代即可完成，表明该方法即便在强关联 Mott 绝缘体区域也没有遭遇数值失稳（见图 2）。
对噪声和粗分辨力的鲁棒性：使用仅 100 个时间点和有限阻尼，仍能捕捉关键的谱学特征，说明该框架对量子硬件上常见的受限资源有较好的弹性。

实验结果分析

论文在四个代表性相互作用强度下进行了数值验证：U=2（弱关联金属）、U=4（中关联）、U=6（临近相变）、U=8（Mott 绝缘体）。

收敛性：图 2 上排显示，所有 U 下 $\epsilon_{RMS}$ 在约 30 次迭代内下降数个量级，体现了算法的高效收敛。下排的浴参数演化还显示，随着 U 增加，杂质-浴耦合 $t_0$ 单调减小，而位能 $\epsilon_1$ 、 $\epsilon_2$ 的绝对值增，大，直观反映了谱权重从费米能级向高能转移的物理过程。

实时格林函数：图 3 给出了插值后的 $G_{imp}^R(t)$ 。弱耦合时缓慢衰减，强耦合时出现快速振荡，对应于高能 Hubbard 带的形成。对数坐标下的拟合质量证明，即便在中间时间区域，离散浴也较好地再现了真实的衰减行为。

谱函数与自能：图 4 是核心结果。在 U=2 时， $A(\omega)$ 呈半圆形金属态，自能虚部在 $\omega=0$ 处极小，指示良定义准粒子。U=4 时开始出现三峰结构，中心准粒子峰收窄，两侧出现隆起的 Hubbard 带。U=6 时中心峰高度大幅抑制，自能实部在费米能级附近变得格外陡峭，标志着进入临界区域。到 U=8，中心峰完全消失，谱权重全部转移到 $\omega \approx \pm 4$ 的 Hubbard 带，自能虚部在 $\omega=0$ 处出现发散——这是 Mott 绝缘体的标志性特征。这些结果都是在只有 5 个浴位点的条件下得到的，充分证明了方法的有效性。

局限与待解决问题

尽管该实时 DMFT 框架展示了巨大潜力，但存在若干显见的局限：

弱关联极限精度不足：作者明确指出，在 U≲1 的极弱相互作用区域，5 个浴位点产生的离散能级间隔会导致谱函数中出现人为能隙，无法重现连续态密度。因此该方法在近自由电子极限下的定量结果不可靠，需增加浴点数才能改善。
频率分辨力受限于有限演化时间：时间窗口 $t_{max}=20$ 和阻尼 $\eta=0.2$ 决定了频率域的分辨力约在 0.05 的量级，不可能分辨更精细的低能激发，这限制了在近费米面区域的分析精度。
当前仅为经典验证：本文用精确对角化作为杂质求解器，尚未在真实的量子硬件上实现。从经典到量子的迁移仍需解决状态制备、Trotter 误差、量子测量噪声等实际问题。
模型简单：工作仅限于单带 Hubbard 模型和 Bethe 晶格，要推广到多轨道、一般晶格或实际材料，还需解决杂化函数拟合的更复杂情形。
非平衡扩展有待实现：虽然作者提出未来可推广到非平衡 DMFT，但目前只在平衡态零温极限下进行了演示，对于泵浦-探测等动力学问题尚需进一步发展。

总体而言，这一工作为在近期量子设备上实现有物理意义的 DMFT 模拟奠定了坚实基础，但在精度、可扩展性和实验实现上仍有大量挑战待克服。