[论文解读] Various integral estimations and screening schemes for extended systems in PySCF

Qiming Sun|arXiv (Cornell University)|Feb 22, 2023

Electromagnetic Scattering and Analysis被引用 5

一句话总结

本论文利用 PySCF 的范围分离和傅里叶变换积分算法，针对扩展体系中的双电子积分，提出了分析导出的上界和筛选参数。该研究确定了确保所需精度的同时避免不必要的计算成本的距离和能量截断值，数值测试结果证实了其在多种晶体体系、基组和参数下的鲁棒性能。

ABSTRACT

In this document, we briefly review the two-electron integral algorithms based on the range-separated algorithms and Fourier transform integral algorithms that are implement ed in the PySCF package. For each integral algorithm, we estimate the upper bound of relevant integrals and derive the necessary conditions for the screening parameters, including distance cutoff and energy cutoff, to reach the desired accuracy. Our numerical tests show that the proposed integral estimators and screening parameters can effectively address the required accuracy while computational efforts are not wasted on unintended accuracy.

研究动机与目标

• 利用 PySCF 开发周期性体系中双电子积分的可靠积分估计与筛选方案。
• 推导出分析上合理的距离和能量截断值，以确保积分误差在用户指定容差范围内。
• 通过避免高估（浪费资源）和低估（损害精度）来确保计算效率。
• 验证筛选参数在不同体系参数（如晶格尺寸、k 点网格、基组紧凑性及库仑衰减）下的鲁棒性。
• 通过可靠且具有定量基础的截断估计器，实现对晶体积分算法更激进的优化。

提出的方法

• 利用改进的 Schwarz 类估计器推导双电子积分的上界，其中引入了电荷密度间距因子。
• 将这些上界应用于推导范围分离和傅里叶变换积分算法中所需的距离截断 $ R_{\text{cut}} $ 和能量截断 $ E_{\text{cut}} $ 的必要条件。
• 使用变换方程 (33) 和 (36)，将实空间晶格求和与平面波展开中长程和短程部分的截断相互映射。
• 对短程 ERI 采用三重晶格求和，对长程部分采用平面波求和，通过调节截断值以满足目标精度 $ \tau $。
• 通过系统性数值测试验证估计器性能，涵盖不同晶格参数、k 点网格、基函数及库仑衰减参数。
• 调节参数如 $ R_{\text{cut}} $、$ E_{\text{cut}} $ 和 $ \omega $，以在晶体积分计算中实现精度与计算成本的平衡。

实验结果

研究问题

• RQ1 在扩展体系中，能否推导出双电子积分的分析上界，以指导高效筛选？
• RQ2 如何系统性地确定距离和能量截断值，以确保在用户定义容差 $ \tau $ 内的精度，同时避免高估？
• RQ3 晶格尺寸、k 点网格、基组紧凑性及库仑衰减等体系参数如何影响积分筛选性能？
• RQ4 原始积分中的数值不确定性在多大程度上限制了可实现的精度，尤其是在小晶胞或使用弥散基函数时？
• RQ5 所推导的截断估计器能否在 PySCF 中跨多种积分算法（如 RSDF、CCDF、RSJK、傅里叶变换）实现一致的精度？

主要发现

• 所推导的截断估计器可可靠地实现目标精度 $ \tau $，误差通常接近或略低于 $ \tau $，表明计算资源得到高效利用。
• 对于小晶胞（如 $ a = 1.0 $ Å），即使 $ \tau = 10^{-11} $，精度仍受限于 $ \sim 10^{-10} $，可能源于原始积分库的数值不确定性。
• 弥散基函数（$ \alpha = 0.2 $）的精度在 $ 10^{-9} $ 左右趋于平台，表明原始积分评估中存在数值不稳定性或精度限制。
• 库仑衰减参数 $ \omega < 0.5 $ 导致较大误差，精度上限约为 $ 10^{-10} $，原因可能是数值噪声；而较大的 $ \omega $ 值则显示出截断值的轻微高估。
• 筛选方案在 $ N_k = 1^3 $ 到 $ 4^3 $ 的 k 点网格上表现一致，所有配置均达到所需精度。
• 角动量和基函数紧凑性对截断估计影响极小，表明该方法在各类基组中均具有鲁棒性。