[論文レビュー] Estimating Spectral Density Functions for Sturm-Liouville problems with two singular endpoints

Charles T. Fulton, D. B. Pearson+1 more|arXiv (Cornell University)|Mar 12, 2013

Spectral Theory in Mathematical Physics参考文献 18被引用数 5

ひとこと要約

本稿では、2つの特異端点をもつSturm-Liouville問題の固有密度関数の計算のための新しい数値的手法を開発した。これは、正則左問題に対して既存のアルゴリズムを拡張したものである。x=0における正則特異点の近傍でFrobenius級数を用い、数値解法に区分的三角双曲線スプラインを適用し、無限遠点における初期値問題としてのTitchmarsh-Weyl m関数の境界値問題定式化を組み合わせることで、径方向水素原子およびBessel方程式において機械精度に近い精度を達成した。SLEDGEソフトウェアに比べ、精度と速度の両面で顕著に優れた性能を示した。

ABSTRACT

In this paper we consider the Sturm-Liouville equation -y"+qy = lambda*y on the half line (0,infinity) under the assumptions that x=0 is a regular singular point and nonoscillatory for all real lambda, and that either (i) q is L_1 near x=infinity, or (ii) q' is L_1 near infinity with q(x) --> 0 as x --> infinity, so that there is absolutely continuous spectrum in (0,infinity). Characterizations of the spectral density function for this doubly singular problem, similar to those obtained in [12] and [13] (when the left endpoint is regular) are established; corresponding approximants from the two algorithms in [12] and [13] are then utilized, along with the Frobenius recurrence relations and piecewise trigonometric - hyperbolic splines, to generate numerical approximations to the spectral density function associated with the doubly singular problem on (0,infinity). In the case of the radial part of the separated hydrogen atom problem, the new algorithms are capable of achieving near machine precision accuracy over the range of lambda from 0.1 to 10000, accuracies which could not be achieved using the SLEDGE software package.

研究の動機と目的

• x=0が正則特異点で、x=∞が非振動的限界点である2つの特異端点をもつSturm-Liouville問題に対する固有密度関数 f(λ) = ρ′(λ) の特徴付けを拡張すること。
• f(λ) を、無限遠点における特定の漸近的条件を満たす一階常微分方程式系 (P,Q,R)ᵀ の解の有理関数として表現すること。
• 正則左問題に対して開発された数値アルゴリズム [12,13] を、x=0に正則特異点、x=∞に限界点をもつ二重特異問題に適応するため、有限の一致点 x₀(λ) における高精度初期値解を用いること。
• 大規模な λ 範囲（径方向水素原子およびBessel方程式を含む）において、高い精度と効率を達成できる新しい数値コード（AutoB）の開発および検証を行うこと。
• SLEDGEソフトウェアの制限、特に大規模な λ 値における精度と計算時間の点で、新手法がその限界を克服できることを示すこと。

提案手法

• 固有密度 f(λ) は、f(λ) = 1 / [π(Pϕ² + Qϕϕ′ + Rϕ′²)] として表現され、(P,Q,R)ᵀ は行列 A が λ と q(x) に依存する一階線形常微分方程式系 dU/dx = A(x,λ)U を満たし、無限遠点における初期条件をもつ。
• 各 λ に対して一意な解を定義するため、境界条件として limₓ→∞ (P,Q,R)ᵀ = (√λ, 0, 1/√λ)ᵀ を用いる。
• x=0 の近傍では、ϕ(x,λ) と ϕ′(x,λ) を正確なFrobenius級数展開により計算し、有限な x₀(λ) > 0 における初期条件の高精度を保証する。
• 区分的三角双曲線スプラインを用いたシューティング法により、Sturm-Liouville方程式および (P,Q,R)ᵀ 系を x₀(λ) から一致点 x まで統合し、f(λ) の高精度計算を実現する。
• Titchmarsh-Kodairaの公式を活用して固有密度と m 関数を関連づけ、数値スキームに既知の解析的性質を組み込む。
• 得られたアルゴリズムは AutoB コードとして実装され、径方向水素原子およびBessel方程式に対してテストされ、SLEDGE との比較を通じて検証および性能評価が行われた。

実験結果

リサーチクエスチョン

• RQ1 正則左問題と同様に、Titchmarsh-Weyl m 関数の境界値問題定式化を用いて、2つの特異端点をもつSturm-Liouville問題の固有密度関数を特徴付けられるか？
• RQ2 正則左問題用に設計された数値アルゴリズムを、x=0に正則特異点、x=∞に限界点をもつ二重特異問題に適応するにはどうすればよいか？
• RQ3 Frobenius級数、スプライン統合、漸近的境界条件を組み合わせた新規数値手法による f(λ) の達成可能な精度と計算効率はどの程度か？
• RQ4 特に量子力学的問題（径方向水素原子など）において、SLEDGE などの既存ソフトウェアに比べ、新手法が精度と速度の両面で優れているか？
• RQ5 長大な λ 範囲にわたり、固有密度関数に対してどの程度機械精度に近い精度を達成できるか？

主な発見

• AutoBコードは、径方向水素原子の固有密度関数に対して機械精度に近い精度を達成し、λ ∈ [0.1, 10000] の範囲で相対誤差が 10⁻⁸ 未満であった。SLEDGEに比べ顕著に優れた性能を示した。
• ν=1 のBessel方程式に対して、τ=10⁻⁴ 時に相対誤差 1.30×10⁻⁷、τ=10⁻¹⁰ 時に 1.73×10⁻¹³ を達成し、τ=10⁻¹⁰ 時の計算時間は3秒未満であった。
• 水素原子の場合、SLEDGEは λ > 1000 で 10⁻³ の精度を達成できず、λ=1000 の計算に177秒以上を要したが、AutoBは同じ計算を4秒未満で10⁻⁸の精度で実行した。
• 水素原子（ℓ=1）およびBessel方程式を含む複数のポテンシャルに対して、本手法は一貫した高精度を示し、全テスト範囲および許容誤差レベルで安定した性能を発揮した。
• x₀(λ) における初期条件にFrobenius級数を用い、区分的三角双曲線スプラインを用いた統合により、両端点に強い特異性をもつ問題に対しても高精度統合が可能となった。
• 径方向水素原子の固有密度関数は λ=10000 で 144274264.3 に達し、AutoBは正確な値と小数点以下8桁まで一致した。一方、SLEDGEは相対誤差が10⁻³を超えた。