[論文レビュー] Transmutation operators: construction and applications

Vladislav V. Kravchenko, Sergii M. Torba+1 more|arXiv (Cornell University)|Aug 1, 2017

Numerical methods for differential equations参考文献 8被引用数 3

ひとこと要約

本稿では、変数係数を有するシュレーディンガー方程式の解に対する新しい正確な表現を、移行作用素を用いて提示する。この表現は、移行核のフーリエ・レジェンドルおよびエルミート級数展開に依存している。主な貢献は、周波数パラメータ ω に依存しない近似精度を実現することであり、これにより、変数係数方程式の多数の固有データを高精度に計算可能となる。

ABSTRACT

Recent results on the construction and applications of the transmutation (transformation) operators are discussed. Three new representations for solutions of the one-dimensional Schrödinger equation are considered. Due to the fact that they are obtained with the aid of the transmutation operator all the representations possess an important for practice feature. The accuracy of the approximate solution is independent of the real part of the spectral parameter. This makes the representations especially useful in problems requiring computation of large sets of eigendata with a nondeteriorating accuracy. Applications of the exact representations for the transmutation operators to partial differential equations are discussed as well. In particular, it is shown how the methods based on complete families of solutions can be extended onto equations with variable coefficients.

研究の動機と目的

• 変数ポテンシャルを有する1次元シュレーディンガー方程式の解に対して、正確かつ数値的に取り扱いやすい表現を開発すること。
• 直交多項式展開（レジェンドルおよびエルミート）を用いて、移行核 K(x,t) を構築し、計算効率を向上させること。
• 解の表現における近似誤差が、特に大きな値または複素数の ω に対しても、周波数パラメータ ω に依存しないことを保証すること。
• 完全な解の族に基づく手法（例えば、基本解法）を、変数係数を有する偏微分方程式へと拡張すること。
• 移行作用素を用いて、(Δ − q(x))u = 0 のような方程式に対する完全な解の族を体系的に生成するフレームワークを提供すること。

提案手法

• 移行核 K(x,t) をフーリエ・レジェンドル級数として表現する：K(x,t) = Σₙ₌₀^∞ [βₙ(x)/x] Pₙ(t/x) とし、係数 βₙ(x) を形式的べき乗 φₖ(x) およびレジェンドル多項式の係数から導出する。
• f'' − q(x)f = 0 を満たす解 f(x) に対して、f(0)=1、f'(0)=0 を満たし、f²(s) を用いた繰り返し積分を用いて、形式的べき乗 φₖ(x) を再帰的に構築する。
• 解の表現 u(ω,x) = e^{iωx} + Σₙ₌₀^∞ iⁿβₙ(x)jₙ(ωx) を導出する。ここで jₙ は球ベッセル関数であり、一様収束性と ω に依存しない誤差境界を保証する。
• 第二の表現として、エルミート多項式を用いる：K(x,t) を Hₙ(t)e^{-t²} に展開し、u(ω,x) = e^{iωx} + √π e^{-ω²/4} Σₙ₌₀^N cₙ(x)(iω)^n を得る。係数 cₙ(x) は K(x,t) と Hₙ(t) の積分から計算される。
• ω に依存しない誤差境界を確立する：|u(ω,x) − u_N(ω,x)| ≤ ε_N(x) × C とし、C は ω の虚部にのみ依存する。これにより、広い周波数範囲にわたり安定した精度が保証される。
• 調和関数（例：z^m の実部および虚部）および基本解（例：log|x−Z|）に移行作用素 T を適用し、(Δ − q(x))u = 0 の完全な解の族を生成する。

実験結果

リサーチクエスチョン

• RQ1 移行核 K(x,t) は、効率的かつ正確な数値計算を可能にする閉形式の級数展開として表現可能か？
• RQ2 提案されたシュレーディンガー方程式の解の級数表現は、周波数パラメータ ω が変化しても一様に精度を保つのか、特に大きな値または複素数の ω に対してもそうか？
• RQ3 移行法に基づく手法は、変数係数を有する偏微分方程式へと、基本解法をどの程度まで拡張可能か？
• RQ4 (Δ − q(x))u = 0 の完全な解の族を、移行作用素を用いて体系的に生成する方法は何か？
• RQ5 解の切断級数近似における誤差は、ω に依存せず境界づけられるか？これにより、大規模な固有データ計算における頑健性が保証されるか？

主な発見

• 移行核 K(x,t) は正確にフーリエ・レジェンドル級数として表現可能である：K(x,t) = Σₙ₌₀^∞ [βₙ(x)/x] Pₙ(t/x) であり、係数 βₙ(x) は形式的べき乗 φₖ(x) およびレジェンドル多項式の係数から計算可能である。
• シュレーディンガー方程式の解 u(ω,x) は、u(ω,x) = e^{iωx} + Σₙ₌₀^∞ iⁿβₙ(x)jₙ(ωx) として表現可能であり、jₙ は球ベッセル関数であり、x に関して一様収束する。
• 近似誤差 |u(ω,x) − u_N(ω,x)| は、実数の ω に対しては ε_N(x) × √(2x) で、虚部 |Im ω| ≤ C を満たす複素数の ω に対しては ε_N(x) × sinh(Cx)/C で境界づけられ、ω に依存しない精度が証明される。
• 第二のエルミート多項式を用いた表現では、u(ω,x) = e^{iωx} + √π e^{-ω²/4} Σₙ₌₀^N cₙ(x)(iω)^n が得られ、誤差は π^{1/4} e^{(Im ω)^2/2} ε_N(x) で境界づけられ、再び ω に依存しない。
• 調和多項式に移行作用素 T を適用することで、(Δ − q(x))u = 0 の完全な解の族が生成され、φₖ(x) および y に関する単項式を用いて明示的に表される。
• 本手法により、基本解法が変数係数PDEへと拡張可能である：T[log|x−Z|] は、レジェンドル関数 Qₙ(Z/x) を含む級数として表現可能であり、(Δ − q(x))u = 0 の基本解の構築が可能となる。