[论文解读] Scalable and predictive spectra of correlated molecules with moment truncated iterated perturbation theory

Oliver J. Backhouse, Alejandro Santana‐Bonilla+1 more|arXiv (Cornell University)|Jun 7, 2021

Spectroscopy and Quantum Chemical Studies参考文献 94被引用 19

一句话总结

该论文提出了一种可扩展、可预测且自洽的多体微扰方法——矩量截断迭代微扰理论（AGF2），用于计算关联分子的完整单粒子谱函数。通过利用谱矩守恒压缩有效自能动力学，该方法实现了 O[N⁵] 的计算量，支持大规模并行计算（最多约 1000 个轨道），并准确解析了复杂体系（如青蒿素）的电子激发谱，解决了长期存在的其最低能量激发态归属争议。

ABSTRACT

A reliable and efficient computation of the entire single-particle spectrum of correlated molecules is an outstanding challenge in the field of quantum chemistry, with standard density functional theory approaches often giving an inadequate description of excitation energies and gaps. In this work, we expand upon a recently-introduced approach which relies on a fully self-consistent many-body perturbation theory, coupled to a non-perturbative truncation of the effective dynamics at each step. We show that this yields a low-scaling and accurate method across a diverse benchmark test set, capable of treating moderate levels of strong correlation effects, and detail an efficient implementation for applications up to $\sim1000$ orbitals on parallel resources. We then use this method to characterise the spectral properties of the artemisinin anti-malarial drug molecule, resolving discrepancies in previous works concerning the active sites of the lowest energy fundamental excitations of the system.

研究动机与目标

• 开发一种可扩展且精确的方法，用于计算强关联分子的完整单粒子格林函数谱。
• 解决标准 DFT 和微扰方法在描述关联体系中激发能和带隙时的局限性。
• 实现大规模模拟（最多约 1000 个轨道），并具备高效的并行实现，适用于复杂分子体系。
• 解决文献中关于青蒿素抗疟药物分子最低能量激发态空间局域化归属的争议。

提出的方法

• 该方法采用完全自洽的辅助格林函数方法（AGF2），通过辅助自由度将二阶自能重新表述为静态特征值问题。
• 采用矩量截断近似，保留完整二阶自能的一阶谱矩，确保精度的同时降低计算成本。
• 有效自能表示为极点之和的形式：Σpq(ω) = Σk (vpk v†qk)/(ω − Ek)，从而将狄拉克方程转化为矩阵特征值形式。
• 该算法在开源 PySCF 框架中实现，支持高效并行化和大规模应用。
• 该方法避免了将粒子和空穴激发分离，从而可同时获取所有能量尺度上的完整激发谱。
• 该方法被证明等价于 ADC(2) 方法的自洽扩展，由于对动力学自能分量的系统性截断，谱精度得到提升。

实验结果

研究问题

• RQ1 矩量截断的自洽多体微扰方法是否能在关联分子谱中同时实现高精度和低计算量？
• RQ2 与 ADC(2)、GW 和 EOM-CCSD 等成熟方法相比，矩量截断 AGF2 方法在多样化分子体系中的谱精度如何？
• RQ3 该方法在多大程度上能够解决复杂分子（如青蒿素）中激发态特征归属的模糊性？
• RQ4 在并行架构上，该方法在保证可靠谱精度的前提下，可处理的最大体系尺寸（轨道数）是多少？
• RQ5 与完整动力学自能描述相比，矩量截断策略是否能提升关联体系中的谱分辨率？

主要发现

• AGF2 方法实现了 O[N⁵] 的计算量，具有高效的并行性，支持最多约 1000 个轨道的体系模拟。
• 在 GW100 基准集上，AGF2 对带电激发态的计算精度具有竞争力，与 ADC(2) 和 GW 相当，且比 EOM-CCSD 更具可扩展性。
• 该方法成功解决了文献中关于青蒿素最低能量激发态归属的矛盾，确认活性位点为过氧桥结构，与实验光电子能谱结果一致。
• 矩量截断的自能近似保留了关键谱特征，并在精度上优于非截断的自洽方法，尽管其违反了严格的守恒近似。
• 在 PySCF 中的实现支持可重复的大规模计算，并可直接获取完整的能量依赖谱函数。
• 由于完全自洽性，该方法对参考态的选择不敏感，增强了在强关联区域的鲁棒性。