[论文解读] Multi-determinant generalized Hartree-Fock wave functions in Monte Carlo calculations

Chia‐Chen Chang, Miguel A. Morales|arXiv (Cornell University)|Nov 6, 2017

Physics of Superconductivity and Magnetism被引用 4

一句话总结

该论文提出了一种用于辅助场量子蒙特卡罗（AFQMC）模拟的多 determinant 广义 Hartree-Fock（GHF）试探波函数，显著降低了强关联 Hubbard 模型计算中的系统性误差。通过使用具有变分轨道的非正交 Slater 表 determinant 展开，该方法在各种晶格几何结构下实现了接近精确的基态能量，关键情况下的误差低于 0.01%。

ABSTRACT

The quantum Monte Carlo algorithm is arguably one of the most powerful computational many-body methods, enabling accurate calculation of many properties in interacting quantum systems. In the presence of the so-called sign problem, the algorithm typically relies on trial wave functions to eliminate the exponential decay of signal-to-noise ratio, usually at the expense of a bias in the result. The quality of the trial state therefore is critical for accurate simulations. In this work, benchmark results of the ground state auxiliary-field quantum Monte Carlo method are reported for the Hubbard model on several geometries. We demonstrate that when multi-determinant generalized Hartree-Fock states are used as trial wave functions, the systematic errors can be systematically reduced to a low level and the results compare favorably with the exact diagonalization data.

研究动机与目标

• 通过改进试探波函数以实现对强关联体系的精确模拟，解决量子蒙特卡罗中的符号问题。
• 通过采用多 determinant 广义 Hartree-Fock（GHF）态作为试探波函数，减少辅助场量子蒙特卡罗（AFQMC）中的系统性偏差。
• 在不同晶格几何结构和电子填充度下，将多 determinant GHF 试探波函数的性能与精确对角化数据进行基准测试。
• 证明非正交 Slater 表 determinant 展开中的变分自由度可提升 AFQMC 模拟的精度。

提出的方法

• 该方法采用非正交 Slater 表 determinant 展开形式 $|\Phi\rangle = \sum_{i=1}^{n_d} c_i |\varphi_i\rangle$，其中每个 $|\varphi_i\rangle$ 通过参考 determinant 的幺正轨道旋转生成。
• 试探波函数通过广义 Hartree-Fock（GHF）态构建，包含变分参数 $c_i$ 和轨道旋转矩阵 $Z^i$，允许非正交 determinant。
• 采用约束路径辅助场量子蒙特卡罗（AFQMC）方法计算基态能量，试探波函数通过施加约束路径近似以缓解符号问题。
• 变分参数通过两种方法优化：FED（Fiedler 嵌入对角化）和 ResHF（残差-HF），以提升试探态的质量。
• 该方法应用于不同体系尺寸和电子填充度的二维 Hubbard 模型，包括正方、蜂窝、三角和 kagome 晶格。
• 通过与精确对角化数据对比，量化系统性误差并评估精度。

实验结果

研究问题

• RQ1 多 determinant GHF 试探波函数是否能显著降低 AFQMC 模拟 Hubbard 模型时的系统性误差？
• RQ2 与单 determinant 或正交化试探函数相比，使用非正交 Slater 表 determinant 展开如何提升量子蒙特卡罗计算的精度？
• RQ3 GHF 展开中的变分参数和轨道旋转在多大程度上增强了对强关联电子相关性的描述能力？
• RQ4 FED 和 ResHF 优化方案在提升试探波函数质量方面表现如何比较？
• RQ5 在不同晶格类型和电子填充度下，AFQMC 使用多 determinant GHF 试探函数的定量精度如何？

主要发现

• 在 (5,5) 个电子、U=4 的 4×4 正方晶格上，使用 ResHF 优化的 GHF 试探波函数的 AFQMC 能量为 -19.58085(62)，与精确对角化结果 -19.58094 的相对误差仅为 0.0004%。
• 在 (7,7) 个电子、U=4 的 4×4 三角晶格上，AFQMC 能量为 -25.6598(15)，与精确值的偏差仅为 0.0155%。
• 在 (6,6) 个电子、U=4 的 2×2 kagome 晶格上，AFQMC 结果为 -5.62334(11)，与精确能量相比的相对误差仅为 0.0183%。
• 在 (5,5) 个电子、U=4 的 4×3 蜂窝晶格上，AFQMC 能量为 -11.14386(51)，与精确对角化值相比的相对误差为 0.0003%。
• 该方法在多个基准案例中持续将系统性误差降低至 0.01% 以下，优于单 determinant 或正交化试探函数。
• 所有测试构型中，ResHF 优化的 GHF 试探波函数均实现了最低的相对误差，证明其具有更优的变分灵活性和精度。