中高能重离子碰撞——集体流·对称能·QCD相图

把两个金核加速到光速的 99.99%，然后让它们对撞。撞出来的不是碎金子。是一团温度高达几万亿度的火球——比太阳核心热十万倍——里面质子和中子不再存在，它们化成了夸克和胶子的等离子体。这团火球的直径不到十的负十四次方米，寿命不到十的负二十二秒。但它内部的物理，掌握着宇宙诞生后最初几微秒的规则。

这就是中高能重离子碰撞。近五十年，它始终是核物理和粒子物理交汇处最活跃也最昂贵的实验方向。

说两件事。第一件，怎么"看"这团火球。火球在膨胀。膨胀的时候，内部压力把粒子往外推。如果火球不是一个完美的球——它不可能是——某些方向的压力大，某些方向小。这些压力差，会刻在飞出火球的粒子的角分布里。

把一个碰撞事件的所有粒子按方位角统计——它们是各向同性的，还是偏向某个方向的？偏向多少？这就是集体流。两种最基本的流：指向流（v₁）和椭圆流（v₂）。

指向流告诉你，碰撞早期的"反弹"有多强。两个核撞在一起，它们不会正好对心——总是偏一点。碰撞区域的核物质往哪个方向飘？往两个核飞来的方向，还是往垂直于这些方向的面？指向流就量这个偏向。如果核物质状态方程硬——核物质不可压缩——反弹就大，指向流就强。如果软——像海绵，一捏就扁——指向流就弱。

椭圆流告诉你，碰撞区域的"杏仁形"有多大。非对心碰撞产生一个杏仁形的重叠区。这个杏仁的短轴方向压力梯度大、长轴方向压力梯度小。火球膨胀时，粒子更多沿短轴方向喷出来。椭圆流就量这个"长短轴差"。它比指向流更直接地反映早期压力——因为它在碰撞后极短的时间内就被冻结了，后续的演化对它影响小。

这两个流——v₁ 和 v₂——是重离子碰撞物理里最核心的两个观测量。它们被测量了三十年，至今还在产出令人意外的结果。

最令人意外的一个，是指向流的斜率——dv₁/dy——在某个碰撞能量下会出现一个极小值。这个极小值在 1990 年代就被提出了：如果核物质在某个密度和温度下发生一阶相变——从强子气体突然变成夸克胶子等离子体——状态方程会在那个点"软化"——压力突然变小。一阶相变的标志不是升温时突然涨大的压力，而是升温时平坦的压力。更软的状态方程导致更弱的指向流。所以 dv₁/dy 会在相变能区出现一个坑——一个极小值。

STAR 合作组在 RHIC 的束流能量扫描（BES）项目中，系统地追踪了这个坑。BES-II——第二期扫描，统计量比第一期高 3-5 倍——在 2024-2025 年交出了精确到让人屏息的结果。

质子的 dv₁/dy 极小值出现在√s≈11.5-19.6 GeV（质心系每核子对能量），对应的重子化学势 μB≈200-300 MeV。更重的粒子——比如净 kaon——极小值出现在更低的能量：√s≈4.5-7.7 GeV。这个"质量排序"是对的：更重的粒子产生得晚，感受到的压力差来自更晚期、更软的状态方程段。

2024 年 STAR 的一篇新分析（Duckworth 等人）发现，如果把质子的 v₁ 分解成一个"介质生成"分量和一个"迁移质子"分量，后者的标度行为从 200 GeV 一路延续到约 10-14.5 GeV——但在这个能区以下断裂了。断裂意味着：介质的状态变了。低于大约 7.7 GeV——μB 大于约 420 MeV——碰撞产生的物质不再是一个光滑的连续流。可能是强子气体，可能是一个稠密的重子物质，可能和"上面"的夸克胶子等离子体有一个清晰的边界。

这个边界——如果存在——叫 QCD 临界端点。找到它，是 BES-II 最高优先级的目标。

临界点怎么找？用涨落。二阶相变的关键特征是涨落发散——接近临界点时，系统的密度涨落会变得越来越大，因为恢复力（压缩率）趋于零。在实验里，你测量碰撞事件中质子数的分布——不是平均质子数，而是每个事件中质子数偏离平均值的程度。二阶（方差）、三阶（偏度）、四阶（峰度）累积量——特别是净质子数的四阶累积量 C₄/C₂——应该对临界点的存在产生"非单调"的能量依赖。

2025 年，STAR 合作组在 Physical Review Letters 上发表了 BES-II 的最精确净质子涨落测量。7.7 到 27 GeV，六档能量，总计约 9 亿个 Au+Au 碰撞事件。

结果：C₄/C₂ 在 √s≈19.6 GeV（μB≈206 MeV）处出现了一个下沉——一个极小值。偏离所有无需临界点的模型——UrQMD 输运模拟、强子共振气体模型、边缘碰撞数据——达到了 3.2-4.7 个标准差。

这不是发现。4σ 在高能物理中是"提示"，不是"确认"。但方向和形状是对的——理论预言的就是一个下沉，而且这个下沉的精确位置和宽度，可以用来约束临界点的位置。如果当前对 QCD 伪临界线的参数化是对的——只到 μB² 阶——那临界点应该在 T<132 MeV 且 μB>470 MeV 的地方。也就是 √s<7.7 GeV。BES-II 的固定靶模式（FXT）正好覆盖了 √s=3.0-4.5 GeV，这个区域的 C₄/C₂ 目前和强子基线一致——没有看到下沉的回归（从下沉回升上去的"峰"）。但如果临界点在更高或更低的 μB 上，需要更多、更低的能量覆盖。

这是 2025 年最活跃的待解决区之一。

再说第二件事：状态方程和对称能的约束。

把一颗中子星和一锅重离子碰撞数据放在一起，看看能不能用同一个状态方程来描述两者。这是过去两年的核心方法论突破。

中子星的观测——最大质量至少 2 倍太阳质量，半径大约 11-13 公里——给状态方程设定了天文学约束。重离子碰撞——集体流、粒子产额——给状态方程设定了实验约束。这两条约束在 2024-2025 年开始收敛。

Cozma 在 2024 年（PRC 110, 064911）用 dcQMD 输运模型同时对 FOPI 实验的 v₁、v₂ 和 stopping 做多变量拟合，得到：不可压缩系数 K₀=230±10 MeV（中软），对称能斜率 L=63±10 MeV。关键一步是介质内核子-核子截面被压低到真空值的约 40%——如果没有这个压低，你存不了同时拟合 v₁ 和 v₂。

王等人在 2024 年（arXiv:2406.07051）用 IQMD 和高斯过程模拟器对 Au+Au 400 MeV/nucleon 的质子椭圆流做了贝叶斯推断。如果包含动量相关势，压缩系数 K₀=191^{+4}_{-6} MeV——很软。如果不加动量相关势，K₀≈235 MeV——和巨共振值一致。这一差异不是系统误差，而是物理：动量相关性让核子在较远距离上就能感到对方的斥力，等效地软化了状态方程。所以你不能同时忽略动量相关性和测到椭圆流。

Kireyeu 等人在 2024-2025 年（arXiv:2411.04969）用 PHQMD 对比了软静态、硬静态、软动量相关三种势，通过同时拟合 HADES 和 FOPI 的质子及轻核 v₁ 和 v₂ 数据，结论：动量相关软势胜出。而且椭圆流对轻核质量数 A 存在标度行为——v₂(pT)/A 在不同核之间几乎相同——表明轻核主要是通过核子聚合产生的。

对称能——状态方程中区分中子和质子贡献的那部分——正在通过多个实验通道收紧。

2025 年，Ciampi 等人在 Physics Letters B 上发表了目前对对称能量精确的约束。他们用 INDRA-FAZIA 探测器阵列测量了费米能区（约 30-50 MeV/nucleon）重离子碰撞的同位旋输运比——一种比较中子过丰弹丸和缺中子弹丸打在对称靶上的中子-质子比差异的观测量——然后和 BUU 输运计算比较。结果：在饱和密度 ρ/ρ₀≈1.01 处，对称能 S=29.0±0.7 MeV。这个值和从天体物理观测（中子星半径、潮汐形变）独立导出的值高度一致。

Tsang 等人在 2024 年（PLB 853）用 SπRIT 时间投影室的 Sn+Sn 270 MeV/nucleon 数据做了贝叶斯分析，同时约束 S₀、L、同位旋标量和矢量有效质量。一个重要发现：椭圆流对对称能的敏感度被质子-中子有效质量劈裂——另一种核介质效应——强烈调制。你不把有效质量劈裂摸清楚，你就抽不出干净的对称能。

最引人注目的跨学科约束来自 2025 年 PLB 上的一篇论文：用同一个手征平均场模型（Chiral Mean Field Model），同时拟合中子星观测和 HADES 重离子碰撞数据。同一个状态方程，从饱和密度到 5 倍饱和密度的行为，把中子星最大质量（>2M⊙）和 HADES 的亚阈 K⁺ 产额两块拼图拼在了一起。结果消除了历史上"轻的 K⁺ 产额要软 EOS，重的脉冲星要硬 EOS"的矛盾——因为平均场势是动量相关的，它在低动量下软化（有利于 K⁺ 产额），在高动量下硬化（有利于脉冲星质量上限）。

再说第三件事：亚阈产生。

在核子-核子碰撞中，产生一个 K⁺ 介子需要至少约 1.6 GeV 的入射动能——这是"阈值"。如果你的束流能量低于这个值——比如 0.8-1.2 GeV/nucleon——单个核子-核子碰撞的能量不够造 K⁺。但如果很多核子同时参与——三体甚至多体碰撞——它们的能量可以累加起来，造出这个 K⁺。这叫亚阈产生。

亚阈产生是密度计。多体碰撞的概率正比于密度的二次方甚至三次方。压缩越狠——状态方程越软——密度越高，多体碰撞越多，亚阈 K⁺ 越多。这就是为什么 1990 年代 KaoS 实验的 K⁺ 产额指向了一个较软的状态方程。

HADES 探测器在 GSI 把这个逻辑推到了新高度。2024 年，HADES 完成了 Au+Au 在 0.8 GeV/nucleon 的束流实验——深亚阈，入射能量只有自由 NN 阈值的一半。同时，Ag+Ag 在 1.58 GeV/nucleon 的数据分析产出了第一批结果：K⁺ 介子的横向流 v₁ 是负的（被推离碰撞轴），表明 K⁺ 感受到一个排斥性的平均场势。Λ 超子的横向流也被提取出来了——这是这个能区 Λ 流的首次测量——对超子-核子相互作用势的约束，直接和"超子之谜"（中子星里为什么能有 Λ 超子而核力理论不能解释）相关。

最漂亮的 HADES 2024-2025 结果是关于 φ 介子的。φ 是隐藏奇异性——由一对奇异-反奇异夸克组成——它的质量在真空中是 1020 MeV。HADES 在 Au+Au 1.23 GeV/nucleon 中测到了 φ/K⁻ 产额比的增强。标准的统计强子化模型无法解释——φ 的被压低是因为奇异性规范压低。但 φ 的增强可以用"碰撞展宽"解释：介质中的核子不停撞 φ，缩短了它的寿命，按理说应该减少它的产额——实际上是撞的频率太快，φ 还没来得及衰变成 K⁺K⁻ 就被再生的 K⁺ 和 K⁻ 重新"凑"出来了。这个效应的前提是手征对称性在介质中部分恢复——φ 的宽度变宽了。如果它能被定量验证，这会是手征对称性恢复的间接信号。

再说全球装置全景。

2025 年，重离子碰撞物理有四台主力装置在运行或即将运行：

STAR 在 RHIC 上的 BES-II 完成了数据采集，正在做全统计量的最终分析。最精确的涨落数据、流数据、奇异粒子数据都已经或正在发表。固定靶模式延伸到 √s=3.0 GeV，覆盖了 μB≈700 MeV 的区域——这是全球唯一能在这个 μB 区做对撞物理的装置。

HADES 在 GSI 继续研究 0.8-2.0 GeV/nucleon 区。它的核心优势是双轻子探测——测量虚光子衰变成正负电子对的产额——直接探测火球内部的温度，不经过强子末态相互作用的污染。2024-2025 年的重点是可升级的 FAIR 阶段零物理——在 CBM 建成前，HADES 是唯一能在高重子密度区做双轻子测量的探测器。

NA61/SHINE 在 CERN SPS 完成了二维扫描（束流能量 × 核大小），2024-2025 年的重点是 150A GeV Pb+Pb 碰撞中粲强子的产物——这是 SPS 能区粲粒子的首次直接测量。SPS 的 17 GeV 能区正好是"起始解禁闭"的区域——K⁺/π⁺ 的"角峰"指出色禁闭在这个能量附近破裂。粲粒子——因为质量大、产生得早——是这个能区最干净的 QGP 探针之一。

NICA/MPD 在 JINR（杜布纳）在 2025 年进入了关键的一年。超导螺线管冷却到了 4.5K 并保持稳定。MPD 探测器的 TOF 系统——300 个 MRPC 探测器——全部组装完毕并测试通过。2025 年 8-9 月预期在 NICA 对撞环中第一次注入并循环离子束流（氙-124，B=1T）。第一个物理实验——Xe+Bi 碰撞在 √s≈2.5 GeV——预计在 2026 年开始。NICA 填补的是 √s=4-11 GeV 的对撞区——正好是 BES-II 覆盖不到的 μB 适中区。

CBM 在 FAIR 仍然在建设中。2024-2025 年，微顶点探测器的最终原型传感器（MIMOSIS-2.1）在多次束流测试中验证通过。STS 硅追踪器的读出芯片——特制的 Xyter ASIC——80 片晶圆完成批量测试。mCBM 原型在 SIS18 上持续跑束流，验证了"自由流式"无硬件触发读出——每秒 10 MHz 的反应率、每个事件自触发、数据流直写硬盘——在真实束流条件下可行。

中国在重离子碰撞物理中主要角色是理论和 STAR 实验参与。华中师范大学是 STAR 合作组中最大的非美国单位之一。2025 年，施梳苏课题组在 PRL 上发表了一篇重要论文：在 √s≤3.2 GeV 的 Au+Au 碰撞中，强子椭圆流的"组分夸克数标度"完全破裂。组分夸克数标度——把一个核子的 v₂ 除以 3（三个夸克），一个介子的 v₂ 除以 2，发现所有强子的 v₂/nq 落在同一条曲线上——是夸克自由度参与集体流动的经典信号。低于 3.2 GeV，它碎了。在 3.2-4.5 GeV 之间，它开始恢复。这个"恢复阈值"~3.2 GeV，直接对应了核物质从强子自由度主导到夸克自由度主导的过渡。

同年在华中师大成立了华中核理论中心，研究方向包括相对论重离子碰撞、流体力学、硬探针、机器学习在 QCD 中的应用。

再谈输运理论。它是在"第一性原理"和"实验"之间最重要的计算桥梁。

2024-2025 年，输运模型评估项目（TMEP）由 Wolter 等人牵头，系统比较了 BUU 族（BUU、GiBUU、PHSD、SMASH）和 QMD 族（IQMD、UrQMD、AMD、PHQMD）的代码。在完全相同的物理设定下——同样的状态方程、同样的截面、同样的初始化——两种代码对核子观测量仍有不可忽略的差异，π 介子产额的差异可以达到因子 1.5。这些差异不是 bug。它们反映了两种方法在底层假设的不同——BUU 用测试粒子代表分布的平滑演化，QMD 用高斯波包代表全 n 体关联——而这种不同在开放系统（碰撞）中被放大了。

最大的技术突破在耦合输运+流体力学混合模型。SMASH+vHLLE 框架在 2024-2025 年完成了几个关键步骤。动态初始化——在 SMASH 模拟中实时监测局部能量密度，一旦超过 0.5 GeV/fm³ 就切换到 vHLLE 粘滞流体力学——解决了低能区（√s<7.7 GeV）标准"离子穿过时间 τ₀"做初始条件的尴尬。首次在有限密度下引入三个守恒荷（B、Q、S）的流体力学——用 NEOS-4D 的格点 QCD 状态方程——实现了从 SMASH 初始化到 Cooper-Frye 粒子化，再到 SMASH 强子后散射的全链条。2025 年 PRC 112 上的贝叶斯推断用了这个混合框架，发现流体的剪切粘滞系数比之前用标准 Glauber 初始条件推出的值有显著偏移——初始条件对粘滞系数的提取影响巨大。

输运理论的另一个活跃前沿是机器学习。高斯过程模拟器——用少量训练好的输运计算来内插和外推整个参数空间——让贝叶斯推断可处理原本因计算量太大而放弃的多参数问题。Wang 等人的非弹性系数推断、Tsang 等人的对称能五参数同时约束，都依赖这项技术。

最后说中国自己的加速器。2025 年，兰州的 HIRFL-CSR（重离子加速器冷却储存环）在外靶实验上持续产出中能重离子碰撞数据。近物所张亚鹏课题组重点研究超核集体流——测量 Λ 超子的 v₁ 和 v₂——来约束超子-核子的介质内势。华南理工大学的合作团队在 2025 年发表了用 AMPT 模型做的超核产生与集体流研究。HIAF（强流重离子加速器，惠州）2025 年 10 月全线贯通——虽然第一阶段的重点是核结构和原子物理，但它的高温高密核物质研究——用高强度的重离子束流系统地扫描高重子密度区——是第二阶段的核心科学目标之一。

物理的线条拉直了看是这样。三十年前，我们用重离子碰撞"发现"了集体流。二十年前，我们发现指向流在某个能量处有一个坑——它是相变的指纹吗？十年前，BES-I 告诉我们可能是——统计不够，但方向对。2025 年，BES-II 的数据精度比 BES-I 高 4-5 倍，涨落分析达到了 4σ 的提示水平。状态方程、对称能、输运系数正在被贝叶斯方法框进越来越窄的置信区间。中子星和重离子碰撞——天和地——开始用同一个状态方程说话。

但最核心的问题——QCD 临界点在哪里，甚至它是否存在——仍然悬着。临界点可能就在 7.7 GeV 以下——在 μB>420 MeV 的深处，正好是 BES-II 固定靶模式和 NICA/FAIR/CBM 的领地。能不能抓到它，取决于未来三年的装置、统计量、和理论模拟器的精度。一旦抓到——如果抓到——它会是 21 世纪核物理最重大的发现之一。

参考文献

STAR Collaboration, Phys. Rev. Lett. 135, 142301 (2025) — Precision measurement of (net-)proton number fluctuations in Au+Au at BES-II: C₄/C₂ minimum at 19.6 GeV, 3.2-4.7σ deviation from non-CP baselines (BNL/STAR)
STAR Collaboration, Phys. Rev. Lett. (2025) — Emergence of number-of-constituent-quark scaling of hadron elliptic flow at √s>3.2 GeV (施梳苏/华师 STAR group)
E. Duckworth et al. (STAR Collaboration), arXiv: — Proton directed flow BES-II: dv₁/dy minimum at 11.5-19.6 GeV, transport v₁ excess breakdown, Sept. 2024 (BNL/STAR)
M.D. Cozma, Phys. Rev. C 110, 064911 (2024) — EOS from v₁+v₂+stopping: dcQMD+multivariate fit to FOPI data, K₀=230±10 MeV, L=63±10 MeV, in-medium σNN~0.4×vacuum (IFIN-HH Bucharest)
L. Ciampi et al. (INDRA-FAZIA), Phys. Lett. B 868, 139815 (2025) — Symmetry energy S=29.0±0.7 MeV at ρ/ρ₀≈1.01 from BUU+isospin transport ratio (INDRA-FAZIA/GANIL/LPC Caen)
M.B. Tsang et al. (SπRIT Collaboration), Phys. Lett. B 853, 138662 (2024) — Bayesian analysis of v₁+v₂+stopping in Sn+Sn 270 MeV/u: S₀, L, effective mass splitting, in-medium σNN (SπRIT/RIBF/FRIB)
M. Kurata-Nishimura, P. Pawłowski et al., Phys. Lett. B 871, (2025) — Directed and elliptic flow of light clusters in neutron-rich and -deficient Sn+Sn at 270 MeV/u: AMD comparison (RIBF RIKEN)
Y. Wang et al., arXiv:2406.07051 (2024) — Bayesian inference of K₀ from proton v₂ in Au+Au 400 MeV/n: K₀=191 MeV (soft, MD) vs 235 MeV (no MD) (IQMD+GP)
V. Kireyeu et al., arXiv:2411.04969 (2024) — PHQMD: soft-momentum-dependent EOS best reproduces HADES+FOPI v₁+v₂; cluster A-scaling of v₂ (JINR/GSI)
L.A. Tarasovicova et al., SMASH+vHLLE, Phys. Rev. C 112, 014910 (2025) — First Bayesian inference of hybrid transport+hydro model for Au+Au at 0.4-1.23 AGeV (Univ. Frankfurt/GSI)
HADES Collaboration (J. Orliński et al.), EPJ Web Conf. 316, 03012 (2025) — Charged K meson transverse flow v₁+v₂ in Ag+Ag at 1.58 AGeV; Λ v₁ first measurement at this energy (Univ. of Warsaw/HADES)
HADES Collaboration (N. Schild et al.), Nucl. Phys. A 1060, (2025) — Dilepton observables: fireball temperature, radial and anisotropic flow, groundwork for CBM (GSI/HADES)
HADES Collaboration, Phys. Lett. B 867, 139605 (2025) — Simultaneous description of high-density QCD matter from HADES data and neutron star observations via Chiral Mean Field Model+UrQMD
H. Wolter et al. (TMEP), Phys. Rev. C 109, 044609 (2024) — BUU vs QMD code comparison: pion production at 270 A MeV, persistent differences, Bayesian model averaging proposed (Univ. of Munich/TMEP)
J. Aichelin, NUSYM 2024 — "How to explore the hadronic equation of state": BUU vs QMD foundations review, momentum-dependent potentials, clustering capability (Subatech/Nantes)
SMASH Collaboration, Universe 11, 27 (2025) — Elliptic and quadrangular flow v₂+v₄ of protons in high-μB region: v₄/v₂²≈0.5, mean-field essential (SMASH)
Z. Paulínyová et al., Quark Matter 2025 — SMASH+vHLLE dynamical initialization for low √s: continuous energy-density threshold, fluidization at 0.5 GeV/fm³ (Univ. Frankfurt/GSI)
CBM Collaboration (T. Tollefson et al.), arXiv:2506.20545 (2025) — CBM detectors and electronics: free-streaming DAQ, MVD/MIMOSIS-2.1, STS/Xyter wafer testing, mCBM prototyping status (FAIR/GSI)
NICA/MPD Collaboration, 15th MPD Meeting (April 2025) — MPD TOF 100% complete (300 MRPC), SC solenoid at 4.5K, first Xe beam Aug-Sep 2025, full commissioning mid-2026 (JINR Dubna)
NA61/SHINE Collaboration, CERN SPSC-SR-368 (2025) — 2024-2025: Pb+Pb 150A GeV open charm, 2D scan complete, fragmentation for cosmic-ray physics (CERN SPS)
S. Lan et al., Chin. Phys. C (2025) — AMPT: energy dependence of v₂{PP}/v₂{RP} ratio in Au+Au 7.7-200 GeV: saturation above ~39 GeV (CCNU/SINAP)
Y. Zhang et al. (CCNU/STAR), Phys. Rev. Lett. (2025) — NCQ scaling emergence at 3.2-4.5 GeV: transition from hadronic to partonic collectivity (Central China Normal Univ./STAR)
C.M. Ko, NUSYM 2024 — Review of sub-threshold kaon production and its constraint on EOS: soft EOS from KaoS, extension to multi-strange and symmetry energy (Texas A&M)
Song, Aichelin, Bratkovskaya, PoS QNP2024 (2025) — Off-shell PHSD: simultaneous reproduction of KaoS+FOPI+HADES (anti)kaon data with in-medium kaon potentials and φ broadening (Subatech/FIAS)
J. Xu et al. (TMEP), Phys. Rev. C 109, 044609 (2024) — BUU vs QMD for pion/Δ production: code convergence and uncertainty propagation for symmetry energy extraction
H. Elfner et al., EPJ Special Topics (2025) — SMASH transport mini-review: collective flow, dilepton v₂, critical fluctuation survival through hadronic stage (GSI/Goethe Univ. Frankfurt)
R. Wittmann, U. Heinz, arXiv:hep-ph/9509328 (updated 2024) — Subthreshold antiproton spectra: meson-baryon processes dominant, T_eff≈T/2 (Univ. Regensburg/Ohio State)
STAR Collaboration, J. Subatomic Part. Cosmol. 3, 100053 (2025) — Proton cumulants from hydrodynamics: BES-II factorial cumulants, CP location constraints (Vovchenko et al.)
近物所/华南理工, 2025 — 超核产生与集体流: Λ v₁+v₂ from AMPT, constraint on in-medium hyperon-nucleon interaction (IMP/SCUT)