Origin of the electric hexadecapole isomer in $^{93}$Mo

基本信息

• 论文ID: 2510.22352
• 标题: Origin of the electric hexadecapole isomer in $^{93}$Mo
• 作者: B. Maheshwari (GANIL), P. Van Isacker (GANIL), P. M. Walker (University of Surrey)
• 分类: nucl-th (核理论), nucl-ex (核实验)
• 提交时间: 2025年10月25日
• 论文链接: https://arxiv.org/abs/2510.22352

摘要

本文通过壳模型分析研究了$^{93}$Mo核的异常$21/2^+$同质异能态（isomer）的微观起源。该态是核电子俘获激发（NEEC）过程的重要候选者。其独特之处在于衰变由慢速电六十四极（E4）跃迁主导，而通常更快的电四极（E2）衰变路径在能量上被禁戒。研究通过详细分析初态和末态波函数结构及E4跃迁矩阵元，揭示了这一现象的微观机制，并与其粒子-空穴共轭核$^{99}$Cd（不存在E4跃迁）进行了对比分析。

研究背景与动机

研究问题

1. 核心物理问题：$^{93}$Mo中的$21/2^+$同质异能态为何通过罕见的E4跃迁衰变，而非常见的E2跃迁？这一现象的微观核结构起源是什么？
2. 实际应用背景：该同质异能态是唯一报道过核电子俘获激发（NEEC）现象的核态。NEEC过程涉及用5 keV触发从长寿命$21/2^+$态（半衰期6.85小时）到短寿命$17/2^+$态（半衰期3.53纳秒）的跃迁，可能实现至少50倍的能量增益。

问题重要性

1. 基础物理意义：高多极跃迁对波函数的径向和角向重叠极为敏感，是探测核结构的优秀探针
2. 能源应用前景：理解NEEC机制对开发潜在的核能存储和释放技术至关重要
3. 同质异能态物理：长寿命激发态的形成机制是核物理的基本问题

现有研究局限

1. 大规模壳模型计算虽能重现实验特征，但缺乏透明的定量分析
2. 对E4跃迁本身的微观机制（波函数成分、耦合方式）缺乏明确解释
3. 前人对关键的B(E2; $17/2^+ → 21/2^+$)值估计可能存在偏差

研究动机

通过简化的壳模型空间（$\nu1d_{5/2} \otimes \pi0g^2_{9/2}$）进行精细分析，明确量化中子-质子相互作用在能级反转中的作用，并通过与$^{99}$Cd的对比验证机制的普适性。

核心贡献

1. 揭示能级反转机制：定量证明了$21/2^+$和$17/2^+$态的能级反转源于最大自旋对齐组态中强吸引的中子-质子（νπ）相互作用
2. E4跃迁微观分解：首次明确分解了E4跃迁矩阵元的各个贡献成分，发现主要由质子两体矩阵元$\pi0g^2_{9/2}(J_\pi=8) → \pi0g^2_{9/2}(J_\pi=4)$贡献，中子和质子部分呈相长干涉
3. 修正NEEC关键参数：将$17/2^+ → 21/2^+$的B(E2)理论值相比前人估计降低约40%，对NEEC过程评估具有重要意义
4. 粒子-空穴对比验证：通过与$^{99}$Cd的系统对比，证明了Pandya变换下粒子-空穴组态中νπ相互作用变为排斥，导致无能级反转和E4同质异能态，从而验证了提出机制的正确性
5. 波函数精细分析：准确计算了各能级的波函数组成，最低$13/2^+$态约80%为$J_\pi=4$质子耦合，$17/2^+$态主要为$J_\pi=6$耦合

方法详解

任务定义

输入：核壳模型哈密顿量参数（单粒子能量、两体相互作用矩阵元）
输出：能级结构、波函数组成、电磁跃迁概率B(EL)
约束：限定在$\nu1d_{5/2} \otimes \pi0g^2_{9/2}$组态空间

理论框架

1. 壳模型哈密顿量

采用标准形式： $\hat{H} = \sum_i \epsilon_i \hat{n}_i + \frac{1}{4}\sum_{ijkl} V_{ijkl} a^\dagger_i a^\dagger_j a_l a_k + \cdots$

参数确定策略：

• 同类粒子相互作用$V_{\pi\pi}$：从$^{92}$Mo（两价质子在$0g_{9/2}$轨道）的实验能谱经验确定
  • $V_{\pi\pi}(J=0)=0.000$, $V_{\pi\pi}(J=2)=1.510$, $V_{\pi\pi}(J=4)=2.283$, $V_{\pi\pi}(J=6)=2.612$, $V_{\pi\pi}(J=8)=2.761$ MeV
• 中子-质子相互作用$V_{\nu\pi}$：从$^{92}$Nb（$\nu1d_{5/2}$与$\pi0g_{9/2}$耦合）的能谱固定
  • $V_{\nu\pi}(J=2)=0.135$, $V_{\nu\pi}(J=3)=0.286$, $V_{\nu\pi}(J=4)=0.480$, $V_{\nu\pi}(J=5)=0.357$, $V_{\nu\pi}(J=6)=0.501$, $V_{\nu\pi}(J=7)=0.000$ MeV
  • 关键特征：最大对齐态$J=7$具有最强吸引作用（设为零点）

2. 电磁跃迁算符

约化跃迁概率： $B(EL;\alpha_i J_i \to \alpha_f J_f) = \frac{|\langle \alpha_f J_f ||\hat{O}(EL)||\alpha_i J_i\rangle|^2}{2J_i+1}$

电跃迁算符： $\hat{O}(ELM) = \sum_k e_k r^L_k Y^M_L(\theta_k, \phi_k)$

多体态展开： $|\alpha_i J_i M_i\rangle = \sum_k a_k(\alpha_i J_i)|kJ_i M_i\rangle$

导致广义跃迁概率： $B(EL;\alpha_i J_i \to \alpha_f J_f) = \frac{1}{2J_i+1}\left|\sum_{kl} a_k(\alpha_i J_i)b_l(\alpha_f J_f)\langle lJ_f||\hat{O}(EL)||kJ_i\rangle\right|^2$

该表达式可解释中子和质子贡献的干涉效应。

3. 单粒子约化矩阵元

对E4跃迁： $\langle \nu1d_{5/2}||\hat{O}_\nu(E4)||\nu1d_{5/2}\rangle = \frac{513}{4}\sqrt{\frac{1}{7\pi}}e_\nu b^4$$\langle \pi0g_{9/2}||\hat{O}_\pi(E4)||\pi0g_{9/2}\rangle = \frac{9}{4}\sqrt{\frac{715}{\pi}}e_\pi b^4$

其中$b$为谐振子长度参数：$b^2 = 41.46/(45A^{-1/3} - 25A^{-2/3})$ fm²

技术创新点

1. 简化模型空间的有效性：虽然只考虑一个中子轨道和一个质子轨道，但通过精心选择两体相互作用参数，能准确重现实验能谱和跃迁强度
2. 系统的参数确定：不是拟合$^{93}$Mo本身，而是从邻近核的实验数据确定参数，保证了预测的可靠性
3. 粒子-空穴对称性分析：利用Pandya变换建立$^{93}$Mo和$^{99}$Cd的理论联系，通过对比验证机制
4. 跃迁矩阵元的完整分解：明确给出各波函数成分对跃迁强度的贡献及相位关系

实验设置

计算工具

• 使用Mathematica代码shell.m在J-scheme下进行计算
• 采用分数母系数（coefficients of fractional parentage）方法
• 需要数值对角化求解本征值问题

有效电荷选择

E2跃迁：

• 质子有效电荷：$e_\pi = 1.32$（标准值）
• 中子有效电荷：采用三个值
  • $e_\nu = 1.18$（从$^{91}$Zr的$Q(5/2^+)$四极矩确定）
  • $e_\nu = 1.73$（从$^{92}$Zr的$B(E2; 2^+_1 \to 0^+_1)$确定）
  • $e_\nu = 1.48$（平均值）

E4跃迁：

• 质子有效电荷：$e_\pi = 1.20$（从$^{92}$Mo的测量B(E4)值确定）
• 中子有效电荷：沿用E2算符的值

对比方法

主要与Hasegawa等人（2011）的大规模壳模型计算结果对比

实验结果

主要结果

1. 能级结构（表1）

关键发现：成功重现$21/2^+$态（2.437 MeV）低于$17/2^+$态（2.454 MeV）的能级反转，误差仅13 keV，与实验5 keV的能量差相当接近。

2. 电磁跃迁概率（表2，单位：W.u.）

重要结果：

• E4跃迁强度与实验符合优异（1.37 W.u. vs 1.449 W.u.，取平均有效电荷）
• $17/2^+ \to 21/2^+$的B(E2)值比Hasegawa的估计低约40%（2.21 vs 3.5 W.u.）
• $17/2^+ \to 13/2^+$的E2跃迁作为波函数精度的关键检验，理论值3.98 W.u.与实验4.48 W.u.符合良好

波函数分析（图2）

$13/2^+_1$态组成：

• $J_\pi=4$：~80%（主导成分）
• $J_\pi=6$：~15%
• $J_\pi=8$：~5%

$17/2^+_1$态组成：

• $J_\pi=6$：~85%（主导成分，较低对齐）
• $J_\pi=4$和$J_\pi=8$：各约7-8%

$21/2^+$态：

• 100%纯$J_\pi=8$组态（最大对齐）

E4跃迁矩阵元分解（图5）

质子贡献（单位：$b^4$）：

• $\pi0g^2_{9/2}(J_\pi=8) \to \pi0g^2_{9/2}(J_\pi=4)$：最大贡献（~3.0）
• $\pi0g^2_{9/2}(J_\pi=8) \to \pi0g^2_{9/2}(J_\pi=6)$：中等贡献（~1.5）
• $\pi0g^2_{9/2}(J_\pi=8) \to \pi0g^2_{9/2}(J_\pi=8)$：小贡献且相位相反（~-0.5）

中子贡献：$1d_{5/2} \to 1d_{5/2}$，矩阵元为23.32 $b^4$

总跃迁概率： $B(E4; 21/2^+ \to 13/2^+) = (0.83e_\nu + 3.53e_\pi)^2 b^8$

中子和质子贡献符号相同，呈相长干涉，增强了跃迁强度。

与$^{99}$Cd的对比（图1）

$^{93}$Mo：

• $21/2^+$态（2.424 MeV）< $17/2^+$态（2.429 MeV）→ E2禁戒 → E4同质异能态
• νπ相互作用：粒子-粒子，最大对齐时吸引

$^{99}$Cd：

• $21/2^+$态远高于$17/2^+$态 → E2允许 → 无E4同质异能态
• νπ相互作用：粒子-空穴（相对$^{100}$Sn核），最大对齐时排斥（Pandya变换）

这一对比完美验证了νπ相互作用性质决定能级序的机制。

物理图像分析

角动量特性（图3）

展示了4阶球谐函数$Y^M_4(\theta,\phi)$（M=-4至+4）的三维形态：

• 复杂的节面结构和角对称性
• 核子必须经历显著的轨道运动重排才能匹配这些角分布
• 红/蓝色表示正/负相位

径向重叠（图4）

比较了E2和E4跃迁的归一化径向被积函数：

• E4的$r^4$依赖使跃迁对核表面波函数行为高度敏感
• E4贡献峰值在更大半径（~4-5 fm），偏离质子在核内部的最可几位置（~2-3 fm）
• $0g_{9/2}$轨道比$1d_{5/2}$更紧凑，但都属于同一谐振子壳层

径向积分值：

• $\langle \nu1d_{5/2}|r^4|\nu1d_{5/2}\rangle = 42.75 b^4$
• $\langle \pi0g_{9/2}|r^4|\pi0g_{9/2}\rangle = 35.75 b^4$

E4跃迁的实现需要复杂角重排和径向重叠同时满足。

相关工作

核同质异能态研究

1. 经典理论：Auerbach和Talmi（1964）最早提出$^{93}$Mo中νπ相互作用导致能级反转的定性解释
2. 大规模计算：Hasegawa等（2011）用更大组态空间的壳模型计算重现了能谱特征
3. 本文改进：在简化空间内提供透明的定量分析和修正的跃迁强度

NEEC研究

1. 实验报道：Chiara等（2018, Nature）首次报道$^{93}$Mo的NEEC现象
2. 争议验证：Guo等（2022, PRL）对该结果提出质疑，仍在积极研究中
3. 应用前景：Carroll和Chiara（2024）讨论了NEEC的能源应用潜力

高多极跃迁

• 高多极跃迁因Weisskopf估计的强烈抑制而罕见
• E4跃迁对核结构极为敏感，是检验理论模型的严格测试

粒子-空穴对称性

• Pandya变换（1956）建立了粒子-粒子和粒子-空穴相互作用的关系
• 本文是该对称性在实际核结构问题中的成功应用案例

结论与讨论

主要结论

1. 能级反转机制明确：$^{93}$Mo中$21/2^+$和$17/2^+$态的反转完全源于$\nu1d_{5/2} \otimes \pi0g^2_{9/2}$组态中最大自旋对齐（J=7）时强吸引的νπ相互作用
2. E4同质异能态成因：能级反转使快速E2衰变路径在能量上被禁戒，迫使$21/2^+$态通过慢速E4跃迁到$13/2^+$态衰变
3. 跃迁机制定量化：E4跃迁由质子两体矩阵元主导（$J_\pi=8 \to 4$贡献最大），中子和质子部分相长干涉，总强度为$(0.83e_\nu + 3.53e_\pi)^2 b^8$
4. NEEC参数修正：关键的$B(E2; 17/2^+ \to 21/2^+)$值比前人估计降低~40%，对NEEC过程评估有重要影响
5. 机制普适性验证：$^{99}$Cd中粒子-空穴组态使νπ相互作用变排斥，无能级反转和E4同质异能态，完美验证了提出的机制

局限性

1. 模型空间限制：仅考虑$\nu1d_{5/2} \otimes \pi0g^2_{9/2}$组态，忽略了其他轨道的可能贡献（如$\pi0g_{7/2}$、$\nu0g_{7/2}$等）
2. 有效电荷不确定性：中子有效电荷在1.18-1.78范围内变化，虽对E4影响较小，但仍带来系统误差
3. 能量差精度：计算的$21/2^+-17/2^+$能量差为13 keV，实验仅5 keV，相对误差较大（虽绝对误差小）
4. 三体及更高相互作用：壳模型截断在两体层面，可能遗漏多体关联效应
5. 芯激发忽略：假设惰性芯$^{88}$Sr，实际可能有芯极化和芯激发贡献

未来方向

1. M4同质异能态：提出$^{84}$Rb中的M4同质异能态需要类似研究（文中提及正在进行）
2. 更多候选核：在核素图其他区域寻找类似的高多极同质异能态
3. NEEC实验验证：为正在进行的NEEC实验争议提供更精确的核结构输入
4. 同质异能态耗竭：研究诱导同质异能态耗竭的可能途径
5. 扩展模型空间：在更大组态空间验证简化模型的结论

深度评价

优点

1. 物理图像清晰：通过简化模型空间和系统的参数确定，将复杂的核结构问题归结为清晰的物理机制——νπ相互作用在最大对齐组态中的强吸引
2. 定量分析透彻：
   • 明确给出各波函数成分的百分比
   • 完整分解E4跃迁矩阵元的各项贡献
   • 量化中子和质子的干涉效应
3. 实验符合优异：
   • 能级结构误差<100 keV
   • E4跃迁强度与实验符合在10%以内
   • 关键的E2跃迁作为检验也符合良好
4. 对比分析有力：$^{93}$Mo与$^{99}$Cd的系统对比通过粒子-空穴对称性验证了机制，增强了结论的可信度
5. 可视化出色：
   • 图3展示球谐函数的三维形态，直观理解角动量特性
   • 图4的径向分布对比揭示E4跃迁的表面敏感性
   • 图5的矩阵元分解清晰展示各成分贡献
6. 实用价值高：为NEEC研究提供了修正的关键参数，对能源应用研究有直接指导意义

不足

1. 模型空间验证不足：虽然简化模型给出了好的结果，但未系统检验包含更多轨道时结论的稳定性
2. 不确定度分析欠缺：
   • 未给出计算结果的理论误差估计
   • 有效电荷的选择对结果的影响虽有讨论但不够系统
   • 参数拟合的误差传播未量化
3. $^{99}$Cd分析简略：虽提到可重现其能级结构，但未给出详细数据和图表，对比分析可更充分
4. 芯效应讨论缺失：完全忽略芯极化和芯激发，未评估其可能影响
5. E2禁戒程度：虽能级反转导致E2"能量禁戒"，但实际能差仅5 keV，该跃迁虽慢但并非完全禁止，文中对这一细节的讨论不够
6. 实验预测有限：除了修正B(E2)值外，未提出可供实验检验的新预测

影响力

1. 学术贡献：
   • 为核同质异能态物理提供了典型案例的深入分析
   • 展示了简化壳模型在理解复杂核结构中的有效性
   • 量化了νπ相互作用在能级结构中的关键作用
2. 实用价值：
   • 直接服务于NEEC实验争议的解决
   • 为核能存储和释放的新概念提供理论支持
   • 可指导寻找其他高多极同质异能态
3. 可复现性：
   • 模型空间和参数完全明确
   • 计算方法标准且成熟
   • 结果易于重复验证
4. 潜在应用：
   • 核钟开发（利用长寿命同质异能态）
   • 核电池技术
   • 基础核数据完善

适用场景

1. 理论研究：
   • 需要理解核同质异能态微观起源的研究
   • νπ相互作用效应的定量研究
   • 高多极跃迁机制的探索
2. 实验指导：
   • NEEC过程的实验设计和数据分析
   • 同质异能态寻找的候选核选择
   • 精密核谱学实验的理论支持
3. 方法学借鉴：
   • 简化模型空间的有效使用策略
   • 从邻近核确定相互作用参数的方法
   • 粒子-空穴对称性在核结构分析中的应用
4. 教学案例：
   • 壳模型应用的典型例子
   • 核能级结构和跃迁理论的实例
   • 理论与实验结合的范例

参考文献

本文引用的关键文献包括：

1. 同质异能态综述: Walker, Jain, Maheshwari (2024) - 核同质异能态专题
2. NEEC首次报道: Chiara et al., Nature 554, 216 (2018)
3. NEEC争议: Guo et al., PRL 128, 242502 (2022)
4. 经典理论: Auerbach & Talmi, Phys. Lett. 9, 153 (1964)
5. 前人计算: Hasegawa et al., Phys. Lett. B 696, 197 (2011)
6. 作者前期工作: Maheshwari et al., Phys. Rev. C 112, 024308 (2025)
7. 核结构理论: Talmi, "Simple Models of Complex Nuclei" (1993)
8. Pandya变换: Pandya, Phys. Rev. 103, 956 (1956)

总体评价：这是一篇高质量的核理论研究论文，通过简化而精巧的壳模型分析，深入揭示了$^{93}$Mo中罕见E4同质异能态的微观起源。物理图像清晰，定量分析透彻，与实验符合优异，并通过粒子-空穴对比验证了机制的正确性。论文不仅解决了基础核结构问题，还为NEEC这一前沿应用研究提供了重要的修正参数。主要不足在于模型空间的系统性验证和不确定度分析可进一步加强。该工作对核同质异能态物理和NEEC研究具有重要的学术价值和实用意义。