[論文レビュー] Disconnected quark loop contributions to nucleon observables using $N_f=2$ twisted clover fermions at the physical value of the light quark mass

Abdou Abdel-Rehim, Constantia Alexandrou+6 more|arXiv (Cornell University)|Nov 2, 2015

Particle physics theoretical and experimental studies参考文献 21被引用数 4

ひとこと要約

本研究では、物理的パイオン質量（$m_\pi \approx 133$ MeV）における$N_f=2$ねじれクローラーフェルミオンを用いて、核子の物理的観測量への非連結クォークループ寄与を計算した。GPUアクセラレーションされた格子上で、軽いクォークには正確なデフラーションを、 strange クォークには切断ソルバー法（TSM）を組み合わせることで、計算コストを低減しつつ高統計的結果を達成した。$\sigma_s$ および $g_A^s$ について、先行の物理的点への外挿結果と強い一致を示した。本手法により、物理的クォーク質量での直接的かつ高精度な計算が可能となり、核子構造の格子QCD計算が前進した。

ABSTRACT

We compute the disconnected quark loops contributions entering the determination of nucleon observables, by using a $N_f = 2$ ensemble of twisted mass fermions with a clover term at a pion mass $m_π= 133$ MeV. We employ exact deflation and implement all calculations in GPUs, enabling us to achieve large statistics and a good signal.

研究の動機と目的

• 物理的パイオン質量（$m_\pi \approx 133$ MeV）における核子観測量への非連結クォークループ寄与を計算すること。これは、格子QCDにおいて計算コストが非常に高いという主な課題である。
• すべてのクォーク伝播関数を評価する困難を、正確なデフラーションや切断ソルバー法（TSM）などの高度なアルゴリズム的手法を用いることで克服すること。
• GPUアクセラレーションと大統計を活用し、核子のスピン項（$\sigma_{\pi N}$, $\sigma_s$）、軸性結合定数（$g_A$）、テンソル結合定数（$g_T$）について、最小限の系統的誤差で高精度な結果を得ること。
• 結果を先行の物理的点への外挿値と比較することで、本手法の信頼性を検証し、有限体積効果や離散化効果の影響を暗黙的に評価すること。
• 物理的クォーク質量での直接的計算の可能性と正確性を示し、今後の研究におけるチャイral外挿の必要性を回避すること。

提案手法

• 物理的軽いクォーク質量に調整された、$a \approx 0.093$ fm および $m_\pi \approx 133$ MeV の $48^3 \times 96$ 格子上に $N_f=2$ ねじれクローラーフェルミオンを用いた。
• 軽いクォークの低モードを効率的に処理するため、ヒルベルト空間上の演算子 $Q = \gamma_5 M$ に対して PARPACK とランチョス固有値ソルバーを用いた正確なデフラーションを実装し、CGソルバーのコストを無視できるほどに削減した。
• Strange クォークループに切断ソルバー法（TSM）を適用：高精度（HP）63回と低精度（LP）1024回の反復を用い、計算コストを削減しながらも、確率的バイアス補正により精度を維持した。
• 各配置ごとに $N = 2250$ 個の $Z_4$ ノイズ源を用い、$\frac{1}{N}\sum_j |\eta_j\rangle\langle\eta_j| \approx \mathbb{I}$ により不偏推定量を保証し、体積和縮約により信号対ノイズ比を向上させた。
• QUDA ライブラリを用いて、すべての反復、縮約、FFT処理をGPU上で実装し、1配置あたり最大721,600回の測定を可能にした。これにより、効率的かつ大統計的計算が実現した。
• トレース縮約を計算するための恒等式を適用：$\textrm{Tr}[\Gamma(G_u - G_d)] = -\frac{2i\mu a}{N}\sum_j \langle s_j \Gamma \gamma_5 s_j \rangle$ および $\textrm{Tr}[\Gamma(G_u + G_d)] = \frac{2}{N}\sum_j \langle s_j \gamma_5 \Gamma \gamma_5 D_W s_j \rangle$ により、ゲージ不変性を保ちつつ統計を向上させた。

実験結果

リサーチクエスチョン

• RQ1 物理的パイオン質量における核子観測量への非連結クォークループ寄与を、$N_f=2$ ねじれクローラーフェルミオンを用いて高精度かつ十分な統計で直接計算できるか？
• RQ2 正確なデフラーションと切断ソルバー法（TSM）は、それぞれ軽いクォークループと strange クォークループにおいて、計算コストをどの程度低減し、信号対ノイズ比を向上させられるか？
• RQ3 核子スピン項（$\sigma_s$）および軸性結合定数（$g_A^s$）の結果が、先行の物理的点への外挿結果とどの程度一致するか。これは、本手法の信頼性を示す。
• RQ4 励起状態の汚染が非連結行列要素の決定に与える影響は何か。また、和集合法はその汚染をどの程度効果的に抑制できるか？
• RQ5 この直接的な物理的点計算において、有限体積効果や離散化効果は無視できる程度か。先行結果との一致に基づいて評価する。

主な発見

• 1配置あたり2250個の確率的ソースを用い、最大721,600回の測定を実施することで、非連結クォークループの高統計的計算を達成し、信号対ノイズ比を著しく向上させた。
• 軽いクォークループでは、正確なデフラーションにより600個の固有モードを計算することで、CG反復コストを無視できるほどに削減し、連結および非連結相関関数の効率的かつ正確な計算を可能にした。
• TSMは、63回の高精度（HP）および1024回の低精度（LP）ソースを用いて、strange クォークループに成功裏に適用し、計算コストを削減しながらも、確率的バイアス補正により精度を維持した。
• 和集合法は励起状態の汚染を効果的に抑制し、特にスピン項において顕著であった。これにより、大源・ソース間隔での安定したプラトー領域の特定が可能になった。
• 計算された $\sigma_s$ および $g_A^s$ の値は、先行の物理的点への外挿結果の重み付き平均と顕著な一致を示し、有限体積効果や離散化効果が小さいことが示された。
• 本結果は、GPUアクセラレーションされたアルゴリズムを用いた物理的クォーク質量での直接的格子QCD計算の可能性と正確性を示しており、今後の1次導関数型演算子の研究や、複数のエンsemblesにおける系統的誤差解析への道筋を示した。