Kinetic Langevin Diffusion

• Kinetic Langevin Diffusion is a family of stochastic processes that integrates momentum to enhance sampling efficiency compared to overdamped methods.
• Discrete-time implementations like the GJF–2GJ thermostat yield robust configurational and kinetic statistics across a range of physical systems.
• The method’s leap-frog updates and companion velocity ensure accurate thermodynamic averages and time-step independent diffusion metrics.

Kinetic Langevin Diffusion is a family of stochastic dynamical systems that augment classical overdamped (first-order) Langevin dynamics with velocity (momentum-like) variables, yielding a second-order or underdamped diffusion process. This framework is widely used in computational statistical mechanics, MCMC sampling, and molecular dynamics, and is the basis for a number of high-performance algorithms and analytical tools. The kinetic Langevin approach enables improved mixing rates, dimensionally favorable convergence properties, time-step robust discretizations, and more accurate configurational and kinetic statistics compared to classical overdamped schemes. Discrete-time kinetic Langevin methods, such as the Grønbech-Jensen–Farago (GJF) thermostat, yield robust samplers for both equilibrium and transport statistics over a wide class of physical and probabilistic systems (Jensen et al., 2019).

1. Mathematical Formulation of Kinetic Langevin Diffusion

Kinetic Langevin diffusion describes the evolution of the phase-space vector $(r, v)$ (position and velocity) as a solution to the coupled stochastic differential equations: $$\begin{cases} dr_t = v_t\,dt, \ m\,dv_t = f(r_t)\,dt - \alpha\,v_t\,dt + \sqrt{2\alpha k_B T}\,dW_t, \end{cases}$$ where $m$ is mass, $\alpha$ the friction coefficient, $f(r_t)$ a deterministic force field (often $-\nabla V(r_t)$), and $W_t$ standard Brownian motion. The system admits the canonical Boltzmann–Gibbs invariant measure: $$\mu(dr\,dv) = Z^{-1}\exp\left(-\beta[V(r)+(1/2)m|v|^2]\right) dr\,dv,$$ where $\beta = (k_B T)^{-1}$. Assuming sufficient regularity, this SDE enjoys exponential ergodicity and hypocoercivity properties under mild convexity and Poincaré conditions (Stoltz et al., 2016).

2. Discrete-Time Kinetic Langevin: GJF Thermostat and Companion Velocity

Discrete-time schemes are essential for computational implementation. The GJF thermostat (Jensen et al., 2019) is a discrete integrator for Langevin dynamics that maintains robust configurational and kinetic sampling:

Update scheme (with time step $dt$):

$$\begin{aligned} & r^{n+1} = r^n + b \left[ dt\,v^n + \frac{dt^2}{2m} f^n + \frac{dt}{2m} \beta^{n+1} \right], \ & v^{n+1} = a\,v^n + \frac{dt}{2m}(a\,f^n + f^{n+1}) + \frac{b}{m}\, \beta^{n+1} \end{aligned}$$

where $a = \frac{1 - \alpha dt/(2m)}{1 + \alpha dt/(2m)}$, $b = \frac{1}{1 + \alpha dt/(2m)}$, and $\beta^{n+1}$ is Gaussian noise.

Companion half-step velocity:

$$u^{n+\frac{1}{2}} = \frac{r^{n+1} - r^n}{\sqrt{b} dt}$$

This companion velocity yields time-step-independent kinetic statistics, exactly reproducing the Maxwell–Boltzmann average kinetic energy even for large $dt$. The method supports a leap-frog formulation efficiently updating both position and velocity (Jensen et al., 2019).

3. Exact Kinetic and Configurational Statistics, Green–Kubo Diffusion

A hallmark of the GJF–2GJ (GJF with companion velocity) scheme is its statistical accuracy in both configurational and kinetic observables:

• $\langle (u^{n+\frac{1}{2}})^2 \rangle = k_B T / m$ is exact for all $dt \leq 2/\Omega_0$ in harmonic systems.
• The scheme yields exact time-step-independent variance of kinetic energy and correct discrete-time autocorrelations.
• The discrete-time Green–Kubo formula evaluated via the autocorrelation function of the companion velocity yields:

$D = \frac{k_B T}{\alpha}$

independently of $dt$, provided stability constraints hold (Jensen et al., 2019).

4. Extension to Nonlinear and Many-Body Systems

Empirical studies confirm that GJF–2GJ maintains statistical accuracy beyond the harmonic case:

• For nonlinear oscillators (non-quadratic potentials), the method produces kinetic and configurational temperatures and fluctuations accurate to within 1% up to the stability threshold.
• In three-dimensional molecular dynamics with hundreds of particles (Lennard-Jones solids and liquids), the method preserves correct phase-space and thermodynamic observables as $dt$ approaches the maximum stable step, outperforming classical integrators such as BBK and PBS which exhibit systematic drift under similar conditions (Jensen et al., 2019).

5. Algorithmic Summary: Leap-Frog and Stability Constraints

The scheme is efficiently implemented using a leap-frog update:

Step	Update Equation	Notes
Half-step velocity update	$$u^{n+\frac{1}{2}} = a u^{n-\frac{1}{2}} + (\sqrt{b} dt/m) f(r^n) + (\sqrt{b}/2m)(\beta^n+\beta^{n+1})$$
Position update	$r^{n+1} = r^n + \sqrt{b} dt u^{n+\frac{1}{2}}$
Force evaluation	$f(r^n) = -\partial V/\partial r |_{r^n}$

The stability domain is governed by $\Omega_{max} dt < 2$ and $\alpha dt/(2m) < 1$, where $\Omega_{max}$ is the highest oscillator frequency. For maximal efficiency, $dt$ should approach but not reach the stability boundary.

6. Comparison with Continuous-time Langevin and Other Discretizations

Continuous-time Langevin yields exact equipartition at any moment. However, naive discretizations introduce $dt^2$ errors in either configuration or kinetic averages. The GJF–2GJ method, by contrast, ensures exact statistics (energies, fluctuations, diffusion) up to the method's stability limit and aligns with continuous dynamics in the zero time-step limit.

In practice:

• BBK, PBS, and conventional splitting schemes show systematic time-step-dependent deviations in kinetic or configurational averages.
• GJF–2GJ exhibits weak second-order accuracy (error $O(dt^2)$) in generic systems, and error-free sampling in the harmonic limit (Finkelstein et al., 2021).

7. Practical Recommendations and Broader Context

For effective simulation:

• Choose $dt$ to be 5–10% below the strict analytic stability bound.
• Friction $\gamma$ ($\alpha/m$) should be commensurate with the slowest relevant timescale; moderate values avoid overdamping while ensuring good sampling.
• For precise diffusion estimates, use Riemann-summed velocity autocorrelations.
• The method's validity and robustness extend to complex, nonlinear, and many-body systems without additional computational cost.

Kinetic Langevin discretizations, particularly in the GJF–2GJ class, provide a practical, thermodynamically sound standard for phase-space sampling, molecular dynamics, and stochastic optimization where accurate kinetic, configurational, and transport properties are required (Jensen et al., 2019).