The paper introduces a modified version of a standard algorithm used in statistics and physics called Langevin dynamics, which is a method for sampling from complicated probability distributions. Standard Langevin dynamics works by combining random noise with a drift that pushes particles toward regions of high probability. The new method adds an extra ingredient: a nonlocal velocity field computed by solving a related equation called a weighted Poisson equation. This extra term acts as a kind of global transport, moving probability mass around in a coordinated way rather than relying purely on local gradient information. The key feature is that this construction only requires knowing the target distribution up to a constant, which is practically important because in many applications the normalizing constant is unknown.
The authors prove rigorously that this combined system of equations is mathematically well-behaved. Specifically, they show that starting from reasonable initial conditions, the equations have a unique smooth solution that exists for all time. This is a nontrivial result because the system couples two different types of equations together: a parabolic equation governing how the probability distribution evolves over time, and an elliptic equation that must be solved at each moment to determine the extra velocity field. The analysis is carried out on a bounded region with boundary conditions that prevent probability from leaking out.
The most striking result concerns how fast the method converges to the correct target distribution. The authors show that the difference between the current distribution and the target shrinks exponentially fast at a rate that does not require the target distribution to be log-concave, a restrictive geometric assumption that most existing convergence guarantees depend on. They also derive a clean identity describing how information-theoretic quantities like KL divergence and Fisher information decrease over time, attributing the decay to two distinct physical mechanisms corresponding to the two components of the dynamics. This suggests the method could be practically useful for sampling from the multimodal or otherwise complex distributions that are common in Bayesian inference and machine learning.