skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics

Journal Article · · Journal of Chemical Physics
DOI:https://doi.org/10.1063/1.4996204· OSTI ID:1497877
ORCiD logo [1]; ORCiD logo [1];  [1];  [2]; ORCiD logo [1]; ORCiD logo [1]; ORCiD logo [1];  [2]
  1. Peking University, Beijing (China). Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering
  2. Beijing Normal University, Beijing (China). College of Chemistry and Center for Advanced Quantum Studies, Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education

Langevin dynamics has become a popular tool to simulate the Boltzmann equilibrium distribution. When the repartition of the Langevin equation involves the exact realization of the Ornstein-Uhlenbeck noise, in addition to the conventional density evolution, there exists another type of discrete evolution that may not correspond to a continuous, real dynamical counterpart. This virtual dynamics case is also able to produce the desired stationary distribution. Different types of repartition lead to different numerical schemes, of which the accuracy and efficiency are investigated through studying the harmonic oscillator potential, an analytical solvable model. By analyzing the asymptotic distribution and characteristic correlation time that are derived by either directly solving the discrete equations of motion or using the related phase space propagators, it is shown that the optimal friction coefficient resulting in the minimum characteristic correlation time depends on the time interval chosen in the numerical implementation. In conclusion, when the recommended “middle” scheme is employed, both analytical and numerical results demonstrate that, for good numerical performance in efficiency as well as accuracy, one may choose a friction coefficient in a wide range from around the optimal value to the high friction limit.

Research Organization:
Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC)
Sponsoring Organization:
USDOE Office of Science (SC)
Grant/Contract Number:
AC02-05CH11231
OSTI ID:
1497877
Journal Information:
Journal of Chemical Physics, Vol. 147, Issue 18; ISSN 0021-9606
Publisher:
American Institute of Physics (AIP)Copyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 20 works
Citation information provided by
Web of Science

References (39)

Algorithms for brownian dynamics journal February 1982
A fourth-order numerical integrator for stochastic Langevin equations journal June 1998
A unified thermostat scheme for efficient configurational sampling for classical/quantum canonical ensembles via molecular dynamics journal July 2017
Robust and efficient configurational molecular sampling via Langevin dynamics journal May 2013
Symplectic integrators from composite operator factorizations journal March 1997
Brownian dynamics with hydrodynamic interactions journal August 1978
Path integrals for Fokker–Planck dynamics with singular diffusion: Accurate factorization for the time evolution operator journal August 1998
Stochastic boundary conditions for molecular dynamics simulations of ST2 water journal March 1984
Rational Construction of Stochastic Numerical Methods for Molecular Sampling journal June 2012
A Leap-frog Algorithm for Stochastic Dynamics journal March 1988
Quasilinear approximations for the propagator of the Fokker-Planck equation journal June 1993
Efficient stochastic thermostatting of path integral molecular dynamics journal September 2010
Structure of positive decompositions of exponential operators journal January 2005
Algorithms for Brownian dynamics journal June 2003
Hybrid exponential product formulas for unbounded operators with possible applications to Monte Carlo simulations journal June 1995
A simple and effective Verlet-type algorithm for simulating Langevin dynamics journal January 2013
Long-Run Accuracy of Variational Integrators in the Stochastic Context journal January 2010
Construction of higher order symplectic integrators journal November 1990
Path integral Liouville dynamics: Applications to infrared spectra of OH, water, ammonia, and methane journal January 2016
Decomposition formulas of exponential operators and Lie exponentials with some applications to quantum mechanics and statistical physics journal April 1985
The computation of averages from equilibrium and nonequilibrium Langevin molecular dynamics journal January 2015
General theory of fractal path integrals with applications to many‐body theories and statistical physics journal February 1991
The determination of an accurate isotope dependent potential energy surface for water from extensive ab initio calculations and experimental data journal March 1997
Linearized semiclassical initial value time correlation functions using the thermal Gaussian approximation: Applications to condensed phase systems journal September 2007
Numerical integration of the Langevin equation: Monte Carlo simulation journal April 1980
Design of quasisymplectic propagators for Langevin dynamics journal July 2007
A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat journal July 2016
Accurate sampling using Langevin dynamics journal May 2007
Langevin stabilization of molecular dynamics journal February 2001
Monte Carlo Calculation of Quantum Systems. II. Higher Order Correction journal November 1984
Exact numerical simulation of the Ornstein-Uhlenbeck process and its integral journal August 1996
High order three part split symplectic integrators: Efficient techniques for the long time simulation of the disordered discrete nonlinear Schrödinger equation journal May 2014
Accurate sampling using Langevin dynamics text January 2008
Efficient stochastic thermostatting of path integral molecular dynamics text January 2010
A simple and effective Verlet-type algorithm for simulating Langevin dynamics text January 2012
High order three part split symplectic integrators: Efficient techniques for the long time simulation of the disordered discrete nonlinear Schroedinger equation text January 2013
Robust and efficient configurational molecular sampling via Langevin Dynamics text January 2013
The computation of averages from equilibrium and nonequilibrium Langevin molecular dynamics text January 2013
A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat text January 2016

Cited By (6)

Sampling the thermal Wigner density via a generalized Langevin dynamics journal September 2019
Dimension-free path-integral molecular dynamics without preconditioning journal March 2020
Critical role of quantum dynamical effects in the Raman spectroscopy of liquid water journal November 2017
Complete set of stochastic Verlet-type thermostats for correct Langevin simulations journal September 2019
Critical role of quantum dynamical effects in the Raman spectroscopy of liquid water text January 2018
Critical role of quantum dynamical effects in the Raman spectroscopy of liquid water text January 2017

Similar Records

A generalized class of strongly stable and dimension-free T-RPMD integrators
Journal Article · Fri Jan 08 00:00:00 EST 2021 · Journal of Chemical Physics · OSTI ID:1497877

Noise and error analysis and optimization in particle-based kinetic plasma simulations
Journal Article · Wed May 05 00:00:00 EDT 2021 · Journal of Computational Physics · OSTI ID:1497877

Formulation of a moment method for multidimensional Fokker-Planck equations
Journal Article · Thu Jun 01 00:00:00 EDT 1995 · Physical Review. E, Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics · OSTI ID:1497877