Determining the dynamic deformation of $^{140}$Ce by constraining coupled-channels parameters for fusion

技術要約：結合チャネル・パラメータの制約による$^{140}$Ceの動的変形の決定

問題提起
核形状パラメータ（四重極 $\beta_2$、八重極 $\beta_3$、および六重極 $\beta_4$）の抽出は、核物理学における中心的な課題である。融合反応は、結合チャネル（CC）解析を通じて核変形に対する敏感なプローブとして機能するが、現在の手法には重大な限界がある。具体的には、パラメータのチューニングが個々の反応系に限定されることが多く、異なる入射粒子と標的の組み合わせに適用した場合、一貫性のない、あるいは非物理的な結果を招く。さらに、融合障壁分布（$D$）を導出するための従来の二階微分（DD）法は、ステップサイズに対して非常に敏感であり、高エネルギー域での精度が低いため、抽出される構造情報の信頼性を制限している。振動型核として知られる$^{140}$Ceについては、これまでの研究において、使用される実験的プローブや理論モデルによって、得られる変形パラメータ（例：$\beta_2$ が 0.003 から 0.15 まで変動）が著しく異なるという結果が出ている。

手法
本研究では、新しい実験データ、堅牢な障壁分布の解析フレームワーク、および複数の反応系にわたる厳密な統計解析を組み合わせた体系的なアプローチを採用している。

1. 実験データの取得:

   • $^{16}$O + $^{140}$Ce: 新しい融合断面積（$\sigma_{fus}$）の測定を、インド大学間加速器センター（IUAC）の重イオン反応解析器（HIRA）を用いて行った。データは、バス障壁（Bass barrier, 57–76 MeV）の約12.4%下限から約17%上限までの範囲で収集された。
   • $^{36}$S + $^{140}$Ce: 文献から既存の融合データを利用した。
   • $^{28}$Si + $^{140}$Ce: 既存の励起関数を検証用のテストベッドとして使用した。

2. 障壁分布解析:

   • 融合障壁分布（$D$）は、従来のDD法と**ガウス解析レシピ（Gaussian analytic recipe）**の2つの手法を用いて導出された。後者は、ステップサイズへの依存性がなく、$\sigma_{fus}$ の大きな不確かさに対しても耐性があるため、完全で明確に定義されたスペクトルを提供できるとして選択された。

3. 結合チャネル（CC）計算:

   • 計算は、Woods-Saxon 核ポテンシャルを用いた ccfull コードを使用して行われた。
   • 標的モデリング: $^{140}$Ce は、一フォノン的な四重極（$2^+$）および八重極（$3^-$）励起を含む純粋な振動型核としてモデル化された。
   • パラメータ空間: $\beta_2$ および $\beta_3$ を 0.00 から 0.60 まで系統的に変化させた。
   • 検証: $^{28}$Si + $^{140}$Ce 系については、入射粒子の固有変形（振動型および回転型の両方として扱う）と、正の Q 値を持つ二中性子（$2n$）転送チャネルを含む、改訂された結合スキームを用いて解析を行った。

4. 統計解析:

   • $\chi^2$ 最小化: 解析的な障壁分布と理論的な CC 結果との差を最小化することにより、変形パラメータを抽出した。
   • ベイズ解析: $^{16}$O および $^{36}$S の両方の系における $\beta_2$ と $\beta_3$ の事後確率分布を決定するために、マルコフ連鎖モンテカルロ（MCMC）フレームワークを使用した。
   • モデル平均化: 2つの系の結果を組み合わせるために、独立したベイズモデル平均化（NBMA）を適用し、関連する不確かさを伴う最適値を算出した。

主な結果

• 障壁分布: ガウス解析レシピは、$^{36}$S + $^{140}$Ce 系に対して、 $E \simeq 104$ MeV に主要なピークを持ち、広範な高エネルギー構造を示す、明確に定義された障壁分布を生成することに成功した。これは、高エネルギー域で偽の3ピーク構造を生じさせた従来の DD 法とは対照的であり、高エネルギーにおける DD 法の限界を浮き彫りにしている。
• 抽出された変形パラメータ:
  • $^{16}$O + $^{140}$Ce: $\beta_2 = 0.06^{+0.07}_{-0.06}$ および $\beta_3 = 0.18^{+0.02}_{-0.04}$。
  • $^{36}$S + $^{140}$Ce: $\beta_2 = 0.09^{+0.02}_{-0.02}$ および $\beta_3 = 0.18^{+0.01}_{-0.01}$。
  • NBMA による結合結果: 本研究は、$^{140}$Ce の動的変形パラメータを $\beta_2 = 0.09 \pm 0.03$ および $\beta_3 = 0.18 \pm 0.02$ と決定した。これらの値は、両方の入射粒子系において一致している。
• $^{28}$Si + $^{140}$Ce における検証:
  • 抽出されたパラメータ（$\beta_2=0.09, \beta_3=0.18$）を $^{28}$Si + $^{140}$Ce 系に適用したところ、CC 計算は融合励起関数と解析的障壁分布の両方を正常に再現した。
  • 正の Q 値を持つ $2n$ ピックアップ・チャネルの導入が、障壁下の融合増強および障壁分布の形状を再現するために不可欠であることが判明した。
  • $^{28}$Si を振動型または剛体回転型として扱うことによる結果の違いは無視できる程度であり、抽出された標的パラメータの堅牢性を示唆している。
• 文献との比較: 抽出された $\beta_2$ 値は、Pritychenko ら（$\beta_2 = 0.10$）と密接に一致しているが、$\beta_2 \approx 0.15$ を示唆する以前の融合ベースの研究とは異なっている。抽出された $\beta_3$ は、先行文献で報告されている値よりも大きい。

意義および主張
本論文は、複数の反応系によって制約されることで、抽出される変形パラメータの一貫性と堅牢性が示されることを実証しており、単一系解析で見られる曖昧さに対処していると主張している。また、本研究は、融合障壁分布を導出するための効果的かつ敏感なツールとしてガウス解析レシピを検証しており、実験的な不確かさを扱う上で従来の DD 法よりも優れていることを証明している。さらに、単一の $^{140}$Ce 変形パラメータを用いて、3つの異なる入射粒子系（$^{16}$O, $^{36}$S, $^{28}$Si）のデータを成功裏に再現できたことは、採用された手法の信頼性を強調している。著者らは、このようなマルチシステムによる制約が、物理的に正当化された核構造パラメータを導出するために極めて重要であり、ガウス解析アプローチが固有の核変形を探索するための強力なプローブとして機能すると述べている。