MFGSB (ver. 1.0): Computer code for self-consistent mean-field calculations of atomic nuclei using Gaussian expansion method

千葉大学リポジトリに公開された「MFGSB」は、ガウス展開法を用いた原子核の自己無撞着平均場計算を行うためのコンピュータコードです。

要約

この論文は、**「原子核（原子の中心にある小さな世界）」をコンピューターでシミュレーションするための新しい道具箱「MFGSB」**の紹介書です。

まるで、複雑なパズルを解くための「魔法のルーレット」や「万能の設計図」のようなものです。専門用語を噛み砕いて、日常の風景に例えながら解説します。

1. この道具箱（MFGSB）は何をするもの？

原子核は、陽子と中性子という小さな粒がぎっしり詰まった、とても複雑な世界です。これらがどう動いているか、どう形を作っているかを理解するには、膨大な計算が必要です。

この「MFGSB」というプログラムは、**「原子核の形や動きを、自分自身で考えながら（自己無撞着に）計算する」**ためのツールです。

• 従来の方法： 丸いボールや箱のような単純な形を仮定して計算する。
• この道具箱のすごい点： 「ガウス関数（ガウスの鐘の形のような曲線）」を積み重ねる方法（ガウス展開法）を使います。これにより、どんなに複雑に歪んだ形や、不思議な動きをする原子核でも、高い精度で描き出すことができます。

2. 何ができるの？（料理に例えると）

この道具箱を使えば、原子核の「レシピ」をいくつか作ることができます。

• 0. 下ごしらえ（ポテンシャル計算）： 決まった形（リンゴや卵のような形）で、とりあえず計算してみる。
• 1. 基本の味付け（ハートリー・フォック計算）： 粒子同士がどう影響し合うかを計算し、基本の形を決める。
• 2. 仲間の絆（HF+BCS 計算）： 粒子同士が「ペア」になって仲良くする現象（超流動など）も考慮に入れる。
• 3. 究極の融合（HFB 計算）： 粒子がペアになる現象と、基本の形を同時に、最も精密に計算する。

さらに、原子核を「丸い球」のまま計算することも、「ひしゃげた楕円」や「ねじれた形」で計算することも可能です。まるで粘土細工のように、原子核の形を自由に変えてシミュレーションできるのです。

3. この道具箱の「魔法」はどこにあるの？

このプログラムには、他の道具にはない 3 つの「魔法」があります。

1. 万能な「味付け」ができる（相互作用）：
   原子核の中にある粒子同士の力（相互作用）を、どんな複雑なルール（湯川力やテンソル力など）でも扱えます。特に、他のプログラムでは難しい「粒子の向きによる力」まで正確に計算できるのが最大の特徴です。また、計算の邪魔になる「全体の動き（重心運動）」も、無理やり近似せず、正確に処理できます。

   • 例え話： 他の料理人が「塩とコショウ」しか使えないのに対し、この道具箱は「スパイスの王様」まで使いこなせる大料理人のようなものです。

2. エネルギーに合わせた「変身」ができる：
   粒子のエネルギーが変わると、その振る舞い（波の形）も変わります。この道具箱は、その変化を無駄なく、かつ正確に表現できます。

   • 例え話： 水が氷になったり蒸気になったりするように、粒子の状態に合わせて形を柔軟に変えられる「変身能力」を持っています。

3. 一度の「設計図」で全原子核をカバー：
   ここが最も画期的です。通常、原子核の種類（元素）や形が変わるたびに、計算の基礎となるパラメータ（設計図の寸法）を微調整する必要があります。しかし、この道具箱は**「一度設定した設計図」で、周期表にあるほぼ全ての原子核に適用できます。**

   • 例え話： 従来の方法は、家を作るたびに「基礎の深さ」や「壁の厚さ」を一つ一つ調整する必要がありましたが、この道具箱は**「万能の型紙」**を持っています。これを使えば、小さな家から巨大な城まで、型紙を調整せずに作れてしまうのです。

4. 使うにはどんな準備が必要？

• メモリ（記憶装置）： 原子核の形を細かく設定する（ℓcut=7 の場合）と、約 8GB のメモリが必要です。設定を細かくすればするほど、メモリは倍々で増えます（16GB、32GB...）。
• 必要な道具： 計算を助けるための「BLAS」や「LAPACK」という計算用ライブラリが入っている必要があります。
• データの大きさ： 計算に必要な「相互作用データ（レシピ本）」が非常に大きく、約 16.6GB あります。そのため、本体とデータを別々にダウンロードできる仕組みも用意されています。

まとめ

この論文は、「原子核という複雑な世界を、より正確に、より簡単に、そして広範囲にシミュレーションできる新しい強力なツール」の完成と公開を報告するものです。

研究者たちは、この「万能型紙（MFGSB）」を使うことで、これまで解けなかった原子核の謎を解き明かしたり、新しい元素の性質を予測したりできるようになります。まるで、宇宙の奥深くにある小さな世界の地図を、初めて詳細に描き出すためのコンパスを手に入れたようなものです。

技術概要

以下は、提示された論文（arXiv:2602.13356v1）に基づく、MFGSB コードに関する詳細な技術的サマリーです。

MFGSB (ver. 1.0) 技術的サマリー

1. 背景と課題 (Problem)

原子核の構造を記述する「自己無撞着平均場（SCMF）計算」において、従来の手法にはいくつかの限界が存在しました。

• 相互作用の制限: テンソル力を含むヤウカワ（Yukawa）相互作用など、多様な有効核力を用いた計算が困難でした。
• 基底関数の調整: 異なる原子核や変形に対して、基底関数のパラメータを個別に調整（チューニング）する必要があり、計算コストと手間がかかっていました。
• 漸近挙動の記述: 単粒子（s.p.）または準粒子（q.p.）波動関数のエネルギー依存性を持つ漸近挙動を、高精度かつ効率的に記述する手法が限られていました。
• 補正項の近似: 重心運動（c.m. motion）の 2 体項やクーロン交換項を扱う際に、追加の近似が必要になるケースがありました。

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

本論文で紹介されるコードMFGSBは、**ガウス展開法（Gaussian Expansion Method: GEM）**を採用した FORTRAN77 言語で記述された SCMF 計算コードです。

• ガウス展開法（GEM）の活用:
  • 基底関数としてガウス関数を用いることで、波動関数の空間的な広がりや漸近挙動を柔軟に記述します。
  • 基底関数のパラメータは原子核の種類や変形に依存しないため、核図表の広範囲にわたって同一のパラメータセットを適用可能です。
• 対称性の仮定:
  計算は以下の対称性仮定の下で実行されます。
  • A) 球対称（J）、パリティ保存（P）、時間反転対称（T）。
  • B) 軸対称（z 軸周りの回転対称 $J_z$）、P 保存、T 保存、y 軸に関する反射対称（R）。
  • C) $J_z$ 対称、P 保存、$RT$ 対称（R と T の積）。
• 数値解法:
  • 自己無撞着状態を達成するために、反復法（ハミルトニアンを毎回対角化）または共役勾配法（最終段階で対角化）を選択可能です。
  • HF+BCS 計算では、HF 部分と BCS 部分で異なる解法を選択できます。
• 制約条件:
  • 四重極モーメントの拘束項を B) および C) の対称性条件下で追加可能です。
  • 重心運動の処理はオプションで選択可能です。

3. 主要な貢献と特徴 (Key Contributions)

MFGSB コードは、GEM の特性を活かした以下の革新的な機能を提供します。

1. 多様な核力への対応:
   • テンソルチャネルを含むヤウカワ相互作用まで適用可能な、ユニークな SCMF コードです。
   • 2 核子有効相互作用の関数形式に柔軟に対応します。
2. 高精度な物理項の扱い:
   • クーロン交換項および重心運動の 2 体項を、追加の近似なしで正確に扱えます。
   • 単粒子・準粒子波動関数の「エネルギー依存性を持つ漸近挙動」を、適度な精度で効率的に記述します。
3. 汎用性の高い基底パラメータ:
   • 基底関数のパラメータが原子核種や変形に敏感でないため、パラメータの微調整が不要です。
   • 1 つの相互作用行列要素ファイルで、ほぼ全ての原子核をカバーできます。
4. 計算オプションの多様性:
   • 自己無撞着計算（HF, HF+BCS, HFB）に加え、Woods-Saxon ポテンシャルや調和振動子ポテンシャル下での計算（非自己無撞着）も可能です。

4. 結果と性能 (Results & Performance)

• メモリ要件:
  • 基底関数のカットオフ値 $\ell_{cut}$ に依存します。デフォルト値（$\ell_{cut} = 7$）では約 8 GB のメモリを消費します。
  • $\ell_{cut}$ が 1 増えるごとに、必要なメモリは約 2 倍になります。
• 依存ライブラリ:
  • 計算には BLAS および LAPACK ライブラリのインストールが必要です。
• データ容量:
  • 完全版（mfgsb.zip）は約 16.6 GB あり、その大部分は相互作用データ（int_data_l7 フォルダ）を占めています。
  • 必要に応じて、相互作用ファイルを別途ダウンロードする「軽量版（mfgsb_light.zip）」も提供されています。

5. 意義と重要性 (Significance)

MFGSB は、原子核物理における理論計算の枠組みを大きく広げる重要なツールです。

• 理論的厳密性の向上: テンソル力を含むヤウカワ相互作用や、重心運動の厳密な扱いなど、従来近似されていた物理効果を高精度に計算可能にしました。
• 計算効率と汎用性: 基底パラメータの調整が不要であるため、未知の原子核や極端な変形を持つ核種に対する系統的な研究が容易になります。これにより、核図表全体を網羅する大規模な計算が現実的なものとなります。
• 研究コミュニティへの貢献: Chiba University Repository から公開されており、研究者が自らの研究に容易にアクセス・利用できるようになっています。

本コードは、原子核の構造理解を深め、特に非対称核や変形核における微視的な記述を可能にするための強力な基盤技術として位置づけられます。