2025-11-23T21:07:17.387352

Energy dependence of $J/Ï$ production in pp collisions with the PACIAE model

Ye, Wang, Shi et al.

In this work we investigate the $J/Ï$ production in proton-proton collisions at the center-of-mass energy ($\sqrt{s}$) equal to 2.76, 5.02, 7, 8 and 13 TeV with a parton and hadron cascade model PACIAE 2.2a. It is based on PYTHIA but extended considering the partonic and hadronic rescatterings before and after hadronization, respectively. In the PYTHIA sector the $J/Ï$ production quantum chromodynamics processes are selected specially and a bias factor is proposed correspondingly. The calculated total cross sections, the differential cross sections as a function of the transverse momentum and the rapidity of $J/Ï$ in the forward rapidity region reproduce the corresponding experimental measurements reasonably well. In the mid-rapidity region, the double differential cross sections at $\sqrt{s}=$ 5.02, 7 and 13 TeV are also in a good agreement with the experimental data. Moreover, we interpolate the double differential cross section as well as the total cross section of $J/Ï$ in the mid-rapidity region at $\sqrt{s}=$ 8 TeV, which could be validated if the experimental data is available.

academic

Energy dependence of $J/\psi$ production in pp collisions with the PACIAE model

基本信息

论文ID: 2310.12627
标题: Energy dependence of $J/\psi$ production in pp collisions with the PACIAE model
作者: Kai-Fan Ye, Qiang Wang, Jia-Hao Shi, Zhi-Ying Qin, Wen-Chao Zhang, An-Ke Lei, Zhi-Lei She, Yu-Liang Yan, Ben-Hao Sa
分类: hep-ph (高能物理-现象学), nucl-th (核理论)
发表时间: 2024年2月9日
论文链接: https://arxiv.org/abs/2310.12627v2

摘要

本研究使用部分子和强子级联模型PACIAE 2.2a研究了质心能量 $\sqrt{s}$ 为2.76、5.02、7、8和13 TeV的质子-质子对撞中 $J/\psi$ 介子的产生。该模型基于PYTHIA，但扩展考虑了强子化前的部分子散射和强子化后的强子散射。在PYTHIA部分特别选择了 $J/\psi$ 产生的量子色动力学过程，并相应地提出了偏置因子。计算得到的总截面、前向快度区域 $J/\psi$ 的横动量和快度微分截面与相应的实验测量结果符合良好。在中快度区域， $\sqrt{s}=5.02、7和13$ TeV的双微分截面也与实验数据很好吻合。此外，研究还内插了 $\sqrt{s}=8$ TeV中快度区域 $J/\psi$ 的双微分截面和总截面，可在实验数据可用时进行验证。

研究背景与动机

核心问题: $J/\psi$ 介子作为最轻的矢量粲偶素介子，其产生机制的理解对于研究夸克-胶子等离子体(QGP)至关重要。在相对论重离子碰撞中， $J/\psi$ 的抑制被认为是QGP形成的重要探针。
问题重要性: 为了区分热介质效应和冷核物质效应，必须首先理解质子-质子对撞中 $J/\psi$ 的产生机制，这是研究核-核对撞的重要基准。
现有方法局限性:
- 色单态模型、非相对论QCD模型和色蒸发模型等理论方法都无法同时描述 $J/\psi$ 的极化、横动量谱和截面的能量依赖性
- PYTHIA 6.4在描述 $J/\psi$ 产量与带电粒子多重数的关联方面存在不足
- 缺乏考虑部分子和强子散射效应的完整模型
研究动机: 使用PACIAE 2.2a模型系统研究不同能量下pp对撞中 $J/\psi$ 产生的能量依赖性，该模型在PYTHIA基础上增加了部分子和强子散射过程。

核心贡献

模型应用: 首次使用PACIAE 2.2a模型系统研究了多个LHC能量下pp对撞中 $J/\psi$ 产生
偏置因子方法: 提出了处理特殊选择的 $J/\psi$ 产生过程的偏置因子方法
全面对比: 在前向快度区域(2.5 < y < 4)和中快度区域(-0.9 < y < 0.9)与ALICE实验数据进行了系统对比
预测结果: 为 $\sqrt{s}=8$ TeV中快度区域提供了 $J/\psi$ 产生截面的理论预测

部分子初始化阶段: 基于PYTHIA 6.4产生硬散射、初态和末态辐射
部分子散射阶段: 使用领头阶微扰QCD部分子-部分子相互作用截面，引入K因子考虑高阶效应
强子化阶段: 通过弦碎裂或聚合模型将部分子转化为强子
强子散射阶段: 考虑强子间的弹性和非弹性散射过程

$J/\psi$ 产生过程选择:

设置MSEL=0和MSUB(i)=1 (i=86,87,88,89,104,105,106)
对应色单态过程: $gg \to J/\psi g$ , $gg \to \chi_{0c}g$ , $gg \to \chi_{1c}g$ , $gg \to \chi_{2c}g$ , $gg \to \chi_{1c}$ , $gg \to \chi_{2c}$ , $gg \to J/\psi\gamma$

技术创新点

偏置因子方法: $B = \frac{\int_{p_T} \frac{d^2\sigma_{J/\psi}}{dp_T dy}\Big|_{\text{exp}} dp_T}{\int_{p_T} \frac{d^2\sigma_{J/\psi}}{dp_T dy}\Big|_{\text{sim}} dp_T}$
K因子优化: 通过最小化 $\chi^2$ 方法确定最优K因子: $\chi^2 = \sum_{i=1}^{N} \frac{(y_i - f(p_T^i; K))^2}{(\Delta_i)^2}$
散射过程: 考虑了详细的 $J/\psi$ -强子非弹性散射过程，如 $J/\psi + \pi \to D + D^*$ 等

实验设置

数据集

实验数据: ALICE合作组发表的pp对撞 $J/\psi$ 产生数据
能量范围: $\sqrt{s} = 2.76, 5.02, 7, 8, 13$ TeV
快度区间: 前向快度(2.5 < y < 4)和中快度(-0.9 < y < 0.9)
模拟样本: 每个能量点生成 $10^6$ 个事例

评价指标

总截面 $\sigma_{J/\psi}$
横动量微分截面 $d^2\sigma_{J/\psi}/dp_T dy$
快度微分截面 $d\sigma_{J/\psi}/dy$
$\chi^2/\text{ndf}$ 拟合优度

对比方法

与ALICE实验测量结果直接对比，参考了PYTHIA 6.4、PYTHIA 8.2、EPOS3等模型的研究结果

实现细节

K因子范围: 1.1-1.9，步长0.1
最优K因子通过 $\chi^2$ 最小化确定
考虑了 $\chi_c$ 态向 $J/\psi$ 的辐射衰变贡献

实验结果

主要结果

前向快度区域表现:
- 横动量微分截面与实验数据符合良好(图2)
- 快度微分截面准确再现实验测量(图3)
- 总截面随能量增长趋势与实验一致(图5)
中快度区域结果:
- $\sqrt{s} = 5.02, 7, 13$ TeV的双微分截面与实验数据吻合(图7)
- 总截面在误差范围内与实验符合(图8)
K因子优化结果:
- 2.76 TeV: K = 1.6 ( $\chi^2/\text{ndf} = 3.46/6$ )
- 5.02 TeV: K = 1.7 ( $\chi^2/\text{ndf} = 10.92/6$ )
- 7 TeV: K = 1.3 ( $\chi^2/\text{ndf} = 16.62/6$ )
- 8 TeV: K = 1.5 ( $\chi^2/\text{ndf} = 13.40/6$ )
- 13 TeV: K = 1.7 ( $\chi^2/\text{ndf} = 54.86/7$ )