[論文レビュー] Multi-gluon one-loop amplitudes numerically

Achilleas Lazopoulos|ArXiv.org|Dec 16, 2008

Particle physics theoretical and experimental studies参考文献 22被引用数 24

ひとこと要約

本稿では、NLO QCDにおける1ループ多ゲルミオン振幅を数値的に計算するための $D_s$ 次元ユニタリティ切断アルゴリズムのC++実装を提示する。この手法は、5次元におけるゲルミオンおよびスカラーを含む木レベルの色順序付き振幅を構築ブロックとして用いる。本手法により、仮想振幅の有理部および切断構成可能部の両方を正確に評価でき、四重精度のフォールバックにより数値的安定性が保証され、最大22外部ゲルミオンまでに文献結果と完全に一致する。

ABSTRACT

A c++ implementation of the D_s-dimensional unitarity cut algorithm for the numerical evaluation of the virtual contribution to NLO QCD amplitudes is presented. The current version includes an arbitrary number of external gluons with gluonic propagators in the loop. The building blocks are tree level color-ordered amplitudes with gluons and with gluons and two scalars in five dimensions. Numerical stability issues are addressed and agreement has been reached with the results published in the literature.

研究の動機と目的

• 任意の外部ゲルミオン数を想定した、1ループ仮想振幅を数値的に安定して計算できる汎用的なC++ツールの開発。
• $D_s$ 次元ユニタリティ切断法（EGKMアルゴリズム）の実装により、振幅の有理部および切断構成可能部の両方を効率的かつ正確に評価する。
• 四重精度算術をフォールバックメカニズムとして統合することで、高多重度プロセスにおける数値的ロバスト性を確保する。
• 同じフレームワーク内でフェルミオンおよびゲージボソンへの将来的な拡張を可能とし、LHCにおける物理学的関連プロセスをサポートする。

提案手法

• アルゴリズムは、被積分関数レベルでのOPP還元法を用い、1ループ振幅をスカラー四角形、三角形、バブル、タドル積分に分解する。
• 仮想振幅は $A^{D_s} = A^0 + A^1 \cdot D_s$ と表現され、内部粒子の次元 $D_s$ を2つの値で評価することで、$\epsilon$-依存部および有理部を抽出可能となる。
• ゲルミオンおよびゲルミオン＋2スカラーを含む木レベルの色順序付き振幅は、Berends-Giele再帰関係を用いて計算される。
• 数値的安定性を維持するため、QDライブラリを用いた四重精度版のコードを用いて、$D_s$ 次元における複素ユニタリティ切断を実施する。
• アルゴリズムは、任意の外部ヘリシティおよびフェーズスペース点を想定したNゲルミオン振幅に適用され、$10^5$ 個のフェーズスペース点における平均CPU時間を性能指標として測定される。
• 特定のフェーズスペース点に対して文献結果と照合し、数値精度内での一致を確認する。

実験結果

リサーチクエスチョン

• RQ1 NLO QCDにおける多ゲルミオン過程において、有理部を含む完全な仮想1ループ振幅をどのように数値的に計算できるか？
• RQ2 多くの外部粒子を含む高多重度振幅において、数値的安定性を維持するための有効な数値戦略は何か？
• RQ3 $D_s$ 次元ユニタリティ切断法を、汎用的なC++フレームワーク内で効率的かつロバストに実装できる範囲はどの程度か？
• RQ4 アルゴリズムの性能は外部ゲルミオン数にどのように依存するか？また、物理学的に関連するプロセスを実行するための計算リソースはどの程度必要か？

主な発見

• C++実装は、最大22外部ゲルミオンまで完全な1ループ振幅を計算可能であり、Intel Xeon X5450上で4ゲルミオンでは3.4 ms、22ゲルミオンでは1.85秒の典型的なCPU時間を記録した。
• N=6ゲルミオンの場合、スイッチ精度を $10^{-4}$ に設定したとき、約5%のフェーズスペース点で四重精度算術が必要であった。これは性能への影響が管理可能であることを示している。
• N=6の場合の二重極および単一極における相対誤差は対数分布を示し、Nが増加するにつれて誤差のピークがより大きな値にシフトしたが、有限部は大多数のケースで安定していた。
• コードは文献からの結果と完全に一致しており、$D_s = 5$ および $6$ での照合も実施され、実装の正しさが確認された。
• 振幅の $\epsilon$ 展開における $1/\epsilon^2$ および $1/\epsilon$ の極、および有理有限部が正しく回復された。
• フレームワークはフェルミオンおよびゲージボソンへの拡張が可能であり、先行研究で示されたように、今後は完全な物理学的実装への道筋が確立された。