跳到正文
返回核结构
至真探索 · 核结构

核物理第一性原理计算NCSM-IMSRG-QMC-SCGF

2026-05-18

# 核物理第一性原理计算：直接解多体薛定谔方程

假设你手里有一个公式，写出所有质子和中子之间的相互作用力——质量、电荷、自旋、交换项全在里面。理论上，把 A 个核子放进这个公式，解一个 A 体薛定谔方程，就能得到原子核的所有信息：结合能、激发态、能谱、反应截面、形变……物理学家把这叫"从第一性原理出发"。

但 A 等于 1 的时候很容易。A 等于 2 的时候也很容易。从 A 等于 3 开始，麻烦来了。N 个量子粒子构成的体系，它的希尔伯特空间维度随粒子数指数增长。4 个核子的空间就已经大到必须在超级计算机上做取舍。一般原子核里有几十到上百个核子，不能直接解任何超过 A ≈ 5-6 的全空间薛定谔方程——组合爆炸会瞬间吞噬任何可用计算资源。

所以四套方法诞生了。它们面临的是一样的难题——怎么在指数爆炸的多体空间里找到物理上有意义的近似——但走了四条完全不同的路。NCSM、IMSRG、QMC、SCGF，四条路的终点是一致的：给任意原子核（不管它有多重）提供从第一性原理出发的基态和激发态预言。但它们走到终点的路径、擅长什么、卡在哪里，都不同。这四套方法放在一起，才构成核物理第一性原理计算的全景。

四条路要解决同一个核心困难：多体薛定谔方程在核子数 A 增长时，自由度急剧膨胀。这个困难本身是量子多体物理的普遍问题，但核物理有它的特殊性——核力既不是两体力在三体里可以包住的简单形式，也不是纯排斥的费米统计。核力有各种量子数依赖、张量耦合、以及短程硬核，让传统的量子化学方法（比如耦合簇理论对电子气体的处理）难以直接套用。

核物理第一性原理计算在 2010 年代后期经历了一次质的飞跃。2020 年代，它们已经能处理 A ≈ 100 的核，并且在某些同位素链上的结合能计算做到了 2% 以内的精度。这个精度在十年前看来是不可能的。下面逐一展开这四条路线。

---

NCSM（无核芯壳模型） 是最直接、最"暴力"的方法——把所有核子放在一个谐振动子基组里，把哈密顿量对角化，然后解出基态和激发态。它叫"无核芯"，因为它不做平均场近似——所有 A 个核子都被平等对待，没有"冻结"在核芯里的粒子。这个设定使它在对称性上是干净的，结果不受任何平均场假设的污染。

代价是计算量随 A 指数增长。16 个核子（氧-16）的 NCSM 模型空间可以达到 10¹⁰ 维以上。为了把问题压到能算的范围，NCSM 有两个关键部件。一个是截断：把振荡量子数总和超过某个值的基组去掉。另一个是重整化：利用 Lee-Suzuki 变换把原始核力做"软化"，让它在截断空间里也能给出合理的低能物理。软化之前的核力在短程有很强的排斥芯——两个核子靠得太近时排斥力剧烈上升。这个排斥芯在多体计算里很难处理，因为需要海量的高能基组才能描述它的量子态波函数的正确节点行为。Lee-Suzuki 变换在保持低能物理不变的前提下，把高能部分"积分掉"，得到等效的低能哈密顿量——这是 NCSM 成功的关键。

NCSM 的经典成果来自 Navrátil 等人（2007，J Phys G 34, S731）——他们用 NCSM 计算了锂-6、铍-10 等轻核的结合能和激发态能谱，与实验数据在约 1 MeV 以内吻合。2013 年 Barrett、Navrátil 和 Vary（Prog Part Nucl Phys 69, 131）系统回顾了 NCSM 的全部发展，把 NCSM 推到了氧-24 和氟-25。

但 NCSM 的上限——至少对标准 NCSM 而言——是 A ≈ 16 左右。超过这个质量数，模型空间的维数就开始失控。无核芯蒙特卡洛壳模型（NCSM-MCSM）可以拉到 A ≈ 20。Abe 等人（2021，Phys Rev C 104, 054315）用 NCSM-MCSM 计算了从氦-4 到氖-20 的 4n 自共轭核的基态性质——这是 NCSM 族方法系统性地推向中重核的重要一步。但必须指出，NCSM-MCSM 目前对集体转动特征的精确再现仍有距离，在氖-20 上的计算与实验的偏差比轻核更大。

---

IMSRG（核内相似度重整化群） 走了另一条路。它不直接对全空间对角化，而是通过一系列连续么正变换，把哈密顿量"流"到一个预先选定的简化基态附近。想象一个流动的哈密顿量，从一个初始的高维形式开始，随着某个参数的演化，逐渐"冻结"掉基态以外的自由度，最终得到局限在基态附近的有效哈密顿量。这个过程不截断空间——它在演化过程中始终保持哈密顿量的完整结构——但最终只需要在很小的子空间里对角化。

IMSRG 的开创性工作来自 Hergert、Bogner、Schwenk 等人。2016 年他们发表了完整的理论综述（Hergert et al., Phys Rep 621, 165），奠定了 IMSRG 的方法学基础。2017 年 Hergert 又发表了关于闭壳和开壳核的 IMSRG 方法的综合讨论（Hergert, Phys Scr 92, 023002）。

IMSRG 的优点是它随 A 增长的缩放不像 NCSM 那样恐怖——大约是 A² 到 A³，而不是指数级。这使得 IMSRG 可以一直算到 A ≈ 100。代价是在每一步演化中，哈密顿量的形式都是截断在某种算子水平（通常是两体算符）上的——所谓 IMSRG(2)——因为三体算符的存储和计算成本跟 A⁶ 成正比。这个截断引入了系统误差。

Stroberg、Holt、Schwenk 和 Simonis 用 VS-IMSRG+χEFT 核力计算了从氦到铁（Z = 2-26）全同位素链的结合能和分离能。2021 年他们发表了里程碑式的工作：预测了从轻核到中等质量核的滴线（Stroberg et al., Phys Rev Lett 126, 022501），涵盖了接近 700 种核素。他们也在 Annual Review of Nuclear and Particle Science 上撰写了综述（Stroberg et al., Annu Rev Nucl Part Sci 71, 327），系统介绍了 valence-space IMSRG 方法。这些结果的重要意义在于：它没有用到任何在目标核素上调参数的自由度——IMSRG 的全部 LEC 都来自对 NN 散射和氘、氚、氦-3、氦-4 的数据拟合，然后直接外推到中重核区。2-3% 的偏差不是参数拟合的"结果"，而是理论真正预测能力的度量。

IMSRG 的三体算符问题是近年来的一个重要焦点。2024 年 Stroberg、Morris 和 He 完成了 IMSRG 中三体算符完全演化的计算（Stroberg et al., Phys Rev C 110, 044316）：对氧-22 和钙-40 等核，三体算符的完全包含让结合能相对于 IMSRG(2) 提高了大约 1.5-2%，跟实验的偏差降到了 1% 以内。这印证了理论家的一个共识：两步截断（同时截断在二阶流方程和三体算符）中的第二步其实是更重要的误差来源。更系统的 IM-SRG(3) 理论框架由 Tichai 等人（2024，Phys Rep 1054, 1）完成综述。

在开壳核（价核子不在闭壳层上的核）方面，IMSRG 也取得了显著进展。闭壳核（氧-16、钙-40、钙-48）的 IMSRG 计算已经很成熟；开壳核需要先做一个参考态变形破缺、然后投影恢复对称性。多参考态 IMSRG（MR-IMSRG）和多粒子数投影 IMSRG 已经在钙同位素链和氩同位素链上展示了系统性能力（Hergert et al., Phys Rep 621, 165; Hergert, Phys Scr 92, 023002），使 IMSRG 不再只是球形核的工具。2024 年 Xu 等人（Nucl Sci Tech 35, 215）进一步综述了 IMSRG 与耦合道方法及连续谱耦合的进展。

---

QMC（量子蒙特卡洛） 是计算量最大但理论上最精确的方法。它直接对多体波函数做随机抽样，不需要任何有限基组的截断，也不需要有效的哈密顿量——直接使用裸的核力。QMC 在少数约束下可以对波函数做最准确的描述，优势在于它天然"感受"了连续空间，没有任何离散化误差。

但 QMC 有几个致命局限。一个是符号问题——费米子的交换对称性导致波函数在某些区域为负，破坏了概率解释，引起统计误差发散。核子数越多，符号问题越严重。标准 GFMC（格林函数蒙特卡洛）只能处理 A ≤ 12 左右的核（碳-12 已经是极限）。AFDMC（辅助场扩散蒙特卡洛）通过 Hubbard-Stratonovich 变换把自旋-同位旋自由度用辅助场积分掉，减轻了符号问题，可以把边界推到 A ≈ 40 甚至更重。第二个限制是核力的形式——三体力必须被简化为"准二体"形式才能在 QMC 框架内高效计算，但这个简化在重核区的精度未经系统验证。

2015 年 Carlson、Gandolfi、Pederiva、Pieper、Schiavilla、Schmidt 和 Wiringa 发表了全面的 QMC 核物理综述（Rev Mod Phys 87, 1067），详细介绍了 VMC、GFMC 和 AFDMC 三种方法。2020 年 Gandolfi 和 Lonardoni 更新了 QMC 结合 χEFT 核力的新进展（Gandolfi & Lonardoni, Front Phys 8, 117）。

GFMC 方面近年来的进展包括 Lonardoni 等人（2018，Phys Rev C 97, 044313）用 GFMC + 局域手征核力计算了碳-12 和氧-16 的性质，展示了手征 EFT 在 GFMC 框架中的可行性。

AFDMC 方面最重要的进展来自 Hu、Sun、Hagen、Jansen 和 Papenbrock（2024，Phys Lett B 858, 139010）。他们用耦合簇方法（而非 AFDMC——这里需要澄清）计算了从镍-78 沿 N = 50 向钙-70 的核素链。这篇工作首次用 ab initio 方法成功再现了镍-78 的低位转动带，并预测了铁、铬、钛核中的基态形变。AFDMC 本身在重核上的直接计算仍然受限于符号问题和三体力的简化处理——声称 AFDMC 已经算到镍-78 的说法需要审慎对待。目前 AFDMC 在 A ≈ 40 以下的核上结果可靠。

关于费米子符号问题的处理，近年来也有新思路。2022 年 Li 等人提出了微扰量子蒙特卡洛方法（Phys Rev Lett 128, 242501），通过微扰修正缓解约束 QMC 中的符号问题偏差。2025 年进一步的 sign-problem-free QMC 方法也在发展中，但目前仍限于相对较轻的核。

---

SCGF（自洽格林函数法） 跟前面三个方法都不一样。它不直接算波函数，而是算"传播子"（格林函数）——描述一个核子被加入或移出系统时，系统如何响应。从格林函数里可以直接提取谱函数、分离能、能谱加权求和规则。SCGF 的路线很聪明：与其直接面对 A 体波函数的复杂性，不如转去计算单粒子 Green 函数和双粒子 Green 函数，然后用一套自洽方程（Dyson 方程）把它们连起来。

这个方法的优势是可以直接处理"连续谱"——原子核并非孤立存在的，它的激发态天然是连续谱的一部分（在阈值以上），而 SCGF 在处理这方面比 NCSM 和 IMSRG 自然得多。它的缺点是计算三体力在 SCGF 框架内的实现极其复杂——三体力对应的是三粒子-三空穴顶点的自洽处理，这个计算的技术门槛比 IMSRG 还高。

SCGF 的理论基础由 Dickhoff 和 Barbieri 在 2004 年的综述中奠定（Prog Part Nucl Phys 52, 377）。2011-2014 年，Somà、Duguet 和 Barbieri 发展了 Gorkov-SCGF 方法，首次将 SCGF 推广到真正开壳的原子核（Somà et al., Phys Rev C 84, 064317; Phys Rev C 89, 024323）。

SCGF 在中等质量核上的代表作是 Somà 等人（2014，Phys Rev C 89, 024323）对钙和钛同位素链的谱函数计算。他们用 χEFT 核力，计算了这些中等质量核的价核子谱函数——能谱被"展宽"的效果（对应的效应是核子的有限寿命），跟高分辨率 (e,e'p) 实验数据有可比性。

2021 年 Frosini 等人（Eur Phys J A 57, 113）用 Gorkov-SCGF 方法系统计算了从氟到铁的七条完整同位素链的结合能和电荷半径，展示了 SCGF 从单幻数核向真正开壳核过渡的能力。SCGF 的连续综述由 Barbieri 等人完成，收录在 Hjorth-Jensen 等人编辑的教材中（Barbieri & Carbone, "Self-Consistent Green's Function Approaches", in Lecture Notes in Physics 936, Springer 2017）。

---

四条路线现在可以比较了。

NCSM 在 A ≤ 16 范围内提供了精度最高的"准精确解"（对 A ≈ 4-6 核，NCSM 事实上就是精确解），也提供了最干净的激发态谱描述。超过 A ≈ 20 之后，计算成本失控。NCSM-MCSM 在 A ≈ 20 附近提供了有限但有用的延伸。

IMSRG 是"A ≤ 100 全覆盖"的实际冠军——结合能精度 2-3%，在某些核区可达到 1% 以内。它对闭壳核已经非常成熟，对开壳核的 MR-IMSRG 扩展正在迅速进步。缺点是它对连续谱的处理弱（但已有 Gamow-IMSRG 等耦合方案出现）。

QMC 提供了一个原则上精确但实践上受限于符号问题和核力形式的替代方案。GFMC 在轻核区（A ≤ 12）提供了对 γ 跃迁、电弱流等算符的无截断计算——这是它的独有优势。AFDMC 正在向中等质量核推进，但在 A ≈ 40 以上，符号问题和三体力的简化处理仍是瓶颈。

SCGF 在连续谱、开壳、对关联这三个维度上有综合优势。它在中等质量核的谱函数、分离能、以及价核子寿命方面是最全面的信息来源。但 SCGF 的核力实现（特别是三体力）技术门槛最高。

四条路线之间如何交叉验证是一个重要的课题。目前已经有不少方法间的交叉基准测试，比如 NCSM 与 CC 在轻核上的一致性检验（Hagen et al., Phys Rev C 76, 044305），IMSRG 与 CC 在中等质量核上的对比（Hergert et al., Phys Rep 621, 165），以及 QMC 与 NCSM 在少数核子体系上的准精确解比较。但据我所知，截至 2026 年 5 月，尚未有 NCSM、IMSRG、QMC、SCGF 四种方法在同一核力输入下对同一种核的系统性"四路交叉验证"发表在 Nature Communications 等顶级期刊上。此前的文章中声称的"Keblasi 等人 Nature Commun 2025"这一引用，我没有找到可靠对应文献，应当删除。

即使没有四路交叉验证，现有的方法间比较已经表明——至少对闭壳核和中等质量核——不同的多体方法在相同的核力输入下可以给出相当一致的结果，偏差在 1-2% 以内。这暗示多体方法的系统误差在现有精度水平上被合理控制了，剩下的主要偏差来自核力本身。

核力本身的不确定性——特别是在同位旋破坏、三体力参数、以及手征展开截断误差等高阶项上——仍然限制着第一性原理计算的最终精度。目前最先进的 N³LO 核力在 A ≈ 50-100 区域的多体计算中，截断误差大约 2-5%，与多体方法本身的收敛误差相当。这意味着如果不进一步缩小核力的不确定性（推进到 N⁴LO 三体力+N⁵LO 四体力等），单纯改进多体方法的算法将不会显著提升总体精度。

---

交叉验证的结果指向一个问题：多体方法可以怎么继续改进？

NCSM 方向的改进重点是"重要采样 + 高效并行"。NCSM-MCSM 的下一步是使用张量网络（TN）或者量子计算机上的量子模拟来扩展可达的核子数。2022 年 Kiss 等人（Phys Rev C 106, 034325）首次在量子模拟器上实现了锂-6 核的变分计算，使用了 ordered unitary coupled clusters 方法。这是量子计算在核物理中的概念验证，结果与经典 NCSM 计算定性一致。这只是个起步，但暗示了未来的方向。

IMSRG 方向的改进重点是"三体算符的完整演化"和"连续谱耦合"。Stroberg 等人（2024，Phys Rev C 110, 044316）已经完成了三体算符在 IMSRG 中的完整演化计算。完整 IMSRG(3) 的计算复杂度大约是 A⁶，对 A ≈ 40 以下的核是可行的，对更重核则需要新的算法。

QMC 方向的改进重点是"更好的核力表示"和"符号问题的进一步抑制"。AFDMC 在 A ≈ 40 核上的应用仍有挑战——符号问题在 A ≈ 100 附近会重新变得失控。微扰 QMC 等新方法（Li et al., Phys Rev Lett 128, 242501, 2022）提供了缓解符号问题的思路，但距离实用仍有距离。

SCGF 方向的改进重点是"三体力的协调"和"开壳核的系统性扩展"。三体力 SCGF 的完整形式化已经建立（Somà et al., Phys Rev C 84, 064317），但数值上完全自洽地包含三体力仍在发展中。一旦完成，SCGF 应该是四个方法中唯一一个可以统一处理开壳、闭壳、连续谱的框架。

最后要说的是：这四个方法之间的竞争远弱于它们的互补性。NCSM 在高能谱和 γ 跃迁上提供最干净的基准点；IMSRG 在中重核的结合能和基态性质上提供最高效率；QMC 在电弱流算符上提供最准确的计算（因为不截断基组）；SCGF 在连续谱和开壳核上提供最宽的覆盖。把它们合在一起，第一性原理核物理从前只能在 A ≤ 4 的极轻核区立足，现在它的边界已经推到了 A ≈ 100 的区域。

在这个过程中，核物理第一性原理计算不仅仅是在模仿实验——它已经进入了能预测未测量核素的性质的阶段。在 FRIB 和 RIKEN 等放射性核束装置即将产生数千种新核素的时代，第一性原理计算提供了一种不依赖实验标定的理论预测能力。这个能力本身会驱动新的实验设计，也会让核物理的知识边界向外推进更大的一步。

---

参考文献

  1. Barrett B R, Navrátil P, Vary J P. Ab initio no-core shell model. *Prog Part Nucl Phys*, 2013, 69: 131-181.
  2. Hergert H, Bogner S K, Morris T D, Schwenk A, Tsukiyama K. The In-Medium Similarity Renormalization Group: A novel ab initio method for nuclei. *Phys Rep*, 2016, 621: 165-222.
  3. Stroberg S R, Holt J D, Schwenk A, Simonis J. Ab initio limits of atomic nuclei. *Phys Rev Lett*, 2021, 126: 022501.
  4. Stroberg S R, et al. Nonempirical interactions for the nuclear shell model. *Annu Rev Nucl Part Sci*, 2021, 71: 327-359.
  5. Carlson J, Gandolfi S, Pederiva F, Pieper S C, Schiavilla R, Schmidt K E, Wiringa R B. Quantum Monte Carlo methods for nuclear physics. *Rev Mod Phys*, 2015, 87: 1067-1118.
  6. Gandolfi S, Lonardoni D. Atomic nuclei from quantum Monte Carlo calculations with chiral interactions. *Front Phys*, 2020, 8: 117.
  7. Dickhoff W H, Barbieri C. Self-consistent Green's function method for nuclei and nuclear matter. *Prog Part Nucl Phys*, 2004, 52: 377-496.
  8. Somà V, Duguet T, Barbieri C. Ab initio Gorkov-Green's function calculations of open-shell nuclei. *Phys Rev C*, 2011, 84: 064317.
  9. Somà V, Barbieri C, Duguet T. Ab initio self-consistent Gorkov-Green's function calculations of medium-mass nuclei. *Phys Rev C*, 2014, 89: 024323.
  10. Navrátil P, et al. Recent developments in no-core shell-model calculations. *J Phys G*, 2007, 34: S731.
  11. Abe T, Maris P, Otsuka T, Shimizu N, Utsuno Y, Vary J P. Ground-state properties of light 4n self-conjugate nuclei in ab initio no-core Monte Carlo shell model calculations with nonlocal NN interactions. *Phys Rev C*, 2021, 104: 054315.
  12. Stroberg S R, Morris T D, He B C. In-medium similarity renormalization group with flowing 3-body operators, and approximations thereof. *Phys Rev C*, 2024, 110: 044316.
  13. Tichai A, et al. IM-SRG(3) with three-body operators. *Phys Rep*, 2024, 1054: 1.
  14. Hu B S, Sun Z H, Hagen G, Jansen G R, Papenbrock T. Ab initio computations from ⁷⁸Ni towards ⁷⁰Ca along neutron number N=50. *Phys Lett B*, 2024, 858: 139010.
  15. Frosini M, et al. Moving away from singly-magic nuclei with Gorkov Green's function theory. *Eur Phys J A*, 2021, 57: 113.
  16. Hergert H. In-Medium Similarity Renormalization Group for closed and open-shell nuclei. *Phys Scr*, 2017, 92: 023002.
  17. Hjorth-Jensen M, Lombardo M P, van Kolck U (Eds.). An Advanced Course in Computational Nuclear Physics: Bridging the Scales from Quarks to Neutron Stars. *Lecture Notes in Physics*, Vol. 936. Springer, 2017.
  18. Barbieri C, Carbone A. Self-Consistent Green's Function Approaches. In: Hjorth-Jensen M et al. (Eds.), *An Advanced Course in Computational Nuclear Physics*, Springer, 2017, pp. 361-436.
  19. Xu X Y, Fan S Q, Yuan Q, Hu B S, Li J G, Wang S M, Xu F R. Progress in ab initio in-medium similarity renormalization group and coupled-channel method with coupling to the continuum. *Nucl Sci Tech*, 2024, 35: 215.
  20. Lonardoni D, et al. Properties of nuclei up to A=16 with local chiral interactions. *Phys Rev C*, 2018, 97: 044313.
  21. Li Z, et al. Perturbative quantum Monte Carlo method for nuclear physics. *Phys Rev Lett*, 2022, 128: 242501.
  22. Kiss O, Grossi M, Lougovski P, Sanchez F, Vallecorsa S, Papenbrock T. Quantum computing of the ⁶Li nucleus via ordered unitary coupled clusters. *Phys Rev C*, 2022, 106: 034325.