Abstract / Description of output
We provide a framework to analyze the convergence of discretized kinetic Langevin dynamics for M-\nablaLipschitz, m-convex potentials. Our approach gives convergence rates of \scrO(m/M), with explicit step size restrictions, which are of the same order as the stability threshold for Gaussian targets and are valid for a large interval of the friction parameter. We apply this methodology to various integration schemes which are popular in the molecular dynamics and machine learning communities. Further, we introduce the property ``\gamma-limit convergent"" to characterize underdamped Langevin schemes that converge to overdamped dynamics in the high-friction limit and which have step size restrictions that are independent of the friction parameter; we show that this property is not generic by exhibiting methods from both the class and its complement. Finally, we provide asymptotic bias estimates for the BAOAB scheme, which remain accurate in the high-friction limit by comparison to a modified stochastic dynamics which preserves the invariant measure.
Original language | English |
---|---|
Pages (from-to) | 1226-1258 |
Number of pages | 33 |
Journal | Siam journal on numerical analysis |
Volume | 62 |
Issue number | 3 |
Early online date | 22 May 2024 |
DOIs | |
Publication status | Published - 30 Jun 2024 |
Keywords / Materials (for Non-textual outputs)
- kinetic Langevin dynamics
- MCMC sampling
- underdamped Langevin dynamics
- Wasserstein convergence