Interplay of quadrupole and octupole degrees of freedom in the Gd isotopes

技術要約：Gd同位体における四重極および八重極自由度の相互作用

問題提起
希土類領域における原子核は、集団的自由度と単一粒子自由度が相互作用する複雑な景観を呈している。特にガドリニウム（Gd, Z=64）同位体は、Z=64におけるサブシェル閉鎖や、中性子数（N）の増加に伴う球形から回転形への形状変化を含む、多様な構造現象を示す。現在の理論モデルにとっての大きな課題は、全同位体系列にわたって四重極（$\beta_2$）および八重極（$\beta_3$）の集団性を同時に記述することである。近年の実験データは、この領域における最大の八重極集団性が$^{150}\text{Gd}$（N=86）に位置している可能性を示唆しており、これは八重極マジックナンバーをN=88またはN=90付近に置いていた従来の仮定に疑問を投げかけている。さらに、X(5)臨界点対称性（振動型から回転型への形状相転移に関連するもの）は$^{154}\text{Gd}$のような核に対して提案されているが、他の理論的研究では、その最良の実装は$^{150}\text{Gd}$または$^{152}\text{Gd}$である可能性が示唆されている。四重極変形の進化と八重極集団性の傾向を同時に再現することは、既存のアプローチにとって困難な課題である。

手法
著者らは、ボール・ハミルトニアンに基づく四重極・八重極軸対称集団モデルを用いている。理論的枠組みは、集団的部分（$\hat{H}_{QO}$）と回転項（$C\hat{L}^2$）からなるハミルトニアン（式1）を利用している。集団ポテンシャル $U(\beta_2, \beta_3)$ は、新しい変数セットを用いて扱われる。それは、全体的な変形尺度 $\tilde{\beta}$ と、四重極および八重極変形の相対的な寄与を表す角度変数 $\phi$ である。

本手法の主な特徴は以下の通りである：

• ポテンシャルの構築： 角度ポテンシャル $u_L(\phi)$ は、放物線状の内側領域を導入することで、$\pm\phi$（または $\pm\beta_3$）の構成間の相互作用を許容するように修正されており、これにより、八重極の安定性または振動と一致する単一または二重の極小値を可能にしている。
• パラメータ・フィッティング： モデルパラメータ（$|\phi_c|$, $w_0$, $c$）は、質量数 $A = 148$ から $160$ のGd同位体の、交互のパリティを持つイラスト状態および最低正パリティ励起（$\beta$）バンドの実験的エネルギー準位へのフィッティングによって最適化される。
• 電磁的観測量： 遷移確率（E0, E1, E2, E3）は、$\tilde{\beta}$ および $\phi$ の変数に関する特定の積分を用いて計算される。スケーリングパラメータ（$s$ および $x$）は、理論結果を絶対値（ワイスコップ単位または適切な単位）に変換するために、利用可能な実験的遷移強度へのフィッティングによって決定される。

主要な貢献と結果
本研究は、$84 \le N \le 96$ のGd同位体系列の系統的な調査を提供し、以下の結果を得た：

1. エネルギー準位の再現： 本モデルは、全系列にわたる正および負のパリティ状態の両方について、実験的なエネルギー準位を正常に再現している。軽い核（$^{148}\text{Gd}$, $^{150}\text{Gd}$）における交互パリティ配列から、重い同位体において負パリティ状態が高エネルギー側にシフトする構造への遷移を捉えており、これは八重極変形から八重極振動への進化と一致している。
2. 四重極変形の進化： 結果は、$N=90$ での明確なジャンプによって強調される、四重極変形（$\beta_2$）の滑らかな進化を示している。このジャンプはX(5)臨界点対称性と一致しており、球形から回転形への形状相転移をマークしている。計算された変形は、$^{148}\text{Gd}$ の $\approx 0.01$ から $^{158}\text{Gd}$ および $^{160}\text{Gd}$ の $\approx 0.35$ へと増加している。
3. 八重極変形と集団性：
   • モデルは、最も軽い同位体である $^{148}\text{Gd}$ と $^{150}\text{Gd}$ についてのみ、非ゼロの静的な八重極変形（$\beta_3$）を予測する。
   • 重い同位体（$N \ge 88$）については、ポテンシャルはゼロの八重極変形における単一の極小値を示し、これは静的な変形ではなく、八重極ソフトモードまたは振動であることを示している。
   • E3強度： 計算された $B(E3)$ 遷移強度は、軽い同位体において増加傾向を示し、$^{152}\text{Gd}$（$N=88$）でピークに達した後、より重い核では減少する。絶対値には実験との定量的な偏差が見られるものの、モデルは $N=88$ 付近で最大となる八重極集団性の傾向を定性的に再現している。
4. 遷移確率：
   • E1遷移： モデルは、実験的なE1遷移強度に対して著しく高いレベルの一致を実現しており、これは著者らによって、既存の理論モデルではめったに達成されない結果であると述べられている。
   • E2遷移： モデルは、バンド内およびバンド間のE2遷移の一般的な傾向を再現しているが、絶対値を過大評価する傾向がある。また、一部のデータで見られる $N=92$ での $B(E2)$ 値の急激な減少を捉えている。
   • E0遷移： モデルは、E0強度の傾向（$^{154}\text{Gd}$ の臨界性と一致する $N=90$ での極大を含む）を再現しているが、絶対的な大きさについては過大評価している。

意義と主張
本論文は、少数の調整可能なパラメータを用いたフェノメノロジー的モデルを用いて、Gd同位体系列における四重極および八重極の集団性を「統一的な記述」を提供することを主張している。著者らは、本アプローチが形状相転移と八重極相関の複雑な相互作用を同時に再現することに成功したことを強調している。

具体的には、本論文は以下を強調している：

• 「八重極駆動数」の挙動の再現であり、E3強度のピークが $N=88$（$^{152}\text{Gd}$）にあることを予測しており、これは八重極マジックナンバーの理論的予測および最近の実験的示唆と一致している。
• 反転非対称（静的八重極変形）から反転対称（八重極振動）の核への転移を、同一の形式内で記述する能力。
• 通常、波動関数の混合に対して敏感であるためにモデリングが困難とされる、E1遷移強度の再現。

著者らは、より重い同位体や高スピン状態については、本予測を完全に検証するためにさらなる実験データが必要であると考えているが、本モデルは、この領域における構造進化を理解するための堅牢な枠組みを提供するものであると結論付けている。彼らは、同様の系統的な研究が隣接する同位体系列（Nd, Sm, Dy）においても必要であり、本モデルは、以前は概念実証の例としてのみ用いられていた重核の再検討を求めるものであると示唆している。