Convergence of Langevin Monte Carlo in Chi-Squared and Renyi Divergence
This provides incremental theoretical improvements for researchers in computational statistics and machine learning, enhancing convergence guarantees for sampling algorithms in high-dimensional settings.
The paper tackles the problem of sampling from a target distribution using the unadjusted Langevin Monte Carlo algorithm, proving that it achieves improved convergence rates in Chi-squared and Renyi divergence under strong dissipativity and smoothness conditions, with a rate of Õ(dε⁻¹) for strongly convex potentials.
We study sampling from a target distribution $ν_* = e^{-f}$ using the unadjusted Langevin Monte Carlo (LMC) algorithm when the potential $f$ satisfies a strong dissipativity condition and it is first-order smooth with a Lipschitz gradient. We prove that, initialized with a Gaussian random vector that has sufficiently small variance, iterating the LMC algorithm for $\widetilde{\mathcal{O}}(λ^2 dε^{-1})$ steps is sufficient to reach $ε$-neighborhood of the target in both Chi-squared and Renyi divergence, where $λ$ is the logarithmic Sobolev constant of $ν_*$. Our results do not require warm-start to deal with the exponential dimension dependency in Chi-squared divergence at initialization. In particular, for strongly convex and first-order smooth potentials, we show that the LMC algorithm achieves the rate estimate $\widetilde{\mathcal{O}}(dε^{-1})$ which improves the previously known rates in both of these metrics, under the same assumptions. Translating this rate to other metrics, our results also recover the state-of-the-art rate estimates in KL divergence, total variation and $2$-Wasserstein distance in the same setup. Finally, as we rely on the logarithmic Sobolev inequality, our framework covers a range of non-convex potentials that are first-order smooth and exhibit strong convexity outside of a compact region.