A microscopic analysis of sub-barrier photo-induced fission in $^{236}$U$(\gamma,f)$ based on the non-equilibrium Green function method

本論文では非平衡グリーン関数法を用いて$^{236}$U の光誘起核分裂を微視的に解析し、実験データを再現するとともに、分裂確率が第一固有チャネルに支配されていることを示すことで、ボーア・ウィーラーの遷移状態説を微視的観点から裏付けた。

要約

1. 研究の舞台：「不安定な城」と「光のハンマー」

まず、ウラン 236という原子核を想像してください。これは、もともと不安定で、少しの刺激で崩壊しようとしている「巨大な城」のようなものです。

通常、この城を壊すには、壁（核分裂の障壁）を越えるだけの大きな力が必要です。しかし、この研究では、**「壁を越える力よりも少し弱い光（ガンマ線）」**を当てています。

• 壁を越える力：大きな岩を投げて城を壊すようなもの。
• 今回の光：小さな石を投げて、城の壁を「トンネル」のようにすり抜けて壊すようなもの。

この「壁をすり抜ける現象」を**「量子トンネル効果」**と呼びますが、これを正確に計算するのは非常に難しいのです。

2. 使われた道具：「非平衡グリーン関数（NEGF）」という新しい地図

これまでの研究では、この現象を計算する方法にいくつかの限界がありました。

• 古い地図（TDDFT）：城が崩壊する「後の様子」は描けるが、壁をすり抜ける「瞬間」は描けない。
• 別の地図（TDGCM）：壁をすり抜けることは描けるが、光が当たった「入り口」からの動きを正確に描くのが苦手。

そこで、この論文の著者（K. Uzawa 氏）は、電子工学（半導体など）で使われている「非平衡グリーン関数（NEGF）」という手法を、原子核の世界に持ち込みました。

• アナロジー：
  これまで原子核の分裂を計算する人たちは、「川の流れ」をシミュレートすることに熱中していました。しかし、NEGF は**「電子が回路をどう流れるか」を計算するプロの技術です。
  著者は、この「電子の流れる道」を計算する高度な技術を使って、「光が当たった入り口」から「分裂して飛び出す出口」まで、原子核の中をどうやって粒子が移動するかを、「光の道筋」のように精密に追跡**しました。

3. 実験の結果：「壁をすり抜ける」現象の再現

著者は、この新しい計算方法を使って、**「5.0 メV から 6.0 メV のエネルギーを持つ光」**をウランに当てたとき、どれくらい分裂が起きるかを計算しました。

• 実験との比較：
  実際の実験データ（過去の研究者たちが測った値）と比較すると、「壁を越えるエネルギー（5.7 メV）より低い領域」でも、計算結果が実験とよく一致していました。
  特に、エネルギーが低い（壁を越えるのが難しい）領域でも、計算が実験をうまく再現できたことは、この新しい方法の信頼性を示しています。

4. 最大の発見：「一本の道」が支配している

この研究で最も面白い発見は、**「分裂の道筋」**についての分析です。

• アナロジー：
  光が当たって分裂する際、原子核の中には「何万通りもの道筋（経路）」が存在するはずです。しかし、NEGF という方法で分析したところ、「その何万通りもの道筋のうち、たった 1 本の道（第 1 固有チャネル）」が、分裂の 99% を担っていることがわかりました。
  他の道筋はほとんど無視できるほど小さく、**「一本の太い幹道」**だけが、光のエネルギーを分裂に変換しているのです。

• 意味すること：
  これは、昔から言われていた**「ボーア・ホイーラーの遷移状態モデル」**という理論が、ミクロな視点から見ても正しいことを裏付ける結果です。
  「分裂するときは、決まった一本の『通り道』を通る」というイメージが、実は物理的に正しかったのです。

5. まとめ：なぜこの研究が重要なのか？

この論文は、単に「ウランが割れる」ことを計算しただけではありません。

1. 新しい計算手法の確立：電子工学の高度な計算手法を、原子核物理に応用し、成功させました。
2. 壁を越える現象の解明：エネルギーが足りない状態（サブバリア領域）でも、分裂がどう起きるかを説明できました。
3. シンプルさの発見：複雑に見える原子核の分裂も、実は「一本の道」が支配しているという、驚くほどシンプルな構造を持っていることがわかりました。

結論として：
この研究は、**「複雑怪奇に見える原子核の分裂という現象を、新しい計算の『メガネ』で見ると、実は一本の明確な道筋で説明できる」**と示した、画期的な一歩です。

将来、この手法を使えば、実験が難しい「中性子過剰な原子核」の挙動も予測できるようになり、原子力エネルギーや宇宙での元素合成（r プロセス）の理解がさらに進むことが期待されています。

技術概要

以下は、提供された論文「A microscopic analysis of sub-barrier photo-induced fission in 236U(γ, f) based on the non-equilibrium Green function method」に基づく詳細な技術的サマリーです。

論文概要

タイトル: 非平衡グリーン関数法に基づく 236U(γ, f) におけるサブバリア光誘発核分裂の微視的解析
著者: K. Uzawa (日本原子力研究開発機構 原子核データセンター)

---

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

核分裂は、原子核が大振幅変形を起こし、集団運動と単粒子運動が相互作用してエネルギー散逸を伴う複雑な量子多体過程である。この過程を微視的な自由度に基づいて理解することは核物理学における重要な課題である。

既存の手法には以下のような限界があった：

• 時間依存密度汎関数理論 (TDDFT): 大きな変形や 1 体散逸を記述できるが、核分裂障壁を越える「量子トンネル効果」を記述できないため、障壁貫通ダイナミクスや分裂確率の計算には適用できない。
• 時間依存生成座標法 (TDGCM): 障壁貫通を記述可能だが、初期状態の準備（特に反応入口チャネルからの波動パケットとしての記述）に恣意性があり、断面積の計算や入口チャネル効果の解析において一貫した扱いが難しい。

これに対し、非平衡グリーン関数法 (NEGF) は、電子輸送物理学で確立された手法であり、複合核反応理論に基づいて入口チャネルから出口チャネルへの遷移を自然に記述できるため、核分裂断面積の計算に適している。しかし、光誘発核分裂（特にサブバリア領域）への適用と、その微視的メカニズムの解明は十分に行われていなかった。

2. 手法と理論的枠組み (Methodology)

本研究では、NEGF 法を 236U(γ, f) 反応に適用し、サブバリア領域（5 MeV ≦ Eγ ≦ 6 MeV）での核分裂を微視的に解析した。

• モデル空間の構築:

  • 分裂経路に沿った Skyrme ハートリー・フォック (SHF) 波動関数（変形パラメータ $Q$ と粒子 - 空孔励起 $E_\mu$ でラベル付け）を基底状態として重ね合わせ、多体基底を構成した。
  • 変形パラメータ空間は $Q_{20}$-$Q_{30}$ 平面の 2 次元ポテンシャルエネルギー曲面を用い、19 の基準状態と粒子 - 空孔励起状態（全 145,124 状態）を含む大規模なモデル空間を構築した。
  • 対相関は残存相互作用として後から取り入れた。

• ハミルトニアンと残存相互作用:

  • 対角要素は変形ポテンシャルと励起エネルギー。
  • 非対角要素には、以下の 3 種類の残存相互作用を取り入れた：
    1. 単極型対相互作用: 0+ 状態の励起エネルギーを再現するよう調整。
    2. ランダム残存相互作用: 散逸的な形状進化を記述するため、正規分布からサンプリングしたランダム数を用いた（特に陽子 - 中性子チャネルで重要）。
    3. Diabatic 相互作用: 2 次項に基づく結合項。

• 減衰幅行列 (Decay Width Matrices):

  • 光吸収/放出幅: 実験値に基づいた E1 ギャムマ強度関数とレベル密度を用いて決定。Brink-Axel 仮説を適用。
  • 核分裂幅: 経験的な係数を用いず、モデル空間内で収束する値（100 keV）を採用。
  • 非エルミットな Hill-Wheeler 方程式を解き、複素固有エネルギーと固有状態を得た。

• 断面積の計算:

  • 遅延グリーン関数 $G(E)$ を用いて、光吸収チャネルから核分裂チャネルへの透過係数 $T_{\gamma in, fis}$ を計算。
  • 大規模な線形方程式を反復解く計算コストを回避するため、シフト・インバート・アーノルディ法（Shift-invert Arnoldi method）を用いてグリーン関数を評価した。

• 固有チャネル解析 (Eigenchannel Analysis):

  • 透過係数を特異値分解（SVD）し、各固有チャネルの寄与を分離。これにより、どのモードが核分裂を支配しているかを解析した。

3. 主要な結果 (Key Results)

• 核分裂断面積の再現性:

  • 計算された 236U(γ, f) の核分裂断面積は、実験データ（Caldwell et al., Yester et al., Soldatov et al.）と定性的に良好な一致を示した。
  • サブバリア領域 (5.0 - 5.2 MeV): 実験値の約 2 倍以内で一致し、障壁以下の抑制効果をよく再現。
  • 障壁近傍 (5.3 - 5.7 MeV): 若干過大評価されたが、全体的な減少傾向は再現。
  • 障壁以上 (Eγ ≥ 5.7 MeV): 実験値の 5 倍以内で一致。
  • 光吸収断面積も経験値と 1.2 倍以内で一致しており、光吸収後の核分裂確率の微視的計算が妥当であることを示唆。

• 固有チャネル解析によるメカニズムの解明:

  • 核分裂確率は、第 1 固有チャネル（最低次の透過モード）が支配的であることが判明した（$|t_1|^2 \gg |t_n|^2$ for $n \ge 2$）。
  • 第 2 以降のチャネルの寄与はエネルギーとともに増加するが、全エネルギー領域で第 1 チャネルが優勢。
  • この結果は、ボーア・ホイーラーの遷移状態モデル (Bohr-Wheeler transition-state picture) を微視的な観点から支持するものである。サブバリア領域では、量子トンネル効果を通じて最低次の遷移状態を通過することが核分裂の主要な経路であることを示している。

4. 貢献と意義 (Contributions & Significance)

• NEGF 法の核分裂への適用拡大:
  従来の中性子誘発核分裂（235U(n, f)）から、光誘発核分裂（236U(γ, f)）へ手法を拡張し、特に量子トンネル効果が支配的なサブバリア領域での適用可能性を実証した。
• 微視的観点からの遷移状態モデルの裏付け:
  従来の統計モデル（遷移状態モデル）が仮定する「特定の遷移状態を通過する」という描像が、非平衡グリーン関数法による微視的な波動関数の伝搬解析（固有チャネル分解）によって裏付けられた。
• 計算手法の革新:
  大規模な基底空間（10 万超）におけるグリーン関数の評価に、シフト・インバート・アーノルディ法を適用することで、高エネルギー分解能での断面積計算を現実的な計算コストで実現した。
• 将来への展望:
  本研究の枠組みは、実験が困難な中性子過剰核（r 過程に関与する核種）の核分裂断面積や、核分裂生成物の質量分布の予測に応用可能な基礎となる。今後の課題として、より現実的な残存相互作用の導入や、多次元変形空間への拡張が挙げられている。

結論

本論文は、非平衡グリーン関数法を用いて 236U の光誘発核分裂を微視的に記述し、サブバリア領域の実験データを再現することに成功した。特に、核分裂確率が単一の主要な固有チャネル（最低遷移状態）によって支配されているという発見は、核分裂の量子多体ダイナミクスに対する理解を深め、統計モデルの微視的基盤を強化する重要な成果である。