Microcanonical simulated annealing: Massively parallel Monte Carlo simulations with sporadic random-number generation.
Saved in:
| Title: | Microcanonical simulated annealing: Massively parallel Monte Carlo simulations with sporadic random-number generation. |
|---|---|
| Authors: | Bernaschi, M.1 (AUTHOR), Chilin, C.2,3 (AUTHOR), Fernandez, L.A.2 (AUTHOR), González-Adalid Pemartín, I.1 (AUTHOR), Marinari, E.3,4 (AUTHOR), Martin-Mayor, V.2 (AUTHOR), Parisi, G.3,4,5 (AUTHOR), Ricci-Tersenghi, F.3,4 (AUTHOR), Ruiz-Lorenzo, J.J.6,7 (AUTHOR), Yllanes, D.1,8,9,10 (AUTHOR) david.yllanes@unizar.es |
| Source: | Computer Physics Communications. Aug2026, Vol. 325, pN.PAG-N.PAG. 1p. |
| Subjects: | Monte Carlo method, Random number generators, Simulated annealing, Parallel programming, Combinatorial optimization, Graphics processing units, Ising model |
| Abstract: | Numerical simulations of models and theories that describe complex systems such as spin glasses are becoming increasingly important. Beyond fundamental research, these computational methods also find practical applications in fields like combinatorial optimization. However, Monte Carlo simulations, an important subcategory of these methods, are plagued by a major drawback: they are extremely greedy for (pseudo) random numbers. The total fraction of computer time dedicated to random-number generation increases as the hardware grows more sophisticated, and can get prohibitive for special-purpose computing platforms. We propose here a general-purpose microcanonical simulated annealing (MicSA) formalism that dramatically reduces such a burden. The algorithm is fully adapted to a massively parallel computation, as we show in the particularly demanding benchmark of the three-dimensional Ising spin glass. We carry out very stringent numerical tests of the new algorithm by comparing our results, obtained on GPUs, with high-precision standard (i.e. , random-number-greedy) simulations performed on the Janus II custom-built supercomputer. In those cases where thermal equilibrium is reachable (i.e. , in the paramagnetic phase), both simulations reach compatible values. More significantly, barring short-time corrections, a simple time rescaling suffices to map the MicSA off-equilibrium dynamics onto the results obtained with standard simulations. [ABSTRACT FROM AUTHOR] |
| Copyright of Computer Physics Communications is the property of Elsevier B.V. and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.) | |
| Database: | Engineering Source |
Be the first to leave a comment!
