Langevin dynamics with general kinetic energies

Gabriel Stoltz 1, 2 Zofia Trstanova 3
Inria de Paris, CERMICS - Centre d'Enseignement et de Recherche en Mathématiques et Calcul Scientifique
3 NANO-D - Algorithms for Modeling and Simulation of Nanosystems
Inria Grenoble - Rhône-Alpes, LJK - Laboratoire Jean Kuntzmann, INPG - Institut National Polytechnique de Grenoble
Abstract : We study Langevin dynamics with a kinetic energy different from the standard, quadratic one in order to accelerate the sampling of Boltzmann-Gibbs distributions. In particular, this kinetic energy can be nonglobally Lipschitz, which raises issues for the stability of discretizations of the associated Langevin dynamics. We first prove the exponential convergence of the law of the continuous process to the Boltzmann-Gibbs measure by a hypocoercive approach and characterize the asymptotic variance of empirical averages over trajectories. We next develop numerical schemes which are stable and of weak order two by considering splitting strategies where the discretizations of the fluctuation/dissipation are corrected by a Metropolis procedure. We use the newly developed schemes for two applications: optimizing the shape of the kinetic energy for the so-called adaptively restrained Langevin dynamics (which considers perturbations of standard quadratic kinetic energies vanishing around the origin) and reducing the metastability of some toy models using nonglobally Lipschitz kinetic energies.
Complete list of metadatas
Contributor : Gabriel Stoltz <>
Submitted on : Monday, September 12, 2016 - 9:35:40 PM
Last modification on : Friday, April 19, 2019 - 11:38:18 AM

Links full text



Gabriel Stoltz, Zofia Trstanova. Langevin dynamics with general kinetic energies. Multiscale Modeling and Simulation: A SIAM Interdisciplinary Journal, Society for Industrial and Applied Mathematics, 2018, 16 (2), pp.777-806. ⟨10.1137/16M110575X⟩. ⟨hal-01364821⟩



Record views