利用等距张量超压缩的辅助场量子蒙特卡洛高效方法

《Journal of Chemical Theory and Computation》：Efficient Auxiliary-Field Quantum Monte Carlo Using Isometric Tensor Hypercontraction

　　

辅助场量子蒙特卡洛（AFQMC）已成为处理强关联电子体系的有力框架，在计算成本与精度之间取得了良好平衡。本文中，研究人员提出一种新颖的AFQMC方法，利用等距张量超压缩（ITHC）技术，通过引入额外虚构费米模式，对分子电子哈密顿量的双体库仑相互作用进行对角化。传统低秩近似如张量超压缩（THC）分解此前虽用于AFQMC以紧凑表示双电子相互作用从而减少内存占用，而本研究方法不仅降低内存使用，相较于标准AFQMC方法，在传播与局域能量评估上均展现出更低的理论复杂度和更优的实际性能。研究人员通过计算线性H10链与苯分子的基态能量验证了该方法的有效性。结果表明，扩展基组AFQMC以与耦合簇（CC）或密度矩阵重正化群（DMRG）等高阶波函数方法相当精度恢复多体关联，同时显著改善了标度行为。

研究背景与问题提出

求解电子薛定谔方程是现代计算化学与计算材料物理的核心。基于电子密度的密度泛函理论（DFT）虽具O(N³)–O(N⁴)优势标度，却在强电子关联体系下失效；而组态相互作用（CI）与耦合簇（CC）等波函数方法虽精确，但标度高达O(N⁴)–O(N¹⁰)甚至指数级。辅助场量子蒙特卡洛（AFQMC）以随机采样确定基态，标度接近DFT且免自相互作用误差，但双电子排斥积分（ERI）张量的存储构成主要内存瓶颈，既有Cholesky或张量超压缩（THC）表示仅缓解存储，未能加速walker传播。为此，研究人员开展将等距张量超压缩（ITHC）引入AFQMC的研究，通过在扩展空间对角化双体项，获得Hubbard型相互作用以同时降低复杂度与内存。

主要关键技术方法

研究人员采用基于插值可分离密度拟合（ISDF）生成通用THC再经三步变换得ITHC因子；将分子哈密顿量双体项分解为V_pqrs=Σ_αβu_pαu_qαW_αβu_rβu_sβ，引入N_aux–N个虚构费米模b_m构造扩展空间湮灭算符c?_α=Σ_pu_pαa?_p+Σ_mv_mαb?_m，使双体项化为W?=?Σ_α≠βW_αβn?_αn?_β；虚时间演化用二阶Trotter-Suzuki分解结合量子Zeno动力学投影回物理空间；Hubbard-Stratonovich（HS）变换将e^–ΔτW?转为单体外场积分；采用phaseless AFQMC（ph-AFQMC）以力偏置（force bias）与重要性采样克服费米符号问题；基于ipie包实现扩展基传播子，walker以Slater行列式表示，通过等距矩阵u作基变换。

研究结果

1. 2.
   Theory（理论）
   2.1 Notation（符号）：界定N、N_e、N_aux、N_w及ITHC、ISDF、FCI、MBE-FCI、CCSDT、HS等缩略语。
   2.2 Diagonalizing the Interaction Hamiltonian（相互作用哈密顿量对角化）：研究人员借ITHC将ERI分解为等距形式，引入虚构模做正则变换，证明在原物理子空间投影下W?等价于原V?，且ITHC与THC同精度时辅助数N_aux相当。
   2.3 Extended Basis AFQMC（扩展基AFQMC）：虚时间步e^–ΔτH?写为e^–Δτh?/2P e^–ΔτW?P e^–Δτh?/2，Zeno误差O(Δτ)主导；W?经本征/Cholesky分解为Σ_γω?_γ²，HS变换得含N_aux维辅助场x的积分形式，B?(x)为单体算符；ph-AFQMC下walker按|?_i(τ+Δτ)?=B?(x_i–x?_i)|?_i(τ)?更新，权重乘I_ph=|S_ie^{x_i·x?_i–x?_i2/2}|max(0,cos(arg S_i))，力偏置x?_i,γ=–√Δτ Σ_αw_γα?Ψ_trial|n?_α|?_i?/?Ψ_trial|?_i?。
2. 3.
   Implementation Details（实现细节）
   3.1 Enlarged Basis Propagator（扩大基传播子）：单体e^–Δτh?/2在原基作用后，walker矩阵B⁽ⁱ⁾经B?_αk⁽ⁱ⁾=Σ_pu_pαB_pk⁽ⁱ⁾旋至扩展基，成本O(N_wN_auxNN_e)；扩展基中双体传播因仅含数算符，成本降至O(N_wN_auxN_e)。
   3.2 One-Body Green’s Function and Force Bias（单体格林函数与力偏置）：扩展基格林函数G?=(A?*)^T(B?^TA?)^–1B?^T，仅取对角算力偏置，总复杂度O(N_wNN_e²+N_wN_auxN_e)，内存O(N_aux²)，优于标准Cholesky AFQMC的O(N_wN_auxN²)与O(N_auxN²)。
   3.3 Energy Estimator（能量估计）：局域能量E_loc,i=?Ψ_trial|H?|?_i?/?Ψ_trial|?_i?，双体项?V??=?Σ_α≠βW_αβ(G?_ααG?_ββ–G?_αβG?_βα)，成本O(N_wN_aux²N_e)，低于标准法O(N_wN_auxNN_e²)。

讨论与结论翻译

讨论指出，ITHC-AFQMC因双体项在扩展基对角，使传播子与估计器均避开了Cholesky张量收缩，虽N_aux（ITHC秩）常大于Cholesky秩N_chol，但在GPU上线性氢链基准仍获更短运行时间；该方法在保持CC/DMRG级精度的同时改善标度，为强关联分子基态计算提供新路径。研究结论为：研究人员提出的扩展基ITHC-AFQMC通过对双电子相互作用在含虚构费米模的扩展空间做等距对角化，将双体项化为仅含数算符的Hubbard型形式，从而使传播与局域能量评估的理论复杂度分别降至O(N_wN_auxNN_e)与O(N_wN_aux²N_e)、内存降至O(N_aux²)；在线性H₁₀链与苯分子上恢复的基态能量与CCSDT及DMRG参照值一致，证明该法在精度可比前提下显著优于标准AFQMC与THC-AFQMC的实践效率。论文发表于《Journal of Chemical Theory and Computation》。