[論文レビュー] New developments in PJFry

J. Fleischer, T. Riemann+1 more|DESY (CERN, DESY, Fermilab, IHEP, and SLAC)|Oct 15, 2012

Black Holes and Theoretical Physics参考文献 29被引用数 6

ひとこと要約

この論文では、量子場理論における1ループテンソル積分の数値安定な計算を可能にする、オープンソースのC++パッケージPJFryを提示する。次元的再帰と級数展開を用いることで、逆グラム行列式を回避し、ランク6までのペンタゴンおよびヘキサゴン積分の堅牢な評価を可能にし、高エネルギー物理学の振幅計算や自動化された振幅計算に応用可能である。

ABSTRACT

We report on recent progress in numerical evaluation of one loop tensor integrals. A public C++ package PJFry implementing algorithms from arXiv:1009.4436 and its extension to hexagons up to rank 6 are presented.

研究の動機と目的

• 多腿過程における逆グラム行列式に起因する1ループテンソル積分の数値不安定性を解消すること。
• 4点関数を超える高ランクテンソル積分を対象とする、堅牢で公開可能なツールの開発。
• 任意の質量とランクを想定した、ペンタゴンからヘキサゴンまでの高ポイント関数へのテンソル積分還元の適用範囲の拡張。
• 小グラム行列式領域におけるテンソル型係数の正確な評価を可能にする漸近級数展開の利用。
• 自動化された振幅計算やループ振幅ライブラリに利用可能な、高性能で再利用可能なフレームワークの提供。

提案手法

• パッケージは、テンソル積分を、分母のべきをずらした高次元のスカラー積分に表現する次元的再帰関係を用いる。
• 式(6)および(7)の再帰関係を用いて、テンソル型係数を高次元のスカラー積分に還元し、グラム行列式の直接逆算を回避する。
• グラム行列式が小さい領域では、式(10)から導かれる再帰関係に基づく級数展開を適用し、数値的安定性を維持する。
• 組合せ的係数およびべき乗の上昇・低下作用素を介して、テンソル型係数を一般化されたスカラー積分にマッピングする。
• キャッシュ機構により、低次元のテンソル積分を再利用することで、計算を高速化し、重複評価を削減する。
• C、FORTRAN、C++、Mathematica、GoSamなど複数のインターフェースをサポートし、複素質量のための拡張精度をqdライブラリで実現する。

実験結果

リサーチクエスチョン

• RQ1 小またはゼロに近いグラム行列式が存在する状況下でも、1ループテンソル積分をどのように安定して評価できるか。
• RQ2 5点および6点関数の高ランクテンソル積分を効率的かつ数値的に安定して還元する最良の方法は何か。
• RQ3 次元的再帰と級数展開を組み合わせることで、逆行列式を導入せずに小グラム領域での精度を維持できるか。
• RQ4 任意の質量配置を想定した、ランク6までのテンソル積分を扱える汎用的で高性能なライブラリをどのように設計できるか。
• RQ5 テンソル積分還元における中間結果を最適にキャッシュ・再利用する戦略は何か。

主な発見

• PJFryは、小グラム行列式領域を含むすべての相空間領域において、ランク6までの5点および6点関数のテンソル積分を完全に数値安定的に評価できる。
• 再帰に基づく還元と漸近級数展開により、逆グラム行列式を回避することで、あらゆる相空間領域で堅牢性を確保する。
• 低次元テンソル積分の再利用を可能にするキャッシュ機構のおかげで、高い性能を達成している。
• 6外部粒子と1本の質量のある内部線を含むテスト相空間点において、スカラー・ヘキサゴン $F_0$ は $\epsilon=0$ で $-30.3579 - 205.8213i$ に評価され、$\epsilon=0,1,2$ で一貫した結果が得られた。
• 同じフレームワークを用いて、ベクトル・ヘキサゴン成分も正確に計算されており、複素平面における期待される挙動と一致している。
• Mathematicaインターフェースは、LoopToolsの標準的表記をサポートし、記号的振幅計算への直接的な利用を可能にしている。