Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics
Abstract
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.
- Authors:
-
- Peking University, Beijing (China). Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering
- Beijing Normal University, Beijing (China). College of Chemistry and Center for Advanced Quantum Studies, Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education
- Publication Date:
- Research Org.:
- Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC)
- Sponsoring Org.:
- USDOE Office of Science (SC)
- OSTI Identifier:
- 1497877
- Grant/Contract Number:
- AC02-05CH11231
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Journal of Chemical Physics
- Additional Journal Information:
- Journal Volume: 147; Journal Issue: 18; Journal ID: ISSN 0021-9606
- Publisher:
- American Institute of Physics (AIP)
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY
Citation Formats
Li, Dezhang, Han, Xu, Chai, Yichen, Wang, Cong, Zhang, Zhijun, Chen, Zifei, Liu, Jian, and Shao, Jiushu. Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics. United States: N. p., 2017.
Web. doi:10.1063/1.4996204.
Li, Dezhang, Han, Xu, Chai, Yichen, Wang, Cong, Zhang, Zhijun, Chen, Zifei, Liu, Jian, & Shao, Jiushu. Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics. United States. https://doi.org/10.1063/1.4996204
Li, Dezhang, Han, Xu, Chai, Yichen, Wang, Cong, Zhang, Zhijun, Chen, Zifei, Liu, Jian, and Shao, Jiushu. Thu .
"Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics". United States. https://doi.org/10.1063/1.4996204. https://www.osti.gov/servlets/purl/1497877.
@article{osti_1497877,
title = {Stationary state distribution and efficiency analysis of the Langevin equation via real or virtual dynamics},
author = {Li, Dezhang and Han, Xu and Chai, Yichen and Wang, Cong and Zhang, Zhijun and Chen, Zifei and Liu, Jian and Shao, Jiushu},
abstractNote = {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.},
doi = {10.1063/1.4996204},
journal = {Journal of Chemical Physics},
number = 18,
volume = 147,
place = {United States},
year = {Thu Nov 09 00:00:00 EST 2017},
month = {Thu Nov 09 00:00:00 EST 2017}
}
Web of Science
Works referenced in this record:
Algorithms for brownian dynamics
journal, February 1982
- van Gunsteren, W. F.; Berendsen, H. J. C.
- Molecular Physics, Vol. 45, Issue 3
A fourth-order numerical integrator for stochastic Langevin equations
journal, June 1998
- Hershkovitz, Eli
- The Journal of Chemical Physics, Vol. 108, Issue 22
A unified thermostat scheme for efficient configurational sampling for classical/quantum canonical ensembles via molecular dynamics
journal, July 2017
- Zhang, Zhijun; Liu, Xinzijian; Chen, Zifei
- The Journal of Chemical Physics, Vol. 147, Issue 3
Robust and efficient configurational molecular sampling via Langevin dynamics
journal, May 2013
- Leimkuhler, Benedict; Matthews, Charles
- The Journal of Chemical Physics, Vol. 138, Issue 17
Symplectic integrators from composite operator factorizations
journal, March 1997
- Chin, Siu A.
- Physics Letters A, Vol. 226, Issue 6
Brownian dynamics with hydrodynamic interactions
journal, August 1978
- Ermak, Donald L.; McCammon, J. A.
- The Journal of Chemical Physics, Vol. 69, Issue 4
Path integrals for Fokker–Planck dynamics with singular diffusion: Accurate factorization for the time evolution operator
journal, August 1998
- Drozdov, Alexander N.; Talkner, Peter
- The Journal of Chemical Physics, Vol. 109, Issue 6
Stochastic boundary conditions for molecular dynamics simulations of ST2 water
journal, March 1984
- Brünger, Axel; Brooks, Charles L.; Karplus, Martin
- Chemical Physics Letters, Vol. 105, Issue 5
Rational Construction of Stochastic Numerical Methods for Molecular Sampling
journal, June 2012
- Leimkuhler, B.; Matthews, C.
- Applied Mathematics Research eXpress
A Leap-frog Algorithm for Stochastic Dynamics
journal, March 1988
- Van Gunsteren, W. F.; Berendsen, H. J. C.
- Molecular Simulation, Vol. 1, Issue 3
Quasilinear approximations for the propagator of the Fokker-Planck equation
journal, June 1993
- Drozdov, Alexander N.
- Physica A: Statistical Mechanics and its Applications, Vol. 196, Issue 2
Efficient stochastic thermostatting of path integral molecular dynamics
journal, September 2010
- Ceriotti, Michele; Parrinello, Michele; Markland, Thomas E.
- The Journal of Chemical Physics, Vol. 133, Issue 12
Structure of positive decompositions of exponential operators
journal, January 2005
- Chin, Siu A.
- Physical Review E, Vol. 71, Issue 1
Algorithms for Brownian dynamics
journal, June 2003
- Ricci, Andrea; Ciccotti, Giovanni
- Molecular Physics, Vol. 101, Issue 12
Hybrid exponential product formulas for unbounded operators with possible applications to Monte Carlo simulations
journal, June 1995
- Suzuki, Masuo
- Physics Letters A, Vol. 201, Issue 5-6
A simple and effective Verlet-type algorithm for simulating Langevin dynamics
journal, January 2013
- Grønbech-Jensen, Niels; Farago, Oded
- Molecular Physics, Vol. 111, Issue 8
Long-Run Accuracy of Variational Integrators in the Stochastic Context
journal, January 2010
- Bou-Rabee, Nawaf; Owhadi, Houman
- SIAM Journal on Numerical Analysis, Vol. 48, Issue 1
Construction of higher order symplectic integrators
journal, November 1990
- Yoshida, Haruo
- Physics Letters A, Vol. 150, Issue 5-7
Path integral Liouville dynamics: Applications to infrared spectra of OH, water, ammonia, and methane
journal, January 2016
- Liu, Jian; Zhang, Zhijun
- The Journal of Chemical Physics, Vol. 144, Issue 3
Decomposition formulas of exponential operators and Lie exponentials with some applications to quantum mechanics and statistical physics
journal, April 1985
- Suzuki, Masuo
- Journal of Mathematical Physics, Vol. 26, Issue 4
The computation of averages from equilibrium and nonequilibrium Langevin molecular dynamics
journal, January 2015
- Leimkuhler, Benedict; Matthews, Charles; Stoltz, Gabriel
- IMA Journal of Numerical Analysis
General theory of fractal path integrals with applications to many‐body theories and statistical physics
journal, February 1991
- Suzuki, Masuo
- Journal of Mathematical Physics, Vol. 32, Issue 2
The determination of an accurate isotope dependent potential energy surface for water from extensive ab initio calculations and experimental data
journal, March 1997
- Partridge, Harry; Schwenke, David W.
- The Journal of Chemical Physics, Vol. 106, Issue 11
Linearized semiclassical initial value time correlation functions using the thermal Gaussian approximation: Applications to condensed phase systems
journal, September 2007
- Liu, Jian; Miller, William H.
- The Journal of Chemical Physics, Vol. 127, Issue 11
Numerical integration of the Langevin equation: Monte Carlo simulation
journal, April 1980
- Ermak, Donald L.; Buckholz, Helen
- Journal of Computational Physics, Vol. 35, Issue 2
Design of quasisymplectic propagators for Langevin dynamics
journal, July 2007
- Melchionna, Simone
- The Journal of Chemical Physics, Vol. 127, Issue 4
A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat
journal, July 2016
- Liu, Jian; Li, Dezhang; Liu, Xinzijian
- The Journal of Chemical Physics, Vol. 145, Issue 2
Accurate sampling using Langevin dynamics
journal, May 2007
- Bussi, Giovanni; Parrinello, Michele
- Physical Review E, Vol. 75, Issue 5
Langevin stabilization of molecular dynamics
journal, February 2001
- Izaguirre, Jesús A.; Catarello, Daniel P.; Wozniak, Justin M.
- The Journal of Chemical Physics, Vol. 114, Issue 5
Monte Carlo Calculation of Quantum Systems. II. Higher Order Correction
journal, November 1984
- Takahashi, Minoru; Imada, Masatoshi
- Journal of the Physical Society of Japan, Vol. 53, Issue 11
Exact numerical simulation of the Ornstein-Uhlenbeck process and its integral
journal, August 1996
- Gillespie, Daniel T.
- Physical Review E, Vol. 54, Issue 2
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
- Skokos, Ch.; Gerlach, E.; Bodyfelt, J. D.
- Physics Letters A, Vol. 378, Issue 26-27
Accurate sampling using Langevin dynamics
text, January 2008
- Bussi, Giovanni; Parrinello, Michele
- arXiv
Efficient stochastic thermostatting of path integral molecular dynamics
text, January 2010
- Ceriotti, Michele; Parrinello, Michele; Markland, Thomas E.
- arXiv
A simple and effective Verlet-type algorithm for simulating Langevin dynamics
text, January 2012
- Grønbech-Jensen, Niels; Farago, Oded
- arXiv
High order three part split symplectic integrators: Efficient techniques for the long time simulation of the disordered discrete nonlinear Schroedinger equation
text, January 2013
- Skokos, Ch.; Gerlach, E.; Bodyfelt, J. D.
- arXiv
Robust and efficient configurational molecular sampling via Langevin Dynamics
text, January 2013
- Leimkuhler, Benedict; Matthews, Charles
- arXiv
The computation of averages from equilibrium and nonequilibrium Langevin molecular dynamics
text, January 2013
- Leimkuhler, Benedict; Matthews, Charles; Stoltz, Gabriel
- arXiv
A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat
text, January 2016
- Liu, Jian; Li, Dezhang; Liu, Xinzijian
- arXiv
Works referencing / citing this record:
Sampling the thermal Wigner density via a generalized Langevin dynamics
journal, September 2019
- Plé, Thomas; Huppert, Simon; Finocchi, Fabio
- The Journal of Chemical Physics, Vol. 151, Issue 11
Dimension-free path-integral molecular dynamics without preconditioning
journal, March 2020
- Korol, Roman; Rosa-Raíces, Jorge L.; Bou-Rabee, Nawaf
- The Journal of Chemical Physics, Vol. 152, Issue 10
Critical role of quantum dynamical effects in the Raman spectroscopy of liquid water
journal, November 2017
- Liu, Xinzijian; Liu, Jian
- Molecular Physics, Vol. 116, Issue 7-8
Complete set of stochastic Verlet-type thermostats for correct Langevin simulations
journal, September 2019
- Grønbech-Jensen, Niels
- Molecular Physics, Vol. 118, Issue 8
Critical role of quantum dynamical effects in the Raman spectroscopy of liquid water
text, January 2018
- Liu, Xinzijian; Liu, Jian
- Taylor & Francis
Critical role of quantum dynamical effects in the Raman spectroscopy of liquid water
text, January 2017
- Liu, Xinzijian; Liu, Jian
- arXiv