[論文レビュー] Solving QCD evolution equations in rapidity space with Markovian Monte Carlo

K. Golec-Biernat, S. Jadach+2 more|arXiv (Cornell University)|Aug 14, 2007

Particle physics theoretical and experimental studies参考文献 2被引用数 5

ひとこと要約

本稿では、急速性空間におけるDGLAP、CCFM、および関連モデルのQCD進化方程式を解くためのマルコフ連鎖モンテカルロ（MMC）法を提示する。進化変数をパートンの4次元運動量に写像することで、約0.1%の精度を達成した。この手法は非モンテカルロ手法であるAPChebプログラムと比較され、パートンシャワー・モンテカルロ開発における高い正確性と頑健性を示した。

ABSTRACT

This work covers methodology of solving QCD evolution equation of the parton distribution using Markovian Monte Carlo (MMC) algorithms in a class of models ranging from DGLAP to CCFM. One of the purposes of the above MMCs is to test the other more sophisticated Monte Carlo programs, the so-called Constrained Monte Carlo (CMC) programs, which will be used as a building block in the parton shower MC. This is why the mapping of the evolution variables (eikonal variable and evolution time) into four-momenta is also defined and tested. The evolution time is identified with the rapidity variable of the emitted parton. The presented MMCs are tested independently, with ~0.1% precision, against the non-MC program APCheb especially devised for this purpose.

研究の動機と目的

• DGLAP、CCFM、および関連モデルを含む多様なモデルにおけるQCD進化方程式を解くためのマルコフ連鎖モンテカルロ（MMC）フレームワークの開発。
• 統計誤差が0.1%未満となるよう、部分的分布関数（PDF）の正確で安定した数値解を得ること。
• 進化変数（エイコナールおよび急速性）を物理的4次元運動量に写像する手法の確立、これによりパートンシャワー・モンテカルロシミュレーションへの直接応用を可能にする。
• MMC手法を非MC手法であるAPChebプログラムと照合し、ベンチマークとしての妥当性を検証し、将来の制約付きモンテカルロ（CMC）フレームワークへの応用を想定する。
• 急速性空間における進化のための堅牢でスケーラブルな基盤を提供することで、より現実的なマルチパートンシャワーシミュレーションの開発を支援すること。

提案手法

• 本手法は、一般化されたQCD進化方程式 ∂ₜD_f(t,x) = Σ_f' ∫_x¹ du K_ff'(t,x,u) D_f'(t,u) を、マルコフ連鎖モンテカルロによるサンプリングで解く。
• 進化時間tは、放出されたパートンの急速性と同一視され、物理的4次元運動量への直接的写像が可能になる。
• MMCアルゴリズムは、放出歴の確率的サンプリングを用い、DGLAP、CCFM、および拡張モデルに対応した進化カーネルを適応的に使用する。
• 高精度での相違検証のため、非モンテカルロ手法としてのチェビシェフ多項式コロケーション法（APCheb）が用いられる。
• 進化カーネルの対角化を必要とせず、グルーオン、クォーク、反クォークの全パートン種を同時に取り扱える。
• 大規模なMC計算（約10¹⁰イベント）により統計的収束を監視・制御し、APChebとの一致が0.2%未満に保たれることを確認した。

実験結果

リサーチクエスチョン

• RQ1 マルコフ連鎖モンテカルロ技法は、CCFMやDGLAPを含む多様なモデルにおいて、QCD進化方程式の高精度な数値解を得られるか？
• RQ2 MMCは、特に小x領域において、非モンテカルロ手法（APCheb）の結果をどの程度正確に再現できるか？
• RQ3 進化変数（エイコナールおよび急速性）を、パートンシャワー・モンテカルロシミュレーションへの応用を想定して、物理的4次元運動量に一貫して写像できるか？
• RQ4 複雑な進化カーネルを扱う際、従来のグリッドベース手法に比べ、MMCの安定性、スケーラビリティ、適用可能性にどのような利点があるか？
• RQ5 MMCは、より高度な制約付きモンテカルロ（CMC）アルゴリズムのテストフレームワークとして、信頼できる役割を果たせるか？

主な発見

• MMC手法は、非モンテカルロ手法であるAPChebプログラムと約0.1%の精度で一致し、数値的頑健性が裏付けられた。
• α((1−z)Q)を伴う進化タイプ（B）において、グルーオンおよびクォークPDF（xD_f(x)）について、Q = 1, 10, 100, 1000 GeVの全エネルギー領域で、MMCとAPChebの結果が統計誤差約0.2%の範囲内で一致した。
• α(kᵀ)を伴う進化タイプ（C）においても、MMCとAPChebの一致は約0.2%のレベルを維持し、異なる進化カーネル間の整合性が確認された。
• 本手法は、先行研究と比較して統計誤差が小さく、数値的バイアスも低減したLL DGLAP結果を正確に再現した。
• 進化変数から4次元運動量への写像が成功裏に定義され、テストを経て、パートンシャワー・モンテカルロフレームワークへの直接的応用が可能となった。
• MMCフレームワークは、QCD進化シミュレーションにおけるより複雑な制約付きモンテカルロ（CMC）アルゴリズムのテストに、信頼性があり、高精度なツールとして有効であることが示された。