⚛️ nuclear theory

Nuclear level density studied in odd-mass nuclei in the framework of the projected shell model

本文扩展了一种投影壳模型框架，用于计算奇质量核的核能级密度，揭示了单个核子的阻塞效应会削弱配对效应并抑制低能结构变化，从而导致与相邻偶偶核相比，其统计行为的出现更早且具有规则的高斯自旋分布。

原作者： Jiaqi Wang, Saumi Dutta, Cui-Juan Lv, Long-Jun Wang, Yang Sun

发布于 2026-08-13

📖 1 分钟阅读🧠 深度阅读

原作者： Jiaqi Wang, Saumi Dutta, Cui-Juan Lv, Long-Jun Wang, Yang Sun

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一下，原子核并非一个坚实的弹珠，而是一个繁忙的微观舞池。在这个微小的舞厅里，质子和中子就是舞者。在最稳定、最快乐的原子核中，这些舞者完美地配对，手牵手，同步旋转。这种“配对”是核舞的一个关键规则，它让原子核变得稳定且可预测。然而，有时一个原子核拥有奇数个舞者，这会让一个人孤零零地站在场边，无法找到舞伴。这名落单的舞者改变了整个派对的氛围。

科学家们非常感兴趣于计算这些舞者在不同能量层级下有多少种不同的排列方式。这种计数被称为“核能级密度”。这为什么重要呢？因为了解在不同能量水平下舞池有多拥挤，有助于我们理解恒星是如何燃烧的、核反应堆是如何工作的，以及我们如何处理核废料。这就像是了解房间里有多少人，从而预测噪音会有多大或热量会扩散得有多快。长期以来，科学家们对那些“完美配对”的舞池（偶-偶核）已经有了很好的掌握，但那些“奇数人数”的派对（奇质量核）却一直是个谜。他们想知道：那个孤独的舞者是否会让人群的行为变得不同，还是他们最终会融入人群之中？

本论文通过使用一种名为“投影壳模型”的高级计算机模拟，来观察两个特定奇质量核——镝-163（Dysprosium-163）和铽-163（Terbium-163）的舞池。研究人员发现，那名单一未配对核子的存在就像是一个游戏规则改变者。在完美配对的原子核中，舞池在低能级时非常有结构，有着特定的模式和规则。但在这些奇质量核中，这种结构坍塌得快得多。那名未配对的舞者本质上“阻碍”了配对过程，削弱了其他舞者之间手牵手的力量。结果，奇质量核比它们的配对邻居更早地进入了一种统计学上的混沌状态——即舞者们正进行着随机且可预测的运动。

具体而言，研究表明，虽然配对的原子核需要达到更高的能量水平（约 4.0 MeV）才会开始表现得像一群随机的人群，但奇质量核在较低的能量（约 2.25 到 2.75 MeV）时就开始呈现出这种行为。在这种“混沌”区域，奇质量核展现出一种完美的平衡：一半的舞者具有一种类型的自旋（宇称），另一半则具有另一种；它们的自旋遵循平滑的钟形曲线分布。作者认为，这种统计行为的提前到来是那名未配对核子直接导致的结果。他们还计算了一个特定的数值，称为离散度 $\sigma$ ，对于这些原子核而言约为 6.1，描述了自旋分布的离散程度。这一发现意义重大，因为它表明对于奇质量核，我们可以比以前认为的更早地在较低能量下信任简单的统计公式，这有助于改进用于核天体物理学和反应堆技术的模型。本文并不声称解决了所有的谜团，但它提供了一个清晰的模拟图景，展示了那一个额外的舞者是如何重塑整个核派对的。

技术摘要：投影壳模型框架下的奇质量原子核核能级密度

问题陈述
核能级密度（NLD）是豪瑟-费希尔巴赫（Hauser-Feshbach）理论中的关键输入，该理论描述了对天体物理学和核技术至关重要的辐射中子俘获反应。虽然实验技术（如奥斯陆方法）可以提供 NLD 数据，但这些数据通常给出的是自旋和宇称混合的总密度 $\rho(E_x)$ 。为了在反应计算中使用这些数据，必须假设变量分离： $\rho(E_x, I, \pi) = \rho(E_x)g(E_x, I)F(E_x, \pi)$ 。然而，确定自旋分布函数 $g$ 和宇称分布 $F$ 仍然是一个挑战。唯象模型（如费米气体模型）缺乏结构细节，而平均场微观模型往往违反角动量守恒。此外，与偶-偶系统相比，针对奇质量（奇-A）核的现有理论处理较为匮乏，且单个未配对核子对 NLD 中统计行为演化的具体影响尚未得到充分探索。

方法论
作者采用投影壳模型（PSM）来计算 $^{163}\text{Dy}$ （偶-奇）和 $^{163}\text{Tb}$ （奇-偶）这两个奇质量核的 NLD，并将其与它们的偶-偶邻居 $^{164}\text{Dy}$ 进行比较。

哈密顿量与基矢： 该模型使用包含可分离四极-四极相互作用、单极配对和四极配对相互作用的哈密顿量。配对强度经过调整（与偶-偶情况相比降低了约 5%），以解释奇-A 系统中减弱的配对关联。
构型空间： 不同于从 0-准粒子真空开始的偶-偶系统，奇-A 系统构建在 1-准粒子（1-qp）基矢之上。构型空间包括由变形 Nilsson 轨道构建的多准粒子激发（1-qp, 3-qp, 5-qp, 7-qp 等）。
角动量投影： 该方法从变形内禀态投影角动量以恢复旋转对称性，确保所有计算出的本征态都具有确定的自旋（ $I$ ）和宇称（ $\pi$ ）。
计算策略： 作者进行全对角化计算，直至激发能达到 4.5 MeV。为了将计算扩展到更高能量（超过中子分离能）并分析特定准粒子阶数的贡献，他们使用了“未混合”的 NLD。这些 NLD 通过对投影态进行求和而不进行最后的矩阵对角化步骤（忽略构型混合）来计算，这显著降低了计算成本，同时保留了定性的结构和量子数。

核心贡献与结果

结构主导变化的抑制：
研究发现，在偶-偶核中观察到的低能区显著的结构依赖性变化在奇-A 核中得到了极大抑制。在奇-A 核中，由于存在一个未配对的核子阻挡了配对形成，导致配对能隙（ $\Delta$ ）变弱。因此，与偶-偶邻居相比，奇-A 系统在显著较低的激发能下就向统计行为过渡。
机制分类与配对破坏：
奇-A 核的 NLD 演化被分为三个阶段：
- 未配对机制（ $^{163}\text{Dy}$ 为 0 – ~1.25 MeV， $^{163}\text{Tb}$ 为 0 – ~1 MeV）： 由未配对的 1-qp 态主导。由于靠近费米能级的单粒子态的存在，其能级密度本质上比相邻的偶-偶核高出几个数量级（几倍到十倍）。
- 配对破坏机制： 在奇-A 核中，配对破坏仅发生在未被阻挡的核子类型内（例如，奇-中子核中的质子对）。这所需的能量（ $\sim 2\Delta \approx 1.5$ MeV）比偶-偶系统中中子和质子对同时破坏所需的能量（ $\sim 2(\Delta_n + \Delta_p) \approx 3$ MeV）要低。因此，奇-A 核中的配对破坏阶段及随后的统计行为起始点出现在更低的能量处。
- 多准粒子（混沌）机制： 在超过 $\sim 2.25$ MeV ( $^{163}\text{Tb}$ ) 和 $\sim 2.75$ MeV ( $^{163}\text{Dy}$ ) 之后，系统由高阶多准粒子态（5-qp, 7-qp）主导。在此阶段，NLD 表现出统计特性。
宇称分布：
一个关键发现是奇-A 核中宇称分布达到相等的速度很快。虽然偶-偶系统在高能下仍表现出偶宇称态占主导的现象，但在奇-A 核中，奇宇称与偶宇称密度的比值在低至 2 MeV 时就趋于单位值（相等分配）。这是因为未配对核子允许比偶-偶系统中需要通过配对破坏才能进入奇宇称态的情况更早地实现宇称态的平衡分布。
自旋分布与埃里克森（Ericson）公式：
PSM 结果提供了具有确定自旋的本征态，从而允许直接进行自旋分布分析。
- 在低能未配对机制中，自旋分布是不规则且依赖于宇称的。
- 在多准粒子机制中（ $> 2.75$ MeV），两种宇称的自旋分布均呈现规则的高斯形状。
- 作者成功地将这些分布拟合到埃里克森自旋分布公式，并提取了能量依赖的色散参数 $\sigma$ 。对于 $^{163}\text{Dy}$ 中的多准粒子机制，两种宇称的 $\sigma \approx 6.1$ 。
- 研究表明，高自旋态比低自旋态更早地接近统计行为。
验证：
计算得到的 $^{163}\text{Dy}$ 和 $^{163}\text{Tb}$ 的离散能级和 NLD 与现有的能谱数据（NNDC）以及可获得的奥斯陆方法测量结果表现出良好的一致性。“未混合”的 NLD 计算成功重现了中子分离能处的实验数据点（ $^{163}\text{Dy}$ 为 6.271 MeV），验证了该方法用于向更高能量外推 NLD 的有效性。

意义
本文声称，投影壳模型提供了一个鲁棒的微观框架，用于计算具有显式自旋和宇程依赖性的 NLD，而这在唯象或平均场方法中往往是缺失的。其主要意义在于证明，由于单核子阻挡及随后的单类配对破坏机制，奇-A 核比其偶-偶对应物更早地表现出统计行为（即宇称相等分配和高斯自旋分布）。这一发现表明，对于奇-A 核，假设宇称相等分布并使用埃里克森自旋分布公式在相对较低的能量（ $\sim 2$ MeV）开始就是有效的，这对于提高豪瑟-费希尔巴赫形式中反应截面预测的准确性至关重要。作者还指出，他们的方法提供了一套完整的多体波函数，为未来研究核热力学性质和相变提供了途径。

类似论文