Abstract:Sampling a target probability distribution with an unknown normalization constant is a fundamental challenge in computational science and engineering. Recent work shows that algorithms derived by considering gradient flows in the space of probability measures open up new avenues for algorithm development. This paper makes three contributions to this sampling approach by scrutinizing the design components of such gradient flows. Any instantiation of a gradient flow for sampling needs an energy functional and a metric to determine the flow, as well as numerical approximations of the flow to derive algorithms. Our first contribution is to show that the Kullback-Leibler divergence, as an energy functional, has the unique property (among all f-divergences) that gradient flows resulting from it do not depend on the normalization constant of the target distribution. Our second contribution is to study the choice of metric from the perspective of invariance. The Fisher-Rao metric is known as the unique choice (up to scaling) that is diffeomorphism invariant. As a computationally tractable alternative, we introduce a relaxed, affine invariance property for the metrics and gradient flows. In particular, we construct various affine invariant Wasserstein and Stein gradient flows. Affine invariant gradient flows are shown to behave more favorably than their non-affine-invariant counterparts when sampling highly anisotropic distributions, in theory and by using particle methods. Our third contribution is to study, and develop efficient algorithms based on Gaussian approximations of the gradient flows; this leads to an alternative to particle methods. We establish connections between various Gaussian approximate gradient flows, discuss their relation to gradient methods arising from parametric variational inference, and study their convergence properties both theoretically and numerically.

Global Optimization Via Schrödinger-Föllmer Diffusion.

Convergence Analysis of Schr{ö}dinger-F{ö}llmer Sampler Without Convexity

Schr{ö}dinger-F{ö}llmer Sampler: Sampling without Ergodicity

Schr{\"o}dinger-F{\"o}llmer Sampler: Sampling without Ergodicity

Global Optimization Algorithm through High-Resolution Sampling

Convergence Analysis of Schr{\"o}dinger-F{\"o}llmer Sampler without Convexity

Sampling and estimation on manifolds using the Langevin diffusion

Global Convergence of Langevin Dynamics Based Algorithms for Nonconvex Optimization

Projected Langevin dynamics and a gradient flow for entropic optimal transport

Analysis of Primal-Dual Langevin Algorithms

Sampling from the Sherrington-Kirkpatrick Gibbs measure via algorithmic stochastic localization

Sampling via Gradient Flows in the Space of Probability Measures

Sampling from Mean-Field Gibbs Measures via Diffusion Processes

Solving Fredholm Integral Equations of the First Kind via Wasserstein Gradient Flows

Sampling as optimization in the space of measures: The Langevin dynamics as a composite optimization problem

Stein Boltzmann Sampling: A Variational Approach for Global Optimization

Minimizing $f$-Divergences by Interpolating Velocity Fields

A Langevin-like Sampler for Discrete Distributions

Derivative-Free Global Minimization in One Dimension: Relaxation, Monte Carlo, and Sampling

Weak approximation of Schrödinger-Föllmer diffusion

Projected Langevin Monte Carlo algorithms in non-convex and super-linear setting