[論文レビュー] Spectral decomposition of 1D Fokker - Planck differential operator

Igor A. Tanski|arXiv (Cornell University)|May 1, 2006

Fractional Differential Equations Solutions参考文献 4被引用数 4

ひとこと要約

本稿では、分離変数法とフーリエ変換を用いて1次元Fokker-Planck微分作用素の完全なスペクトル分解を提示する。固有関数はエルミート多項式で表され、固有値は λ = −αn − k(2πm/(αa))² として得られる。初期値問題（Cauchy問題）の基本解は固有関数展開により構成され、先行研究の結果と一致することが示され、周期的境界条件および全空間におけるFokker-Planck方程式に対して、きめ細かなスペクトル的枠組みを提供する。

ABSTRACT

We construct spectral decomposition of 1D Fokker - Planck differential operator. This reveal solution of Cauchy problem. We develop fundamental solution of Cauchy problem and compare it with one obtained by other means in our former work [5].

研究の動機と目的

• 周期的境界条件の下で1次元Fokker-Planck微分作用素の完全なスペクトル分解を導出すること。
• 分離変数法とフーリエ変換技術を用いて、Fokker-Planck作用素およびその随伴作用素の固有関数と固有値を構成すること。
• スペクトル展開を用いて初期値問題の基本解を導出し、先行研究の結果と整合性を確認すること。
• フーリエ積分の極限を用いて、周期的境界条件からのスペクトル解を全空間に拡張すること。

提案手法

• 固有値方程式 L(φ) = λφ に分離変数法を適用し、φ(x,v) = X(x)V(v) を仮定することで、x と v に関する常微分方程式が得られる。
• 周期的境界条件を満たす x に依存する常微分方程式を解き、Xm = exp(2πimx/a) を得て、離散的 m モード解が得られる。
• v に依存する常微分方程式をフーリエ領域の式に変換し、F(ω) についての1階常微分方程式を解き、F(ω) = exp(−k/(2α)(ω+ω*)²)(ω−ω*)ⁿ を得る。
• ガウス関数の変換、シフト定理、微分定理などのフーリエ変換の性質を用いて、F(ω) を逆変換し、固有関数をエルミート多項式で表現する。
• 変数変換と複素数シフトを施し、最終的な固有関数を平面波と複素引数をもつエルミート関数の積として表現する。
• 作用素およびその随伴作用素の固有関数を用いたスペクトル展開により基本解を構成し、a→∞ の極限をとることで全空間上の解が得られる。

実験結果

リサーチクエスチョン

• RQ1 周期的境界条件の下での1次元Fokker-Planck微分作用素の完全な固有関数と固有値は何か？
• RQ2 フーリエ解析と特殊関数を用いて、非自己随伴なFokker-Planck作用素のスペクトル分解をどのように構成できるか？
• RQ3 1次元Fokker-Planck方程式の初期値問題に対する基本解の明示的形は何か？ また、先行研究の結果とどのように比較できるか？
• RQ4 フーリエ積分の極限を用いて、周期的境界条件におけるスペクトル解が全空間にどのように拡張されるか？

主な発見

• Fokker-Planck作用素の固有値は、n が非負整数で m が整数モードインデックスであるとき、λ = −αn − k(2πm/(αa))² で与えられる。
• 固有関数は φₘₙ(x,v) = exp(2πim(x−v/α)/a) × exp(−αv²/(2k)) × Hₙ(√(α/(2k))(v + 4πimk/(α²a))) で表され、Hₙ はエルミート多項式を表す。
• 全空間における基本解は G∞(t,x,ξ,v,u) = (1/(4πk√D*)) × exp(−Q(ˆx, ˆv)/(4kD*)) として導出され、Q は変数 ˆx と ˆv の二次形式である。
• この全空間の解は、著者の先行研究 [5] における経路積分法および特性関数法で得られた結果と一致し、一貫性が確認された。
• 周期的境界条件ではスペクトルが離散的であり、全空間では連続的であり、すべての固有値は実数かつ非正である。
• ゼロ固有値は、一様な空間密度とゼロ平均速度をもつマクスウェル速度分布に対応し、平衡状態を表す。