[論文レビュー] Stochastic evaluation of four-component relativistic second-order many-body perturbation energies: A potentially quadratic-scaling correlation method

Julio César Cruz, Jorge Garza+2 more|arXiv (Cornell University)|Mar 10, 2022

Nuclear physics research studies参考文献 101被引用数 2

ひとこと要約

本論文は、実時間および虚時間の両方におけるモンテカルロ積分を用いて、計算コストの高い四添い積分変換を回避する確率的で4成分の相対論的2次多体摂動（MP2）手法を導入する。この手法は、系のサイズに対して潜在的に2次式スケーリング（O(n²)）を達成しており、従来のO(n⁵)スケーリングと比べて顕著に改善され、高い精度を維持しながら4096プロセッサで92％の強スケーリングを実現する。

ABSTRACT

A second-order many-body perturbation correction to the relativistic Dirac-Hartree-Fock energy is evaluated stochastically by integrating 13-dimensional products of four-component spinors and Coulomb potentials. The integration in the real space of electron coordinates is carried out by the Monte Carlo (MC) method with the Metropolis sampling, whereas the MC integration in the imaginary-time domain is performed by the inverse-CDF (cumulative distribution function) method. The computational cost to reach a given relative statistical error for spatially compact but heavy molecules is observed to be no worse than cubic and possibly quadratic with the number of electrons or basis functions. This is a vast improvement over the quintic scaling of the conventional, deterministic second-order many-body perturbation method. The algorithm is also easily and efficiently parallelized with demonstrated 92% strong scalability going from 64 to 4096 processors for a fixed job size.

研究の動機と目的

• 強い電子相関および相対論的効果を示す重元素化合物に対して、スケーラブルな電子相関法を開発すること。
• 従来の4成分相対論的MP2におけるO(n⁵)計算コストのボトルネックを解消するため、四添い積分変換を排除すること。
• 確率的積分が、局所相関法が失敗する空間的に凝縮した重い系において、スケーリングを低減（潜在的にO(n²)）できることを示すこと。
• 高性能コンピューティングを活用して、多数の電子および基底関数を有する系に対して大規模な相対論的相関計算を可能にすること。

提案手法

• 4成分スピンダーとクーロンポテンシャルを含む13次元積分を、モンテカルロ（MC）積分による確率的評価に用いる。
• 実空間の電子座標統合にはメトロポリスサンプリングを、虚時間の統合には逆CDFサンプリングを用いる。
• 決定的四添い積分変換を確率的サンプリングに置き換え、大規模な積分の保存および変換を回避する。
• 独立したMCサンプリングに内在する並列性を活用し、数千のプロセッサに効率的に分散可能である。
• 統計誤差の低減と収束の向上を図るため、重み関数を用いた重要度サンプリングを用いる。
• 相対論的効果を初期段階から組み込むために、ディラック–ハートリー–フォック基底法と統合する。

実験結果

リサーチクエスチョン

• RQ1 確率的積分は、重元素系において4成分相対論的MP2のスケーリングをO(n⁵)から立方未満に低減できるか？
• RQ2 空間的に凝縮した重い分子において、確率的アプローチが著しく低い計算コストで高い精度を維持できるか？
• RQ3 数千のCPUやGPUなどの大規模並列アーキテクチャにおいて、この手法は効率的にスケーリングできるか？
• RQ4 さまざまな系のサイズおよび基底関数サイズにおいて、統計誤差は制御可能で予測可能か？
• RQ5 局所相関近似が失敗する系、例えばランタノイドやアクチノイドにこの手法を適用できるか？

主な発見

• 計算コストは、電子数や基底関数数に対して立方未満、あるいは潜在的にO(n²)まで低減しており、従来のO(n⁵)MP2と比べて顕著な改善である。
• 64から4096プロセッサにスケーリングする際、92％の強スケーリングを達成しており、優れた並列効率を示している。
• 統計誤差は制御可能であり、サンプリングを増やすことで減少し、高精度な結果が得られる。
• 確率的アプローチにより、四添い積分の保存や変換の必要性を完全に回避でき、相対論的相関計算における主要なボトルネックを解消した。
• 従来の局所相関法が失敗する空間的に凝縮した重い系、例えばアクチノイドやランタノイドに対しても、本手法は適用可能である。
• 数値的証拠により、決定的アプローチと比較して著しく低い計算コストで高い精度を達成できることが示された。