A Discrete Occupation-Flow Diagnostic for Orbital-Filling Crossovers in sd-Shell Nuclei: A Multi-Interaction Benchmark in Ne, Mg, and Si

技术摘要：用于 sd 壳核系中轨道填充交叉的离散占据-流向诊断法

问题陈述
在核壳模型中，轨道占据数传统上被报告为特定同位素的静态快照（如 $n(1d_{5/2})$, $n(2s_{1/2})$, $n(1d_{3/2})$ 表格）。虽然这种表示法是准确的，但它未能明确表征随着中子增加，原子核如何重新组织其价层构型。现有的壳演化诊断工具——如激发能 $E(2_1^+)$、跃迁强度 $B(E2)$ 和有效单粒子能 (ESPE)——要么对全多体波函数进行积分，要么仅提供单一 $N$ 的快照。目前缺乏一种诊断工具，能够在不要求标准壳模型输出之外的额外计算机制的情况下，显式地追踪沿同位素链的轨道填充的定向、步进式重组过程。

方法论
作者引入了一种被称为轨道占据-流向向量 $\vec{F}(N)$ 的离散动力学量，它定义为具有中子数 $N$ 和 $N+2$ 的同位素之间三个 $sd$壳轨道（$1d_{5/2}$, $2s_{1/2}$, $1d_{3/2}$）的中子占据数的有限差分：
$$\vec{F}(N) = (\Delta n_{d5/2}, \Delta n_{s1/2}, \Delta n_{d3/2})$$
其中 $\Delta n_j = n_j(N+2) - n_j(N)$。

从该向量中可以导出两个标量指标：

1. 曲率 ($\kappa(N)$)： 相邻流向向量 $\vec{F}(N)$ 与 $\vec{F}(N+2)$ 之间的夹角（以度为单位）。它衡量了主导填充轨道的改变的剧烈程度。
2. 散度 ($\nabla \cdot \vec{F}(N)$)： 步进之间流向向量幅度的变化，指示了重组速率的加速或减速。

本研究利用 NuShellX 代码对 $sd$价层空间（以$^{16}$O 为核）内的氖 (Ne)、镁 (Mg) 和硅 (Si) 的偶偶同位素链进行壳模型计算，直至$N=20$ 壳闭合。研究采用了三种独立的 USD 系列相互作用（USDA、USDB 和 USDC）来测试该诊断法的鲁棒性。作者还计算了直接源自相互作用哈密顿量的基于单极修正的有效单粒子能 (ESPE)，以将流向行为锚定在底层物理机制上。

关键结果

• 曲率的鲁棒性： 曲率 $\kappa(N)$ 在所有三种相互作用（USDA、USDB、USDC）中均表现出高度可复现的模式，其变异系数 (CV) 通常低于 9%，在主峰处低于 2%。这证实了 $sd$ 壳结构中轨道填充的顺序交叉是一个鲁棒的特征，而非特定相互作用参数化的产物。
• 识别交叉点： $\kappa(N)$ 识别出了主导填充轨道发生切换的具体中子数。在 Ne 和 Si 中，在 $N=14$ 处出现一个主峰（标志着 $1d_{5/2} \to 2s_{1/2}$ 的交叉），而在 Mg 中，峰值出现在 $N=12$。观察到次级交叉（$2s_{1/2} \to 1d_{3/2}$）发生在 $N=16-18$ 附近。
• 与集体性的独立性： 与最初的假设相反，统计检验（皮尔逊和斯皮尔曼相关分析）显示，$\kappa(N)$ 与 $B(E2; 0_1^+ \to 2_1^+)$ 之间不存在系统性相关。作者得出结论，$\kappa$ 追踪的是接收中子的轨道的身份（单粒子流向），而 $B(E2)$ 追踪的是由此产生的四极矩形变的相干性（集体相干性）。这两者是不同的物理层面。
• 相互作用敏感性： 基于三种相互作用下占据数标准差之和的相互作用敏感性图谱显示，在所有三个链中，$N=16$ 都是最大相互作用依赖点。这对应于 $2s_{1/2} \to 1d_{3/2}$ 交叉正在进行的区域，以及相互作用的相对单粒子能差异最显著的区域。
• 验证与局限性： 与实验数据的对比表明，在远离 $N=20$ 的区域，$sd$壳计算在$E(2_1^+)$和$B(E2)$方面是定量可靠的。然而，正如预期，该模型在$^{30}$Ne 和$^{32}$Mg（“反转岛”）处无法定性地描述其结构，因为它无法重现由$pf$壳入侵构型引起的低$E(2_1^+)$和高$B(E2)$。该模型之所以失效，是因为$sd$ 壳空间排除了这些入侵态。
• 哈密顿量联系： 基于单极子的 ESPE 计算证实，$\kappa(N)$ 的峰值对应于由中子-中子单极相互作用驱动的 $1d_{3/2} - 2s_{1/2}$ 能隙的缩小，从而将流向诊断法建立在哈密顿量之上。

意义与主张
本文声称，占据-流向向量 $\vec{F}(N)$ 及其曲率 $\kappa(N)$ 提供了一种新型的、低成本的壳演化诊断工具。其主要意义在于：

1. 动态洞察： 它将静态的占据数表格转化为定向的、步进式的轨道重组图谱，揭示了静态数据所掩盖的填充“轨迹”。
2. 物理区分： 它确立了轨道填充交叉（由 $\kappa$ 追踪）与集体结构变化（由 $B(E2)$ 追踪）并不一定耦合，挑战了“填充发生剧烈变化必然同时导致集体性改变”的假设。
3. 诊断效用： 它识别出 $N=16$ 是 $sd$壳中模型敏感性的关键区域，在此区域，即使衍生观测值（如$B(E2)$）保持稳定，占据成分的预测也会高度依赖于相互作用的选择。
4. 通用性： 该框架被呈现为一种通用的方法论，适用于任何固定壳模型空间内的同位素链，且无需超出标准对角化输出之外的额外计算成本。

作者明确指出，$\kappa$ 不是占据数的参数化，而是一个捕捉轨迹信息的独立量。他们还注意到，虽然 $\kappa$ 在 $sd$壳内是鲁棒的，但其在反转岛上的应用需要跨壳模型空间（$sd$-$pf$），这超出了目前的基准测试范围。未来的工作拟将此扩展到静态四极矩$Q(2_1^+)$，并将质子-中子单极项纳入 ESPE 计算，以进行更具预测性的测试。