原子核团簇与α凝聚-从霍伊尔态到中子星
弗雷德·霍伊尔在1953年走进加州理工的凯洛格实验室,对一个叫沃德·怀廷的实验物理学家说:碳-12的核里一定有一个激发态,在7.68兆电子伏左右。如果你不打到这个态,宇宙里的碳就不够。不够造恒星,不够造行星,不够造你和我。
怀廷觉得这个英国人在胡扯——一个天文学家跑来告诉核物理学家原子核能级应该在哪?但他还是去做了。1953年9月,他们在物理评论上宣布找到了这个态:实际位置7.65 MeV,与霍伊尔的预测差了0.03 MeV[1]。
这件事后来被物理学史学家反复争论。霍伊尔是不是用了"人择原理"?他是不是说"因为碳基生命存在,所以碳-12必然有这个共振态"?历史学家Helge Kragh的研究给出了更平实的答案:霍伊尔只是在解一个具体的恒星核合成问题[2]。三个α粒子(氦-4核)需要同时碰到一起才能合成碳-12,但中间体⁸Be寿命只有6.7×10⁻¹⁷秒——比三体碰头的窗口窄太多。如果不共振,三α过程的反应率不足。碳不会生成。恒星会在氦燃烧阶段卡住。但如果碳-12恰好有一个0⁺激发态,能量刚好在三个α粒子的总静止能量之上一点点,α粒子就可以通过共振隧道被捕获——反应率飙升一千万到一亿倍。
这个7.65 MeV的态——今天被称为霍伊尔态——在当时只是核天体物理的一个聪明答案。没有人意识到它提出了一个更深的问题:什么样的核结构能让三个α粒子以这么低的总能量聚在一堆,却不塌成一个紧密的球?
答案等了将近五十年。
2001年,京都大学的藤崎昭博(Akihiro Tohsaki)、堀内昶(Hisashi Horiuchi)、德国的Peter Schuck和Gerd Röpke——四个人,两个国家——发表了一种全新的波函数[3]。他们把它叫做THSR波函数(四个姓氏的首字母)。它的数学形式像一个高斯函数——就像你把几个α粒子放进同一个最低能量的轨道,像玻色-爱因斯坦凝聚(BEC)里的原子一样。公式写出来很简单:几个α粒子占据同一个0s轨道,轨道宽度由一个参数B控制。当B远大于每个α粒子本身的尺寸时,这些α粒子互不挤压,像一团稀薄的气体。
这跟传统的核团簇模型完全不同。传统模型——比如Brink的α团簇模型——把α粒子看作硬球,放在特定的几何位置上。藤崎他们不做这个假设。他们只要求α粒子在同一条轨道上,具体位置由量子力学自己决定。结果呢?这个THSR波函数和霍伊尔态的传统团簇波函数(经过完全微观计算的)之间的重叠——不是80%、不是90%——是接近100%。
这在量子多体问题中极为罕见。一个如此简单的解析波函数,几乎完美地捕获了一个极其复杂的关联态的本质。这意味着霍伊尔态的本质就是"α粒子的BEC"——三个α粒子凝聚在同一个量子态上[4]。
这个图像解释了好几个过去说不清的事。为什么霍伊尔态的体积比碳-12基态大那么多——基态的α粒子间距离大约2.5 fm,霍伊尔态膨胀到大约6.8 fm,体积扩大了近十倍。为什么壳模型算不出霍伊尔态的低能量——因为壳模型用的是单核子自由度,而霍伊尔态需要三个α粒子作为基本单元。霍伊尔态的BEC解释还连通了一条更大的物理线索:池田清美(Kiyomi Ikeda)在1968年提出的"阈值规则"——团簇态总是出现在粒子发射阈值附近。霍伊尔态恰好在3α阈值以上约300 keV。这不是巧合。在阈值附近,核系统在"束缚"和"解离"的边缘——这正是团簇自由度被激活的条件。
但直到2023-2026年,这个图像才经历了两次严酷的检验:一次来自第一性原理计算,一次来自实验直接观测。
2026年2月,东京大学和理研的大塚孝治(Takaharu Otsuka)、佛罗里达州立的Alexander Volya和大阪公立大学的板垣直之(Naoyuki Itagaki)在《欧洲物理杂志A》发表了一篇综述[5]——但不止是综述。他们用了蒙特卡罗壳模型(MCSM)加Daejeon16和JISP16这两种真实的核子-核子相互作用——从头算,没有任何团簇预设——去算¹²C的基态和激发态。α团簇自己出现了。不是人为放进去的,是算出来的。
这不止是团簇的存在得到了验证。大塚他们还发现了一个至今未在其他任何核上被同时观察到的现象:¹²C内部存在两种完全不同的旋转模式。一种是"致密物体旋转"——整个碳-12像一个压扁的橄榄球绕着轴转,这是基态转动带的模式。另一种是"分离物体旋转"——三个α粒子彼此不挤压,像一组三体星绕着公共的质心转,这是霍伊尔态转动带的模式。同一堆核子,两种转法。大塚把¹²C称为"在团簇描述和平均场描述等级跃迁边界上的核"——再往前一步,核子就彻底分家了。这也解释了为什么碳-12的基态实际上也混了一小部分(约0.8 MeV的关联能来自)α团簇成分——团簇不只存在于激发态,它在基态中也有痕迹,只是被压缩到了最小限度。
实验侧的检验来自另一条更激进的路线。如果α粒子可以凝聚,那为什么中子对不能凝聚?
2023年12月,北京大学叶沿林教授和杨再宏副教授领导的团队——包括博士生陈嘉豪、魏康和合作单位日本理研、北航、复旦、韩国IBS等——在《物理评论快报》上发表了一篇Editors' Suggestion论文[6]。他们在世界上最丰中子的稳定核——⁸He——中找到了一个奇特的0₂⁺激发态。⁸He只有两个质子和六个中子,处于中子滴线的边缘。传统壳模型说它的低激发态应该是单粒子激发模式。叶团队用理研RIBF的次级束流做敲出反应实验,重建了衰变产物的四动量。他们发现了一个0⁺态——不是壳模型预言的那种,而是展现出三个团簇特征:一个α粒子核心,加两个双中子对,彼此松散耦合。
三个证据:第一,自旋-宇称被独立确定为0⁺。第二,这个态的同位旋标量单极跃迁强度(ISMR)异常大——这是BEC态的标志信号,说明核子在集体地、相干地重新分布。第三,它发射出的中子不是随机的——是强关联的中子对。α+2n+2n微观模型计算与实验数据完美吻合[7]。
这是人类第一次在原子核中看到"中子对的凝聚"——不是α粒子凝聚,而是中子对(双中子)凝聚。在标准的BEC类比中,⁸He的0₂⁺态是α粒子凝聚向中子-rich系统的第一次跨越。叶沿林在2025年5月的INPC大会上报告了这个系列的后续:¹⁴C中初步发现了α+α+⁶He的凝聚态候选,这意味着BEC类比不仅适用轻核,可以在整个核素图向中子过剩方向推移[8]。
与此同时,霍伊尔态向更重核的推广也在加速。
2023年6月,叶沿林团队的陈嘉豪、魏康等人在《科学通报》上发表了¹⁶O中霍伊尔类结构的新证据[9]。他们在北京中国原子能科学研究院的HI-13串列加速器上,用96 MeV的¹⁶O束流轰击¹²C靶。八组硅条带电粒子望远镜组成了高精度阵列。在4α衰变道(¹⁶O* → 4α)中,他们在这个霍伊尔态衰变道上空看到了一组窄共振——四个态,挤在4α分离阈值(14.44 MeV)以上不到2 MeV的窗口里。每一个态优先衰变到¹²C的霍伊尔态(7.65 MeV)再加一个α——就像霍伊尔态是"三α凝聚",这些态是"四α凝聚"的候选者。
魏康、陈嘉豪和叶沿林在2024年的《原子核物理评论》和《核科学与技术》综述[10,11]中总结了这条实验线索:从¹²C的三α(霍伊尔态)→ ¹⁶O的四α → 再到¹⁴C的α+α+⁶He(突破了α共轭限制,迈入中子-rich)→ 甚至⁸He的α+2n+2n(中子对凝聚)。每一步都是把BEC团簇的版图向外推一格。关键实验瓶颈是多重符合探测——要同时抓到四个或五个带电粒子并准确重建它们的动量,这需要覆盖大立体角、高粒度、同时具备粒子鉴别能力的大型阵列。
理论端也在追。2025年1月,一个解析可解模型被用于系统计算¹²C、¹⁶O和²⁰Ne中的α凝聚态能量[12]。模型预测:¹²C中的霍伊尔态凝聚能符合实验值;¹⁶O中预测的凝聚态能量接近实验候选态位置;²⁰Ne中预测的5α凝聚能量大约在1 MeV左右——远低于目前实验看到的候选态(约3 MeV)。这意味着要么实验态不是纯5α凝聚(混了壳模型成分),要么模型的范德华修正不够。²⁰Ne的5α凝聚态——五个α粒子共用同一个0s轨道——还没被确认。
这引出一个贯穿始终的物理问题:凝聚态到底有多"纯"?
答案随着中子和质子的加入在变。在¹²C中,霍伊尔态中α粒子占据0s轨道的概率约为70%——不是100%,因为泡利不相容原理留下了30%的"被禁轨道"影响。在¹⁶O的候选4α态中,这个比例可能更低——因为多了一个α,Pauli压制更强。在¹⁴C和⁸He的中子-rich候选态中,"凝聚态"这个标签可能需要更松弛:中子可能不是以纯α粒子核心的形式参与凝聚,而是以配对中子——双中子——的形式形成玻色子,再与α粒子一起占据同一个模式。这是从一个概念(α凝聚)向一个更广的概念(核内玻色子凝聚)的延伸。
这个延伸最重要的是它指向的方向:中子星。
中子星内部——内壳层那大约一公里厚的区域——是大自然能提供的最大尺度的核物质实验室。根据标准图像,内壳由中子-rich的原子核(像浸在超流中子海里的岛屿)和均匀中子物质之间的复杂共存构成。"核面食"(nuclear pasta)——棒状、板状的极端形变核结构——就出现在内壳底部,密度接近核饱和密度的一半。
α团簇和BEC类比在这里扮演了一个意想不到的角色。内壳中那些嵌入中子海的原子核——主要是铁峰和更重的核——在极高中子密度下会变得比地球上更"松散"。中子海把核子波函数往外推,降低了α粒子发射阈值——就像碳-12中霍伊尔态靠近3α阈值一样。在足够高的中子密度下,核可能会完全分解成α粒子加自由中子的混合BEC态。2024-2025年的理论探索——包括Argonne的Lovato团队用神经量子态(NQS)研究稀薄核物质[13],以及印第安纳大学等组用量子蒙特卡罗研究中子物质中的团簇涌现[14]——都指向一个结论:在从核物质到中子物质的密度渡越区,α团簇的重组可以剧烈改变状态方程。
这不是一个边缘假设。如果内壳存在α凝聚成分,超流中子海的黏滞系数会变。脉冲星的自转突变——glitch——的恢复时间取决于内壳的超流涡旋动力学,而涡旋-核散射截面取决于核的尺寸和内部结构。如果核是固态小球,截面小,恢复慢。如果核是蓬松的、带α凝聚的松散态,截面大,恢复快。Vela脉冲星和Crab脉冲星的不同glitch行为,可能部分源自其内壳核团簇态的差异。
还有一个更远、更美的连接。中子星并合——2017年GW170817让人类第一次"听到"的事件——产生的千新星光变曲线取决于r过程核合成的输出,而r过程的起点取决于中子星潮汐瓦解时抛出的物质是以α粒子、重核还是自由中子为主。如果内壳在剥离时先分解成α粒子,后续核合成的路径会根本不同。这是目前天体物理模拟几乎完全没有考虑的因素。
2026年的核团簇物理处在一个辉煌又痛苦的位置。
辉煌:从霍伊尔1953年的天才预言,到THSR波函数揭示BEC本质,到MCSM从头算确认团簇涌现和双旋转模式,到⁸He/¹⁶O/¹⁴C实验证据链——七十三年间,一个解释"碳为什么存在"的问题变成了一个解释"核物质多种存在形式"的问题。霍伊尔态不再是核物理边缘的奇珍异宝——它是核物质相图上一个连续区域的原型。
痛苦:5α凝聚态还没找到。中子-rich凝聚态的纯净度还没定量。中子星内壳中的α凝聚分数到底是0.1%还是10%——没人知道。探测多重团簇的符合实验需要的探测器阵列还不够大、不够快。理论端,团簇模型和壳模型仍然各自为政——只有很少的框架(大塚的MCSM、CNCI)能在壳模型的基础上自然产生团簇。
北京大学叶沿林团队正在建新一代探测器阵列,瞄准14C和更重中子-rich体系的凝聚态指纹。日本的板垣直之和堀内昶正在把THSR框架推向更重的核[15]。大塚的MCSM正在攻击¹⁶O和²⁰Ne的从头算团簇结构。三个方向汇聚到同一个科学问题:核物质在"凝聚"和"独立粒子"之间到底是怎么过渡的——而这个过渡决定了中子星内部是一锅均匀汤还是一盘有料的面食。
霍伊尔如果活到今天,大概会说:我只要一个7.65 MeV的态,你们给我推出来一个宇宙。
---
参考文献
[1] F. Hoyle, D.N.F. Dunbar, W.A. Wensel, W. Whaling (Kellogg Radiation Laboratory, California Institute of Technology, Pasadena, USA). "A State in C¹² Predicted from Astrophysical Evidence." *Physical Review* 92, 649 (1953). F. Hoyle. "On Nuclear Reactions Occurring in Very Hot Stars. I. The Synthesis of Elements from Carbon to Nickel." *Astrophysical Journal Supplement* 1, 121 (1954).
[2] H. Kragh. "An anthropic myth: Fred Hoyle's carbon-12 resonance level." *Archive for History of Exact Sciences* 64, 721–751 (2010). 历史学家视角: 霍伊尔的预测源自恒星核合成问题, 非人择推理.
[3] A. Tohsaki (Kyoto University, Japan), H. Horiuchi (Kyoto University), P. Schuck (Université Paris-Sud, France / GSI), G. Röpke (Universität Rostock, Germany). "Alpha Cluster Condensation in ¹²C and ¹⁶O." *Physical Review Letters* 87, 192501 (2001). THSR波函数首次提出.
[4] H. Horiuchi (Kyoto University / RCNP, Osaka University, Japan), D. Blaschke (University of Wrocław, Poland / JINR Dubna). "THSR wave function and non-localized clustering." *European Physical Journal A* 60, 174 (2024). THSR形式体系综述: 扩展至非气体态团簇; ²⁰Ne中>99%重叠于Brink-GCM波函数.
[5] T. Otsuka (University of Tokyo / RIKEN Nishina Center, Japan), A. Volya (Florida State University, USA), N. Itagaki (Osaka Metropolitan University, Japan). "Theoretical studies of α clustering in nuclei and beyond." *European Physical Journal A* 62, 108 (2026). arXiv:2602.24175. 第一性原理MCSM展示α团簇自发涌现; 双旋转模式(致密vs分离)首次提出.
[6] Z.H. Yang, Y.L. Ye (Peking University, Beijing, China), et al. (RIKEN Nishina Center; Beihang University; Fudan University; IBS Korea). "Observation of the Exotic 0₂⁺ Cluster State in ⁸He." *Physical Review Letters* 131, 242501 (2023). Editors' Suggestion & Featured in Physics. α+2n+2n类BEC结构; 双中子对凝聚首次实验证据.
[7] 北京大学物理学院核物理与核技术国家重点实验室. "成果:杨再宏、叶沿林课题组实验发现丰中子滴线核氦-8中的中子对凝聚态." 国家自然科学基金委员会 (NSFC) 报道, 2023年12月.
[8] Y.L. Ye (Peking University). "BEC-analog states in light neutron-rich nuclei." INPC 2025 (The 29th International Nuclear Physics Conference), Daejeon, Korea, May 2025. ¹⁴C中的α+α+⁶He候选态初步结果.
[9] J.H. Chen, Y.L. Ye, K. Ma, J.X. Han, D.X. Wang, C.J. Lin, H.M. Jia, L. Yang, G. Li, L.S. Yang, Z.Y. Hu, Z.W. Tan, K. Wei, W.L. Pu, Y. Chen, J.L. Lou, X.F. Yang, Q.T. Li, Z.H. Yang, T.P. Luo, D.H. Huang, S.H. Zhong, Z.H. Li, J.Y. Xu (Peking University; China Institute of Atomic Energy; Guangxi Normal University). "New evidence of the Hoyle-like structure in ¹⁶O." *Science Bulletin* 68, 1119–1126 (2023). 4α窄共振在α+¹²C(Hoyle)衰变道中首次明确识别.
[10] K. Wei, J.H. Chen, Y.L. Ye (Peking University). "Progress and perspective on the study of the BEC-analog structure in nuclei." *Nuclear Physics Review* 42, 2023091 (2025). ChinaXiv:202402.00174 (2024). THSR理论+实验BEC态系统综述.
[11] K. Wei, J.H. Chen, Y.L. Ye (Peking University). "Clustering in nuclei: progress and perspectives." *Nuclear Science and Techniques* 35, 216 (2024). 团簇物理多前沿综述.
[12] X. Sun, et al. "The α condensate states of atomic nuclei ¹²C, ¹⁶O and ²⁰Ne in an analytical solvable model." arXiv:2501.09406 (2025). 解析可解模型系统预测α凝聚能: ²⁰Ne 5α凝聚~1 MeV.
[13] J. Kim, B. Fore, A. Lovato (Argonne National Laboratory, USA), et al. "Neural quantum states for cold Fermi gases in the BCS-BEC crossover." *Communications Physics* 7, 148 (2024). NQS首次用于核多体问题的BEC渡越区.
[14] A. Lovato (Argonne National Laboratory) et al. "Neural Quantum States for dilute nucleonic matter and self-emerging clusters." APS DNP Meeting / ECT* Workshop, 2025. NQS展示稀薄核物质中团簇自发涌现.
[15] N. Itagaki (Osaka Metropolitan University), H. Horiuchi (RCNP, Osaka University). "Container evolution of cluster structure with excitation energy." In: Horiuchi & Blaschke, *EPJ A* 60, 174 (2024). THSR框架向更重核和容器演化图像的推广.
[16] P. Schuck (Université Paris-Sud / GSI), Y. Funaki (RIKEN Nishina Center), H. Horiuchi (RCNP), G. Röpke (Universität Rostock), A. Tohsaki (Kyoto University), T. Yamada (Kanto Gakuin University). "Alpha-particle condensation in nuclei: A review." *Progress in Particle and Nuclear Physics* 89, 101 (2016). BEC类比奠基性综述.
[17] S. Shen, T.A. Lähde, D. Lee, U.-G. Meißner (Forschungszentrum Jülich / Universität Bonn). "Emergent geometry and duality in the carbon-12 nucleus." *Nature Communications* 14, 2777 (2023). 核格点EFT: 霍伊尔态为钝角三角形3α构型; 等边三角形构型存在壳模型-团簇双描述.
[18] C. Drischler, J.W. Holt, C. Wellenhofer (Ohio University / Texas A&M / Technical University of Munich). "Chiral effective field theory and the high-density nuclear equation of state." *Annual Review of Nuclear and Particle Science* 71, 403 (2021). 手征EFT对中子星内壳EoS的约束; α团簇对低密度不稳定性的调制.