[论文解读] Solutions of the Two-Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

J. P. F. LeBlanc, Andrey E. Antipov+23 more|arXiv (Cornell University)|Jan 1, 2015

Physics of Superconductivity and Magnetism参考文献 137被引用 581

一句话总结

本文通过广泛运用多种数值方法（包括量子蒙特卡洛、图解蒙特卡洛、DMRG 和动力学平均场方法）对二维 Hubbard 模型进行了全面的基准测试研究。在弱耦合、半满及远离半满的参数区域，研究建立了可靠的结果；然而，在接近半满及中等相互作用强度区域，由于量子相的竞争，仍存在持续的不确定性。

ABSTRACT

Numerical results for ground-state and excited-state properties (energies, double occupancies, and Matsubara-axis self-energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary-field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed-node approximation, unrestricted coupled cluster theory, and multireference projected Hartree-Fock methods. Comparison of results obtained by different methods allows for the identification of uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic-limit values is emphasized. Cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.

研究动机与目标

• 评估不同数值方法在求解二维 Hubbard 模型时的准确性和一致性。
• 识别不同方法在误差范围内收敛的参数区域，为未来理论与算法的发展建立可靠的基准。
• 量化由有限体系尺寸效应、图解截断和统计采样引起的数值模拟中的不确定性。
• 突出参数空间中——尤其是接近半满和中等 U 值——系统性误差与相之间竞争导致方法间分歧的区域。

提出的方法

• 采用辅助场量子蒙特卡洛（AFQMC）方法，获得热力学极限下的数值精确结果。
• 应用原始与粗线图解蒙特卡洛（DiagMC）方法，系统性地包含费曼图并施加受控近似。
• 采用双 fermion（DF）方法，研究超越动力学平均场理论的长波长关联。
• 应用密度矩阵嵌入理论（DMET）与密度矩阵变分组态重整化组（DMRG），在有限体系中结合自能修正，研究强关联效应。
• 采用固定节点扩散蒙特卡洛（FN-DMC）与非限制性耦合簇（UCCSD）方法，对比变分与微扰方法。
• 结合多参考投影哈特ree-福克（MRPHF）与簇动力学平均场理论（DCA），探测对称性破缺与非均匀相。

实验结果

研究问题

• RQ1 在热力学极限下，不同数值方法在二维 Hubbard 模型的基态能量与双占据数上收敛到何种程度？
• RQ2 各种数值方法中的主要系统性误差来源是什么？这些误差如何随相互作用强度与电子填充率变化？
• RQ3 在哪些参数区域（U, n）内，多种方法在误差范围内达成一致，从而可作为可靠的基准？
• RQ4 为何在接近半满与中等 U 值区域仍存在显著差异？这是由于数值伪影还是物理相之间的竞争？
• RQ5 数值精确方法（如 AFQMC、DiagMC）之间的一致性是否可用于验证 DMRG 或 UCCSD 等近似方法？

主要发现

• 在所有相互作用强度下，AFQMC、DiagMC 及其他数值精确方法在半满时表现出极佳的一致性，确立了莫特绝缘相的可靠基准。
• 在弱耦合区域及远离半满的载流子浓度下，多种方法在基态能量与双占据数上收敛一致，不确定性极小。
• 在接近半满及中等相互作用强度区域（U ≈ 8–12），方法间出现显著差异，表明存在较大不确定性，且可能存在相之间的竞争。
• 在弱耦合与半满条件下，AFQMC 与 DiagMC 在 Matsubara 轴上的自能表现出一致行为，验证了这些方法在动力学性质研究中的适用性。
• 本研究发现，DMRG 与 UCCSD 的系统性误差在中等 U 区域更为显著，尤其是在接近半满时，可能无法捕捉竞争序。
• 作者得出结论：中等参数区域中剩余的不确定性更可能源于物理原因（如接近量子相变点），而非纯粹的数值伪影。