[論文レビュー] Exact solutions of the QCD evolution equations using Monte Carlo method

S. Jadach, M. Skrzypek|ArXiv.org|Dec 26, 2003

High-Energy Particle Collisions Research参考文献 2被引用数 4

ひとこと要約

本稿では、現代のCPU性能を活用して、Markov連鎖モンテカルロ法を用いてQCDの発展方程式を正確かつ高精度に数値的に解く手法を提示している。この手法により、$Q$ と $x$ の広い範囲で約0.1%の精度が達成された。本研究は、現代のCPU性能が、一部の分布関数の正確なモンテカルロベースの発展を現実可能としていることを示し、初期状態放射の完全なQCD発展を組み込んだ非Markov連鎖モンテカルロアルゴリズムの開発を挑戦課題として提示している。この手法は、事前に計算されたPDFテーブルに依存しない。

ABSTRACT

We present the exact and precise (~0.1%) numerical solution of the QCD evolution equations for the parton distributions in a wide range of $Q$ and $x$ using Monte Carlo (MC) method, which relies on the so-called Markovian algorithm. We point out certain advantages of such a method with respect to the existing non-MC methods. We also formulate a challenge of constructing non-Markovian MC algorithm for the evolution equations for the initial state QCD radiation with tagging the type and $x$ of the exiting parton. This seems to be within the reach of the presently available computer CPUs and the sophistication of the MC techniques.

研究の動機と目的

• 現代のCPU性能を活用して、モンテカルロ法を用いてQCD発展方程式を正確かつ高精度に解く実用的妥当性を検証すること。
• Markov連鎖モンテカルロ手法が、広い$Q$および$x$の範囲で発展方程式を約0.1%の精度で解くことができることを示すこと。
• 初期状態放射の完全なQCD発展を組み込んだ非Markov連鎖モンテカルロアルゴリズムの開発を促進する挑戦を提示すること。この際、一部の分布関数タイプと$x$-割合を特定する。
• 今後のパートオンシャワー生成器に正確なモンテカルロ発展を統合することで、LHCおよびDISのデータを同時にフィッティング可能にするよう提言すること。
• 事前に計算されたPDFテーブルの代わりに、動的発展された分布関数をイベントジェネレータに組み込む可能性を検討すること。

提案手法

• QCD発展方程式を$t = \ln Q$における確率過程として解釈することで、Markov連鎖モンテカルロ法を用いて発展方程式を解く。$t$方向へのランダムなステップと一部の分布関数の崩壊を含む。
• 発展方程式の積分形を、正の被積分関数を持つ繰り返し積分の系列に書き直し、適切に正規化された遷移確率を持つMarkov過程として解釈できる。
• 遷移確率$\omega(t_i, K_i | t_{i-1}, K_{i-1})$は、分裂関数と崩壊率を用いて定義され、正規化および確率的整合性を保証する。
• アルゴリズムは、$t$と一部の分布関数状態空間におけるランダムウォークにより、シャワーをシミュレートし、最終的な$Q$スケールで停止ルールを適用する。
• グリッドベースの離散化を避けて、発展経路の確率的サンプリングに依存することで、$D_I(t)$を正確に計算する。
• 初期状態放射の非Markov連鎖アルゴリズムの拡張も行い、高エネルギー$Q$での分布関数タイプと$x$を特定する。ブートストラップまたは最適化されたサンプリングを用いる。

実験結果

リサーチクエスチョン

• RQ1 Markov連鎖モンテカルロ法は、$Q$ と $x$ の広い範囲で、QCD発展方程式を正確かつ高精度（約0.1%）に数値的に解くことができるか？
• RQ2 QCD発展のためのモンテカルロ法と従来の非モンテカルログリッドベースのソルバーとを比較した場合、実用的な利点は何か？
• RQ3 今日の技術的水準では、初期状態放射の完全なQCD発展を組み込み、分布関数タイプと$x$-割合を特定する非Markov連鎖モンテカルロアルゴリズムを構築するのは可能か？
• RQ4 現代のCPU性能とモンテカルロ技術を活用することで、事前に計算されたテーブルに依存せずに、低$Q$から高$Q$へPDFを発展させる非Markov連鎖モンテカルロアルゴリズムを構築できるか？
• RQ5 正確なモンテカルロ発展をパートオンシャワー生成器に統合することで、LHCおよびDISのデータを同時にフィッティングする際にどのような影響が生じるか？

主な発見

• Markov連鎖モンテカルロ法は、$Q$ が1 GeVから1 TeVの範囲で、QCD発展方程式を約0.1%の精度で解くことに成功した。
• $Q = 1$ GeVから$Q = 1$ TeVまでのプロトンシャワーにおける平均的な分布関数発生回数は約20回であり、顕著な放射が発生していることを示している。
• $q,\bar{q} \to G$および$G \to q,\bar{q}$の遷移の平均回数は低く、0.5から1の範囲にとどまり、これらの過程は事前に表形式で保存できる可能性がある。
• 本研究では、標準的なMarkov連鎖モンテカルロ法が初期状態放射に不適切であるという低効率性が、過去15年間でCPU性能が10,000倍に向上したことで、もはや根本的な障壁ではないと示している。
• 著者らは、内蔵された発展機能を持つ非Markov連鎖モンテカルロアルゴリズムが今や実現可能であり、より正確で自己整合性のあるイベントジェネレータの実現が可能になると主張している。
• 著者らは、$Q_{\min}$から$Q_{\max}$までNLL $\overline{\text{MS}}$カーネルを用いてPDFを発展させるモンテカルロアルゴリズムの開発を挑戦課題として提示している。この際、$Q_{\max}$で固定された分布関数タイプと$x$を保持し、事前に計算されたPDFテーブルに依存しない。