Static Fission Properties of Even-Even Actinides within the Warsaw Macroscopic-Microscopic Model Using Fourier-over-Spheroid Parameterization

ワルシャワ・マクロミクロモデルとフーリエ過球面パラメータ化を用いた大規模計算により、トリウムからカリホルニウムまでの偶々アクチノイド核の分裂障壁高さを高精度に評価し、IAEA RIPL-4 データセットとよく一致する結果を得るとともに、トリウム同位体では第 3 極小の存在を確認したがウランやプルトニウムでは確認されなかったことを報告しています。

要約

🌋 原子核の分裂：巨大な「粘土の山」の物語

まず、原子核（原子の中心にある粒）を想像してください。これは、**「電気的な反発力」と「表面張力」という 2 つの力がせめぎ合っている、「電気を持った粘土の球」**のようなものです。

通常、この粘土の球は丸くて安定しています。しかし、何らかのきっかけで伸び縮みし始めると、やがて「くびれ」ができ、2 つの塊に分かれて分裂します。これを核分裂と呼びます。

このとき、粘土が分裂するまでに越えなければならない**「高い壁（山）」のことを「核分裂の障壁（バリア）」**と呼びます。この壁の高さが正確にわかれば、原子炉の設計や、新しい元素の発見に役立ちます。

🗺️ 研究の目的：「完璧な地形図」を描く

これまでの研究では、この「粘土の山」の地形図（エネルギーの地図）を描くとき、いくつかの**「仮定」や「単純化」**をしていました。それは、地図を作るのが大変だからです。

しかし、今回の研究（ワルシャワのチーム）は、**「もっと細かく、もっと立体的に」**地図を描くことにしました。

1. 新しい道具：「 Fourier-over-Spheroid（フォス）」

彼らは、新しい**「粘土の形を表現する言語」**を使いました。

• 昔の方法： 粘土を「球」に「いくつかの突起」を足して表現していました。でも、これだと、分裂の直前にあるような「くびれた奇妙な形」を正確に表現するのが難しかったのです。
• 今回の方法（フォス）： 粘土の表面を、**「 Fourier（フーリエ）変換」**という数学的な波の重ね合わせで表現しました。
  • イメージ： 昔は「丸いボールに、少しの凹凸」しか描けなかったのが、今回は**「しなやかな波のように、くびれや膨らみを自由自在に描ける」**ようになったのです。これにより、分裂の瞬間の「くびれた首」のような複雑な形も、くまなく捉えられるようになりました。

2. 巨大な計算：「1 億 3000 万の地点」

彼らは、この新しい地図を描くために、**「1 億 3000 万」**もの地点を計算しました。

• イメージ： 山岳地帯のハイキングで、これまで「主要な道」しかチェックしていなかったのが、**「森の中の草むら、岩の隙間、すべての小道」**まで、1 メートルごとに調べるようなものです。
• これにより、「地図にない隠れた谷（エネルギーの低い場所）」や「見落としがちな小さな山（障壁）」を、見逃さずに発見できるのです。

🔍 発見された驚きの事実

この精密な地図を使って、**「アクチノイド（ウランやプルトニウムなど、重い元素）」**の核分裂を調べたところ、いくつかの面白いことがわかりました。

① 「第 3 の谷」の謎

長い間、物理学者たちは「重い原子核の分裂の途中に、**『第 3 の谷（超変形した安定した状態）』**という、深い穴があるのではないか？」と議論していました。

• トロン（Th）という元素： この研究では、トロンには**「浅いけれど、確かに存在する第 3 の谷」**が見つかりました。まるで、山を越える途中に、一休みできる小さな窪みがあるようなものです。
• ウラン（U）やプルトニウム（Pu）： しかし、これらより重い元素では、「その第 3 の谷は全く見当たりませんでした」。
  • なぜ？ 以前の実験データでは「ウランにも深い谷があるはずだ」と言われていましたが、今回の「超精密な地図」では、そこはただの「なだらかな坂」だったのです。
  • 意味： これは、「実験データの解釈」か「古い地図の描き方」に問題があった可能性を示唆しています。

② 壁の高さの予測精度

計算された「壁の高さ（障壁）」は、実験で測定された値と非常に良く一致していました（平均して 1 メガ電子ボルト以下の誤差）。

• これは、今回の「新しい地図の描き方（フォス）」が、原子核の振る舞いを非常に正確に再現できていることを証明しています。

🏁 まとめ：なぜこれが重要なのか？

この研究は、単に数字を計算しただけではありません。

1. より正確な「未来の地図」： 新しい数学的な方法（フォス）を使うことで、原子核の分裂の道筋を、これまで以上に正確に予測できるようになりました。
2. 謎の解決： 「第 3 の谷」のような長年の議論に対し、「元素の種類によって、その谷の有無や深さは違う」という答えを出しました。
3. 実用への貢献： 原子力発電所の安全性向上や、新しい超重い元素の合成、さらには宇宙での元素の生成（r 過程）を理解する上で、この「正確な地形図」は不可欠なツールです。

一言で言えば：
「これまでの地図では見落としていた、原子核の分裂という『複雑な山岳地帯』の、くまなく詳細な地形図を描き上げ、そこで隠れていた『小さな谷』の正体に迫った研究」です。

彼らの努力により、人類は原子核という「見えない世界」の風景を、これまで以上に鮮明に捉えることができるようになりました。

技術概要

以下は、提示された論文「Static Fission Properties of Even–Even Actinides within the Warsaw Macroscopic–Microscopic Model Using Fourier-over-Spheroid Parameterization（ワルシャワ・マクロ・ミクロモデルにおけるフーリエ・オーバー・スフェロイドパラメータ化を用いた偶数 - 偶数アクチノイド核の静的核分裂特性）」の技術的な要約です。

1. 研究の背景と課題 (Problem)

原子核の分裂障壁（fission barrier）は、原子核の安定性、超重元素の合成、原子炉設計、および天体物理学的 r 過程における元素合成において決定的な役割を果たす重要な物理量です。

• 課題: 従来のマクロ・ミクロモデルにおける核分裂障壁の計算では、核形状のパラメータ化（記述方法）が結果に大きく影響します。特に、核分裂過程における複雑な形状変化（極端な伸長や首のくびれなど）を正確に記述するために、高次元の形状空間を網羅的に探索する必要があります。
• 既存の限界: 従来の球面調和関数（多極モーメント）による展開では、非物理的な形状が含まれたり、計算コストが高すぎたり、あるいは特定の形状領域（特に分裂直前や超変形領域）での精度が不足する可能性があります。また、第三の極小（超変形極小）の存在とその深さについては、理論モデルと実験データ（中性子誘起反応の断面積など）の間で長年の議論が続いています。

2. 手法とアプローチ (Methodology)

本研究では、ワルシャワ・マクロ・ミクロモデル（Warsaw macroscopic–microscopic model）を採用し、以下の革新的手法を用いて、Th（トリウム）から Cf（カリホルニウム）までの偶数 - 偶数アクチノイド核の静的核分裂特性を系統的に研究しました。

• 形状パラメータ化（FoS）:
  • 従来の多極展開の代わりに、「フーリエ・オーバー・スフェロイド（Fourier-over-Spheroid: FoS）」パラメータ化を採用しました。
  • 軸対称な核形状を、球体からの偏差をフーリエ級数で記述する 5 次元パラメータ空間 $(c, a_3, a_4, a_5, a_6)$ で定義します。
    • $c$: 伸長パラメータ
    • $a_3$: 反射非対称変形（質量非対称性）
    • $a_4$: 首（neck）の発達
    • $a_5, a_6$: より高次の微細構造
  • このパラメータ化は、基底状態から分裂断片に至るまでの形状進化を物理的に直感的かつ柔軟に記述できます。
• 網羅的なポテンシャルエネルギー面（PES）の探索:
  • 各核種に対して、約 $1.3 \times 10^8$点（$101 \times 26 \times 51 \times 31 \times 31$）という極めて高密度な 5 次元変形グリッドを計算しました。
  • 補間を一切行わず、離散的なグリッド点から直接極小点と鞍点を特定することで、数値的不確実性を排除しました。
• 障壁高さの抽出:
  • 「浸水法（Immersion Water Flow: IWF）」というアルゴリズムを改良・適用し、5 次元空間における真の鞍点（saddle point）を特定し、内側障壁と外側障壁の高さを抽出しました。
• モデル構成:
  • マクロ部分：Yukawa-plus-exponential モデル。
  • ミクロ部分：Strutinsky 殻補正法（変形 Woods-Saxon ポテンシャル）および BCS 理論による対相関。
  • 非軸対称変形（triaxiality）は今回は考慮していませんが、将来的な課題として挙げられています。

3. 主要な成果と結果 (Key Contributions & Results)

• 基底状態質量と障壁高さの精度:

  • 計算された基底状態の質量過剰（mass excess）は、実験値（AME2020）と非常に良く一致しており、平均偏差は +0.114 MeV、RMS 誤差は 0.146 MeV でした。
  • 核分裂障壁高さも、IAEA の RIPL-3 および RIPL-4（未発表だがレビュー済み）データセットや Smirenkin による評価値と比較され、平均偏差が 1 MeV 未満という良好な一致を示しました。
    • 内側障壁: 実験値に対して平均 +0.46 MeV 過大評価（RMS 1.53 MeV）。これは非軸対称変形の欠如が主な原因と考えられます。
    • 外側障壁: 実験値に対して平均 +0.39 MeV 過大評価（RMS 0.64 MeV）。内側障壁よりも実験値との一致が良く、非軸対称性の影響が小さい領域であることを示唆しています。

• 第三極小（超変形極小）の存在に関する知見:

  • トリウム（Th）同位体: 浅いながらも明確な「第三の極小（hyperdeformed minimum）」が観測されました。特に $^{230}\text{Th}$ では、$c \approx 1.6, a_4 \approx 0.2$ 付近に深さ約 0.5 MeV の極小が存在しました。
  • ウラン（U）およびプルトニウム（Pu）同位体: 同様の FoS パラメータ化を用いた計算では、U や Pu の重核種において、この第三極小は完全に消失し、外側鞍点を過ぎるとポテンシャルエネルギーは単調に減少する（分裂へ向かう）ことが示されました。
  • パラメータ化の影響: この結果は、第三極小の存在や深さが、マクロ・ミクロモデルの枠組みそのものよりも、核表面を記述する幾何学的なパラメータ化の完全性に敏感に依存することを示しています。

• 高次変形項の重要性:

  • 3 次元パラメータ空間 $(c, a_3, a_4)$ と 5 次元空間 $(c, a_3, a_4, a_5, a_6)$ を比較した結果、高次項（$a_5, a_6$）の導入は、基底状態や鞍点近傍のエネルギーを 1 MeV 以上低下させる重要な効果を持つことが確認されました。特に重アクチノイド（Cm など）では、基底状態の正確な記述にも高次項が不可欠であることが示されました。

4. 意義と結論 (Significance & Conclusion)

• 理論的妥当性の確認: 高密度な 5 次元グリッドと FoS パラメータ化を用いることで、パラメータの追加調整なしに、基底状態から分裂経路に至るまで一貫性のある高精度な記述が可能であることを実証しました。
• 第三極小論争への寄与: 実験的な核分裂断面積の解析（特に軽ウラン同位体）では第三極小の存在が仮定されることが多い一方、本研究のような厳密な静的計算では U 核で第三極小が見られないという矛盾を浮き彫りにしました。これは、第三極小の存在が形状パラメータ化の選択に依存するか、あるいは動的効果（遠心力障壁の変形依存性など）が静的ポテンシャルとは異なる振る舞いを生んでいる可能性を示唆しています。
• 将来への展望:
  • 本研究で得られた静的ポテンシャルエネルギー面は、多次元ランジュバン方程式やブラウン運動シミュレーションによる動的核分裂過程（断片質量分布、分裂時間スケールなど）の計算への入力として利用可能です。
  • 将来的には、奇数核への拡張、非軸対称変形（トライアキシアル性）の導入、およびマクロエネルギー関数（LSD モデルや ISOLDA）の比較検討を通じて、さらに精度の高いモデル構築が進められる予定です。

総じて、本研究は、核分裂障壁の予測において形状パラメータ化の重要性を再確認し、特に重アクチノイド領域における第三極小の性質について、実験データとの整合性を高めるための新たな視点を提供する重要な成果です。