[论文解读] Associated Legendre Polynomials and Spherical Harmonics Computation for Chemistry Applications

Taweetham Limpanuparb, Josh Milthorpe|arXiv (Cornell University)|Oct 7, 2014

Scientific Research and Discoveries参考文献 2被引用 8

一句话总结

本文提出了一种高效且高精度的C++算法，用于在量子化学应用中计算连带勒让德多项式（ALPs）和球谐函数（SHs）。通过归一化ALPs，使用稳定的一阶递推关系计算三角函数，并预计算系数，该方法在度数和阶数高达1000时，实现了绝对误差或相对误差低于$10^{-10}$，单个ALP的平均计算时间少于4个周期。

ABSTRACT

Associated Legendre polynomials and spherical harmonics are central to calculations in many fields of science and mathematics - not only chemistry but computer graphics, magnetic, seismology and geodesy. There are a number of algorithms for these functions published since 1960 but none of them satisfy our requirements. In this paper, we present a comprehensive review of algorithms in the literature and, based on them, propose an efficient and accurate code for quantum chemistry. Our requirements are to efficiently calculate these functions for all non-negative integer degrees and orders up to a given number (<=1000) and the absolute or the relative error of each calculated value should not exceed 10E-10. We achieve this by normalizing the polynomials, employing efficient and stable recurrence relations, and precomputing coefficients. The algorithm presented here is straightforward and may be used in other areas of science.

研究动机与目标

• 开发一种在量子化学中计算连带勒让德多项式（ALPs）和球谐函数（SHs）的数值稳定且高效的算法。
• 确保在度数和阶数高达1000时，ALP和SH值的绝对或相对误差均低于$10^{-10}$。
• 通过使用三角函数的一阶递推关系和预计算系数来优化性能。
• 提供一种简单、可重用的实现，适用于超出量子化学领域的更广泛科学应用。

提出的方法

• 对ALPs进行归一化，以在递推计算中提高数值稳定性。
• 使用$\cos(m\phi)$和$\sin(m\phi)$的一阶递推关系，其数值稳定且相对误差低于$10^{-13}$。
• 预计算系数以避免重复计算，提升性能。
• 通过沿$\ell$递增方向的递推关系计算$P_{\ell}^{m}$值，起始于归一化的$P_{\ell}^{\ell}$。
• 应用反射公式，将$P_{\ell}^{-m}$与$P_{\ell}^{m}$关联，以处理负阶数。
• 最终通过归一化的$P_{\ell}^{m}$和预计算相位的三角函数计算球谐函数。

实验结果

研究问题

• RQ1 基于递推的ALPs和SHs算法能否在度数和阶数高达1000时实现低于$10^{-10}$的绝对或相对误差？
• RQ2 在此上下文中，使用一阶递推关系计算三角函数是否比两阶替代方法更具稳定性和准确性？
• RQ3 预计算系数和归一化能否显著提升ALP和SH评估的精度与性能？
• RQ4 当利用向量化和流水线技术时，该算法在$\ell = 1000$以上是否仍保持高效与准确？

主要发现

• 该算法在所有$\ell \leq 1000$和$|m| \leq \ell$的ALP和SH值计算中，实现了绝对或相对误差低于$10^{-10}$。
• $\cos(m\phi)$和$\sin(m\phi)$的一阶递推关系产生的相对误差低于$10^{-13}$，确保了SH计算的高精度。
• 现代处理器上单个ALP的平均计算时间少于4个周期，且随着$\ell$增大，由于有效循环向量化和流水线化，性能进一步提升。
• 当多次调用时，初始化时间可忽略不计（$L=100$时约100 μs），使该方法在重复使用中具有实际可行性。
• 球谐函数$Y_{\ell}^{m}$的误差受$P_{\ell}^{m}$和三角函数误差之和的限制，乘积项$\Delta_P \Delta_c$可忽略不计。