[论文解读] Uniform minorization condition and convergence bounds for discretizations of kinetic Langevin dynamics

Alain Durmus, Aurélien Enfroy+2 more|arXiv (Cornell University)|Jul 30, 2021

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.

研究动机与目标

• 为一大类时间离散化动能朗之万动力学数值格式建立统一的极小化条件。
• 推导出在时间步长上一致的几何福斯特–莱普诺夫漂移条件，以确保稳定性和收敛性。
• 提供在总变差和 $V$-范数下的定量收敛界，其缩放与步长 $\gamma$ 呈线性关系。
• 通过确保统一遍历性，支持统计推断和分子动力学中出现的泊松方程解的分析。
• 将现有收敛理论扩展至欠阻尼（动能）朗之万动力学，包括非梯度漂移和非平衡情形。

提出的方法

• 引入一个结合二次项与势能项的柯尔莫哥洛夫函数 $\overline{\mathpzc{W}}_{\gamma,\varpi}$，以捕捉类似能量的行为。
• 通过耦合论证与半椭圆性，建立对所有 $\gamma \leq \bar{\gamma}_{U}$ 成立的统一极小化条件，且与步长无关。
• 推导出几何漂移条件 $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)}$ 确保了对非几何项的统一控制。
• 通过 $\mathfrak{C}_\phi$ 与 $c_{\mathpzc{W}}$ 证明了柯尔莫哥洛夫函数 $\overline{\mathpzc{W}}_{\gamma,\varpi}$ 的非负性与一致有界性。
• 证明了 $\phi_\gamma$ 与 $\Gamma_\gamma$ 的利普希茨连续性，并获得一致常数，从而支持在扰动下的稳定性分析。