[论文解读] The PHOBOS Glauber Monte Carlo

B. Alver, M. Baker+2 more|ArXiv.org|May 28, 2008

High-Energy Particle Collisions Research被引用 136

一句话总结

本文提出了一种用于模拟重离子碰撞初始态几何结构的PHOBOS蒙特卡罗Glauber模型，通过随机核子定位与碰撞参数采样，计算$N_{\text{part}}$、$N_{\text{coll}}$和椭圆度。其主要贡献在于提供了一个公开可用、用户友好的实现版本，包含详细的逐事件可观测量，支持对RHIC与LHC能量下相对论性重离子物理中由涨落驱动效应的精确研究。

ABSTRACT

``Glauber'' models are used to calculate geometric quantities in the initial state of heavy ion collisions, such as impact parameter, number of participating nucleons and initial eccentricity. The four RHIC experiments have different methods for Glauber Model calculations, leading to similar results for various geometric observables. In this document, we describe an implementation of the Monte Carlo based Glauber Model calculation used by the PHOBOS experiment. The assumptions that go in the calculation are described. A user's guide is provided for running various calculations.

研究动机与目标

• 提供重离子碰撞初始态几何结构蒙特卡罗Glauber模型的详细、可复现实现。
• 支持对$N_{\text{part}}$、$N_{\text{coll}}$和椭圆度（$\epsilon_{\text{part}}$、$\epsilon_{\text{RP}}$）等逐事件几何可观测量的精确计算，以服务于PHOBOS分析。
• 提供一个用户友好、可扩展的代码库，并附带教程与文档，适用于实验与理论应用。
• 通过逐事件采样参与核子构型，支持对初始态涨落的研究。
• 对PHOBOS特定的Glauber模型实现进行标准化与文档化，确保分析间的一致性，并促进与其他实验的合作。

提出的方法

• 利用电子散射数据导出的径向概率分布，随机定位核子在原子核内的位置，电荷密度采用费米或高斯形式。
• 引入最小核子间距（$d_{\text{min}}$）以防止物理上不合理的聚集，缺省值为0.4 fm。
• 从$b$加权分布中随机采样碰撞参数$b$，上限为$b_{\text{max}} \approx 20$ fm，两核位于$(-b/2, 0, 0)$与$(b/2, 0, 0)$。
• 通过横向距离判断二元碰撞：若两核子间距小于“球直径”$D = \sqrt{\sigma_{\text{NN}} / \pi}$，则视为发生碰撞。
• 采用eikonal近似：核子沿束流方向直线运动，不考虑纵向动力学。
• 逐事件计算几何可观测量，包括$N_{\text{part}}$、$N_{\text{coll}}$以及参与核子位置的矩，用于椭圆度计算。

实验结果

研究问题

• RQ1 PHOBOS蒙特卡罗Glauber模型如何在重离子碰撞中逐事件生成初始态几何结构？
• RQ2 核子位置的随机性与最小间距对$N_{\text{part}}$与$N_{\text{coll}}$分布有何影响？
• RQ3 参与核子椭圆度$\epsilon_{\text{part}}$与反应平面椭圆度$\epsilon_{\text{RP}}$如何随$N_{\text{part}}$在不同碰撞体系中变化？
• RQ4 在RHIC与LHC能量下，Au+Au、Cu+Cu与Pb+Pb碰撞中，几何结构的逐事件涨落对初始态的影响程度如何？
• RQ5 标准化的、公开可用的Glauber模型实现如何支持重离子实验间的一致性分析？

主要发现

• 该模型成功再现了RHIC与LHC能量下Cu+Cu、Au+Au与Pb+Pb碰撞中$N_{\text{part}}$与$N_{\text{coll}}$的预期分布，如图2所示，基于10,000个事件。
• 在低$N_{\text{part}}$区域，参与核子椭圆度$\epsilon_{\text{part}}$显著大于$\epsilon_{\text{RP}}$，表明初始态存在强烈的非高斯涨落。
• 在所有体系中，$\epsilon_{\text{part}}$随$N_{\text{part}}$增加而减小，反映了中央碰撞中几何结构的平滑化。
• 该代码实现可通过公式$\epsilon_{\text{part}} = \sqrt{(\text{VarY}-\text{VarX})^2 + 4\text{VarXY}^2}/(\text{VarY}+\text{VarX})$精确计算$\epsilon_{\text{part}}$，事件级矩存储于输出ntuple中。
• 总核-核截面在Run()函数内部计算，并可通过GetTotXSect()访问，各事件样本间结果一致。
• 通过与预期分布对比验证了模型的正确性，且代码已公开，可用于独立或嵌入式分析工作流。