arxiv: 2604.10068 · v1 · submitted 2026-04-11 · 🧮 math.AP

Recognition: unknown

Sharp hypocoercive convergence estimates for underdamped Langevin dynamics via the modified L² method

Zexi Fan , Bowen Li , Jianfeng Lu

Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:32 UTC · model grok-4.3

classification 🧮 math.AP

keywords hypocoercivityunderdamped Langevin dynamicsL² convergence ratesPoincaré inequalitymodified L² methodhypoelliptic diffusion

0 comments

The pith

The underdamped Langevin dynamics converges in L² at the explicit rate 1/(6(√(2+K/(2m))+√(4+K/(2m)))) √m under a Poincaré inequality and bounded-below Hessian.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes a sharp, explicit hypocoercive L² convergence rate for the underdamped Langevin dynamics whose invariant measure is the Gibbs measure with potential U(x) + |v|²/2. It assumes only that the position marginal satisfies a Poincaré inequality with constant m and that the Hessian of U is bounded from below by -K. The proof revisits the modified L² method by inserting a gap-shifted corrector built from the overdamped generator and the Hamiltonian part of the dynamics. When U is convex the bound simplifies to the optimal order √m scaling. The result gives a concrete quantitative guarantee on how fast the process mixes in L².

Core claim

Under the stated assumptions on U, the underdamped Langevin semigroup satisfies ||f - ∫f dμ||_{L²(μ)} ≤ C e^{-Λ t} ||f||_{L²(μ)} with the explicit constant Λ = 1/(6(√(2 + K/(2m)) + √(4 + K/(2m)))) √m, recovered by constructing a modified L² norm whose dissipation is controlled via the gap-shifted corrector A_m = (m - L_o)^{-1}(L_a Π_v)^*.

What carries the argument

The gap-shifted corrector A_m = (m - L_o)^{-1}(L_a Π_v)^* that augments the standard L² inner product to capture the hypocoercive dissipation between the overdamped generator L_o and the antisymmetric Hamiltonian part L_a.

If this is right

The L² distance to equilibrium decays exponentially with an explicitly computable rate that depends only on m and K.
When the potential is convex (K=0) the rate is optimal and scales exactly as √m.
The same corrector construction yields a concrete hypocoercive constant that can be inserted into any Lyapunov analysis relying on the modified L² norm.
The bound holds uniformly for all initial data in L²(μ) and does not require higher regularity or moment assumptions beyond the Poincaré and Hessian conditions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The explicit dependence on m and K makes the rate immediately usable for dimension-free or high-dimensional sampling analyses that track only the Poincaré constant.
The gap-shifted corrector technique may extend to other hypoelliptic diffusions whose generator splits into an elliptic part and a transport part.
For non-convex potentials the formula still supplies a positive lower bound on the rate that degrades continuously with K/m.

Load-bearing premise

The position marginal must satisfy a Poincaré inequality with positive constant m while the Hessian of U is bounded from below by -K.

What would settle it

A numerical simulation of the underdamped Langevin dynamics for a convex quadratic potential U(x) = m|x|²/2 that fails to exhibit decay at rate proportional to √m would falsify the claimed optimality.

read the original abstract

In this note, we consider the underdamped Langevin dynamics with invariant measure $\mu(\mathrm{d}x\,\mathrm{d}v) \propto e^{-U(x)-|v|^2/2}\,\mathrm{d}x\,\mathrm{d}v$. Assume that the position marginal $\mu_x(\mathrm{d}x)\propto e^{-U(x)}\,\mathrm{d}x$ satisfies a Poincar\'{e} inequality with constant $m>0$, and that $\nabla^2 U\ge -K\,\mathrm{Id}$ for some $K\ge 0$. We revisit the modified $L^2$ method of Dolbeault--Mouhot--Schmeiser, employing a gap-shifted corrector \begin{equation*} A_m=(m- L_{\mathrm{o}})^{-1}(L_a\Pi_v)^*, \end{equation*} where $L_{\mathrm{o}}=\Delta_x-\nabla U\cdot\nabla_x$ is the overdamped generator, $L_a$ is the generator of the Hamiltonian flow, and $\Pi_v$ denotes averaging over the velocity variable. We establish an explicit hypocoercive $L^2$-convergence rate \begin{equation*} \Lambda=\frac{1}{6\Bigl(\sqrt{2+\frac{K}{2m}}+\sqrt{4+\frac{K}{2m}}\Bigr)}\sqrt{m}. \end{equation*} In particular, for convex $U$, this recovers the optimal $O(\sqrt{m})$ rate.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives an explicit hypocoercive L2 rate for underdamped Langevin via a gap-shifted corrector that recovers the optimal sqrt(m) scaling when the potential is convex.

read the letter

The main takeaway is a closed-form convergence rate Lambda = 1 over 6 times (sqrt(2 + K over 2m) plus sqrt(4 + K over 2m)) times sqrt(m), derived from the modified L2 method with the specific corrector A_m = (m - L_o)^{-1} (L_a Pi_v)^*. This is new in the sense that the gap shift produces an explicit dependence on the Poincaré constant m and the semi-convexity bound K without leaving the constant implicit or requiring further optimization. For K=0 it matches the known optimal order, which is a solid check on the construction. The operator estimates follow from standard commutator bounds once the spectral gap of the overdamped generator is at least m, so the argument is not circular. The assumptions are the usual ones: Poincaré on the position marginal and a lower bound on the Hessian. The 1/6 factor looks like it comes from a few loose norm inequalities in the hypocoercivity estimate, but the paper does not claim sharpness on the prefactor, only on the scaling. That keeps the result usable even if the constant could be tightened later. The derivation stays within the Dolbeault-Mouhot-Schmeiser framework and does not introduce hidden parameters or fitted quantities. This note is aimed at people who need concrete rates for kinetic Langevin or related hypocoercive SDEs, especially in sampling or stochastic PDE contexts. The math is direct and the result is easy to check against the operator setup. It is worth sending to peer review because the explicit formula is a modest but concrete improvement over prior work that left the rate less transparent.

Referee Report

0 major / 2 minor

Summary. The paper considers the underdamped Langevin dynamics with invariant measure proportional to exp(-U(x) - |v|^2/2) dx dv. Under the assumptions that the position marginal satisfies a Poincaré inequality with constant m > 0 and that ∇²U ≥ -K Id for K ≥ 0, the authors apply the modified L² method of Dolbeault-Mouhot-Schmeiser with the gap-shifted corrector A_m = (m - L_o)^{-1} (L_a Π_v)^*, where L_o is the overdamped generator, L_a the Hamiltonian generator, and Π_v the velocity average. They derive the explicit hypocoercive L²-convergence rate Λ = 1/(6(√(2 + K/(2m)) + √(4 + K/(2m)))) √m, which recovers the optimal O(√m) scaling when U is convex (K=0).

Significance. If the central derivation holds, the result is significant because it supplies fully explicit constants in the hypocoercivity estimate, depending only on the Poincaré constant m and the semi-convexity parameter K, without additional fitted quantities. This improves quantitative understanding of convergence rates for underdamped Langevin sampling and recovers the expected √m scaling for convex potentials via a direct operator construction that avoids circularity.

minor comments (2)

[Abstract] Abstract, equation for A_m: the notation L_o, L_a and Π_v is introduced without a one-line reminder of their explicit action on test functions (e.g., L_o f = Δ_x f - ∇U · ∇_x f); adding this would improve immediate readability for readers outside the immediate hypocoercivity literature.
[Abstract] Abstract, final sentence: the assertion that the rate is 'optimal O(√m)' for convex U is stated without a brief comparison to known lower bounds or prior explicit constants; a single sentence referencing the relevant literature would strengthen the claim.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We are grateful to the referee for the positive assessment of our work on explicit hypocoercive convergence rates for underdamped Langevin dynamics. The referee correctly summarizes the main result and its significance in providing fully explicit constants depending only on m and K. We note the recommendation for minor revision and will address any editorial suggestions in the revised manuscript. Since no major comments were provided, we have no specific points to respond to in this rebuttal.

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained from assumptions

full rationale

The paper applies the modified L² method (citing external Dolbeault-Mouhot-Schmeiser) to the underdamped Langevin generator, defining the gap-shifted corrector A_m = (m - L_o)^{-1}(L_a Π_v)^* explicitly from the Poincaré gap m and the overdamped operator L_o. The hypocoercivity inequality and explicit rate Λ are then obtained via standard commutator estimates and norm bounds under the given semi-convexity ∇²U ≥ -K Id; these steps do not reduce to self-definition, fitted inputs renamed as predictions, or load-bearing self-citations. The construction is independent of the target rate and holds by direct operator analysis on the orthogonal complement.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The claim rests on two standard domain assumptions from hypocoercivity theory; no free parameters are fitted and no new entities are introduced.

axioms (2)

domain assumption Position marginal satisfies Poincaré inequality with constant m>0
Invoked to control overdamped convergence and define the gap shift in A_m.
domain assumption ∇²U ≥ -K Id for K≥0
Controls the growth of the potential and appears in the explicit constants inside Λ.

pith-pipeline@v0.9.0 · 5583 in / 1201 out tokens · 48623 ms · 2026-05-10T16:32:23.900217+00:00 · methodology

discussion (0)

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

A sharp hypocoercive entropy decay estimate for underdamped Langevin dynamics
math.AP 2026-05 unverdicted novelty 8.0

Under LSI and convexity assumptions on the target, the relative entropy of underdamped Langevin dynamics decays at an explicit rate proportional to sqrt(ρ) with constant depending on a scaled friction parameter.
Nesterov acceleration for the Wasserstein minimization of displacement-convex free energies
math.AP 2026-05 unverdicted novelty 7.0

Mean-field underdamped Langevin dynamics achieve Nesterov acceleration for Wasserstein gradient flows of displacement-convex free energies.
A sharp hypocoercive entropy decay estimate for underdamped Langevin dynamics
math.AP 2026-05 unverdicted novelty 7.0

Relative entropy of underdamped Langevin dynamics decays at explicit rate Lambda = theta sqrt(rho) / (2(1+theta)) with theta = min(Gamma/12, 1/(4Gamma)) for friction scaled as Gamma sqrt(rho), achieving optimal sqrt(r...

Reference graph

Works this paper leans on

15 extracted references · 3 canonical work pages · cited by 2 Pith papers

[1]

Albritton, S

D. Albritton, S. Armstrong, J.-C. Mourrat, and M. Novack,Variational methods for the kinetic Fokker–Planck equation, Anal. PDE17(2024), no. 6, 1953–2010

2024
[2]

Bakry, I

D. Bakry, I. Gentil, and M. Ledoux,Analysis and geometry of Markov diffusion operators, Grundlehren der mathema- tischen Wissenschaften, vol. 348, Springer, Cham, 2014

2014
[3]

Brigati, F

G. Brigati, F . Lörler, and L. Wang,Hypocoercivity meets lifts, 2025. arXiv:2412.10890

work page arXiv 2025
[4]

Y. Cao, J. Lu, and L. Wang,On explicit L 2-convergence rate estimate for underdamped Langevin dynamics, Arch. Ra- tion. Mech. Anal.247(2023), no. 5, 90

2023
[5]

1, 321–351

Loïs Delande,Sharp spectral gap of adaptive Langevin dynamics, Asymptotic Analysis142(2025), no. 1, 321–351

2025
[6]

Dolbeault, C

J. Dolbeault, C. Mouhot, and C. Schmeiser,Hypocoercivity for kinetic equations with linear relaxation terms, C. R. Math. Acad. Sci. Paris347(2009), 511–516

2009
[7]

,Hypocoercivity for linear kinetic equations conserving mass, Trans. Amer. Math. Soc.367(2015), 3807–3828

2015
[8]

Convergence of non-reversible Markov processes via lifting and flow Poincar

A. Eberle, A. Guillin, L. Hahn, F . Lörler, and M. Michel,Convergence of non-reversible Markov processes via lifting and flow Poincaré inequality, 2025. arXiv:2503.04238

work page arXiv 2025
[9]

Theory Relat

Andreas Eberle and Fabian Lörler,Non-reversible lifts of reversible diffusion processes and relaxation times, Probab. Theory Relat. Fields194(2026), 173–203

2026
[10]

Di Fang, Jianfeng Lu, and Yu Tong,Mixing time of open quantum systems via hypocoercivity, Phys. Rev. Lett.134 (2025), no. 14, 140405

2025
[11]

3-4, 349–359

Frédéric Hérau,Hypocoercivity and exponential time decay for the linear inhomogeneous relaxation Boltzmann equation, Asymptotic Analysis46(2006), no. 3-4, 349–359

2006
[12]

Bowen Li and Jianfeng Lu,Speeding up quantum markov processes through lifting, arXiv preprint arXiv:2505.12187 (2025)

work page arXiv 2025
[13]

G. A. Pavliotis,Stochastic processes and applications: Diffusion processes, the fokker-planck and langevin equations, Texts in Applied Mathematics, vol. 60, Springer, New York, 2014

2014
[14]

Roussel and G

J. Roussel and G. Stoltz,Spectral methods for Langevin dynamics and associated error estimates, ESAIM Math. Model. Numer. Anal.52(2018), no. 3, 1051–1083

2018
[15]

Villani,Hypocoercivity, Mem

C. Villani,Hypocoercivity, Mem. Amer. Math. Soc., vol. 202, 2009. SCHOOL OFMATHEMATICALSCIENCES, PEKINGUNIVERSITY Email address:2200010816@stu.pku.edu.cn DEPARTMENT OFMATHEMATICS, CITYUNIVERSITY OFHONGKONG Email address:boweli4@cityu.edu.hk DEPARTMENT OFMATHEMATICS, DEPARTMENT OFPHYSICS,ANDDEPARTMENT OFCHEMISTRY, DUKEUNIVERSITY Email address:jianfeng@math...

2009