Variational free complement method with Gaussian-expanded complement functions: convergence with fixed Gaussian expansion length

本論文では、ガウス展開された補関数を用い、基底関数の数が無限大に近づく一方でガウス展開長（$n_\mathrm{G}$）が固定されている極限における、変分自由補完法のエネルギー収束について調査する。

要約

あなたは、水素原子（宇宙で最も単純な原子）の完璧な肖像画を描こうとしていると想像してください。このとき、あなたは**変分自由補完法（Variational Free Complement Method）**と呼ばれる特別なデジタル筆を使っています。この筆は、層を重ねるごとに詳細を加えていくことで、「真の」絵（原子の正確なエネルギー）に限りなく近づけるように設計されています。

この論文の中で、著者であるCong Wang氏は、**ガウス関数（Gaussian functions）**を使用する特定のバージョンのこの筆について検証しています。ガウス関数を「柔らかく、ふわふわとした雲」のような塗料だと考えてください。これらは数学的に扱いやすいのですが、特定の形状を持っています。つまり、滑らかで、急速に消えていく性質があります。

著者が行った核心的な実験を、簡単に説明します。

2つの実験

著者は、この「ふわふわした雲」の筆が、最終的に完璧な絵を描き上げることができるのかどうかを調べようとしました。ただし、固定された限られた数の雲の形（これを $n_G$ と呼びます）を使うという制約がある場合でも、です。彼はこう問いかけました。「もし、これらの特定の雲の層を永遠に増やし続けたとしたら、最終的に完璧なエネルギー値に到達できるのだろうか？」

彼は2つの異なるシナリオを実行しました。

シナリオ1：「単一の雲」の限界（固定された $n_G = 1$）

• 設定: 著者は基本的な「スレーター型（Slater-type）」の波からスタートし、補正のためにたった一つのガウス型の雲のみを使用して、それを改良しようと試みました。彼は、この同じ単一の雲の形状を何度も何度も、層として重ね続けました。
• 問題点: ガウス型の雲は「頑固」です。実際の原子に比べて、消えるのが早すぎます。もし雲の種類が一つしかなければ、原子の非常に「広がりのある（拡散した）」部分（diffuse parts）を正しく描くことは決してできません。
• 結果: 著者は1,200層まで計算を行いました。絵はどんどん良くなっていきましたが、途中で止まってしまいました。完璧なエネルギー（-0.5）には非常に近づいたものの、約-0.4998で停滞してしまったのです。それは、底に小さな穴が開いているカップを使ってバケツを満たそうとするようなものでした。何度注いでも、決して満タンにはなりません。
• 結論: 固定された少数の雲の形状を用いている場合、この手法は完璧な答えへと収束しません。どうしても超えられない「天井」に突き当たってしまうのです。

シナリオ2：「無限の雲」の限界（増加する $n_G$）

• 設定: 2番目の実験では、最初から「ガウス型」の初期波（最初から雲であるもの）を使い、雲の形（$n_G$）が無限に増えていくことを許容しました。
• 結果: 今回は、絵は完璧になりました。より多くの異なる種類の雲の形を加えていくことで、エネルギー値は真の答え（-0.5）に正確に収束しました。
• 結論: 雲の多様性を増やすことができれば、この手法は完璧に機能します。

大きな教訓

この論文は、ある特定の問いに答えています。「もし、固定された少数のガウス形状に縛られてしまった場合、ただ永遠に計算を続ければ、いつか方法は機能するのだろうか？」

答えは、**「いいえ」**です。

著者は、なぜそうなるのかを説明するために、**ミュンツ・サーツの定理（Müntz–Szász theorem）**という数学的概念（これは、ある形の集合が、あらゆる可能な曲線を作り出せるかどうかを判定するルールブックのようなものです）を用いています。彼は、固定された数のガウス形状に縛られているとき、原子の「拡散した（広がりを持った）」部分（遠くまで伸びている部分）が欠落していることを示しています。それらの特定の形状をどれほど積み重ねても、欠けているピースを作ることはできないのです。

これが意味すること、そして意味しないこと

• 意味すること: もし、この特定の手法を固定された少数のガウス関数を用いて使用する場合、どれほど計算能力を投入したとしても、数学的に完璧なエネルギーを得ることは決してできません。常にわずかな誤差が残ります。
• 意味しないこと: 著者は、この手法が役に立たないと言っているわけではありません。現実世界の化学において、科学者たちは通常、多くの異なるガウス形状（大きな $n_G$）と、妥当な数の層を使用します。そのような実用的なケースにおいて、この手法は非常によく機能し、高速です。この論文は、もしあなたが「雲」に対してあまりにケチになりすぎた場合（$n_G$ を固定して小さく保った場合）、その手法には決して超えられない限界があるということを警告しているのです。

要約すると: たった一種のレンガだけを使って完璧な家を建てることはできません。いくら積み重ねても無理なのです。隙間を埋めるためには、さまざまなサイズのレンガ（拡散関数）が必要です。この論文は、もしあなたがもっと多くの種類のレンガを使うことを拒むなら、あなたの家には常に、修復不可能な小さな隙間が残ることになる、ということを証明しています。

技術概要

技術要約：ガウス関数展開による補完関数を用いた変分自由補完法（Variational Free Complement Method）

問題提起
本研究は、ガウス関数展開を用いた変分自由補完（Free Complement, FC）法の収束特性を調査するものである。具体的には、自由補完法の次数（$n$）および補完関数の数（$M_n$）が無限大に近づく際、固定された有限のガウス展開長（$n_G$）という制約の下で、手法がハミルトニアンの厳密な固有エネルギーに収束するかどうかという問いに取り組む。

この研究は、他の計算手法や機械学習データセットのための参照値として、定常シュレディンガー方程式の数値的に厳密な解を得る必要性に端を発している。自由補完理論（最近「free complete-element」と改称された）は、こうした解を提供するための枠組みを提供するが、ガウス展開を用いた従来の実装では、積分の困難さや指数関数的な複雑さに直面してきた。これらは階層的デコンストラクション（hierarchical decontraction）によって緩和されているものの、理論的な疑問が残っている。すなわち、固定された $n_G$（STO-$n$G 展開におけるガウス関数の数）は、$n \to \infty$ の際に基底状態エネルギーへの収束に十分であるか、という点である。

手法
著者は、ベンチマークシステムとして非相対論的な電子水素基底状態を用い、以下の2つの異なるシナリオにおける収束性をテストする。

1. $n_G = 1$ で固定し、$n \to \infty$ とするケース：

   • 初期波動関数： スレーター型軌道 $\psi_0 = e^{-\zeta r}$ （ここで $\zeta = 0.5$）。
   • 補完関数： $g = 1 - e^{-\gamma r}$ （ここで $\gamma = 1 - \zeta$）という形式の $g$ 関数を用いて生成。
   • 展開： 補完関数は、単一のガウス関数（$n_G=1$）へと展開される。ガウス関数の指数は、STO-$n$G 基底関数セットから導出されたスケーリング関係に従い、数列 $\alpha(k) \propto (1 + k\gamma/\zeta)^2$ を形成する。
   • 理論的解析： 収束性は、$L^2(\mathbb{R}^3)$ および ソボレフ空間 $W^{(1)}_2(\mathbb{R}^3)$ における基底の完全性に関するミュンツ・シャス（Müntz–Szász）の定理およびその一般化を用いて解析される。完全性のための総和条件は、 $n_G=1$ の制約によって生成される指数の特定の数列に対して評価される。
   • 数値実装： 積分の評価とフィッティングには Python を、一般化固有値問題には Julia を用いて、レイリー・リッツ（Rayleigh-Ritz）変分計算を行う。$n$ の増加に伴うオーバーラップ行列の極端な悪条件（ill-conditioning）に対処するため、高精度算術（1200〜1500桁）を利用する。

2. $n_G \to \infty$ かつ $n \to \infty$ とするケース：

   • 初期波動関数： ガウス型軌道 $\psi_0 = e^{-\alpha r^2}$。
   • 補完関数： 同様の $g$ 関数の形式を用いるが、ガウス展開を行わない（実質的に $n_G \to \infty$）。これにより、関数は $e^{-k\gamma r}\psi_0$ という形式をとる。
   • 比較： このシナリオは、制御群（コントロール）として $n_G=1$ のケースと比較され、単一の初期ガウス関数によって指数の下限が制限されない場合の収束挙動を検証するために用いられる。

主要な結果

• ケース 1 ($n_G = 1$):

  • 理論的知見： 固定された $n_G=1$ の展開によって生成されるガウス指数の数列は、$L^2(\mathbb{R}^3)$ および $W^{(1)}_2(\mathbb{R}^3)$ における完全性のための十分条件を満たさない。完全性条件に関連する総和項は、無限大に発散するのではなく有限の値に収束しており、これは基底が不完全であることを示唆している。
  • 数値的知見： 変分計算の結果、エネルギーは約 $-0.4997995$a.u. へと緩やかに収束するが、これは厳密な基底状態エネルギーである$-0.5$a.u. よりも厳密に高い値である。外挿フィッティング（$E = A + B/M_n^C$）により、厳密な固有値に到達しない極限が確認された。
  • 結論： $n_G=1$ の場合、手法は $n \to \infty$ において厳密な固有エネルギーには収束しない。

• ケース 2 ($n_G \to \infty$):

  • 数値的知見： 初期波動関数がガウス型であり、任意の数の指数を許容する場合（$n_G \to \infty$ をシミュレート）、エネルギーは $n$ の増加に伴い厳密な値である $-0.5$ a.u. に収束する。
  • 結論： ケース 1 における収束の欠如は、ガウス指数の制限（具体的には、拡散関数（diffuse functions）の欠如と、指数の集積点の特異な振る舞い）に起因している。

意義と主張
本論文は、中心的な問いに対して明確な否定的な回答を提示している。すなわち、固定された $n_G$ の場合、ガウス展開された自由補完法は、$n$ が無限大に近づく際、厳密な固有エネルギーには収束しない。

著者は、これらの知見に関するいくつかの控えめながらも重要な主張を行っている：

1. 固定された $n_G$ の限界： ガウス展開された FC 法は、固定された $n_G$（本研究では $n_G=1$）を用いる場合、ミュンツ・シャスの完全性条件を満たす基底セットと比較して、水素原子の基底状態エネルギーを任意の精度で捉える計算能力を欠いている。
2. 拡散関数の役割： 収束の失敗は、拡散関数の不在と指数の集積点の特定の振る舞いに結びついている。ただし、拡散関数の欠如が必ずしもあらゆる文脈での非収束を意味するわけではなく、適切な指数の集積点（$n_G \to \infty$ のケースで見られるもの）によって補完可能であることも著者は指摘している。
3. 実用性： 理論的極限（固定された $n_G$ での $n \to \infty$）で見られる非収束は、有限の $n > 0$ および $n_G > 1$ を選択する実用的な計算におけるデメリットを意味するものではない。
4. 潜在的な救済策： 固定された $n_G$ のシナリオにおける収束問題は、ゼロ次を超えた指数の最適化（明示的に相関したガウス関数法と同様）によって解決できる可能性があるが、これは計算コストが高い。あるいは、ガウス型の $g$ 関数を使用するか、運動演算子（$p+k$ アプローチ）を導入することで、より効果的に完全性条件を満たす基底セットを生成できる可能性がある。

要約すると、本研究はガウス展開された自由補完法の理論的境界を明らかにし、有限の次数では高い精度を提供できるものの、固定されたガウス展開長は漸近的な厳密な固有エネルギーへの収束に対して根本的な限界を課すことを示している。