Mean-value identities as an opportunity for Monte Carlo error reduction

Physical Review E - Statistical, Nonlinear, and Soft Matter Physics

Volumn 79, Issue 5, 2009, Pages

  Fernandez, L A ^a,b   Martin Mayor, V ^a,b  

^a UNIVERSIDAD COMPLUTENSE DE MADRID   (Spain)

^b UNIVERSITY OF ZARAGOZA   (Spain)

Author keywords

[No Author keywords available]

Indexed keywords

CONTROL VARIATES; CPU TIME; DYSON EQUATIONS; ERROR REDUCTION; FIELD THEORY; GAIN FACTORS; MEAN VALUES; MONTE CARLO; MONTE CARLO SIMULATION; NUMBER GENERATOR; PSEUDO RANDOM; STATISTICAL ERRORS; THERMALIZATION; TWO-DIMENSIONAL ISING MODEL;

ISING MODEL; MONTE CARLO METHODS; NUMBER THEORY; STATISTICAL MECHANICS;

STATISTICAL TESTS;

EID: 66749174537     PISSN: 15393755     EISSN: 15502376     Source Type: Journal    
DOI: 10.1103/PhysRevE.79.051109     Document Type: Article

References (32)

1. D. P. Landau and K. Binder, A Guide to Monte Carlo Simulations in Statistical Physics (Cambridge University Press, Cambridge, 2000).
2. A. D. Sokal, in Functional Integration: Basics and Applications (1996 Cargèse School), edited by, C. DeWitt-Morette, P. Cartier, and, A. Folacci, (Plenum, New York, 1997).
3. G. Parisi, R. Petronzio, and F. Rapuano, Phys. Lett. B 128, 418 (1983). 10.1016/0370-2693(83)90930-9
4. M. Lüscher and P. Weisz, J. High Energy Phys. 09 (2001) 010. 10.1088/1126-6708/2001/09/010
5. See, e.g., F. Belletti (the Janus Collaboration), Phys. Rev. Lett. 101, 157201 (2008). 10.1103/PhysRevLett.101.157201
6. J. M. Hammersley and D. C. Handscomb, Monte Carlo Methods (Fletcher and Son, Methuen, London, 1964).
7. R. Y. Rubinstein, Simulation and the Monte Carlo Method (Wiley and Sons, New York, 1981).
8. J. S. Schwinger, Proc. Natl. Acad. Sci. U.S.A. 37, 452 (1951) 10.1073/pnas.37.7.452
9. F. J. Dyson, Phys. Rev. 75, 1736 (1949). 10.1103/PhysRev.75.1736
10. H. B. Callen, Phys. Lett. 4, 161 (1963). 10.1016/0031-9163(63)90344-5
11. F. Guerra, Int. J. Mod. Phys. B 10, 1675 (1996). 10.1142/ S0217979296000751
12. R. H. Swendsen and J. S. Wang, Phys. Rev. Lett. 58, 86 (1987). 10.1103/PhysRevLett.58.86
13. J. Salas and A. D. Sokal, J. Stat. Phys. 98, 551 (2000). 10.1023/A:1018611122166
14. E. Marinari, G. Parisi, and J. J. Ruiz-Lorenzo, Phys. Rev. B 58, 14852 (1998). 10.1103/PhysRevB.58.14852
15. H. G. Katzgraber, M. Palassini, and A. P. Young, Phys. Rev. B 63, 184422 (2001). 10.1103/PhysRevB.63.184422
16. See, e.g., H. G. Ballesteros and V. Martin-Mayor, Phys. Rev. E 58, 6787 (1998), and references therein. 10.1103/PhysRevE.58.6787
17. H. G. Ballesteros, L. A. Fernandez, V. Martin-Mayor, and A. Munoz Sudupe, Phys. Lett. B 387, 125 (1996) 10.1016/0370-2693(96)00984-7
18. H. G. Ballesteros, L. A. Fernandez, V. Martin-Mayor, A. Munoz Sudupe, G. Parisi, and J. J. Ruiz-Lorenzo, Phys. Rev. B 58, 2740 (1998). 10.1103/PhysRevB.58.2740
19. M. Weigel and W. Janke, Phys. Rev. Lett. 102, 100601 (2009). 10.1103/PhysRevLett.102.100601
20. For disordered systems, the time average in Eq. 9 is followed by a disorder average, which strongly diminishes the influence of the particular Monte Carlo dynamics.
21. D. Amit and V. Martin-Mayor, Field Theory: The Renormalization Group and Critical Phenomena, 3rd ed. (World-Scientific, Singapore, 2005).
22. S. S. Lavenberg and P. D. Welch, IBM Corporation, Yorktown Heights, New York Research Report No. RC8161, 1980 (unpublished).
23. A. E. Ferdinand and M. E. Fischer, Phys. Rev. 185, 832 (1969). 10.1103/PhysRev.185.832
24. F. Cooper, B. Freedman, and D. Preston, Nucl. Phys. B 210, 210 (1982). 10.1016/0550-3213(82)90240-1
25. L. A. Fernandez, V. Martin-Mayor, and D. Yllanes, Nucl. Phys. B 807, 424 (2009). 10.1016/j.nuclphysb.2008.07.009
26. L. A. Fernandez, V. Martin-Mayor, S. Perez-Gaviro, A. Tarancon, and A. P. Young e-print arXiv:0905.0322.
27. We shall need to recall the main features of the Swendsen-Wang dynamics, as formulated in Ref. The configuration space is expanded; to the original spins { S x } in the lattice sites, we add bond variables b x y =0,1 in the lattice bonds joining nearest-neighbor sites x and y (the marginal probability for the spins is the Boltzmann weight of the original Ising model). The b=1 bonds are said to be occupied. Two lattice sites connected by a chain of occupied bonds are said to belong to the same cluster. The clusters partition the lattice into connected components. A crucial role is played by the conditional probabilities. That of the spins, given the bonds is: all the spins in a cluster take the same value, the spin values in different clusters are statistically independent and equal ±1 with 50% probability. For that of the bonds, given the spins, different bonds are statistically independent, and a given bond can be one only if it connects spins of equal sign and, in that case, with probability 1- e-2κ. The two Schrödinger-type operators, Pspin and Pbond belong to the heat-bath category; Pspin (Pbond) leave the bonds (spins) unchanged and choose the new spin (bond) configuration according to the conditional probability of the spins (bonds), given the bonds (spins). As all heat-bath dynamics, both Pbond and Pspin verify detailed balance. Their product PSW = Pspin Pbond verifies only the softer balance condition, Eq. 10.
28. X. J. Li and A. D. Sokal, Phys. Rev. Lett. 63, 827 (1989). 10.1103/PhysRevLett.63.827
29. Combining ηCt ηSt □ =0 with Eqs. 3 8 9 10, we recover the well-known cluster estimator for U4: U4 = 3 (□c nc2) 2 -2 □c nc4 □ □c nc2 □ 2.
30. For finite L, gC >0 while the large- L limit for U4 is 3.
31. For instance, M2 □ = Pspin M2 □ = C□. A telegraphic proof of CC BSW (0) =0 is also obtained by combining C M2 □ = Pspin C M2 □ with Pspin C M2 =C Pspin M2 = C2 and with Eq. 18.
32. R. G. Edwards and A. D. Sokal, Phys. Rev. D 38, 2009 (1988). 10.1103/PhysRevD.38.2009