← Back to arXiv
arXivProbabilityarXiv:2608.17022

Sharp hypocoercive convergence estimates for underdamped Langevin dynamics with specular reflection

The paper analyzes a particle simulation technique called underdamped Langevin dynamics, which is used to draw random samples from complicated probability distributions while respecting physical boundaries. Imagine a ball rolling around inside a convex container, bouncing off the walls like a perfect mirror reflection, while also being subject to random kicks and friction. This setup is useful in statistics and machine learning when you want to sample from a distribution but your variables are constrained to a bounded region. The central question is: how quickly does the distribution of the particle's position and momentum settle down to the desired target distribution?

The main result is the first explicit, mathematically rigorous convergence rate for this bouncing-ball sampler. The authors show that the distribution converges exponentially fast to the target, and they pin down exactly how the speed depends on the geometry of the problem. Specifically, the convergence rate scales as the square root of a quantity called the Poincare constant, which measures how "spread out" the target distribution is. This square-root scaling turns out to be twice as fast as a simpler, non-momentum-based version of the sampler in the regime where convergence is slowest. This mirrors a well-known speedup that momentum methods provide in unconstrained settings, confirming that adding momentum helps even when walls are present.

The technical challenge comes from the boundary. Standard mathematical tools for analyzing these dynamics assume the particle moves freely through all of space, and they break down when mirror reflections at a boundary are introduced. The authors adapt a sophisticated framework called hypocoercivity, which handles the interplay between the random noise and the deterministic motion, and replace a key algebraic identity used in the wall-free proof with a geometric version called a weighted Reilly formula. This formula brings in the curvature of the boundary wall, and the fact that the container is convex guarantees that this curvature term contributes favorably, keeping the whole argument consistent.

Read original →