Long-time L^p Wasserstein contraction for diffusion processes without global dissipativity
Pith reviewed 2026-05-15 18:00 UTC · model grok-4.3
The pith
Diffusion processes contract L^p Wasserstein distances in long time even without global dissipativity when the maximal eigenvalue of an associated Feynman-Kac operator is negative.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For a Markov diffusion semigroup on R^d the L^p Wasserstein distance contracts in the long-time limit along suitable synchronous couplings without any global dissipativity assumption on the drift, provided the coefficients satisfy conditions that make the maximal eigenvalue of the associated Feynman-Kac operator negative; the result covers non-elliptic cases and supplies a complete characterization at least in dimension one.
What carries the argument
Synchronous coupling of two copies of the diffusion, whose distance process is controlled by the sign of the maximal eigenvalue of the Feynman-Kac operator built from the drift and diffusion coefficients.
If this is right
- Uniform-in-time stability estimates for numerical discretizations of the SDE hold under these weaker coefficient conditions.
- Long-time behavior of the process remains stable in Wasserstein distance without requiring strong global dissipativity.
- Negative results delimit the precise boundary between contraction and failure, including explicit lower bounds on the contraction rate.
- In dimension one the eigenvalue condition is necessary as well as sufficient.
Where Pith is reading between the lines
- The same eigenvalue test could be used to decide long-time stability for SDEs arising in applications where global dissipativity is unrealistic.
- The approach might extend to other metrics or to processes with jumps once the corresponding Feynman-Kac operator is identified.
- Numerical computation of the eigenvalue for concrete drifts would give immediate, verifiable predictions for contraction or its absence.
Load-bearing premise
Suitable synchronous couplings exist and the drift-diffusion coefficients satisfy the stated conditions that let the contraction be governed by the sign of the Feynman-Kac eigenvalue.
What would settle it
A concrete one-dimensional diffusion whose associated Feynman-Kac operator has positive maximal eigenvalue yet whose L^p Wasserstein distance still contracts (or conversely, negative eigenvalue yet no contraction).
Figures
read the original abstract
The fact that a Markov diffusion semi-group on $\mathbb R^d$ contracts the $L^p$ Wasserstein distance, which has been extensively used to establish uniform-in-time stability estimates (e.g. with respect to numerical discretization errors), is a well-studied question in the case where the distances are in fact deterministically contracted by the drift (global dissipativity condition) or in the case $p=1$ (with reflection couplings). This work focuses on the non-globally dissipative case with $p>1$. This situation was previously considered in \cite{MonmarcheBruit}, but only for elliptic processes, and with a restriction on the diffusivity coefficient (which had to be large enough). Here, we extend this analysis to non-elliptic processes and provide sharper conditions to get contractions along synchronous coupling, including negative results, lower bounds and a characterization (at least in dimension 1) in terms of the maximal eigenvalue of a Feynman-Kac operator.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript establishes long-time L^p Wasserstein contractions for Markov diffusion semigroups on R^d in the absence of global dissipativity of the drift. It extends prior elliptic results to non-elliptic processes via synchronous couplings, supplies sharper coefficient conditions for contraction, includes negative results and lower bounds, and provides (at least in dimension 1) a characterization of the contraction threshold in terms of the maximal eigenvalue of an associated Feynman-Kac operator.
Significance. If the derivations hold, the work supplies sharper, more applicable criteria for uniform-in-time stability estimates of diffusions beyond the globally dissipative regime, directly relevant to numerical discretization error bounds and ergodicity analysis. The inclusion of negative results together with the eigenvalue characterization strengthens the sharpness of the theory and supplies falsifiable predictions.
major comments (2)
- [§3] §3 (synchronous coupling construction): the argument that the quadratic variation term remains controlled when the diffusion matrix degenerates on a positive-measure set relies on the stated drift-diffusion conditions bounding the generator applied to |x-y|^p; however, this does not automatically guarantee that a pathwise synchronous coupling can be realized on the degeneracy set without an additional uniform controllability assumption on the support of the diffusion. If this gap is not closed, the long-time contraction claim for non-elliptic cases fails to hold pathwise.
- [§5] §5 (1D Feynman-Kac characterization): the maximal-eigenvalue criterion is derived under non-degenerate assumptions on the diffusion coefficient; the manuscript must verify whether the same spectral characterization transfers directly to the degenerate case or requires a modified operator, since the negative-result examples in higher dimensions may indicate that the 1D reduction does not capture the full degeneracy effect.
minor comments (2)
- [§2] Notation for the degeneracy set of the diffusion matrix should be introduced explicitly in §2 and used consistently when stating the coefficient conditions.
- [Abstract] The abstract claims 'sharper conditions' relative to MonmarcheBruit; a one-sentence comparison of the new threshold with the previous diffusivity lower bound would improve readability.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below with point-by-point responses, providing clarifications where needed and indicating revisions to strengthen the presentation.
read point-by-point responses
-
Referee: [§3] §3 (synchronous coupling construction): the argument that the quadratic variation term remains controlled when the diffusion matrix degenerates on a positive-measure set relies on the stated drift-diffusion conditions bounding the generator applied to |x-y|^p; however, this does not automatically guarantee that a pathwise synchronous coupling can be realized on the degeneracy set without an additional uniform controllability assumption on the support of the diffusion. If this gap is not closed, the long-time contraction claim for non-elliptic cases fails to hold pathwise.
Authors: We appreciate the referee's observation on the synchronous coupling construction. The L^p-Wasserstein distance is controlled through the evolution of the expectation E[|X_t - Y_t|^p]. Applying Itô's formula to |x-y|^p yields a drift term bounded by the generator under our coefficient assumptions, while the quadratic variation term is nonnegative and can be discarded to obtain an upper bound. The synchronous coupling is realized by driving both processes with the same Brownian motion; degeneracy of the diffusion matrix on a set of positive measure does not invalidate this construction or the resulting differential inequality for the expectation, as the processes remain well-defined under the local Lipschitz and linear growth conditions stated in the paper. Since the contraction result concerns the Wasserstein distance (an expectation), no pathwise controllability on the degeneracy set is required. We will add a clarifying remark in §3 emphasizing that the argument relies on the generator bound for the expectation and does not claim deterministic pathwise contraction. revision: partial
-
Referee: [§5] §5 (1D Feynman-Kac characterization): the maximal-eigenvalue criterion is derived under non-degenerate assumptions on the diffusion coefficient; the manuscript must verify whether the same spectral characterization transfers directly to the degenerate case or requires a modified operator, since the negative-result examples in higher dimensions may indicate that the 1D reduction does not capture the full degeneracy effect.
Authors: The derivation of the maximal-eigenvalue criterion in §5 is carried out under the assumption that the diffusion coefficient is bounded away from zero, ensuring the Feynman-Kac operator is uniformly elliptic and the standard spectral theory applies directly. In the presence of degeneracy, the operator would require modification (for instance, by restricting the domain or incorporating a weighted space that accounts for vanishing diffusion). Our negative results are formulated in dimensions d ≥ 2 and illustrate the failure of contraction under certain local dissipativity conditions; they do not contradict the one-dimensional non-degenerate characterization. The coupling-based contraction results in the paper hold for non-elliptic processes without relying on the spectral criterion. We will revise §5 to state the non-degeneracy assumption explicitly and add a short discussion noting that a full extension of the eigenvalue characterization to degenerate one-dimensional diffusions lies outside the present scope. revision: yes
Circularity Check
No significant circularity; extensions rely on independent analysis of couplings and generators
full rationale
The paper extends prior elliptic results from self-cited work to non-elliptic diffusions by deriving sharper contraction conditions along synchronous couplings, including negative results and a 1D characterization via the maximal eigenvalue of a Feynman-Kac operator. No step reduces a claimed prediction or uniqueness result to a fitted parameter, self-referential definition, or load-bearing self-citation chain; the generator bounds and coupling constructions are developed from the stated drift/diffusion assumptions without tautological reduction. The derivation remains self-contained against external benchmarks such as explicit 1D examples and lower-bound constructions.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Existence and uniqueness of strong solutions to the underlying SDE
- domain assumption The process is a Markov diffusion
Reference graph
Works this paper leans on
-
[1]
Springer Science & Business Media, 2005
Luigi Ambrosio, Nicola Gigli, and Giuseppe Savar´ e.Gradient flows: in metric spaces and in the space of probability measures. Springer Science & Business Media, 2005
work page 2005
-
[2]
Fran¸ cois Bolley, Arnaud Guillin, and Florent Malrieu. Trend to equilibrium and parti- cle approximation for a weakly selfconsistent Vlasov-Fokker-Planck equation.ESAIM: Mathematical Modelling and Numerical Analysis, 44(5):867–884, 2010
work page 2010
-
[3]
Ivan Burenev, Dani¨ el Cloete, Vansh Kharbanda, and Hugo Touchette. An introduction to large deviations with applications in physics.SciPost Physics Lecture Notes, page 104, 2025
work page 2025
-
[4]
Patrick Cattiaux, Max Fathi, and Arnaud Guillin. Self-improvement of the Bakry-Emery criterion for Poincar´ e inequalities and Wasserstein contraction using variable curvature bounds.Journal de Math´ ematiques Pures et Appliqu´ ees, 166:1–29, 2022
work page 2022
-
[5]
Martin Chak and Pierre Monmarch´ e. Reflection coupling for unadjusted generalized Hamiltonian Monte Carlo in the nonconvex stochastic gradient case.IMA Journal of Numerical Analysis, page draf045, 2025
work page 2025
-
[6]
Chatterji, Yasin Abbasi-Yadkori, Peter L
Xiang Cheng, Niladri S Chatterji, Yasin Abbasi-Yadkori, Peter L Bartlett, and Michael I Jordan. Sharp convergence rates for Langevin dynamics in the nonconvex setting.arXiv preprint arXiv:1805.01648, 2018
-
[7]
Underdamped Langevin MCMC: A non-asymptotic analysis
Xiang Cheng, Niladri S Chatterji, Peter L Bartlett, and Michael I Jordan. Underdamped Langevin MCMC: A non-asymptotic analysis. InConference on learning theory, pages 300–323. PMLR, 2018
work page 2018
-
[8]
Arnak S Dalalyan, Avetik Karagulyan, and Lionel Riou-Durand. Bounding the error of discretized Langevin algorithms for non-strongly log-concave targets.Journal of Machine Learning Research, 23(235):1–38, 2022. 35
work page 2022
-
[9]
On the Kantorovich contraction of Markov semi- groups.arXiv preprint arXiv:2511.08111, 2025
Pierre Del Moral and Mathieu Gerber. On the Kantorovich contraction of Markov semi- groups.arXiv preprint arXiv:2511.08111, 2025
-
[10]
Hacene Djellout, Arnaud Guillin, and Liming Wu. Transportation cost-information in- equalities and applications to random dynamical systems and diffusions.The Annals of Probability, 32(3B):2702 – 2732, 2004
work page 2004
-
[11]
Alain Durmus and Andreas Eberle. Asymptotic bias of inexact Markov chain Monte Carlo methods in high dimension.The Annals of Applied Probability, 34(4):3435 – 3468, 2024
work page 2024
-
[12]
Alain Durmus, Andreas Eberle, Aur´ elien Enfroy, Arnaud Guillin, and Pierre Monmarch´ e. Discrete sticky couplings of functional autoregressive processes.The Annals of Applied Probability, 34(6):5032–5075, 2024
work page 2024
-
[13]
Alain Durmus, Andreas Eberle, Arnaud Guillin, and Raphael Zimmer. An elementary approach to uniform in time propagation of chaos.Proceedings of the American Mathe- matical Society, 148(12):5387–5398, 2020
work page 2020
-
[14]
Reflection coupling and Wasserstein contractivity without convexity
Andreas Eberle. Reflection coupling and Wasserstein contractivity without convexity. Comptes Rendus Mathematique, 349(19):1101–1104, 2011
work page 2011
-
[15]
Reflection couplings and contraction rates for diffusions.Probab
Andreas Eberle. Reflection couplings and contraction rates for diffusions.Probab. Theory Relat. Fields, 166:851–886, 2016
work page 2016
-
[16]
Andreas Eberle, Arnaud Guillin, and Raphael Zimmer. Couplings and quantitative con- traction rates for Langevin dynamics.The Annals of Probability, 47(4):1982 – 2010, 2019
work page 1982
-
[17]
Andreas Eberle and Mateusz B. Majka. Quantitative contraction rates for Markov chains on general state spaces.Electronic Journal of Probability, 24(none):1 – 36, 2019
work page 2019
-
[18]
Gr´ egoire Ferr´ e and Gabriel Stoltz. Large deviations of empirical measures of diffusions in weighted topologies.Electronic Journal of Probability, 25(none):1 – 52, 2020
work page 2020
-
[19]
Arnaud Guillin and Pierre Monmarch´ e. Uniform long-time and propagation of chaos estimates for mean field kinetic particles in non-convex landscapes.J Stat Phys, 185, 2021
work page 2021
-
[20]
Arnaud Guillin and Feng-Yu Wang. Degenerate Fokker–Planck equations: Bismut for- mula, gradient estimate and Harnack inequality.Journal of Differential Equations, 253(1):20–40, 2012
work page 2012
- [21]
-
[22]
Martin Hanke. Stochastic modeling of stationary scalar gaussian processes in continuous time from autocorrelation data.Advances in Computational Mathematics, 50(4):60, 2024
work page 2024
-
[23]
Maxime Herda, Pierre Monmarch´ e, and Benoˆ ıt Perthame. Wasserstein contraction for the stochastic Morris-Lecar neuron model.Kinetic and Related Models, 18(1):1–18, 2025. 36
work page 2025
-
[24]
Xing Huang, Eva Kopfer, Pierre Monmarch´ e, and Panpan Ren. Exponential Ergodicity in Relative Entropy andL 2-Wasserstein Distance for non-equilibrium partially dissipative Kinetic SDEs.arXiv e-prints, page arXiv:2507.02518, July 2025
-
[25]
Ald´ eric Joulin and Yann Ollivier. Curvature, concentration and error estimates for Markov chain Monte Carlo.The Annals of Probability, 38(6):2418 – 2442, 2010
work page 2010
-
[26]
Duality on gradient estimates and Wasserstein controls.J
Kazumasa Kuwada. Duality on gradient estimates and Wasserstein controls.J. Funct. Anal., 258(11):3758–3774, 2010
work page 2010
-
[27]
Tony Leli` evre, Xuyang Lin, and Pierre Monmarch´ e. Convergence rates for an Adaptive Biasing Potential scheme from a Wasserstein optimization perspective.arXiv preprint arXiv:2501.17979, 2025
-
[28]
Using Witten Laplacians to locate index-1 saddle points
Tony Leli` evre and Panos Parpas. Using Witten Laplacians to locate index-1 saddle points. SIAM Journal on Scientific Computing, 46(2):A770–A797, 2024
work page 2024
-
[29]
Linshan Liu, Matteus B. Majka, and Pierre Monmarch´ e. L2-Wasserstein contraction for Euler schemes of elliptic diffusions and interacting particle systems.Stochastic Processes and their Applications, 179:104504, 2025
work page 2025
-
[30]
A Bakry-´Emery Approach to Lipschitz Transportation on Manifolds
Pablo L´ opez-Rivera. A Bakry-´Emery Approach to Lipschitz Transportation on Manifolds. Potential Analysis, 62(2):331–353, 2025
work page 2025
-
[31]
Dejun Luo and Jian Wang. Exponential convergence in-Wasserstein distance for diffu- sion processes without uniformly dissipative drift.Mathematische Nachrichten, 289(14- 15):1909–1926, 2016
work page 1909
-
[32]
Florent Malrieu. Convergence to equilibrium for granular media equations and their Euler schemes.The Annals of Applied Probability, 13(2):540–560, 2003
work page 2003
-
[33]
Long-time behaviour and propagation of chaos for mean field kinetic particles.Stochastic Process
Pierre Monmarch´ e. Long-time behaviour and propagation of chaos for mean field kinetic particles.Stochastic Process. Appl., 127(6):1721–1737, 2017
work page 2017
-
[34]
Pierre Monmarch´ e. High-dimensional MCMC with a standard splitting scheme for the underdamped Langevin diffusion.Electronic Journal of Statistics, 15(2):4117 – 4166, 2021
work page 2021
-
[35]
Almost sure contraction for diffusions onR d
Pierre Monmarch´ e. Almost sure contraction for diffusions onR d. Application to gener- alised Langevin diffusions.Stochastic Processes and their Applications, 2023
work page 2023
-
[36]
Pierre Monmarch´ e. Wasserstein contraction and Poincar´ e inequalities for elliptic diffu- sions with high diffusivity.Annales Henri Lebesgue, 6:941–973, 2023
work page 2023
-
[37]
Pierre Monmarch´ e. Long-time propagation of chaos and exit times for metastable mean- field particle systems.arXiv e-prints, page arXiv:2503.00157, February 2025
-
[38]
Pierre Monmarch´ e, Zhenjie Ren, and Songbo Wang. Time-uniform log-Sobolev inequali- ties and applications to propagation of chaos.Electronic Journal of Probability, 29(none):1 – 38, 2024
work page 2024
-
[39]
Pierre Monmarch´ e and Songbo Wang. Logarithmic Sobolev Inequalities for Non- equilibrium Steady States.Potential Analysis, 63(3):1381–1405, 2025. 37
work page 2025
-
[40]
Felix Otto and C´ edric Villani. Generalization of an inequality by Talagrand and links with the logarithmic Sobolev inequality.Journal of Functional Analysis, 173(2):361–400, 2000
work page 2000
-
[41]
Panpan Ren and Feng-Yu Wang. Exponential convergence in entropy and Wasserstein for McKean–Vlasov SDEs.Nonlinear Analysis, 206:112259, 2021
work page 2021
-
[42]
R´ emi Rhodes and Vincent Vargas. Gaussian multiplicative chaos and applications: A review.Probability Surveys, 11(none):315 – 392, 2014
work page 2014
-
[43]
Michael R¨ ockner and Feng-Yu Wang. Log-Harnack inequality for stochastic differential equations in Hilbert spaces and its consequences.Infinite Dimensional Analysis, Quantum Probability and Related Topics, 13(01):27–37, 2010
work page 2010
-
[44]
Justin Salez. Modern aspects of Markov chains: entropy, curvature and the cutoff phe- nomenon.arXiv preprint arXiv:2508.21055, 2025
-
[45]
Katharina Schuh. Global contractivity for Langevin dynamics with distribution- dependent forces and uniform in time propagation of chaos.Annales de l’Institut Henri Poincare (B) Probabilites et statistiques, 60(2):753–789, 2024
work page 2024
-
[46]
von Renesse and Karl-Theodor Sturm
Max-K. von Renesse and Karl-Theodor Sturm. Transport inequalities, gradient esti- mates, entropy and Ricci curvature.Communications on Pure and Applied Mathematics, 58(7):923–940, 2005
work page 2005
-
[47]
Application of coupling methods to the Neumann eigenvalue problem
Feng-Yu Wang. Application of coupling methods to the Neumann eigenvalue problem. Probability Theory and Related Fields, 98:299–306, 1994
work page 1994
-
[48]
Liming Wu. A deviation inequality for non-reversible Markov processes.Annales de l’Institut Henri Poincare (B) Probability and Statistics, 36(4):435–445, 2000. 38
work page 2000
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.