[论文解读] Repartition of the quasi-stationary distribution and first exit point density for a double-well potential

Dorian Le Peutrec, Boris Nectoux|arXiv (Cornell University)|Feb 17, 2019

Quantum chaos and dynamical systems参考文献 11被引用 9

一句话总结

本文研究了在温度 $ h \to 0^+ $ 时，一个在双势阱中具有退化势垒的过阻尼朗之万过程的拟平稳分布与首达点密度。结果表明，由于亚稳态性，拟平稳分布通常集中在单个势阱中，但在对称情况下可能在两个势阱间分裂，揭示了隧穿效应，该效应通过半经典分析技术对特征函数和边界通量的精确渐近展开得到证实。

ABSTRACT

Let f : R d $\ ightarrow$ R be a smooth function and (Xt) t$\\ge$0 be the stochastic process solution to the overdamped Langevin dynamics dXt = ----f (Xt)dt + $\\sqrt$ h dBt. Let $\\Omega$ $\\subset$ R d be a smooth bounded domain and assume that f | $\\Omega$ is a double-well potential with degenerate barriers. In this work, we study in the small temperature regime, i.e. when h $\ ightarrow$ 0 + , the asymptotic repartition of the quasi-stationary distribution of (Xt) t$\\ge$0 in $\\Omega$ within the two wells of f | $\\Omega$. We show that this distribution generically concentrates in precisely one well of f | $\\Omega$ when h $\ ightarrow$ 0 + but can nevertheless concentrate in both wells when f | $\\Omega$ admits sufficient symmetries. This phenomenon corresponds to the so-called tunneling effect in semiclassical analysis. We also investigate in this setting the asymptotic behaviour when h $\ ightarrow$ 0 + of the first exit point distribution from $\\Omega$ of (Xt) t$\\ge$0 when X0 is distributed according to the quasi-stationary distribution. 1 Setting and results 1.1 Quasi-stationary distribution and purpose of this work Let (X t) t$\\ge$0 be the stochastic process solution to the overdamped Langevin dynamics in R d : dX t = ----f (X t)dt + $\\sqrt$ h dB t , (1) where f : R d $\ ightarrow$ R is the potential (chosen C $\\infty$ in all this work), h > 0 is the temperature and (B t) t$\\ge$0 is a standard d-dimensional Brownian motion. Let $\\Omega$ be a C $\\infty$ bounded open and connected subset of R d and introduce $\ au$ $\\Omega$ = inf{t $\\ge$ 0 | X t / $\\in$ $\\Omega$} the first exit time from $\\Omega$. A quasi-stationary distribution for the process (1) on $\\Omega$ is a probability measure $\\mu$ h on $\\Omega$ such that, when X 0 $\\sim$ $\\mu$ h , it holds for any time t > 0 and any Borel set A $\\subset$ $\\Omega$, P(X t $\\in$ A | t < $\ au$ $\\Omega$) = $\\mu$ h (A).

研究动机与目标

• 研究当温度 $ h \to 0^+ $ 时，双势阱中两个势阱内拟平稳测度的渐近分布。
• 刻画当初始分布为拟平稳分布时，从有界区域 $ \Omega $ 中的首达点分布。
• 识别拟平稳分布集中在单个势阱还是两个势阱的条件，特别是在对称或近似对称势阱中。
• 利用半经典与威滕拉普拉斯算子技术，建立生成元主特征函数及其在边界上法向导数的精确渐近展开。
• 将双峰集中现象与量子和随机亚稳态中的隧穿效应联系起来。

提出的方法

• 使用过阻尼朗之万 SDE $ dX_t = -\nabla f(X_t)dt + \sqrt{h}\, dB_t $，其中 $ f $ 为有界区域 $ \Omega \subset \mathbb{R}^d $ 上的光滑双势阱势能函数。
• 将拟平稳分布 $ \nu_h $ 定义为生成元 $ L_{f,h}^{(0)} = -\frac{h}{2}\Delta + \nabla f \cdot \nabla $ 的狄利克雷实现的归一化基态。
• 对威滕拉普拉斯算子 $ \Delta_{f,h}^{(0)} = -h^2\Delta + |\nabla f|^2 - h\Delta f $ 应用半经典分析，重点关注低能级特征值与特征函数。
• 利用拉普拉斯方法与谐振子近似，对局部极小点与边界临界点附近的特征函数进行渐近展开。
• 推导特征函数梯度与边界层函数 $ \psi_{i_j} $ 的 $ L^2 $ 内积的精确估计，以捕捉通量贡献。
• 通过等距性论证与对称性假设表明，当 $ f $ 近似为偶函数时，拟平稳分布等分于两个势阱之间，且在展开式中 $ \mu(h) = 0 $。

实验结果

研究问题

• RQ1 在具有退化势垒的双势阱中，拟平稳分布何时集中在单个势阱，何时在两个势阱间分裂？
• RQ2 当初始分布为拟平稳测度时，首达点分布的渐近行为如何，随着 $ h \to 0^+ $？
• RQ3 对称性（例如存在交换两个势阱的等距变换）在拟平稳分布的双峰集中现象中起什么作用？
• RQ4 主特征函数及其在边界上的法向导数的渐近行为如何依赖于势能与区域的几何结构？
• RQ5 当势能对称或近似对称时，隧穿效应如何在拟平稳分布中体现？

主要发现

• 当 $ h \to 0^+ $ 时，拟平稳分布 $ \nu_h $ 通常集中在单个势阱中，每个势阱中的概率质量以 $ h^{-1} $ 的指数形式衰减，除非对称性阻止了这种集中。
• 当势能 $ f|_{\overline{\Omega}} $ 在等距变换 $ \Phi $ 下（近似）对称时，拟平稳分布等分于两个势阱之间，且在渐近展开中 $ \mu(h) = 0 $。
• 首达点密度集中在 $ \partial\Omega $ 上使势能最小的边界点 $ z_{p,\ell} $ 附近，渐近通量与 $ h^{\frac{d-6}{4}} e^{-H'/h} $ 成正比，其中 $ H' = \min_{\partial\Omega}f + H $。
• 在 $ \Sigma \subset \partial\Omega $ 上，特征函数的法向导数 $ \partial_n u_h $ 满足：若 $ \Sigma $ 不包含临界边界点，则 $ \int_\Sigma F \partial_n u_h e^{-2f/h} = \mathcal{O}(e^{-H'/h}) $；当 $ \Sigma $ 包含此类点时，其值与 $ z_{p,\ell} $ 处的通量匹配。
• 在对称条件下，特征函数 $ u_h $ 满足 $ \int_{\Omega} u_h e^{-2f/h} = (1 + \mathcal{O}(e^{-c/h})) \sqrt{2} \int_{\mathcal{O}_i} \widetilde{u}_i e^{-2f/h} $，从而证实了双峰集中。
• 特征函数展开中主项系数包含 $ \left( \det \text{Hess}\, f(x_i) \right)^{-1/4} $，表明极小点处曲率的影响。