[論文レビュー] Eigenvalue problems, spectral parameter power series, and modern applications

Kira V. Khmelnytskaya, Vladislav V. Kravchenko+1 more|arXiv (Cornell University)|Dec 7, 2011

Quantum Mechanics and Non-Hermitian Physics参考文献 70被引用数 7

ひとこと要約

本稿では、変数係数を伴うスツルム＝リウヴィル方程式の固有値問題を解くための新規な解析的・数値的手法として、スペクトルパラメータのべき級数（SPPS）法を導入する。解をスペクトルパラメータλに関するべき級数として表現することにより、テイラー級数の形で明示的な特徴方程式が得られ、解析関数の零点を求める手法により、固有値を高精度かつ高効率に計算可能となる。数値例では14桁の精度が達成されたことが示された。

ABSTRACT

In the present review we deal with the recently introduced method of spectral parameter power series (SPPS) and show how its application leads to an explicit form of the characteristic equation for different eigenvalue problems involving Sturm-Liouville equations with variable coefficients. We consider Sturm-Liouville problems on finite intervals; problems with periodic potentials involving the construction of Hill's discriminant and Floquet-Bloch solutions; quantum-mechanical spectral and transmission problems as well as the eigenvalue problems for the Zakharov-Shabat system. In all these cases we obtain a characteristic equation of the problem which in fact reduces to finding zeros of an analytic function given by its Taylor series. We illustrate the application of the method with several numerical examples which show that at present the SPPS method is the easiest in the implementation, the most accurate and efficient. We emphasize that the SPPS method is not a purely numerical technique. It gives an analytical representation both for the solution and for the characteristic equation of the problem. This representation can be approximated by different numerical techniques and for practical purposes constitutes a powerful numerical method but most important it offers additional insight into the spectral and transmission problems.

研究の動機と目的

• 変数係数を伴う固有値問題に対する明示的な分散関係を求めるという長年の課題に取り組む。
• 特に境界条件に複素固有値やスペクトルパラメータを含む場合に限界を示す従来の数値的手法の制限を克服する。
• 理論的洞察と高精度な数値計算を両立できる統一的で解析的根拠を持つ手法を開発する。
• 周期的ポテンシャルやザハロフ＝シャバト系を含む多様な物理的問題に、SPPS法の有効性を示す。
• 有限差分法や射撃法が失敗する状況においても、計算効率が高く実装可能な代替手法を提供する。

提案手法

• SPPS法は、スツルム＝リウヴィル方程式の線形独立な2つの解を、スペクトルパラメータλに関するべき級数として構成する：u₁(x) = Σ aₖ(x)λᵏ および u₂(x) = Σ bₖ(x)λᵏ。
• 係数aₖとbₖは、p, rおよび方程式 (pv′)′ + qv = 0 の特解vを用いて再帰的に計算され、複雑な反復積分を回避する。
• 微分作用素のポリア型因数分解を活用することで、係数の体系的かつ安定的な生成が可能になる。
• 境界条件を適用することで特徴方程式が導かれる。その結果、固有値は関数κ(λ) = Σ aₙλⁿ の零点として得られる。
• 数値解は級数を多項式κₙ(λ)に切り詰め、標準的なソルバーを用いてその根を求める。
• MATLABなどの科学計算環境を用いた実装により、最小限のプログラミング作業で高精度な数値的安定性を達成できる。

実験結果

リサーチクエスチョン

• RQ1 SPPS法は、変数係数を伴うスツルム＝リウヴィル問題の特徴方程式に対して、明示的な解析的表現を提供できるか？
• RQ2 SPPS法は、有限差分法や射撃法といった従来の数値技術と比較して、精度および効率において優れているか？
• RQ3 SPPS法は、複素固有値や境界条件にスペクトルパラメータを含む固有値問題を効果的に扱えるか？
• RQ4 SPPS法は、周期的ポテンシャルやザハロフ＝シャバト系のような問題へ、どの程度一般化可能か？
• RQ5 ベンチマークケースにおいて、SPPS法は正確な解や基準解と比較して、どの程度の精度で固有値を計算できるか？

主な発見

• A=1の長方形ポテンシャル井戸に対して、m=2000およびN=120を用いたSPPS法は、Wolfram MathematicaのNSolveと12桁まで一致した。
• パラメータを増大させた（m=4000, N=180）結果、同じケースで14桁の一致精度が達成された。
• A=4の場合、SPPS法は3つの固有値を正確に計算し、特に高分解能下でNSolveの結果と同等またはそれを上回る精度を示した。
• ザハロフ＝シャバト系の固有値は、テイラー級数で定義された解析関数の零点を求める問題に還元することで、SPPS法により正しく計算された。
• ザハロフ＝シャバト問題の特徴方程式は、λに関するべき級数として明示的に導出され、係数は再帰的公式により計算可能であった。
• 特に非実固有値に対して複素平面での順序付けが欠落する場合に射撃法が失敗する状況において、SPPS法は優れた性能を示した。