胶子饱和与色玻璃凝聚

如果你拿一把能量越来越高的"电子显微镜"去看一个质子，你看到的不是一个越来越清晰的三夸克小球。你看到的是越来越多、越来越细碎的小点——胶子。它们的密度随能量指数增长，在某一刻——能量够高、x够小——不再增加。不是胶子停止了增殖。是新胶子被已有胶子"吸收"了——一个胶子和另一个胶子合并成一个。这就是胶子饱和。描述这个状态的理论，叫色玻璃凝聚（CGC）。

二十年来，CGC是QCD最没被实验验明正身的漂亮框架。它有一个优美到不真实的核心预言——几何标度：所有小x物理截面，不管你独立变化Q²和x，只依赖一个组合变量τ=Q²/Q_s²(x)。Q_s是饱和动量——胶子密度的"天花板"。HERA的H1和ZEUS数据在1990年代末就显示了这个标度，但没人能说它是饱和的独特信号——线性偶极子模型也能拟合同样的数据。

2024-2026年，这个领域翻到了新的一章。

理论端，CGC终于从领头阶精度走到了次领头阶。这不是小升级。领头阶的CGC只是一个定性框架——它告诉你"质子在高能下像一堵胶子墙"，但它算不出墙上具体有多少胶子、墙有多密实。次领头阶的BK方程（Balitsky-Kovchegov演化方程）2008年就推导出来了，但此后十四年没人能稳定地求解它——方程里的双对数项让散射截面在数值上变成负数。一套话来说是"共线不稳定性"——在Q²远大于Q_s²的过渡区，NLO修正出现了不物理的振荡。

2025年，三个独立方案几乎同时解了这个问题。捷克技术大学的Cepila等人给出了首个碰撞参数依赖的NLO BK稳定数值解——这在技术上意味着可以处理质子"边缘"（大b）vs"中心"（小b）的饱和差异了。BNL的Boussarie、Caucal和Mehtar-Tani从第一性原理出发，在CGC算符空间里做了一个基变换——用"复合偶极子算符"替换原来的BK振幅，干脆利落地消掉了所有共线双对数。芬兰Jyvaskyla大学的Casuga、Hänninen和Mäntysaari把NLO BK方程和NLO DIS冲击因子结合起来，做了首个贝叶斯全局分析——不只解方程，而是把所有参数的不确定性用一个概率分布描述出来。他们从HERA数据提取的初始条件中，反常维数γ≈1——比此前需要γ≈1.1-1.2的拟合更"自然"。

与此同时，同阶层的JIMWLK重整化群——比BK更完整的、不依赖大Nc近似的基本演化方程——在2024年遭遇了一场概念地震。Ben-Gurion大学的Kovner、NCSU的Lublinsky、Skokov和Zhao在JHEP上发了一篇论文，标题就说明了一切："不是所有β₀都来自β函数——NLO JIMWLK中的DGLAP重求和与跑动耦合"。简单说：过去二十年里，NLO方程中那些跟β₀有关的对数项，被一致地当成了跑动耦合修正。但这篇论文证明，其中相当一部分根本不是跑动耦合，而是弹丸波函数的DGLAP演化——一个完全不同的物理过程。这个发现动摇了传统跑动耦合方案的理论根基。

理论精度在飞跃，但实验信号呢？

最诚实的回答来自ERC 2024年启动的GlueSatLight项目（居韦斯屈莱大学，199万欧元），它的项目摘要第一句就说："迄今为止尚未观察到胶子饱和的明确信号。"

不是没有信号。是这些信号都有替代解释。

HERA小x的F₂数据呈现完美的几何标度——τ=Q²/Q_s²把所有数据点压缩到一条曲线上。CGC预言了这个标度。但线性BFKL演化（不含饱和）也可产生类似行为——它只是强耦合跑动导致的巧合式标度。区别在哪？CGC说标度来自非线性（胶子重组的截断），BFKL说标度来自线性（对数重求和）。两种机制产生几乎一样的观测量——这是二十年来这领域的核心困境。

RHIC和LHC的前向快度强子压低——d+Au或p+Pb碰撞中，快度大于2的区域强子产额比标度后的pp少30%-70%——一度被认为是饱和的第二个有力证据。CGC说，前向快度探针的是靶核中小于10⁻⁴的x区间，那里的胶子密度饱和了——少产强子是因为部分子"不够"了。但2025年，Perepelitsa在PRC上做了一个关键演示：现代核PDF组（EPPS21、nNNPDF）——不含非线性饱和、只用线性DGLAP演化——能完全描述前向双强子压低的数据。饱和不是必须的。nPDF同样可以产生压低——它只是用不同的参数化方式吸收了"胶子不够"的效应。

不过，有一组数据开始倾斜天平。LHC的铅-铅超周边碰撞（UPC）中的遍举J/ψ光产生，在Royon 2025年的分析中明确偏好BK演化（有饱和）而非BFKL（纯线性）。光子打到一个铅核上产生J/ψ的过程，截面正比于小x胶子密度的平方——它对饱和的敏感度比单举观测量高。γ+p靶中BK和BFKL区别不大——质子饱和标度太低了。但γ+Pb靶中区别变得显著——铅的饱和标度比质子大了约A^{1/3}≈5.5倍，进入了非线性区域。

欧洲Jyvaskyla的Mäntysaari团队在2025年PRL上把JIMWLK重整化群方程直接集成到了IP-Glasma初态模型中。IP-Glasma是重离子碰撞初态的标准框架——它用CGC有效场论计算碰撞后一瞬间（τ≈0.2 fm/c）的能量-动量张量，喂给后续的粘滞流体动力学。此前IP-Glasma里的饱和动量Q_s是参数化的——从HERA数据拟合一个函数，再往碰撞能区外推。现在用JIMWLK直接演化——不需要参数化了。关键发现：RHIC能区（200 GeV）两者几乎没区别；但LHC能区（5.02 TeV）JIMWLK让核子几何更平滑、压力梯度更小、平均横动量更低——小系统p+Pb、O+O的效应最显著。

O-O碰撞——2025年LHC首次运行的轻离子碰撞——是系统尺寸扫描的关键棋子。从pp到pPb到O-O到Pb-Pb，饱和标度Q_s²≈A^{1/3}连续增长。O核的质量数是16，铅是208——中间隔着近两个数量级的核厚度。Cepila等人在2025年计算了O-O和Ne-Ne的UPC矢量介子光产生，预言非相干截面是一个比单举测量更清晰的饱和信号。

说到小系统——pp和pPb碰撞中的"ridge"（长程方位角关联）——这是当前最尖锐的张力点。Carrington和Mrówczyński在2024年的Acta Physica Polonica上做了一个关键的解析计算：他们用proper time展开计算了Glasma（碰撞后经典胶子场阶段，在热化之前）的能-动张量T^{μν}，发现横向Poynting矢量T^{0i}几乎完美地跟随能量密度梯度——跟流体动力学流的行为一模一样。但这纯粹是经典胶子场的效应，跟热化、跟流体动力学无关。空间几何的偏心率被经典场直接转化成了动量空间的各向异性。他们的结论："Glasma模拟了流体动力学行为。"

这意味着什么？你在pp碰撞中看到粒子按方位角不对称地出来（v₂>0），你不一定需要QGP流体动力学——经典Yang-Mills场自己就能做到。但ALICE 2025年3月的预印本又测到了|Δη|>6.5的超长程关联——纯CGC初态关联在快度上退关联很快（|Δη|≈2-3就没了），解释不了超长程。所以当前共识是"混合"：CGC初态产生短程关联，流体动力学或部分热化拉长长程关联。

Grosse-Oetringhaus和Wiedemann在2024年7月发了一篇81页的十年综述，结尾是："十年来没有任何唯象框架令人满意地描述所有系统尺寸下的集体性。"这不是宣告CGC失败——是在说"初态+末态"的拼图还缺关键的一块。

---

现在说中国。

华中师范大学是中国CGC研究的核心枢纽。2025年5月新成立的华中核理论中心（C3NT）计划在2027年做一个大型国际项目——"精密饱和物理时代"。王新年在这个领域做了一件漂亮的活儿。

2025年，王新年、华南师范大学邢宏喜、UCLA康忠波和Farid Salazar在PRL和PRD上同时发表了两篇论文——CGC有效场论和高扭度共线因子化理论在twist-4层次上的一致性证明。这两个框架分别来自完全不同的理论传统：CGC源自小x重求和的非线性动力学（强场极限），高扭度共线因子化源自大Q²的算符乘积展开（稀疏场极限）。它们应该在某个中间区域重合——三十年来大家都这么说，但从来没人证明过。王新年他们在直接光子产生中做了这件事——超冲击波近似、LPM干涉效应、sub-eikonal相位全部包含——然后展示了两个框架在twist-4层次的对应关系。这是一个"统一场论"式的成果：从稀疏胶子气到凝聚态胶子墙的完全相图。

肖博文不在华中师大了。他2020年调到了香港中文大学（深圳）。但他的CGC框架下的Sudakov重求和——2013年与Mueller和Yuan合作的PRL——是过去十年被引用最频繁的CGC基础理论工作之一。Caucal和Salazar 2024年那篇NLO双强子TMD因子化的PRL，核心引用之一就是肖博文2013年的Sudakov方案。肖博文团队2024年还在EPJ Web of Conferences上发表了CGC NLO前向强子产生的阈值重求和——解决了NLO出现负截面的长期问题，让RHIC和LHC的所有相关数据可以在一个统一的公式下描述。

山东大学的周剑在做CGC和TMD的衔接。他在2024年的JHEP上重新分析了衍射双喷注的方位角不对称度。他更早（2017年）就和肖博文、袁烽合作建立了Collins-Soper（TMD）+BK（小x）联合重求和的框架——这是现在Caucal和Iancu在PRD 2025上那篇JIMWLK+DGLAP+CSS三种演化统一框架的先驱之一。

中国科学技术大学的唐泽波、查王妹团队在做EicC的饱和探测方案设计——碰撞参数操控、中子标记、胶子层析。他们2025年在Nuclear Science and Techniques上发表了碰撞参数操控的详细方案。EicC还在概念设计阶段，比美国EIC晚一个身位——但作为EIC的互补设施，它有自己的独特能区覆盖。

差距也是明显的。JIMWLK数值模拟、IP-Glasma完整事件模拟——这两个方向的中国团队几乎空白。Schenke（BNL）和Mäntysaari（Jyvaskyla）是IP-Glasma的核心开发者；Lappi是JIMWLK数值求解的领军者。中国没有自研的对标代码。华师大的CLVisc是3+1D流体力学模型——但它不做CGC初态，初态还是从IP-Glasma导入的。

人才流失也是问题。肖博文从华师大调到港中深，廖劲峰（他不在华师大——更正一下——在印第安纳大学，而且2024-2025年主攻手征磁效应而非CGC），CGC中国核心理论人才处于散落状态。华师大的C3NT正在试图重新聚合这个群体——2025年7月办了国内首届LHC小x物理研讨会，理论家和实验家（CMS/LHCb/ATLAS/ALICE）第一次聚在一起讨论饱和物理。

---

2024到2026年，胶子饱和研究的进展可以总结为三件事。

理论精度从"定性对"走向"定量对"——NLO BK的共线不稳定性被两个独立方案解决，JIMWLK+DGLAP+CSS三种演化的统一框架建立，衍射DIS的完整NLO计算完成，三圈（NNLO）BK哈密顿量的非共形部分已给出。这是CGC从领头阶到次领头阶的成人礼。

但实验仍然没有"smoking gun"。HERA几何标度可以被线性BFKL解释，LHC前向压低可以被nPDF解释。唯一开始倾斜天平的是UPC遍举矢量介子光产生——铅靶中BK明确优于BFKL。EIC——2024-2025年的理论准备显示——F_L（纵结构函数）和F_2^{cĉ}（粲结构函数）是两个黄金测量道，信号差异达到20%以上，远超EIC的预期系统误差。但EIC最早也要到2030年代中期才能出数据。

CGC作为初态理论在重离子碰撞中的地位反而在上升。IP-Glasma+JIMWLK让初态变成真正无参数的第一性原理计算；Glasma阶段被证明可以"模拟"流体动力学行为——这对"流=QGP"的简单标识提出了根本挑战；KøMPøST-3D建立了从CGC初态经预平衡到3+1D流体动力学的完整链路。换句话说，CGC作为"初态产生器"的地位比作为"新物质态探针"的地位站得更稳。

与此同时，最根本的问题——胶子饱和在质子/核里到底存不存在——二十年了还摆在那里。但理论工具和实验设计已经走到了一道分界线之前。EIC将第一次给我们一个足够大的Q_s（核靶的A^{1/3}倍增强）、足够清洁的初态（光子代替强子）、和足够的精度（3%系统误差）去面对这个问题。如果连EIC都分不出饱和和非线性之外的解释——那就算看到了，也算没看到。

---

参考文献

[1] Cepila J., Contreras J.G., Matas M., Vaculciak M., "Collinearly-improved BK equation at NLO with impact parameter dependence," *Physical Review D* 111, 096015 (2025); arXiv:2412.08571. (Czech Technical University in Prague)

[2] Casuga C., Hänninen H., Mäntysaari H., "Confronting CGC at NLO with HERA data using Bayesian inference," *Physical Review D* 112, 034003 (2025) — Editors' Suggestion; arXiv:2604.22332. (University of Jyväskylä)

[3] Boussarie R., Caucal P., Mehtar-Tani Y., "Fixing small-x: first-principles resolution of collinear instabilities in NLO BK evolution," arXiv:2509.20236 (2025). (Brookhaven National Laboratory)

[4] Brunello F., Caron-Huot S., Crisanti G., Giroux M., Smith D., "Three-loop BK Hamiltonian from space-like to time-like correspondence," *Journal of High Energy Physics* 11, 055 (2025); arXiv:2508.03794. (McGill University)

[5] Kovner A., Lublinsky M., Skokov V., Zhao Z., "Not all that is β₀ is β-function: the DGLAP resummation and the running coupling in NLO JIMWLK," *Journal of High Energy Physics* 07, 148 (2024); arXiv:2308.15545. (Ben-Gurion University / North Carolina State University)

[6] Altinoluk T., Beuf G., Lublinsky M., Skokov V., "On running coupling in the JIMWLK evolution and its Langevin formulation," *Journal of High Energy Physics* 03, 131 (2024). (NCBJ Warsaw / Ben-Gurion University / NCSU)

[7] Beuf G., Lappi T., Mäntysaari H., Paatelainen R., Penttala J., "Diffractive deep inelastic scattering at NLO in the dipole picture," *Journal of High Energy Physics* 2024, 024 (2024); arXiv:2401.17251. (University of Jyväskylä / NCBJ Warsaw)

[8] Caucal P., Iancu E., "Evolution of the transverse-momentum dependent gluon distribution at small x: JIMWLK + DGLAP + CSS unified framework," *Physical Review D* 111, 074008 (2025); arXiv:2406.04238. (Université Paris-Saclay / CEA Saclay)

[9] Mäntysaari H., Schenke B., Shen C., Zhao W., "IP-Glasma with JIMWLK evolution: from RHIC to LHC," *Physical Review Letters* 135, 022302 (2025); arXiv:2502.05138. (University of Jyväskylä / Brookhaven National Laboratory)

[10] Carrington M.E., Mrówczyński S., "The earliest phase of relativistic heavy-ion collisions: the Glasma," *Acta Physica Polonica B* 55, 4-A3 (2024); arXiv:2401.13705. (Brandon University / Jan Kochanowski University)

[11] Grosse-Oetringhaus J.F., Wiedemann U.A., "A decade of collectivity in small systems," arXiv:2407.07484 (2024). (CERN)

[12] Perepelitsa D.V., "Description of di-hadron saturation signals within a universal nuclear parton distribution function approach," *Physical Review C* 111, 054901 (2025); arXiv:2501.18347. (University of Colorado Boulder)

[13] Garrido J., Guevara A., Sanhueza M., "Confronting NLO CGC/Saturation approach with combined HERA data," *PoS DIS2024* 092 (2025); arXiv:2412.15234. (Universidad Técnica Federico Santa María, Chile)

[14] Fu Y., Kang Z.-B., Salazar F., Wang X.-N., Xing H., "Correspondence between Color Glass Condensate and High-Twist formalism," *Physical Review Letters* 135, 032301 (2025); *Physical Review D* 112, 014029 (2025); arXiv:2406.01684. (Central China Normal University / South China Normal University / UCLA)

[15] Shi Y., Wang L., Wei S.-Y., Xiao B.-W., "Pursuing the precision study for CGC in forward single hadron productions," *EPJ Web of Conferences* 296, 10008 (2024). (CUHK-Shenzhen / Shandong University / CCNU)

[16] Tong X.-B., Xiao B.-W., Zhang Y.-Y., "Harmonics of parton saturation in lepton-jet correlation at EIC," *Physical Review D* 109, 054004 (2024). (CUHK-Shenzhen)

[17] Xing H. et al., "eHIJING: an event generator for jet tomography in electron-ion collisions," *Physical Review D* (2024). (South China Normal University)

[18] Tang Z.-B., Zha W.-M., Wu X., Li X.-B., Wang K.-Y., "Impact parameter manipulation in exclusive photoproduction in Electron-Ion Collisions," *Nuclear Science and Techniques* 36 (2025). (University of Science and Technology of China)

[19] Liu H.-Y., Liu X.-H., Pan J.-C., Yuan F., Zhu H.-X., "Nucleon Energy Correlators for the Color Glass Condensate," arXiv:2301.01788 (2023-2024). (Beijing University of Chemical Technology / Beijing Normal University / IHEP / LBNL / Zhejiang University)

[20] Royon C., "UPC J/ψ production at the LHC: BK vs BFKL," *Acta Physica Polonica B Proceedings Supplement* 18, 1-A48 (2025). (University of Kansas)

[21] Hauksson S., Iancu E., Mueller A.H., Triantafyllopoulos D.N., Wei S.-Y., "TMD factorization for diffractive jets in DIS at small x," *Journal of High Energy Physics* 05, 024 (2024). (CEA Saclay / Columbia / ECT* / Shandong University)

[22] Du X., Schlichting S., Zhu Y., "KøMPøST-3D: pre-equilibrium transport in full 3+1D," *Physical Review D* 113, 096009 (2026); arXiv:2512.08007. (Universität Bielefeld)

[23] Ipp A., Leuthner M., Müller D., Schlichting S., Schmidt P., Singh A., "Energy-momentum tensor of the dilute (3+1)D Glasma," *Physical Review D* 109, 094040 (2024); arXiv:2401.10320. (TU Wien / Universität Bielefeld)

[24] Boussarie R., Fucilla M., Szymanowski L., Wallon S., "Probing gluonic saturation in deeply virtual meson production beyond leading power," *Physical Review Letters* (2024); arXiv:2407.18203. (BNL / NCBJ Warsaw / IJCLab Orsay)

[25] LHCb Collaboration, "Nuclear modification factor of D⁰ mesons in pPb collisions at √s_NN = 8.16 TeV," *Physical Review Letters* 131, 102301 (2023). (CERN)

[26] ALICE Collaboration, "First observation of ultra-long-range correlations in pp and pPb collisions," CERN-EP-2025-080 (2025). (CERN)

[27] Caucal P., Salazar F., "NLO semi-inclusive dihadron production in the CGC," *Journal of High Energy Physics* 2024, 130 (2024); arXiv:2401.01934. (BNL / UCLA)

[28] Lu Z., Tang Z.-B., Wu X., Zha W.-M., "Bayesian reweighting of nuclear PDFs with coherent J/ψ photoproduction data at RHIC and LHC," arXiv:2502.09063 (2025). (USTC)

[29] GlueSatLight ERC Project, Mäntysaari H. et al., University of Jyväskylä, 2024-2029. ERC Consolidator Grant, €1.99M.

[30] Boer D., Hagiwara Y., Zhou J., Zhou Y., "Energy evolution of T-odd gluon TMDs at small x," arXiv:2203.00267 (updated 2024). (Shandong University)