Cluster heat bath method on a quasi-one-dimensional ising model

We have proposed a cluster heat bath (CHB) method of a quasi-one-dimensionaI Ising model and demonstrated that it reproduces the exact results of the two-dimensional Ising model with axially anisotropic exchange interactions. The key point of this method is to select one of the equilibrium spin configurations of each of the chains which are subjected to effective fields of surrounding spins. Thus, the chains are always in their equilibrium state during the simulation, and the d-dimensional model effectively becomes a (d - 1)-dimensional model of fictitious spins with their freedom of 2(N), where N is the number of the spins on the chain, and the simulation time is drastically reduced.