Physics-informed interpretable forms for dipole resonances and quantities relevant for astrophysics

技术摘要：用于偶极共振及天体物理相关量的物理启发式可解释形式

问题陈述
同位旋巨偶极共振（IVGDR）是一种关键的集体激发模式，为理解核结构、动力学及核对称能提供了重要见解。虽然微观模型（如 TDHF、RPA）能够较为准确地预测 IVGDR 中心能量（$E_G$），但它们往往难以精确重现共振宽度（$\Gamma_G$），且在大规模应用中计算成本过高。相反，经典流体动力学模型虽提供了简单的解析表达式，但缺乏捕捉同位旋不对称性、核形变及壳效应等结构效应的精度。此外，现有的数据驱动机器学习方法通常缺乏物理可解释性。因此，需要改进的解析公式，既能保持清晰的物理可解释性，又能准确重现整个核图谱上的 IVGDR 系统性规律，以支持大规模天体物理计算（如辐射俘获过程）。

研究方法
作者采用贝叶斯推断框架从实验光吸收数据中提取并优化 IVGDR 参数。

1. 数据来源： 分析利用了来自 EXFOR 数据库的 81 个球形和 60 个轴向形变核的 $\gamma$ 射线光谱。
2. 建模： 光吸收截面使用简化洛伦兹（SLO）分布进行建模。作者提出了具有物理启发性的新函数形式来描述中心能量（$E_G$）和宽度（$\Gamma_G$），其中明确引入了：
   • 同位旋不对称性： 在中心能量标度中加入了一个依赖于 $I = |N-Z|/A$ 的项，扩展了传统的 Goldhaber-Teller ($A^{-1/6}$) 和 Steinwedel-Jensen ($A^{-1/3}$) 模型。
   • 核形变： 对于轴向形变核，中心能量和宽度根据由实验 IVGDR 强度分布分裂得出的四极形变参数 $\beta_2$ 进行调整。
   • 壳效应： 宽度公式中的阻尼系数 ($k_s$) 经过精细化处理，纳入了对基态壳修正能 ($\delta W$) 的依赖性，从而解决了幻数附近的局部波动问题。
3. 优化： 使用贝叶斯推断来确定模型参数（对于能量为 $a_0, a_1, a_I$；对于宽度为 $k_s, \delta$）的后验概率分布，并量化它们之间的相关性。
4. 天体物理应用： 提取的参数用于约束体积对称能（$K_v$）和表面对称能（$K_s$）系数。随后，这些参数被用于预测核偶极极化率（$\alpha_D$），并利用各种 Skyrme 能量密度泛函（EDFs）研究其与中子皮厚度（$r_{skin}$）的相关性。

核心贡献

• 可解释参数化： 本文提出了紧凑的解析表达式，用于描述 IVGDR 中心能量和宽度，成功重现了球形及形变核的实验系统性。
• 结构集成： 与以往的全局拟合不同，这些形式显式地集成了同位旋不对称性、形变和壳效应，从而提供了更稳健的核图谱描述。
• 贝叶斯约束： 研究直接从 IVGDR 数据出发，为核对称能系数（$K_v, K_s$）提供了严谨的贝叶斯约束。
• 中子皮预测工具： 作者建立了偶极极化率与对称能之积（$\alpha_D J$）与中子皮厚度（$r_{skin}$）之间的线性相关性，为估算重核的中子皮厚度提供了一个实用的公式。

结果

• 中心能量： 提出的表达式 $E_G(A, I) = (a_0 A^{-1/3} + a_1 A^{-1/6})(1 - a_I I^2)$ 捕捉到了中心能量微妙的同位旋依赖性，在 $A < 50$ 的核素中比传统形式表现出更好的实验数据符合度。分析显示，$A^{-1/3}$ 与 $A^{-1/6}$ 系数之间存在强负相关（$r \approx -0.99$），表明存在共同的底层物理机制。
• 宽度与壳效应： 在阻尼系数 $k_s$ 中引入壳修正，显著改善了对实验宽度的重现，特别是在 $k_s$ 被抑制的幻数附近。对于形变核，研究发现必须使用不同的幂律指数（$\delta_\parallel \approx 3.66, \delta_\perp \approx -1.06$）来描述宽度各向异性，其表现优于单指数模型。
• 对称能与极化率： 优化的对称能系数（$K_v \approx 36.26$ MeV, $K_s \approx 57.64$ MeV）与近期测量值吻合良好。预测的偶极极化率（$\alpha_D$）与实验数据匹配良好，其中基于方程 (18) 的模型提供上限，而液滴模型（方程 20）提供下限。
• 中子皮厚度： 在重核（$A > 150$）中，观察到 $\alpha_D J$ 与 $r_{skin}$ 之间存在强线性相关性（方程 18 的 $r > 0.80$，方程 20 的 $r > 0.90$）。将此关系应用于 $^{208}\text{Pb}$ 时，得到的中子皮厚度与 PREX 实验协作组及电偶极响应研究的结果重叠。

意义与主张
作者声称，其研究结果为“大规模计算与核天体物理相关的辐射俘获过程有关的参数提供了有价值的输入”。具体而言，导出的解析形式提供了一种计算高效且具有物理动机的 IVGDR 强度函数描述，这对于模拟天体物理环境（如中子星合并）中的光诱导及辐射俘获反应速率至关重要。

论文强调，提取的对称能参数与独立的偶极极化率测量结果一致，从而加强了所推导的核状态方程的可靠性。此外，建立的 $\alpha_D J$ 与 $r_{skin}$ 之间的相关性为跨广泛质量范围估算中子皮厚度提供了一条“实用路径”，这对于约束对称能的密度依赖性以及理解中子星结构至关重要。这项工作被视为连接 IVGDR 可观测物理量与基本同位旋性质的桥梁，提升了对极端条件下天体物理情景的理解。