[論文レビュー] New Developments in FormCalc 4.1

Thomas Hahn|ArXiv.org|Jun 21, 2005

Parallel Computing and Optimization Techniques被引用数 8

ひとこと要約

FormCalc 4.1 は、素粒子物理学における自動的一次ループ行列要素計算のための主要な向上を実装した。具体的には、フェルミオン振幅計算の効率化を図るウェイル・ヴァン・ダール・ヴァーデンスピンル形式の導入、多次元位相空間積分のためのCubaライブラリ連携、最適化されたFortranコード生成、およびMSSMにおけるFeynHiggsおよび非最小なフラバー違反の新インターフェースを含む。これらの改善により、ヒッグス粒子生成やトップクォーク崩壊のような過程における精度計算の性能、数値的安定性、および使いやすさが顕著に向上した。

ABSTRACT

FormCalc is a matrix-element generator that turns FeynArts amplitudes up to one loop into a Fortran code for computing the squared matrix element. The generated code can be run with FormCalc's own driver programs or used with other `frontends', e.g. Monte Carlos. Major new or enhanced features in Version 4.1 are: treatment of external fermions, phase-space integration, code-generation functions, extensions for the MSSM, the HadCalc frontend.

研究の動機と目的

• 外部フェルミオンを含む一回りの行列要素計算の効率性と数値的安定性の向上を図ること。
• FeynHiggs 2.2 との連携により、MSSMにおける高精度なヒッグス質量計算を可能にすること。
• 上・下クォーク型スクォークの完全な6×6混合行列を用いたMSSMにおける非最小なフラバー違反（NMFV）のサポートを可能にすること。
• 公開されたFortranコード生成関数と新規フロントエンドを通じて、コード生成および統合ワークフローの簡素化を図ること。
• Mersenne Twisterサンプリングとグリッドメモリゼーションを備えたCuba 1.2ライブラリを用いた、高度なモンテカルロおよび決定的アルゴリズムによる位相空間積分能力の強化。

提案手法

• FormCalcにウェイル・ヴァン・ダール・ヴァーデンスピンル形式を実装し、トレース手法を用いずに振幅を直接計算することで、フェルミオンチェーン数nに対して計算複雑度をO(n²)からO(n)に低減する。
• 多次元位相空間積分のためのCuba 1.2ライブラリを統合し、Vegas、Suave、Divonne、Cuhreアルゴリズムをサポート。Mersenne Twisterサンプリングとグリッドメモリゼーションの機能を備える。
• 最適化されモジュラーなFortranコードを生成するための公開されたMathematica関数を導入。共通部分式の削減と、大きな式を複数ファイルに分割する機能を含む。
• プリプロセッサフラグを介したFeynHiggs 2.2との連携により、完全な一次ループ精度で放射修正を含むヒッグス質量を計算可能にする。
• モデルファイルと数値評価の拡張により、上・下クォーク型スクォークの完全な6×6混合行列を用いたMSSMにおけるNMFVのサポートを実現。
• PDF畳み込み、運動量のカット、インバリアント質量・ラピディティ・横運動量に関する微分分布を備えた、自動的ハドロン断面積計算を可能にする新規HadCalcフロントエンドの提供。

実験結果

リサーチクエスチョン

• RQ1 外部フェルミオンを含む行列要素の二乗の計算コストを、正確性を損なわずにどのように低減できるか？
• RQ2 高度なモンテカルロおよび決定的アルゴリズムを用いることで、位相空間積分の性能と安定性にどのような向上が達成できるか？
• RQ3 高エネルギー物理学の計算において、大きな式やモジュラーコンパイルを想定した自動Fortranコード生成をどのように最適化できるか？
• RQ4 FeynHiggsを介した完全な一次ループヒッグス質量補正の組み込みが、MSSMにおける精度予測に与える影響は何か？
• RQ5 非最小なフラバー違反を自動行列要素生成器に体系的に組み込むにはどうすればよいか？

主な発見

• ウェイル・ヴァン・ダール・ヴァーデン形式により、外部フェルミオンを含む振幅の縮約回数がO(n²)からO(n)に削減され、多数の外部フェルミオンを持つ振幅の効率的計算が可能になった。
• Mersenne Twisterサンプリングとグリッドメモリゼーションを備えたCuba 1.2連携により、特に多次元位相空間において収束性とロバスト性が向上した。
• 新規のFortranコード生成関数により、共通部分式の削減と複数ファイルへの分割をサポートする最適化されモジュラーなコードが生成可能になった。
• FeynHiggs 2.2との連携により、正確な一次ループヒッグス質量計算が可能となり、FeynHiggsFast近似の最高度の精度または完全な一次ループ結果の両方の選択が可能になった。
• プリプロセッサフラグと拡張されたモデルファイルを通じて、MSSMにおけるNMFVは完全な6×6スクォーク混合行列を用いてサポートされた。
• HadCalcフロントエンドにより、PDF畳み込み、運動量カット、およびインバリアント質量・ラピディティ・横運動量に関する微分分布を備えた、自動的かつ高精度なハドロン断面積計算が可能になった（要請に応じて提供可能）。