Extracting nuclear charge radii from binding energies: a single-parameter empirical formula with structural corrections

技术摘要：从结合能中提取核电荷半径

问题陈述
核电荷半径（$r_{ch}$）是表征质子电荷空间分布的基本可观测物理量，而每核结合能（$\epsilon = B/A$）则量化了有限原子核的平均结合强度。尽管这两个物理量从互补的角度描述了同一个多体系统（空间延伸度 vs. 能量结合能），但经验半径公式通常依赖于几何缩放（$A^{1/3}$）并辅以人为修正，而微观模型（如 HFB、RMF、DFT）则高度依赖于特定的有效相互作用，且计算强度极大。此外，与质量数据相比，实验电荷半径数据在富中子和富质子区域较为稀缺。作者指出，由于 $\epsilon$ 编码了关于饱和、同位旋和壳结构的信息，它应当作为一种具有物理动机的输入，用于提取电荷半径尺度，从而为半径估计提供一条新的经验途径。

方法论
作者提出了一种单参数经验公式 BECR1p，该公式直接从每核结合能（$\epsilon$）而非仅仅通过质量数缩放来推导领先的电荷半径尺度。该方法分为三个阶段：

1. 宏观项 (BECR$_{mac}$)：基于体积-表面耦合图景构建了一个唯象项。有效长度尺度被推导为与 $(N/\epsilon)^{1/5}$ 成正比，其中 $N$ 是中子数。为了解释 $r_{ch}$ 探测的是质子而非总物质，减去了一个与 $- \frac{\sqrt{3}}{2} \frac{N}{Z}$ 成正比的同位旋修正项。仅靠该宏观项在 893 个实验半径上的均方根偏差（RMSD）为 0.0345 fm。
2. 局部结构修正：为了恢复被宏观项平滑掉的局部结构细节，作者引入了以下修正：
   • 壳效应与奇偶效应：通过基于接近幻数（magic numbers）的 Casten 类因子（$e_{sh}$）和按 $A^{-5/6}$ 缩放的配对因子（$e_{pair}$）进行建模。
   • 有限尺寸效应：一个对于轻核至关重要的 $1/N$ 项。
   • 形变：一个涉及形变参数 $\beta_2$ 和 $\beta_4$（取自 WS3.3 质量表）的乘性因子。
3. 参数简化：完整的公式最初包含四个拟合参数（$a_1, a_2, c_0, r_\beta$）。通过全局优化，作者发现这些参数之间存在强烈的经验相关性，允许将它们表示为一个单一的有效参数 $\kappa$。这使得公式简化为 BECR1p 形式。
4. 残差修正 (AKRR)：为了探测残差中剩余的系统信息（例如高阶壳效应或形状共存），作者应用各向异性核岭回归（AKRR）作为后处理诊断工具。这并非独立的预测器，而是利用 $(N, Z)$ 作为输入对 BECR1p 的残差进行的修正。

主要贡献与结果

• 单参数精度：仅含一个拟合参数（$\kappa = 0.7190$）的 BECR1p 公式在 893 个精选实验电荷半径集上的 RMSD 达到了 0.0138 fm。其表现优于代表性的经验公式，如 WS3.3-R (0.0170 fm) 和 Sheng (0.0247 fm)。
• 结构再现：该公式成功再现了整个核图谱中复杂的局部特征，包括：
  • 幻数处的壳层转折（kinks）（例如 Ca, Pb）。
  • 奇偶交替现象。
  • 与形状转变相关的剧烈半径跳变（例如 $N \approx 60$ 附近的 Sr 链）。
  • 形变区域（Cs, Gd）的趋势。
• 预测能力：利用该公式，作者预测了整个核图谱中 11,205 个原子核 的电荷半径。
• 残差分析：将 AKRR 应用于 BECR1p 残差，使留一法交叉验证（LOOCV）的 RMSD 进一步降低至 0.0081 fm，其中 98.4% 的点落在 $\pm 0.02$ fm 之内。核函数的各向异性（$\sigma_1 > \sigma_2$）表明，残差相关性沿中子数（$N$）的变化比沿质子数（$Z$）的变化更为平滑。
• 局限性：在对九个被排除在拟合集之外的同位素链（Tb, Tm, Lu, Re, Po, Rn, Fr, Ra, Cm）进行测试时，该模型能够再现总体趋势，但表现出系统性偏差（RMSD $\approx$ 0.04 fm）。这些偏差归因于这些链中的实验归一化歧义以及复杂的理论结构（形状共存、八极形变），这些是全局经验公式难以捕捉的。

意义与主张
本文主张，每核结合能（$\epsilon$）包含了足以建立核电荷半径领先尺度的内在结构信息。通过利用这种相关性，BECR1p 公式提供了一种互补且物理透明的电荷半径建模方法，而不依赖于微观模型所需的复杂有效相互作用。

作者对其方法的适用范围保持了谦逊的态度：

• BECR1p 被呈现为一种深化对结合能-半径相关性理解的经验工具，而非微观模型的替代品。
• 对于向滴线核或远离稳定性的奇异区域的外推，作者指出微观模型（如 HFB-29）仍然是“更具物理鲁棒性的指导”。
• AKRR 修正被明确界定为一种用于已知区域内插的“残差分析器”；由于其依赖于附近的训练覆盖范围，应谨慎将其应用于数据稀疏的区域。

总之，这项工作证明了基于结合能相关性的单参数公式可以实现极高的精度，从而在实验数据匮乏的区域为半径估计提供可靠的基准，同时也承认了全局系统规律在面对复杂核结构时的局限性。