Long-Run Accuracy of Variational Integrators in the Stochastic Context
arXiv:0712.4123 · doi:10.1137/090758842
Abstract
This paper presents a Lie-Trotter splitting for inertial Langevin equations (Geometric Langevin Algorithm) and analyzes its long-time statistical properties. The splitting is defined as a composition of a variational integrator with an Ornstein-Uhlenbeck flow. Assuming the exact solution and the splitting are geometrically ergodic, the paper proves the discrete invariant measure of the splitting approximates the invariant measure of inertial Langevin to within the accuracy of the variational integrator in representing the Hamiltonian. In particular, if the variational integrator admits no energy error, then the method samples the invariant measure of inertial Langevin without error. Numerical validation is provided using explicit variational integrators with first, second, and fourth order accuracy.
30 pages
References in corpus (5)
Cited by in corpus (35)
- Robust and efficient configurational molecular sampling via Langevin Dynamics
- Rational Construction of Stochastic Numerical Methods for Molecular Sampling
- The computation of averages from equilibrium and nonequilibrium Langevin molecular dynamics
- Explicit symplectic approximation of nonseparable Hamiltonians: algorithm and long time performance
- Randomized Hamiltonian Monte Carlo
- On the numerical treatment of dissipative particle dynamics and related systems
- Non-intrusive and structure preserving multiscale integration of stiff ODEs, SDEs and Hamiltonian systems with hidden slow dynamics via flow averaging
- Time step rescaling recovers continuous-time dynamical properties for discrete-time Langevin integration of nonequilibrium systems
- Dimension-free path-integral molecular dynamics without preconditioning
- Using nonequilibrium fluctuation theorems to understand and correct errors in equilibrium and nonequilibrium discrete Langevin dynamics simulations
- Stochastic Discrete Hamiltonian Variational Integrators
- A unified operator splitting approach for multi-scale fluid-particle coupling in the lattice Boltzmann method
- High-order integrator for sampling the invariant distribution of a class of parabolic SPDEs with additive space-time noise
- Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics
- Sampling the isothermal-isobaric ensemble by Langevin dynamics
- Pairwise adaptive thermostats for improved accuracy and stability in dissipative particle dynamics
- Convergence of unadjusted Hamiltonian Monte Carlo for mean-field models
- Order conditions for sampling the invariant measure of ergodic stochastic differential equations on manifolds
- Pressure-induced recovery of Fourier's law in one dimensional momentum-conserving systems
- Postprocessed integrators for the high order integration of ergodic SDEs
- Langevin dynamics with general kinetic energies
- Irreversible Samplers from Jump and Continuous Markov Processes
- Stable schemes for dissipative particle dynamics with conserved energy
- Stochastic variational principles for the collisional Vlasov-Maxwell and Vlasov-Poisson equations
- Variational integrators for the dynamics of thermo-elastic solids with finite speed thermal waves
- Parameter Estimation in Nonlinear Multivariate Stochastic Differential Equations Based on Splitting Schemes
- A symmetry-adapted numerical scheme for SDEs
- Variational integrators for interconnected Lagrange-Dirac systems
- The Hamiltonian Mechanics of Stochastic Acceleration
- A uniformly accurate scheme for the numerical integration of penalized Langevin dynamics
- Coupling a reactive potential with a harmonic approximation for atomistic simulations of material failure
- Strang Splitting for Parametric Inference in Second-order Stochastic Differential Equations
- Hopf algebra structures for the backward error analysis of ergodic stochastic differential equations
- GROMACS Stochastic Dynamics and BAOAB are equivalent configurational sampling algorithms
- NySALT: Nyström-type inference-based schemes adaptive to large time-stepping