Convergence of Stochastic Approximation Monte Carlo and modified Wang-Landau algorithms: Tests for the Ising model

Abstract We investigate the behavior of the deviation of the estimator for the density of states (DOS) with respect to the exact solution in the course of Wang–Landau and Stochastic Approximation Monte Carlo (SAMC) simulations of the two-dimensional Ising model. We find that the deviation saturates in the Wang–Landau case. This can be cured by adjusting the refinement scheme. To this end, the 1 ∕ t -modification of the Wang–Landau algorithm has been suggested. A similar choice of refinement scheme is employed in the SAMC algorithm. The convergence behavior of all three algorithms is examined. It turns out that the convergence of the SAMC algorithm is very sensitive to the onset of the refinement. Finally, the internal energy and specific heat of the Ising model are calculated from the SAMC DOS and compared to exact values.

[1]  R. Belardinelli,et al.  Fast algorithm to calculate density of states. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[2]  Wang-Landau algorithm: a theoretical analysis of the saturation of the error. , 2007, The Journal of chemical physics.

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

[4]  L. Onsager.  Crystal statistics. I. A two-dimensional model with an order-disorder transition , 1944 .

[5]  E. Ising.  Beitrag zur Theorie des Ferromagnetismus , 1925 .

[6]  Adam D. Swetnam,et al.  Improving the Wang–Landau algorithm for polymers and proteins , 2011, J. Comput. Chem..

[7]  Thomas Wüst,et al.  Versatile approach to access the low temperature thermodynamics of lattice polymers and proteins. , 2009, Physical review letters.

[8]  Berg,et al.  Multicanonical ensemble: A new approach to simulate first-order phase transitions. , 1992, Physical review letters.

[9]  D. Huse,et al.  Optimizing the ensemble for equilibration in broad-histogram Monte Carlo simulations. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  Wolfhard Janke,et al.  MULTICANONICAL SIMULATION OF THE TWO-DIMENSIONAL 7-STATE POTTS MODEL , 1992 .

[11]  Hesselbo,et al.  Monte Carlo simulation and global optimization without parameters. , 1995, Physical review letters.

[12]  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.

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

[14]  A Theory on Flat Histogram Monte Carlo Algorithms , 2006 .

[15]  Chenggang Zhou,et al.  Optimal modification factor and convergence of the Wang-Landau algorithm. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[17]  Eric Moulines,et al.  Stability of Stochastic Approximation under Verifiable Conditions , 2005, Proceedings of the 44th IEEE Conference on Decision and Control.

[18]  R. Carroll,et al.  Stochastic Approximation in Monte Carlo Computation , 2007 .

[19]  T. Shakirov,et al.  Stochastic approximation Monte Carlo and Wang-Landau Monte Carlo applied to a continuum polymer model , 2015, Comput. Phys. Commun..

[20]  B. Kaufman.  Crystal Statistics: II. Partition Function Evaluated by Spinor Analysis. III. Short-Range Order in a Binary Ising Lattice. , 1949 .

[21]  Wolfhard Janke,et al.  Thermodynamics and structure of macromolecules from flat-histogram Monte Carlo simulations. , 2016, Soft matter.

[22]  Beale.  Exact distribution of energies in the two-dimensional ising model. , 1996, Physical review letters.

[23]  Thomas Wüst,et al.  Optimized Wang-Landau sampling of lattice polymers: ground state search and folding thermodynamics of HP model proteins. , 2012, The Journal of chemical physics.

[24]  Q. Yan,et al.  Fast calculation of the density of states of a fluid by Monte Carlo simulations. , 2003, Physical review letters.

[25]  R. Belardinelli,et al.  Analysis of the convergence of the 1t and Wang-Landau algorithms in the calculation of multidimensional integrals. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.