[논문 리뷰] Uniform minorization condition and convergence bounds for discretizations of kinetic Langevin dynamics

Alain Durmus, Aurélien Enfroy외 2명|arXiv (Cornell University)|2021. 07. 30.

Markov Chains and Monte Carlo Methods참고 문헌 53인용 수 4

한 줄 요약

이 논문은 넓은 범위의 시간 이산화된 운동학적 랑주뱅 동역학에 대해 균일한 미니멀라이제이션 조건과 기하학적 드리프트 조건을 수립하여, 총 변동과 $V$-노름에서 시간단계에 비례하여 선형적으로 스케일링되는 수렴 경계를 가능하게 한다. 주요 기여는 시간단계에 관계없이 균일한 기하학적 에르고딕성 결과를 도출한 것으로, 분자 동역학 및 베이지안 추론에서 신뢰할 수 있는 샘플링을 보장하고 포아송 방정식 해석에 필수적이다.

ABSTRACT

We study the convergence in total variation and $V$-norm of discretization schemes of the underdamped Langevin dynamics. Such algorithms are very popular and commonly used in molecular dynamics and computational statistics to approximatively sample from a target distribution of interest. We show first that, for a very large class of schemes, a minorization condition uniform in the stepsize holds. This class encompasses popular methods such as the Euler-Maruyama scheme and the schemes based on splitting strategies. Second, we provide mild conditions ensuring that the class of schemes that we consider satisfies a geometric Foster--Lyapunov drift condition, again uniform in the stepsize. This allows us to derive geometric convergence bounds, with a convergence rate scaling linearly with the stepsize. This kind of result is of prime interest to obtain estimates on norms of solutions to Poisson equations associated with a given numerical method.

연구 동기 및 목표

• 넓은 범위의 수치적 방법에 대해 시간 이산화된 운동학적 랑주뱅 동역학에 대해 균일한 미니멀라이제이션 조건을 수립하기 위해.
• 시간단계에 대해 균일한 기하학적 포스터-라이아프노프 드리프트 조건을 도출하여 안정성과 수렴성을 보장하기 위해.
• 시간단계 $\gamma$ 에 대해 선형적으로 스케일링되는 총 변동과 $V$-노름에서의 정량적 수렴 경계를 제공하기 위해.
• 통계적 추론과 분자 동역학에서 발생하는 포아송 방정식 해석 분석을 지원하기 위해 균일한 에르고딕성을 보장하기 위해.
• 기존 수렴 이론을 기울기 없는 드리프트와 비평형 설정을 포함한 과도한(운동학적) 랑주뱅 동역학으로 확장하기 위해.

제안 방법

• 에너지 유사 행동을 포착하기 위해 이차형 및 위치 에너지 항을 조합한 리아프노프 함수 $\overline{\mathpzc{W}}_{\gamma,\varpi}$ 를 도입한다.
• 모든 $\gamma \leq \bar{\gamma}_{U}$ 에 대해 시간단계에 무관하게 유효한, 커플링 방법과 초타원성(hypoellipticity)을 통한 균일한 미니멀라이제이션 조건을 수립한다.
• 기하학적 드리프트 조건 $R_{\gamma}\overline{\mathpzc{W}}_{\gamma,\varpi} \leq \lambda_{U}^{\gamma}\overline{\mathpzc{W}}_{\gamma,\varpi} + \gamma b_U$ 를 유도하며, 이는 $\lambda_U^\gamma < 1$ 를 만족하여 기하학적 에르고딕성을 보장한다.
• 전이 밀도를 근사하고 미니멀라이제이션 조건을 뒷받침하기 위해 말리아빈 미적분학과 파라메트릭스 방법을 활용한다.
• 이산화에서의 변형을 제어하기 위해 $\phi_\gamma$ 와 $\Gamma_\gamma$ 의 리프시츠 연속성과 균일 유계성을 적용한다.
• 시간 평균 커플링 전략과 지수 모멘트 경계를 활용하여 균일한 수렴 속도를 도출한다.

실험 결과

연구 질문

• RQ1 운동학적 랑주뱅 이산화에 대해 넓은 범위의 이산화 방법에 대해 시간단계에 대해 균일한 미니멀라이제이션 조건을 확립할 수 있는가?
• RQ2 이러한 방법들이 시간단계 $\gamma$ 에 대해 균일한 기하학적 포스터-라이아프노프 드리프트 조건을 만족하는가?
• RQ3 이러한 방법에 대해 수렴 속도가 시간단계 $\gamma$ 에 따라 어떻게 의존하는가?
• RQ4 유도된 수렴 경계를 통계적 응용에서 발생하는 포아송 방정식 해석에 활용할 수 있는가?
• RQ5 결과는 기울기 없는 드리프트와 비평형 역학으로 어떻게 일반화되는가?

주요 결과

• 모든 $\gamma \leq \bar{\gamma}_{U}$ 에 대해, 올리버-마루야마 및 스플리팅 기반 방법을 포함한 광범위한 방법 클래스에 대해 균일한 미니멀라이제이션 조건이 성립한다.
• 기하학적 드리프트 조건은 시간단계 $\gamma$ 에 대해 균일하게 확립되며, $R_{\gamma}\overline{\mathpzc{W}}_{\gamma,\varpi} \leq \lambda_{U}^{\gamma}\overline{\mathpzc{W}}_{\gamma,\varpi} + \gamma b_U$ 를 만족하며, 여기서 $\lambda_U^\gamma < 1$ 이다.
• 총 변동과 $V$-노름에서의 수렴 속도는 시간단계에 대해 선형적으로 스케일링되며, $\log(\rho_\gamma) \propto \gamma$ 를 만족한다.
• 비기하학적 항을 균일하게 제어하기 위해 경계 $b_U = (\mathscr{K} - \log \lambda_U) \mathrm{e}^{\mathscr{K} \bar{\gamma}_U + \varpi_U(1 + \mathfrak{C}_\phi K_U)}$ 가 보장된다.
• 리아프노프 함수 $\overline{\mathpzc{W}}_{\gamma,\varpi}$ 는 $\mathfrak{C}_\phi$ 와 $c_{\mathpzc{W}}$ 를 통해 음이 아닌 균일 유계성을 입증한다.
• $\phi_\gamma$ 와 $\Gamma_\gamma$ 의 리프시츠 연속성이 균일한 상수를 갖는 것으로 증명되어, 변형에 대한 안정성 분석이 가능해진다.