Fast vectorized algorithm for the Monte Carlo Simulation of the Random Field Ising Model

The simulation of the two-dimensional Ising model is used as a benchmark to show the computational capabilities of Graphic Processing Units (GPUs). The rich programming environment now available on GPUs and flexible hardware capabilities allowed us to quickly experiment with several implementation ideas: a simple stencil-based algorithm, recasting the stencil operations into matrix multiplies to take advantage of Tensor Cores available on NVIDIA GPUs, and a highly optimized m...

Phase transitions in the three-dimensional diluted Ising antiferromagnet in an applied magnetic field are analyzed numerically. It is found that random magnetic field in a system with spin concentration below a certain threshold induces a crossover from second-order phase transition to first-order transition to a new phase characterized by a spin-glass ground state and metastable energy states at finite temperatures.

We have developed an efficient Monte Carlo algorithm, which accelerates slow Monte Carlo dynamics in quasi-one-dimensional Ising spin systems. The loop algorithm of the quantum Monte Carlo method is applied to the classical spin models with highly anisotropic exchange interactions. Both correlation time and real CPU time are reduced drastically. The algorithm is demonstrated in the layered triangular-lattice antiferromagnetic Ising model. We have obtained the relation between...

Graphics processing units (GPUs) are recently being used to an increasing degree for general computational purposes. This development is motivated by their theoretical peak performance, which significantly exceeds that of broadly available CPUs. For practical purposes, however, it is far from clear how much of this theoretical performance can be realized in actual scientific applications. As is discussed here for the case of studying classical spin models of statistical mecha...

This discussion serves as an introduction to the use of Monte Carlo simulations as a useful way to evaluate the observables of a ferromagnet. Key background is given about the relevance and effectiveness of this stochastic approach and in particular the applicability of the Metropolis-Hastings algorithm. Importantly the potentially devastating effects of spontaneous magnetization are highlighted and a means to avert this is examined. An Ising model is introduced and used to...

A definition for a class of asynchronous cellular arrays is proposed. An example of such asynchrony would be independent Poisson arrivals of cell iterations. The Ising model in the continuous time formulation of Glauber falls into this class. Also proposed are efficient parallel algorithms for simulating these asynchronous cellular arrays. In the algorithms, one or several cells are assigned to a processing element (PE), local times for different PEs can be different. Althoug...

We propose a kinetic Ising model to study phase separation driven by surface diffusion. This model is referred to as "Model S", and consists of the usual Kawasaki spin-exchange kinetics ("Model B") in conjunction with a kinetic constraint. We use novel multi-spin coding techniques to develop fast algorithms for Monte Carlo simulations of Models B and S. We use these algorithms to study the late stages of pattern dynamics in these systems.

Lattice spin models are useful for studying critical phenomena and allow the extraction of equilibrium and dynamical properties. Simulations of such systems are usually based on Monte Carlo (MC) techniques, and the main difficulty is often the large computational effort needed when approaching critical points. In this work, it is shown how such simulations can be accelerated with the use of NVIDIA graphics processing units (GPUs) using the CUDA programming architecture. We ha...

The multi-spin coding of the Monte Carlo simulation of the three-state Potts model on the simple cubic lattice is presented. The ferromagnetic (F) model, the antiferromagnetic (AF) model, and the random mixture of the F and AF couplings are treated. The multi-spin coding technique is also applied to the block-spin transformation. The block-spin transformation of the F Potts model is simply realized by the majority rule, whereas the AF three-state Potts model is transformed to...

We present a surprisingly simple approach to high-accuracy calculations of critical properties of the three-dimensional Ising model. The method uses a modified block-spin transformation with a tunable parameter to improve convergence in Monte Carlo renormalization group. The block-spin parameter must be tuned differently for different exponents to produce optimal convergence.