[論文レビュー] Sharp estimate of the mean exit time of a bounded domain in the zero white noise limit

Boris Nectoux|arXiv (Cornell University)|Oct 20, 2017

Markov Chains and Monte Carlo Methods参考文献 16被引用数 5

ひとこと要約

本稿は、領域内に非退化した最小値をもつポテンシャルと境界上で正の法線微分を持つ条件下で、ゼロノイズ極限（$\varepsilon \to 0$）における過減衰ランジュバン拡散の平均抜出自時刻に関する鋭い漸近公式の厳密な証明を提供する。主な結果として、出口時刻の正確な指数的推定が確立され、先行研究における形式的公式を裏付け、関連するフォッカー＝プランク作用素の最小固有値に対する鋭い漸近推定を可能にする。

ABSTRACT

We prove a sharp asymptotic formula for the mean exit time from a bounded domain $D\subset \mathbb R^d$ for the overdamped Langevin dynamics $$d X_t = - abla f(X_t) d t + \sqrt{2\ve} \ d B_t$$ when $\ve o 0$ and in the case when $D$ contains a unique non degenerate minimum of $f$ and $\pa_{\mbf n}f>0$ on $\pa D$. This formula was actually first derived in~\cite{matkowsky-schuss-77} using formal computations and we thus provide, in the reversible case, the first proof of it. As a direct consequence, we obtain when $\ve o 0$, a sharp asymptotic estimate of the smallest eigenvalue of the operator $$L_{\ve}=-\ve Δ+ abla f\cdot abla$$ associated with Dirichlet boundary conditions on $\pa D$. The approach does not require $f|_{\partial D}$ to be a Morse function. The proof is based on results from~\cite{Day2,Day4} and a formula for the mean exit time from $D$ introduced in~\cite{BEGK, BGK}.

研究の動機と目的

• 過減衰ランジュバン過程の有界領域からの平均抜出自時刻に関する鋭い漸近公式の数学的厳密な証明を提供すること。この公式は、以前は形式的計算でのみ得られていた。
• 同じ漸近的状況下で、ディリクレ境界条件を満たすフォッカー＝プランク作用素 $L_\varepsilon = -\varepsilon\Delta + \nabla f \cdot \nabla$ の最小固有値に対する鋭い推定を確立すること。
• 境界上に非退化した臨界点を含まない「退化」した境界設定に対しても、アイring-Kramers型の公式を適用可能にするために、境界上での非退化臨界点の必要性を排除すること。
• 特に $\partial_\mathbf{n}f > 0$ on $\partial D$ および $D$ 内に一意の非退化最小値があるという最小限の幾何的仮定の下で、既存の遷移的出口時刻に関する結果を統一的かつ一般化して証明すること。

提案手法

• 証明は [3,4] で得られた平均抜出自時刻の公式と、[7,6] における容量推定とスペクトル漸近に関する結果を組み合わせる。
• 境界付近での局所座標と変数変換を用いて、問題を法線方向の1次元ディリクレ問題に還元する。
• 容量定式化におけるテスト関数の構成と、Lax-Milgram定理の適用により、法線微分とポテンシャルを含む1次元楕円型方程式の弱解を得る。
• 適切な関数クラス上で定義されたディリクレ型エネルギー汎関数の最小化を用いて、容量 $\mathrm{cap}_C(D^c)$ の下界を導出する。
• 境界積分 $\int_{\partial D} \partial_\mathbf{n}f(\sigma) e^{-f(\sigma)/\varepsilon} d\sigma$ におけるラプラス法を用いて、$\varepsilon \to 0$ の極限における主要寄与を分析する。
• 最終的な推定は、容量の下界と、平均抜出自時刻と容量との既知の関係を組み合わせることで得られ、鋭い漸近公式が得られる。

実験結果

リサーチクエスチョン

• RQ1 過減衰ランジュバン過程の有界領域からの平均抜出自時刻の正確な漸近的挙動は、ノイズ強度 $\varepsilon$ が 0 に近づく際にどのように振る舞うか？
• RQ2 境界にポテンシャルの非退化臨界点を含まない場合に、形式的 Eyring-Kramers型の平均抜出自時刻公式を、数学的に厳密に証明することは可能か？
• RQ3 ディリクレ境界条件を満たすフォッカー＝プランク作用素 $L_\varepsilon = -\varepsilon\Delta + \nabla f \cdot \nabla$ の最小固有値 $\lambda_\varepsilon$ は、小ノイズ極限においてどのように振る舞い、平均抜出自時刻とどのように関連するか？
• RQ4 境界上での法線微分 $\partial_\mathbf{n}f$ が出口時刻に果たす役割は何か？また、Morse型の仮定を必要としない状況で、鋭い推定を得るためにどのように利用できるか？

主な発見

• 本稿では、平均抜出自時刻の鋭い漸近公式を確立した：$\mathbb{E}_x[\tau_{D^c}] = \frac{(2\pi\varepsilon)^{d/2}}{\sqrt{\det \mathrm{Hess}\, f(x_0)} \int_{\partial D} \partial_\mathbf{n}f(\sigma) e^{-f(\sigma)/\varepsilon} d\sigma} e^{-f(x_0)/\varepsilon} (1 + O(\varepsilon))$ として $\varepsilon \to 0$ のとき、$D$ のコンパクト部分集合上のすべての $x$ に対して一様に成立する。
• ディリクレ境界条件を満たす作用素 $L_\varepsilon = -\varepsilon\Delta + \nabla f \cdot \nabla$ の最小固有値 $\lambda_\varepsilon$ は、$\varepsilon \to 0$ のとき $\lambda_\varepsilon \sim \frac{\int_{\partial D} \partial_\mathbf{n}f(\sigma) e^{-f(\sigma)/\varepsilon} d\sigma}{(2\pi\varepsilon)^{d/2} \sqrt{\det \mathrm{Hess}\, f(x_0)}} e^{f(x_0)/\varepsilon} (1 + O(\varepsilon))$ を満たす。
• 容量 $\mathrm{cap}_C(D^c)$ に対する下界が得られた：$\mathrm{cap}_C(D^c) \geq \int_{\partial D} \partial_\mathbf{n}f(\sigma) e^{-f(\sigma)/\varepsilon} d\sigma (1 + O(\varepsilon))$。これは主結果の証明に不可欠である。
• 証明は $f|_{\partial D}$ がMorse関数である必要がないため、Eyring-Kramers公式の適用範囲が退化境界幾何にまで拡張される。
• 余剰項 $O(\varepsilon)$ の完全な漸近展開を避けるが、そのような展開が存在することは予想される。