[論文レビュー] MBR Monte Carlo Simulation in PYTHIA8

R. Ciesielski, K. Goulianos|arXiv (Cornell University)|May 7, 2012

Particle physics theoretical and experimental studies被引用数 10

ひとこと要約

この論文は、50 TeVまでにまで全エネルギー領域で全断面積、弾性断面積、非弾性断面積を予測する、正規化されたレッジ理論に基づく現象論的モデルであるmbr（最小バイアス・ロックフェラー）モンテカルロシミュレーションをPYTHIA8に実装した。中央フラグメンテーション（二重ポメロン交換）プロセスのPYTHIA8における初の実装を含み、単一・二重・中央フラグメンテーション分解の完全なシミュレーションが行われ、CDFデータと照合され、'Diffraction:PomFlux = 5'フラグにより有効化されている。

ABSTRACT

We present the MBR (Minimum Bias Rockefeller) Monte Carlo simulation of (anti)proton-proton interactions and its implementation in the PYTHIA8 event generator. We discuss the total, elastic, and total-inelastic cross sections, and three contributions from diffraction dissociation processes that contribute to the latter: single diffraction, double diffraction, and central diffraction or double-Pomeron exchange. The event generation follows a renormalized-Regge-theory model, successfully tested using CDF data. Based on the MBR-enhanced PYTHIA8 simulation, we present cross-section predictions for the LHC and beyond, up to collision energies of 50 TeV.

研究の動機と目的

• LHCエネルギーダイナミクス範囲における最小バイアスpp相互作用の包括的モンテカルロシミュレーションを開発すること。
• PYTHIA8フレームワーク内に、単一・二重・中央フラグメンテーションの3つのフラグメンテーションプロセスのモデルを実装および検証すること。
• PYTHIA8に、二重ポメロン交換（中央フラグメンテーション）プロセスの完全なシミュレーションを初導入すること。
• CDFデータ制約に基づき、50 TeVまでに及ぶエネルギー依存の全断面積・弾性断面積・非弾性断面積の予測を提供すること。
• 標準化されたPYTHIA8ステアリングカードを用いて、正確なフラグメンテーションダイナミクスを有する現実的な最小バイアスイベントをシミュレートできるようにすること。

提案手法

• モデルは、ポメロン軌道 α(t) = 1.104 + 0.25 GeV⁻²·t と、急速性ギャップ分布を補正するためのフラックス正規化係数 N_gap(s) を用いた正規化された-Regge理論フレームワークに依拠している。
• フラグメンテーション断面積は、β(t)（ポメロン-陽子結合係数）、κ（三重ポメロン結合係数比）、および s₀ と s_F パラメータに基づくエネルギー依存項を含む微分方程式を用いて計算される。
• イベント生成は、t分布を e^{2α′Δyt} を用いてサンプリングし、陽子の運動量移動およびハドロン系の4元運動量に物理的制約を課すことによって行われる。
• 中央フラグメンテーション（DPE）では、急速性ギャップがガウス型抑制関数を用いて生成され、中央系の急速性が一様にサンプリングされ、t₁ と t₂ が ξ₁ と ξ₂ から計算される。
• PYTHIA8における実装は 'Diffraction:PomFlux = 5' により有効化され、中央フラグメンテーション（SoftQCD:centralDiffractive = on）およびプロセス固有パラメータの専用フラグが用意されている。
• イベント生成には、非-ppビームのバイトチェックが含まれており、β(0) = 6.566 GeV⁻¹、σ₀ = 2.82 mb、m²₀ = 1.5 GeV² といった事前定義済みのデフォルトパラメータが使用されている。

実験結果

リサーチクエスチョン

• RQ1 全断面積・弾性断面積・非弾性断面積のpp断面積は、50 TeVまでに中心系エネルギーの関数としてどのようにスケーリングされるか？
• RQ2 LHCエネルギーダイナミクス範囲において、単一・二重・中央フラグメンテーションが全非弾性断面積に与える相対的寄与はどの程度か？
• RQ3 正規化された-Regge理論モデルは、1.8 TeVにおけるCDFデータの全断面積および弾性断面積をどの程度正確に再現できるか？
• RQ4 二重ポメロン交換（中央フラグメンテーション）を含めることで、PYTHIA8における最小バイアスイベント生成にどのような影響が生じるか？
• RQ5 急速性ギャップ分布および運動的制約は、高エネルギーpp衝突におけるフラグメンテーションプロセスのシミュレーションにどのような影響を及えるか？

主な発見

• mbrモデルは、√s = 1.8 TeV における全断面積を 80.03 ± 2.24 mb と予測し、CDF測定値と一致し、正規化基準点として機能する。
• 弾性断面積/全断面積比は、LHC以前のデータへのグローバルな-Reggeフィットに基づき、0.3 TeVで約0.2、14 TeVで約0.1に減少する。
• √s = 14 TeV における非弾性断面積は σ_inel ≈ 63.2 mb と予測され、SD（≈15.5 mb）、DD（≈10.8 mb）、CD（≈10.5 mb）の寄与に加え、非フラグメンテーション成分が ≈26.4 mb である。
• 二重ポメロン交換（DPE）プロセスは、PYTHIA8で初めて完全にシミュレートされ、専用のプロセスID（106）と、中央系に rho_diff0 を使用するイベントレコード構造が採用された。
• 低急速性ギャップ（Δy < 2.0）における非フラグメンテーション汚染を低減するために、σ_S = 0.5 のガウス型抑制関数が導入されている。
• モデルは、50 TeVまでのエネルギー依存性を正確に再現できており、√s < 1.8 TeV では s^0.104 のように増加し、高エネルギー領域ではFroissart限界に近づく。