[논문 리뷰] Sharp estimate of the mean exit time of a bounded domain in the zero white noise limit

Boris Nectoux|arXiv (Cornell University)|2017. 10. 20.

Markov Chains and Monte Carlo Methods참고 문헌 16인용 수 5

한 줄 요약

이 논문은 잠재력이 영역 내에 고유한 비퇴화 최소값을 가지며 경계에서 잠재력의 법선 도함수가 양수인 조건 하에, 영역 외부로의 평균 이탈 시간에 대한 날카로운 점근적 공식을 영소음 한계($\varepsilon \to 0$)에서 엄밀하게 증명한다. 주요 결과는 이탈 시간에 대한 정확한 지수적 추정을 수립하여 이전 연구에서 형식적으로 유도된 공식을 확인하고, 관련된 Fokker-Planck 연산자의 최소 고유값에 대한 날카로운 점근적 추정을 가능하게 한다.

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

연구 동기 및 목표

• 소음 강도 $\varepsilon \to 0$의 극한에서, 과다진동 랭지비안 과정이 유계 영역에서의 평균 이탈 시간에 대한 날카로운 점근적 공식에 대한 엄밀한 수학적 증명을 제공하는 것. 이 공식은 이전에 형식적 계산을 통해 유도된 바 있다.
• 동일한 점근적 극한에서, 딜리클레 경계 조건을 갖는 Fokker-Planck 연산자 $L_\varepsilon = -\varepsilon\Delta + \nabla f \cdot \nabla$의 최소 고유값에 대한 날카로운 추정을 수립하는 것.
• 경계에서 $f|_{\partial D}$가 모어스 함수가 아니더라도, Eyring-Kramers 유형의 공식을 고려할 수 있도록, 경계에서 비퇴화 임계점이 필요 없도록 하는 것.
• 특히 $\partial_\mathbf{n}f > 0$ on $\partial D$ 및 $D$ 내에 고유한 비퇴화 최소값이 존재하는 조건 하에, 최소한의 기하학적 가정으로 기존의 비정상적 이탈 시간 결과를 통합하고 일반화하는 것.

제안 방법

• 증명은 [3,4]에서 유도된 평균 이탈 시간 공식과 [7,6]에서의 용량 추정 및 스펙트럼 점근적 성질에 관한 결과를 결합한다.
• 경계 근처에서의 좌표 변환과 국소 좌표계를 사용하여 문제를 법선 방향의 일차원 딜리클레 문제로 축소한다.
• 용량 공식에서 시험 함수를 구성하고, 라크스-밀그램 정리를 적용하여 법선 도함수와 잠재력이 포함된 일차원 타원형 미분방정식의 약한 해를 도출한다.
• 적절한 함수 클래스 위에서 딜리클레 유형의 에너지 함수를 최소화하는 변분 원리를 사용하여, 용량 $\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 소음 강도가 작아질 때, Fokker-Planck 연산자 $L_\varepsilon = -\varepsilon\Delta + \nabla f \cdot \nabla$의 최소 고유값 $\lambda_\varepsilon$는 어떻게 행동하는가? 이는 평균 이탈 시간과 어떻게 연결될 수 있는가?
• RQ4 법선 도함수 $\partial_\mathbf{n}f$가 경계에서 이탈 시간에 미치는 역할은 무엇이며, 모어스 유형의 가정 없이도 날카로운 추정을 유도하는 데 어떻게 활용될 수 있는가?

주요 결과

• 논문은 평균 이탈 시간에 대한 날카로운 점근적 공식을 수립한다: $\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))$ as $\varepsilon \to 0$, $D$의 컴acts부분집합에 대해 균일하게 성립한다.
• 디리클레 경계 조건을 갖는 연산자 $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}$가 모어스 함수일 필요가 없으므로, Eyring-Kramers 공식의 적용 범위를 퇴화된 경계 기하학에까지 확장한다.
• 잔여항 $O(\varepsilon)$의 전체 점근 전개를 피하는 방법을 사용하지만, 그러한 전개가 존재할 것으로 예상된다.