从头计算的 E1 强度如何改变镍同位素链上的中子俘获率？

E1 强度补上了仅用 M1 时缺失的部分

富中子核的中子俘获率是 r-process 核合成模拟的输入，而其中许多核无法直接测量。标准工具是 Hauser-Feshbach 统计模型，本文用反应程序 TALYS 实现，它需要光子强度函数和核能级密度作为输入。本文只隔离其中一个输入：中子阈值附近及以上的电偶极（E1）强度，作者用第一性原理计算它，而不是取自唯象模型。

承载结果的比较被刻意限定得很窄。作者先只用唯象 M1 强度计算反应率，再在同一 M1 基准上叠加他们从头计算的 E1 强度。因此报告的增强是加入 E1 强度带来的效应，而不是对全部光子强度输入的完整替换。这一基准选择很重要：下文百分比描述的是相对仅含 M1 的反应率的增长，而不是相对完整反应率计算的增长。

LIT-CC 借助轴对称参考态进入开壳镍核

E1 强度由洛伦兹积分变换耦合簇方法（LIT-CC）得到。耦合簇理论在 Hartree-Fock 参考态之上构建关联，洛伦兹积分变换则避免显式构造构成粒子阈以上连续区响应的散射态。响应改为通过数值反演一个变换来恢复，而所报不确定度的很大一部分正来自这里。

闭壳核允许使用球对称参考态，但处于壳层中部的镍同位素不行。作者把方法推广到轴对称参考态，此时好量子数是角动量在对称轴上的投影，而不是总角动量。E1 算符随后按分量分别处理，响应由三个球分量取平均来近似。零投影分量与正负一投影分量之间的差异，被用作对称性破缺效应的估计，因为在旋转不变的计算中这一差异应当消失。

计算采用手征有效场论给出的两核子和三核子相互作用，主要是 Delta NNLO GO(394) 相互作用；对 68Ni 和 78Ni 还使用 N3LO Texas 相互作用，以评估对哈密顿量的依赖。覆盖范围是 58Ni 以及偶偶同位素 68Ni 至 80Ni。

两个数据锚点先检验方法再外推

在沿链应用方法之前，作者先用测量量检验它。计算的光吸收截面与 58Ni 的实测截面比较，计算的 E1 强度与 68Ni 的实测 E1 强度比较。正因为有这两个锚点，后续富中子预言才被呈现为经过检验的框架的延伸，而不是无约束的模型输出。

58Ni 的情形还说明了一个会改变其反应率行为的细节。对复合核 58Ni，第一个粒子阈值是 8.17 MeV 的质子分离能，比 12.22 MeV 的中子分离能低数 MeV。因此 E1 强度由跃迁到中子分离能加入射中子能再减质子分离能以下末态的跃迁来探测，这也是 57Ni(n,gamma)58Ni 反应率在所考察的整个温度区间都对 E1 强度敏感的原因。

低能 E1 强度随中子数增加而增强

沿链观察，作者在 68Ni 至 80Ni 的所有偶偶同位素中都发现了中子阈以上的低能 E1 强度。其幅度随 neutron number 增加，这一增长正是反应率增强也随 neutron number 增加的结构性原因。阈值附近的强度与中子俘获最相关，因为它恰好位于光子强度函数为俘获计算提供输入的位置。

这一低能强度的幅度目前受积分变换反演的限制。作者刻意保留了一族较宽的可接受反演解，改变非线性参数和基函数数目，并把其包络作为不确定度带报告。他们指出，更严格的接受判据或不同的反演策略未来可以减小这一不确定度。对 78Ni，模型空间收敛检验显示，在变换变量低于约 30 MeV 的区间，两个模型空间大小的变换相差不到 1%，而这正是偶极响应相关的区域。

生成 68-78Ni 的反应率在 3 至 10 GK 上升

把从头计算的 E1 强度输入 TALYS 后，生成偶偶同位素 68Ni 至 78Ni 的中子俘获率主要在 3 至 10 GK 之间被提高。效应随温度和 neutron number 增加。在不确定度带的上边缘，5 GK 时相对仅含 M1 基准的增强从 67Ni(n,gamma)68Ni 的 16% 升至 77Ni(n,gamma)78Ni 的 24%。在 10 GK，相应的总反应率达到各自仅含 M1 基准的 1.8 倍和 3 倍。

这些数字应被读作反演不确定度范围内的上边缘估计，而不是中心值。它们描述的也是相对仅含 M1 计算的比值，因此本身并不能说明完整反应率与测量值或数据库值相差多远。作者指出，他们的反应率与 JINA REACLIB 和 BRUSLIB 的表列反应率在一个数量级内一致，这一差距与该区域统计模型反应率的典型不确定度相当；但他们把这一比较称为指示性的，因为三套反应率都没有完整的 uncertainity 量化。

80Ni 把增强提前到约 0.5 GK

最明显偏离整链规律的情形出现在 80Ni。那里紧邻阈值之上存在强 pygmy 偶极强度，79Ni(n,gamma)80Ni 的反应率在约 0.5 GK 就已偏离仅含 M1 的结果。这一温度区间远低于支配链中较轻同位素的 3 至 10 GK 窗口，它直接来自低能强度相对于阈值所处的位置。

作者还把他们的趋势与早先一项基于能量密度泛函的研究比较。那项研究用含与不含 pygmy 偶极共振的反应率之比来量化该共振，同时两者都保留巨偶极共振。作者提醒，他们自己的比较方式不同，因为他们比较的是阈以上全部 E1 贡献与不含 E1 强度的计算，所以两组比值不能直接比较。他们报告的共同点是 E1 效应随 neutron number 增加的趋势。

从头计算输入止于粒子阈值

计算的 E1 强度覆盖第一个粒子发射阈以上的响应。在阈值以下，作者对闭壳的 68Ni 和 78Ni 各找到一个约 1 MeV 的低能 1 负态，它只携带极小部分 E1 强度；作者判断这在反演不确定度的尺度上可忽略，对本文讨论的俘获率也不重要。他们没有在开壳同位素中寻找束缚态，因为在对称性破缺框架中准确描述这些态需要角动量投影，而这留待未来工作。

这一边界界定了本文交付了什么、没有交付什么。它给出了一条镍同位素链在阈值以上的从头计算 E1 贡献，而不是完整的光子强度函数。在这些反应率计算中，M1 部分仍是唯象的；作者把将从头计算描述扩展到第一个粒子阈值以下的 E1 和 M1 强度（包括 M1 驱动的低能上弯）列为下一步，这需要角动量投影。

80Ni 的低温增强是最值得关注的情形

本文最重要的区分不是同位素之间的一般差异，而是 80Ni 与链中其余核的差异。对 68Ni 至 78Ni，E1 贡献在 3 至 10 GK 窗口起作用并随 neutron number 增大，在不确定度带上边缘、10 GK 时最多达到仅含 M1 基准的 3 倍。对 80Ni，紧邻阈值的 pygmy 强度把起始温度下移到约 0.5 GK，这是另一个天体物理区间，而不是同一效应的更大版本。

作者据此提出后续问题：本文发现的 E1 效应，尤其是 80Ni 的低温增强，是否会在质量数 100 以下的丰度上留下印记。这正是这份从头计算输入所支持的具体检验，它取决于阈值附近的低能强度、目前约束它的反演不确定度，以及每个同位素反应率实际响应的温度。

本文引用的文献