Comparison of Modern Langevin Integrators for Simulations of Coarse-Grained Polymer Melts

Abstract

For a wide range of phenomena, current computational ability does not always allow for fully atomistic simulations of high-dimensional molecular systems to reach time scales of interest. Coarse-graining (CG) is an established approach to alleviate the impact of computational limits while retaining the same algorithms used in atomistic simulations. It is of importance to understand how algorithms such as Langevin integrators perform on non-trivial CG molecular systems, and in particular how large of an integration time step can be used without introducing unacceptable amounts of error into averaged quantities of interest. To investigate this, we examined three different Langevin integrators on a CG polymer melt: the recently developed BAOAB method by Leimkuhler and Matthews, the Gronbech-Jensen and Farago method, or G-JF, and the frequently used Brunger-Brooks-Karplus integrator, also known as BBK. We compute and analyze key statistical properties for each. Our results indicate that the three integrators perform similarly when using a small friction parameter; however, outside of this regime the use of large integration steps produces significant deviations from the predicted diffusivity and steady-state distributions for all integration methods examined with the exception of G-JF.

20 pages, 7 figures

References in corpus (4)

• Langevin thermostat for robust configurational and kinetic sampling
• Accurate Configurational and Kinetic Statistics in Discrete-Time Langevin Systems
• Assessing numerical methods for molecular and particle simulation
• General order conditions for stochastic partitioned Runge-Kutta methods

Cited by in corpus (9)

• Algorithms for Brownian dynamics across discontinuities
• On the Application of Non-Gaussian Noise in Stochastic Langevin Simulations
• The Challenge of Stochastic Størmer-Verlet Thermostats Generating Correct Statistics
• A Langevin dynamics approach for multi-layer mass transfer problems
• Defining velocities for accurate kinetic statistics in the GJF thermostat
• Bringing discrete-time Langevin splitting methods into agreement with thermodynamics
• On the Definition of Velocity in Discrete-Time, Stochastic Langevin Simulations
• Linear Analysis of Stochastic Verlet-Type Integrators for Langevin Equations
• Self-averaging parameter estimation for coarse-grained particle models