Event-by-event fluctuations of elliptic flow in ultrarelativistic O+O collisions

技术摘要：超相对论性 O+O 碰撞中椭圆流的逐事件涨落

问题陈述
本文研究了在 $\sqrt{s_{NN}} = 5.36$ TeV 的氧-氧（O+O）碰撞中，各向异性流（特别是椭圆流 $v_2$）的本质。对于大核（如 Pb+Pb）而言，主流观点认为椭圆流主要源于由碰撞参数决定的重叠区域的几何不对称性；然而，作者提出 O+O 碰撞可能代表了一个不同的机制。在较小系统中，由于核子数量有限，能量沉积的逐事件涨落可能会主导初始几何结构，这使得 O+O 碰撞在流的起源方面可能更接近于小系统（如 p+p），而非重离子碰撞。本研究旨在量化这些涨落，并确定一组缩减后的涨落模态是否能够准确重现椭圆流的统计特性。

方法论
作者采用了一个结合了 McDipper（初始态）、MUSIC（流体动力学演化）和 iSS（粒子化）的全三维数值框架。

• 初始态： 使用核晶格有效场论（NLEFT）生成 $^{16}$O 核子的构型。这些构型经过随机旋转并与随机的碰撞参数进行碰撞。基于颜色玻璃凝聚有效场论（GBW 模式）的 McDipper 代码生成了三维能量密度分布 $\Phi(x, y, \eta_s)$，其中包含了亚核热点和厚度涨落。
• 流体动力学： 这些分布作为 MUSIC 的初始条件，以恒定的比剪切粘度（$\eta/s = 0.16$）和随温度变化的体粘度进行演化，直到冻结温度为 155 MeV。
• 分解框架： 本研究的一个关键方法论创新是将三维初始态能量密度分解为一个平均分布 $\bar{\Psi}$ 和一组不相关的涨落模态 $\{\Psi_l\}$（及其系数 $c_l$）。这是通过对能量密度涨落的自相关矩阵进行特征分解实现的。该方法首次被应用于非玻尔兹曼不变（non-boost-invariant）的三维初始态。
• 线性响应： 作者利用线性响应框架，将观测物理量（离心率 $\varepsilon_2$ 和流 $v_2$）近似为涨落模态的线性函数。他们将“逐事件”模拟（每个中心度类 5000 个事件）与通过采样高斯分布系数生成的“线性化”观测值进行了对比。

核心贡献与结果

1. 涨落的主导作用： 研究表明，在 O+O 碰撞中，初始椭圆离心率（$\varepsilon_2$）主要由能量沉积的逐事件涨落而非有限的碰撞参数驱动。因此，即使在非中心碰撞（20–30% 中心度）中，平均初始态 $\bar{\Psi}$ 也几乎是圆形的，其平均离心率非常小（$\varepsilon_2 \approx 0.087$）。这与重离子碰撞形成鲜明对比，在后者中，平均态呈现显著的杏仁形。
2. 涨落权重： 研究发现 O+O 中平均态的相对权重（$\bar{w}$）仅约为 $2\%$，显著低于 Pb+Pb 中观察到的 $\sim 13\%$。这从定量角度证实了在 O+O 中，涨落是初始几何结构的主要贡献来源。
3. 模态分解精度： 作者展示了使用一小组涨落模态（具体为 10 个模态）可以较为准确地重现 $\varepsilon_2$ 和 $v_2$ 的联合概率分布。“线性模型”和“解析高斯涨落近似”（G.f.a.）在 $\varepsilon_2 \lesssim 0.6$ 时，与全逐事件模拟的结果高度一致。
4. 非线性局限性： 在大离心率情况下（$\varepsilon_2 > 0.6$），由于非线性流体动力学响应（特别是关系式 $v_2 \approx \kappa_{2,2}\varepsilon_2 + \kappa_{2,222}\varepsilon_2^3$ 中的三次项），线性近似偏离了逐事件数据。线性模型无法捕捉分布中的“香蕉”形状以及由这些非线性效应导致的更宽的 $v_2$ 尾部。
5. 纵向去相关： 论文指出，虽然积分流可以通过少量模态来捕捉，但对纵向结构敏感的观测物理量（如离心率的纵向去相关 $r_2$）需要多得多的模态（数百个）才能重现逐事件结果。这归因于快度奇（rapidity-odd）涨落模态的重要性，这类模态出现在模态谱的高阶项中。

意义与主张
本文声称，O+O 碰撞不仅仅是重离子碰撞的“缩放版”。相反，O+O 中椭圆流的起源在本质上是不同的，它由初始态涨落而非全局碰撞参数几何结构所主导。这表明 O+O 碰撞为研究涨落作为主要驱动力的集体现象起始阶段提供了一个独特的实验室。

作者谦虚地声称，他们的逐模态分解框架是一种可行的、成本更低的替代方案，可以用于描述椭圆流的整体统计行为，从而替代全逐事件三维流体动力学模拟。他们证明了，只要考虑到涨落相对于平均几何结构的支配地位，有限的一组涨落模态足以重现 $v_2$ 和 $\varepsilon_2$ 的平均值及概率分布。然而，他们也承认线性近似在处理高离心率下的非线性响应时存在局限性，且更复杂的观测物理量（如纵向去相关）需要显著增大的模态基数。