[論文レビュー] Electric polarizability of hadrons with overlap fermions on multi-GPUs

Michael Lujan, Andrei Alexandru+1 more|arXiv (Cornell University)|Nov 27, 2011

Quantum Chromodynamics and Particle Interactions参考文献 17被引用数 3

ひとこと要約

本稿では、格子QCDを用いてハドロンの電気感受率、特に中性π中間子の電気感受率を計算するための高性能GPU実装であるオーバーラップフェルミオン形式を提示する。マルチGPUアーキテクチャとデフラクションを施した適応的共役勾配ソルバーを活用することで、1GPUあたり約30CPUコア相当の性能を達成し、クリティカルなクリスタル的ダイナミクスが働く低π中間子質量領域でのシミュレーションが可能になった。予備的な結果では、500 MeV近辺で中性π中間子の感受率の符号が反転することが示され、理論的予想と整合的である。

ABSTRACT

Electric polarizability is an important parameter for the internal structure of hadrons. Previous studies of polarizabilities have been done at relatively heavy pion masses, leaving the chiral region largely unexplored. In this report, we use overlap fermions which are known to be computationally demanding to properly capture the chiral dynamics. We present an implementation strategy to construct overlap on multi-GPUs. We find that our GPU code has an equivalent of \sim30 CPU cores to 1 GPU. We also present preliminary results for the polarizability of the neutral pion.

研究の動機と目的

• チャーミカル対称性を保存し、例外的配置を回避するオーバーラップフェルミオンを用いてハドロンの電気感受率を計算すること。
• オーバーラップフェルミオンの計算ボトルネックを、マルチGPUアーキテクチャへの効率的実装によって克服すること。
• 感受率が非自明な振る舞いを示すと予想される低π中間子質量（約240 MeVまで）でのシミュレーションにより、QCDのクリスタル的領域を探索すること。
• 大規模な格子QCDシミュレーションに適した、スケーラブルでメモリ効率の良いGPU戦略を開発すること。この戦略は、適応的マルチマスCGソルバーと固有値系デフラクションを組み合わせる。

提案手法

• マトリクスの符号関数をZolotarevの有理近似により近似し、多項式近似を小さな固有値部分空間と大きな固有値部分空間に分割することで、オーバーラップフェルミオン作用を実装する。
• 固有値系ソルバーは、$ H_w $ および $ D_{\text{ov}} $ の低エネルギー固有モードを計算するための暗黙的再起動Arnoldi因子分解を用い、収束を加速するデフラクションを可能にする。
• マルチマスで適応的な共役勾配（CG）法を用いてオーバーラップディラック作用素を逆行列化し、多項式の次数と精度（単精度算術を含む）を動的に調整して性能最適化を行う。
• バックグラウンド場法を適用する際には、ゲージリンクに位相因子を導入して定常電場を導入し、共変微分を $ D_\mu = \partial_\mu - iG_\mu - iqA_\mu $ に修正する。
• 5つのクォーク伝播関数（$ u $ および $ d $ を別々に含む）を1配置あたり用い、ソース位置を変更することで統計的精度を向上させ、電場あり・なしの相関関数を計算する。
• メモリ制限に対処するため、データを複数のGPUに分散させ、GPUメモリが不足した場合にはハイブリッドCPU-GPUメモリ戦略を用いる。

実験結果

リサーチクエスチョン

• RQ1 中性π中間子の電気感受率は、低π中間子質量、特にクリスタル極限近辺でどのように変化するか？
• RQ2 計算コストが非常に高いオーバーラップフェルミオン形式を、マルチGPUアーキテクチャに効率的に実装することで、物理的π中間子質量でのシミュレーションを可能にする方法は何か？
• RQ3 チャーミカル有効場理論の予測通り、$ m_\pi \approx 500 \, \text{MeV} $ の領域で中性π中間子感受率の符号が反転するか？
• RQ4 デフラクションと適応的精度CG法が、GPU上でのオーバーラップフェルミオン伝播関数の計算コストをどの程度低減できるか？

主な発見

• GPU実装により、1GPUあたり約30CPUコア相当の性能を達成し、オーバーラップフェルミオン計算が顕著に高速化された。
• 適応的マルチマスCGに動的精度チューニングを適用した結果、$ m_\pi \approx 200 \, \text{MeV} $ での1伝播関数計算時間は、非適応的の場合の1.6時間から0.8時間に短縮された。
• 中性π中間子感受率は $ m_\pi \approx 500 \, \text{MeV} $ で符号を反転し、理論的予測および先行研究と整合的である。
• 統計的精度が限られているため誤差の範囲が大きい。低質量領域での信号対雑音比を向上させるために、4倍の配置を生成中である。
• 16GPUに固有値系および伝播関数計算を分散させ、必要に応じてハイブリッドCPU-GPUメモリ使用戦略を採用することで、メモリ制限を効果的に処理した。
• 本手法により、π中間子質量が300 MeV未満の領域で、オーバーラップフェルミオンを用いた中性π中間子感受率の格子QCD計算が初めて実現された。物理的点に近い領域での計算が可能になった。