VLSI computing architectures for Ising model simulation

Abstract Efficient VLSI architectures for parallel computations involved in the Monte Carlo simulation of the Ising model are presented. The architecture provides an extremely area-efficient, spatially-distributed set of pseudorandom numbers, which are required in the local nondeterministic decisions at the various sites in the lattice, using pseudorandom number generators based upon cellular automata (CA). This provides a highly parallel fine-grained VLSI architecture specifically aimed at Ising model computations. It is shown that the time-intensive task of sampling the Ising configurations is optimised by this approach. Therefore, these architectures can be considered to form the basis for Ising computers which could act as coprocessors for computational statistical mechanics in much the same manner as a floating-point accelerator is used to speed arithmetic computation. These hardware experts can be used to report pertinent information such as the magnetisation to a host computer. It is demonstrated that this architecture can provide a speedup of several orders of magnitude over conventional Monte Carlo simulation and are also an improvement over previous parallel computations. Measurements from a prototype constructed using a custom VLSI chip implementation indicate that a 1000 × 1000 Ising lattice can be completely updated and lattice energy and magnetisation calculated in less than 1 μsec. The validity of this approach is verified by computer simulation of the new architecture which yielded the correct Ising model critical exponents and lattice energy and magnetisation for the two-dimensional square lattice. This provides a stringent test against long term correlation effects in the CA-based distributed random numbers over and above our earlier verification using standard empirical random number tests.

[1]  A. Compagner,et al.  A special-purpose processor for the Monte Carlo simulation of ising spin systems , 1983 .

[2]  PETER D. HORTENSIUS,et al.  Importance Sampling for Ising Computers Using One-Dimensional Cellular Automata , 1989, IEEE Trans. Computers.

[3]  Peter Dirk Hortensius Parallel computation of non-deterministic algorithms in vlsi , 1987 .

[4]  G. Vichniac Simulating physics with cellular automata , 1984 .

[5]  B. McCoy,et al.  The Two-Dimensional Ising Model , 1973 .

[6]  A. Verhagen,et al.  An exactly soluble case of the triangular ising model in a magnetic field , 1976 .

[7]  K. Binder Monte Carlo methods in statistical physics , 1979 .

[8]  L. Onsager Crystal statistics. I. A two-dimensional model with an order-disorder transition , 1944 .

[9]  John L. Richardson,et al.  A Fast Processor for Monte-Carlo Simulation , 1983 .

[10]  Stephen Wolfram,et al.  Universality and complexity in cellular automata , 1983 .

[11]  Ian G. Enting,et al.  Crystal growth models and Ising models: Disorder points , 1977 .

[12]  Kurt Binder,et al.  Monte Carlo Investigation of Dynamic Critical Phenomena in the Two-Dimensional Kinetic Ising Model , 1973 .

[13]  D. Landau Phase transitions in the Ising square lattice with next-nearest-neighbor interactions , 1980 .

[14]  D. Landau,et al.  Finite-size behavior of the simple-cubic Ising lattice , 1976 .

[15]  S. F. Reddaway,et al.  A very high speed Monte Carlo simulation on DAP , 1985 .

[16]  Robert H. Swendsen,et al.  Statistical mechanics and disordered systems , 1985, CACM.

[17]  Kurt Binder,et al.  Phase diagrams and critical behavior in Ising square lattices with nearest- and next-nearest-neighbor interactions , 1980 .

[18]  Eytan Domany,et al.  Equivalence of Cellular Automata to Ising Models and Directed Percolation , 1984 .

[19]  S. Wolfram Statistical mechanics of cellular automata , 1983 .

[20]  S. Wolfram Random sequence generation by cellular automata , 1986 .

[21]  Pearson,et al.  Finite-size scaling in the three-dimensional Ising model. , 1985, Physical review. B, Condensed matter.

[22]  Howard C. Card,et al.  Parallel Random Number Generation for VLSI Systems Using Cellular Automata , 1989, IEEE Trans. Computers.

[23]  David P. Landau,et al.  Finite-size behavior of the Ising square lattice , 1976 .

[24]  K. Wilson Problems in Physics with many Scales of Length , 1979 .

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

[26]  J. Hammersley,et al.  Monte Carlo Methods , 1965 .