[논문 리뷰] Repartition of the quasi-stationary distribution and first exit point density for a double-well potential

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

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 $ 에서의 첫 번째 탈출점 분포를 기술하는 것.
• 준정적분포가 한쪽 함몰지에 집중되는지, 또는 두 함몰지에 걸쳐 분포되는지의 조건을 규명하는 것 — 특히 대칭적이거나 거의 대칭적인 퍼텐셜에서의 경우.
• 주요 고유함수와 그 경계에서의 법선 도함수에 대한 정밀한 점근 전개를 반고전적 및 위튼 라플라스 연산자 기법을 사용하여 수립하는 것.
• 이중모드 집중 현상이 양자역학적 및 확률적 메타안정성에서의 터널링 효과와 어떻게 연결되는지 밝혀내는 것.

제안 방법

• 유한 도메인 $ \Omega \subset \mathbb{R}^d $ 에서 매끄러운 이중함몰 퍼텐셜 $ f $ 를 가진 과다진 랑주뱅 방정식 $ dX_t = -\nabla f(X_t)dt + \sqrt{h}\, dB_t $ 를 사용한다.
• 준정적분포 $ \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 퍼텐셜이 대칭적이거나 거의 대칭일 경우, 터널링 효과는 준정적분포에 어떻게 나타나는가?

주요 결과

• 준정적분포 $ \nu_h $ 는 일반적으로 $ h \to 0^+ $ 일 때 한 함몰지에 집중되며, 각 함몰지의 확률 질량은 $ 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} $ 를 포함하며, 이는 최솟값에서의 곡률의 영향을 보여준다.