No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
The full text of this article is not currently available.
s
A new approach to Monte Carlo simulations in statistical physics: Wang-Landau sampling
Abstract
We describe a Monte Carlo algorithm for doing simulations in classical statistical physics in a different way. Instead of sampling the probability distribution at a fixed temperature, a random walk is performed in energy space to extract an estimate for the density of states. The probability can be computed at any temperature by weighting the density of states by the appropriate Boltzmann factor. Thermodynamic properties can be determined from suitable derivatives of the partition function and, unlike “standard” methods, the free energy and entropy can also be computed directly. To demonstrate the simplicity and power of the algorithm, we apply it to models exhibiting first-order or second-order phase transitions.
© 2004 American Association of Physics Teachers
Received 15 December 2003
Accepted 20 February 2004
Published online 13 September 2004
/content/aapt/journal/ajp/72/10/10.1119/1.1707017
1.
1.D. P. Landau and K. Binder, A Guide to Monte Carlo Methods in Statistical Physics (Cambridge U. P., Cambridge, 2000).
2.
2.D. P. Landau and R. Alben, “Monte Carlo calculations as an aid in teaching statistical mechanics,” Am. J. Phys. 41, 394–400 (1973).
3.
3.N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, “Equation of state calculations by fast computing machines,” J. Chem. Phys. 21, 1087–1092 (1953).
4.
4.R. H. Swendsen and J.-S. Wang, “Nonuniversal critical dynamics in Monte Carlo simulations,” Phys. Rev. Lett. 58, 86–88 (1987).
5.
5.U. Wolff, “Collective Monte Carlo updating for spin systems,” Phys. Rev. Lett. 62, 361–364 (1989).
6.
6.B. A. Berg and T. Neuhaus, “Multicanonical ensemble: A new approach to simulate first-order phase transitions,” Phys. Rev. Lett. 68, 9–12 (1992).
7.
7.W. Janke and S. Kappler, “Multibondic cluster algorithm for Monte Carlo simulations of first-order phase transitions,” Phys. Rev. Lett. 74, 212–215 (1995).
8.
8.W. Janke, “Multicanonical simulation of the two-dimensional 7-state Potts model,” Int. J. Mod. Phys. C 3, 1137–1146 (1992).
9.
9.B. A. Berg, U. Hansmann, and T. Neuhaus, “Simulation of an ensemble with varying magnetic field: A numerical determination of the order-order interface tension in the D=2 Ising model,” Phys. Rev. B 47, 497–500 (1993).
10.
10.W. Janke, “Multicanonical Monte Carlo simulations,” Physica A 254, 164–178 (1998).
11.
11.B. A. Berg and T. Celik, “New approach to spin-glass simulations,” Phys. Rev. Lett. 69, 2292–2295 (1992).
12.
12.N. A. Alves and U. H. E. Hansmann, “Partition function zeros and finite size scaling of helix-coil transitions in a polypeptide,” Phys. Rev. Lett. 84, 1836–1839 (2000).
13.
13.A. M. Ferrenberg and R. H. Swendsen, “New Monte Carlo technique for studying phase transitions,” Phys. Rev. Lett. 61, 2635–2638 (1988);
13.A. M. Ferrenberg and R. H. Swendsen, “Optimized Monte Carlo data analysis,” Phys. Rev. Lett. 63, 1195–1198 (1989).
14.
14.P. D. Beale, “Exact distribution of energies in the two-dimensional Ising Model,” Phys. Rev. Lett. 76, 78–81 (1996).
15.
15.J. Lee, “New Monte Carlo algorithm: Entropic sampling,” Phys. Rev. Lett. 71, 211–214 (1993).
16.
16.P. M. C. de Oliveira, T. J. P. Penna, and H. J. Herrmann, “Broad histogram method,” Braz. J. Phys. 26, 677–683 (1996);
16.P. M. C. de Oliveira, T. J. P. Penna, and H. J. Herrmann, “Broad histogram Monte Carlo,” Eur. Phys. J. B 1, 205–208 (1998).
17.
17.J.-S. Wang and L. W. Lee, “Monte Carlo algorithms based on the number of potential moves,” Comput. Phys. Commun. 127, 131–136 (2000).
18.
18.A. R. Lima, P. M. C. de Oliveira, and T. J. P. Penna, “A comparison between broad histogram and multicanonical methods,” J. Stat. Phys. 99, 691–705 (2000).
19.
19.F. Wang and D. P. Landau, “Efficient, multiple-range random walk algorithm to calculate the density of states,” Phys. Rev. Lett. 86, 2050–2053 (2001);
19.F. Wang and D. P. Landau, “Determining the density of states for classical statistical models: A random walk algorithm to produce a flat histogram,” Phys. Rev. E 64, 056101 (2001).
20.
20.C. Yamaguchi and Y. Okabe, “Three-dimensional antiferromagnetic q-state Potts models: application of the Wang-Landau algorithm,” J. Phys. A 34, 8781–8794 (2001).
21.
21.Y. Okabe, Y. Tomita, and C. Yamaguchi, “Application of new Monte Carlo algorithms to random spin systems,” Comput. Phys. Commun. 146, 63–68 (2002).
22.
22.M. Troyer, S. Wessel, and F. Alet, “Flat histogram methods for quantum systems: Algorithms to overcome tunneling problems and calculate the free energy,” Phys. Rev. Lett. 90, 120201 (2003).
23.
23.P. N. Vorontsov-Velyaminov and A. P. Lyubartsev, “Entropic sampling in the path integral Monte Carlo method,” J. Phys. A 36, 685–693 (2003).
24.
24.W. Koller, A. Prüll, H. G. Evertz, and W. von der Linden, “Uniform hopping approach to the ferromagnetic Kondo model at finite temperature,” Phys. Rev. B 67, 104432 (2003).
25.
25.T. S. Jain and J. J. de Pablo, “Calculation of interfacial tension from density of states,” J. Chem. Phys. 118, 4226–4229 (2003).
26.
26.Q. L. Yan, R. Faller, and J. J. de Pablo, “Density-of-states Monte Carlo method for simulation of fluids,” J. Chem. Phys. 116, 8745–8749 (2002).
27.
27.R. Faller and J. J. de Pablo, “Density of states of a binary Lennard-Jones glass,” J. Chem. Phys. 119, 4405–4408 (2003).
28.
28.E. B. Kim, R. Faller, Q. Yan, N. L. Abbott, and J. J. de Pablo, “Potential of mean force between a spherical particle suspended in a nematic liquid crystal and a substrate,” J. Chem. Phys. 117, 7781–7787 (2002).
29.
29.T. S. Jain and J. J. de Pablo, “A biased Monte Carlo technique for calculation of the density of states of polymer films,” J. Chem. Phys. 116, 7238–7243 (2002).
30.
30.N. Rathore and J. J. de Pablo, “Monte Carlo simulation of proteins through a random walk in energy space,” J. Chem. Phys. 116, 7225–7230 (2002).
31.
31.N. Rathore, T. A. Knotts, and J. J. de Pablo, “Density of states simulations of proteins,” J. Chem. Phys. 118, 4285–4290 (2003).
32.
32.T. J. H. Vlugt, “Measurement of chemical potentials of systems with strong excluded volume interactions by computing the density of states,” Mol. Phys. 100, 2763–2771 (2002).
33.
33.F. Calvo, “Sampling along reaction coordinates with the Wang-Landau method,” Mol. Phys. 100, 3421–3427 (2002).
34.
34.F. Calvo and P. Parneix, “Statistical evaporation of rotating clusters. I. Kinetic energy released,” J. Chem. Phys. 119, 256–264 (2003).
35.
35.M. A. de Menezes and A. R. Lima, “Using entropy-based methods to study general constrained parameter optimization problems,” Physica A 323, 428–434 (2003).
36.
36.V. Mustonen and R. Rajesh, “Numerical estimation of the asymptotic behaviour of solid partitions of an integer,” J. Phys. A 36, 6651–6659 (2003).
37.
37.C. Yamaguchi and N. Kawashima, “Combination of improved multibondic method and the Wang-Landau method,” Phys. Rev. E 65, 056710 (2002).
38.
38.B. J. Schulz, K. Binder, and M. Müller, “Flat histogram method of Wang-Landau and N-fold way,” Int. J. Mod. Phys. C 13, 477–494 (2002).
39.
39.A. Hüller and M. Pleimling, “Microcanonical determination of the order parameter critical exponent,” Int. J. Mod. Phys. C 13, 947–956 (2002).
40.
40.M. S. Shell, P. G. Debenedetti, and A. Z. Panagiotopoulos, “Generalization of the Wang-Landau method for off-lattice simulations,” Phys. Rev. E 66, 056703 (2002).
41.
41.B. J. Schulz, K. Binder, M. Müller, and D. P. Landau, “Avoiding boundary effects in Wang-Landau sampling,” Phys. Rev. E 67, 067102 (2003).
42.
42.See, for example, L. Onsager, “Crystal statistics. I. A two-dimensional model with an order-disorder transition,” Phys. Rev. 65, 117–149 (1944).
43.
43.F. Y. Wu, “The Potts model,” Rev. Mod. Phys. 54, 235–268 (1982).
44.
44.B. A. Berg and W. Janke, “Multioverlap simulations of the 3D Edwards-Anderson Ising spin glass,” Phys. Rev. Lett. 80, 4771–4774 (1998).
45.
45.B. A. Berg, “Algorithmic aspects of multicanonical simulations,” Nucl. Phys. B 63, 982–984 (1998).
46.
46.B. A. Berg, T. Celik, and U. Hansmann, “Multicanonical study of the 3D Ising spin-glass,” Europhys. Lett. 22, 63–68 (1993).
47.
47.N. Hatano and J. E. Gubernatis, “Bivariate multicanonical Monte Carlo of the 3D ±J spin glass,” in Computer Simulation Studies in Condensed Matter Physics XII, edited by D. P. Landau, S. P. Lewis, and H.-B. Schüttler (Springer, Berlin, 2000), pp. 149–161.
48.
48.B. A. Berg and U. H. E. Hansmann, “Configuration space for random walk dynamics,” Eur. Phys. J. B 6, 395–398 (1998).
49.
49.U. H. E. Hansmann, “Effective way for determination of multicanonical weights,” Phys. Rev. E 56, 6200–6203 (1997).
50.
50.U. H. E. Hansmann and Y. Okamoto, “Monte Carlo simulations in generalized ensemble: Multicanonical algorithm versus simulated tempering,” Phys. Rev. E 54, 5863–5865 (1996).
51.
51.W. Janke, B. A. Berg, and A. Billoire, “Multi-overlap simulations of free-energy barriers in the 3D Edwards-Anderson Ising spin glass,” Comput. Phys. Commun. 121/122, 176–179 (1999).
52.
52.J.-S. Wang, T. K. Tay, and R. H. Swendsen, “Transition matrix Monte Carlo reweighting and dynamics,” Phys. Rev. Lett. 82, 476–479 (1999).
53.
53.D. P. Landau, “Finite-size behavior of the Ising square lattice,” Phys. Rev. B 13, 2997–3011 (1976).
54.
54.A. E. Ferdinand and M. E. Fisher, “Bounded and inhomogeneous Ising models. I. Specific-heat anomaly of a finite lattice,” Phys. Rev. 185, 832–846 (1969).
55.
55.K. Binder, K. Vollmayr, H. P. Deutsch, J. D. Reger, M. Scheucher, and D. P. Landau, “Monte Carlo methods for first-order phase transitions: some recent progress,” Int. J. Mod. Phys. C 3, 1025–1058 (1992).
56.
56.M. S. S. Challa, D. P. Landau, and K. Binder, “Finite-size effects at temperature-driven first-order transitions,” Phys. Rev. B 34, 1841–1852 (1986).
57.
57.See EPAPS Document No.
E-AJPIAS-72-006406 for a sample code of the Wang-Landau algorithm for the 2D Ising model. This document may also be retrieved via the EPAPS homepage (http://www.aip.org/pubservs/epaps.html) or from ftp.aip.org in the directory /epaps. See the EPAPS homepage for more information.
[Supplementary Material]
http://aip.metastore.ingenta.com/content/aapt/journal/ajp/72/10/10.1119/1.1707017
Article metrics loading...
/content/aapt/journal/ajp/72/10/10.1119/1.1707017
2004-09-13
2016-04-07
Most read this month
Article
content/aapt/journal/ajp
Journal
5
3
true
Commenting has been disabled for this content