Wasserstein Stability of Contracting Flows: Effective Rates, Euler Self-Correction, and Noise Tightening

Abstract

Contraction theory guaranties exponential convergence between trajectories of a stable nonlinear system. When initial conditions are uncertain and represented as probability distributions, as in ensemble control, Bayesian estimation, and generative modeling, this guaranty extends to the distributional level via Wasserstein distance. However, the classical distributional bound is tight only for linear systems; for nonlinear dynamics, it can be significantly conservative because it collapses the spatially varying local contraction rate to a single worst-case constant, discarding distributional information entirely. We address three concrete consequences of this conservatism. First, we derive a tighter Wasserstein bound by replacing the worst-case rate with a displacement-weighted distributional average of the local contraction rate, which strictly improves upon the classical bound for every nonlinear contracting system. Second, we provide the first theoretical characterization of the self-correcting Euler discretization error under contraction: the error profile is non-monotone, peaks at a universal time that depends only on the contraction rate, and then decays exponentially, a behavior absent in non-contracting dynamics. Third, we prove that nonlinear contracting drifts always achieve strictly smaller stationary variance than a linear system sharing the same worst-case contraction rate, formally establishing the noise-rejection advantage of nonlinear controllers. All results are validated on a representative suite of one- and two-dimensional vector fields.