[論文レビュー] Langevin Approach to Levy Flights in Fixed Potentials: Exact Results for Stationary Probability Distributions

A. A. Dubkov, Bernardo Spagnolo|ArXiv.org|Oct 5, 2008

Complex Systems and Time Series Analysis参考文献 2被引用数 14

ひとこと要約

本稿では、コーシー安定レヴィノイズを伴うランジュバン方程式から、関数的アプローチを用いて分数階のフォッカー・プランク方程式を導出し、対称的単阱ポテンシャル内での定常確率分布（SPD）の正確な解析的計算を可能にする。主な結果として、$ m \geq 2 $ のポテンシャルプロファイルに対して、SPDは有限分散を有する二峰性を示し、超拡산的レヴィ飛行にもかかわらず粒子の閉じ込めが生じることを示している。これはノイズの強度とポテンシャルの傾きの間の相互作用によるものである。

ABSTRACT

The functional method to derive the fractional Fokker-Planck equation for probability distribution from the Langevin equation with Levy stable noise is proposed. For the Cauchy stable noise we obtain the exact stationary probability density function of Levy flights in different smooth potential profiles. We find confinement of the particle in the superdiffusion motion with a bimodal stationary distribution for all the anharmonic symmetric monostable potentials investigated. The stationary probability density functions show power-law tails, which ensure finiteness of the variance. By reviewing recent results on these statistical characteristics, the peculiarities of Levy flights in comparison with ordinary Brownian motion are discussed.

研究の動機と目的

• 関数解析を用いて、対称的 $\alpha$-安定レヴィノイズを伴うランジュバル方程式から分数階フォッカー・プランク方程式を導出すること。
• コーシーノイズ下での対称的単阱ポテンシャル内における定常確率分布（SPD）の正確な解析的表現を求めること。
• 超拡散的レヴィ飛行が閉じ込められ、有限分散を示す条件を調査すること。
• ポテンシャルの傾きとノイズ強度が、SPDの二峰性構造に与える役割を明確にすること。

提案手法

• 無限に可除な分布および特性関数に基づく関数的アプローチを用い、対称的 $\alpha$-安定レヴィノイズを伴うランジュバル方程式から、分数階フォッカー・プランク方程式を導出する。
• このアプローチにより、任意のレヴィノイズに対して有効な空間分数階微分を含む一般化されたフォッカー・プランク方程式に確率的力学を写像する。
• コーシーノイズ（$ \alpha = 1 $）の場合、ノイズ強度 $ D $ とポテンシャルの傾き $ \gamma $ の比 $ \beta $ を含む有理関数として、SPDの正確な閉形式表現が得られる。
• レヴィ安定過程の特性関数を活用し、逆フーリエ変換を適用して実空間におけるSPDを導出する。
• 有限な二階モーメントを保証するため、$ U(x) = \gamma x^{2m}/(2m) $ という形のポテンシャル $ m \geq 2 $ を対象とする。
• 得られたSPDは有理関数の積として表現され、$ m = 3, 4, 5 $ に対して明示的な形が導出され、奇数および偶数の $ m $ に対して一般化される。

実験結果

リサーチクエスチョン

• RQ1 固定ポテンシャル内での超拡散的レヴィ運動が、有限分散を有する閉じ込められた定常状態を示す条件は何か？
• RQ2 レヴィ駆動系において、ポテンシャルの傾きとノイズ強度を変化させたとき、定常確率分布（SPD）の形状はどのように変化するか？
• RQ3 コーシーノイズを伴う対称的単阱ポテンシャル内でのレヴィ飛行の定常分布の解析的形は何か？
• RQ4 なぜレヴィ飛行では二峰性のSPDが出現するのか？また、ブラウン運動における単峰性のボルツマン分布とはどのように異なるか？
• RQ5 関数的アプローチが、レヴィノイズを伴うランジュバル方程式から分数階フォッカー・プランク方程式を正確に導出する仕組みは何か？

主な発見

• ポテンシャルプロファイルが $ m \geq 2 $ の場合、粒子座標の二階モーメントが有限であり、超拡散的レヴィ飛行にもかかわらず閉じ込められていることを示している。
• すべての調査されたポテンシャル（$ m \geq 3 $）に対して、定常確率分布は二峰性を示し、ノイズ強度とポテンシャルの傾きの制御パラメータ $ \beta $ にかかわらず成立する。
• $ \beta $ の増加に対しても二峰性構造が維持されることから、ノイズの増加やポテンシャルの傾きの低下が二峰性を破壊しないことが示された。
• $ m = 3, 4, 5 $ のSPDは明示的に有理関数として導出され、$ W_{st}(x) = \frac{\beta^5}{\pi(x^2 + \beta^2)(x^4 - 2\beta^2 x^2 \cos\pi/5 + \beta^4)} $ および同様の形の高次 $ m $ に対する式が得られた。
• 奇数 $ m = 2n+1 $ および偶数 $ m = 2n $ に対して一般化された解析的表現が得られ、SPDが $ \pi $ の有理倍角のコサイン項を含む有理関数の積として表現されることを示した。
• SPDのべき則的尾部は有限分散を保証し、ポテンシャルの傾きが増加するかノイズ強度が低下するにつれて分布はより狭くなる。