Comparing Monte Carlo methods for finding ground states of Ising spin glasses: Population annealing, simulated annealing, and parallel tempering.

Population annealing is a Monte Carlo algorithm that marries features from simulated-annealing and parallel-tempering Monte Carlo. As such, it is ideal to overcome large energy barriers in the free-energy landscape while minimizing a Hamiltonian. Thus, population-annealing Monte Carlo can be used as a heuristic to solve combinatorial optimization problems. We illustrate the capabilities of population-annealing Monte Carlo by computing ground states of the three-dimensional Ising spin glass with Gaussian disorder, while comparing to simulated-annealing and parallel-tempering Monte Carlo. Our results suggest that population annealing Monte Carlo is significantly more efficient than simulated annealing but comparable to parallel-tempering Monte Carlo for finding spin-glass ground states.

[1]  Alexander K. Hartmann,et al.  SCALING OF STIFFNESS ENERGY FOR THREE-DIMENSIONAL J ISING SPIN GLASSES , 1999 .

[2]  F. Barahona.  On the computational complexity of Ising spin glass models , 1982 .

[3]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[4]  Cristopher Moore,et al.  The Nature of Computation , 2011 .

[5]  A. Young,et al.  Triviality of the Ground State Structure in Ising Spin Glasses , 1999, cond-mat/9906323.

[6]  A. Middleton,et al.  Improved extremal optimization for the Ising spin glass. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[7]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[8]  James Edward Gubernatis,et al.  The Monte Carlo Method in the Physical Sciences: Celebrating the 50th Anniversary of the Metropolis Algorithm , 2003 .

[9]  J. Machta,et al.  Population annealing with weighted averages: a Monte Carlo method for rough free-energy landscapes. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  Heiko Rieger,et al.  Optimization algorithms in physics , 2001 .

[11]  Helmut G. Katzgraber,et al.  Evidence against a mean-field description of short-range spin glasses revealed through thermal boundary conditions , 2014, 1408.0438.

[12]  H. Rieger,et al.  New Optimization Algorithms in Physics , 2004 .

[13]  K Hukushima.  Domain-wall free energy of spin-glass models: numerical method and boundary conditions. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[14]  Elena Bautu,et al.  Searching Ground States of Ising Spin Glasses with Genetic Algorithms and Binary Particle Swarm Optimization , 2007, NICSO.

[15]  M. Troyer,et al.  Optimized parallel tempering simulations of proteins. , 2006, The Journal of chemical physics.

[16]  F. Nieto,et al.  The ground state energy of the Edwards–Anderson spin glass model with a parallel tempering Monte Carlo algorithm , 2008, 0806.1054.

[17]  G. Rinaldi,et al.  Exact ground states of Ising spin glasses: New experimental results with a branch-and-cut algorithm , 1995 .

[18]  R. Ellis,et al.  Monte Carlo Methods for Rough Free Energy Landscapes: Population Annealing and Parallel Tempering , 2011, 1104.1138.

[19]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[20]  Proceedings of the 2010 Winter Simulation Conference, WSC 2010, Baltimore, Maryland, USA, 5-8 December 2010 , 2010, Winter Simulation Conference.

[21]  Finding Low-Temperature States With Parallel Tempering, Simulated Annealing And Simple Monte Carlo , 2002, cond-mat/0209248.

[22]  Timothy J. Robinson,et al.  Sequential Monte Carlo Methods in Practice , 2003 .