Optimal Displacement Parameters in Monte Carlo Simulations.

An adaptive algorithm optimizing single-particle translational displacement parameters in Metropolis Monte Carlo simulations is presented. The optimization is based on maximizing the mean square displacement of a trial move. It is shown that a large mean square displacement is strongly correlated with a high precision of average potential energy. The method is here demonstrated on model systems representing a Lennard-Jones fluid and a dilute polymer solution at poor solvent conditions. Our adaptive algorithm removes the need to provide values of displacement parameters in simulations, and it is easily extendable to optimize parameters of other types of trial moves.

[1]  P. Linse.  Interaction between colloids with grafted diblock polyampholytes. , 2007, The Journal of chemical physics.

[2]  P. Linse.  On the convergence of simulation of asymmetric electrolytes with charge asymmetry 60:1 , 1999 .

[3]  Rescic Jurij,et al.  MOLSIM: A modular molecular simulation software , 2015, J. Comput. Chem..

[4]  Raymond D. Mountain,et al.  Quantative measure of efficiency of Monte Carlo simulations , 1994 .

[5]  Robert H. Swendsen.  How the maximum step size in Monte Carlo simulations should be adjusted , 2011 .

[6]  Kumar,et al.  Efficient Monte Carlo methods for the computer simulation of biological molecules. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[7]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[8]  Optimum Monte Carlo simulations: some exact results , 2003, cond-mat/0302366.

[9]  M. Rao,et al.  On the force bias Monte Carlo simulation of water: methodology, optimization and comparison with molecular dynamics , 1979 .

[10]  Andrew J Schultz,et al.  Quantifying Computational Effort Required for Stochastic Averages. , 2014, Journal of chemical theory and computation.

[11]  S. Schneider,et al.  Effect of the Molecular Architecture on the Internal Complexation Behavior of Linear Copolymers and Miktoarm Star Polymers , 2015 .

[12]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[13]  J. Mccammon,et al.  Simulation methods for protein structure fluctuations , 1980, Biopolymers.

[14]  P. Linse,et al.  Diblock polyampholytes grafted onto spherical particles: Monte Carlo simulation and lattice mean-field theory. , 2004, Langmuir : the ACS journal of surfaces and colloids.

[15]  Mark A. Miller,et al.  Should one adjust the maximum step size in a Metropolis Monte Carlo simulation , 2000 .

[16]  P. Linse.  Simulation of charged colloids in solution , 2005 .