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

让原子核动起来的那套工具箱

2026-09-25

*引用格式说明：* 正文中的 [N] 与文末参考文献列表第 N 条一一对应；文末另附「核实与纠错记录」，逐条列出被纠正的误引与陷阱。

一、先说这篇站在哪儿

前面几篇讲的是「原子核长什么样」。这篇讲的是「原子核动起来是什么样」。

这两个问题是两套家伙事。算基态、算能级、算谱因子，用的是对角化、迭代、微扰那一类办法——本质上都是把一个静态的东西解出来。可原子核实际上一半以上的时间在动：它在振、在裂、在撞、在把中子从一个核搬到另一个核。你要算这些，得换一套完全不同的机器。

这套机器叫「含时方法」。核心就一句话：别去解能量最低的那个状态，去解状态怎么随时间变。

听起来只是换个解法，其实不是。它换掉的是整个思路的立足点，换来的东西很值钱，代价也很疼。这篇文章要讲的，就是这套工具箱里到底有哪几件扳手、每件能拧什么螺丝、以及每件分别在哪儿会打滑。

和之前几篇的分工，我先说清楚，省得重复。

20260702T141003 那篇讲的是不含时的方法总论——壳模型、耦合簇、IMSRG、格点 EFT 那一整排，那是「静止的原子核」的地图。20260624T195611、20260702T140026、20260922T210910 这三篇讲裂变这件事本身——裂变机理、朗之万方程、亚铅区的不对称岛、碎片纠缠。20260512T215230 那篇讲用机器学习加速计算。

这篇不碰上面任何一块。这篇只干一件事：把「含时」这两个字单独拎出来，问它买到了什么、要付什么账、以及为什么这个家族里会有十几套看起来互相打架的方法。裂变、熔合、转移只是这套工具箱的试验场，不是主角。

二、「含时」这两个字，到底买到了什么

先说清楚它的卖点，不然没法理解为什么这么多人前赴后继。

第一件：它天然守规矩。 一个含时的方程解出来，粒子数、总能量、角动量会自己守恒——因为这些量都跟哈密顿量对易，方程一写出来就自动带上了。你不用去证明，它是白送的。这一点看着不起眼，但后面会看到，很多先进方法恰恰在这里栽跟头。

第二件：它能走出小振幅。 传统的线性响应（RPA 那一套）只能算「轻轻推一下、核小幅度晃一晃」。可裂变是核被拉成两个碎片，熔合是两个核并成一个——这都不是小晃。含时方法不需要这个前提，它直接让密度矩阵往前走，走多远都行。

第三件：一体耗散是白送的。 核撞核的时候，动能会变成内部的激发。平均场框架里，这个过程不是加上去的，是自然而然出来的：核子各自在自己的平均场里跑，场又由所有核子共同决定，这个互相拖拽本身就构成耗散。

第四件：时间本身成了分辨率。 前面几篇讨论裂变时反复出现一个麻烦——很多东西发生在「一瞬间」，可你只看初态和末态，中间过程是黑箱。含时方法把中间过程摊开了：哪一步密度开始变化、哪一步断裂、哪一步中子开始流动，都能看。这篇后面会看到，正是靠这个，才有人发现原子核在短时间尺度上有一种带随机性的密度涨落 [35]。

代价呢？下面几节一件件说。我先把最疼的一句放这儿：含时方法最擅长描述的，是「最可能发生的那条路」；它不擅长描述的，是「事情有多分散」。 这句话几乎解释了本文后面所有的方法分支。

三、最小的一套：从 TDHF 走到 TDHFB

最早的那一件扳手叫 TDHF（含时哈特里-福克）。方程就一行：

``iℏ ∂ρ/∂t = [ h[ρ], ρ ]``

ρ 是一体密度矩阵，h[ρ] 是平均场哈密顿量，它自己又是 ρ 的函数——所以这是个非线性方程，得一步一步往前推。

历史上第一脚是 1976 年迈出来的。Bonche、Koonin、Negele 三个人做了个一维的计算 [1]，把两个「核板」撞在一起，看它们会不会并起来。那不是真实三维，但足以证明这条路走得通。六年后 Negele 写了一篇很厚的综述 [2]，把平均场怎么同时管结构和动力学讲了一遍。到 1979 年，Błocki 和 Flocard 用 TDHF 算出了巨共振 [3]——这是第一次有人用含时方法算集体振动。这一支后来的整体框架，被 Nakatsukasa、Matsuyanagi、Matsuo、Yabana 四个人整理成了一篇标准综述 [4]，想系统入门从那里进最省事。

TDHF 有一个很显著的性质，值得单独说：它不允许量子隧穿。 原因很直白——它推的是一条经典式的轨迹，没有「同时走两条路再叠加」这回事。所以垒下熔合（能量不够、靠隧穿过去）它算不了；裂变鞍点之前那段（核还没决定要裂，但已经在「穿」了）它也算不了。

这个缺陷催生了大量后续工作，本文第六、第九节都要回到这里。

第二件扳手是配对。原子核里有配对关联——核子喜欢两两抱成对。要给 TDHF 加上这个，就得到 TDHFB（含时哈特里-福克-博戈留博夫）。多了一组 U、V 变换的方程，密度的概念从「占据数」变成「准粒子」。

在往下走之前，先记一个已经有人算过、而且挺有意思的结果。Loebl、Maruhn、Reinhard 用一个新工具——维格纳分布函数——去查 TDHF 到底有没有把系统推到平衡 [53]。他们算 ¹⁶O+¹⁶O 和 ⁹⁶Zr+¹³²Sn，看出两件事：一是初始相空间里那个「体」几乎不参与熔合，熔合靠的是表面那一层；二是熔合反应里所有观测量的均值守恒得很好，这被解读成存在记忆效应——系统记得自己是从哪儿来的。结论摆得很清楚：看到了强耗散，但拿不到完全平衡的证据。

但配对加进来，麻烦也跟着来了。这一节先记住一件事：TDHF 描述「最可能的路」，加上配对也不改变这一点。 真正的账在下一节。

四、配对那一环：TDHF+BCS 为什么会漏粒子

这是这个家族里最诚实、也最值得看的一段技术史，出自 Scamps、Lacroix、Bertsch、Washiyama 四人的一篇对比研究 [10]。

配对有一套正规的实现叫正则基 TDHFB（Cb-TDHFB），Ebata、Nakatsukasa、Inakura、Yoshida、Hashimoto、Yabana 六个人搭起来，并且做了系统的线性响应计算 [9]。但更多人用的是下面这条便宜得多的路。

事情是这样的。给 TDHF 加配对，最省事的做法是走 BCS 近似——不引入准粒子，只在原来那些单粒子态上给个占据概率。便宜、好写、跑得快。很多人就这么干。

然后他们发现：这么做会违反连续性方程。 连续性方程是「粒子数不会凭空产生或消失」的数学形式。违反了会怎样？你会看到粒子密度出现非物理的振荡——密度在某些地方自己涨起来又自己掉下去，好像在无中生有。

那怎么补救？一个土办法是把占据数冻住，演化过程中不让它变。振荡是压下去了，但作者原话说得挺重：BCS 近似还有别的毛病，「使它用于反应动力学是可疑的」。

然后他们做了件更狠的事：拿 TDHFB 的结果去跟两粒子薛定谔方程的数值精确解对照。结论分两半：

  • 短时间：发射率差不多，能对上。
  • 长时间：总发射量可以差很多。

差别的来源作者讲得很清楚——不是数值误差，是 HFB 理论本身的平均场性质。平均场再精确也只是一个平均，长时间累积下来，平均和真实就走岔了。

同一篇里还有个细节：TDHFB 的 U、V 方程组不唯一，因为博戈留博夫变换有规范自由度。他们试了几种写法，结论是「常见的那种写法恰好也是最高效的」，算是个安慰。

再补一件相关的活。Avez、Simenel、Chomaz 用 TDHFB 去找「配对振动」——把配对强度当成一个集体自由度，用双中子转移算符去激发它 [11]。他们在 ¹⁸O、²⁰O、²²O 和 ⁴²Ca、⁴⁴Ca、⁴⁶Ca 上做，用 SLy4 泛函。结果有个干净的对比：巨配对振动在氧同位素里看到了，在钙同位素里没有。 这是个能拿来做判据的结论。

配对这件事到这里还没完。第五节的四条局限里，有两条跟它有关。

五、TDDFT 的四个硬伤

现在把话说圆。上面讲的 TDHF、TDHF+BCS、TDHFB，其实都是同一个东西的不同实现，统称 TDDFT（含时密度泛函理论）。它的边界在哪里，最清楚的一份清单来自北大那组人 [31]。我把它拆成四条，用自己的话讲。

硬伤一：小振幅的时候，它其实就是 RPA。

你把幅度调到很小，TDDFT 会退化成 RPA。RPA 是什么？是把核的激发当成「一个核子从下面跳上去」这种最简单的组态（专业上写 1p-1h）。问题是真实的巨共振里，一个核子跳上去之后还会带动旁边一堆核子重新排布——也就是 2p-2h、3p-3h 甚至更多。这些 TDDFT 里没有。

后果很具体：巨共振的强度分布，算出来会比实验窄。 缺的那些高阶组态就是让峰会「胖」起来的东西。这个「胖出来的宽度」有个专门名字叫扩展宽度（spreading width），TDDFT 系统性低估它 [31]。

硬伤二：大振幅的时候，它算不出「分散」。

这个是本文反复出现的那句话。TDDFT 能算出一体耗散，能算出平均的动能损失，但它没有集体空间里的量子涨落。结果就是：裂变碎片的质量分布宽度、电荷分布宽度，它算不出来 [31][8]。它能告诉你「大概裂成什么样」，不能告诉你「有多少种裂法、各占多少」。

硬伤三：隧穿算不了。

前面讲过 TDHF 没有量子隧穿。这条在裂变里的具体表现是：鞍点之前那段——核已经变形了但还没越过势垒——它处理不了。垒下熔合同理 [48][49]。

绕开这条的办法叫 DC-TDHF（密度约束 TDHF）：在演化的每个瞬间，把质子密度和中子密度「钉」在瞬时值上，做一次静态能量最小化，这样就能从 TDHF 轨迹里反推出一条核-核相互作用势，再拿这个势去做垒下外推 [45][46]。它在 Ca+Ca [47] 和 ¹⁶O+²⁰⁸Pb [46] 上都能对上数据，还顺带给出一个纯 TDHF 自己给不出的信息：势垒高度会随碰撞能量变——能量越高，势垒越接近「冻结密度」极限 [48]。

硬伤四：长期传播的稳定性。

非线性方程一步一步往前推，误差会累积。这条路在数值上比解一个静态问题要脆弱得多。这一条不是理论上的硬伤，是工程上的现实，但它决定了另一件事：为什么含时方法算得慢，以及为什么很多工作只能跑很短的时间。

四条摆在一起，你会看到一个很清楚的分岔：

  • 想修硬伤一（宽度），得把高阶组态加回来。
  • 想修硬伤二（分散），得给集体空间加量子涨落。
  • 想修硬伤三（隧穿），得跳出「一条轨迹」的框架。
  • 硬伤四没有理论解，只能靠算力和算法。

第六节开始，就是四条不同的路。

六、超越平均场的四条路

这一节的骨架直接来自北大的那份分类 [31]（他们的方法就是在这四条路之外又开了一条，见第八节）。我按顺序讲。

第一条：TDDM，含时密度矩阵理论。

思路是「把 BBGKY 层级往上截一段」。什么叫 BBGKY？简单说，一体密度矩阵的演化方程里会出现二体量，二体方程里会出现三体量，三体里出现四体，无穷往上。截断在最省事的那一级——只算一体密度矩阵和二体关联矩阵，把三体及以上全扔掉——就是 TDDM [12][13]。

它真的算出了新东西。Tohyama 用 TDDM 算 ⁴⁰Ca 的巨四极共振，发现强度分成了两个成分，而主成分的扩展宽度比实验小[12]。这个「还是偏小」很关键——说明往前走了一步，但没走够。

TDDM 的代价是维度。二体关联矩阵的元素数按单粒子态数的四次方涨。所以实际用的时候还得再截一刀，最常见的是只留下时间反演配对的那些关联 [31]。

第二条：TDRPA，含时随机相位近似。

这条路的源头是 Balian 和 Vénéroni 的变分原理 [14][15]。他们的想法很巧：不要求你算得准，要求你要算的那个观测量准——按观测量去优化状态的演化。用它来算一体观测量的涨落，就是 TDRPA。

Simenel 在 2011 年把它用到了两个 ⁴⁰Ca 的碰撞上 [16]。结果是：在深部非弹碰撞里，BV 变分原理给出的涨落远大于标准 TDHF 的结果，而且更接近实验。他同时还算出了质子和中子数的关联——据这篇的说法，这是第一次在量子微观框架里算出这个关联。

但 TDRPA 有个硬约束：它现在只能处理对称碰撞[20]。非对称的算不了。

第三条：随机 TDHF / SMF。

Ayik 在 2008 年提了个很朴素的招 [17]：既然一条确定性轨迹算不出涨落，那我就跑一群轨迹，每条给一个带随机波动的初始密度矩阵，跑完之后看这批结果的分散程度。这就是随机平均场（SMF）。

这条路走得很远，从熔合的涨落-耗散 [18]、到核子交换的量子描述 [21]、一路做到了好几个真实的转移反应体系 [20]。这套方法的系统阐述见 [19]；如果你想把更早的随机方法全景也看一遍，Abe、Ayik、Reinhard、Suraud 四个人 1996 年那篇两百页的长篇综述是绕不过去的 [22]。它的强项是能处理非对称碰撞——正好补上 TDRPA 的短板。

代价也很明确，而且这一点必须讲清楚：出来的轨迹是统计平均的，不是相干平均[31]。换句话说，它算的是「一群互不相干的可能性」，不是「一个量子态的所有分支叠加在一起」。这两件事在量子力学上不一样。另外它的预测力取决于你怎么采样——换一套采样方案，结果会变 [31]。

第四条：gd-TDGCM，广义含时生成坐标法。

这条最讲究，单独用一整节讲。

先给出这四条路的账本：

方法修的是哪条硬伤主要代价
TDDM扩展宽度（硬伤一）二体矩阵维度太大，还得再截
TDRPA集体空间涨落（硬伤二）只能算对称碰撞；仍属小振幅
随机 TDHF集体空间涨落（硬伤二）统计平均而非相干平均；依赖采样方案
gd-TDGCM涨落 + 保留 TDDFT 的耗散基态反馈被忽略，严格能量守恒失守

七、gd-TDGCM 内部：两个范式，一笔没还的账

先说生成坐标法（GCM）是什么。它是这么个想法：你别去解一个复杂的大方程，你先生成一堆「典型形状」，然后把真实的核态写成这些形状的叠加，叠加系数由变分确定。对裂变来说，「典型形状」就是核拉长的各种程度。

传统 GCM 是静态的——形状是固定的几何配置。gd-TDGCM 把「含时」加了进去：生成的那些状态本身也在随时间演化，权重函数也依赖时间[26][27]。结果就是一堆 TDDFT 轨迹被相干地叠加起来。

这件事的意义在哪？它同时抓住了两头：TD-DFT 带来的一体耗散（轨迹本身是耗散的），和GCM 带来的集体空间量子涨落（轨迹之间是相干叠加的）。这是它比前面三条路都更完整的地方。

它已经被推到了非相对论版和相对论版 [26][27]，也在真实体系上出过一批结果：质量参数 [25]、耗散效应 [23]、总动能分布 [24]。

现在说那笔没还的账。这是 [31] 里点出来的，我觉得是整篇最诚实的一句批评：

基态是沿着 TDDFT 轨迹各自独立演化的，所以组态混合对基态的反馈被忽略了；结果是严格的能量守恒不再有保证。

拆开讲。本来是「一堆组态混合在一起」决定了基态长什么样。但因为每条轨迹都自己独立往前走，谁也不回头看别人，混合带来的反作用就丢了。丢了反作用，能量就守不严。

这也解释了为什么第八节那个新方法要把「能量守恒」当成一个卖点专门写进摘要——因为在这个家族里，它真的不是自动的。

顺便说一个容易混的地方。这几年有一批工作是把 TDGCM 用在裂变上、署名里都有 Vretenar 和 Nikšić，容易被当成同一批人同一件事。实际上作者组合在变：Zhao、Nikšić、Vretenar 三人做过一组 [23][24]，其中一篇还多一个 Zhou [25]；Li 牵头的相对论版是五人 [26]，后来加了 Zhao 变成六人 [27]。引用的时候得看清是哪一篇。

还有一条独立的路子值得记一笔：Regnier、Dubray、Schunck 那组人走的是「用有限元解 TDGCM 积分方程」的工程路线 [29][30]，并且在镄同位素上做了从不对称裂变到对称裂变的演化 [29]。他们的题目和「寿命」无关，是分布形状的演化——这个别记混。

以及，Verriere 和 Regnier 写过一篇专门讲 TDGCM 的综述 [28]，想快速入门的话从那里进最省事。

八、CI-TDDFT：把基态和激发放在同一块地基上

前面六条路有一个共同点——它们都是在 TDDFT 外面打补丁。基态还是那个平均场基态，补丁加在激发或者涨落那一头。

2026 年 5 月挂出来的一篇，走的是另一条路：换地基[31]。

作者是北大核物理与核技术全国重点实验室的 Y. P. Wang、B. Li、P. W. Zhao、J. Meng，加上萨格勒布大学的 D. Vretenar 和 T. Nikšić，一共 6 人（这一条我特意去原文数过，二手摘要里有写成 5 人的）。

他们干的事叫「组态相互作用含时密度泛函」（CI-TDDFT）。关键改动是：不再用一条平均场波函数，而是把多体波函数展开在一组正交的「多粒子组态」上——从一个共同的正交单粒子基出发，把 np-nh 激发（一个核子跳上去、两个跳上去……）都摆进来当基矢。然后让单粒子基和展开系数在同一个变分原理（Dirac-Frenkel 含时变分原理）下一起演化。

为什么这一步重要？回到硬伤一：TDDFT 缺的就是 2p-2h 及以上的组态。CI-TDDFT 把它们显式放进基里了。而且因为基态和激发共享同一套基，「基态反馈」这件事不再是问题。

他们拿 ⁵⁸Ni 和 ⁶⁰Ni 的等标量巨单极共振做了验证，用的是相对论点耦合泛函 PC-PK1 加一个单极对相互作用。三个结果值得记：

  1. 能量和粒子数都守恒，演化过程中相对偏差在 4×10⁻⁴ 以内。
  2. 跟常规 TDDFT 比，强度分布更宽。
  3. 同时，主峰位置跟 TDDFT 保持接近。

第 2、3 条合起来看才有意思：不是把峰挪歪了才换来宽度，是峰没动、宽度长出来了。作者把宽度归给「价空间的组态混合」，并且说这暗示单极振荡跟别的集体自由度耦合上了。

九、精确 TDGCM 加投影：绕开那个近似

2026 年 8 月还有一篇，处理的是 TDGCM 内部一个很技术、但影响很实在的近似 [32]。

TDGCM 传统上有两个省事的做法：一是高斯重叠近似（GOA），把两个不同形状的核态之间的重叠积分近似成一个高斯；二是半经典 WKB，用来处理隧穿。这两个近似在裂变寿命计算里都会带误差。

Lau、Scamps、Tan、Bernard 四个人把这俩都绕开了，做了一套「带投影的精确 TDGCM」，用在自发裂变上。这篇目前是预印本（arXiv:2608.02220v2），还没查到正式期刊版，引用时得注意。

十、从平均场一步跨到从头算：含时耦合簇

前面所有方法有个共同的天花板：它们都建立在能量密度泛函之上。泛函是个有效理论，参数是拟合出来的，所以你没法问你算的结果「离真实答案有多远」。

含时耦合簇（TDCC）走的是另一条垂直方向：不换水平方向的底座，换垂直方向的高度。

它的做法是把含时多体薛定谔方程拿耦合簇的方式解——用指数型的关联算符作用在参考态上，然后往前演化。这跟静态耦合簇（前面 20260702T141003 那篇讲过）是同一套关联的哲学，只是加了时间。

这条路的历史比很多人想的早。2012 年 Pigg、Hagen、Nam、Papenbrock 四个人就把框架搭起来了 [33]。

真正让我觉得这一块值得单独写一节的是 2026 年的两篇。

第一篇是 Bonaiti、Balos、Godbey、Hagen、Papenbrock、Woodward 六个作者，发表在 *Phys. Rev. C* 2026 年第 113 卷第 2 期，文号 024312 [34]。他们的思路很直白：解含时的 A 体薛定谔方程，把含时跃迁矩阵元记录下来，然后做傅里叶变换——谱信息就从频域里出来了。

验证用的是 ⁴He 和 ¹⁶O 的电偶极跃迁：把响应函数分布的矩跟等价的静态框架对照，差别可以忽略。

然后他们看¹⁶O 和 ²⁴O 里的质子和中子密度怎么随时间变，看到软偶极共振和巨偶极共振作为质子/中子的集体振荡从计算里自己长出来了——注意，这是从头算出来的，不是拟合出来的图像。

再然后他们做了一件挺有意思的事：在强电场下看这个核会怎样。 结论是系统的行为变混沌（chaotic）。这个结果跟另一条完全不同的路——含时相对论平均场算巨共振——能够呼应，后者早就发现同位旋振荡是混沌的 [51]。两套完全不同的理论、相隔二十多年、在「集体运动在某个条件下会失稳」这件事上对上了，我觉得这是这一块最耐看的一个点。这个「同位旋振荡不稳」的图像在相对论 RPA 那边也能找到呼应——Ma、Van Giai、Wandelt、Vretenar、Ring 五个人的计算显示，Dirac 海态对等标量单极模的影响相当可观，而相对论 RPA 给出的单极能量跟含时相对论平均场的结果是一致的 [52]。

第二篇更值得说，因为它直接改写了第二节里那句「时间成了分辨率」的具体内容。

作者是同一个组的三人：Bonaiti、Hagen、Papenbrock [35]。他们从哈特里-福克密度出发做含时耦合簇演化，去看密度涨落的动力学。发现的东西用他们自己的话说挺兴奋：存在短时、小振幅、短程、带随机性的涨落。 而且这些涨落在几十飞米每光速（fm/c）这个尺度上的振荡，跟含时平均场「电影」里看到的定性上相似。

这一条的意义在于：它第一次在从头算层面确认了平均场方法里那些看起来像数值噪声的东西，有一部分是真的。

顺带说一句时间步长的量级，帮读者建立尺度感：早期那篇 TDCC 的工作里，步长取 0.05 fm/c 的时候能量漂移约 10⁻¹⁰ MeV 量级，放大到 1 fm/c 就掉到 10⁻² MeV 量级 [33]。含时演化的精度是被步长死死管着的，这也是硬伤四的具体样子。

十一、双时格林函数与 Kadanoff-Baym：另一条完全不同的路

到这里我要换一个坐标系。前面十节讲的，全都是「波函数/密度矩阵跟着时间走」。现在讲一条从根上就不一样的路。

非平衡格林函数 / Kadanoff-Baym 方程。

区别在哪？看那个函数长什么样。前面所有方法里，你追踪的是一个单时对象——t 时刻的密度矩阵 ρ(t)。KB 方法追踪的是双时对象：G(t, t′)，两个时间都在变。

这带来一个直接后果：变量数翻倍。空间本来就有两个坐标，现在时间也有两个。计算成本上去了，但换来三样东西：

  1. 关联是内建的。 双时格林函数本身就能携带关联信息，不像平均场那样必须额外加。
  2. 总能量守恒——注意，是包括相互作用势能在内的总能量。这一点比前面那些「能量守不守得住要看情况」的方法干净。
  3. 频谱是有限宽的。 密度分布不再是「在某个能量上的一个尖峰」，而是有宽度的谱函数的频率积分。

它还自带一个漂亮的检验：取半经典极限，KB 方程会退回玻尔兹曼方程。 也就是说，它不是一个孤立的新理论，而是把玻尔兹曼那套输运理论当成了自己的一个极限情形。

Danielewicz 在 1984 年最早做了数值应用 [35]。后来 Köhler 和 Morawetz 写了一篇双时格林函数的关联理论 [37]——注意这篇文章的正式标题是「多体系统里用双时格林函数的关联」，而不是常被引成的「Kadanoff-Baym 方程」，引的时候别改题名。

这条路为什么没变成主流？看那个「变量数翻倍」就知道。核物理里真实的三维、非均匀体系，双时再加上两个空间坐标，算起来极其昂贵。所以这一支的成果主要集中在核物质、一维模型、以及半导体和等离子体这些能忍受成本的体系里。

它现在的价值更像是一把标尺：给那些便宜的方法提供一个「上限在哪」的参照。

十二、洛伦兹积分变换：用复能量换束缚态

这一节讲一个特别巧的招，值得单独说，因为它的思路是「绕过问题」而不是「解决问题」。

先说它要绕开什么。你想算一个响应函数 R(ω)——比如核被光子打了一下，各个能量上吸收多少。这个东西本质上是连续的：光子可以把核打到任意激发能，包括把它打散。连续谱的东西最难算，因为你没法把无穷多个末态一个个列出来。

洛伦兹积分变换（LIT）的招数是：别直接算 R(ω)，去算它的一个积分变换。

``L(ω₀, Γ) = ⟨Ψ_L | Ψ_L⟩, 其中 (H − z)|Ψ_L⟩ = Ô|Ψ₀⟩, z = E₀ + ω₀ + iΓ``

关键就在那个 z 是复数。因为带了虚部 iΓ，这个方程里的 |Ψ_L⟩ 是束缚的——可以用正常的束缚态方法算出来。而 Γ 就是那个「分辨率参数」。

也就是说：用一个复能量的技巧，把一个连续谱问题变成了一堆束缚态问题。 积分变换算完之后，再通过反演把 R(ω) 恢复出来。

代价当然有：Γ 不能无限小。Γ 越小分辨率越高，但 |Ψ_L⟩ 越接近非束缚、越算不出来。 所以这套方法天生带着一个有意的平滑，反演过程本身也会引入不确定度。这也意味着「分辨率」和「可算性」之间有一个必须自己选的权衡点。

这套方法的系统综述是 Efros、Leidemann、Orlandini、Barnea 四人写的 [38]，原始文献是 1994 年那篇 [39]。

2025–2026 年这条路有两个新进展，值得记。

一个是把 LIT 和神经网络量子态接起来 [40]。作者是 Parnes、Barnea、Carleo、Lovato、Rocco、Zhang 六人。他们的亮点是能给出理论不确定度的精确上界——这在从头算响应函数里是很稀缺的东西，因为大多数方法只能给个结果，说不出离真值多远。

另一个是把它和量子-经典混合计算接起来 [41]。作者是 Du、Yang、Liu、Yang、Vary 五人，发表在 *Phys. Lett. B* 2026 年第 878 卷、文号 140538。署名单位是中国科学院近代物理研究所（兰州）、先进能源科学与技术广东省实验室（惠州）、劳伦斯伯克利国家实验室应用数学与计算研究部、爱荷华州立大学物理与天文系。他们做的是「LIT 加哈密顿量输入方案加三种提取协议」，在 ¹⁹O 上演示了束缚态谱和响应函数。

十三、连续谱与晕核：NCSMC 那条线

前面十一节的路线大致是一条：从平均场往上加东西。这里要讲的是完全另一条——从头算方法自己往连续谱方向扩。

代表是 NCSMC（无核壳模型含连续谱）。它把「核子在一个束缚的壳模型基里」这件事松开，允许核子跑到连续谱里去，于是可以用来描述晕核——那种有一个或两个中子在核外飘着、几乎要散掉的核。

2026 年有一篇系统综述 [42]，作者是 Navrátil、Quaglioni、Hupin、Gennari、Kravvaris 五人，发表在 *Particles* 第 9 卷第 2 期、文号 57。他们的输入是手征 2N+3N 力，没有别的拟合参数。综述里覆盖了 ⁶He、⁸B、¹¹Be、¹⁵C 这些弱束缚态和共振，也讨论了 ¹⁴C(n,γ)¹⁵C 怎么产生 ¹⁵C。

这篇要不要算「含时」？我的判断是：不算，但必须放进来。

为什么不算？它算的是静态的散射态和共振——波函数不随时间演化，只是不再束缚而已。为什么又要放进来？因为它和含时方法在解同一类问题：核快要散架的时候是什么样。TDHF 的「没有隧穿」和 NCSMC 的「有连续谱」是同一堵墙的两面。读者如果把这两条线放在一起看，会更容易明白为什么含时方法在这些体系上格外吃力。

写这篇的时候要提醒一句：这类方法容易被笼统归到「含时理论」里，这是错的，引用时要分清。

十四、量子计算上的实时演化

最后一条新路：把时间演化本身交给量子计算机。

思路不难理解。经典计算机算含时演化，成本随体系大小指数涨；量子计算机本身就是个量子系统，让它照着演一遍就行。难的是资源。

2026 年 9 月发表在 *Quantum* 期刊上的一篇 [43]，作者是 Spagnoli、Lissoni、Roggero 三人，单位在特伦托大学、ECT*、FBK 和 INFN 特伦托。他们做的是一阶量子化下的资源估计：拿完整的领头阶无π EFT 哈密顿量，比较「乘积公式」和「量子信号处理」两种策略。

给出的结论挺清楚：这个哈密顿量的时间演化，资源随粒子数是多项式的，随单粒子基态数是取对数的。跟同一个哈密顿量在二阶量子化下的早期工作比，是指数级的改善 [43]。落到具体数字：低能核散射的模拟「几十万个 T 门 + 几百个逻辑量子比特」，因此作者判断简单的核反应有可能落在早期容错量子平台的能力范围内。

还有一条配套的路子，是先把波包准备好 [44]——Rule 和 Stetcu 两个人提的「用含时波包处理反应」的方案，给量子计算做输入。

对这条路的判断我放在第十七节。

十五、中国这一边的账

这个家族里中国团队的份量不轻，而且集中在最近三四年。我按单位捋一遍，因为这几家的贡献很容易被混起来。

北京大学（核物理与核技术全国重点实验室）——P. W. Zhao 与 J. Meng 组，是这篇里出现频率最高的。

  • 线性链 α 团簇的含时相对论密度泛函：[54]，Ren、Zhao、Meng 三人。
  • 用含时协变密度泛函算多核子转移：[55]，Zhang、Vretenar、Nikšić、Zhao、Meng 五人，*PRC* 109 卷 2 期 024614。顺带一个可用的结论：他们用 PC-PK1 和 DD-ME2 两个相对论泛函、以及前人用 SLy5 的 Skyrme 泛函，得到的转移总截面差不多；但跟半经验的 GRAZING 模型比，微观模型跟数据吻合得好得多。
  • gd-TDGCM 的相对论版：[26][27]。
  • CI-TDDFT：[31]。
  • 多核子转移里的退激发与纠缠：[56]，Yang、Zhang、Vretenar、Li、Nikšić、Zhao、Meng 七人，*PRC* 113 卷 4 期 044626。这篇做了一件挺新的事：把含时协变密度泛函跟统计退激发模型接起来，然后用互信息量化退激发过程怎么削弱碎片之间的初始量子纠缠。

中国科学院理论物理研究所（ITP-CAS）——出现在两处。一处是 Zhao、Nikšić、Vretenar、Zhou 的 gd-TDGCM 质量参数那篇 [25]，Shan-Gui Zhou 在这里。另一处是第 56 号文献里 Y. C. Yang 的署名含 ITP-CAS。

中国科学院近代物理研究所（IMP-CAS，兰州）——第 41 号那篇 [41]，LIT 与量子-经典框架。

广东这一边——先进能源科学与技术广东省实验室（惠州）出现在 [41] 的署名里。

注意一个坑： 第 55 号 [55]、第 56 号 [56] 和 [26][27] 的作者组合互不相同（分别是 5 人、7 人、5 人、6 人），但都以北大署名打头。做检索的时候如果把「北大 + Vretenar + 含时」当一个特征去搜，会捞出一大串，必须逐篇对作者数。这是我在核实阶段花时间最多的一个地方。

十六、这些方法各自的病：一张对照表

把前面散落的东西收进一张表。

方法时间处理守不守恒主要病灶
TDHF单时、一条轨迹粒子数、能量守恒无隧穿；无集体涨落；宽度偏低
TDHF+BCS单时连续性方程被违反（TU 密度振荡）冻住占据数可缓解；作者自己说「可疑」
TDHFB单时、准粒子守恒（方程唯一）长时发射量跟精确解差很多
TDDM单时、一体+二体需再截断二体矩阵维度爆炸
TDRPA单时、按观测量优化依观测量而定只能对称碰撞；本质仍小振幅
随机 TDHF单时、系综统计守恒统计平均 ≠ 相干平均；依赖采样方案
gd-TDGCM生成态与权重都含时严格能量守恒不保证基态反馈被忽略
CI-TDDFT正交多体基含时能量与粒子数守恒（4×10⁻⁴）基空间大小；尚在轻核验证阶段
精确 TDGCM+投影生成态含时—计算成本；仍为预印本
TDCC单时、A 体薛定谔方程依截断和要求而定步长敏感；算力
Kadanoff-Baym双时总能量守恒（含势能）变量翻倍，成本极高
LIT用复能量换束缚态—Γ 与分辨率/可算性的权衡；反演不确定度

这张表我建议倒着读——最后一行的「变量翻倍换守恒」，恰恰是前面所有便宜方法都想绕开的那件事。 这个家族几十年的技术史，本质上就是「用多少守恒和多少涨落，去换多少可算性」的一串折中。

十七、我的判断

三条。

第一，这个家族最值得记的不是某一套方程，而是那条分界线。

「能算最可能的路」和「能算事情有多分散」是两件事，而且后者难得多。TDHF 到 TDHFB 是加了个物理自由度，但那条分界线没动；TDDM、TDRPA、随机 TDHF、gd-TDGCM、CI-TDDFT 这五条路，全都是为了跨过这条线才出现的，而且各有各的过法。你看裂变碎片质量分布的宽度、看深部非弹的电荷分散、看巨共振的胖瘦——这些都是同一堵墙在不同场景下的样子。

第二，中国团队在这个家族里的位置，跟前面几篇不一样。

在裂变微观刻画那篇里，中国团队主要在应用和验证端。但在含时方法这个方法论层面，北大那一组的角色是提供方法本身——gd-TDGCM 的相对论推广、CI-TDDFT，都是往外走的通用工具，不是给某个体系做的专用计算。这个位置更难拿到，也更值得留意。

第三，关于量子计算，我想说句不那么热的。

第 44 号那篇给出的资源估计（几十万 T 门、几百逻辑比特）确实比二阶量子化好了指数级，但它对应的是领头阶无π EFT 下的简单核反应。这个哈密顿量本身已经把核力简到不能再简了。从「简单的核反应」到「真实的裂变动力学」，中间隔着的是三个核子力和连续谱，不是几个数量级的问题。

我判断这条路在中期内的真正价值，是给经典含时方法当参照——告诉你哪些结果是算法的局限、哪些是物理的真相，而不是替代经典计算去算裂变。如果哪天有人拿量子计算机算出 ²⁴⁰Pu 的裂变碎片分布，那才是真的换代。在那之前，这一支的成绩应该按「资源估计」而不是「物理结果」来读。

这套工具箱里没有一件是「更好」的——它们只是各自选择了在守恒、涨落、可算性这三样里放弃哪一个。谁把这三个字同时拿住，谁就赢了这个方向。

参考文献

正文中的引用标注 [N] 与本列表第 N 条一一对应。

  1. Bonche P., Koonin S. E., Negele J. W., "One-dimensional nuclear dynamics in the time-dependent Hartree-Fock approximation", *Physical Review C* 13(3), 1226–1258 (1976). DOI 10.1103/PhysRevC.13.1226. 3 位作者。TDHF 早期奠基性数值工作之一（一维核板碰撞）。
  2. Negele J. W., "The mean-field theory of nuclear structure and dynamics", *Reviews of Modern Physics* 54(4), 913–1015 (1982). DOI 10.1103/RevModPhys.54.913. 单一作者，署名麻省理工学院核科学实验室（MIT LNS）。〔注意：常被引成 1983 年，实际是 1982 年。〕
  3. Błocki J., Flocard H., "Calculations of giant resonances with time-dependent Hartree-Fock", *Physics Letters B* 85(3), 163–166 (1979). DOI 10.1016/0370-2693(79)90568-9. 2 位作者。首次用含时方法算集体振动。
  4. Nakatsukasa T., Matsuyanagi K., Matsuo M., Yabana K., "Time-dependent density-functional description of nuclear dynamics", *Reviews of Modern Physics* 88(4), 045004 (2016). DOI 10.1103/RevModPhys.88.045004；arXiv:1606.04717. 共 4 位作者。含时密度泛函在核物理里的标准综述。
  5. Simenel C., Umar A. S., "Heavy-ion collisions and fission dynamics with the time-dependent Hartree-Fock theory and its extensions", *Progress in Particle and Nuclear Physics* 103, 19–66 (2018). DOI 10.1016/j.ppnp.2018.07.002；arXiv:1807.01859. 共 2 位作者（C. Simenel；A. S. Umar）。
  6. Simenel C., "Nuclear quantum many-body dynamics (2nd edition)", *The European Physical Journal A* 61(7), 181 (2025). DOI 10.1140/epja/s10050-025-01624-3；arXiv:2506.04261. 单一作者。题献 Paul Bonche（TDHF 应用于核系统的先驱）。〔注意：这篇是 Simenel 一个人写的；常被与第 5 条那篇 Simenel & Umar 的 PPNP 综述混为一谈。〕
  7. Simenel C., "Nuclear quantum many-body dynamics", *The European Physical Journal A* 48(11), 152 (2012). DOI 10.1140/epja/i2012-12152-0. 单一作者。第 6 条的第 1 版，两版不要互引为同一篇。
  8. Schunck N., Regnier D., "Theory of Nuclear Fission", *Progress in Particle and Nuclear Physics* 125, 103963 (2022). DOI 10.1016/j.ppnp.2022.103963；arXiv:2201.02719. 共 2 位作者。
  9. Ebata S., Nakatsukasa T., Inakura T., Yoshida K., Hashimoto Y., Yabana K., "Canonical-basis time-dependent Hartree-Fock-Bogoliubov theory and linear-response calculations", *Physical Review C* 82(3), 034306 (2010). DOI 10.1103/PhysRevC.82.034306；arXiv:1007.0785. 共 6 位作者（正则基 TDHFB，Cb-TDHFB）。
  10. Scamps G., Lacroix D., Bertsch G. F., Washiyama K., "Pairing dynamics in particle transport", *Physical Review C* 85(3), 034328 (2012). DOI 10.1103/PhysRevC.85.034328；arXiv:1202.2681. 共 4 位作者。署名单位：GANIL（Scamps、Lacroix）；华盛顿大学西雅图分校（Bertsch）；布鲁塞尔自由大学 + GANIL（Washiyama）。〔注意：「连续性方程被违反」针对的是 TDHF+BCS 的常用表述；原文明确指出 HFB 的方程是唯一的、守恒律是满足的。别把两件事并成一句。〕
  11. Avez B., Simenel C., Chomaz P., "Pairing Vibrations Study with the Time-Dependent Hartree-Fock-Bogoliubov theory", *Physical Review C* 78(4), 044318 (2008). DOI 10.1103/PhysRevC.78.044318；arXiv:0808.3507. 共 3 位作者。研究体系 ¹⁸,²⁰,²²O 与 ⁴²,⁴⁴,⁴⁶Ca，SLy4 泛函。〔注意结论方向：巨配对振动在氧同位素中观察到了，在钙同位素中没有。〕
  12. Tohyama M., "Giant quadrupole resonances in time-dependent density-matrix theory", *Nuclear Physics A* 657(4), 343–354 (1999). DOI 10.1016/S0375-9474(99)00327-9. 单一作者，署名杏林大学医学部（Kyorin University School of Medicine）。结论：⁴⁰Ca 的 GQR 强度分成两个成分，主成分的扩展宽度比实验小。
  13. Tohyama M., "Double giant dipole resonances in time-dependent density-matrix theory", *Physics Letters B* 484(3-4), 231–235 (2000). DOI 10.1016/0370-2693(00)00685-7；arXiv:nucl-th/0003034. 单一作者。结论：¹⁶O 的 DGDR 因朗道阻尼而有很大宽度；⁴⁰Ca 的 DGDR 更接近谐和。
  14. Balian R., Vénéroni M., "Extension Through Time Smoothing of the Time Dependent Mean Field Theory", *Annals of Physics* 135(2), 270–311 (1981). DOI 10.1016/0003-4916(81)90156-1. 共 2 位作者。署名：萨克莱（Saclay）与奥赛核物理研究所（IPN Orsay）。
  15. Balian R., Vénéroni M., "Correlations and fluctuations in static and dynamic mean field approaches", *Annals of Physics* 216(2), 351–430 (1992). DOI 10.1016/0003-4916(92)90181-K. 共 2 位作者。BV 变分原理的系统阐述，TDRPA 的直接来源。
  16. Simenel C., "Particle-number fluctuations and correlations in transfer reactions obtained using the Balian-Vénéroni variational principle", *Physical Review Letters* 106(11), 112502 (2011). DOI 10.1103/PhysRevLett.106.112502；arXiv:1011.2293. 单一作者。发表时署名单位为澳大利亚国立大学（ANU）核物理系与法国 CEA Saclay IRFU/SPN 双署名。结论：深部非弹碰撞中 BV 涨落远大于标准 TDHF、且更接近实验；首次在量子微观框架内给出质子-中子数关联。
  17. Ayik S., "A stochastic mean-field approach for nuclear dynamics", *Physics Letters B* 658(4), 174–179 (2008). DOI 10.1016/j.physletb.2007.09.072；arXiv:0709.2686. 单一作者，署名田纳西理工大学（Tennessee Technological University）。随机平均场（SMF）的原始文献。
  18. Ayik S., Washiyama K., Lacroix D., "Fluctuation and dissipation dynamics in fusion reactions from a stochastic mean-field approach", *Physical Review C* 79(5), 054606 (2009). DOI 10.1103/PhysRevC.79.054606. 共 3 位作者。
  19. Lacroix D., Ayik S., "Stochastic quantum dynamics beyond mean field", *The European Physical Journal A* 50(6), 95 (2014). DOI 10.1140/epja/i2014-14095-8. 共 2 位作者。SMF 的综述性阐述。
  20. Ayik S., Arik M., Yilmaz O., Umar A. S., "Stochastic Mean-Field Theory and Applications to Multinucleon Transfer and Kinetic Energy Dissipation Processes in Heavy-Ion Collisions", arXiv:2509.18914 (2025). 共 4 位作者。〔纠错：这一篇常被误记为含 Lacroix——SMF 系列早期确实多为 Lacroix 与 Ayik 合作，但本篇署名是 Ayik、Arik、Yilmaz、Umar 四人。〕截至本文核实日未见正式期刊版。文中明确指出的两条限制：TDRPA 现行形式不能用于非对称碰撞；SMF 的轨迹是统计平均（statistically averaged）而非相干平均。
  21. Ayik S., Yilmaz O., Yilmaz B., Umar A. S., Gokalp A., Turan G., Lacroix D., "Quantal description of nucleon exchange in a stochastic mean-field approach", *Physical Review C* 91(5), 054601 (2015). DOI 10.1103/PhysRevC.91.054601. 共 7 位作者。
  22. Abe Y., Ayik S., Reinhard P.-G., Suraud E., "On stochastic approaches of nuclear dynamics", *Physics Reports* 275(2-3), 49–196 (1996). DOI 10.1016/0370-1573(96)00003-8. 共 4 位作者。随机方法的经典长篇综述。
  23. Zhao J., Nikšić T., Vretenar D., "Time-dependent generator coordinate method study of fission: Dissipation effects", *Physical Review C* 105(5), 054604 (2022). DOI 10.1103/PhysRevC.105.054604. 共 3 位作者。
  24. Zhao J., Nikšić T., Vretenar D., "Time-dependent generator coordinate method study of fission. II. Total kinetic energy distribution", *Physical Review C* 106(5), 054609 (2022). DOI 10.1103/PhysRevC.106.054609. 共 3 位作者。
  25. Zhao J., Nikšić T., Vretenar D., Zhou S.-G., "Time-dependent generator coordinate method study of fission: Mass parameters", *Physical Review C* 101(6), 064605 (2020). DOI 10.1103/PhysRevC.101.064605. 共 4 位作者（多出 S.-G. Zhou，署名中国科学院理论物理研究所）。
  26. Li B., Vretenar D., Nikšić T., Zhao P. W., Meng J., "Generalized time-dependent generator coordinate method for small and large amplitude collective motion", *Physical Review C* 108(1), 014321 (2023). DOI 10.1103/PhysRevC.108.014321；arXiv:2304.13369. 共 5 位作者。
  27. Li B., Vretenar D., Nikšić T., Zhao J., Zhao P. W., Meng J., "Generalized time-dependent generator coordinate method for induced fission dynamics", *Frontiers of Physics* 19(4), 44201 (2024). DOI 10.1007/s11467-023-1381-4；arXiv:2309.12564. 共 6 位作者（比第 26 条多出 J. Zhao）。〔注意：第 26 与第 27 条作者数不同，别并成一句引。〕
  28. Verriere M., Regnier D., "The Time-Dependent Generator Coordinate Method in Nuclear Physics", *Frontiers in Physics* 8, 233 (2020). DOI 10.3389/fphy.2020.00233. 共 2 位作者。TDGCM 的入门综述。
  29. Regnier D., Dubray N., Schunck N., "From Asymmetric to Symmetric Fission in the Fermium Isotopes within the Time-Dependent GCM Formalism", *Physical Review C* 99(2), 024611 (2019). DOI 10.1103/PhysRevC.99.024611；arXiv:1810.08402. 共 3 位作者。署名单位：奥赛核物理研究所（IPN Orsay）、法国 CEA DAM DIF、劳伦斯利弗莫尔国家实验室（LLNL）。〔注意：这篇讲的是分布形状的演化，不是自发裂变寿命。〕
  30. Regnier D., Verrière M., Dubray N., Schunck N., "FELIX-1.0: A finite element solver for the time dependent generator coordinate method with the Gaussian overlap approximation", *Computer Physics Communications* 200, 350–363 (2016). DOI 10.1016/j.cpc.2015.11.013. 共 4 位作者。TDGCM 的有限元求解器。
  31. Wang Y. P., Li B., Vretenar D., Nikšić T., Zhao P. W., Meng J., "Configuration-interaction time-dependent density functional theory for nuclear dynamics", arXiv:2605.21037 (2026). 共 6 位作者。署名单位：北京大学物理学院核物理与核技术全国重点实验室（Wang、Li、Zhao、Meng）；克罗地亚萨格勒布大学理学院物理系（Vretenar、Nikšić）。〔纠错：二手摘要中有把作者数写成 5 人的，arXiv 原文署名是 6 人。〕应用：PC-PK1 相对论点耦合泛函 + 单极对相互作用，⁵⁸Ni 与 ⁶⁰Ni 的等标量巨单极共振；能量与粒子数守恒，相对偏差在 4×10⁻⁴ 以内；强度分布比常规 TDDFT 更宽，主峰位置保持接近。截至本文核实日未见正式期刊版。
  32. Lau N.-W. T., Scamps G., Tan P., Bernard R. N., "An exact time-dependent generator coordinate method with projection for spontaneous fission reactions", arXiv:2608.02220 (2026). 共 4 位作者。避开高斯重叠近似（GOA）与半经典 WKB。截至本文核实日为预印本（v2）。
  33. Pigg D. A., Hagen G., Nam H., Papenbrock T., "Time-dependent coupled-cluster method for atomic nuclei", *Physical Review C* 86(1), 014308 (2012). DOI 10.1103/PhysRevC.86.014308；arXiv:1205.4002. 共 4 位作者。TDCC 的框架奠基。
  34. Bonaiti F., Balos C., Godbey K., Hagen G., Papenbrock T., Woodward C. S., "Computing nuclear response functions with time-dependent coupled-cluster theory", *Physical Review C* 113(2), 024312 (2026). DOI 10.1103/scn2-dnyg；arXiv:2510.19940. 共 6 位作者。报告号 LLNL-JRNL-2012532。正式发表日期 2026-02-09。〔纠错：部分数据库条目的作者列表缺 C. S. Woodward，只列 5 人；arXiv 原文与橡树岭方面页面均为 6 人。〕方法：解含时 A 体薛定谔方程 → 记录含时跃迁矩 → 傅里叶变换提取谱。验证：⁴He 与 ¹⁶O 的电偶极跃迁。结果：¹⁶O 与 ²⁴O 中出现软偶极与巨偶极共振的质子/中子集体振荡图像；在强电场下系统行为变混沌（chaotic）。
  35. Bonaiti F., Hagen G., Papenbrock T., "Dynamics of Density Fluctuations in Atomic Nuclei", *Physical Review Letters* 137(7), 072502 (2026). DOI 10.1103/xdr8-4tp3；arXiv:2604.12089. 共 3 位作者。〔纠错：多处二手来源仍按「arXiv 预印本」引用，该工作已正式发表于 PRL 137 卷 7 期 072502。〕结论：存在短时、小振幅、短程、带随机性的密度涨落；其在几十 fm/c 尺度上的振荡与含时平均场结果定性相似。
  36. Danielewicz P., "Quantum theory of nonequilibrium processes, 1", *Annals of Physics* 152(2), 239–304 (1984). DOI 10.1016/0003-4916(84)90092-7. 单一作者，署名劳伦斯伯克利实验室核科学部。Kadanoff-Baym 方程的早期系统性数值应用。〔注意：期刊是 *Annals of Physics*，不是 *Nuclear Physics*；同题分 part 1 与 part 2 两篇。〕
  37. Köhler H. S., Morawetz K., "Correlations in Many-Body Systems with Two-time Green's Functions", *Physical Review C* 64(2), 024613 (2001). DOI 10.1103/PhysRevC.64.024613；arXiv:nucl-th/0102059. 共 2 位作者。署名：亚利桑那大学（Köhler）、德累斯顿马克斯·普朗克复杂系统物理研究所（Morawetz）。〔注意：正式题名是「多体系统里用双时格林函数的关联」，不是常被引成的「Kadanoff-Baym 方程」，引时勿改题。〕
  38. Efros V. D., Leidemann W., Orlandini G., Barnea N., "The Lorentz integral transform (LIT) method and its applications to perturbation-induced reactions", *Journal of Physics G: Nuclear and Particle Physics* 34(12), R459 (2007). DOI 10.1088/0954-3899/34/12/R02；arXiv:0708.2803. 共 4 位作者。署名：库尔恰托夫研究所（莫斯科）；特伦托大学 + INFN 特伦托；希伯来大学拉卡物理研究所 + 华盛顿大学核理论研究所（INT）。〔纠错：部分数据库条目把页码记作「R459-2456」，尾段是录入错误；正式页码自 R459 起。〕
  39. Efros V. D., Leidemann W., Orlandini G., "Response functions from integral transforms with a Lorentz kernel", *Physics Letters B* 338(2-3), 130–133 (1994). DOI 10.1016/0370-2693(94)91355-2；arXiv:nucl-th/9409004. 共 3 位作者。LIT 方法的原始文献。
  40. Parnes E., Barnea N., Carleo G., Lovato A., Rocco N., Zhang X., "Nuclear responses with neural-network quantum states", arXiv:2504.20195 (2025). 共 6 位作者。用神经网络量子态配合 LIT 计算偶极响应与光电吸收截面，并可给出理论不确定度的精确上界。截至本文核实日未见正式期刊版。
  41. Du W., Yang Y., Liu Z., Yang C., Vary J. P., "Quantum-classical computational framework for many-fermion response and structure", *Physics Letters B* 878, 140538 (2026). DOI 10.1016/j.physletb.2026.140538；arXiv:2602.08357. 共 5 位作者。署名单位：中国科学院近代物理研究所（兰州）；先进能源科学与技术广东省实验室（惠州）；劳伦斯伯克利国家实验室应用数学与计算研究部；爱荷华州立大学物理与天文系。内容：LIT 加哈密顿量输入方案加三种提取协议，在 ¹⁹O 上演示束缚态谱与响应函数。
  42. Navrátil P., Quaglioni S., Hupin G., Gennari M., Kravvaris K., "Halo Nuclei from Ab Initio Nuclear Theory", *Particles* 9(2), 57 (2026). DOI 10.3390/particles9020057；arXiv:2604.02612. 共 5 位作者。署名单位：TRIUMF + 维多利亚大学（Navrátil）；劳伦斯利弗莫尔国家实验室（Quaglioni、Kravvaris）；IJCLab 奥赛（Hupin）；美因茨大学（Gennari）。综述 ⁶He、⁸B、¹¹Be、¹⁵C；讨论 ¹⁴C(n,γ)¹⁵C。〔注意：这是含连续谱的静态方法，不属于含时方法，不可归入本文的含时家族。〕
  43. Spagnoli L., Lissoni C., Roggero A., "Quantum Simulation of Nuclear Dynamics in First Quantization", *Quantum* 10, 2200 (2026). DOI 10.22331/q-2026-09-02-2200；arXiv:2507.22814. 共 3 位作者。署名单位：特伦托大学、ECT*、FBK CIT、INFN 特伦托。〔纠错：常被按「arXiv 预印本」引用，该工作已正式发表于 Quantum 期刊 10 卷 2200（2026-09-02）。〕结论：领头阶无π EFT 哈密顿量的时间演化，资源随粒子数为多项式、随单粒子基态数为对数；相比同模型在二阶量子化下的早期工作为指数级改善；低能核散射估计需要「几十万个 T 门 + 几百个逻辑量子比特」。
  44. Rule E., Stetcu I., "A time-dependent wave-packet approach to reactions for quantum computation", arXiv:2603.26881 (2026). 共 2 位作者。为量子计算准备含时波包的方案。
  45. Umar A. S., Oberacker V. E., "Heavy-ion interaction potential deduced from density-constrained time-dependent Hartree-Fock calculation", *Physical Review C* 74(2), 021601 (2006). DOI 10.1103/PhysRevC.74.021601. 共 2 位作者，署名范德堡大学。DC-TDHF 方法原始文献。
  46. Umar A. S., Oberacker V. E., "Density-constrained time-dependent Hartree-Fock calculation of ¹⁶O+²⁰⁸Pb fusion cross-sections", *The European Physical Journal A* 39(2), 243–247 (2009). DOI 10.1140/epja/i2008-10712-5. 共 2 位作者。
  47. Keser R., Umar A. S., Oberacker V. E., "Microscopic study of Ca + Ca fusion", *Physical Review C* 85(4), 044606 (2012). DOI 10.1103/PhysRevC.85.044606；arXiv:1202.3932. 共 3 位作者。署名：里泽大学 + 范德堡大学（Keser）；范德堡大学（Umar、Oberacker）。
  48. Umar A. S., Simenel C., Oberacker V. E., "Energy dependence of potential barriers and its effect on fusion cross-sections", *Physical Review C* 89(3), 034611 (2014). DOI 10.1103/PhysRevC.89.034611；arXiv:1401.7230. 共 3 位作者。文中明确指出 TDHF「只能用于垒上能量的熔合，因为该理论不允许多体隧穿」，但其熔合阈值能很好地预言动态绝热势垒顶端能量。
  49. Umar A. S., Oberacker V. E., Horowitz C. J., "Microscopic sub-barrier fusion calculations for the neutron star crust", *Physical Review C* 85(5), 055801 (2012). DOI 10.1103/PhysRevC.85.055801；arXiv:1203.5565. 共 3 位作者。署名：范德堡大学（Umar、Oberacker）；印第安纳大学 + 田纳西大学 + 橡树岭国家实验室（Horowitz）。
  50. Vretenar D., Lalazissis G. A., Behnsch R., Pöschl W., Ring P., "Monopole giant resonances and nuclear compressibility in relativistic mean field theory", *Nuclear Physics A* 621(4), 853–878 (1997). DOI 10.1016/S0375-9474(97)00192-9；arXiv:nucl-th/9612042. 共 5 位作者。署名：萨格勒布大学；亚里士多德大学塞萨洛尼基分校；慕尼黑工业大学物理系。〔注意：该文结论是核物质压缩模量 K_nm ≈ 250–270 MeV；文中另有「早期 GCM 工作曾给出 ≈300 MeV」的历史值，两者不可混引。〕
  51. Vretenar D., Paar N., Ring P., Lalazissis G. A., "Nonlinear dynamics of giant resonances in atomic nuclei", *Physical Review E* 60(1), 308–319 (1999). DOI 10.1103/PhysRevE.60.308；arXiv:nucl-th/9809003. 共 4 位作者。发表于 *Physical Review E*（非核物理期刊，易被误搜）。结论：等标量单极振荡为规则运动，同位旋振荡为混沌动力学。
  52. Ma Z.-Y., Van Giai N., Wandelt A., Vretenar D., Ring P., "Isoscalar compression modes in relativistic random phase approximation", *Nuclear Physics A* 686(1-4), 173–186 (2001). DOI 10.1016/S0375-9474(00)00523-6. 共 5 位作者。结论：Dirac 海态对等标量单极模影响可观；相对论 RPA 的单极能量与含时相对论平均场结果吻合。
  53. Loebl N., Maruhn J. A., Reinhard P.-G., "Equilibration in the time-dependent Hartree-Fock approach probed with the Wigner distribution function", *Physical Review C* 84(3), 034608 (2011). DOI 10.1103/PhysRevC.84.034608；arXiv:1104.1941. 共 3 位作者。研究体系 ¹⁶O+¹⁶O 与 ⁹⁶Zr+¹³²Sn。结论：初始相空间体几乎不参与熔合；熔合反应中观测量均值守恒，指示存在记忆效应；观察到强耗散，但没有完全平衡的证据。
  54. Ren Z. X., Zhao P. W., Meng J., "Dynamics of the linear-chain alpha cluster in microscopic time-dependent relativistic density functional theory", *Physics Letters B* 801, 135194 (2020). DOI 10.1016/j.physletb.2019.135194；arXiv:2001.02834. 共 3 位作者，署名北京大学。
  55. Zhang D. D., Vretenar D., Nikšić T., Zhao P. W., Meng J., "Multinucleon transfer with time-dependent covariant density functional theory", *Physical Review C* 109(2), 024614 (2024). DOI 10.1103/PhysRevC.109.024614；arXiv:2401.06539. 共 5 位作者。结论：PC-PK1、DD-ME2 两个相对论泛函与 SLy5 Skyrme 泛函给出的多核子转移总截面相近；相比半经验的 GRAZING 模型，微观模型与数据吻合好得多。
  56. Yang Y. C., Zhang D. D., Vretenar D., Li B., Nikšić T., Zhao P. W., Meng J., "De-excitation effects on entanglement in multi-nucleon transfer reactions", *Physical Review C* 113(4), 044626 (2026). DOI 10.1103/fl3l-fb88. 共 7 位作者。署名单位：北京大学核物理与核技术全国重点实验室（Yang、Li、Zhao、Meng）；中国科学院理论物理研究所（Zhang、Yang）；萨格勒布大学（Vretenar、Nikšić）。正式发表 2026-04-27。内容：把含时协变密度泛函与统计退激发模型结合，用 Shannon 熵与互信息量化退激发对碎片间初始量子纠缠的削弱。

核实与纠错记录

#常见说法实际情况
1CI-TDDFT（arXiv:2605.21037）是 5 位作者arXiv 原文署名 6 位：Y. P. Wang、B. Li、D. Vretenar、T. Nikšić、P. W. Zhao、J. Meng。部分摘要式数据库条目只列前 5 位。[31]
2CI-TDDFT 正文里引的「YangYC2026PRC」指向一篇小作者数的 PRC 论文指的是 Yang Y. C. 等 7 位作者、*PRC* 113(4), 044626 (2026)、DOI 10.1103/fl3l-fb88，即多核子转移退激发与纠缠那篇。[56]
3含时耦合簇响应函数那篇是 5 位作者arXiv 原文与橡树岭方面页面均为 6 位，包含 C. S. Woodward。部分数据库条目漏掉此人。[34]
4「Dynamics of density fluctuations in atomic nuclei」只是 arXiv 预印本该工作已正式发表：*PRL* 137(7), 072502 (2026)，DOI 10.1103/xdr8-4tp3。二手来源大量滞后。[35]
5随机平均场（SMF）那篇综述有 Lacroix 参加该综述（arXiv:2509.18914）署名是 S. Ayik、M. Arik、O. Yilmaz、A. S. Umar 四人，没有 Lacroix。SMF 系列早期确实多为 Lacroix 与 Ayik 合作，容易顺手续错。[20]
6「含时核多体动力学」综述是 Simenel 与 Umar 合写的那是 2018 年 *PPNP* 103, 19 那篇 [5]。2025 年的第 2 版是 Simenel 单一作者，*EPJA* 61(7), 181。[6][7]
7Simenel 的 BV 涨落 PRL 有多位合作者*PRL* 106(11), 112502 (2011) 是单一作者；发表时署名单位为澳大利亚国立大学核物理系与法国 CEA Saclay IRFU/SPN 双署名。[16]
8LIT 综述页码是「R459–R2456」页码尾段是数据库录入错误。正式页码自 R459 起。[38]
9Negele 的综述发表于 1983 年实为 1982 年，*RMP* 54(4), 913–1015。[2]
10含时相对论平均场算巨单极给出的核物质压缩模量是 300 MeV该文结论是 K_nm ≈ 250–270 MeV。文中提到的「≈300 MeV」是早期 GCM 工作的历史结论，不是该文的结论。[50]
11巨共振非线性动力学那篇发表在核物理期刊上发表在 *Physical Review E* 60(1), 308 (1999)，按核物理期刊关键词检索会漏掉。结论是「等标量规则、同位旋混沌」，方向别写反。[51]
12gd-TDGCM 那两篇（PRC 108,014321 与 Front. Phys. 19,44201）是同一批作者作者数不同：前者 5 位（Li、Vretenar、Nikšić、P. W. Zhao、Meng），后者 6 位（多出 J. Zhao）。引两条时不能并成一句「同一作者组」。[26][27]
13gd-TDGCM 的质量参数那篇是 3 位作者是 4 位，多出 S.-G. Zhou（中国科学院理论物理研究所）。[25]
14「TDHFB 违反连续性方程」违反连续性方程的是 TDHF+BCS 的常用表述，会造成非物理的粒子密度振荡。原文明确说 HFB 的方程是唯一的、守恒律是满足的。[10]
15巨配对振动在氧和钙同位素里都看到了只在氧同位素里看到，钙同位素里没有。[11]
16镄同位素那篇 TDGCM 算的是自发裂变寿命算的是从不对称裂变到对称裂变的分布演化，期刊版题名即为「From Asymmetric to Symmetric Fission in the Fermium Isotopes」。[29]
17双时格林函数那篇的题名叫「Kadanoff-Baym equations」正式题名是「Correlations in Many-Body Systems with Two-time Green's Functions」，作者 2 位（Köhler、Morawetz）。引时勿改题。[37]
18Danielewicz 1984 那篇发在 *Nuclear Physics*发在 *Annals of Physics* 152(2), 239 (1984)，且同题分 part 1 / part 2 两篇。[36]
19量子计算核动力学那篇只是预印本已正式发表于 *Quantum* 10, 2200 (2026)，DOI 10.22331/q-2026-09-02-2200（2026-09-02）。[43]
20NCSMC 晕核综述属于含时方法不属于。它是含连续谱的静态方法（束缚态+散射态，波函数不随时间演化）。本文只把它作为「同一堵墙的另一面」对照，不归入含时家族。[42]
21光锥积分变换（LIT）里 Γ 越小越好Γ 越小分辨率越高，但那道复能量方程解出来的 Ψ_L 就越接近非束缚态、越算不出来。分辨率与可算性之间存在必须显式选择的权衡。[38][39]
22北大＋Vretenar＋含时的工作可以用一套作者列表概括不能。同一署名结构下至少存在 5 人 [55]、6 人 [31]、7 人 [56] 三套不同的作者列表；且 [26] 是 5 人、[27] 是 6 人。必须逐篇对作者数。
23搜索结果里的「August 24, 2026」是论文日期那是 ar5iv 等镜像站页面的渲染日期，不是论文发表或提交日期。必须回 arXiv 原文或期刊页取日期。
24以下站点可以直接当文献来源引用不可以：tomesphere.com、ar5iv.arxiv.org、ar5iv.labs.arxiv.org、arxiv-vanity.com、xr-scholar.cn、scixplorer.org、exaly.com、handwiki.org、llmpedia.net、indico.ectstar.eu。前几类是 arXiv 的聚合/镜像渲染站，后几类是学术社交或会议幻灯片页；其中 llmpedia.net 关于 Baym–Kadanoff 的条目通篇无可核查处，明显由模型生成。这些一律只能当线索，必须回到 arXiv 原文、期刊官网或 CrossRef 落地核实。
25用 arXiv 检索接口按标题批量反查很方便在本次核实环境中，arXiv API 的检索型查询（search_query=ti:...）因 Accept 头不匹配返回 HTTP 406 Not Acceptable；用 id_list= 的按编号直查则正常。标题未知时的批量反查需退回 CrossRef 书目检索。