---
title: Underdamped Langevin Diffusion
url: https://www.emergentmind.com/topics/underdamped-langevin-diffusion
type: topic
---

# Underdamped Langevin Diffusion

Underdamped Langevin diffusion is a phase-space diffusion that evolves both a configuration variable and a momentum or velocity variable, thereby retaining inertia, friction, and thermal noise. In its standard sampling form,
$$
dx_t = p_t\,dt,\qquad dp_t = -\nabla f(x_t)\,dt - \gamma p_t\,dt + \sqrt{2\gamma}\,dB_t,
$$
or, more generally,
$$
dx_t = v_t\,dt,\qquad dv_t = -\gamma v_t\,dt - u\nabla f(x_t)\,dt + \sqrt{2\gamma u}\,dB_t,
$$
and it admits the Gibbs stationary law
$$
\Pi(x,p)\propto \exp\!\left(-f(x)-\frac{1}{2}\|p\|^2\right),
$$
whose \(x\)-marginal is the target distribution \(\pi(x)\propto e^{-f(x)}\) [2410.08987], [1707.03663], [2101.03446]. Because noise enters only through momentum, the process is hypoelliptic rather than elliptic, and this structural feature governs its convergence theory, transport asymptotics, thermodynamic properties, and numerical discretizations [2503.12610].

## 1. Phase-space formulation and invariant structure

The canonical physical form of the underdamped Langevin equation writes
$$
d\boldsymbol{x} = \frac{\boldsymbol{p}}{m}\,dt,\qquad
d\boldsymbol{p} = -\gamma \frac{\boldsymbol{p}}{m}\,dt + \boldsymbol{F}(\boldsymbol{x},t)\,dt + d\boldsymbol{B}_t,
$$
with Brownian forcing of covariance \(2D\,dt\), \(D=\gamma T\), and corresponding Kramers dynamics for the phase-space density \(f(x,p,t)\) [2405.12305]. In the sampling literature, the same structure is commonly rewritten with \(f\) as the negative log target density and with unit mass, yielding a Hamiltonian
$$
H(x,p)=f(x)+\frac{1}{2}\|p\|^2
$$
and a compact drift-diffusion representation in \((x,p)\) variables [2410.08987]. The invariant measure is then Gibbsian on phase space and Gaussian in the momentum coordinate, a fact used both in MCMC analysis and in nonequilibrium thermodynamics [2410.08987], [2102.09115].

The explicit momentum variable is the essential distinction from overdamped dynamics. It introduces inertial transport through \(dx_t=p_t\,dt\), suppresses purely diffusive random-walk motion, and places underdamped Langevin methods in close conceptual proximity to Hamiltonian Monte Carlo; one non-asymptotic analysis states that the underdamped Langevin MCMC scheme “can be viewed as a version of Hamiltonian Monte Carlo (HMC)” [1707.03663]. A plausible implication is that phase-space geometry, rather than position-space drift alone, is the correct organizing principle for understanding acceleration, metastability, and response in these systems.

## 2. Equilibration, transport coefficients, and singular limits

In generalized Langevin settings with colored noise and auxiliary variables, the same underdamped degeneracy persists in a stronger form: dissipation is degenerate and standard coercive arguments fail. For the quasi-Markovian GLE in a periodic potential, sharp longtime equilibration estimates are obtained by hypocoercivity, with exponential decay in a twisted Sobolev norm equivalent to a standard Sobolev norm:
$$
\|e^{t\mathcal{L}} f\|_{H^1_0(\mu),\gamma,\nu}
\le
\exp\!\left(-C_1 \min\!\left(\gamma,\frac{1}{\gamma},\frac{\gamma}{\nu^4}\right)t\right)
\|f\|_{H^1_0(\mu),\gamma,\nu},
$$
for mean-zero \(f\) [2007.16087]. The same work formulates effective diffusion through the Poisson problem
$$
D=\int \phi\,p\,d\mu,\qquad -\mathcal{L}\phi=p,
$$
which is the basic homogenization object for long-time transport in periodic media [2007.16087].

In the underdamped limit of that GLE, corresponding to \(\gamma\to 0\), the effective diffusion coefficient obeys
$$
D_{\gamma,\nu}\sim \frac{D_\nu^*}{\gamma},
$$
so diffusion diverges as friction vanishes, while the prefactor retains explicit dependence on the memory kernel and the periodic potential [2007.16087]. In the white-noise limit \(\nu\to 0\), the prefactor converges to the classical underdamped Langevin value \(D_{\mathrm{und}}\), and numerical results further report
$$
|D_\nu^*-D_{\mathrm{und}}|=\mathcal{O}(\nu^4),
$$
confirming the theoretical expansion [2007.16087]. This suggests a robust separation between universal scaling in friction and non-universal dependence on temporal correlations of the noise.

## 3. Discretization, MCMC, and algorithmic variants

For strongly log-concave targets with smooth gradients, underdamped Langevin diffusion admits a sharper non-asymptotic MCMC theory than overdamped Langevin. One foundational result proves that a discretization of the underdamped process achieves \(2\)-Wasserstein error \(\varepsilon\) in
\[
\tilde{\mathcal{O}}\!\left(\frac{\sqrt{d}\,\kappa^2}{\varepsilon}\right)
\]
steps, compared with
\[
\mathcal{O}\!\left(\frac{d\,\kappa^2}{\varepsilon^2}\right)
\]
for overdamped Langevin under the same smoothness and strong convexity assumptions [1707.03663]. The analysis is based on contraction of the continuous-time semigroup,
$$
W_2(\Phi_t q_0,q^*)\le e^{-t/(2\kappa)}\,W_2(q_0,q^*),
$$
followed by explicit control of discretization error [1707.03663].

A parallel line of work studies geometric splitting methods for the discrete chain itself. The OBABO integrator is a symmetric Strang splitting for the kinetic Langevin SDE and requires only one gradient evaluation per iteration [2007.05455]. For convex and smooth targets, the discrete-time chain enjoys dimension-free Wasserstein contraction and explicit total variation regularization, and the unadjusted chain has Wasserstein efficiency bounds of order \(\sqrt d/\varepsilon\) in the general case, \(\sqrt{d/\varepsilon}\) if the Hessian is Lipschitz, and \(d^{1/4}/\sqrt{\varepsilon}\) for separable targets [2007.05455]. The same paper establishes non-asymptotic bounds for a Metropolis-adjusted version as well [2007.05455].

Higher-order gradient-only schemes exploit additional smoothness of the target. The shifted ODE method matches the leading stochastic Taylor terms of underdamped Langevin, uses only gradients, and becomes practical through a third-order Runge–Kutta implementation (SORT) or a fourth-order splitting implementation (SOFA); the reported per-step overhead for SORT is two additional gradient evaluations, and logistic-regression experiments show faster convergence than competing unadjusted Langevin algorithms [2101.03446]. More recently, QUICSORT is presented as the first gradient-only ULD method with third-order convergence, with step complexities
\[
\mathcal{O}(\sqrt d/\varepsilon),\qquad
\mathcal{O}(\sqrt d/\sqrt{\varepsilon}),\qquad
\mathcal{O}(\sqrt d/\varepsilon^{1/3})
\]
under, respectively, Lipschitz gradient, Lipschitz Hessian, and Lipschitz third derivative assumptions [2508.16485].

Several variants modify the basic dynamics while preserving the same invariant law. Gradient-adjusted underdamped Langevin dynamics add a term \(-a\nabla f(x_t)\,dt\) to the \(dx_t\) equation together with matching noise, retain the Gibbs stationary distribution, and for Gaussian targets yield discrete mixing time toward a biased target with square-root dependence on the condition number of the covariance matrix [2410.08987]. Skew-symmetric perturbations in both position and momentum likewise preserve the invariant Gibbs measure and can reduce asymptotic variance, especially for Gaussian targets, provided the perturbations are properly aligned with the observable and covariance structure [1705.00170]. Underdamped transitions have also been incorporated into variational inference: Langevin Diffusion Variational Inference numerically simulates the underdamped forward diffusion and its time reversal, augments backward transitions with a score network, and reports consistently stronger ELBOs than ULA, MCD, and UHA across logistic-regression and time-series tasks [2208.07743].

## 4. Anomalous transport and heterogeneous environments

Underdamped Langevin diffusion has also served as a vehicle for anomalous transport in heterogeneous media. In an isothermal one-dimensional medium with spatially dependent diffusion coefficient \(D(x)\sim |x|^c\), analytical diffusion-equation solutions and numerical simulations of the underdamped Langevin equation agree at long times for \(c<1\), yielding subdiffusion for \(c<0\), normal diffusion for \(c=0\), and superdiffusion for \(0<c<1\) [1606.08632]. The same study shows that for \(c>1\) the diffusion equation predicts faster-than-ballistic spreading, but this is unphysical: in the full underdamped dynamics, motion cannot exceed the ballistic limit \(\langle x^2\rangle\sim t^2\), because particles eventually coast inertially once friction becomes negligible [1606.08632]. A recurrent misconception is therefore that position-only diffusion equations remain valid when spatially varying friction destroys the time-scale separation needed for overdamped closure.

A different failure of the overdamped description appears in underdamped scaled Brownian motion, where temperature, damping, and diffusion coefficients depend explicitly on time,
$$
T(t)=T_0\left(1+\frac{t}{\tau_0}\right)^{2\alpha-2},\qquad
\gamma(t)=\gamma_0\left(1+\frac{t}{\tau_0}\right)^{\alpha-1},\qquad
D(t)=D_0\left(1+\frac{t}{\tau_0}\right)^{\alpha-1}.
$$
For this model, the overdamped limit can fail to describe the long-time behavior, and for small \(\alpha\) or the ultraslow case \(\alpha=0\), inertial corrections remain significant over practical observation windows [1609.07250]. The paper connects this behavior to free cooling granular gases and emphasizes persistent inertia in both ensemble and time-averaged statistics [1609.07250].

Non-Markovian timing mechanisms generate further anomalies. A weakly damped Langevin system coupled to an \(\alpha\)-dependent subordinator with \(1<\alpha<2\) leads to a fractional Klein–Kramers equation, and its mean-squared displacement depends crucially on the two-point distribution of the inverse subordinator [1810.04332]. For \(\gamma=0\), the model yields super-ballistic diffusion,
\[
\langle x^2(t)\rangle\propto t^3,
\]
whereas for \(\gamma\neq 0\) it gives sub-ballistic superdiffusion,
\[
\langle x^2(t)\rangle\propto t^{3-\alpha},\qquad 1<\alpha<2,
\]
together with an ergodicity-breaking parameter
\[
\mathcal{EB}=\frac{1}{\alpha-1}.
\]
These features motivate the description “Lévy-walk-like” [1810.04332].

Randomization of medium parameters produces another class of anomalies. In an underdamped Langevin system with time-dependent diffusivity \(\nu t^{\beta-1}\) and a random relaxation timescale \(\tau\), a heavy-tailed distribution of \(\tau\) with finite mean suppresses the decay of the velocity correlation function and promotes diffusion [2110.05846]. The ensemble MSD scales as \(t^\beta\) for fixed \(\tau\) and \(v_0\), as \(t^{1-\alpha+\beta}\) after averaging over heavy-tailed \(\tau\), and as \(t^{2-\alpha}\) in the special random-initialization regime \(v_0^2=\nu\tau\) [2110.05846]. The same model exhibits aging and broken ergodicity in the time-averaged MSD [2110.05846]. This suggests that underdamped diffusion is particularly effective at disentangling inertia, quenched heterogeneity, and temporal nonstationarity.

## 5. Heat, entropy production, and periodically driven states

In stochastic thermodynamics, underdamped Langevin diffusion differs from the overdamped theory because heat includes explicit kinetic-energy boundary terms. For the free particle, linear potential, and harmonic potential in arbitrary dimension, the heat over \([0,\tau]\) can be written as
$$
Q[\vec r]=\frac{1}{2}m(v_\tau^2-v_0^2)+V(\vec r_\tau)-V(\vec r_0),
$$
and exact characteristic functions and probability distributions are available in each case [2102.09115]. The free-particle heat law is Besselian, the free and harmonic cases are symmetric in \(Q\), and the linear potential produces an asymmetric distribution because the external field shifts the mean heat [2102.09115]. The same analysis emphasizes leptokurtic tails and “heat leakage” from kinetic fluctuations that are absent in overdamped reductions [2102.09115].

Precision bounds for fluctuating currents must also be modified in the finite-mass regime. For steady-state underdamped Langevin dynamics, the overdamped thermodynamic uncertainty relation fails in general, and the relevant bound becomes
$$
\eta_J \le \frac{1}{2}\Sigma,
$$
where \(\Sigma\) contains the overdamped entropy-production term together with two additional positive contributions associated with local mean acceleration and velocity fluctuations, plus a transient term \(\Omega_\nu\) [2202.10696]. The theory extends to magnetic fields, anisotropic temperature, and joint bounds for several observables, and in periodic potentials the correlation TUR can make the bound tight when correlations with auxiliary observables are incorporated [2202.10696].

A separate modified TUR addresses entropy-production inference directly. For underdamped Langevin systems, one may define a cumulant current \(J(w)\) and a stochastic current \(J_S(w)\) so that
$$
\frac{\Sigma}{2}\ge \frac{\langle J(w)\rangle^2}{\operatorname{var}(J_S(w))}.
$$
The bound is saturable when \(w\) is proportional to the irreversible current \(u(x,p,t)\), and the method requires only the damping-to-mass ratio \(\gamma/m\) and the diffusion constant \(D\), rather than full knowledge of the force field or stationary density [2405.12305]. The same work states that the construction applies to both overdamped and underdamped Langevin dynamics, for arbitrary time-dependent protocols and nonsteady scenarios [2405.12305].

Periodic driving produces explicitly time-dependent nonequilibrium asymptotic states. For a harmonic underdamped Langevin system with \(T\)-periodic friction \(\gamma(t)\), stiffness \(k(t)\), and diffusion \(D(t)\), the asymptotic distribution is Gaussian with \(T\)-periodic covariance matrix \(\Sigma(t)\),
$$
P^{\mathrm{os}}(x,v,t)=\frac{1}{2\pi\sqrt{|\Sigma(t)|}}
\exp\!\left[-\frac{1}{2}(x,v)\Sigma^{-1}(t)(x,v)^T\right],
$$
and the full nonperturbative dynamics can be reduced via an \(SL_2\) symmetry to a Hill equation with periodic coefficients [1911.03106]. This framework yields exact or near-exact asymptotic distributions, two-time correlation functions, and linear responses to drift and diffusion perturbations [1911.03106].

## 6. Metastability, applications, and kinetic mean-field extensions

At low temperature, underdamped Langevin diffusion exhibits metastable transitions between basins of attraction. For a double-well Morse potential with metastable minimum \(m\), global minimum \(s\), and saddle \(\sigma\), the sharp Eyring–Kramers asymptotic for the mean transition time from \((m,0)\) to a neighborhood of \(s\) is
$$
\mathbb{E}_{(m,0)}[\tau_{\mathcal S_\epsilon}]
=
[1+o_\epsilon(1)]
\frac{2\pi}{\mu^\sigma}
\sqrt{\frac{-\det \mathbb H_U^\sigma}{\det \mathbb H_U^m}}
\exp\!\left(\frac{U(\sigma)-U(m)}{\epsilon}\right),
\qquad \epsilon\to 0,
$$
with
$$
\mu^\sigma=\frac{-\gamma+\sqrt{\gamma^2+4\lambda^\sigma}}{2}>0.
$$
The derivation circumvents traditional potential-theoretic tools, which are ill-defined for this hypoelliptic process [2503.12610].

Outside equilibrium theory, underdamped simulations have been used as direct computational tools for diffusion in layered or discontinuous media. In a one-dimensional two-layer model of a drug-eluting stent, the motion is simulated with the underdamped Langevin equation
$$
m\frac{dv}{dt}=-\alpha(x)v+\beta(x(t)),
$$
using the G-JF integrator, the inertial convention for spatially varying friction, and a fictitious-mass device to enlarge the ballistic time without altering the long-time diffusive solution [1805.06697]. A semi-permeable membrane at the interface is implemented by an exponential waiting time
$$
p(t)=\frac{1}{T}e^{-t/T},
$$
and the reported finding is that the membrane has virtually no delay effect on the overall rate of delivery for realistic parameter values [1805.06697].

Underdamped Langevin dynamics also underlies kinetic reaction-diffusion limits. In particle-based reactive Langevin dynamics, single-particle motion is governed by
$$
dQ_t=V_t\,dt,\qquad dV_t=-b_A V_t\,dt+b_A\sqrt{2D_A}\,dW_t,
$$
while reactions are driven by stochastic kernels and placement densities on phase space [2606.03079]. In the large-population limit, the empirical phase-space measures converge to deterministic nonlocal kinetic reaction-diffusion partial integro-differential equations with hypoelliptic transport operators and reaction terms that retain both spatial and velocity structure [2606.03079]. Relative to overdamped PBSRD models, the limiting equations are kinetic rather than purely configurational, and they preserve momentum-level information absent from position-only closures [2606.03079].

A final connection runs toward stochastic optimization. For stochastic gradient descent with momentum, quantitative error estimates to the underdamped Langevin diffusion are available in both \(1\)-Wasserstein and total variation distance [2410.17297]. Under the stated smoothness, moment, and learning-rate assumptions, the discrepancies scale as \(\sqrt{n_n}\) in \(W_1\) and \(n_n^{1/4}\) in total variation, with an additional \(1/\sqrt N\) term reflecting minibatch noise [2410.17297]. This suggests that underdamped Langevin diffusion functions not only as a sampling process and a physical transport model, but also as a mathematically controlled continuum surrogate for momentum-based stochastic optimization.

Source: https://www.emergentmind.com/topics/underdamped-langevin-diffusion