Wasserstein Contractivity

• Wasserstein contractivity is a property where dynamical evolution reduces Wasserstein distances between probability measures, ensuring stability and convergence.
• Key methodologies include coupling techniques and tailored metrics, like concave and Lyapunov-weighted distances, to achieve exponential or polynomial decay.
• Applications span kinetic equations, filtering, and MCMC, where contractivity underpins uniqueness, mixing, and quantitative convergence to equilibrium.

Wasserstein contractivity is a fundamental property of dynamical systems, Markov processes, and partial differential equations wherein the evolution in time contracts Wasserstein distances between probability measures, often implying mixing, stability, uniqueness, and quantitative convergence to equilibrium. Contractivity in Wasserstein metric can manifest as strict exponential decay, monotonicity, or, in weaker settings, as dimension-free or polynomial decay estimates. The mechanisms underlying contractivity are diverse: synchronous or reflection coupling, convexity conditions, functional inequalities, regularizing effects of noise, and structural properties of the drift and interaction. The literature provides both sharp positive results and geometric obstructions, relating contractivity to curvature (Bakry–Émery), uniform log-concavity, spectral properties, or special metric choices.

1. Definitions and General Principles

The $p$-Wasserstein distance between probability measures $\mu$, $\nu$ on a metric space $(X,d)$,

$$W_p(\mu,\nu) = \Bigl( \inf_{\gamma \in \Gamma(\mu,\nu)} \int_{X \times X} d(x,y)^p\, d\gamma(x,y) \Bigr)^{1/p},$$

quantifies the minimal $p$-th moment of the transport cost required to move $\mu$ to $\nu$. When a semigroup $(P_t)_{t \ge 0}$ (or a Markov kernel $P$) acts on measures, one says it is contractive in $\mu$0 with rate $\mu$1 if

$\mu$2

for all $\mu$3 in a suitable class. For discrete-time kernels, analogous statements involve a contraction factor $\mu$4 per step. In many applications, contractivity can be proved for modified metrics—weighted, concave, or adapted to the system's geometry—to handle weak or nonuniform dissipativity.

Contractivity is tightly coupled with the system's structural properties: monotonicity of the drift, convexity conditions, spectral properties, or the geometry of the noise. Notably, in one dimension, the $\mu$5-Wasserstein distance $\mu$6 often provides the sharpest possible mode of control, while higher dimensions may require additional regularity or strict convexity.

2. Criteria and Mechanisms for Contractivity

2.1. Synchronous and Reflection Coupling

For stochastic differential equations in $\mu$7,

$\mu$8

a classic approach is synchronous coupling (driving two solutions with the same Brownian motion). If the drift satisfies

$\mu$9

for all $\nu$0 and some positive-definite $\nu$1, parallel coupling yields

$\nu$2

and even almost sure contraction of distances along coupled trajectories (Monmarché, 2020). This criterion extends to hypoelliptic diffusions and to scenarios with nonconstant diffusion matrices, with suitable adaptation to $\nu$3 (Bouvrie et al., 2019).

Reflection coupling, as developed in (Eberle, 2013), handles cases where the drift is only contractive at infinity or not globally one-sided Lipschitz. By constructing a coupling that forces sticky encounters and exploits local concavity of a distance function, one obtains exponential contractivity in modified (concave) Kantorovich distances—interpolating between total variation and standard $\nu$4.

2.2. Modified, Concave, and Lyapunov-Weighted Metrics

Wasserstein contractivity sometimes fails for the standard $\nu$5; however, contractivity can be restored by using a distance adapted to the drift's structure, non-uniform dissipation, or domain constraints. Notable constructions include:

• Concave distances: Define $\nu$6 with $\nu$7 strictly concave and monotone, tailored to match local dissipativity and containing the necessary curvature to guarantee exponential decay (Eberle, 2013).
• Lyapunov-weighted metrics: On domains or for processes with degenerate diffusion, distances involving Lyapunov functions and convexity weights provide contraction in stochastic systems with reflection or boundary constraints (Lekang et al., 2021).

Properties such as $\nu$8 for an effective curvature $\nu$9 ensure positive contractivity rates.

2.3. Convolution and Heat Semigroup Contractivity

On Euclidean space, convolution with mollifiers (heat kernel, Gaussian smoothing) produces classical monotonicity: $(X,d)$0 for any $(X,d)$1 and all mollifiers $(X,d)$2 (Fathi et al., 4 Dec 2025). However, this is not a strict contraction in general. Only under further structure (e.g., log-concavity, positive curvature) does one recover exponential decay in $(X,d)$3 or $(X,d)$4 (Chen et al., 2020). On Riemannian manifolds with positive Ricci curvature, the heat semigroup contracts $(X,d)$5 at an exponential rate prescribed by the curvature tensor.

2.4. Markov Chains and Stein–Steinsaltz Theory

For discrete-time Markov kernels, the Kantorovich coefficient $(X,d)$6 quantifies the contraction: $(X,d)$7 If $(X,d)$8, then $(X,d)$9 is contractive, and the Wiener–Hopf (or Neumann) series directly constructs solutions to functional equations and yields exponential mixing (Hofstadler, 22 Feb 2026, Madras et al., 2011). Reflection and other drift functions enable geometric ergodicity proofs in Wasserstein metric (Madras et al., 2011) for a wide class of stochastic maps and Markov chains.

3. Applications and Model Classes

3.1. Langevin, McKean–Vlasov, Kolmogorov–Fokker–Planck, Landau Dynamics

In kinetic and mean-field models, Wasserstein contractivity underpins propagation of chaos, uniqueness, stability, and long-time behavior:

• Hypoelliptic and kinetic Fokker–Planck equations: Exponential contractivity in $$W_p(\mu,\nu) = \Bigl( \inf_{\gamma \in \Gamma(\mu,\nu)} \int_{X \times X} d(x,y)^p\, d\gamma(x,y) \Bigr)^{1/p},$$0 or $$W_p(\mu,\nu) = \Bigl( \inf_{\gamma \in \Gamma(\mu,\nu)} \int_{X \times X} d(x,y)^p\, d\gamma(x,y) \Bigr)^{1/p},$$1 is proved via duality (oscillation or coupling variables) and short-time hypocoercive smoothing (Forcillo et al., 13 Oct 2025). Construction of Lyapunov functions is often essential.
• Interacting particle systems with Lévy noise or weak interaction: Contractivity is achieved using coupled stochastic processes and metric constructions blending reflection and synchronous coupling (Schuh, 2022, Liu et al., 2024).
• Landau equation: For Maxwellian molecules, monotonicity of $$W_p(\mu,\nu) = \Bigl( \inf_{\gamma \in \Gamma(\mu,\nu)} \int_{X \times X} d(x,y)^p\, d\gamma(x,y) \Bigr)^{1/p},$$2 is proved by lifting to a higher-dimensional convexity framework and exploiting displacement convexity (Caja et al., 18 Apr 2025).
• Kac's model: Uniform contractivity in $$W_p(\mu,\nu) = \Bigl( \inf_{\gamma \in \Gamma(\mu,\nu)} \int_{X \times X} d(x,y)^p\, d\gamma(x,y) \Bigr)^{1/p},$$3 on the energy simplex, with rate equal to the spectral gap, is established by energy-coupling analysis (Hauray, 2015).

3.2. Filtering and Inference

The nonlinear filtering map for partially observed diffusions with affine drift and log-concave likelihoods is contractive in $$W_p(\mu,\nu) = \Bigl( \inf_{\gamma \in \Gamma(\mu,\nu)} \int_{X \times X} d(x,y)^p\, d\gamma(x,y) \Bigr)^{1/p},$$4, with explicit and dimension-free rates—as long as log-concavity conditions are strong enough to overcome lack of ergodicity in the signal (Whiteley, 2017). This underpins stability and robustness in high-dimensional filtering and Bayesian inference.

3.3. Markov Chain Monte Carlo and Sampling

Wasserstein contractivity is directly linked to efficient convergence and mixing in Markov chain Monte Carlo algorithms, including random walks, slice sampling, and Metropolis–Hastings. Explicit dimension-dependent rates and coupling arguments are crucial for practical tuning and for understanding the effect of interactions and drift/non-reversibility (Hofstadler, 22 Feb 2026).

4. Quantitative Rates and Rigidity, Limitations, and Obstructions

4.1. Rigidity and Limiting Rates

Contractivity factors in convolution settings cannot generally be made strict unless further symmetry, convexity, or log-concavity assumptions hold. Quantitative rigidity results for Wasserstein contraction under convolution show that near-equality forces measures to be nearly translates (for $$W_p(\mu,\nu) = \Bigl( \inf_{\gamma \in \Gamma(\mu,\nu)} \int_{X \times X} d(x,y)^p\, d\gamma(x,y) \Bigr)^{1/p},$$5) or otherwise structurally close in their marginals (Fathi et al., 4 Dec 2025). For the heat semigroup on $$W_p(\mu,\nu) = \Bigl( \inf_{\gamma \in \Gamma(\mu,\nu)} \int_{X \times X} d(x,y)^p\, d\gamma(x,y) \Bigr)^{1/p},$$6, in the absence of curvature, contraction is always polynomial, with the decay rate dictated by moment matching of the input measures (Chen et al., 2020).

4.2. Obstructions and Non-contraction

On non-Euclidean Minkowski spaces or Finsler manifolds (i.e., normed spaces where the norm is not derived from an inner product), the heat semigroup (linear diffusion) is not contractive in any $$W_p(\mu,\nu) = \Bigl( \inf_{\gamma \in \Gamma(\mu,\nu)} \int_{X \times X} d(x,y)^p\, d\gamma(x,y) \Bigr)^{1/p},$$7 unless the norm is Euclidean (Ohta et al., 2010). The fundamental obstruction is rooted in the absence of skew-convexity, or nonpositive curvature, in the displacement geometry.

For stochastic systems with only local dissipativity or non-globally convex potential, $$W_p(\mu,\nu) = \Bigl( \inf_{\gamma \in \Gamma(\mu,\nu)} \int_{X \times X} d(x,y)^p\, d\gamma(x,y) \Bigr)^{1/p},$$8 contractivity may fail, especially for $$W_p(\mu,\nu) = \Bigl( \inf_{\gamma \in \Gamma(\mu,\nu)} \int_{X \times X} d(x,y)^p\, d\gamma(x,y) \Bigr)^{1/p},$$9 unless weighted or concave-metric adaptations are used (Monmarché, 28 Feb 2026). Sharp Lyapunov exponents and Feynman–Kac operators provide exact conditions for contractivity in these regimes.

5. Uniqueness, Propagation of Chaos, and Functional Inequalities

Exponential contractivity in Wasserstein metrics typically implies uniqueness of invariant measures, mixing, and strong ergodicity. In mean-field particle systems, uniform contractivity ensures uniform-in-time propagation of chaos and quantifies finite-$p$0 convergence rates of empirical measures (Schuh, 2022, Liu et al., 2024). Contractivity is also a key ingredient in Poincaré and log-Sobolev inequalities, leading to concentration of measure and non-asymptotic variance bounds for ergodic averages (Liu et al., 2023, Forcillo et al., 13 Oct 2025).

Uniqueness and convergence are particularly robust in settings where the contractivity rate is uniform over time and particle number, and where the metric is adapted to system structure (Choi et al., 2020, Carrillo et al., 2013).

6. Methods of Proof and Technical Innovations

PDE and probabilistic proofs of Wasserstein contractivity are deeply interconnected:

• Coupling by reflection and synchronous dynamics: Directly constructs sample-path couplings realizing contraction.
• Pseudo-inverse and monotonicity arguments: For 1D dynamics or phase-coupled models (Kuramoto, Cucker–Smale), the evolution of pseudo-inverse quantile functions yields explicit decay rates (Choi et al., 2020, Carrillo et al., 2013).
• Oscillation/doubling variables techniques: In kinetic PDEs, estimate contraction by controlling weighted seminorms and exploiting adjoint equations (Forcillo et al., 13 Oct 2025).
• Concave and Lyapunov metrics: Adapt the geometry to the system’s actual dissipation.
• Rigidity theorems via duality and uniform convexity: Turn qualitative contraction bounds into quantitative stability and uniqueness results under convolution (Fathi et al., 4 Dec 2025).

7. Connections with Functional Inequalities and Information Theory

Contractivity in Wasserstein distances is closely related to Poincaré and log-Sobolev inequalities, displacement convexity of entropy, and information-theoretical decay (relative entropy, $p$1, and total variation). Decay of Fisher information in $p$2, as seen in the Fisher infinitesimal model, offers a non-classical perspective: contractivity at the level of sup-norm derivatives of log-densities, not just quadratic quantities (Calvez et al., 2023). In the rigor of Markov semigroups, decay in $p$3 frequently implies, or is implied by, functional inequalities for the generator and its invariant measure.

---

In summary, Wasserstein contractivity is a central organizing principle underlying stability, ergodicity, synchronization, and robust statistical convergence in modern stochastic analysis. Its realization and quantification require careful geometric, analytic, and probabilistic constructions, and its obstructions shed light on the structural limits of stability in high- or infinite-dimensional systems. The interplay between contractivity, curvature, coupling, and metric design continues to drive advances in kinetic theory, stochastic dynamics, sampling, and functional analysis.