Path-length dependence of parton energy loss across collision systems: a Bayesian analysis of charged-particle RAA, consistent with a universal exponent from O+O to Pb+Pb

技术摘要：跨碰撞系统的部分子能量损失的路径长度依赖性

问题陈述
夸克-胶子等离子体（QGP）中部分子能量损失机制的控制规律，编码在平均能量损失 $\langle \Delta E \rangle$ 如何随中介路径长度 $L$ 变化的标度关系中。理论框架预测了不同的标度行为：碰撞能量损失呈线性标度（$\langle \Delta E \rangle \propto L$），介质诱导辐射损失呈二次方标度（$\langle \Delta E \rangle \propto L^2$），而强耦合（全息）模型则预测立方标度（$\langle \Delta E \rangle \propto L^3$）。从历史上看，提取这种标度指数的复杂之处在于，通过分析单一大型系统（如 Pb+Pb）内的中心度变化来提取，这会引入来自中心度选择和基于 Glauber 归一化（$T_{AA}$）在边缘碰撞中的偏差。LHC 上轻离子碰撞数据（O+O, Ne+Ne）的出现提供了一个互补的杠杆：通过改变系统大小（$A$）并在最小偏置（minimum bias）条件下进行研究，可以探测与重离子边缘碰撞相当的路径长度，从而隔离出系统大小依赖性。

方法论
作者对四种碰撞系统（O+O, Ne+Ne, Xe+Xe, 和 Pb+Pb）的带电粒子核修正因子（$R_{AA}$）进行了贝叶斯分析，这些系统的质量数 $A$ 从 16 跨越到 208。

1. 前向模型（Forward Model）： 一个参数化模型将 $R_{AA}$ 与分数能量损失 $\Delta p_T$ 联系起来。该模型假设 $\Delta p_T \propto \rho G^n$，其中 $\rho$ 是相对介质密度，$G$ 是无量纲几何比，$n$ 是待确定的系统大小指数。谱指数 $a(p_T)$ 使用来自 $pp$碰撞的数据驱动测量值进行固定，并应用$x_T$ 标度以适用于不同的碰撞能量。
2. 几何结构： 碰撞几何使用两种独立的 Glauber 实现进行计算：一种是光学 Glauber 模型，另一种是自定义的蒙特卡洛 Glauber 模型（通过具有硬核排斥的核子位置采样）。关键几何代理量包括 $A^{1/3}$、$\langle N_{part} \rangle^{1/3}$ 以及平均退出路径长度 $\langle L \rangle$。
3. 统计框架：
   • 似然函数： 使用了包含分箱间协方差矩阵的相关高斯似然函数。系统误差通过相关长度 $\xi = 4$ 个分箱进行建模，以处理部分相关性，避免因完全相关的系统误差导致过度约束。
   • 推断： 使用仿射不变系综采样器（EMCEE）对参数（$\kappa, n, \beta$）的后验分布进行采样。
   • 模型选择： 通过嵌套采样（DYNESTY）计算贝叶斯证据，以比较固定指数（$n=1, 2, 3$）的模型。
   • 验证： 该框架包含一个使用正规化流神经后验估计器的“基于模拟的推断”（SBI）支柱，并通过基于模拟的校准（SBC）和留一系统法（LOSO）交叉验证进行验证。此外还训练了一个高斯过程模拟器来预测被留出的系统。
4. 鲁棒性： 分析通过 160 种变体进行了测试，涉及不同的几何代理、协方差长度、谱基线和先验分布。

关键结果

• 有效指数： 对于 $\langle N_{part} \rangle^{1/3}$ 几何结构，分析提取出的有效系统大小指数为 $n = 1.78 \pm 0.15 \text{ (stat)} \pm 0.05 \text{ (syst)}$。
• 模型选择： 贝叶斯证据果断支持辐射标度（$n=2$），优于碰撞标度（$n=1$）和强耦合（$n=3$）假设。对数证据差异为 $2\Delta \ln Z \approx -29$（对于 $n=1$）和 $-48$（对于$n=3$），两者均超过了“决定性”阈值。
• 微观界限： 通过分析在存在能量损失涨落的情况下微观指数与有效指数之间的映射，作者确定了涨落只能使有效指数低于其微观指数。因此，该测量结果为微观指数设定了一个下界（$n_{micro} \ge 1.78$），有效地排除了在固定几何下的纯碰撞能量损失（$n_{micro}=1$）。
• 普适性： 贝叶斯因子检验发现，在小系统（O+O, Ne+Ne）和大型系统（Xe+Xe, Pb+Pb）之间不存在能量损失机制发生变化的证据，支持了全范围内的统一指数。
• O+O 抑制： 该框架为 O+O 中非零的部分子能量损失提供了决定性的贝叶斯证据（$2\Delta \ln Z \approx 110$），量化了最小系统中的抑制起始点。
• 输运系数： 提取的能量损失量对应于有效喷注输运系数 $\hat{q}/T^3 \approx 2\text{--}5$，这与独立的 JETSCAPE 测定结果一致。

意义与主张
本文声称提供了首次经过统计校准的、基于系统大小的路径长度指数的部分子能量损失提取。其主要贡献包括：

1. 一个鲁棒的有效指数： 它确立了 $n \approx 1.78$ 作为一个鲁棒的、几何层面的可观测物理量，该物理量独立于特定的输运实现，可作为未来多阶段能量损失计算的基准。
2. 排除碰撞机制： 通过涨落界限，它证明了即使考虑到涨落，纯碰撞能量损失也与观察到的系统大小趋势不相容。
3. 方法论的严谨性： 该工作引入了一个全面的贝叶斯框架，用于传播相关不确定性，通过闭合测试验证区间覆盖率，并利用基于模拟的推断来交叉检查结果，确保了提取过程在统计上的可辩护性。
4. 预测能力： 该框架为未来的轻离子测量（Ar+Ar 和 Kr+Kr）提供了可证伪的预测，预测 $R_{AA}(10 \text{ GeV})$ 分别为 $0.55 \pm 0.03$ 和 $0.41 \pm 0.03$。

局限性与范围
作者明确指出，所得结果是一个有效指数，而非对微观路径长度依赖性的直接测量。该分析受限于在最小偏置系统大小扫描中介质密度与路径长度之间的简并性（两者都随 $\langle N_{part} \rangle^{1/3}$ 缩放），这意味着提取的 $n$ 结合了这两者的效应。本文并不声称在没有进一步输入（如中心度微分数据或方位角各向异性 $v_2$ 测量）的情况下解决了微观机制（例如，区分 BDMPS-Z 与有限不透明度修正）。其结论严格受限于所分析的四个系统以及所采用的具体前向模型假设。