[论文解读] 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

一句话总结

本文在零白噪声极限（$ε \to 0$）下，针对势能函数在区域内部具有唯一非退化极小值、且在边界上法向导数为正的条件下，严格证明了过阻尼朗之万扩散过程从有界区域中退出的平均首次 hitting 时间的精确渐近公式。关键结果建立了退出时间的精确指数估计，验证了先前工作中形式推导的公式，并为相关福克-普朗克算子的最小特征值提供了精确渐近估计。

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}.

研究动机与目标

• 为过阻尼朗之万过程在噪声强度 $ε \to 0$ 极限下从有界区域退出的平均首次 hitting 时间的精确渐近公式提供严格的数学证明，该公式此前仅通过形式计算得出。
• 在相同渐近极限下，建立福克-普朗克算子 $L_\varepsilon = -\varepsilon\Delta + \nabla f \cdot \nabla$ 在 Dirichlet 边界条件下的最小特征值的精确估计。
• 通过消除对边界上非退化临界点的依赖，将 Eyring-Kramers 型公式的适用范围扩展至退化边界情形，即当 $f|_{\partial D}$ 不是 Morse 函数时。
• 通过在最小几何假设下统一并推广现有亚稳态退出时间结果，特别是当 $\partial_\mathbf{n}f > 0$ 在 $\partial D$ 上成立且 $D$ 内部存在唯一非退化极小值时。

提出的方法

• 证明依赖于文献 [3,4] 中推导的平均首次 hitting 时间公式，结合文献 [7,6] 关于容量估计与谱渐近的成果。
• 通过在边界附近引入变量变换与局部坐标系，将问题约化为法向方向的一维 Dirichlet 问题。
• 该方法涉及在容量公式中构造一个测试函数，并应用 Lax-Milgram 定理，推导出包含法向导数与势能的一维椭圆方程的弱解。
• 利用变分原理，通过对某一类合适函数最小化狄利克雷型能量泛函，推导出容量 $\mathrm{cap}_C(D^c)$ 的下界。
• 通过在边界积分 $\int_{\partial D} \partial_\mathbf{n}f(\sigma) e^{-f(\sigma)/\varepsilon} d\sigma$ 上应用 Laplace 方法，分析其渐近行为，该积分在 $\varepsilon \to 0$ 时主导了整体贡献。
• 最终通过结合容量的下界估计与平均首次 hitting 时间和容量之间的已知关系，得到精确渐近公式。

实验结果

研究问题

• RQ1 当噪声强度 $\varepsilon$ 趋近于零时，过阻尼朗之万过程从有界区域退出的平均首次 hitting 时间的精确渐近行为是什么？
• RQ2 当区域边界不包含势能函数的非退化临界点时，能否严格证明此前形式推导的 Eyring-Kramers 型平均首次 hitting 时间公式？
• RQ3 在小噪声极限下，福克-普朗克算子 $L_\varepsilon = -\varepsilon\Delta + \nabla f \cdot \nabla$ 的最小特征值如何变化？其与平均首次 hitting 时间之间是否存在关联？
• RQ4 边界上的法向导数 $\partial_\mathbf{n}f$ 在决定退出时间中起什么作用？如何在不依赖 Morse 型假设的前提下，利用其推导出精确估计？

主要发现

• 本文建立了平均首次 hitting 时间的精确渐近公式：$\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$ 在 Dirichlet 边界条件下的最小特征值 $\lambda_\varepsilon$ 满足：$\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))$，当 $\varepsilon \to 0$ 时成立。
• 导出了容量 $\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)$ 的完整渐近展开，尽管预计此类展开存在。