[论文解读] MBR Monte Carlo Simulation in PYTHIA8

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

Particle physics theoretical and experimental studies被引用 10

一句话总结

本论文介绍了在 PYTHIA8 中实现的 mbr（最小偏差 Rockefeller）蒙特卡洛模拟，这是一种基于重整化 Regge 理论的唯象模型，可预测高达 50 TeV 的质子-质子总截面、弹性截面和非弹性截面。该模型首次在 PYTHIA8 中实现了双 Pomeron 交换（中心衍射）过程，完整模拟了单侧、双侧和中心衍射 dissociation，经 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.

研究动机与目标

• 开发一种基于唯象 Regge 理论模型的完整蒙特卡洛模拟，用于 LHC 能量范围内的最小偏差质子-质子相互作用。
• 在 PYTHIA8 框架内实现并验证三种衍射过程——单侧、双侧和中心衍射的模型。
• 扩展 PYTHIA8，首次实现对双 Pomeron 交换（中心衍射）过程的完整模拟。
• 基于 CDF 数据约束，提供高达 50 TeV 的总截面、弹性截面和非弹性截面的能量依赖性预测。
• 使研究人员能够使用标准化的 PYTHIA8 控制卡片，模拟具有精确衍射动力学的真实最小偏差事件。

提出的方法

• 该模型采用重整化 Regge 理论框架，其中 Pomeron 轨迹为 α(t) = 1.104 + 0.25 GeV⁻²·t，且通过通量重整化因子 N_gap(s) 来描述快度间隙分布。
• 衍射截面通过涉及 β(t)（Pomeron-质子耦合）、κ（三 Pomeron 耦合比）以及基于 s₀ 和 s_F 参数的能量依赖项的微分方程计算得出。
• 通过采样 e^{2α′Δyt} 形式的 t 分布生成事件，并对质子动量转移和强子系统四momentum 施加运动学约束。
• 对于中心衍射（DPE），使用高斯抑制函数生成快度间隙，并对中心系统的快度进行均匀采样，通过 ξ₁ 和 ξ₂ 计算 t₁ 和 t₂。
• 在 PYTHIA8 中通过 'Diffraction:PomFlux = 5' 激活该实现，设有专门的标志用于中心衍射（SoftQCD:centralDiffractive = on）和过程特异性参数。
• 事件生成包括对非 pp 束流的屏蔽检查，并使用预定义的默认参数，如 β(0) = 6.566 GeV⁻¹、σ₀ = 2.82 mb 和 m²₀ = 1.5 GeV² 作为最小质量阈值。

实验结果

研究问题

• RQ1 在高达 50 TeV 的质心系能量范围内，总截面、弹性截面和非弹性截面如何随能量变化？
• RQ2 在 LHC 能量范围内，单侧、双侧和中心衍射对总非弹性截面的相对贡献是什么？
• RQ3 基于重整化 Regge 理论的模型在 1.8 TeV 能量下，对 CDF 测量的总截面和弹性截面的重现精度如何？
• RQ4 包含双 Pomeron 交换（中心衍射）对 PYTHIA8 中最小偏差事件生成有何影响？
• RQ5 快度间隙分布和运动学约束如何影响高能 pp 碰撞中衍射过程的模拟？

主要发现

• 在 √s = 1.8 TeV 时，mbr 模型预测总截面为 80.03 ± 2.24 mb，与 CDF 测量值一致，作为归一化参考点。
• 弹性截面与总截面之比由 LHC 前数据的全局 Regge 拟合得出，σ_el/σ_tot 从约 0.2（0.3 TeV 时）下降至约 0.1（14 TeV 时）。
• 在 √s = 14 TeV 时，模型预测 σ_inel ≈ 63.2 mb，其中 SD（≈15.5 mb）、DD（≈10.8 mb）和 CD（≈10.5 mb）贡献，非衍射分量约为 26.4 mb。
• 首次在 PYTHIA8 中模拟了双 Pomeron 交换（DPE）过程，使用专用过程 ID（106）和使用 rho_diff0 表示中心系统的事件记录结构。
• 通过使用 σ_S = 0.5 的高斯形状抑制函数，在低快度间隙（Δy < 2.0）区域引入抑制因子，减少非衍射污染。
• 该模型成功再现了高达 50 TeV 的截面能量依赖性，σ_tot 在 √s < 1.8 TeV 时随 s^0.104 增长，而在更高能量下趋于 Froissart 界限。