Monte Carlo and mean field slow cooling simulations for spin glasses: relation to NP-completeness

We discuss the dependence of the ground state energy as the Monte Carlo (MC) cooling rate in a variety of spin glass and other optimization problems. Our work is motivated by recent interest in the concept of simulated annealing and the need to find efficient numerical optimization schemes. We find that the ground state energy for six spin-glass models and the optimal tour length for the N-city traveling salesman problem are a sensitive function of the cooling rate. For the one dimensional nearest neighbor Gaussian and two two dimensional spin-glass models, the nearest neighbor ± J and Gaussian models, we find E(r) = E0 + C1 rx, while for the three dimensional t J model, a two layer t J model, the infinite range model and the N-city traveling salesman problem, the dependence on r is slower, E(r) = Eo - C2 (Inr)-1.Here C1 and C2 are constants.We believe that this difference is related to the fact that finding the ground state energy for the latter four problems is an intractable problem, i.e., it is Np-complete. We assert that a logarithmic cooling rate dependence for finding the ground state is a necessary consequence of Np-completeness. Non Np-complete problems may depend on the cooling rate either as a power-law or logarithmically. However, we have found no example for the latter.

[1]  Fisher,et al.  Residual energies after slow cooling of disordered systems. , 1986, Physical review letters.

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

[3]  Scott Kirkpatrick,et al.  Optimization by simulated annealing: Quantitative studies , 1984 .

[4]  Paolo Carnevali,et al.  Image Processing by Simulated Annealing , 1985, IBM J. Res. Dev..

[5]  Young,et al.  Search for a transition in the three-dimensional , 1985, Physical review letters.

[6]  A. P. Young,et al.  Lack of Ergodicity in the Infinite-Range Ising Spin-Glass , 1982 .

[7]  Levin,et al.  Cooling-rate dependence for the spin-glass ground-state energy: Implications for optimization by simulated annealing. , 1986, Physical review letters.

[8]  Daniel H. Rothman,et al.  Automatic estimation of large residual statics corrections , 1986 .

[9]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[10]  D. M. Nicholson,et al.  Monte Carlo optimization of pair distribution functions: Application to the electronic structure of disordered metals , 1984 .

[11]  N-city traveling salesman problem: Optimization by simulated annealings , 1986 .

[12]  David S. Johnson,et al.  Computers and Intractability: A Guide to the Theory of NP-Completeness , 1978 .

[13]  A. B. Bortz,et al.  A new algorithm for Monte Carlo simulation of Ising spin systems , 1975 .

[14]  Giorgio Parisi,et al.  Infinite Number of Order Parameters for Spin-Glasses , 1979 .

[15]  E. Bonomi,et al.  The N-City Travelling Salesman Problem: Statistical Mechanics and the Metropolis Algorithm , 1984 .

[16]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[17]  S. Kirkpatrick,et al.  Infinite-ranged models of spin-glasses , 1978 .

[18]  Morgenstern,et al.  Critical behavior of three-dimensional Ising spin-glass model. , 1985, Physical review letters.

[19]  I. Morgenstern,et al.  Evidence Against Spin-Glass Order in the Two-Dimensional Random-Bond Ising Model , 1979 .

[20]  Gary S. Grest,et al.  Irreversibility and metastability in spin-glasses. I. Ising model , 1983 .

[21]  Scott Kirkpatrick,et al.  Global Wiring by Simulated Annealing , 1983, IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems.

[22]  S. Kirkpatrick,et al.  Solvable Model of a Spin-Glass , 1975 .

[23]  Rammile Ettelaie,et al.  Residual entropy and simulated annealing , 1985 .

[24]  P. Carnevali,et al.  Image processing by stimulated annealing , 1985 .

[25]  H H Barrett,et al.  Reconstruction of objects from coded images by simulated annealing. , 1983, Optics letters.

[26]  Steven G. Louie,et al.  A Monte carlo simulated annealing approach to optimization over continuous variables , 1984 .