[论文解读] Alchemical perturbation density functional theory (APDFT)

von Rudorff, Guido Falk, von Lilienfeld, O. Anatole|arXiv (Cornell University)|Sep 5, 2018

Advanced Chemical Physics Studies被引用 32

一句话总结

本文提出了一种无轨域的密度泛函理论方法——化学势扰动密度泛函理论（APDFT），该方法通过在化学势扰动（如核电荷变化）下对电子密度导数进行泰勒级数展开，预测等电子分子的电子性质。该方法仅需一个参考电子密度及其扰动项，即可在极低计算成本下实现与PBE0等高阶方法相当的精度。

ABSTRACT

We introduce an orbital free electron density functional approximation based on alchemical perturbation theory. Given convergent perturbations of a suitable reference system, the accuracy of popular self-consistent Kohn-Sham density functional estimates of properties of new molecules can be systematically surpassed---at negligible cost. The associated energy functional is an approximation to the integrated energy derivative, requiring only perturbed reference electron densities: No self-consistent field equations are necessary to estimate energies and electron densities. Electronic ground state properties considered include covalent bonding potentials, atomic forces, as well as dipole and quadropole moments.

研究动机与目标

• 开发一种计算高效的电子性质预测方法，用于等电子分子。
• 仅利用参考电子密度及其导数，超越标准DFT的精度。
• 在性质预测中消除自洽场迭代的需要。
• 通过高阶扰动项实现系统性精度提升。
• 通过一次高精度参考计算，实现对大规模化学空间的筛选。

提出的方法

• 使用化学势扰动理论，将核电荷变化视为连续参数λ。
• 在参考体系（λ = 0）附近对能量和电子密度进行泰勒级数展开。
• 应用Hellmann-Feynman定理，将能量导数表示为外部势差与密度导数乘积的积分。
• 通过密度导数的求和构建λ平均电子密度˜ρ，从而实现能量预测：∫∆v·˜ρ dr。
• 利用有限差分或参考密度的解析导数，推导密度和能量导数。
• 截断级数至少数项，对小ΔZ（如ΔZ = 1）系统在固定几何结构下实现收敛。

实验结果

研究问题

• RQ1 是否可利用化学势扰动理论在无需自洽场迭代的情况下预测电子性质？
• RQ2 是否可仅通过参考密度导数系统性地提升DFT类性质的精度？
• RQ3 λ平均密度˜ρ是否在等电子体系中收敛并能提供准确的能量与性质？
• RQ4 该方法能否以高保真度预测力、偶极矩和四极矩？
• RQ5 一次高精度参考计算是否可实现对广阔化学空间的可靠预测？

主要发现

• 对于以N2为参考的CO和BF体系，APDFT在3–4阶扰动下，能量、力和电荷性质预测精度与PBE0相当。
• 对于He体系，使用HF或DFT泛函时，n = 4阶展开下能量误差低于1 mH。
• 该方法预测的偶极矩和四极矩精度与参考方法一致。
• 在类氢原子和自由原子体系中观察到收敛性，且小ΔZ体系的数值证据支持收敛性。
• 该方法可通过化学势梯度实现几何优化和化学组成扫描，更高阶项可将适用范围扩展至更大的化学空间区域。
• 仅需参考电子密度及其导数，因此可应用于任何输出电子密度的参考方法。