[论文解读] Fast simulation of flow effects in central and semi-central heavy ion collisions at LHC

I. Lokhtin, Alexander Snigirev|ArXiv.org|Dec 16, 2003

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

一句话总结

本论文提出了一种用于模拟LHC能量下中心及半中心重离子碰撞中集体流效应（纵向、横向及椭圆流）的快速蒙特卡洛事件生成器。该方法采用简化的流体动力学模型，结合热辐射与通过快速度依赖的洛伦兹提升变换实现的集体流，无需复杂调参即可获得真实的粒子谱和流各向异性，关键输出结果包括 $\langle p_T \rangle = 0.55$ GeV/c，以及与实验趋势相符的粒子比值。

ABSTRACT

The simple method for simulation of ``thermal'' hadron spectra in ultrarelativistic heavy ion collisions including longitudinal, transverse and elliptic flow is developed. The model is realized as fast Monte-Carlo event generator.

研究动机与目标

• 开发一种计算高效的重离子碰撞集体流模拟方法，适用于LHC能量范围。
• 再现中心及半中心Pb-Pb碰撞中实验观测到的横向与椭圆流模式。
• 为LHC重离子物理提供一种快速、用户友好的事件生成器，避免全微观模型带来的高计算成本。
• 通过一种基于流体动力学的快速模拟工具，实现LHC探测器软件与物理模型的验证和交叉比较。

提出的方法

• 模型在局部静止系中使用具有玻色-爱因斯坦或费米-狄拉克统计的相对论性玻尔兹曼分布生成强子四动量。
• 通过由横向与纵向快度定义的空间可变四维速度场实施洛伦兹提升变换，其中 $u^r = (r/R_f) \sinh Y_T^{\text{max}}$ 且 $u^z = \sinh \eta$。
• 通过基于初始重叠偏心度 $\epsilon_0 = b/(2R_A)$ 推导出的角度依赖有效冻结-out半径 $R_f(b, \Phi)$ 引入椭圆流。
• 模型使用泊松多重度分布，并按核重叠函数缩放以确定粒子产额。
• 最终粒子谱通过在 $T_f = 140$ MeV 的冻结-out超曲面上积分获得，其中 $Y_L^{\text{max}} = 5$ 且 $Y_T^{\text{max}} = 1$。
• 该生成器以Fortran实现，并通过JETSET的LUJETS公共块接口输出事件，支持最多150,000个粒子。

实验结果

研究问题

• RQ1 简化的流体动力学模型能否准确再现LHC重离子碰撞中的横向与椭圆流模式？
• RQ2 如何在不使用完整微观事件生成器的前提下高效模拟集体流效应？
• RQ3 在LHC能量尺度下，热模型结合集体流会产生怎样的粒子谱与比值？
• RQ4 该模型在稀疏碰撞及高 $p_T$ 或前向快度区域的适用性如何，其有效性在多大程度上仍成立？

主要发现

• 在 $T_f = 140$ MeV、$Y_T^{\text{max}} = 1$ 且 $Y_L^{\text{max}} = 5$ 条件下，模型生成的平均横向动量为 $\langle p_T^h \rangle = 0.55$ GeV/c。
• 粒子比值 $\pi^\pm : K^\pm : p^\pm = 24:6:1$ 与中快速度区域的实验趋势相符良好。
• 模型准确再现了 $\pi^\pm : \pi^0 = 2:1$、$K^\pm : K^0 = 1:1$ 以及 $p:n = 1:1$ 的比值，与统计模型预期一致。
• 粒子的方位各向异性由椭圆流形式良好描述，$v_2$ 随 $p_T$ 增加至 $\sim 2$ GeV/c，与RHIC观测结果一致。
• 该生成器计算效率高，适用于中心及半中心碰撞，且无需对内部参数进行额外调校。
• 该模型不适用于极稀疏碰撞（$b \sim 2R_A$）、极高前向快度区域（$|y| \gtrsim 3$）或高 $p_T$ 区域（$\gg 2-5$ GeV/c），因为在这些区域流体动力学假设不再成立。