Angular and Kinetic Properties of Scission Neutrons within Time-dependent Density Functional Theory

本研究は、時間依存密度汎関数理論を用いて、核分裂片由来の中性子が$^{235}\mathrm{U}$、$^{239}\mathrm{Pu}$、および$^{252}\mathrm{Cf}$の核分裂におけるプロンプト核分裂中性子スペクトルの重要な高エネルギー成分であることを示し、それらの包含が従来の蒸発のみを考慮したモデルで見られる高エネルギー収率の系統的な過小評価を解決することを直接的な証拠として提示するものである。

要約

巨大で不安定な風船（原子核）が、突然二つに弾け飛ぶ場面を想像してみてください。これが核分裂です。数十年来、科学者たちは、この現象が起きると、結果として生じる二つの破片（断片と呼ばれます）が猛烈な速度で飛び散り、中性子と呼ばれる小さな粒子を吐き出すことを知っていました。

長い間、物理学者たちは、これらの中性子のすべては、熱いコーヒーから立ち昇る湯気のように、後になってから「蒸発」したものだと考えてきました。彼らは、断片が完全に形成され、安定して動いている状態で、中性子を放出し始めるのだと想定していました。

しかし、この新しい論文は、これらの中性子のいくつかは、まさに風船が弾ける瞬間に「蹴り出されている」ものであることを示唆しています。これらは**切断中中性子（scission neutrons）**と呼ばれます。これらは、分裂の後の穏やかに冷却される断片からではなく、分裂の瞬間の激しい混沌の中から生まれるのです。

研究者たちが、これら「スナップタイム（破裂の瞬間）」の中性子の存在を証明した方法は、以下のように簡単に説明できます。

1. スーパーコンピュータによるシミュレーション

分裂の間に何が起きているのかを見るために、科学者たちは顕微鏡を使ったのではなく、**時間依存密度汎関数理論（TDDFT）**という理論を用いて、スーパーコンピュータでその出来事の「映画」を再生しました。

これは、原子がどのように踊り、壊れていくかをシミュレートする、高精度な3Dビデオゲームのようなものです。以前のバージョンのこの「ゲーム」では、仮想世界があまりに小さすぎました。中性子が画面の端に当たってしまうため、科学者たちは中性子がどれくらいの速さで、どの方向に飛んでいるのかを正確に把握することができませんでした。

今回の研究では、彼らははるかに大きな仮想世界（以前よりも約3倍大きいもの）を構築しました。これにより、中性子が十分に飛び出し、落ち着くためのスペースが確保され、シミュレーションの「壁」によってデータが乱されることなく、正確に測定できるようになりました。

2. 「速度制限」の発見

視界をクリアにした後、彼らは特定の角度（主に分裂に対して横方向、あるいはやや後方への角度）に飛んでいく中性子を観察しました。すると、驚くべきことが分かりました。

• 「進入禁止」ゾーン： 低エネルギー（約150万から200万電子ボルト未満）の切断中中性子は存在しません。まるで速度制限があるかのようです。低速なものは「スナップタイム」の中性子にはなれないのです。
• 高速の集団： その代わりに、これらの中性子はすべて高速です。これらの中性子は特定の高い速度（3〜3.5 MeV）付近に集中し、さらに高速な粒子へと続く長い裾野を描いています。

これは、ダイビングボードから飛び降りる人々の群れのようなものです。「蒸発した」中性子は、後でプールサイドからカジュアルに歩いて出てくる人々です。一方、「切断面（スナップタイム）」の中性子は、ボードが壊れたまさにその瞬間に、激しく弾き飛ばされる人々です。ボードから弾き飛ばされた人々は、常に速く動いています。その特定のイベントから、ゆっくりとした動きの人は決して現れません。

3. 「失われたエネルギー」の謎の解決

科学者たちは長年、コンピュータモデルを現実世界の実験結果と一致させようと苦心してきました。そこには一つの問題がありました。

• 旧モデル： もし「蒸気（蒸発した中性子）」だけをカウントすると、コンピュータモデルが予測する高エネルギー中性子の数は少なくなってしまいます。これは、小さなカップでバケツを満たそうとしているようなもので、バケツが必要とする水の量に対して、カップが提供できる量が足りない状態です。
• 新モデル： 研究者が、先ほどのシミュレーションで見つけた「スナップタイム」の中性子を「蒸発した」中性子に加えたところ、計算がようやく成立しました。この組み合わせたモデルは、ウランやカリホルニウムの実験で測定された高エネルギーデータの数値と完璧に一致したのです。

4. なぜこれが重要なのか

これは大きなニュースです。なぜなら、純粋に微視的な理論（単に何かが存在すると仮定したり推測したりするのではなく）によって、これらの「スナップタイム」の中性子が予測され、その実在が証明された初めてのケースだからです。

• 以前： 科学者たちは、計算が合わないために、これらの中性子が存在すると「推測」するしかありませんでした。
• 今： コンピュータ・シミュレーションは、指示されることなく、自然にこれらの中性子を生み出しました。それは、天気予報に基づいて嵐を予測するのではなく、雲の動きを観察することで嵐を予測するようなものです。

まとめ

この論文は、原子が分裂するとき、中性子の総量のうち小さくも重要な一部（約6%から10%）が、破壊の激しい瞬間に生まれると結論付けています。これらの中性子は、特定の角度において常に高速であり、決して低速ではないという点で、他の中性子とは明確に区別されます。

データの中にこの「指紋」を見つけ出すことで、研究者たちはついに「スナップタイム」の中性子を「蒸発した」中性子から分離することに成功し、核分裂が実際にどのように機能しているのかについて、より鮮明で正確な姿を明らかにしました。これは、物質を繋ぎ止め、そして引き裂く根本的な力の理解を深める助けとなります。

技術概要

技術要約：時間依存密度汎関数理論における切断中性子（Scission Neutrons）の角度および運動学的特性

問題提起
核分裂の発見から約一世紀が経過した現在も、そのプロセスの完全な微視的記述は依然として困難であり、特に外側サドルポイントから切断（scission）に至る過程、およびその後の一次核分裂片（FFs）の形成に関する記述が課題となっている。現在の統計的蒸発モデルは、すべての即発中性子が完全に加速された核分裂片から等方的に放出されると仮定している。しかし、この枠組みでは、ネックの破断として知られる高度に非平衡な動的プロセス中に放出される中性子、すなわち切断中性子（SNs）を説明することができない。SNsは1939年以来、様々なモデルを通じて推測・研究されてきたが、即発中性子の観測量のみでは、切断時に放出された中性子と後に蒸発した中性子を区別できないため、実験的には未解決のままである。TDDFT（時間依存密度汎関数理論）を用いたこれまでの予測試みは、シミュレーション領域が小さかったために、中性子が計算境界に到達する前に角度および運動学的特性を同時に決定することができなかった。

手法
本研究では、$^{235}\text{U}(n_{th}, f)$、$^{239}\text{Pu}(n_{th}, f)$、および$^{252}\text{Cf}(sf)$におけるSN放出を調査するため、超流動系へと拡張されたTDDFTに基づく完全な微視的枠組みを採用している。

• シミュレーション設定: 著者らは、SkM* エネルギー密度汎関数（EDF）を用いた軸対称調和振動子ソルバー AxialHOHFB を使用した。重要な進展は、シミュレーション領域を大幅に拡大し（半軸 $Z_{max} = 70$ fm および $R_{max} = 55$ fm の楕円体）、これは従来の研究よりも約3.2倍大きいものである点である。
• 二段階の積分: 計算上の制約を管理しつつ精度を維持するために、二段階の時間積分手順が採用された。系を、より小さな基底（$B_1$）において静的な拘束ヘッセ・ボルディン（HFB）解から切断の瞬間まで発展させ、その後、中性子が境界に近づくまで発展を継続するために、より大きな基底（$B_2$）へと変換した。
• SN特性の抽出: 集団流エネルギー密度と固有運動エネルギー密度の比が閾値 $\eta = 3$ を超える領域 $V^{(SN)}_\eta$ を定義することで、SNを特定した。この基準により、外向きの集団流によって支配される中性子を分離している。
• 角度および運動学的解析: この領域内において、局所的な速度場および中性子あたりの集団流エネルギーを算出した。解析は、境界反射が発生する前にSNの収率が飽和する角度範囲 $110^\circ \leq \theta \leq 150^\circ$（軽核分裂片の方向に対する角度）に焦点を当てた。
• 実験との比較: 算出されたSNの運動エネルギースペクトルを、蒸発中性子（EN）成分に対するマクスウェル型蒸発モデルと組み合わせた。蒸発モデルのパラメータ（断片温度 $T$ および中性子多重度比 $\bar{\nu}_L/\bar{\nu}_H$）は、低エネルギー実験データ（$E \leq E_{thr}$）のみによって制約されている。この結合モデルを、$^{239}\text{Pu}$ および $^{252}\text{Cf}$ の高エネルギー即発中性子スペクトル（PFNS）データと比較した。

主な結果

• SNスペクトルの特性: 抽出された3つの系すべてのSNスペクトルは、明確な閾値挙動を示す。SNは特定のエネルギー閾値（$E_{thr}$）以下では存在せず、この閾値は約1.5–2 MeV（誘導分裂では約2 MeV、自発核分裂では約1.5 MeV）である。スペクトルは3–3.5 MeV付近でピークに達し、約18 MeVまで続く指数関数的な高エネルギー・テイルを特徴とする。
• 収率の推定: $^{239}\text{Pu}(n_{th}, f)$ について、全SN収率は $0.19 \pm 0.01$（全即発中性子の約6.6%）と算出された。$^{252}\text{Cf}(sf)$ については、収率は $0.39 \pm 0.02$（約10.4%）であった。
• モデルの妥当性: 算出されたSN成分を（低エネルギーデータのみで制約された）蒸発モデルに加えると、得られた全PFNSは、$^{239}\text{Pu}$ および $^{252}\text{Cf}$ の両方において測定された高エネルギー収率を正確に再現する。対照的に、蒸発のみのモデルは、95%信頼区間内であっても、高エネルギー収率を系統的に過小評価する。
• 角度の安定性: 本研究は、信頼できる漸近的特性が角度範囲 $110^\circ \leq \theta \leq 150^\circ$ 内でのみ抽出可能であることを確認している。この範囲外では、過渡的な干渉効果と境界反射により、現在のシミュレーション領域内で収率が飽和することはない。

意義および主張
本論文は、既存の測定された即発中性子スペクトルの中に、完全な微視的枠組みを用いて初めて切断中性子のシグネチャを特定したと主張している。TDDFTが事前の仮定なしに自然にSNを生成すること、およびこれらのSNが高エネルギー実験データと蒸発のみのモデルとの間の乖離を解消することを示すことで、本研究は、即発核分裂において無視できない切断中性子成分が存在することの直接的な証拠を提示している。

著者らは、本研究が標準的な統計的蒸発の図式を超えて、核分裂中性子スペクトルを解釈するための微視的な基礎を確立するものと考えている。さらに、特定の実験的シグネチャを特定した。それは、エネルギー閾値（$E_{thr}$）以下におけるSNとENの寄与の分離である。これは、指定された角度範囲における低エネルギー測定が蒸発成分に対する最も堅牢な制約を提供し、高エネルギー・テイルがSN寄与の主要な指標となることを意味している。本研究は、現在のシミュレーションが特定の角度範囲に限定されているものの、本手法が、一世紀近く理論的に推測されてきた切断中性子の測定可能なシグネチャを、見事に孤立させたものであると結論付けている。