量子多体理论与计算核物理方法
2024年，Nature发了一篇论文，标题很朴素——"波函数匹配"。
Elhatisari、Lee和Meißner做了一件事：把手征有效场论搬到时空格点上，用蒙特卡洛方法直接模拟原子核。从氘核到铁-56——两个核子到五十六个核子——束缚能和电荷半径的误差只有几个百分点。在此之前，格点方法只能处理轻核，重一点的核就会出现符号问题，计算直接崩溃。他们绕开的办法很巧妙：先在一个Wigner SU(4)对称性保护的无符号问题模型上做计算，再把波函数"匹配"回现实的手征核力。这个trick让格点EFT从A≤16一下子跳到了A=56。
差不多同一时间，橡树岭国家实验室的Hagen组在开发一套叫N3LO_Texas的手征核力。跟以往的做法不同，他们把氧-16的束缚能和电荷半径放进了核力拟合程序里。这听起来像是细节——但过去没有人这样做，因为氧-16太大了，从头算一次要消耗太多计算资源。Hagen组用耦合簇方法硬算，把氧-16做进了核力的低能常数校准中。结果是第一套可以从氘核到铅-208同时正确描述束缚能和电荷半径的手征核力。
这两件事加在一起，标志着一个转折点。从头算核结构——从基本核力出发，不做任何唯象近似——终于走到重核区了。
走到，但不等于走稳了。
---
先退一步看这个领域的整体结构。
核物理有一个独特的问题：核力不是基本力。基本力是量子色动力学——夸克和胶子之间的强相互作用。但在核的尺度上（10⁻¹⁵米），夸克被禁闭在核子里，核子之间的剩余强相互作用才是主角。这个剩余相互作用太复杂了，没法从QCD直接推导出来——除非你愿意在格点上模拟整个核，而这个计算量在可预见的未来都是天文数字。
所以核物理的"基本理论"其实是一层一层的。最底层是手征有效场论——它从QCD的对称性出发，用π介子交换构造核子之间的力，包含两体力、三体力、四体力……逐阶展开。手征EFT给了我们一个系统可改进的核力框架。往上是量子多体方法——怎么用这个核力去解A个核子的薛定谔方程。再往上是唯象模型——密度泛函、壳模型——它们不一定从手征核力出发，但可以覆盖整个核素图。
2024到2026年，这个三层结构正在打通。三种过去不太对话的方法——从头算、密度泛函、壳模型——开始激烈地交互。
先看从头算这边。
方法上有七条路。NCSM（无核芯壳模型）在谐振子基里精确对角化全部A个核子的哈密顿量——最干净，但只能算到A=16左右。耦合簇方法（CC）用指数参数化波函数，自动满足广延性，已经算到了铅-208。IMSRG（介质中相似重整化群）通过连续幺正变换把关联吸收进有效算符——变体最丰富，VS-IMSRG退耦到价空间接壳模型对角化，D-IMSRG处理形变，G-IMSRG处理连续谱。量子蒙特卡罗（GFMC/AFDMC）在坐标空间中随机行走求精确解——精度最高但受限于A≤12-20。自洽格林函数（SCGF）直接解戴森方程，自然输出谱函数和碎裂信息。格点EFT把时空离散化后在格点上模拟。多体微扰论（MBPT）从HF参考态出发逐阶加微扰——HF基下三阶就能和CC一致。
这些方法在重叠区域互相验证。CCSD(T)和IMSRG(2)在钙-48的结合能上差异不到1%——大概2-5 MeV。MBPT三阶就能给出跟CC一致的结果。但也有出现裂缝的地方。比如NCSM的谱因子——2024年Xu等人的一项研究发现了让人不安的事：即使能谱收敛得非常漂亮，谱因子也没有确定的收敛趋势。这对于过去二十年用NCSM谱因子计算敲出反应截面的做法是个直接挑战。
不过最戏剧性的张力不出在方法之间，出在核力本身。
2025年7月，Cirigliano、Dawid、Dekens和Reddy在PRL上发了一篇论文，声称某类三核子力被"增强"了——超出朴素量纲分析的预期。如果他们的论证成立，这类三体力应该从N⁵LO提升到N³LO，对核物质产生非常大的效应。这捅了马蜂窝。
2026年4月，Epelbaum、Gasparyan、Gegelia、Hog和Krebs在PRC上发表了反驳：增强标度对应的是KSW重整化方案，不是核结构计算实际采用的Weinberg方案。在Weinberg方案中，这些项就是O(1)，跟量纲分析一致。核心反驳词是"方案依赖性"——重整化群论证的结论取决于你选什么方案。
2025年3月，INT在华盛顿大学举办了一周的研讨会"手征EFT：新视角"。会议承认了Weinberg方法的根本问题：它产生的核振幅是不可重整化的，LO在α粒子以上就不稳定了，需要靠短程相互作用的精细调节来描述轻核，而且很难同时描述轻核和中重核。
这不是参数调一调的问题。这是框架层面的裂缝。二十年来手征EFT一直是核物理的"标准模型"——现在有人问这个标准模型本身的内部一致性有没有被严格证明过。
---
密度泛函和壳模型这边，进展的基调不太一样。不是"裂缝"，而是"打通"。
DRHBc——形变相对论连续谱Hartree-Bogoliubov理论——完成了偶Z核全局质量表。4829个束缚核素，覆盖Z=8到120。用PC-PK1点耦合泛函，RMS偏差约1.433 MeV——在所有DFT类方法中排名第一。更重要的是，2025年Qu等人用296个2021-2024年间新测量的质量数据做了基准检验，PC-PK1的RMS进一步降到~1.35 MeV。奇Z核的表还在计算中，三轴版本（TRHBc）也已经在发展。
DRHBc不是唯一的突破。广义TDGCM——含时生成坐标方法——被扩展到了包含动力学对关联。李志攀、赵鹏巍、孟杰团队在PRC 111上发表了²⁴⁰Pu裂变的电荷产额和总动能释放的首次无自由参数同时描述。这对裂变理论来说是件大事：之前没有任何方法能从同一个微观框架同时给出这两个量。
壳模型这边也在经历质变。Cao和Jiao在PLB 871上第一次把VS-IMSRG加上投影生成坐标方法（PGCM）扩展到完整的sdfp两个主壳层——这是一个极其庞大的价空间。结果很清楚：扩大价空间后，N=20反转岛的2⁺激发能被显著压低，B(E2)跃迁增强——跨壳的多粒子-多空穴激发对打破N=20幻数是不可或缺的。这在物理上是预期的，但之前从未被实际计算证实过。
而东京大学CNS的大塚孝治等人继续推进蒙特卡罗壳模型——MCSM在镍同位素中发现了球形、扁椭、长椭和gamma不稳定四种形状的共存。特苏诺达和清水等人发展的准粒子真空壳模型（QVSM）已经能处理钕-钐区域的球形到形变形状相变。更重要的是VS-DMRG——把IMSRG退耦后的有效哈密顿量用密度矩阵重整化群对角化。在镍-78上，DMRG给出了稳健收敛的激发能3.007 MeV，而传统组态相互作用方法在约2200亿个Slater行列式的空间中挣扎，给出的值是3.141 MeV，而且误差棒大得多。
在中国这边。2026年2月，中山大学尧江明团队的PRL 136——幻数起源的微观机制。这是一篇Editors' Suggestion加Featured in Physics。第一作者是直博生丁晨蓉。她们的发现很优雅：当你降低核力的动量分辨率——相当于从高分辨率显微镜切换到粗粒化视角——核子的壳结构经历一次相变，从自旋对称性变成赝自旋对称性。幻数28、50、82在这个过程中自然涌现。三体力在这个过程中起了关键作用——增强了自旋-轨道劈裂。
北京大学许甫荣、裴俊琛团队在PLB 862上首次把IMSRG用于致密核物质计算——这是国内中子星物态方程研究首次有了自主的从头算平台。甄昕是第一作者。同一年他们还实现了基于NCSM加手征核力的氚核量子计算——在国产16量子比特超导芯片上，GC分组测量加对称性投影让效率提升了三个量级。
清华龙桂鲁团队的QCSH——全量子壳模型——用LCU形式在纯量子计算机上计算了从氘到钙-40共12个核素的结合能。这是2023年的工作，到2025年持续被国际引用。跟VQE这种量子-经典混合算法不同，QCSH是纯量子算法，避免了经典和量子之间慢速切换的瓶颈。
福州大学吴鑫辉、北大任中洲和赵鹏巍的合作——机器学习无轨道密度泛函，2025年发表在Communications Physics上——是个真正的方法论突破。他们用核岭回归在24000个数据点上训练了一个直接映射：核子密度进去，动能加自旋轨道能出来。不需要显式的单粒子轨道。在球形氧-16和形变氖-20上都成功重现了结合能、半径、四极形变和势能曲线。计算复杂度从Kohn-Sham的O(N³)降到O(N)。这是一个教科书级别的挑战被攻克了——Hohenberg-Kohn定理说了可以做，但从1964年到现在，没人能在形变核上真的做出来。
---
量子计算这一块，氛围不太一样。它还不是"工具"，而是"赌注"。
VQE（变分量子本征求解器）已经从轻核推进到了中等质量核。Bhoy和Stevenson 2024年用12个量子比特模拟了镍-58的fp壳——哈密顿量有1691个Pauli项。Sarma和Stevenson 2025年在IBM真实硬件上模拟了锂到铅-210，经零噪声外推降噪后偏差不到4%。Singh等人引入了Gray-code编码，把量子比特需求大幅降低。Costa等人用配对基编码方案把CNOT门数降了两到三个数量级。自适应ansatz（ADAPT-VQE）、对称性编码、新量子比特映射策略——这是三大创新方向。
但经典方法也在进步。VS-DMRG在单个DGX-H100节点上达到了约250 TFLOPS。2024年，中国团队在神威超算上用4190万核验证了Google Sycamore的量子霸权实验——经典计算并未停滞。量子优势的起跑线是一个移动靶。
美国能源部2025年的评估说，供应商路线图和科学应用资源需求预计在5到10年内交汇重叠——也就是2030到2035年。到那时，核物理可能需要1000到2000个逻辑量子比特。NISQ时代能做的目前还很有限——目前达到p壳和sd壳的轻核，经典方法在这些系统上根本不费力。
张量网络方法可能是一个更现实的路径。Rausch等人在Annalen der Physik上证明了原子核在配对哈密顿量视角下是低纠缠度系统——可以用一维MPS/DMRG高效描述，计算成本和BCS差不多但没有平均场的伪影。他们计算了铅-178到铅-220所有同位素的奇偶质量差，铅-208的前一百个激发态，铅-210的双中子移除谱函数——都是机器精度。
神经网络量子态（NQS）加变分蒙特卡洛是另一条路。Gnech等2024年在PRL上计算了A≤20核素的结合能和电荷半径，达到百分之一精度——而且通过一个外磁场方案，NQS自发涌现了核壳层结构。Fore等人把NQS用到了中子星壳层相关低密度核物质——在低密度区域显著优于辅助场扩散蒙特卡洛，还观测到了Cooper配对的自发涌现。
冷原子量子模拟这边，Halimeh等2025年在Nature Physics上综述了格点规范理论的量子模拟进展。冷原子光晶格已经成为模拟Z₂、U(1)、SU(N)格点规范理论的领先平台。NCSU和MIT的Schäfer、König、Zwierlein团队有一个2023-2026的DOE项目，用超冷原子直接模拟稀薄中子气体的动力学和输运性质。
---
如果我们退一步看全景，有三条线索正在交织。
第一条：三种方法层级正在互相渗透。VS-IMSRG产出壳模型的有效相互作用——这是从头算输送到壳模型的管道。手征EFT的π交换项被嵌入Skyrme泛函的接触项中——这是核力输送到DFT的管道。SCGF和AFDMC的从头算结果被用作局域密度近似来构建EDF——这是微观约束唯象理论的管道。它们不是三座孤岛，而是同一个三角洲上的三条河道。
第二条：计算和理解的张力正在升级。ML无轨道DFT证明了Hohenberg-Kohn定理在形变核中的实际可行性——壳效应不需要显式的单粒子轨道。这套直接挑战了"壳效应需要组态混合"的传统观念。NuCLR——Nolte等人的多任务深度学习模型——在没有明确物理指令的情况下自发涌现了壳结构的内部表征。Trifinopoulos等人解剖了NuCLR的潜空间，发现它学会了Jaffe在1969年发现的能量分解函数。ML不是在做"记忆"，它在做"压缩"——从数据中提炼物理规律。
第三条：中国正在从"参与"转向"引领"。尧江明团队的PRL 136幻数起源是领域顶峰级工作。DRHBc全局质量表是国际唯一的自洽引入形变加对关联加连续谱的协变DFT质量表。ML无轨道DFT是方法论的原创贡献。QCSH是国际上最早的系统性全量子壳模型方案之一。C3NT——华中核理论中心——2025年5月正式成立，标志着系统性加强核理论研究的战略布局。
差距也在：耦合簇理论的重核应用中国几乎没有团队在做。裂变动力学全微观模拟基本空白。致密核物质的微扰QCD约束完全依赖国际团队。核天体物理r过程全局网络的系统性研究不足。手征核力本身的开发——中国的参与度很低。
---
回到那个根本问题：我们能用第一性原理算清楚一个原子核吗？
2024到2026年的答案是：能——但只是某些核，用某些方法，在某个精度内。
耦合簇算到了铅-208。格点EFT算到了铁-56。IMSRG加GCM可以处理形变核的形状共存。密度泛函可以覆盖整个核素图——但同位旋矢量道可能有个系统性缺陷还没解决。壳模型可以做精确谱学——但重核的价空间大到连蒙特卡罗都吃力。量子计算做出了原理演示——但离"比经典快"还有至少五到十年。ML代理模型有惊人的精度——但外推到训练域以外的可靠性没人敢打包票。
最诚实的评估可能是：多条路径正在交汇。从头算往下渗透到壳模型和DFT，DFT和壳模型往上吸收从头算的约束，量子计算和张量网络从侧面开辟新维度的计算资源。全面的"量子多体统一理论"还不存在——但各条路径的交叉点正在变得越来越密集。
而那个最终的问题——核力到底是什么？——Cirigliano和Epelbaum的争论给出了一个诚实的答案：我们用了二十年的标准框架，内部一致性还没有被严格证明。
---
参考文献
[1] Elhatisari S., Lee D., Meißner U.-G. et al., "Wave function matching for the quantum many-body problem," *Nature* 630, 59-64 (2024). (KFUPM / FRIB-MSU / Universität Bonn)
[2] Sun Z.H., Ekström A., Forssén C., Hagen G., Jansen G.R., Papenbrock T., "Multiscale physics from first principles: Coupled cluster with angular momentum projection," *Physical Review X* 15, 011028 (2025). (ORNL / UT Knoxville / Chalmers University)
[3] Hagen G. et al., "N3LO_Texas chiral interaction for global ab initio nuclear structure," to be published. Presented at PAINT 2025, INT Workshop. (ORNL)
[4] Heinz M., Miyagi T., Stroberg S.R., Tichai A., Hebeler K., Schwenk A., "IMSRG(3)-N^7: Systematic quantification of many-body truncation uncertainties," *Physical Review C* 111, 034311 (2025). (TU Darmstadt / ORNL / University of Tsukuba / Notre Dame)
[5] Belley A., Stroberg S.R., Hergert H. et al., "Ab initio uncertainty quantification of the neutrinoless double-beta decay nuclear matrix element in ^76Ge," *Physical Review Letters* 132, 182502 (2024). (FRIB / MSU / TRIUMF / Notre Dame)
[6] Yaghi A., Hupin G., Navrátil P., "Complex scaling within the similarity renormalization group framework: Application to ab initio no-core shell model," arXiv:2507.01595 (2025). (TRIUMF / IJCLab / LLNL)
[7] Cirigliano V., Dawid M., Dekens W., Reddy S., "Enhanced three-nucleon forces from renormalization group," *Physical Review Letters* 135, 022501 (2025). (INT / UW / LANL)
[8] Epelbaum E., Gasparyan J., Gegelia J., Hog H., Krebs H., "On the renormalization of pion-exchange contact three-nucleon forces," *Physical Review C* 113, 044005 (2026); arXiv:2512.14117. (Ruhr-Universität Bochum / JGU Mainz / Tbilisi State University)
[9] 丁晨蓉, Wang C.C., Hergert H., Bogner S.K., Liang H.Z., 尧江明, "Origin of nuclear magic numbers: A phase transition from spin to pseudospin symmetry from chiral nuclear forces," *Physical Review Letters* 136, 052501 (2026). (Sun Yat-sen University / CAEP / MSU-FRIB / University of Tokyo)
[10] Zhou E.F., 丁晨蓉, Luo Q.Y., 尧江明, Hergert H., "Shape coexistence in neutron-rich Mg isotopes from ab initio IM-GCM," *Physics Letters B* 865, 139464 (2025). (Sun Yat-sen University / FRIB-MSU)
[11] 甄昕, 裴俊琛, 许甫荣, "Nonperturbative ab initio calculation of dense nuclear matter with IMSRG and chiral nuclear forces," *Physics Letters B* 862, 139350 (2025). (Peking University)
[12] Guo P. et al. (DRHBc Mass Table Collaboration), "Deformed relativistic Hartree-Bogoliubov theory in continuum: The DRHBc mass table for even-Z nuclei," *Atomic Data and Nuclear Data Tables* 158, 101661 (2024). (~30 institutions worldwide, initiated by Peking University / IMP)
[13] Qu X.B., Chen Y., Pan C., Yu K.Y., Zhang K.Y., "Benchmarking global nuclear mass models with 296 new experimental masses," *Nuclear Science and Techniques* 36 (2025); arXiv:2505.09914. (IMP / PKU / et al.)
[14] 赵鹏巍 et al., "Rotational anomalies in transfermium nuclei: Octupole deformation in ^252,254No," *Physical Review Letters* 133, 022501 (2024). (Peking University)
[15] Li Z.P., Vretenar D., Nikšić T., 赵鹏巍, 孟杰, "Generalized time-dependent generator coordinate method with dynamical pairing for fission dynamics," *Physical Review C* 111, L051302 (2025). (Southwest University / Zagreb / Peking University)
[16] Cao L.J., Jiao J.J., "Ab initio VS-IMSRG+PGCM in two complete major harmonic oscillator shells for N=20 island of inversion," *Physics Letters B* 871, 140034 (2025). (Zhengzhou University)
[17] Wu X.H., Ren Z.X., 赵鹏巍, "Machine learning orbital-free density functional theory for deformed nuclei," *Communications Physics* 8, 316 (2025). (Fuzhou University / Peking University / FZ Jülich / Universität Bonn)
[18] Bhoy B., Stevenson P.D., "Variational quantum eigensolver for ^58Ni shell-model calculation," *New Journal of Physics* 26, 075001 (2024). (University of Surrey)
[19] Rausch R., Plorin S., Peschke M., Karrasch C., "DMRG solves nuclear pairing with machine precision at BCS cost," *Annalen der Physik* 536, 2300436 (2024). (TU Braunschweig)
[20] Gnech A., Lovato A., Carleo G. et al., "Neural-network quantum states for ab initio nuclear structure," *Physical Review Letters* 133, 142501 (2024). (ANL / EPFL / TRIUMF)
[21] Fore B., Kim J.M., Hjorth-Jensen M., Lovato A., "Neural-network quantum states for low-density neutron matter," *Communications Physics* 8, 1 (2025). (ANL / MSU / University of Oslo)
[22] Halimeh J.C. et al., "Quantum simulation of lattice gauge theories with cold atoms," *Nature Physics* 21, 25-36 (2025). (Heidelberg University / JHU / University of Maryland)
[23] Lv P., Wei S.J., Xie H.N., 龙桂鲁, "QCSH: A full quantum computer nuclear shell model package," *Science China Physics, Mechanics & Astronomy* 66, 240311 (2023); arXiv:2205.12087. (Tsinghua University / BAQIS)
[24] Sarma C., Stevenson P.D., "Slater determinant to qubit mappings for VQE on real quantum hardware," arXiv:2510.02124 (2025). (University of Surrey)
[25] Trifinopoulos S., Kitouni O., Nolte N., Williams M. et al., "From neurons to neutrons: Mechanistic interpretability of neural network mass models," ICML 2024; arXiv:2508.08370 (2025). (MIT/IAIFI / Meta FAIR / Harvard)
[26] Lee D., "Nuclear lattice effective field theory: A systematic framework for nuclear structure and reactions," *Annual Review of Nuclear and Particle Science* 75, 109-128 (2025). (FRIB / MSU)
[27] Nolte N. et al., "NuCLR: Nuclear Co-Learned Representations — Self-supervised discovery of nuclear concepts from neural networks," ICML 2023 Workshop; ongoing work 2024-2025. (TU Darmstadt / MIT)
[28] Tichai A., Legeza Ö., Knecht S. et al., "VS-DMRG: Valence-space density matrix renormalization group for ab initio shell model," *Physics Letters B* 845, 138139 (2023); extended to ^78Ni and N=50 isotones (2024). (TU Darmstadt / Wigner RCP)
[29] 许甫荣, 孟子颜 et al., "手征核力与核结构第一性原理计算综述," *中国科学:物理学 力学 天文学* (2025); *科学通报* 71 (2025). (北京大学)
[30] INT Workshop, "Chiral EFT: New Perspectives," Institute for Nuclear Theory, University of Washington, March 10-14, 2025.
[31] Singh A., Siwach P., Arumugam P., "Gray-code encoding and qubit-ADAPT-VQE for nuclear shell model," *Physical Review C* 112, 034320 (2025). (IIT Roorkee)
[32] 潘琮 et al. (DRHBc Mass Table Collaboration), "Pseudo-DRHBc mass table for odd-Z nuclei and r-process applications," *European Physical Journal A* 61, 275 (2025); arXiv:2503.09324. (~30 institutions)