Monte Carlo study of the Yukawa coupled two-spin Ising model

We consider a particular four state spin system composed of two Ising spins (s_x, σ_x) with independent hopping parameters κ₁ κ₂, coupled by a bilinear Yukawa term, ys_xσ_x. The Yukawa term is solely responsible for breaking the global Z₂ × Z₂ symmetry down to Z2. This model is intended as an illustration of general coupled Higgs system where scalars can arise both as composite and elementary excitations. For the Ising example in 2d, we give convincing numerical evidence of the universality of the two-spin system with the one-spin Ising model, by Monte Carlo simulations and finite size scaling analysis. We also show that as we approach the phase transition, universality arises by having a single correlation length that diverges.     

Monte Carlo study of surface roughening in the three-dimensional Ising model

Above a roughening temperature of about 0.56T _c the thickness of the two-dimensional interface separating domain on a simple cubic lattice is found to increase roughly logarithmically with system size. Also, the interface tension is determined for temperatures belowT _c.     

Single-cluster Monte Carlo dynamics for the Ising model

We present an extensive study of a new Monte Carlo acceleration algorithm introduced by Wolff for the Ising model. It differs from the Swendsen-Wang algorithm by growing and flipping single clusters at a random seed. In general, it is more efficient than Swendsen-Wang dynamics ford>2, giving zero critical slowing down in the upper critical dimension. Monte Carlo simulations give dynamical critical exponentsz _w=0.33±0.05 and 0.44+0.10 ind=2 and 3, respectively, and numbers consistent withz _w=0 ind=4 and mean-field theory. We present scaling arguments which indicate that the Wolff mechanism for decorrelation differs substantially from Swendsen-Wang despite the apparent similarities of the two methods.     

Critical dynamics of the three-dimensional Ising model: A Monte Carlo study

Extensive Monte Carlo simulations have been performed to analyze the dynamical behavior of the three-dimensional Ising model with local dynamics. We have studied the equilibrium correlation functions and the power spectral densities of odd and even observables. The exponential relaxation times have been calculated in the asymptotic one-exponential time region. We find that the critical exponentz=2.09 ±0.02 characterizes the algebraic divergence with lattice size for all observables. The influence of scaling corrections has been analyzed. We have determined integrated relaxation times as well. Their dynamical exponentz _int agrees withz for correlations of the magnetization and its absolute value, but it is different for energy correlations. We have applied a scaling method to analyze the behavior of the correlation functions. This method verifies excellent scaling behavior and yields a dynamical exponentz _scal which perfectly agrees withz.     

Monte Carlo studies of the square Ising model with next-nearest-neighbor interactions

We apply a new entropic scheme to study the critical behavior of the square-lattice Ising model with nearest- and next-nearest-neighbor antiferromagnetic interactions. Estimates of the present scheme are compared with those of the Metropolis algorithm. We consider interactions in the range where superantiferromagnetic (SAF) order appears at low temperatures. A recent prediction of a first-order transition along a certain range (0.5–1.2) of the interaction ratio (R=J_nnn/J_nn) is examined by generating accurate data for large lattices at a particular value of the ratio (R=1). Our study does not support a first-order transition and a convincing finite-size scaling analysis of the model is presented, yielding accurate estimates for all critical exponents for R=1. The magnetic exponents are found to obey “weak universality” in accordance with a previous conjecture.     

Monte Carlo simulation of three-dimensional dilute Ising model

A multispin coding program for site-diluted Ising models on large simple cubic lattices is described in detail. The spontaneous magnetization is computed as a function of temperature, and the critical temperature as a function of concentration is found to agree well with the data of Marro et al.⁽⁴⁾ and Landau⁽³⁾ for smaller systems.The first successful epsilon expansion seems to be by D. E. Khmelnitskii,ZhETF 68:1960 (1975), English translationSov. Phys. JETP 41:981 (1975); for numerical estimates see K. E. Newman and E. K. Riedel,Phys. Rev. H25:264 (1982), for experiments see R. J. Birgenau, R. A. Cowley, G. Shirane and H. Yoshizawa,J. Stat. Phys. 34:817 (1984).     

Monte Carlo study of the antiferromagnetical Ising model on a centred honeycomb lattice

We use the Monte Carlo method to study an antiferromagnetical Ising spin system on a centred honeycomb lattice, which is composed of two kinds of 1/2 spin particles A and B. There exist two different bond energies J_A-A and J_A-B in this lattice. Our study is focused on how the ratio of J_A-B to J_A-A influences the critical behaviour of this system by analysing the physical quantities, such as the energy, the order parameter, the specific heat, susceptibility, {etc} each as a function of temperature for a given ratio of J_A-B to J_A-A. Using these results together with the finite-size scaling method, we obtain a phase diagram for the ratio J_A-B / J_A-A. This work is helpful for studying the phase transition problem of crystals composed of compounds.     

Monte Carlo study of an Ising gauge system

The phase diagram of the 3-d Ising gauge model with additional ferromagnetic nearest neighbour Ising coupling is explored by Monte Carlo simulations.     

Bond dilution in the 3D Ising model: a Monte Carlo study

We study by Monte Carlo simulations the influence of bond dilution on the three-dimensional Ising model. This paradigmatic model in its pure version displays a second-order phase transition with a positive specific heat critical exponent . According to the Harris criterion disorder should hence lead to a new fixed point characterized by new critical exponents. We have determined the phase diagram of the diluted model, starting from the pure model limit down to the neighbourhood of the percolation threshold. For the estimation of critical exponents, we have first performed a finite-size scaling study, where we concentrated on three different dilutions to check the stability of the disorder fixed point. We emphasize in this work the great influence of the cross-over phenomena between the pure, disorder and percolation fixed points which lead to effective critical exponents dependent on the concentration. In a second set of simulations, the temperature behaviour of physical quantities has been studied in order to characterize the disorder fixed point more accurately. In particular this allowed us to estimate ratios of some critical amplitudes. In accord with previous observations for other models this provides stronger evidence for the existence of the disorder fixed point since the amplitude ratios are more sensitive to the universality class than the critical exponents. Moreover, the question of non-self-averaging at the disorder fixed point is investigated and compared with recent results for the bond-diluted q = 4 Potts model. Overall our numerical results provide evidence that, as expected on theoretical grounds, the critical behaviour of the bond-diluted model is indeed governed by the same universality class as the site-diluted model.Received: 24 February 2004, Published online: 28 May 2004PACS: 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 64.60.Fr Equilibrium properties near critical points, critical exponents - 75.10.Hk Classical spin models     

Monte Carlo simulation of two-dimensional Ising model with logarithmic exchange interaction

Monte Carlo simulation was done to a 2-D Ising model with exchange interaction depending logarithmically on distance. This model is equivalent to a 2-D Coulomb gas.A configuration corresponding to bound pairs of opposite charges appears below T_c, while above T_c pairs are dissolved into free charges.     

A Markov-property Monte Carlo method: One-dimensional Ising model

In this paper we introduce a new Monte Carlo procedure based on the Markov property. This procedure is particularly well suited to massively parallel computation. We illustrate the method on the critical phenomena of the well known one-dimensional Ising model. In the course of this work we found that the autocorrelation time for the Metropolis Monte Carlo algorithm is closely given by the square of the correlation length. We find speedup factors of the order of 1 million for the method as implemented on the CM2 relative to a serial machine. Our procedure gives error estimates which are quite consistent with the observed deviations from the analytically known exact results.     

Series and Monte Carlo study of high-dimensional Ising models

Ising models in high dimensions are used to compare high-temperature series expansions with Monte Carlo simulations. Simulations of the magnetization on four-, six-, and seven-dimensional hypercubic lattices give consistent values of the critical temperature from both equilibrium and nonequilibrium data ford=6 and 7. We tabulate 15 terms of series expansions for the susceptibility for generald and giveJ/k _B T _c =0.092295 (3) and 0.077706 (2) ford=6 and 7. In contrast to five dimensions, where earlier series found nonanalytic scaling corrections, for d=6 and 7 the leading scaling correction may be analytic inT-T _c . In most cases these expansions gave more accurate results than these simulations.     

Monte Carlo simulation of cluster kinetics in three-dimensional Ising model

Simulations in lattices of size 100³ at the critical point of the three-dimensional Glauber kinetic Ising model indicate that the 1935 Becker-Doring equation no longer works there: The growth rates decay in time. These conclusions confirm those in two dimensions.     

Monte Carlo simulation of Ising model on decagonal covering structure

We study the Ising model on a two-dimensional quasilattice developed from the decagonal covering structure. The periodic boundary conditions are applied to a patch of rhombus-like covering pattern. By means of the Monte Carlo simulation and the finite-size scaling analysis the critical temperature is estimated as 2.317±0.002. Two critical exponents are obtained being 1/v=0.992±0.003 and η=0.247±0.002, which are close to the values of the two-dimensional regular lattices as well as the Penrose tilings.     

Monte Carlo simulation of a kinetic Ising model of the glass transition

Results of Monte Carlo experiments for the two-spin facilitated kinetic Ising model on a cubic lattice are presented and compared with a theoretical prediction.