Monte Carlo simulation of ice models

Physical Review E - Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics

Volumn 57, Issue 1, 1998, Pages 1155-1166

  Barkema, G T ^a   Newman, M E J ^b  

^a FORSCHUNGSZENTRUM JÜLICH   (Germany)

^b SANTA FE INSTITUTE   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

EID: 0000227735     PISSN: 1063651X     EISSN: None     Source Type: Journal    
DOI: 10.1103/PhysRevE.57.1155     Document Type: Article

References (27)

1. J. D. Bernal and R. H. Fowler, J. Chem. Phys. 1, 515 (1933).
2. L. Pauling, J. Am. Chem. Soc. 57, 2680 (1935)
3. E. H. Lieb, Phys. Rev. 162, 162 (1967)
4. E. H. Lieb, Phys. Rev. Lett. 18, 1046 (1967); PRLTAO
5. E. H. Lieb Phys. Rev. Lett. 19, 108 (1967).
6. R. J. Baxter, Exactly Solved Models in Statistical Mechanics (Academic, London, 1982).
7. A. Lenard (unpublished).
8. R. Raghavan, C. L. Henley, and S. L. Arouh, J. Stat. Phys. 86, 517 (1997)
9. B. Nienhuis, in Phase Transitions and Critical Phenomena, edited by C. Domb and J. L. Lebowitz (Academic, London, 1987), Vol. 11.
10. J. Kondev, Phys. Rev. Lett. 78, 4320 (1997)
11. A. Rahman and F. H. Stillinger, J. Chem. Phys. 57, 4009 (1972).
12. A. Yanagawa and J. F. Nagle, Chem. Phys. 43, 329 (1979).
13. N. Bjerrum, K. Dan. Vidensk. Selsk. Mat. Fys. Skr. 27, 1 (1951); KVMFAS
14. N. BjerrumScience 115, 385 (1952).
15. P. V. Hobbs, Ice Physics (Clarendon, Oxford, 1974).
16. It is not too hard to show that the loops which wrap around the periodic boundary conditions change the polarization, Eq. (9), of the system, whereas the ones which do not conserve polarization. Thus if we do not allow the loops to wrap around in this way, the polarization would never change.
17. U. Wolff, Phys. Rev. Lett. 62, 361 (1989)
18. R. H. Swendsen and J.-S. Wang, Phys. Rev. Lett. 58, 86 (1987)
19. J.-S. Wang, R. H. Swendsen, and R. Kotecky, Phys. Rev. B 42, 2465 (1990)
20. In fact, we can show that the loop size scales exactly as L5/3. This result follows from calculations using the Coulomb gas representation of the model. See H. Saleur, Nucl. Phys. B 360, 219 (1991)
21. M. P. Nightingale and H. W. J. Blöte, Phys. Rev. Lett. 76, 4548 (1996)
22. J. C. Slater, J. Chem. Phys. 9, 16 (1941).
23. F. Rys, Helv. Phys. Acta 36, 537 (1963)
24. H. J. F. Knops, Phys. Rev. B 20, 4670 (1979)
25. E. Carlon, G. Mazzeo, and H. van Beijeren, Phys. Rev. B 55, 757 (1997)
26. Another widely studied loop algorithm, the algorithm of Evertz, Lana, and Marco, possesses a much higher dynamic exponent than this, around z=0.7 [H. G. Evertz, G. Lana, and M. Marcu, Phys. Rev. Lett. 70, 875 (1993)]. PRLTAO
27. Since this algorithm is also quite complicated to implement, we recommend the short loop algorithm for simulations of the F model, or better still, the full-lattice three-color algorithm. However, the algorithm of Evertz, Lana, and Marcu has proved useful for the simulation of certain one-dimensional quantum systems [H. G. Evertz, in Numerical Methods for Lattice Quantum Many-Body Problems, edited by D. J. Scalapino (Addison-Wesley, Reading, MA, in press)].