A worm algorithm for the fully-packed loop model
Abstract
We present a Markov-chain Monte Carlo algorithm of worm type that correctly simulates the fully-packed loop model with n=1 on the honeycomb lattice, and we prove that it is ergodic and has uniform stationary distribution. The honeycomb-lattice fully-packed loop model with n=1 is equivalent to the zero-temperature triangular-lattice antiferromagnetic Ising model, which is fully frustrated and notoriously difficult to simulate. We test this worm algorithm numerically and estimate the dynamic exponent z=0.515(8). We also measure several static quantities of interest, including loop-length and face-size moments. It appears numerically that the face-size moments are governed by the magnetic dimension for percolation.
- Publication:
-
Nuclear Physics B
- Pub Date:
- June 2009
- DOI:
- 10.1016/j.nuclphysb.2009.01.007
- arXiv:
- arXiv:0811.2042
- Bibcode:
- 2009NuPhB.814..461Z
- Keywords:
-
- 02.70.Tt;
- 05.10.Ln;
- 64.60.De;
- 64.60.F-;
- Justifications or modifications of Monte Carlo methods;
- Monte Carlo methods;
- Statistical mechanics of model systems;
- Equilibrium properties near critical points critical exponents;
- Condensed Matter - Statistical Mechanics
- E-Print:
- 31 pages, 10 figures. Several new figures added, and some minor typos corrected. Uses the following latex packages: algorithm.sty, algorithmic.sty, elsart-num.bst, elsart1p.cls, elsart.cls