Far Eastern Mathematical Journal

To content of the issue


An enumeration decomposition algorithm for solving the Ising mode


Trukhin V.O., Lobanova E.A., Anisich A.I., Makarova K.V., Makarov A.G., Nefedev K.V.

2025, issue 1, P. 102-112
DOI: https://doi.org/10.47910/FEMJ202509


Abstract
The article presents a novel algorithm for the complete enumeration of spin configurations in the Ising model on a square lattice. Particular attention is given to the parallel algorithmization of computations on central processing units (CPUs) using OpenMP and on graphics processing units (GPUs) using CUDA. The structure of the algorithm is described, including its main implementation steps, as well as its application to solving problems in statistical thermodynamics, specifically the calculation of the density of states. A performance comparison is conducted between the proposed algorithm and sequential enumeration algorithms implemented in Python and C. The results demonstrate that the proposed approach significantly accelerates computations and enables efficient analysis of square-lattice spin systems in the Ising model with sizes up to $10 \times 10$ nodes (100 spins) and with an arbitrary distributions of exchange constants.

Keywords:
decomposition algorithm, Ising model, parallel computing.

Download the article (PDF-file)

References

[1] Bekster R., Tochno reshaemye modeli v statisticheskoi mekhanike, Mir, M., 1985.
[2] Makarova K., Makarov A., Strongin V., Titovets Iu., “Canonical Monte Carlo multispin cluster method”, Journal of Computational and Applied Mathematics, 427, (2023), 115153.
[3] Ja lowiecki K., Rams M., Gardas B., “Brute-forcing spin-glass problems with CUDA”, Computer Physics Communications, 260, (2021), 107728.
[4] Shevchenko Yu., Makarov A., Andriushchenko P., Nefedev K., “Multicanonical sampling of the space of states of H(2,n)-vector models”, Journal of Experimental and Theoretical Physics, 124, (2017), 982–993.
[5] Trukhin T.V., Strongin V.S., Chesnokov M.A., Makarov A.G., “Thermodynamic equilibrium of ±J Ising model on square lattice”, Physica A: Statistical Mechanics and its Applications, 655, (2024), 130172.
[6] Onsager L., “Crystal statistics. I. A two-dimensional model with an order-disorder transition”, Physical Review, 65, (1944), 117.
[7] Romero J., Bisson M., Fatica M., Bernaschi M., “High performance implementations of the 2D Ising model on GPUs”, Computer Physics Communications, 256, (2020).
[8] Maren A.J., “A logical topology of neural networks”, Proceedings of the Second Workshop on Neural Networks, 1991.
[9] Grant E.K., Humble T.S., “Adiabatic quantum computing and quantum annealing”, Oxford Research Encyclopedia of Physics, 2020.
[10] Korol' A.O., Kapitan V.YU., “Neironnaya set' dlya opredeleniya temperatury Kyuri dvumernoi modeli Izinga”, Dal'nevost. matem. zhurn., 21:1, (2021), 51–60.
[11] Janke W., “Monte Carlo methods in classical statistical physics”, Computational manyparticle physics, 2008, 79–140.
[12] Markovich L.A., “Parallel algorithm based on the Ising model for solving combinatorial optimization problems”, Information Technology and Systems, 2019, 350–358.
[13] Papadimitriou C. H., “The Euclidean travelling salesman problem is NP-complete”, Theoretical computer science, 4:3, (1977), 237–244.
[14] Karp R.M., “Reducibility among combinatorial problems”, Springer, 2010.
[15] Roma F., Risau-Gusman S., Ramirez-Pastor A. J., Nieto F., “Ground-state topology of the Edwards-Anderson ±J spin glass model”, Physical Review B — Condensed Matter and Materials Physics, 82:21, (2010), 214–401.
[16] Katzgraber H.G., Lee L.W., “Correlation length of the two-dimensional Ising spin glass with bimodal interactions”, Physical Review B—Condensed Matter and Materials Physics, 71:13, (2005), 134–404.
[17] Katz V.J., “The Mathematics of Egypt, Mesopotamia, China, India, and Islam: A Sourcebook”, Physics Letters A, 2007.
[18] Panchenko T.V., Tarasevich Y.Y., “Comparative analysis of the efficiency of application of genetic algorithms and Metropolis algorithm in problems of solid state physics”, Computational methods and programming, 8, (2007), 77–87.

To content of the issue