Systematic study of superheavy nuclei within a microscopic collective Hamiltonian: Impact of quantum shape fluctuations

この論文は、PC-PK1 密度汎関数に基づく拘束三軸相対論的ハートリー・ボゴリューボフ計算から導出された微視的 5 次元集団ハミルトニアンを用いて超重元素の量子形状揺らぎを系統的に研究し、球形殻近傍での動的相関エネルギーが核の形状転移や崩壊エネルギーの予測に決定的な影響を与えることを示しています。

要約

1. 研究の舞台：「超巨大なお城」の建設現場

まず、原子核は「陽子」と「中性子」という小さなレンガが積み重なったお城のようなものです。
通常の元素（金や鉛など）は、このレンガがそこそこ積み上がっただけのお城ですが、この研究では**「陽子 104 個〜126 個」**という、とてつもなく巨大なレンガの山（超巨大原子核）に注目しています。

• これまでの常識（平均場近似）：
  これまでの研究では、「レンガを静かに積み上げれば、一番安定した形（お城の形）が決まる」と考えていました。まるで、静かな湖に石を落とすようなイメージです。
• 今回の新発見（量子形状揺らぎ）：
  しかし、この研究は**「実はレンガは常に微かに震えていて、形がゆらゆら揺れている」**という視点を取り入れました。
  想像してみてください。静かに積み上げたお城も、実は地震（量子の揺らぎ）で常にガタガタ震えていて、そのせいで「安定しているはずの形」が実は崩れやすくなったり、逆に「不安定なはずの形」が少しだけ安定したりするのです。

2. 使った道具：「5 次元の揺れ動く地図」

研究者たちは、この「ゆらゆら揺れるお城」を正確にシミュレーションするために、**「5 次元集団ハミルトニアン（5DCH）」**という高度な計算ツールを使いました。

• アナロジー：
  普通の地図は「経度と緯度（2 次元）」で場所を表しますが、このツールは**「お城の形（丸いか細長いか）、その傾き、そして揺れの大きさ」まで含めた「5 次元の地図」を描くことができます。
  これにより、単に「一番高い山（安定した状態）」を見つけるだけでなく、「山と山の間の谷が、実は揺れすぎていて、本当にそこに住めるかどうか」**まで判断できるようになりました。

3. 驚きの発見：「住めないお城」の存在

この研究で最も衝撃的だったのは、**「理論上は存在するはずのお城が、実は住めない（崩壊してしまう）」**という発見です。

• 発見の内容：
  従来の計算（静かな湖のイメージ）では、「中性子 184 個」や「258 個」のところで、魔法のように安定したお城が作れると予想されていました。
  しかし、**「揺れ（量子形状揺らぎ）」を考慮すると、「その場所にあるお城の壁（分裂の壁）が、揺れに耐えきれずに崩れてしまう」**ことがわかりました。
  • 結果： 中性子が 184 個や 258 個のあたりには、**「実は住めない（存在しない）お城」**が大量にあることが判明しました。研究者はこれを「束縛されていない（unbound）」状態と呼んでいます。
  • 比喩： 「一番高い山頂に家を作ろうとしたら、実はその山頂は常に激しく揺れていて、家が建つ前に転落してしまう場所だった」という感じです。

4. 形の変化：「細長い卵」から「丸い玉」へ

原子核の形は、中性子の数によって変わります。

• 軽い元素： 細長い「卵（プロレート）」の形。
• 中間の元素： 形が定まらず、ぐにゃぐにゃ揺れる「柔らかい状態」。
• 重い元素： 丸い「玉（球体）」の形。

この研究では、**「揺れを考慮すると、形の変化がもっと滑らかになる」こともわかりました。
従来の計算では「ここから急に丸くなる！」とギザギザした変化でしたが、揺れを考慮すると「なめらかに丸くなっていく」**ことが示されました。これは、実際の実験データと非常に良く一致しています。

5. 未来への示唆：「新しい元素」はどこにある？

• 悲観的な側面：
  中性子が 184 個や 258 個の「魔法の数字」の周りは、実は不安定で、新しい元素を作るのが非常に難しい場所かもしれません。
• 楽観的な側面：
  しかし、**「陽子 119 番や 120 番」という、人類がまだ作れていない新しい元素は、揺れを考慮しても「まだ住める（安定している）」ことがわかりました。
  これは、将来の研究所（ロシアのドゥブナや日本の理研など）にとって、「ここなら新しい元素を作れる可能性が高い！」**という重要な地図になります。

まとめ：この研究が教えてくれたこと

1. 原子核は「静か」ではない： 常に揺れており、その揺れが安定性を大きく変える。
2. 「魔法の数字」は崩れた： 昔から言われていた「中性子 184 個で安定」という説は、揺れを考えると少しズレている（実際は 182 個あたりが安定しやすい）。
3. 住めない場所の発見： 多くの超重元素は、実は「分裂して消えてしまう」運命にあるかもしれない。
4. 次の目標： 陽子 119 番・120 番の元素は、まだ作れるチャンスがある！

この論文は、「原子核というお城の建築計画」において、単なる設計図だけでなく、「地震（揺れ）」まで考慮した新しい建築基準を提案した画期的な研究だと言えます。

技術概要

この論文は、超重核（Superheavy Nuclei, SHN）の構造を研究するための微視的集団ハミルトニアン（5DCH）アプローチを用いた体系的な調査に関するものです。以下に、問題提起、手法、主要な貢献、結果、そして意義について詳細な技術的サマリーを日本語で記述します。

1. 問題提起 (Problem)

超重核領域（特に $Z \approx 118$ 以降）では、実験的なデータが限られており、核の安定性や構造を理解することは核物理学の最前線における重要な課題です。従来の平均場近似（Mean-Field Approximation）に基づく密度汎関数理論（DFT）やマクロ・ミクロモデルは、結合エネルギーや変形などの静的性質をある程度記述できますが、**量子形状揺らぎ（Quantum Shape Fluctuations: QSFs）**を自一貫的に取り込むことは困難でした。
QSFs は、核の基底状態エネルギーや励起スペクトル、特に核の安定性（自発核分裂半減期など）に大きな影響を与える動的相関エネルギーを伴います。特に、超重核領域における球形殻閉殻（$N=184, 258$）付近や、変形と球形の遷移領域において、QSFs の効果を体系的に評価し、それが核の「束縛状態（bound state）」の存在にどのような影響を与えるかという点は、未解決の重要な問題でした。

2. 手法 (Methodology)

本研究では、以下の手法を組み合わせて超重核を体系的に解析しました。

• 対象核: 陽子数 $104 \le Z \le 126$、中性子数$N \le 258$ の偶偶核（計 585 核）。
• 基礎理論: 制約付き三方対称相対論的ハートリー・ボゴリューボフ（TRHB）計算。
  • 密度汎関数：相対論的密度汎関数 PC-PK1 を使用。
  • 基底関数：3 次元調和振動子基底（18 主要殻まで）。
• 集団ハミルトニアン: 微視的 5 次元集団ハミルトニアン（5DCH）モデル。
  • 変形変数：四極子変形 $\beta$ と $\gamma$、およびオイラー角。
  • 計算プロセス：TRHB 計算から得られた集団ポテンシャルと質量パラメータ（慣性モーメントなど）を 5DCH に入力し、QSFs を含むシュレーディンガー方程式を解くことで、基底状態および励起状態のエネルギーと波動関数を求めました。
  • 零点振動エネルギー（ZPE）の補正を施し、集団ポテンシャルを再定義しました。
• 解析対象: 結合エネルギー、二核子分離エネルギー（$S_{2n}, S_{2p}$）、$\alpha$ 崩壊エネルギー（$Q_\alpha$）、基底状態の平均変形、$2^+_1$励起エネルギー、$B(E2)$ 遷移確率、および核分裂障壁の高さ。

3. 主要な貢献と結果 (Key Contributions and Results)

A. 結合エネルギーと予測精度の向上

• 平均場近似（TRHB）のみでは、実験値に対して結合エネルギーを過小評価する傾向があり（RMS 偏差 2.30 MeV）、特に超重核領域での精度が課題でした。
• 5DCH による QSFs の取り込みにより、RMS 偏差が 0.57 MeV まで劇的に改善されました（実験値 10 点のみでは 0.35 MeV）。
• このモデルは、同位体鎖全体にわたる核の結合エネルギーのアイソスピン依存性を効果的に捉え、原子番号 $Z$ が増加しても一貫した精度を維持することを示しました。

B. 形状進化と相転移

• 平均場予測: 軽い元素（$Z=104 \sim 118$）では $N=140 \sim 160$ や $N=188 \sim 220$ で大きな偏平（prolate）変形が見られ、$N \approx 178, 250$ で中程度の変形や偏平（oblate）へ遷移し、$N=184, 258$ で球形に近づくという複雑な形状進化が示されました。
• QSFs の効果: 5DCH による計算では、QSFs が形状の急激な変化を平滑化し、より漸進的な形状進化をもたらすことが確認されました。
• 形状転移: 偏平変形から中程度の変形・$\gamma$-軟らかい形状、そして球形への転移が、同位体鎖に沿って明確に観測されました。特に $Z \ge 120$ の同位体では、$N \approx 178$ 付近で偏平（oblate）変形が有利になることが示されました。

C. 「束縛状態」の欠如と核分裂障壁の重要性（最も重要な発見）

• 未束縛状態（Unbound States）の予測: 本研究の最も注目すべき発見は、$N \gtrsim 184$ および $N \gtrsim 240$ の遷移領域において、核分裂鞍点（fission saddle）に束縛された $0^+$ 基底状態が存在しないという結果です。
• メカニズム: これらの核では、平均場ポテンシャルの極小値が存在するものの、ポテンシャルの井戸が非常に浅く、QSFs による零点振動エネルギーが核分裂障壁の高さを上回ってしまいます。その結果、理論的に「束縛された」基底状態は存在せず、核は即座に核分裂するか、共鳴状態としてしか存在しない可能性があります。
• 対象核: 計算された 585 核のうち、175 核がこの「未束縛」領域に分類されました。
• 例外: 融合反応で合成が期待される $Z=119, 120$ ($N \approx 178-180$) の新しい元素は、QSFs を考慮しても依然として「束縛状態」であることが確認されました。

D. 殻閉殻のシフトと動的相関エネルギー

• 平均場計算では $N=184$ と $N=258$ で鋭く変化する二中性子分離エネルギー（$S_{2n}$）や $\alpha$ 崩壊エネルギー（$Q_\alpha$）が、5DCH 計算では$N=182$ と $N=256$ にシフトし、その急峻さが緩和されました。
• これは、球形殻付近での動的相関エネルギー（QSFs に起因）の急速な進化によるもので、実験データとの整合性を高めています。

E. 集団的観測量

• $E(2^+_1)$、$R_{4/2} = E(4^+_1)/E(2^+_1)$、$B(E2; 2^+_1 \to 0^+_1)$ の系統的な計算により、核の集団的性質（回転体、球形、遷移核）が明確に分類されました。
• $N \lesssim 170$ や $196 \lesssim N \lesssim 240$の核は回転体特性（$R_{4/2} \gtrsim 3.2$）を示し、$N \approx 184, 258$付近の核は球形特性（$R_{4/2} \lesssim 2.0$）を示すことが確認されました。

4. 意義 (Significance)

1. 理論的予測能力の確立: 相対論的 DFT と 5DCH を組み合わせたアプローチが、超重核領域の結合エネルギーや構造を高い精度で記述できることを実証し、未発見の超重核の探索指針を提供しました。
2. 核の安定性限界の再定義: 従来の平均場モデルが「安定な島」として予測していた領域の一部（特に $N \gtrsim 184$ の遷移領域）が、QSFs を考慮すると実際には「未束縛」であり、実験的に観測が極めて困難であることを示唆しました。これは、超重核の安定性限界を再評価する上で極めて重要です。
3. 実験への示唆: $Z=119, 120$ 付近の元素は合成可能である一方、より中性子豊富な領域（$N \gtrsim 184$）の核を合成・観測するには、多核子移動反応などの新しい手法が必要であり、その際 QSFs の影響を慎重に評価する必要があることを示しました。
4. 今後の展開: 本研究では八極子変形（octupole deformation）は完全には考慮されていませんが、八極子自由度が核分裂障壁をさらに低下させる可能性があるため、今後は 7 次元集団ハミルトニアン（7DCH）を用いた更なる研究が必要であると結論付けています。

総じて、この論文は、量子形状揺らぎが超重核の構造と安定性に決定的な影響を与えることを体系的に証明し、超重核研究の新たなパラダイムを提供した重要な業績です。