Effects of the centrality determination method for the equation of state and nucleonic observables from Au+Au collisions at $\sqrt{s_{NN}}$ = 2.4 GeV

本研究利用 $\sqrt{s_{NN}}$ = 2.4 GeV 下 Au+Au 碰撞的 UrQMD 模拟表明，中心度确定方法的选择显著影响核物态方程的提取，揭示了基于 Glauber Monte Carlo 的方法引入了比物态方程本身更大的不确定性，而几何解释则与动力学多重数选择保持一致。

要点

想象一下，你正试图通过观察两个巨大的碗互相撞击，来理解一种非常稠密、高温的汤（核物质）的配方。物理学家称之为“重离子碰撞”。当这些碗（金原子）相撞时，它们会挤压内部的物质，创造出仅在中子星核心或早期宇宙中才能发现的条件。

这篇论文的目标是弄清楚这种“核汤”是“硬”还是“软”。这种特性被称为状态方程（Equation of State, EoS）。你可以把它想象成在问：这碗汤是果冻做的（软），还是混凝土做的（硬）？

问题所在：我们如何测量“撞击力度”？

在真实的实验中，你无法直接看到碰撞的中心。它发生得太快，且规模太小。因此，科学家必须通过计算飞出的“碎片”（粒子）来推测撞击到底有多重。

这篇论文的作者提出了一个简单的问题：“我们猜测撞击强度的‘方式’重要吗？我们的猜测会改变我们对汤是果冻还是混凝土的结论吗？”

他们测试了三种不同的猜测撞击“强度”（中心度）的方法：

1. 碎片计数法 (Mch)： “我们假设看到的粒子越多，撞击就越重。”
2. 几何猜测法 (bf)： “我们根据球体重叠程度的一个简单图解来假设撞击强度。”
3. 计算机模型猜测法 (br)： “我们使用一个著名的计算机模拟（Glauber 模型），根据参与其中的人数（粒子数）来猜测撞击强度。”

实验：一场虚拟的碰撞测试

研究人员并没有撞击真实的原子。他们使用了一个名为 UrQMD 的超级计算机模拟程序，将金原子以特定的能量（2.4 GeV）撞击在一起。他们运行了两次碰撞：

• 一次是“果冻汤”（软 EoS）。
• 一次是“混凝土汤”（硬 EOs）。

然后，他们使用上述三种不同的猜测方法观察了结果。

研究结果：“猜测游戏”至关重要

以下是他们的发现，使用了简单的类比：

1. “碎片计数法”对比“几何猜测法” (Mch vs. bf)
当他们使用碎片计数法并将其与简单的几何猜测法进行比较时，结果出奇地相似。

• 类比： 想象一下，你试图通过数破碎的车灯来猜测车祸的严重程度。如果你数了车灯，你会对车祸的严重程度有一个相当好的把握，而且这与车辆重叠程度的简单图解非常吻合。
• 结果： “猜测”带来的不确定性并没有破坏对汤之硬度的测量。 “果冻”和“混凝土”的结果依然界限分明。

2. “碎片计数法”对比“计算机模型” (Mch vs. br)
这是情况变得混乱的地方。当他们使用碎片计数法并将其与计算机模型 (Glauber) 进行比较时，结果大相径庭。

• 类比： 想象一下，你使用一个高科技天气应用来猜测车祸的严重程度。该应用假设“破碎部件越多 = 撞击越重”，这种逻辑在处理大型高速公路车祸时有效，但在处理小型低速追尾事故时却失效了。
• 结果： 在这个特定的能量水平（2.4 GeV）下，计算机模型是错误的。它误判了碰撞的强度。由于这个错误的猜测，导致“果冻汤”和“混凝土汤”之间的差异变得模糊不清。事实上，由于使用错误的猜测方法所造成的误差，竟然比“果冻”和“混凝土”之间的实际差异还要大！

核心结论

论文得出结论，如果你想要准确测量核物质的“硬度”：

• 不要信任旧的计算机模型 (Glauber) 来处理这类特定类型的碰撞。这就像是用城市地图在森林里导航；规则并不适用。
• 坚持使用几何逻辑或直接的碎片计数法。 这些方法在这些能量水平下更符合现实。

简而言之： 如果你用错误的尺子去测量碰撞，你可能会误以为汤是“果冻”，而它实际上是“混凝土”，反之亦然。定义碰撞的方式与你试图测量的物理学本身同样重要。作者警告说，为了得到正确的配方，我们需要一种一致的方法，将“碎片计数”与“实际碰撞几何结构”匹配起来，尤其是在较低的能量水平下。

技术摘要

技术摘要：$\sqrt{s_{NN}} = 2.4$ GeV Au+Au 碰撞中中心度确定对核子观测量的影响

问题陈述
在中等能量重离子碰撞（HIC）分析中，碰撞中心度的确定仍然是主要的系统误差来源。在探测超饱和密度下的核物质状态方程（EoS）时，这一点尤为关键。虽然理论模拟使用碰撞参数（$b$）作为直接输入，但实验中的中心度是通过最终态观测量（例如带电粒子数 $M_{ch}$）利用理论模型推导出来的。不存在将碰撞参数与这些观测量进行严格一一映射的关系。以往的研究表明，所选事件的真实碰撞参数分布所具有的宽度会显著影响集体流。此外，标准中心度确定方法（如 Glauber 蒙特卡洛 MC 模型）在低能（如 2.4 GeV）下的适用性值得商榷，因为该模型假设粒子数与参与者数量成正比——这一假设在高能条件下更为成立。本研究旨在定量评估不同中心度确定方法所带来的不确定性，并评估其对最终态 EoS 敏感观测量的影响。

方法论
本研究利用 UrQMD 模型模拟 $\sqrt{s_{NN}} = 2.4$ GeV（对应于 1.23A GeV 束流能量）的 Au+Au 碰撞。

• 状态方程 (EoS)： 采用了两种动量依赖型 EoS 情景：软 EoS（$K_0 = 200$ MeV，标记为 SM）和硬 EoS（$K_0 = 380$ MeV，标记为 HM）。
• 中心度确定方法： 对三种不同的方法进行了比较，以在特定的中心度区间（0–10%、10–20%、20–30% 和 30–40%）内筛选事件样本：
  1. $M_{ch}$： 基于所有带电粒子的多重数进行筛选。
  2. $b_f$： 基于输入碰撞参数范围的几何解释。
  3. $b_r$： 一种源自 Glauber MC 模型（基于参考文献 [7]）的方法，该方法将多重数映射到碰撞参数范围。
• 观测量： 研究分析了自由质子的产额及其快度分布、定向流（$v_1$）和椭圆流（$v_2$）系数（具体为斜率 $v_{11}$ 和中快度值 $v_{20}$）、横向动量（$p_T$）分布以及平均横向动量（$\langle p_T \rangle$）。

主要结果

• 碰撞参数分布： 通过 $M_{ch}$、$b_f$ 和 $b_r$ 筛选出的事件样本在真实碰撞参数分布上存在显著差异。虽然 $M_{ch}$ 和 $b_f$ 的分布相对一致，但 $M_{ch}$ 与 $b_r$ 之间的偏差非常大。具体而言，对于外围碰撞（例如 30–40%），$b_r$ 样本中处于预期范围之外的事件比例与 $M_{ch}$ 样本相比，偏差高达 60%。
• 质子产额： 自由质子产额对 EoS 硬度和中心度方法均具有敏感性。
  • 当使用几何方法（$b_f$）时，由中心度选择引入的不确定性弱于由 EoS 引起的效应。
  • 当使用基于 Glauber MC 的方法（$b_r$）时，中心度相关不确定性的影响变得更加显著，甚至超过了 EoS 的效应，尤其是在外围碰撞中。$M_{ch}$ 与 $b_r$ 选择之间的产额比随中心度增加而显著上升。
• 集体流： 定向流斜率（$v_{11}$）和椭圆流（$v_{20}$）主要受 EoS 硬度的支配。虽然 $v_{11}$ 对中心度确定方法基本不敏感，但 $v_{20}$ 显示出轻微的敏感性，$b_r$ 选取的事件产生的 $v_{20}$ 值略大于 $M_{ch}$ 选取的事件。
• 横向动量： $p_T$ 分布和平均横向动量遵循与产额相似的趋势。当使用 $b_r$ 时，$p_T$ 分布对中心度定义具有强依赖性，而平均横向动量 $\langle p_T \rangle$ 则被发现对 EoS 和中心度确定方法均相对不敏感。

意义与结论
本文得出结论，建立带电粒子多重数（$M_{ch}$）与碰撞参数之间严谨且一致的映射关系，对于对高密度核 EoS 施加定量约束至关重要。

• 方法的有效性： 在 SIS18 能量（$\sqrt{s_{NN}} = 2.4$ GeV）下，中心度的几何解释（$b_f$）保持有效，并与动力学多重数选择保持一致。相比之下，Glauber MC 基的中心度确定方法（$b_r$）在该能量区间变得不可靠，这可能是由于“粒子多重数严格正比于参与者数量”这一假设的失效所致。
• 对 EoS 约束的影响： 本研究表明，使用不恰当的中心度确定方法（特别是低能下的 Glauber MC）引入的系统误差可能与 EoS 本身的物理效应相当甚至更高。因此，未来的 EoS 约束必须考虑所使用的特定中心度选择方法，确保实验观测量与几何中心度之间的映射在所研究的能量区间内具有物理上的合理性。

未来展望
作者建议，未来的工作应采用贝叶斯推断技术来定量评估模拟与实验数据之间的系统偏差，重点关注不同中心度类别的涨落与相关性。此外，有必要对 Glauber MC 方法的能量依赖性和模型假设进行更详细的研究，以提高低能及中等能量下 EoS 研究的可靠性。