[論文レビュー] New Developments in WHIZARD Version 2.6

Wolfgang Kilian, Brass, S.+5 more|arXiv (Cornell University)|Jan 24, 2018

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

ひとこと要約

WHIZARD 2.6 は、$W$ および $Z$ レゾナントピークの形状を保つために、レゾナントとコンtinuumのパートンシャーリングの滑らかな遷移を実装するなど、レプトン衝突機におけるイベント生成に大きな改善をもたらします。MPI および OpenMP を用いた並列化により、位相空間統合が高速化され、半硬化初期状態放射（ISR）の取り扱いも洗練されています。これらの進歩により、将来の $e^+e^-$ 衝突機、例えば ILC や CLIC における標準模型過程のより正確で、効率的かつスケーラブルなシミュレーションが可能になります。

ABSTRACT

We describe recent additions to the WHIZARD 2 Monte-Carlo event generator which improve the physics description of lepton-collider event samples and speed up the calculation time required for cross sections and event generation.

研究の動機と目的

• 将来のレプトン衝突機における $e^+e^-$ 衝突のシミュレーションにおいて、PYTHIA 6 のパートンシャーの近似によって生じる運動量分布の歪みを是正すること。
• パートンシャー入力におけるレゾナントとコンtinuum仮説の間の滑らかな遷移を実装することで、より正確で効率的なイベント生成を実現すること。
• VAMP 位相空間統合モジュールにおける MPI および OpenMP による完全な並列化を導入することで、断面積計算およびイベント生成の計算時間を大幅に短縮すること。
• ISR とレゾナント崩壊チェーンの運動量マッチングを洗練させることで、半硬化初期状態放射（ISR）の取り扱いを改善すること。
• ILC や CLIC の物理的・検出器研究における高精度シミュレーションの互換性とパフォーマンス向上を確保すること。

提案手法

• レゾナントとコンtinuum仮説の間の滑らかな遷移を実現するために、ガウス関数による重み付けを施した因子分解行列要素を用い、排他的分布における急激な運動量不連続性を解消します。
• VAMP の適応的モンテカルロ統合エンジンを再実装し、MPI を完全にサポートすることで、複数ノードにまたがる分散並列化による位相空間点評価を可能にしました。
• 各計算ノードにおけるマルチコアプロセッサを活用するため、OpenMP 指令を用いて、1つのノード内でスピン状態の並列評価を実現しました。
• 新しい VAMP 実装では、ノード間のグリッド交換中にブロッキングを最小限に抑えるために、非同期 MPI 通信を採用し、通信オーバーヘッドを低減しました。
• レゾナント歴史の選択は、因子分解行列要素と全行列要素の比に基づく確率に従い、レゾナントピークに近いインバリアント質量に近づくようにガイドされます。
• O’Mega と統合することで、オンザフライでの行列要素生成が可能となり、PYTHIA 6 や WHIZARD の解析的シャワー、外部イベントファイル形式とも完全に互換性を保ちます。

実験結果

リサーチクエスチョン

• RQ1 どのようにしてパートンシャーリングを改善すれば、$e^+e^-$ 衝突過程における $W$ および $Z$ レゾナントピークの形状を人工的な不連続性を導入せずに保てるか？
• RQ2 パートンシャワーにおけるレゾナントとコンtinuum放射の間の最適な遷移戦略は何か？ これは行列要素計算と整合性を保つために必要である。
• RQ3 MPI および OpenMP を用いた VAMP 位相空間統合エンジンの並列化は、どのようにして複雑な多粒子最終状態のイベント生成を高速化できるか？
• RQ4 半硬化 ISR 光子の取り扱いが、$e^+e^-$ 衝突機シミュレーションにおけるレゾナントピーク再構築の正確さにどのように影響するか？
• RQ5 新しい並列化 VAMP 実装は、さまざまな高多重度プロセスにおいて、どの程度のパフォーマンス向上とスケーラビリティの限界を達成できるか？

主な発見

• ガウス関数による再重み付けを用いた、レゾナントとコンtinuum仮説の滑らかな遷移により、排他的分布における人工的な段階が完全に解消され、NLO/NNLO マッチングを必要とせず、真のレゾナント構造が保持されました。
• MPI および OpenMP を完全にサポートする新しい VAMP 実装により、スケーラブルな並列化が実現され、プロセスに応じて最大 $O(10\dots 100)$ 個のプロセッサでスピードアップが達成され、標準的なマルチコア CPU においても強いスケーリングが確認されました。
• O’Mega からのプロセス構造を活用した位相空間パrameter化構築アルゴリズムの再実装により、従来のシリアルバッテリーブレークが解消され、将来のバージョンでさらなるパフォーマンス向上が見込まれます。
• $e^+e^- \to q\bar{q}q\bar{q}$ プロセスにおいて、新しい手法は $W$ と $Z$ レゾナントピーク（それぞれ 80 および 91 GeV）を正確に保持するのに対し、レゾナント歴史を考慮しない標準的な PYTHIA 6 シャワーでは、これらのピークは完全に消えてしまいました。
• 半硬化 ISR の改善により、放射の運動量がレゾナント崩壊と正確に一致し、イベントサンプルにおける非物理的な歪みが低減されました。
• 新しい WHIZARD 2.6 フレームワークは、レゾナントおよび非レゾナント寄与を含む完全な行列要素を用いた完全排他的なイベント生成を、LO およびそれ以上の高次の精度で可能にし、ILC や CLIC の高精度研究に適しています。