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$）の関数形式を提案している：
   • アイソスピン非対称性： 従来のGoldhaber-Tellerモデル（$A^{-1/6}$）やSteinwedel-Jensenモデル（$A^{-1/3}$）を超えて、重心エネルギーのスケーリングに$I = |N-Z|/A$に依存する項を追加した。
   • 核変形： 軸対称変形核については、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}$）との間の線形相関を確立し、重い核の$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の観測量と基本的なアイソベクター特性との間の架け橋として提示されており、極限状態における天体物理学的シナリオの理解を深めるものである。