Theoretical calculations on half-lives of spontaneous one-proton radioactivity

技術要約：一陽子放出型放射性崩壊の半減期に関する理論計算

問題提起
中性子ドリップラインを超えた不安定核の研究は、極限的な系における核相互作用および核構造を理解するために不可欠である。中性子ドリップラインは依然として多くが未探索であるが、陽子ドリップラインについては広範な実験的および理論的な調査が行われてきた。陽子放射性崩壊は、負の一陽子分離エネルギーを持つ原子核で発生する自発的崩壊過程であり、開いた量子系、殻構造の進化、およびペアリング相関を調べるためのプローブとして機能する。様々な陽子放出核（$^{53}\text{Co}^m$ から $^{185}\text{Bi}$ まで）に関する実験データが存在するものの、核の変形を系統的に組み込み、半減期計算に対する特定の核入力（崩壊エネルギーおよび分光因子）への感度を評価する包括的な理論的枠組みが必要とされている。さらに、将来の実験的探索を導くために、$fpg$シェル（$30 < Z < 50$）における軽い陽子放出核に対する理論的予測も求められている。

手法
著者らは、陽子放出の半減期を計算するために、微視的なガモフ状態理論に基づく量子トンネルモデルを用いている。この手法の核心は、放出される陽子と娘核との間の二体ポテンシャルを構築することであり、そこでは核の変形を明示的に考慮している。

1. ポテンシャルの構築: 相互作用ポテンシャルは、変形したウッドサックス核ポテンシャル、トーマス型スピン軌道ポテンシャル、および拡張されたクーロンポテンシャルの組み合わせである。

   • 核ポテンシャル: ウッドサックス項およびスピン軌道項の半径パラメータは、四重極（$\beta_2$）および六重極（$\beta_4$）変形パラメータを用いた多重極展開を介して、極角 $\theta$ に依存する。ドリップライン核の強いアイソスピン非対称性を考慮するため、ウッドサックスポテンシャルにはアイソスピン依存の半径公式が用いられている。
   • クーロンポテンシャル: 変形核の場合、クーロンポテンシャルは単純な球近似ではなく、多重極展開法を用いて計算される。
   • ガモフ状態: 放出される陽子は、遠距離において外向きのクーロン波動関数境界条件を満たす、時間独立シュレディンガー方程式を満たすガモフ状態として扱われる。ポテンシャルの深さ（$V_0$）は、実験的な崩壊エネルギー（$Q$値）を再現するように調整される。

2. 半減期の計算:

   • 崩壊幅 $\Gamma(\theta)$ は、ガモフ状態の波動関数の外部領域における振幅を用いて計算される。これは、半古典的なWKB近似よりも微視的な記述を提供する。
   • 全崩壊幅は、方位角にわたる部分幅の平均によって得られる。
   • 半減期（$T_{1/2}$）は、$T_{1/2} = \hbar \ln 2 / (S_p \Gamma)$ という関係式から導出される。ここで $S_p$ は分光因子である。

3. データソース:

   • 崩壊エネルギー: 計算には、実験的な一陽子分離エネルギー（NNDCおよびAME 2020より）と、有限範囲ドロップレットモデル（FRDM）からの理論値の両方が用いられる。
   • 分光因子: 二種類の因子が用いられる。一つは相対論的平均場（RMF）理論から導かれたもの、もう一つは変形依存のフィッティング公式から得られたものである。
   • 変形パラメータ: 四重極および六重極パラメータは、FRDM計算から取得されている。

主な貢献

• 系統的な変形モデル: 本研究は、親核および娘核の両方に対して四重極および六重極変形を含む量子トンネルモデルを拡張し、変形下における陽子放出の半減期の微視的な計算を提供している。
• 包括的なデータセット: 著者らは、実験的に確認された陽子放出核（$50 < Z < 84$）および理論的に予測された放出核（$30 < Z < 50$）の統一された半減期データセットを編纂した。
• 感度分析: 崩壊エネルギー（実験値 vs FRDM）および分光因子（RMF vs フィット関数）の選択に対する計算された半減期の依存性を定量的に分析し、精度を達成する上で実験的な崩壊エネルギーが決定的な役割を果たすことを示した。
• 軽い陽子放出核への予測: 本論文は、$fpg$シェル（$Z=30$から$50$）における陽子放出核の理論的予測を提供しており、これには既に実験的に同定されているものの、この文脈では半減期の詳細が未測定である$^{89}\text{Rh}$ に対する具体的な予測が含まれている。

結果

• 検証: モデルは、代表的な事例として $^{151}\text{Lu}$ を用いて検証された。変形（$\beta_2 = -0.167, \beta_4 = -0.035$）とRMF分光因子を組み込むことで、計算された半減期（70.2 ms）は実験値（127 ms）と合理的な一致を示した。
• 精度およびRMS偏差:
  • 実験的な崩壊エネルギー（$Q_{\text{exp}}$）およびRMF分光因子（$S_p^{\text{RMF}}$）を用いた場合、計算された対数半減期と実験値の間の平方根平均二乗（RMS）偏差は $\sigma = 0.46$ である。
  • FRDMの崩壊エネルギーを使用すると、偏差は著しく増大し（$\sigma \approx 4.90$）、これはモデルが崩壊エネルギーの入力に対して非常に敏感であることを浮き彫りにしている。
  • モデルの性能は、解析的な半古典的手法や拡張されたガイガー・ヌッテルの法則などの他の理論的アプローチと同等である。
• より軽い陽子放出核（$Z < 50$）:
  • $^{89}\text{Rh}$ について、モデルはRMF分光因子を用いて274 ns、フィット関数を用いて972 nsの半減期を予測している。
  • 他の核（例：$^{60}\text{Ga}$, $^{63}\text{As}$, $^{72}\text{Rb}$）の予測は、重い放出核と比較して、一般に短い半減期（ナノ秒からマイクロ秒の範囲）を示している。
  • 結果は、固定された角運動量に対して、崩壊エネルギーが増加するにつれて半減期が減少することを裏付けており、これはガイガー・ヌッテルの法則（$\log_{10} T_{1/2} = a Q^{-1/2} + b$）と一致している。

意義
本論文は、その主要な意義が「陽子ドリップラインに関連する将来の実験的研究のための有用なリファレンスを提供すること」にあると主張している。$^{60}\text{Ga}$ から $^{185}\text{Bi}$ までの原子核の計算および予測された半減期の包括的なデータセットを編纂することで、著者らは実験学者が陽子過剰核を特定し、特徴付けるのを支援することを目指している。本研究は、理論モデルが実験的な崩壊エネルギーを用いることで精密な予測を達成できる一方で、これらの手法をより軽い未観測の核へと拡張することが、ナノ秒からピコ秒の時間スケールを探索できる高度な技術（インフライト崩壊測定など）を用いる実験に対して重要な指針を与えることを強調している。さらなる研究を促進するため、このデータセットはオープンに公開されている。