Inclusive $J/ψ$ productions in pp collisions at $\sqrt{s}=$ 5.02, 7, and 13 TeV with the PACIAE model

本文利用包含部分子与强子再散射机制的 PACIAE 4.0 模型，结合 NRQCD 框架下的色单态与色八重态贡献、簇塌缩及$b$强子弱衰变，系统研究了 5.02、7 和 13 TeV 下质子 - 质子碰撞中$J/\psi$的包容产生，其模拟结果与实验数据吻合良好，并定量分析了不同产生机制及其对能量和快度的依赖关系。

要点

这篇论文就像是在LHC（大型强子对撞机）的“粒子游乐场”里，给科学家们拍的一部高清纪录片。

简单来说，科学家们想搞清楚：当两个质子（原子的核心）以接近光速猛烈相撞时，是如何产生一种叫 $J/\psi$（J/psi） 的微小粒子的？

为了让你更容易理解，我们可以把整个过程想象成一场**“粒子乐高”的搭建与破坏游戏**。

1. 核心任务：寻找“乐高城堡” ($J/\psi$)

• 什么是 $J/\psi$？ 想象它是一个由两个特殊积木（一个“魅”夸克和一个“反魅”夸克）紧紧抱在一起组成的小城堡。它是物理学中非常重要的“探针”，就像侦探一样，能帮我们看清在极端高温高压下（比如大爆炸初期或重离子碰撞中）物质变成了什么样。
• 为什么要研究它？ 以前科学家发现，在两个大原子核（像两个大乐高盒子）相撞时，这个“小城堡”会变少（被抑制）。但这可能是因为“大盒子”本身太硬了（冷核物质效应），也可能是因为里面太热了（夸克 - 胶子等离子体）。为了分清这两种情况，科学家必须先搞清楚：如果只有两个小质子（小乐高块）相撞，这个“小城堡”到底是怎么造出来的？ 这就是这篇论文要做的基础工作。

2. 他们的工具：PACIAE 4.0 模型（超级模拟器）

科学家没有真的去造无数个粒子对撞机，而是用电脑写了一个超级复杂的**“粒子模拟器”**，叫 PACIAE 4.0。

• 它是怎么工作的？ 这个模拟器基于一个著名的程序（PYTHIA 8.3），但作者给它加了两个“超能力”：
  1. 碰撞前的“推挤” (部分子再散射 PRS)： 在积木正式拼合前，里面的小零件（夸克和胶子）会互相推挤、碰撞。
  2. 拼好后的“打闹” (强子再散射 HRS)： 当积木拼成“城堡”后，周围的其他积木（其他粒子）可能会撞上来，把城堡撞散或者改变它的方向。

3. 城堡是怎么造出来的？（三种主要途径）

在模拟器里，这个“小城堡” ($J/\psi$) 有三种主要的出生方式：

1. 正规军 (NRQCD 机制)： 这是最主要的来源。就像两个工匠（胶子）直接按照精密图纸，把两个特殊积木拼成城堡。这包括“直接拼好”和“先拼个半成品（如 $\chi_c$ 或 $\psi(2S)$），再经过几步变形变成城堡”。
   • 比喻： 就像工厂流水线，按标准流程生产。
2. 意外惊喜 (Cluster Collapse)： 有时候，两个特殊积木在混乱中靠得太近，还没来得及找图纸，就直接“粘”在一起变成了城堡。
   • 比喻： 就像两个迷路的人突然撞在一起，手拉手决定不走了，直接组了个队。
3. 外来客 (非 Prompt 过程)： 这个城堡不是直接造出来的，而是由一个更重的“大怪兽”（含有底夸克的粒子）衰变（解体）后留下的残骸。
   • 比喻： 就像一只大鸟飞走了，留下了一颗蛋，这颗蛋孵化出了城堡。

4. 他们发现了什么？（模拟结果 vs 现实数据）

• 模拟很准： 科学家在电脑上模拟了 5.02、7 和 13 TeV（能量单位，相当于把粒子加速到极快）的碰撞。结果发现，模拟出来的“城堡”数量、速度分布，和现实中 ALICE 等实验组测到的数据非常吻合。这说明他们的“模拟器”升级得很成功！

• 谁贡献最大？

  • 正规军 (NRQCD) 是绝对的主力，贡献了大部分（约 75%-86%）。
  • 意外惊喜 (Cluster Collapse) 和 外来客 (非 Prompt) 是配角，但随着能量越高，“外来客”的比例在增加（因为高能更容易产生那些重的大怪兽）。

• 位置的影响 (快度)：

  • 在正中间（中快度）和两边（前快度），城堡的出生方式略有不同。
  • 在两边，由于那里充满了大量“小能量”的零件，正规军和意外惊喜更容易出现；而外来客（需要高能量怪兽）就变少了。

• 积木的拆解 (细分来源)：

  • 在正规军里，直接拼好的城堡和从 $\chi_{c2}$ 变形来的城堡，在低速时最多。
  • 但在高速时，从 $\chi_{c1}$ 变形来的城堡变多了。这就像高速公路上，某种特定型号的变形车跑得更快。

5. 最有趣的发现：碰撞后的“打闹” (再散射效应)

这是论文的一个亮点。科学家分别关掉了模拟器里的“碰撞前推挤”和“碰撞后打闹”功能，看看会发生什么：

• 碰撞前的推挤 (PRS)： 对城堡的数量几乎没有影响。就像在工地开工前，工人们互相推了一下，并不影响最后盖房子的数量。
• 碰撞后的打闹 (HRS)： 这是一个破坏者！当城堡拼好后，周围乱飞的积木会撞它。
  • 结果：大约 8% 的城堡被撞散了（被抑制了）。
  • 但是！那些“外来客”（由大怪兽衰变来的）因为是在很远的地方才出现的，没受到这次打闹的影响，所以数量没变。

总结

这篇论文就像给物理学家提供了一本**《质子对撞造城堡完全指南》**。

它告诉我们：

1. 我们的模拟工具（PACIAE 4.0）非常靠谱，能精准预测不同能量下的结果。
2. 绝大多数城堡是靠“正规流程”造出来的，但能量越高，“意外”和“外来客”的比例会上升。
3. 最关键的发现是：城堡拼好后，周围环境的“打闹”会毁掉大约 8% 的城堡。

这项研究非常重要，因为它帮我们理清了“背景噪音”。只有彻底搞懂了两个质子相撞时城堡是怎么造出来的，未来当两个大原子核相撞时，如果我们发现城堡变少了，才能确定：“看！这肯定是因为里面产生了高温的夸克 - 胶子等离子体，把城堡融化了！” 而不是因为别的什么原因。

技术摘要

这是一份关于论文《Inclusive J/ψ productions in pp collisions at √s = 5.02, 7, and 13 TeV with the PACIAE model》（使用 PACIAE 模型研究 5.02、7 和 13 TeV 质子 - 质子碰撞中的包含性 J/ψ 产生）的详细技术总结。

1. 研究背景与问题 (Problem)

• 物理背景：$J/\psi$ 介子是重离子物理中的关键探针。其在核 - 核碰撞中的压低最初被视为夸克 - 胶子等离子体（QGP）形成的信号，但也可能源于冷核物质效应（如核部分子分布函数的修正）。
• 核心问题：为了区分热核物质（QGP）和冷核物质效应，必须精确掌握无核初始态（即质子 - 质子，pp 碰撞）下的 $J/\psi$ 产生机制。
• 现有局限：
  • 先前的研究（如作者之前的工作 [22]）主要基于 PACIAE 2.2a 和 PYTHIA 6.4，仅考虑了色单态（color-singlet）过程，忽略了色八重态（color-octet）机制。
  • 未充分包含来自 $b$-强子弱衰变的非 prompt 贡献，以及来自 $\psi(2S)$ 和 $\chi_c$ 等激发态的 feed-down 贡献。
  • 缺乏对部分子级和强子级再散射（rescattering）对 $J/\psi$ 产额影响的定量评估。

2. 方法论 (Methodology)

本研究采用了 PACIAE 4.0 模型，该模型基于 PYTHIA 8.3 进行了扩展，主要特点包括：

• 模型架构：模拟过程分为四个阶段：
  1. 部分子初始化：基于 PYTHIA 8.3 生成初始部分子态。
  2. 部分子再散射 (PRS)：在强子化之前，引入部分子间的 $2 \to 2$ 硬散射（基于微扰 QCD）。
  3. 强子化：采用 Lund 弦碎裂模型（Lund string fragmentation）。在此阶段，除了标准的 NRQCD 机制外，还引入了**团簇坍缩（cluster collapse）**机制（当 $c\bar{c}$ 在相空间中足够接近且弦势能不足以产生轻夸克对时，直接形成夸克偶素）。
  4. 强子再散射 (HRS)：在强子化之后，模拟强子间的弹性及非弹性散射（包括 $J/\psi$ 与核子、介子的非弹性反应，如 $J/\psi + n \to \Lambda_c + D^-$ 等）。
• 物理机制：
  • Prompt $J/\psi$：包含非相对论性量子色动力学（NRQCD）框架下的色单态和色八重态贡献，以及团簇坍缩机制。
  • Non-prompt $J/\psi$：来自 $b$-强子的弱衰变。
  • Feed-down：考虑了来自 $\psi(2S)$、$\chi_{c0}$、$\chi_{c1}$ 和 $\chi_{c2}$ 衰变的贡献。
• 参数设置：
  • 模拟能量：$\sqrt{s} = 5.02, 7, 13$ TeV。
  • 事件数：4 亿个非弹性非衍射事件。
  • K 因子调整：为了与实验数据对齐，针对不同能量调整了 K 因子（分别为 0.68, 0.71, 0.75）。
  • 包含多重部分子相互作用（MPIs）和颜色重连（CR）机制。

3. 主要贡献 (Key Contributions)

1. 模型升级与机制完善：将 PACIAE 从仅考虑色单态升级为包含完整的 NRQCD（色单态 + 色八重态）框架，并首次在该模型中系统性地量化了团簇坍缩、非 prompt 过程以及各激发态 feed-down 的贡献。
2. 全面的运动学覆盖：同时研究了**中快度（mid-rapidity）和前向快度（forward rapidity）**区域，覆盖了 LHC 的三个主要能量点。
3. 定量分解：
   • 对 NRQCD 内部各通道（直接产生 vs. 不同激发态衰变）的贡献进行了精细分解。
   • 定量评估了部分子再散射（PRS）和强子再散射（HRS）对 $J/\psi$ 产额的独立及联合影响。
4. 新物理洞察：揭示了不同产生机制随能量和快度的演化规律，特别是色八重态碎裂机制在不同能量下的主导地位变化。

4. 关键结果 (Results)

A. 与实验数据的对比

• PACIAE 4.0 模拟的包含性 $J/\psi$ 横向动量（$p_T$）微分截面在中快度和前向快度区域均与 ALICE 等实验数据符合良好（大部分数据点偏差在 30% 以内）。
• 模型成功复现了 $J/\psi$ 产额随带电粒子多重数增加的趋势。

B. 不同产生机制的相对贡献

• 主导机制：NRQCD 过程是 $J/\psi$ 产生的主要来源。
• 能量依赖性：随着碰撞能量增加（从 5.02 TeV 到 13 TeV）：
  • NRQCD 的相对份额下降（例如中快度从 ~85% 降至 ~74%）。
  • 团簇坍缩（Cluster collapse）和非 prompt（来自 $b$-强子）过程的相对份额上升。
  • 原因：$b\bar{b}$ 产生的截面随能量增长快于 $c\bar{c}$，导致非 prompt 贡献增加；同时高能下色八重态碎裂机制增强。
• 快度依赖性：
  • 前向快度区域：NRQCD 和团簇坍缩的贡献增强，而非 prompt 贡献被抑制。
  • 原因：前向区域小 $x$ 部分子密度极高，有利于 NRQCD 和团簇坍缩；而大 $x$ 部分子稀缺，抑制了产生 $b$ 夸克所需的高能阈值过程。

C. NRQCD 内部通道分解

• 低 $p_T$：直接产生（Direct production）和 $\chi_{c2}$ 的衰变贡献最大。
• 高 $p_T$：直接产生和 $\chi_{c1}$ 的衰变贡献占主导。
• 演化趋势：
  • 直接产生份额随能量增加并趋于饱和。
  • $\psi(2S)$ 和 $\chi_{c1}$ 的 feed-down 份额随能量增加（受益于色八重态碎裂机制的增强）。
  • $\chi_{c0}$ 和 $\chi_{c2}$ 的份额随能量下降（$\chi_{c0}$ 受限于极小的色八重态矩阵元，$\chi_{c2}$ 在高能碎裂中处于自旋劣势）。
• 快度效应：前向快度下，直接产生、$\psi(2S)$ 和 $\chi_{c1}$ 的贡献略有抑制，而 $\chi_{c0}$ 和 $\chi_{c2}$ 的贡献增强，反映了不同机制对部分子动量分数 $x$ 的敏感性差异。

D. 再散射效应 (Rescattering Effects)

• 部分子再散射 (PRS)：对 $J/\psi$ 产额的影响微乎其微（变化小于 1%）。
• 强子再散射 (HRS)：
  • 显著抑制了来自 NRQCD 和团簇坍缩机制的 $J/\psi$ 产额（总抑制效应约为 8%）。
  • 对非 prompt $J/\psi$（来自 $b$-强子衰变）几乎没有影响（因为模型未考虑 $b$-强子的再散射，且 $b$-强子寿命较长，在强子化前未发生显著相互作用）。
  • 物理图像：HRS 中的非弹性散射（如 $J/\psi + \pi \to D + \bar{D}^*$）导致了 $J/\psi$ 的解离。

5. 意义 (Significance)

1. 理论验证：该研究证实了 PACIAE 4.0 模型结合 NRQCD 框架能够准确描述 LHC 能区 pp 碰撞中的 $J/\psi$ 产生，为理解强相互作用非微扰区域提供了可靠工具。
2. 基准建立：提供了高精度的 pp 碰撞基准数据，对于未来在核 - 核碰撞中分离 QGP 效应（热核物质）和冷核物质效应至关重要。
3. 机制解析：通过能量扫描和快度依赖分析，为实验上分离和约束 NRQCD 中的长距离矩阵元（LDMEs），特别是区分不同 $\chi_{cJ}$ 态的贡献提供了强有力的理论依据。
4. 动力学理解：明确了强子再散射是压低 $J/\psi$ 产额的关键因素之一，而部分子再散射在此过程中作用较小，深化了对重夸克偶素在强子介质中演化的理解。

综上所述，这项工作通过引入更完备的物理机制和更先进的模型版本，系统性地解构了高能 pp 碰撞中 $J/\psi$ 产生的复杂动力学，是理解重夸克偶素产生机制的重要一步。