Abstract.The Binder cumulant at the phase transition of Ising models on square lattices with various ferromagnetic nearest and next-nearest neighbour couplings is determined using mainly Monte Carlo techniques. We discuss the possibility to relate the value of the critical cumulant in the isotropic, nearest neighbour and in the anisotropic cases to each other by means of a scale transformation in rectangular geometry, to pinpoint universal and nonuniversal features.
We present a progress report on the Cluster Processor, a special-purpose computer system for the Wolff simulation of the three-dimensional Ising model, including an analysis of simulation results obtained thus far. These results allow, within narrow error margins, a determination of the parameters describing the phase transition of the simple-cubic Ising model and its universality class. For an improved determination of the correctionto-scaling exponent, we include Monte Carlo data for systems with nearest-neighbor and third-neighbor interactions in the analysis.
Population annealing is a promising recent approach for Monte Carlo simulations in statistical physics, in particular for the simulation of systems with complex free-energy landscapes. It is a hybrid method, combining importance sampling through Markov chains with elements of sequential Monte Carlo in the form of population control. While it appears to provide algorithmic capabilities for the simulation of such systems that are roughly comparable to those of more established approaches such as parallel tempering, it is intrinsically much more suitable for massively parallel computing. Here, we tap into this structural advantage and present a highly optimized implementation of the population annealing algorithm on GPUs that promises speed-ups of several orders of magnitude as compared to a serial implementation on CPUs. While the sample code is for simulations of the 2D ferromagnetic Ising model, it should be easily adapted for simulations of other spin models, including disordered systems. Our code includes implementations of some advanced algorithmic features that have only recently been suggested, namely the automatic adaptation of temperature steps and a multi-histogram analysis of the data at different temperatures. The program calculates the internal energy, specific heat, several magnetization moments, entropy and free energy of the 2D Ising model on square lattices of edge length L with periodic boundary conditions as a function of inverse temperature β. Solution method: The code uses population annealing, a hybrid method combining Markov chain updates with population control. The code is implemented for NVIDIA GPUs using the CUDA language and employs advanced techniques such as multi-spin coding, adaptive temperature steps and multi-histogram reweighting. Restrictions: The system size and size of the population of replicas are limited depending on the memory of the GPU device used. Unusual features: Additional comments: Code repository at https://github.com/LevBarash/PAising. Running time: For the default parameter values used in the sample programs, L = 64, θ = 100, β0 = 0, β f = 1, ∆β = 0.005, R = 20 000, a typical run time on an NVIDIA Tesla K80 GPU is 151 seconds for the single spin coded (SSC) and 17 seconds for the multi-spin coded (MSC) program (see Sec. 2 for a description of these parameters).
PROGRAM SUMMARY
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.