Addressing uncertainties of model predictions for extensive air showers initiated by high energy cosmic rays

问题陈述
对超高能宇宙线（UHECR）的实验研究依赖于从广延大气簇射（EAS）特征中重建初级粒子属性，具体包括平均簇射极大值深度（$X_{max}$）和地面处的缪子数（$N_\mu$）。这些重建过程高度依赖于模拟大气中强子相互作用的蒙特卡洛（MC）生成器。量化这些生成器不确定性的过程面临重大挑战。虽然比较现有的不同模型可以提供一些见解，但目前的模型可能存在共同的缺陷或缺失关键物理机制，这可能导致不确定性被低估或产生与现实的系统性偏差。例如，皮埃尔·奥杰天文台（Pierre Auger Observatory）指出，当前的模型可能需要显著更大的 $X_{max}$ 值，才能与观测到的质量组成保持一致；而“缪子谜题”（muon puzzle）则表明，模型低估了缪子在 UHECR 诱发的簇射中的数量。

方法论
作者使用了一种新的强子碰撞 MC 生成器 QGSb，该生成器基于透明的理论形式（夸克-胶子弦模型），但具有足够的参数自由度来探索模型的不确定性。研究通过系统地修改特定模型参数来改变对 $X_{max}$ 和 $N_\mu$ 的预测，同时严格遵守加速器数据（如 LHC、TOTEM、ALICE、NA61/SHINE、EHS-NA22）所施加的约束。

调查重点在于两个主要领域：

1. 修改 $X_{max}$： 作者探索了增加或减少预测的簇射极大值深度的机制。这涉及调整质子-空气非弹性截面的能量依赖性、相互作用的非弹性（$K_{inel}$）以及组成部分子的动量分布。具体的修改包括：
   • 调节“半硬 Pomeron”的过临界性（$\Delta_{sh}$）。
   • 改变海（sea）(反)夸子的动量分布指数（$\alpha^{(sh)}_{sea}$）。
   • 实现一种“双夸克破碎”（diquark breaking）机制，即在多次散射过程中破坏核子中价夸克的相干性，从而可能增加非弹性。
2. 修改 $N_\mu$： 作者通过改变 $\pi^0$ 与稳定次级强子（核子、卡介子）之间的能量平衡，研究了增强缪子含量的方法。由于 $N_\mu$ 主要由稳定强子的前向产生决定，研究重点在于：
   • 增强 $t$ 道 $\pi$ 交换，以有利于带电 $\pi$ 产生而非中性 $\pi$。
   • 通过修改弦碎裂参数（例如，从真空中产生奇异夸克对或双夸克-反双夸克对的概率），增加投影碎片区的前向（反）重子和卡介子的产生。

核心贡献与结果

• 增加 $X_{max}$： 研究发现，增加 $X_{max}$ 最有效的方法是假设附着在颜色场弦上的组成部分子具有较软的动量分布。通过将 $\alpha^{(sh)}_{sea}$ 从默认的 0.5 增加到 0.8 和 0.95（Tune 3 和 4），质子-氮碰撞的非弹性降低了高达 11%，从而使 $X_{max}$ 在 $10^{20}$ eV 时增加了高达 $\approx 40$ g/cm$^2$。

  • 约束： 这些极端的修改被最近的 ALICE $\sqrt{s}=13$ TeV 数据证伪。ALICE 实验观察到中心带电强子多重数与前向中子产额之间存在强烈的反相关关系。而 Tune 3 和 4 预测的反相关性要弱得多，这与数据相矛盾。此外，这些修改导致的极大值缪子产生深度（$X^\mu_{max}$）将需要比铁更重的初级原子核才能匹配皮埃尔·奥杰数据，这在物理上是不太可能的。
  • 替代方案： 适度降低 Pomeron 的过临界性（$\Delta_{sh} = 0.21$ 对比 0.22，Tune 2）会产生较小的 $X_{max}$ 增幅（$\approx 7$ g/cm$^2$），并与 ALICE 相关性保持一致，尽管它略微低估了领先中子的能量损失。

• 减少 $X_{max}$： 作者测试了“双夸克破碎”机制（Tune 5），即多次散射破坏双夸克相干性，理论上会增加非弹性。尽管该机制预测在高能下非弹性增加 3%，但它对平均 $X_{max}$ 没有产生明显的改变（在 1–2 g/cm$^2$ 的统计不确定性范围内）。作者得出结论，默认模型的 $X_{max}$ 预测可能已经接近下限，因为边缘相互作用（主导平均值的相互作用）不受该机制的影响。

• 增强 $N_\mu$： 研究证实 $N_\mu$ 对稳定强子的前向产生高度敏感。

  • 增强前向卡介子产生（Tune 7）使 $N_\mu$ 增加了高达 $\approx 2\%$。然而，该 Tune 在其他碰撞系统（如 $\pi p$ 碰撞中的 $K^+$）中高估了前向卡介子产额，且无法匹配某些实验能谱。
  • 通过修改“残余”夸克的强子化以包含更多重子成分，增强前向（反）核子产生（Tune 8）使 $N_\mu$ 增加了高达 $\approx 5\%$。虽然该 Tune 改善了与某些 NA61/SHINE 数据的符合度，但导致其在 360 GeV/c 下的质子/反质子产额比 LEBC-EHS 数据高出一个数量级。
  • 研究强调，目前 $N_\mu$ 的不确定性受限于不同加速器测量结果之间关于碎片区（反）质子产额的严重矛盾。

意义与主张
本文声称提供了一个详细的“参数空间”映射，用于在当前加速器数据的限制内对强子相互作用模型进行调优。其主要意义在于：

1. 排除极端情况： 它证明了通过（软）部分子分布对 $X_{max}$ 进行的最剧烈修改与 LHC 前向物理数据（ALICE）以及奥杰缪子深度数据是不相容的。
2. 确定 $X_{max}$ 降低的极限： 它挑战了文献中关于双夸克分裂显著降低 $X_{max}$ 的说法，表明在受限于定标靶质子能谱的情况下，此类机制对平均簇射极大值的影响微乎其微。
3. 量化“缪子谜题”的不确定性： 它确立了在 QGSb 框架内增加 5% 的 $N_\mu$ 是理论上可行的，但这需要关于前向（反）核子产生的假设，而这些假设目前与现有的实验数据相冲突。

作者总结道，虽然 QGSb 生成器允许探索这些不确定性，但来自加速器数据（特别是关于前向粒子产生和相关性）的当前约束，显著限制了在不引入尚未考虑的新物理机制的情况下，解决 EAS 观测与模型预测之间差异的能力。