Novel method for evaluating the eigenvalues of the Heun differential equation with an application to the Breit equation

技術要約：Heun微分方程式の固有値を評価するための新規手法、およびブライト方程式への応用

問題提起
静的なクーロンポテンシャルを介して二つのフェルミオン（例：電子と陽子）の相互作用を記述するブライト方程式は、4つの特異点を持つ二階微分方程式を導く。これは、この方程式が標準的な合流型超幾何微分方程式よりも特異であることを示しており、Heun微分方程式に分類される。従来の研究（Yamaguchi、Tursunovらによるものなど）では、変分法や級数解のマッチングを用いてブライト方程式を分析してきたが、ブライト方程式から導出される二階Heun方程式に直接基づいてエネルギー準位を決定した研究は現在存在しない。さらに、微細構造定数の二乗（$\alpha^2$）のべき乗におけるエネルギー準位の展開について、一次の項を超えた明示的な表現も存在しない。

手法
著者らは、スピン一重項（$S=0$）および特定のスピン三重項（$S=1$）状態に対してブライト方程式を解くための理論的枠組みを開発した。その手法は、以下の2つの異なる、かつ相互に検証可能なアプローチを通じて進行する。

1. 連分数（CF）アプローチ：

   • 径方向のブライト方程式を、合流型Heun方程式の標準形式へと変換する。
   • 固有関数を随伴ラゲール多項式を用いて展開する。
   • この展開により、展開係数に関する三項漸化式が得られる。
   • 漸化式を連分数へと変換する。固有値は、この連分数の零点を見つけることによって決定される。
   • 数値評価には、波動関数の正規化のためにミラーの逆向漸化アルデリズムを用い、固有値付近での連分数の収束の遅さを処理するためにピーター・ウィンの$\rho$-アルゴリズム（非線形数列加速法）を利用する。

2. グリーン関数および半無限行列式アプローチ：

   • 微分方程式を、随伴ラゲール演算子の固有関数から構成されるグリーン関数を用いた、第二種の斉次フレドホルム積分方程式へと再定式化する。
   • この積分方程式への解を随伴ラゲール多項式で展開することで、半無限行列式（$\det \Delta_{m,n} = 0$）の条件を導出する。
   • 著者らは、行列式の要素を簡略化するために、特定の特殊関数（$I_{m,n}, S^m_n, \hat{S}^m_n$）を導入する。
   • 数値的な安定性と収束性を確保するため、行列式の対角要素を修正する「収束因子」を導入する。
   • 行列式は、直接計算（ドッジソンの縮退法を使用）、前向き漸化、および相互検証のためのCF法という3つの方法を用いて数値的に解かれる。

主要な貢献

• Heun方程式の導出： 本論文は、スピン一重項の場合（$m=M$）の二階Heun微分方程式をブライト方程式から明示的に導出し、様々なスピン三重項状態に対するパラメータ化を一般化している。
• 新規な行列式法： 著者らは、グリーン関数積分方程式から導かれた半無限行列式に基づく、固有値を計算するための新しい手法を提示している。これは連分数法に対する代替手段となる。
• 高精度な数値結果： 本研究では、スピン一重項およびスピン三重項状態（$J=0$ および $J=L$）について、最初の5つの主量子数（$N=1$ から $5$）に対するエネルギー準位を算出している。得られた結果は25〜30桁の精度に達している。
• 整合性の確認： 2つの異なる数学的アプローチ（CFおよび行列式法）が25桁の精度で一致することが示されており、これにより数値的手法が検証されている。
• 既知の物理の再現： 本手法は、既存の文献（Scottら、Yamaguchiら）の結果を正確に再現し、Pattersonによって以前に議論された$\rho \sim 1/y$付近の「偽状態（spurious state）」の挙動を正しく特定している。
• ディラック極限： 変数の適切なマッピングを通じて、ディラック方程式の解が行列式の条件から回収できることが本フレームワークによって示されており、行列式の構造に対する洞察を提供している。

結果

• 固有値： 論文には、最初の5つの準位に対する固有値（$z_n$）およびエネルギー準位（$E_n$）の表付きの値が記載されている。スピン一重項の場合、結果は18〜25桁の精度で文献値と一致している。
• 波動関数： 径方向の波動関数（$F, K, G$）は、ミラーの逆向漸和アルゴリズムを用いて効率的に計算される。論文ではこれらの成分が可視化されており、期待される漸近挙動が確認されている。
• 収束解析： 級数展開の収束性が分析され、$\rho \ge 0$ における絶対収束が証明されている。また、物理的なパラメータ $y$（微細構造定数に関連）が大きい場合における行列式法の数値的な課題に対処しており、高精度な領域では前向き漸化が直接計算よりも優れていることを示している。
• 解析的構造： 行列式の分析は、固有値の級数展開に非解析的な項、具体的には微細構造定数の対数（$\ln(1/\alpha)$）が含まれる可能性を示唆している。これは既存の文献とは対照的であるが、ポジトロニウムに関する最近の高次計算とは一致している。

意義および主張
著者らは、自らの研究が精度を数桁向上させ、ブライト方程式を解くための堅牢な代替理論的枠組みを提供していると主張している。主な意義は以下の通りである：

1. 境界における級数解のマッチングに依存することなく、Heun方程式に基づいた固有値を直接計算する手法を提供すること。
2. 連分数法の厳密な整合性チェックとして機能する、半無限行列式の条件を提供すること。
3. 今後の研究において、$\alpha^2$ のべき乗（一次を超える）におけるエネルギー準位の明示的な展開を得ること、および、状態の順序問題（例：$2^1P_1$ と $2^1S_0$ の逆転）を解決するために、完全なブライト相互作用（光子交換）を含めるための理論的基礎を築くこと。

本論文は、解析的な展開に関しては控えめな姿勢を保っており、行列式の構造は $\ln(1/\alpha)$ 項の存在を示唆しているものの、本研究において $\alpha^2$ の明示的な展開がまだ成功裏に導出されていないことを認め、それを将来の研究目標としている。