Efficient Langevin and Monte Carlo sampling algorithms: The case of field-theoretic simulations
1School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln LN6 7TS, United Kingdom.
Abstract:
We introduce Langevin sampling algorithms to field-theoretic simulations (FTSs) of polymers that, for the same accuracy, are ∼10× more efficient than a previously used Brownian dynamics algorithm that used predictor corrector for such simulations, over 10× more efficient than the smart Monte Carlo (SMC) algorithm, and typically over 1000× more efficient than a simple Monte Carlo (MC) algorithm. These algorithms are known as the Leimkuhler-Matthews (the BAOAB-limited) method and the BAOAB method. Furthermore, the FTS allows for an improved MC algorithm based on the Ornstein-Uhlenbeck process (OU MC), which is 2× more efficient than SMC. The system-size dependence of the efficiency for the sampling algorithms is presented, and it is shown that the aforementioned MC algorithms do not scale well with system sizes. Hence, for larger sizes, the efficiency difference between the Langevin and MC algorithms is even greater, although, for SMC and OU MC, the scaling is less unfavorable than for the simple MC.
More Related Videos
10:52Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics
Published on: April 12, 2019
11:03An Analog Macroscopic Technique for Studying Molecular Hydrodynamic Processes in Dense Gases and Liquids
Published on: December 4, 2017
Related Concept Videos
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
Maxwell-Boltzmann Distribution: Problem Solving
This distribution function f(v) is defined by saying that the expected number N (v1,v2) of particles with speeds between v1 and v2 is given by
One-Compartment Open Model: Wagner-Nelson and Loo Riegelman Method for ka Estimation
On...
Sampling Theorem
Sampling Methods: Overview
In analytical chemistry, the choice of...
Linear Approximation in Time Domain
For a simple pendulum with a mass evenly distributed along its length and the center of mass located at half the pendulum's length,...
