Probing the onset of hydrodynamization in peripheral p-Pb collisions at $\sqrt{s_{NN}} =$ 5.02 TeV

要点

想象一下，你试图找出能够表现得像流体一样的最小水坑。如果你有一片浩瀚的海洋，它很容易流动。如果你只有一滴水，它可能只是静止不动或散开。但界限在哪里？在什么尺寸下，水分子的集合不再表现得像流体，而是开始表现为一个个混乱的粒子？

本文旨在寻找夸克 - 胶子等离子体（QGP）的确切“临界点”。

什么是 QGP？

将 QGP 想象成宇宙的“原始汤”。它是宇宙大爆炸后几分之一秒内存在的一种物质状态。在这种状态下，构成原子的基本单元（夸克和胶子）被熔化在一起，像一种超热、超密的液体一样自由流动。

通常，科学家通过将两个重原子（如铅）以接近光速的速度相互撞击来制造这种“汤”。但最近，科学家们注意到一个令人费解的现象：即使他们撞击的是更小的物体——例如单个质子撞击铅核（p-Pb 碰撞）——也会出现这种“液态汤”的迹象。

关键问题是：它真的是一种液体，还是仅仅是一堆混乱反弹的粒子？

实验：将质子撞击铅核

本文的作者希望找出这种“汤”能够用流体力学（描述流动液体的数学）来描述的最小尺寸。

他们使用了一个名为JETSCAPE的大型计算机模拟。可以将这个模拟想象成一个高科技视频游戏引擎，它分四个步骤重现整个碰撞过程：

1. 设置（TRENTo）：他们搭建舞台，将质子和铅核放置在起始位置。
2. 赛前（Freestreaming）：在“液体”形成之前，粒子在极短的时间内自由飞行。
3. 流动（MUSIC）：这是流体力学部分。模拟尝试将粒子视为流动的流体来处理。
4. 后果（iSS + SMASH）：随着“汤”冷却，粒子凝固成探测器可以观测到的实际质子、π介子和其他粒子。

测试：这种“汤”有多“液态”？

为了测试这种“汤”是否真的像流体一样行为，科学家们观察了一种称为椭圆流的现象。

类比：想象两辆汽车迎头相撞。如果它们完美圆形且正中相撞，碎片会呈圆形飞出。但如果它们稍微偏离中心撞击（擦边碰撞），碎片会更呈椭圆形飞出（像橄榄球）。

• 如果内部物质表现为完美流体，它会在该椭圆形方向上强烈地挤压流出。
• 如果物质仅仅是混乱粒子在四处反弹，那么椭圆形特征将很弱或根本不存在。

科学家们对“非中心”碰撞（即质子与铅核重叠部分较小的擦边碰撞）进行了模拟。他们问道：这种重叠可以小到什么程度，流体行为才会崩溃？

转折：“弛豫时间”旋钮

在真实流体中，当你推动流体与其做出响应之间存在延迟。在物理学中，这被称为剪切弛豫时间。

作者玩了一个花招：他们将这个“弛豫时间”旋钮调至极端设置。

• 他们问道：“如果流体响应非常迟钝会怎样？如果它非常迅速会怎样？”
• 他们在这些极端条件下观察椭圆流（即椭圆形特征）。

发现：临界点

随着他们模拟越来越“擦边”的碰撞（意味着涉及的物质量，即dN/dy，变得更小），他们观察了流体行为。

• 结果：当物质量下降到约**每单位快度 7 个粒子（dN/dy ≈ 7）**时，流体行为突然开始波动并崩溃。
• 隐喻：想象一群人试图像流体一样移动。如果你有 100 个人，他们流动顺畅。如果你有 10 个人，他们可能仍然流动。但如果减少到 7 个人，他们开始各自相互碰撞，平滑的“流动”就消失了。

论文得出结论，对于他们研究的能量下的质子 - 铅碰撞，当系统尺寸小于约 7 个粒子时，流体力学就不再适用。低于这个数值，“汤”太小而无法像液体一样行为；它只是一堆独立的粒子。

这为什么重要？

这有助于科学家理解自然界的基本极限。它告诉我们，物质的“液态”状态并非魔法；它有最小尺寸要求。如果系统太小，流体动力学的规则就不再适用，我们必须转而研究单个粒子。

作者还指出，他们的结果与他们之前关于更大规模碰撞（如铅 - 铅碰撞）的研究略有不同，这可能是因为他们这次使用的计算机模型更稳定，并且对“赛前”阶段的处理方式不同。

简而言之：他们找到了仍可被称为“流体”的最小夸克 - 胶子等离子体“水坑”，结果发现，这个水坑必须包含至少约 7 个粒子才能维持其整体性。

技术摘要

技术摘要：探测外围 p-Pb 碰撞中流体动力学化的起始

问题陈述
尽管相对论性流体动力学已成功建模了大型重离子碰撞（Pb-Pb、Au-Au）中的夸克 - 胶子等离子体（QGP），但其在小型碰撞系统（高多重数质子 - 质子和质子 - 原子核）中的适用性仍是激烈辩论的焦点。核心问题在于，这些小型系统中观测到的集体流特征究竟源于类 QGP 的热核物质效应，还是冷核物质效应。具体而言，作者旨在确定流体动力学描述失效的最小系统尺寸（以粒子多重数 $dN/dy$ 量化）。作者先前的工作已确立外围 Pb-Pb 和 Au-Au 碰撞中流体动力学化的起始点位于 $dN/dy \approx 10$；本研究将这一探究扩展至 $\sqrt{s_{NN}} = 5.02$ TeV 下更小的 p-Pb 系统。

方法论
本研究利用 JETSCAPE 框架（版本 3.6.4），这是一个旨在一致建模高能碰撞软区和硬区的模块化模拟工具。软区的模拟链包括：

1. 初始条件：TRENTo 模型基于核厚度函数的广义平均值生成初始能量密度，并纳入核子涨落。
2. 预平衡：基于无碰撞玻尔兹曼方程的 Freestreaming 模块填补了初始态与流体动力学阶段之间的空白。
3. 流体动力学：MUSIC 求解器实现了二阶粘性流体动力学（BRSSS 形式体系）。剪切弛豫时间（$\tau_\pi$）作为调节非流体动力学模式的二阶输运系数，是主要关注变量。
4. 强子后燃烧器：iSS（粒子化）和 SMASH（强子输运）模块处理冻结出及随后的强子再散射。

为了探测流体动力学化的起始，作者将剪切弛豫时间（$\tau_\pi$）调整至极端值。由于 $\tau_\pi$ 与剪切粘度与熵密度之比（$\eta/s$）的关系为 $\tau_\pi(T) = b_\pi \frac{\eta}{s(T)} \frac{1}{T}$，他们采用了两种不同的变化策略：

1. 改变归一化常数 $b_\pi$（从 1 到 20），同时保持温度依赖的 $\eta/s$ 固定。
2. 改变温度依赖的 $\eta/s$ 界限（使用下限 KSS 极限和上限贝叶斯极限），同时保持 $b_\pi$ 固定。

流体动力学行为的失效是通过监测椭圆流（$v_2$）及其涨落来识别的。理论论证表明，如果非流体动力学模式（由 $\tau_\pi$ 调节）的贡献超过流体动力学模式的贡献，流体描述将失效，导致流观测量出现急剧偏差。

主要贡献与结果

• 基准测试：该模型成功复现了 ALICE、ATLAS 和 CMS 关于带电粒子快度谱、横向动量（$p_T$）谱以及作为中心度（由 $N_{track}$ 定义）函数的平均横向动量的实验数据。
• 涨落分析：通过模拟外围 p-Pb 碰撞并改变 $\tau_\pi$，作者观察到椭圆流涨落随着系统从中心碰撞向外围碰撞移动而增加。
• 起始点确定：
  • 当改变 $b_\pi$ 时，在 75–85% 至 70–80% 的中心度范围内观察到椭圆流涨落的急剧增加。
  • 当改变 $\eta/s$ 时，涨落在 70–80% 中心度范围以下开始增加。
  • 这些涨落对应于带电粒子多重数 $dN/dy \approx 7$。
• 与先前工作的比较：与作者先前关于更大系统（起始点发现于 $dN/dy \approx 10$）的研究不同，p-Pb 中的转变出现在更低的多重数下。作者指出，低于推断起始点的偏差不如先前研究那样急剧，这可能是由于使用了比先前使用的 ECHO-QGP 更稳定的 MUSIC 流体动力学求解器。

意义与主张
本文声称提供了 5.02 TeV 下 p-Pb 碰撞中可由低阶流体动力学满意建模的最小 QGP 体积的估计值。主要结论是，流体动力学行为在大约 $dN/dy \approx 7$ 处失效。

作者对这一转变的性质保持了谦逊的态度，承认这是否代表一个尖锐的起始点还是一个平滑的转变区域仍然是一个未决问题。他们还强调，他们的研究隔离了软区；对起始点的明确确认理想情况下需要同时研究硬区（例如喷注淬火），以观察硬观测量和软观测量如何相互影响流体动力学化的确定。这项工作是将非流体动力学模式分析应用于约束小型碰撞系统中流体动力学模型有效性的具体实例。