反应计算共用的一块带 30% 误差的地基
*引用格式说明:* 正文中的 [N] 与文末参考文献列表第 N 条一一对应;文末另附「核实与纠错记录」,逐条列出被纠正的误引与陷阱。
一、先说这篇站在哪儿
一个 (d,p) 实验做完了。测出一串角分布点。然后把每个点除以理论算出的截面,取平均,得到一个数,叫谱因子 S。二十多年来,大家拿这个数讨论一件事:核子里那些"单粒子成分",到底还剩几成没被关联搅散。
现在换个角度问一句:如果那个"理论截面"本身带 30% 的不确定度呢?
那 S 少掉的量里,有多少是核真在淬灭,有多少是算出来的?这个问题不好答,而且麻烦的是它不只是算得准不准的问题。它意味着你不可能用一把自己都量不准的尺子去判断东西短了多少。
这话不是我编的。这篇文章要讲的,就是这个 30% 从哪来、为什么二十年过去还没解决、以及最近五年有一批人拿贝叶斯方法在干什么。
先把跟知识库里既有文章的分工说清楚,省得重复。
《弹性散射与光学势核数据基础》(20260825T190342)讲的是最近又有人拟合出了哪几个新的全局势——⁷Li 的、α 核的、氘核 100 MeV/u 的。它的落点是"谁做了什么"。本篇不列这个清单。
《核反应三探针》(20260903T110000)讲的是转移、敲出、库仑激发三种探针怎么配合。《从头算走到 A≈100》(20260919T100434)讲的是从头算结构计算里误差棒怎么进化——那是结构侧。
本篇讲的是反应侧的一块地基。它问的是三个别的问题:
其一,光学势这件事为什么天生就不准——不是算得糙,是数据本身约束不住。
其二,把"一个数"换成"一片云"之后,账变成了什么样——贝叶斯方法到底改了什么,以及它自己新添了哪些麻烦。
其三,这块地基的误差往下游传了多远——传到了谱因子、淬灭因子、中子俘获率,最后传到"这个实验值不值得做"。
二、光学势是什么:一张"我不打算算的部分"的欠条
先说清楚这东西是干嘛的。
一个质子打到一个核上。严格说,这要把核里两百来个核子和这个质子一起放进薛定谔方程解——没人算得动,也没必要算。
光学模型的做法是:只盯住"质子还待在弹道位置上"这个通道,把其余所有可能发生的事——核子被撞飞、核被激发、各种多步过程——打包成一个复数势的虚部。实部管"它受到的推拉",虚部管"流走出去多少"。这就是光学势 [1]。
它是有效势,所以必然依赖靶的质量数、电荷数、和束流能量。写下来大约是这样:
``U(r) = V(r) + i W(r) + V_Coulomb(r)实部 虚部 库仑项``
三个部分各自用几个参数描述形状和深浅,加在一起几十个参数。这些参数不是从核力推出来的,是拿弹性散射数据拟合出来的 [2][3]。
形状用的是 Woods-Saxon 形式——一个在核表面附近平滑过渡的"碗",虚部则通常写成体积项加表面项两项,意思是吸收主要发生在核表面 [4]。局域势要配上能量依赖,才能把实验对上;为了省事,大家习惯再套一层 Perey-Buck 型的高斯非定域修正 [5]。
这套写法用了七十多年,很好用。问题也在这里:这套形式本身没有理论推导做依据,是照着"核是个有表面的东西"这个直觉写的。第十一节会讲,有人从微观计算里发现这个直觉有一处是错的。
三、那 30% 是怎么拆出来的
现在来看那个 30%。
出处是一篇中国团队的文章,《从弹性转移反应提取谱因子》,发在《中国科学:物理学 力学 天文学》上 [6]。作者十个人:李志宏、李云居、苏俊、颜胜权、王友宝、郭冰、南丁、李二涛、甘林、柳卫平。这是一支做低能核反应测量和核天体反应的队伍。
他们做的事很朴素:把从 (d,p)、(d,n) 这类反应里提取谱因子的各条误差来源,一项一项列出来。列出来是这样:
| 误差来源 | 量级 |
|---|---|
| 测量统计 | 3% |
| 靶厚 | 5% |
| 入射道光学势参数 | 15% |
| 出射道光学势参数 | 15% |
| 束缚态几何参数 | 20% |
| 氘核内部的谱因子 | 15% |
结论一句话:提取出的谱因子,典型误差大约 30%。作者还加了一句评价——这个值"非常大",会影响后续的核结构与核天体分析,所以很有必要发展误差更小的新方法 [6]。
把这张表看仔细一点。这里面有两项都叫"光学势",一个管进来的路,一个管出去的路,各 15%。单看它们俩,按平方和加起来是 21%。
但真正要命的不在这儿。 要命的是"典型误差约 30%"这个数,和"淬灭因子通常在 0.5 到 0.7 之间"这个数,是同一个量级。你在讨论一个连自己都带三成误差的量,去判断它少了两三成——这个判断本来就悬着。
我不是说这些实验白做了。我是说,这个误差账一旦摆到桌面上,后面所有关于"淬灭"的结论都得先减掉它。至于能不能减得掉,往下看。
四、弹性散射数据为什么约束不住光学势
三十年前就有人知道这事,有个专门的名字,叫 Igo 模糊。
意思是:在所谓"灵敏半径"那个位置上,不同的参数组给出的实部深度和虚部深度会交在同一点上,取值几乎一样。于是很多组完全不同的参数,都能同样好地拟合同一批角分布数据 [7]。
2019 到 2020 年之间,中国原子能科学研究院核物理所的杨磊、林承键、张英逊等人做了一件挺狠的事:把这套模糊性直接摆到贝叶斯框架下量了一遍 [7]。
他们选的体系是弱束缚核 ⁶Li 打在 ²⁰⁹Bi 上,一共十一个人署名。做法是用三种完全不同的先验分布——阈异常先验、反常阈异常先验、均匀先验——各自跑一遍贝叶斯分析,看后验分布跑到哪儿去。
结果:三种先验都能很好地重现实验的弹性散射角分布。但导出的光学势随能量变化的曲线,走法截然不同。
这句话得慢慢读。它的意思是:弹性散射数据在垒下能区根本没有足够的信息,去约束光学势参数。你要么用不带约束的传统频率方法,要么就用均匀(无偏见)先验——总之不能让先验替你决定结论 [7]。
他们还顺手用转移反应 ²⁰⁸Pb(⁷Li,⁶He)²⁰⁹Bi 去碰同一个问题,看 ⁶He+²⁰⁹Bi 这个中子晕体系在出射道的光学势。在合适先验下,他们确认了"反常阈异常"确实存在——这等于给"色散关系还适用不适用"这个老问题又添了一个问号 [7]。
这里有个细节值得单独说。这篇文章本身是给贝叶斯方法提警告的。 同一批人在国内报道里的原话是:贝叶斯是机器学习的核心算法之一,它基于已有知识预测信息,结果带有先验的偏见;对"新现象、新规律"的探索,无偏见的先验才更合理。
一个用贝叶斯的人来写贝叶斯的局限,比外面的人来写可信得多。第十三节会看到这条警告后来被印证了不止一次。
五、贝叶斯进场:把"一条曲线"换成"一片云"
同一时期,密歇根州立大学那一边走的是反方向的路:既然约束不住,那就把"约束不住"这件事本身算出来。
2018 年,Lovell 和 Nunes 那篇《用贝叶斯定理约束转移截面》是这个方向的起点 [8]。思路不复杂:不再用最小二乘找一组最优参数,而是用马尔可夫链蒙特卡罗跑出参数的后验分布——一组参数变成了一朵云。
然后把这朵云里的每一组参数都喂进反应计算,看输出散成什么样。
他们算的是氘束在 A=48 到 208 的靶上、束能 10 到 25 MeV/u 的转移反应,重点讲了 ⁴⁸Ca(d,p)⁴⁹Ca 基态这个例子。第一个结论是:转移截面的不确定度在 25% 到 100% 之间,随反应不同差很多 [8]。
还有两个副产品,比主结论更有意思。
一是绝热波近似(ADWA)比扭曲波玻恩近似(DWBA)给出的不确定度小。也就是说,模型选得越"完整",误差棒越细——这跟直觉一致,但这是第一次被量化。
二是把实验数据的误差棒缩小,不确定度确实变小,但变小的幅度跟误差缩小的幅度不成比例 [8]。这一条后来被反复验证,第七节再讲。
同年,King、Lovell、Nunes 又补了一篇,专门算 (d,p) 模型里的光学势不确定度 [9]。两年后,Lovell 等人把这一串工作汇总成一篇方法性综述,发在《Journal of Physics G》上 [10]——那一篇是这一整个方向的入门读物。
六、频率学派和贝叶斯吵的那一架,吵出了什么
到这里有个问题绕不开:传统的最小二乘(频率学派)和贝叶斯,到底该信谁?
2019 年,King、Lovell、Neufcourt、Nunes 四个人在《物理评论快报》上给了个很干脆的答案 [11]。他们做了个正面对撞:同一批弹性散射数据,同一套光学模型,两种方法各跑一遍,然后摆在一起看。
先看置信区间。贝叶斯的 95% 区间比频率方法的宽。 这跟大家心里的预期是反的——大家觉得贝叶斯带先验,应该更"保守"才对。
但接下来他们做了一件事,这一件事才是这篇论文真正的分量:经验覆盖率检验。
做法是:把置信区间想象成一套从窄到宽的刻度。然后数一数,实验数据点实际落在区间里的比例,是不是跟名义比例对得上。
结果写得很清楚:对于较大的置信区间,贝叶斯方法的覆盖率是准的;频率方法偏低。偏低意味着什么?意味着频率方法给出的 95% 区间,实际上没有那么宽——它把自己的确定性说高了 [11]。
这是这一整节里最该记住的一句话。误差棒不只是数据的函数,它还是方法的函数。同一批数据,换一套统计学,别人对它的信任程度就变了一次。
七、数据侧的天花板:角度、邻近能量、反应截面
既然不确定度这么大,能不能靠"多测点东西"压下去?Catacora-Rios 那一组人系统地试了一年多。
第一篇(2019)试了四件事 [12] ——数据改成角网格的疏密、加入邻近束能的数据、缩小实验误差棒、以及把总反应截面也塞进拟合。靶选 ⁴⁸Ca 和 ²⁰⁸Pb,然后把不确定度一路传到 ⁴⁸Ca(d,p)⁴⁹Ca 基态和 ²⁰⁸Pb(d,p)²⁰⁹Pb 基态。
结论分三句:角网格换法几乎没影响;加入邻近能量的数据能把不确定度改善最多两倍;而缩小误差棒的收益,"往往比天真预期的小得多"[12]。
第二篇(2021)换了个思路:不比方法,比哪些观测量更有信息量。他们摆出三个统计工具——主成分分析、基于导数的灵敏度分析、贝叶斯证据 [13]。
出现了一个挺尴尬的结果。先看玩具模型:灵敏度和贝叶斯证据都能明确区分开,体积吸收项在更高能量上更有用。但一换到真实的 ⁴⁸Ca 弹性散射:角分布和极化数据约束力大致相当,参数的可变性强烈依赖角度;贝叶斯证据虽然有些起伏,但算出来的贝叶斯因子不足以区分角分布与极化数据 [13]。
翻译一下:我们原以为极化数据能多给点信息,其实没有。
还有一句从实验设计角度很实用的话:能把不确定度压下来的不是"测得准一点",而是"换个角度测、或者换个能区测"——最优的角度区间随观测量变化很大 [13]。
再往后,2023 年 Hebborn、Whitehead、Lovell、Nunes 把这套方法从转移反应搬到了单中子敲出反应上,算 ¹¹Be(晕核)和 ¹²C(紧束缚核)[14]。结论是个能直接引用的门槛:从敲出截面比值反推出的结构信息,理论不确定度至少 20%(晕核)、至少 40%(紧束缚核)。同一篇里他们特别点了一句:¹¹Be 的情形下理论误差跟实验误差是同一个量级,¹²C 的情形下理论误差明显更大 [14]。
也就是说,在紧束缚核上讨论敲出得到的结构信息,理论误差比实验误差还大。这不是技术细节,这是判决力的问题。
八、从"一个势"到"一套基础设施":KDUQ、CHUQ、ELM
上面这些工作有一个共同的形状:每一个反应、每一个靶、每一个能点,单独拟合一次。这在小范围里好用,但没法当公共设施用。
2023 年,Pruitt、Escher、Rahman 三个人把两个被用得最广的全局势——Koning-Delaroche (KD) 和 Chapel Hill '89 (CH) ——补上了标定好的不确定度 [15]。
他们用的是 MCMC 做参数推断。产出的东西不是"一组最优参数",而是参数的完整后验分布 [15]。带上不确定度的版本,在文献里分别叫 KDUQ(Koning-Delaroche, Uncertainty Quantified)和 CHUQ。
这一步为什么是分水岭?因为它把光学势从"一个数"变成了"一个可以被抽样、被传播、被第三方复用的分布"。后面所有的应用都不用再自己拟合一遍——直接把后验样本喂进反应计算就行。
另一条路是 2025–2026 年出现的 ELM(East Lansing Model,东兰辛模型)。这是 Beyer 和 Nunes 两个人的工作 [16],署名在密歇根州立大学的稀有同位素束流装置(FRIB)和物理与天文系。
它跟前面几个全局势有三处不同 [16]。
第一处是形式。 ELM 用了新的中子-质子不对称度项。实部中心项写成同位旋标量项加上一个带 (N−Z)/A 的项,而且同位旋标量与同位旋矢量各有自己独立的半径和弥散度。全模型 21 个自由参数;另有一个简化版 ELM0,18 个参数,把两个分量的几何设成一样。
第二处是数据。 语料是从 EXFOR 里整出来的:17 条 (p,n) 条目共 710 个数据点,57 条弹性散射条目共 5965 个数据点,全是微分截面,束能 10 到 100 MeV,球形靶 A ≥ 40,把 (n,n)、(p,p)、(p,n) 放在同一个框架里。
第三处是统计模型,这个下面单独说。
先说一个反直觉的结果。把 (p,n) 数据单独加进去,并不显著改变参数化 [16]。这一点让作者专门写了一段来解释:信息量不是"加数据"就能用上的,得先有一个能承住这份信息的参数化形式。只有当同位旋矢量项获得独立的几何自由度——也就是允许把中子皮信息编码进光学势——(p,n) 数据才开始干活。加上之后,数据描述变好,不确定度下降。
至于外推能力,ELM 的原话是:往质子滴线和中子滴线两边推的时候,不确定性比目前在用的其他全局势更小 [16]。
但这里有一个很好的对照实验。ELM 跟 KDUQ、CHUQ 比 ¹⁰⁰Sn 到 ¹⁷⁰Sn 这一串锡同位素的质子反应截面,三家的预测差别很明显——而这些差别"并不来自实验约束或物理模型形式,来自统计模型,具体说就是似然函数的重新标度"[16]。
同一批数据、同一个物理形式,换个似然函数,结论就分出岔来。这一句值得抄下来。
九、手征 EFT 那条线:不拟合散射数据,从核力直接算
前面讲的全部是唯象路线——形状是编的,参数是拟合的。这条路线有个天生弱点:你没法把它外推到没测过的核上去,因为你不知道它在数据少的地方会走成什么样。
另一条路是微观路线:从核子-核子相互作用出发,直接算光学势。
2021 年,Whitehead、Lim、Holt 三个人在《物理评论快报》上给了这条路线上最重要的一个结果 [17]。他们用从头算多体微扰论,配上从手征有效场论(手征 EFT)来的真实核力,把核子-核散射的完整微观全局描述做了出来,覆盖 1800 个核素。而且关键在于:这套势里没有任何一个参数是拟合到核子-核散射数据上去的 [17][18]。同时它还给参数的量化不确定度。
这个性质值钱在哪?值钱在它没有被数据"喂"过。你在它身上看到的误差,是理论自己的误差,不是拟合出来的残差。
2025 年,王兵、白栋、徐毅把 WLH 势第一次用到天体物理中子俘获上 [18]。这是 r 过程需要的核心输入。做法是用 Hauser-Feshbach 统计模型算截面和反应率,然后对 26 ≤ Z ≤ 84、从稳定谷到中子滴线之间所有核素算一遍,用后向-前向蒙特卡罗方法传播不确定度,用 f_rms 偏差当 χ² 估计量。
结果里有一个很好看的结构:在 f_rms ≤ 1.56 的约束下,反应率的不确定度在同位旋不对称度 0.28 附近出现明显分界 [18]。作者对这条分界的解读很直接——这说明同位旋依赖是光学势的关键,WLH 势下一步就该重点打理这块。
这是本文里最"应用"的一处。它说明光学势的同位旋依赖不是一个可以晚点再管的技术细节,它直接决定 r 过程网络里哪一半核的反应率你根本没法算准。
同一批搞微观势的人还有另一支。Vorabbi、Gennari、Finelli、Giusti、Navrátil 用 Watson 多重散射理论的一阶展开,把核子-核密度和核子-核子 t 矩阵折起来算光学势——t 矩阵用和密度一致的手征相互作用,密度用无核壳模型(NCSM)算 [19][20][21]。他们往几个方向推:中等质量同位素(钙和镍同位素链,实验室能量 65–200 MeV)[20];核-核弹性碰撞(α 打 ¹²C 和 ¹⁶O,100–300 MeV)[19];用同一条路线反推 ⁷⁸Kr 的物质半径和中子皮 [22]。
其中最诚实的一句出现在 2026 年那篇的标题里:《微观光学势的一个唯象修正》[21]。他们直说,这类微观方法忽略了 Pauli 阻塞以外的介质修正,也没抓住高阶散射贡献,结果是系统地低估了吸收,在低能和中能区角分布描述得不好。于是加了一个能量依赖的唯象修正项进去。
"纯微观"这条路自己承认需要打补丁,这是一个很值得记住的位置。
十、非定域性:到底是"空间"还是"时间"
现在换到更根本的一层。
光学势天生是非定域的:一个核子在 r' 处受到的势,跟它现在在 r 处的波函数有关,两者不是同一点。非定域性的来源主要是反对称化和通道耦合 [23]。
但"非定域"这个词一直在混用两件完全不同的事 [24]:
空间非定域——势把波函数在一个空间范围内抹平。时间非定域——势跟能量有关,也就是跟波函数在时间上的历史有关。
七十多年来,大家是把这两件事一起塞进"能量依赖"里的。局域势要配上能量依赖才拟合得上数据;Perey-Buck 的做法则是反过来,用一个高斯型的空间非定域因子去解释那部分能量依赖 [5]。谁也没有真正分开量过。
2025 年 9 月,Stahulak 和 Holt 把这件事分开量了一次 [24]。两个人,署名在德州农工大学回旋加速器研究所和物理与天文系。
做法是拿 Perey-Buck 的假设当工具——这个假设本来就建立了"局域势的能量依赖"和"高斯型空间非定域"之间的对应关系。然后去问:从手征 EFT 算出来的微观光学势,实部和虚部的那部分能量依赖,各自对应哪一种?
答案是干净的:
微观实部光学势的能量依赖,主要来自空间非定域性。微观虚部光学势的能量依赖,是真正的时间非定域性。
他们给出的是钙同位素链上的非定域核子-核光学势,并且把 Woods-Saxon 形状参数随同位素数的变化规律看了一遍 [24]。
这里要说一句版本上的事。这篇的预印本从 v1(2025 年 9 月)到 v3(2026 年 9 月 7 日)改了很多——文件从 183 KB 涨到 3,804 KB,图从 12 张变成 13 张,摘要里的措辞也从"是真正的时间非定域性"改成了"主要是真正的时间非定域性" [24]。引这句话的时候要认版本,别引 v1。
这篇为什么重要?因为它说明把两种非定域性笼统地记为"能量依赖",在物理上是把两笔不同的账记成了一笔。而这一笔记法,正是所有唯象势的标准做法。第十一节讲的节点结构,是同一个方向上的另一刀。
十一、节点结构:Woods-Saxon 那个"表面吸收峰"可能是错的
上一节是说形式里两类东西混着记。这一节更狠:形式本身画错了。
Woods-Saxon 那套写法里,虚部通常给两项:体积项(吸收发生在核内部)加上表面项(吸收在核表面有个峰)。这个"表面吸收峰"的写法是照直觉来的——核表面密度梯度最大,那里最容易被激发出东西来 [4]。
2024 年,Arellano 和 Blanchon 在《物理评论 C》上给了这个直觉一记重击 [25]。
他们用的是动量空间的介质内折叠模型——先算一个微观的核子-核光学势,然后去看它的结构。扫描范围是束能 40 到 400 MeV、靶质量数 40 ≤ A ≤ 208 [25]。
他们发现的东西有三层。
第一层:光学势具有普适的可分离结构——径向部分和非定域部分可以分开写,而且这个分法在那么宽的能区和质量区里都成立。
第二层:非定域形状因子本质上是复数的,而且是类氢的(形状像氢原子波函数那一类),同时作用于中心项和自旋轨道项。
第三层是真正的发现:虚部的径向形状因子里,一致地出现一个节点。这个节点起的作用是显著压制表面吸收峰 [25]。
把它跟 Woods-Saxon 摆在一起看:Woods-Saxon 假设"核表面有一个吸收峰",而微观计算说"那里恰好有个节点,把峰压掉"。这两句话是互相矛盾的,不是精修级别的差别 [25]。
他们还说,那个复数径向形状因子能写成"均匀球分布与高斯形状因子及汤川项的卷积"这样的组合形式 [25]。
这篇的作者署名要弄清楚。H. F. Arellano 在智利大学物理系(FCFM,圣地亚哥);G. Blanchon 在法国 CEA DAM DIF(Arpajon),同时挂巴黎-萨克雷大学 CEA 的 LMCE 实验室(Bruyères-le-Châtel)[25]。这是智利-法国两家,不是美国组。
那这篇跟"不确定度"有什么关系?关系很直接,而且是最讨厌的那种关系。
前面所有的贝叶斯工作,量的是参数不确定度——你的形式对,只是参数不知道。节点结构这件事是模型不确定度——你的形式本身可能就错了。参数误差可以用更多数据压下去,模型误差压不下去。
Hebborn 那篇 2023 年的综述里专门有一节谈这件事,标题就是"超越两体有效势的模型不确定度"[1]。Arellano 和 Blanchon 的结果,是那一节最具体的注脚。
十二、下游到底被卡住多少
讲了这么多上游的事,该收账了。光学势的毛病到底传下去多少?
先说非定域性这一项。
2016 年,Titus、Nunes、Potel 把有限程 ADWA 推广到能直接吃非定域核子光学势的形式 [23]。他们算的是 ¹⁶O、⁴⁰Ca、⁴⁸Ca、¹²⁶Sn、¹³²Sn、²⁰⁸Pb 上的 (d,p),束能 10、20、50 MeV。
物理图像是很清楚的:非定域性让波函数在核内部振幅变小、往外移。因为转移反应的截面主要靠核外那部分波函数,往外移的结果就是——多数情况下第一个峰的转移截面变大,靶越重、能量越高,这个效应越明显 [23]。
然后是结论那句话,请慢慢读:
如果分析实验转移截面时只用局域光学势,提取出的谱因子可能因为局域近似而错最多 40% [23]。
40%。对照第三节那张表——入出射道光子势在拟合参数上各贡献 15%。这里还有一个和参数无关的 40%,来自"你把非定域的东西按局域处理了"这个动作本身。
再说一个具体的判决例子。
2025 年,蒋英驰、陈杰、B. P. Kay 等 27 人(含中方的袁岑溪等)在《Physics Letters B》上发表了一篇测 ¹⁵C 的工作 [26]。¹⁵C 是个特别干净的体系:一个价中子在 ²s₁/₂ 轨道上,¹⁴C 芯很惰,按独立粒子模型算谱因子应该接近 1。
他们用两个道测:(p,d) 和 (d,t)。
先看光学势带来的差别。用全局势和用拟合过的局域势得到的谱因子,差 18%((p,d))、差 9%((d,t))[26]。注意这里 (d,t) 的差别更小,而 (d,t) 那道本身的误差更大——原因是跟氚有关的反应里光学势的模糊性更大。整篇里谱因子或淬灭因子的标准差,约 14%((p,d))、24%((d,t)),作者明说更大的那一个来自光学势的模糊性 [26]。
他们还试了束缚态半径这个旋钮:把 r₀ 从 1.25 fm 换成推荐的 1.34 fm,结果只小约 3%;而束缚态弥散度 a₀ 的影响远小于不同光学模型参数的影响[26]。
这是个很有用的排序:在这一类反应里,光学势是最粗的那个旋钮,束缚态几何是细旋钮。
最后的结果:两个道平均,¹⁵C 基态的淬灭因子 Rs = 0.72 ± 0.14。而重离子诱导敲出给出的值是 0.96 ± 0.04——两者互相矛盾[26]。
一个说淬灭掉了三成,一个说几乎没淬灭。同一个核,同一个轨道。这篇论文的立场是:转移反应的结果更可信,因为敲出用的是包容式分析而转移是独家式的。
我不打算在这里判谁对。我要指出的是那个 0.14 的误差棒——它占了 0.72 的两成。而这个误差棒里最大的一项,正是光学势。
第三条线:把反应和结构接起来。
圣路易斯华盛顿大学 Dickhoff 那一组走的是另一条路,叫色散光学模型(DOM)——用一条色散关系,把正能区的散射数据和负能区的单粒子性质绑在一起,让同一个势同时解释两边 [27][28]。
这条路对谱因子的含义很直接。2011 年 Nguyen、Waldecker、Nunes、Charity、Dickhoff 把 DOM 用到 (d,p) 上,结论是:DOM 给出的谱因子一般比标准全局势给出的低,而且对束能的依赖更弱,跟 (e,e'p) 的结果更一致 [27]。
也就是说,你换一套更"自洽"的光学势,谱因子的数值就系统性地往下走一截。这正是"谱因子的判决力取决于你用了哪个势"的最干净的例证。
同一条线上还有几个后续:2022 年 Yoshida、Atkinson、Ogata、Dickhoff 第一次把 DOM 用到 (p,2p) 分析里 [29],发现 ⁴⁰Ca(p,2p)³⁹K 给出质子 0d₃/₂ 谱因子 0.560,而与 (e,e'p) 一致的那个 DOM 值是 0.71——两个数差得明显,而且用单粒子波函数、光学势、p-p 相互作用的选择都解释不了这个差距 [29]。他们由此说,核内 p-p 散射的描述急需改进,这对逆运动学下的分析影响很大。2025 年 Atkinson 和 Dickhoff 又用 DOM 从敲出反应里反过来学东西 [30]。
把这些摆在一起看,形状就很清楚了:上游换个势,下游的谱因子就动几成;上游换套统计,下游的误差棒就动几倍。而"淬灭"这个结论,就长在这几成和几倍之间。
十三、一个不好解释的裂缝:5–10% 对 40–100%
到这儿必须停下来处理一个矛盾。这是本文里最难写的一节,因为它没有干净的答案。
按第五节和第十四节,转移截面的不确定度是 25% 到 100%,敲出的结构信息至少 20% 到 40%。
但 2025 年 Hebborn 和 Nunes 用 KDUQ 做的那篇系统研究,给出的数完全不是这个量级 [31]。
他们算 ⁴⁸Ca 上的 (d,p),扫不同束能、不同末态单粒子性质。结论是:参数不确定度很小——68% 可信区间的半宽大约 5–10%,跟实验误差是同一个量级 [31]。
差一个数量级。这怎么解释?
那篇论文自己给了两个理由,我认为都是对的,而且都很难堪 [31]。
第一个理由:似然函数的选择。 早期那些研究用的似然写法,跟这一篇不一样。同一批数据,似然函数写法不同,后验分布的宽窄就不同。
第二个理由:相关性没算进去。 早期做法是把每个光学势独立地传播不确定度——入射道的中子势和质子势之间、它们跟出射道的质子势之间的相关性,都没算。而这一篇发现,束缚态势和光学势的不确定度之间存在相关性,而且这些相关性对正确判断转移截面的不确定度很重要;在他们的例子里,这些项不是按平方和相加的[31]。
还有一个背景差别也值得摆出来:早期那些研究是逐核逐能点单独拟合,一个靶一个能点一次;这一篇用的是已经在上百条数据集上标定过的全局势。约束的信息量根本不在一个量级。
那到底该信谁?
我的判断是这样:两个数都在说真话,但它们量的不是同一件事。 一个是"我给你一个核、一个能量,你现拟合一个势,那这个势有多不确定";另一个是"你已经有了一个全局标定过的势,把它用在 ⁴⁸Ca 上,那这个预测有多不确定"。前者是外推的代价,后者是内插的代价。
所以这篇文章里所有的百分比,都必须在"哪个势、哪套似然、算不算相关性"这三个坐标下读。脱离这三个坐标引用任何一个数,都是在散布一个没有意义的数字。
这一点在第八节那个锡同位素的对照里被直接印证了:ELM、KDUQ、CHUQ 三个都用全局标定、都描述得住弹性散射,但预测起的截面差别明显,差别来源是统计模型而不是物理。
还有一件事要记着。Hebborn 和 Nunes 那篇里有一句提醒:相对不确定度随束能升高而增大,而且不强烈依赖末态性质[31]。这跟第七节那个"换个能区测"的建议合到一起,说明能量这个坐标在不确定度管理上的分量比角度还重。
十四、往 β 不稳定核那边走:ELM 和 DNP 2025 的现场
现在回到用户提到的那条线:光学模型往 β 不稳定核的扩展。
前面讲过 ELM [16]。它要解决的就是这个问题:能不能把不确定度"控制着"往外推。它的做法是给光学势加一个能编码中子皮的不对称度项,然后用 (p,n) 数据去约束它——因为往中子-rich 那边走,同位旋矢量项才是主导 [16]。
ELM 的第一次公开亮相是在 2025 年 10 月的 APS 核物理分会(DNP)秋季会议上。Beyer 和 Nunes 的报告题目是《走向一个面向稀有同位素的全局、不确定度量化光学势》,时间是 2025 年 10 月 18 日 [32]。会议在芝加哥的希尔顿酒店办,会期 10 月 17 到 20 日。
报告里把 ELM 的四步流程讲得挺清楚:从 EXFOR 里整一份实验证据语料((n,n)、(p,p)、(p,n) 到同位旋相似态的微分截面);给物理模型设计一个唯象参数化;给似然函数设计一个能描述"模型怎么描述证据"的形式;然后做贝叶斯标定。最后拿它跟在用的最优势比,重点讲它的同位旋矢量分量有什么独到之处 [32]。
这里有个程序上的事得说清楚。 DNP 2025 那份是会议报告摘要,没有同行评议。ELM 的完整论文到 2026 年 3 月才挂到 arXiv 上 [16]。截至本文写作时(2026 年 9 月),我逐条查了 CrossRef,这篇还没有正式期刊的卷期和页码,所以文末只写 arXiv 号,不写卷期。
顺便说一句这次会议的整体气氛。DNP 2025 一共 887 个报告、160 个分会。其中有一个专门的邀请制工作坊,题目是《核实验与核理论中的系统不确定度与模型不确定度:新工具、新技术、以及合作机会》,分上下两场(编号 1WJ、2WJ),都排在 10 月 17 日上午 [33]。同一个上午还有一场关于弱束缚效应的邀请制工作坊(1WB、2WB)[33]。
把这三件事放在同一个上午:弱束缚核的物理、系统与模型不确定度的方法、以及稀有同位素的计算——这不像是巧合。这说明"光学势的不确定度"这件事,已经从某一组人的私活,变成了整个领域要在会议日程上专门排时段的事。
十五、中国这一边的账
这一节单独写,因为几条关键线索都在中国这一侧。
第一条:30% 那张表。 第三节讲的误差拆账,出自中国科学那篇 [6]。十个人署名,是低能核反应测量和核天体反应这条线上的队伍。他们给的不是抱怨,是一张可以拿去逐项改进的清单——靶厚能压、统计量能加、束缚态几何能用 Hartree-Fock 标定,只有入出射道光学势那两项最难动。他们自己给的出路也是这个方向:发展新的、误差更小的提法。
第二条:贝叶斯的第一脚。 第四节讲的那篇,是把贝叶斯方法第一次用到奇特核体系相互作用势上的工作 [7]。这件事在国内是有报道的,讲法是"首次将贝叶斯方法应用于奇特核体系相互作用势的研究"。
这篇的价值不在于它给出了多准的势,而在于它给出的结论是负面的——三种先验都拟合得上数据,说明弹性散射约束不住;贝叶斯结果强烈依赖先验,所以用它必须极度小心。国内报道把这一点总结得很到位:贝叶斯是机器学习的核心算法之一,结果带先验偏见,"对于新现象的探索,无偏见的先验知识是更为合理的先验分布"。
一个方法刚被引进来,第一批工作就在说它的边界在哪。这种自我设限在文献里不多见。
第三条:DOM 上的一次外推尝试。
2026 年 8 月 20 日,Zhao 和 Du(渤海大学)与 R. Capote(Suncoast Data Evaluation)在《物理评论 C》上发了一篇锆同位素的色散光学模型 [34]。他们做的是核子打在锆同位素上的散射,用 DOM 框架。
这篇的意义在于它接的是另一条线:把自洽的色散光学模型从钙、铅这些经典靶,推到一整条同位素链上去看系统学。而 Capote 的参与让这件事带上了核数据评价的色彩——光学势的误差要真正进入评价库,走的就是这条路 [34]。
十六、收尾:S 本来就不是一个可观测量
讲到这里,我得把一件更根本的事摆上来。前面十四节都在讲"光学势的误差有多大",但还有一个更早的问题:我们凭什么觉得那个"谱因子"是个能测的量?
2011 年,TRIUMF 的 Jennings 写了一篇《谱因子的不可观测性》[35]。
他的论证分两层。
理论层:可以构造出幺正变换,让谱因子在允许的范围内任意变化。也就是说,一个核的谱因子是几,取决于你站在哪个表象上看——它不是一个与表象无关的量。
唯象层:在现象学框架里,"场重定义"起同样的作用。结论是一句很硬的话:从实验里提取出的谱因子,更多取决于你假设的势有多强的能量依赖,而不是取决于你测到的截面 [35]。
他还指出了一个技术条件:规范不变性和 Wegmann 定理合起来,意味着在纯 DWBA 里必然给出 S=1。所以想看到 S 偏离 1,你必须走出纯 DWBA——而一走进去,势的选择就回来了 [35]。
这就是我在开头说的那件事的正式版:那把尺子不光量不准,尺子这个身份本身都成问题。
按 Jennings 的说法,谱因子只在"导出它的那个模型之内"有意义 [35]。这不是说 (d,p) 实验没价值——它当然有价值,它是我们在滴线附近少有的直接探针。而是说,把 S 当成一个核的固有属性去跨实验、跨模型地比较,是一件在方法上就不成立的事。
那出路在哪?
第十二节那三条线其实各自给了一个方向。非定域性那条说:别把非定域的东西按局域处理,那 40% 是白扔的 [23]。DOM 那条说:让反应和结构用同一个势、同一条色散关系绑在一起,这样至少自洽 [27][28]。KDUQ/ELM 那条说:别只报一个数,报一片云,让下游的人知道自己在用什么 [15][16]。
三条路都在做同一件事:把"我不知道的部分"从隐含变成显式。
这篇东西真正想说的判断只有一个。反应计算这套地基,短期内不会变成一个"准"的东西。它更可能变成一个诚实的东西——每个下游结论都带着自己那份误差,而且那份误差明确标着它来自哪里:数据不够、先验太强、似然选错、相关性没算、还是形式本身就写错了。
第十三节那段 5–10% 对 40–100% 的对照,其实就是这件事的预演。同一批数据,两个团队,两组数,一个说误差可以接受,一个说误差大到没法用。差别不在数据,在方法。 而承认这件事,比假装它不存在有用得多。
最后回到开头那句话。那把尺子量不准,这件事本身不是坏消息。坏消息是有人还在拿它的读数去做判决,而不说尺子的事。
参考文献
正文中的引用标注 [N] 与本列表第 N 条一一对应。
- Hebborn C., Nunes F. M., Potel G., Dickhoff W. H., Holt J. W., Atkinson M. C., Baker R. B., Barbieri C., Blanchon G., Burrows M., Capote R., Danielewicz P., Dupuis M., Elster Ch., Escher J. E., Hlophe L., Idini A., Jayatissa H., Kay B. P., Kravvaris K., Manfredi J. J., Mercenne A., Morillon B., Perdikakis G., Pruitt C. D., Sargsyan G. H., Thompson I. J., Vorabbi M., Whitehead T. R., "Optical potentials for the rare-isotope beam era", *Journal of Physics G: Nuclear and Particle Physics* 50(6), 060501 (2023). DOI 10.1088/1361-6471/acc348;arXiv:2210.07293. 共 29 位作者。署名单位(前五位):Hebborn 密歇根州立大学(主)+ 劳伦斯利弗莫尔国家实验室(LLNL);Nunes 密歇根州立大学 + 该校国家超导回旋加速器实验室(NSCL);Potel LLNL;Dickhoff 圣路易斯华盛顿大学;Holt 德州农工大学。本文是 2022 年 3 月 FRIB 理论联盟(FRIB-TA)专题项目「Optical Potentials in Nuclear Physics」的成果。报告号 LLNL-JRNL-840487。
- Koning A. J., Delaroche J. P., "Local and global nucleon optical models from 1 keV to 200 MeV", *Nuclear Physics A* 713(3-4), 231–310 (2003). DOI 10.1016/S0375-9474(02)01321-0. 本文简称 KD。署名:荷兰 NRG Petten;法国 CEA DAM。
- Varner R. L., Thompson W. J., McAbee T. L., Ludwig E. J., Clegg T. B., "A global nucleon optical model potential", *Physics Reports* 201(2), 57–119 (1991). DOI 10.1016/0370-1573(91)90039-O. 共 5 位作者,均属北卡罗来纳大学教堂山分校三角大学核实验室(TUNL)体系。
- Woods R. D., Saxon D. S., "Diffuse Surface Optical Model for Nucleon-Nuclei Scattering", *Physical Review* 95(2), 577–578 (1954). DOI 10.1103/PhysRev.95.577. 两位作者署名加州大学洛杉矶分校(UCLA)。Woods-Saxon 形状因子的原始出处。
- Perey F., Buck B., "A non-local potential model for the scattering of neutrons by nuclei", *Nuclear Physics* 32, 353–380 (1962). DOI 10.1016/0029-5582(62)90345-0. 〔注意:CrossRef 与 INSPIRE 的作者字段把第一作者录成了「F. Ferey」,这是旧库的录入错误,正确写法是 F. Perey;同篇第二作者是 B. Buck。Perey 效应的出处就在这篇。〕
- Li Z. H., Li Y. J., Su J., Yan S. Q., Wang Y. B., Guo B., Nan D., Li E. T., Gan L., Liu W. P., "Extracting spectroscopic factors from elastic transfer reactions", *Science China Physics, Mechanics & Astronomy* 62(3), 32021 (2019). DOI 10.1007/s11433-018-9287-9. 在线日期 2018-10-18,正式卷 2019 年 3 月。共 10 位作者。本文第三节那张 30% 拆账表的出处。
- Yang L., Lin C. J., Zhang Y. X., Wen P. W., Jia H. M., Wang D. X., Ma N. R., Yang F., Zhong F. P., Zhong S. H., Luo T. P., "Bayesian analysis on interactions of exotic nuclear systems", *Physics Letters B* 807, 135540 (2020). DOI 10.1016/j.physletb.2020.135540. 共 11 位作者。第一作者署名单位为中国原子能科学研究院核物理研究所。研究体系:⁶Li+²⁰⁹Bi(弹性散射)与 ²⁰⁸Pb(⁷Li,⁶He)²⁰⁹Bi(转移反应,考察出射道 ⁶He+²⁰⁹Bi)。
- Lovell A. E., Nunes F. M., "Constraining Transfer Cross Sections Using Bayes' Theorem", *Physical Review C* 97(6), 064612 (2018). DOI 10.1103/PhysRevC.97.064612;arXiv:1801.06096. 共 2 位作者。署名:密歇根州立大学国家超导回旋加速器实验室(NSCL)+ 该校物理与天文系。
- King G. B., Lovell A. E., Nunes F. M., "Uncertainty quantification due to optical potentials in models for (d,p) reactions", *Physical Review C* 98(4), 044623 (2018). DOI 10.1103/PhysRevC.98.044623;arXiv:1810.06129. 共 3 位作者。
- Lovell A. E., Nunes F. M., Catacora-Rios M., King G. B., "Recent advances in the quantification of uncertainties in reaction theory", *Journal of Physics G: Nuclear and Particle Physics* 48(1), 014001 (2020). DOI 10.1088/1361-6471/abba72;arXiv:2012.09012. 共 4 位作者。
- King G. B., Lovell A. E., Neufcourt L., Nunes F. M., "Direct comparison between Bayesian and frequentist uncertainty quantification for nuclear reactions", *Physical Review Letters* 122(23), 232502 (2019). DOI 10.1103/PhysRevLett.122.232502;arXiv:1905.05072. 共 4 位作者。研究体系:⁴⁸Ca(n,n) 12 MeV、⁴⁸Ca(p,p) 14 MeV 与 25 MeV。
- Catacora-Rios M., King G. B., Lovell A. E., Nunes F. M., "Exploring experimental conditions to reduce uncertainties in the optical potential", *Physical Review C* 100(6), 064615 (2019). DOI 10.1103/PhysRevC.100.064615. 共 4 位作者。注意:该文的 arXiv 预印本号是 arXiv:2112.00873,2021 年 12 月才挂出,比正式期刊版晚两年——这是作者迟挂,不是两篇不同的文章。
- Catacora-Rios M., King G. B., Lovell A. E., Nunes F. M., "Statistical tools for a better optical model", *Physical Review C* 104(6), 064611 (2021). DOI 10.1103/PhysRevC.104.064611;arXiv:2012.06653. 共 4 位作者。三个统计工具:主成分分析、基于导数的灵敏度分析、贝叶斯证据。
- Hebborn C., Whitehead T. R., Lovell A. E., Nunes F. M., "Quantifying uncertainties due to optical potentials in one-neutron knockout reactions", *Physical Review C* 108(1), 014601 (2023). DOI 10.1103/PhysRevC.108.014601;arXiv:2212.06056. 共 4 位作者。研究体系:¹¹Be(晕核)与 ¹²C(紧束缚核)。
- Pruitt C. D., Escher J. E., Rahman R., "Uncertainty-quantified phenomenological optical potentials for single-nucleon scattering", *Physical Review C* 107(1), 014602 (2023). DOI 10.1103/PhysRevC.107.014602;arXiv:2211.07741. 共 3 位作者。署名:劳伦斯利弗莫尔国家实验室(LLNL)。Rahman 的现址为田纳西大学诺克斯维尔分校,但本文按其发表时署名处理。用 MCMC 给 KD 与 Chapel Hill '89 两个全局势标定不确定度,即文献中的 KDUQ 与 CHUQ。
- Beyer K., Nunes F. M., "The East Lansing Model: a Bayesian uncertainty quantified optical potential for rare isotopes", arXiv:2603.28599 (2026)。〔预印本;截至本文写作时逐条查 CrossRef 未见正式期刊卷期与页码,故只写 arXiv 号。〕提交 2026-03-30,修订版 2026-03-31。署名:密歇根州立大学稀有同位素束流装置(FRIB)+ 该校物理与天文系。ELM 全模型 21 个自由参数,简化版 ELM0 为 18 个。标定语料:17 条 (p,n) 条目共 710 个数据点,57 条弹性散射条目共 5965 个数据点,球形靶 A ≥ 40,束能 10–100 MeV。
- Whitehead T. R., Lim Y., Holt J. W., "Global Microscopic Description of Nucleon-Nucleus Scattering with Quantified Uncertainties", *Physical Review Letters* 127(18), 182502 (2021). DOI 10.1103/PhysRevLett.127.182502;arXiv:2009.08436. 共 3 位作者。文献中简称 WLH 势。注意:arXiv 上的题名是「A global microscopic description of nucleon-nucleus scattering with quantified uncertainties」(多一个 a、措辞略异),正式发表题名以本条目为准。
- Wang B., Bai D., Xu Y., "Application of the microscopic optical potential of chiral effective field theory in astrophysical neutron-capture reactions", *Physics Letters B* 872, 140056 (2026). DOI 10.1016/j.physletb.2025.140056;arXiv:2507.02418. 共 3 位作者。覆盖 26 ≤ Z ≤ 84,用后向-前向蒙特卡罗传播不确定度,f_rms 偏差作 χ² 估计量,约束 f_rms ≤ 1.56。
- Vorabbi M., Gennari M., Finelli P., Giusti C., Navrátil P., "Toward a Microscopic Description of Nucleus-Nucleus Collisions", *Physical Review Letters* 135(17), 172501 (2025). DOI 10.1103/qxtf-5b4y;arXiv:2506.15406. 共 5 位作者。署名:英国萨里大学物理系;德国美因茨约翰内斯古腾堡大学核物理研究所与 PRISMA+ 卓越集群;意大利博洛尼亚大学物理与天文系与 INFN 博洛尼亚分部;INFN 帕维亚分部;加拿大 TRIUMF 与维多利亚大学。研究体系:α 打 ¹²C 与 ¹⁶O,100–300 MeV。
- Vorabbi M., Barbieri C., Somà V., Finelli P., Giusti C., "Microscopic optical potentials for medium-mass isotopes derived at the first order of Watson multiple-scattering theory", *Physical Review C* 109(3), 034613 (2024). DOI 10.1103/PhysRevC.109.034613. 共 5 位作者。密度用从头算自洽格林函数理论算,t 矩阵用与之一致的手征相互作用算。研究体系:钙与镍同位素链,实验室能量 65–200 MeV。
- Vorabbi M., Gennari M., Finelli P., Giusti C., Navrátil P., "A Phenomenological Extension for Microscopic Optical Potentials", arXiv:2607.05080 (2026)。〔预印本;未核实到正式期刊卷期。〕共 5 位作者。文中明确指出这类微观方法忽略了 Pauli 阻塞以外的介质修正,也未充分抓住高阶散射贡献,导致系统性低估吸收。
- Vorabbi M., Finelli P., Giusti C., "Determination of Nuclear Matter Radii by Means of Microscopic Optical Potentials: The Case of ⁷⁸Kr", *Few-Body Systems* 65(2), 47 (2024). DOI 10.1007/s00601-024-01919-z. 共 3 位作者。得出 ⁷⁸Kr 的物质半径 R(rms)_m = 4.12 fm,中子皮 R_np − 0.1 fm。
- Titus L. J., Nunes F. M., Potel G., "Explicit inclusion of nonlocality in (d,p) transfer reactions", *Physical Review C* 93(1), 014604 (2016). DOI 10.1103/PhysRevC.93.014604;arXiv:1601.01338. 共 3 位作者。署名:密歇根州立大学国家超导回旋加速器实验室 + 该校物理与天文系;劳伦斯利弗莫尔国家实验室(LLNL)。研究体系:¹⁶O、⁴⁰Ca、⁴⁸Ca、¹²⁶Sn、¹³²Sn、²⁰⁸Pb,束能 10、20、50 MeV。
- Stahulak L. M., Holt J. W., "Nonlocal nucleon-nucleus optical potentials from chiral effective field theory", *Physical Review C* (2026). DOI 10.1103/n6q9-1www;arXiv:2509.04665. 共 2 位作者。署名:德州农工大学回旋加速器研究所与物理与天文系。DOI 直查确认在线日期为 2026-08-27;截至本文写作时 CrossRef 尚未给出卷期与页码,故只写刊名与年份,不填卷期。研究体系:钙同位素链。版本提示:预印本 v1 为 2025-09-04(183 KB,12 图),v3 为 2026-09-07(3,804 KB,13 图),摘要措辞由「is a genuine time nonlocality」改为「is primarily a genuine time nonlocality」。引用时须认版本。
- Arellano H. F., Blanchon G., "Universal separable structure of the optical potential", *Physical Review C* 109(6), 064609 (2024). DOI 10.1103/PhysRevC.109.064609;arXiv:2405.13301. 共 2 位作者。署名:智利大学物理与数学科学学院物理系(FCFM,圣地亚哥,Av. Blanco Encalada 2008);法国 CEA DAM DIF(F-91297 Arpajon)+ 巴黎-萨克雷大学 CEA 的 LMCE 实验室(91680 Bruyères-le-Châtel)。扫描范围:束能 40–400 MeV,靶质量数 40 ≤ A ≤ 208。
- Jiang Y. C., Chen J., Kay B. P., Hoffman C. R., Tang T. L., Tolstukhin I. A., Xie M. R., Li J. G., Michel N., Avila M. L., Ayyad Y., Bazin D., Bennett S., Clark J. A., Freeman S. J., Jayatissa H., Li G., Liu W. P., Lou J. L., Munoz-Ramos A., Müller-Gatermann C., Nathan T., Santiago-Gonzalez D., Sharp D. K., Shen Y. P., Wuosmaa A. H., Yuan C. X., "Quenching of single-particle strength inferred from nucleon-removal transfer reactions on ¹⁵C", *Physics Letters B* 868, 139789 (2025). DOI 10.1016/j.physletb.2025.139789;arXiv:2506.03685. 共 27 位作者。数据采集于美国阿贡国家实验室的 ATLAS 装置(作者中含 B. P. Kay 等阿贡研究人员)与密歇根州立大学 FRIB(含 D. Bazin 等),合作方还涉及圣路易斯华盛顿大学、密歇根州立大学国家超导回旋加速器实验室(NSCL)等。作者名单与所属机构以 CrossRef 与 arXiv 反查结果为准。研究道:(p,d) 与 (d,t),结论 Rs = 0.72 ± 0.14。
- Nguyen N. B., Waldecker S. J., Nunes F. M., Charity R. J., Dickhoff W. H., "Transfer reactions and the dispersive optical model", *Physical Review C* 84(4), 044611 (2011). DOI 10.1103/PhysRevC.84.044611. 共 5 位作者。署名:密歇根州立大学国家超导回旋加速器实验室 + 该校物理与天文系;圣路易斯华盛顿大学物理系;圣路易斯华盛顿大学化学系(Charity)。
- Dickhoff W. H., Atkinson M. C., "Combining nuclear reactions and structure with the dispersive optical model", *Journal of Physics: Conference Series* 1643(1), 012082 (2020). DOI 10.1088/1742-6596/1643/1/012082. 共 2 位作者。第 27 届国际核物理大会(INPC2019)会议文集。署名:圣路易斯华盛顿大学物理系;TRIUMF。
- Yoshida K., Atkinson M. C., Ogata K., Dickhoff W. H., "First application of the dispersive optical model to (p,2p) analysis within the distorted wave impulse approximation framework", *Physical Review C* 105(1), 014622 (2022). DOI 10.1103/PhysRevC.105.014622;arXiv:2109.02211. 共 4 位作者。署名:日本原子能研究开发机构(JAEA)先端基础研究中心;TRIUMF;大阪大学核物理研究中心(RCNP)与大阪公立大学;圣路易斯华盛顿大学物理系。
- Atkinson M. C., Dickhoff W. H., "Learning from knockout reactions using a dispersive optical model", *Frontiers in Physics* 12, 1505982 (2025). DOI 10.3389/fphy.2024.1505982. 共 2 位作者。
- Hebborn C., Nunes F. M., "Systematic study of the propagation of uncertainties to transfer observables", *Frontiers in Physics* 13, 1525170 (2025). DOI 10.3389/fphy.2025.1525170;arXiv:2507.13063. 共 2 位作者。署名:巴黎-萨克雷大学 CNRS/IN2P3 的 IJCLab(法国奥赛);密歇根州立大学 FRIB;密歇根州立大学物理与天文系。用 KDUQ,研究 ⁴⁸Ca 上的 (d,p)。注意:该文的 68% 可信区间「半宽约 5–10%」与早期逐核拟合工作给出的「25–100%」量级相差一个数量级,两者量的是不同的东西(见正文第十三节)。
- Beyer K., Nunes F. M., "Towards a global, uncertainty-quantified optical potential for rare isotopes",2025 年 APS 核物理分会(DNP)秋季会议分会报告,报告日期 2025-10-18;会议地点美国伊利诺伊州芝加哥希尔顿酒店,会期 2025-10-17 至 20。〔会议报告摘要,未经同行评议。该报告的完整论文见条目 16。〕会议总计 887 个报告、160 个分会。
- 同上会议的工作坊(邀请制):《Workshop: Systematic and Model Uncertainties in Nuclear Experiment and Theory: New Tools, Techniques, and Opportunities for Collaboration》I 与 II,编号 1WJ、2WJ,时间 2025-10-17 上午;以及《Workshop: Aspects of Weak Binding Effects in Nuclear Structure & Reactions》I 与 II,编号 1WB、2WB,同为 2025-10-17 上午。〔会议日程信息,取自 APS 会议存档。〕
- Zhao X., Du W., Capote R., "Dispersive optical model of nucleon scattering on zirconium isotopes", *Physical Review C* 114(2), 024609 (2026). DOI 10.1103/2n39-9g7c. 共 3 位作者。CrossRef 给出的署名单位:渤海大学;渤海大学;Suncoast Data Evaluation。〔Capote 长期任职于国际原子能机构核数据科,本条按其在本论文中的署名填写,不作现职回溯。〕
- Jennings B. K., "Non-observability of Spectroscopic Factors", arXiv:1102.3721 (2011)。〔预印本;未核实到正式期刊卷期。署名:加拿大 TRIUMF。〕重要提示:另有一篇题名相近的综述《Spectroscopic Factors: Observability and Measurability》,载 *Physics of Particles and Nuclei* 53(2), 482–486 (2022),DOI 10.1134/s1063779622020204,其作者是 L. D. Blokhintsev,不是 Jennings,两者不可混引。
*说明:* 参考文献中的单位信息以各论文发表时的署名为准,不作现职回溯。作者人数取自 CrossRef 或 arXiv 作者数组长度,正式期刊版人数可能略多于预印本版。每条文献的核实状态标注如下:已 DOI 直查确认者见各条 DOI 字段;仅有 arXiv 者已在条目内标明「预印本」。
核实与纠错记录
| # | 常见说法 | 实际情况 |
|---|---|---|
| 1 | 光学势带来「约 30% 的系统不确定度」 | 30% 是从 (d,p)/(d,n) 提取谱因子的整套误差合计,不是光学势单独一项。拆账是:统计 3%、靶厚 5%、入射道光学势 15%、出射道光学势 15%、束缚态几何 20%、氘内核谱因子 15%。光学势(进+出)是其中最大的一项,但它不等于 30%。出自 *Sci. China Phys. Mech. Astron.* 62,32021 (2019)[6]。 |
| 2 | 贝叶斯方法一定能给出「更小、更可靠」的误差 | 贝叶斯的 95% 区间比频率方法宽,而经验覆盖率检验显示频率方法的区间被低估(把自己说得太确定),贝叶斯在较宽区间上覆盖率准确 [11]。误差棒是方法的函数,不只是数据的函数。 |
| 3 | 贝叶斯是纯粹的数据驱动、没有偏见 | 同一批 ⁶Li+²⁰⁹Bi 数据,三种先验(阈异常/反常阈异常/均匀)都能很好地重现弹性散射角分布,但导出的势的能量依赖走法截然不同 [7]。提出这一警告的就是使用贝叶斯方法的那批人自己。 |
| 4 | 「¹¹Be 的敲出实验理论误差和实验误差同量级」可以推广到所有核 | 不能。同一篇里 ¹¹Be(晕核)理论误差与实验同量级,¹²C(紧束缚核)的理论误差明显更大。给出的门槛是:至少 20%(晕核)、至少 40%(紧束缚核)[14]。 |
| 5 | 转移截面的不确定度就是 40–100% | 那是逐核逐能点单独拟合的情形 [8][11]。用已在上百条数据集上标定的全局势(KDUQ)算 ⁴⁸Ca(d,p),68% 可信区间半宽只有 5–10% [31]。两组数不矛盾,量的不是同一件事。引用任何一个都必须写清「哪个势、哪套似然、算不算相关性」。 |
| 6 | 谱因子误差里最大的旋钮是束缚态几何 | 在 ¹⁵C 的实测比较里,把束缚态半径 r₀ 从 1.25 换成 1.34 fm,结果只小约 3%;束缚态弥散度的影响远小于不同光学模型参数之间的差异 [26]。光学势才是最粗的那个旋钮。 |
| 7 | Woods-Saxon 的「核表面有吸收峰」是微观计算支持的 | 恰恰相反。基于动量空间介质内折叠模型的微观计算发现,虚部径向形状因子里一致地出现一个节点,显著压制表面吸收峰,与 Woods-Saxon 的假设明显相悖 [25]。 |
| 8 | 「非定域性」是一件东西 | 是两件。空间非定域与时间非定域。手征 EFT 微观光学势里,实部的能量依赖主要来自空间非定域,虚部的能量依赖是真正的时间非定域 [24]。把它们笼统记为「能量依赖」是把两笔账记成一笔。 |
| 9 | 只用局域光学势分析 (d,p),谱因子误差可以忽略 | 可能错最多 40%,而且这 40% 来自「把非定域的东西按局域处理」这个动作本身,与参数拟合得准不准无关 [23]。 |
| 10 | ¹⁵C 的淬灭因子已经被测定了 | 两个道给出的平均是 Rs = 0.72 ± 0.14[26],而重离子诱导敲出给的是 0.96 ± 0.04[26]。两者互相矛盾,尚未解决。0.14 的误差棒里有最大的一项是光学势。 |
| 11 | Perey-Buck 那篇的第一作者是「F. Ferey」 | 是 F. Perey。CrossRef 与 INSPIRE 的作者字段都有这个录入错误,照抄会把一个不存在的人写进参考文献 [5]。 |
| 12 | 引用 Jennings 的「谱因子不可观测」时,可以引 *Phys. Part. Nucl.* 53(2), 482 (2022) | 不行。那一篇的作者是 L. D. Blokhintsev,不是 Jennings。Jennings 的那篇没有核实到正式期刊卷期,只能按 arXiv:1102.3721 引用 [35]。 |
| 13 | ELM 已经有正式期刊论文了 | 截至 2026 年 9 月,逐条 CrossRef 查证未见正式卷期页码,只有 arXiv:2603.28599 [16];它在 DNP 2025 上的报告(2025-10-18)是会议摘要,未经同行评议[32]。两者都要标明性质。 |
| 14 | Stahulak & Holt 那篇说的是「虚部能量依赖是纯时间非定域」 | 预印本 v1 是这么写的,v3(2026-09-07)已改成「primarily a genuine time nonlocality」(主要是)[24],且文件从 183 KB 涨到 3,804 KB、图从 12 张变 13 张。引这句话必须认版本。 |
| 15 | Arellano & Blanchon 是美国组 | 是智利-法国两家。Arellano 在智利大学物理与数学科学学院物理系(圣地亚哥);Blanchon 在法国 CEA DAM DIF(Arpajon)并挂巴黎-萨克雷大学 CEA 的 LMCE 实验室(Bruyères-le-Châtel)[25]。 |
| 16 | ELM、KDUQ、CHUQ 都是「全局标定过的势」,所以预测应当一致 | 不一致。在 ¹⁰⁰Sn–¹⁷⁰Sn 的质子反应截面上三家差别明显,而 ELM 那篇明确指出差别不来自实验约束或物理模型形式,来自统计模型(似然函数的重新标度)[16]。 |
| 17 | 弹性散射数据能唯一确定光学势参数 | 不能。Igo 模糊现象下,不同参数组在灵敏半径处的实虚部深度交于近似同一点 [7];垒下能区尤其如此,三种先验都能拟合得上同一批角分布 [7]。 |
| 18 | 「纯微观」光学势可以直接替代唯象势 | 还不行。Vorabbi 等 2026 年的工作标题就叫《微观光学势的一个唯象修正》,文中直说这类方法忽略了 Pauli 阻塞以外的介质修正、也没抓住高阶散射贡献,导致系统性低估吸收,需要加能量依赖的唯象项打补丁 [21]。 |
| 19 | 极化数据比角分布更能约束光学势 | 未必。贝叶斯证据算出的贝叶斯因子不足以区分 ⁴⁸Ca 弹性散射的角分布与极化数据,两者约束力大致相当 [13]。 |
| 20 | 想减小光学势不确定度,把测量精度提高就行 | 收益「往往比天真预期的小得多」[12]。更有效的两招是换个角度测和加邻近能量的数据(后者可把不确定度改善最多两倍)[12]。 |
| 21 | 谱因子是核的固有属性,可以跨实验比较 | 按 Jennings 的论证,它不是可观测量:存在幺正变换使其在允许范围内任意变化;现象学上它「更多取决于你假设的势的能量依赖,而不是你测到的截面」[35]。它只在导出它的模型之内有意义。 |
| 22 | 光学势的同位旋依赖是个可以延后处理的技术细节 | 不是。WLH 势用于天体物理中子俘获时,反应率的不确定度在同位旋不对称度约 0.28 处出现明显分界,作者据此指出同位旋依赖是下一步的关键 [18]。这直接决定 r 过程网络里大片核的反应率能否算准。 |