[論文レビュー] 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

ひとこと要約

この論文では、量子化学的応用における関連ルジャンドル多項式（ALPs）および球面調和関数（SHs）を計算する、非常に精度が高く効率的なC++アルゴリズムを提示する。ALPsの正規化、三角関数の安定な1項再帰関係の使用、および係数の事前計算により、次数および階数が1000未満の場合に絶対誤差または相対誤差が$10^{-10}$未満の精度を達成し、ALP1つあたりの平均計算時間は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}$未満に保たれることの保証。
• 三角関数の1項再帰関係の使用と係数の事前計算により、パフォーマンスの最適化。
• 量子化学を越えた広範な科学的応用に適した、直感的で再利用可能な実装の提供。

提案手法

• 再帰計算中の数値的安定性を向上させるために、ALPsの正規化。
• $\cos(m\phi)$および$\sin(m\phi)$の1項再帰関係の使用。これらは数値的に安定しており、相対誤差が$10^{-13}$未塔であることが示されている。
• 繰り返し計算を避けるために係数を事前計算し、パフォーマンスを向上。
• 正規化された$P_{\ell}^{\ell}$から始まり、$\ell$を増加する方向への再帰により$P_{\ell}^{m}$の値を計算。
• 負の階数のための反射公式を用いて、$P_{\ell}^{-m}$を$P_{\ell}^{m}$に関連付ける。
• 事前計算された位相を用いた正規化された$P_{\ell}^{m}$および三角関数を用いて、球面調和関数を最終的に計算。

実験結果

リサーチクエスチョン

• RQ1 ALPsおよびSHsの再帰ベースのアルゴリズムは、次数および階数が1000未満の場合に、$10^{-10}$未塔の絶対誤差または相対誤差を達成できるか？
• RQ2 この文脈において、2項再帰関係の代替として、三角関数の1項再帰関係がより安定で正確であるか？
• RQ3 事前計算された係数と正規化は、ALPおよびSHの評価における精度とパフォーマンスを顕著に向上させられるか？
• RQ4 ベクトル化およびパイプライン処理を活用した場合、$\ell = 1000$を超えた範囲でもアルゴリズムは効率的かつ正確に保たれるか？

主な発見

• すべてのALPおよびSHの値について、$\ell \leq 1000$および$|m| \leq \ell$の範囲で、絶対誤差または相対誤差が$10^{-10}$未塔であることが達成された。
• $\cos(m\phi)$および$\sin(m\phi)$の1項再帰関係は、相対誤差が$10^{-13}$未塔を示し、SH計算における高い正確性を保証している。
• 現代のプロセッサ上では、1つのALPを計算する平均時間は4サイクル未塔であり、ループのベクトル化およびパイプライン処理のおかげで、$\ell$が大きくなるほどパフォーマンスが向上する。
• 複数回の呼び出しにわたって平均化した場合、初期化時間はわずか（$L=100$の場合約100 μs）であり、繰り返し使用に実用的である。
• 球面調和関数$Y_{\ell}^{m}$の誤差は、$P_{\ell}^{m}$と三角関数の誤差の和によって上限が定まり、積項$\Delta_P \Delta_c$は無視できるほど小さい。