Histogram-free multicanonical Monte Carlo sampling to calculate the density of states

Abstract We report a new multicanonical Monte Carlo algorithm to obtain the density of states for physical systems with continuous state variables in statistical mechanics. Our algorithm is able to obtain a closed-form expression for the density of states expressed in a chosen basis set, instead of a numerical array of finite resolution as in previous variants of this class of MC methods such as the multicanonical sampling and Wang–Landau sampling. This is enabled by storing the visited states directly and avoiding the explicit collection of a histogram. This practice also has the advantage of avoiding undesirable artificial errors caused by the discretization and binning of continuous state variables. Our results show that this scheme is capable of obtaining converged results with a much reduced number of Monte Carlo steps, leading to a significant speedup over existing algorithms. Program summary Program Title: HistogramFreeMUCA Program Files doi: http://dx.doi.org/10.17632/mtgyvsfzyg.1 Licensing provisions: BSD 3-clause Programming language: Python Nature of problem: This program implements a novel algorithm to obtain the density of states of a physical system, expanded in a chosen basis set. Unlike existing algorithms that return the density of states as a numerical array, this algorithm avoids binning of a continuous variable and is able to express the density of states as a closed-form expression. It is thus suitable for the study of the statistical mechanics and thermodynamic properties of physical systems where the density of a continuous state variable is of interest. Solution method: The new algorithm presented here is a special reweighting method in classical Monte Carlo approaches. In particular, it can be regarded as a descendant and a hybrid method closely related to the multicanonical method and Wang–Landau sampling. Additional comments: Most updated source code can be found at: https://github.com/yingwaili/HistogramFreeMUCA

[1]  D. Landau,et al.  Efficient, multiple-range random walk algorithm to calculate the density of states. , 2000, Physical review letters.

[2]  Thomas Vogel,et al.  Scalable replica-exchange framework for Wang-Landau sampling. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[3]  David P. Landau,et al.  Massively parallel Wang-Landau sampling on multiple GPUs , 2012, Comput. Phys. Commun..

[4]  Bernd A. Berg,et al.  From data to probability densities without histograms , 2007, Comput. Phys. Commun..

[5]  Bernd A. Berg.  Markov Chain Monte Carlo Simulations and Their Statistical Analysis , 2004 .

[6]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[7]  D. V. Lindley,et al.  Handbook of Statistical Tables. , 1962 .

[8]  P. Hohenberg,et al.  Theory of Dynamic Critical Phenomena , 1977 .

[9]  K. Binder,et al.  A Guide to Monte Carlo Simulations in Statistical Physics: Preface , 2005 .

[10]  Proceedings of the Platform for Advanced Scientific Computing Conference , 2018, PASC.

[11]  D. Landau,et al.  Determining the density of states for classical statistical models: a random walk algorithm to produce a flat histogram. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[12]  Johannes Zierenberg,et al.  Scaling properties of a parallel implementation of the multicanonical algorithm , 2013, Comput. Phys. Commun..

[13]  Johannes Zierenberg,et al.  Massively parallel multicanonical simulations , 2017, Comput. Phys. Commun..

[14]  Markus Eisenbach,et al.  Density-functional Monte-Carlo simulation of CuZn order-disorder transition , 2015, 1510.01543.

[15]  Thomas Vogel,et al.  Generic, hierarchical framework for massively parallel Wang-Landau sampling. , 2013, Physical review letters.

[16]  B. Berg,et al.  Multicanonical algorithms for first order phase transitions , 1991 .

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

[18]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[19]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[20]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[21]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[22]  David P. Landau,et al.  Numerical integration using Wang-Landau sampling , 2007, Comput. Phys. Commun..