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.     

The Markov property method applied to Ising model calculations

An efficient method of computation for models possessing the Markov property is set out. We apply this method to the two-dimensional ising model and report exact computations for up to 10 by 10 models with periodic boundary conditions. We find that critical-point, finite-size rounding is quite large in the renormalized coupling constant, which is not divergent at the critical point, in contrast to the energy, which is also not divergent and has no rounding there. The difference is traced to the continuity of the energy and the discontinuity of the renormalized coupling constant at the critical point.     

A new multispin coding algorithm for Monte Carlo simulation of the Ising model

We present a new algorithm for Monte Carlo simulation of the Ising model. The usual serial architecture of a computer is exploited in a novel way, enabling parallel but independent calculations to be carried out on as many spins as there are bits in a computer word in each fundamental move. The algorithm enjoys a substantial increase in execution speed over more usual multispin coding algorithms. By its very nature, the algorithm constitutes a design for a special-purpose processor.     

Monte Carlo study of the critical behavior of pure and site-diluted Ising ferro-and ferrimagnets

Monte Carlo simulations are performed for pure and site-diluted Ising ferro- and ferrimagnets on a simple cubic lattice with up to 40³ sites and with impurity concentrationx. For the diluted ferromagnet (x=0.2) the exponent= 0.392±0.03 is definitely larger than the pure model value of=0.304±0.03. In contrast, for ferrimagnetic systems (x=0, 0.1, 0.2) the values appear to be independent ofx and within the error limits consistent with the value for the pure ferromagnet, possibly because the width of the asymptotic random critical regime (or of the crossover regime) is even smaller than in the case of ferromagnets.     

Tests for Monte Carlo renormalization studies on Ising models

In this expanded version of an earlier letter, we consider many computational details that were omitted for want of space. Ford = 2 Ising spins with up to 13 different short-range interactions, we construct the critical surface in the vicinity of (Onsager's) nearest-neighbor (nn) critical point by using the body of the available information on the solvable nn case. We then see if the Monte Carlo renormalization group flows generated from this point do indeed lie on this surface and quantify the errors if they do not.     

Monte Carlo evidence for a cusp in the uniform susceptibility of the site-diluted simple cubic Ising antiferromagnet

Monte Carlo simulations are performed for pure and site-diluted Ising antiferromagnets on a simple cubic lattice with up to 40³ sites and impurity concentrationx=0, 0.2. and 0.5. Forx=0.5 a cusp emerges for the temperature dependence of the uniform susceptibility at the critical temperature which is contrasted with the smooth behavior for a pure antiferromagnet, in agreement with the theoretical prediction of Fishman and Aharony.     

Fast vectorized algorithm for the Monte Carlo simulation of the random field Ising model

An algorithm for the simulation of the 3-dimensional random field Ising model with a binary distribution of the random fields is presented. It uses multi-spin coding and simulates 64 physically different systems simultaneously. On one processor of a Cray YMP it reaches a speed of 184 million spin updates per second. For smaller field strength we present a version of the algorithm that can perform 242 million spin updates per second on the same machine.     

The roughening transition of the three-dimensional Ising interface: A Monte Carlo study

We study the roughening transition of an interface in an Ising system on a 3D simple cubic lattice using a finite-size scaling method. The particular method has recently been proposed and successfully tested for various solid-on-solid models. The basic idea is the matching of the renormalization-groupflow of the interface with that of the exactly solvable body-centered cubic solid-on-solid model. We unambiguously confirm the Kosterlitz-Thouless nature of the roughening transition of the Ising interface. Our result for the inverse transition temperatureK _r=0.40754(5) is almost two orders of magnitude more accurate than the estimate of Mon, Landau, and Stauffer.     

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 simulation of the compensation and critical behaviors of a ferrimagnetic core/shell nanoparticle Ising model

A Monte Carlo simulation has been used to study the magnetic properties and the critical behaviors of a single spherical nanoparticle, consisting of a ferromagnetic core of spins surrounded by a ferromagnetic shell of S=±1, 0 or , spins with antiferromagnetic interface coupling, located on a simple cubic lattice. A number of characteristic phenomena has been found. In particular, the effects of the shell coupling and the interface coupling on both the critical and compensation temperatures are investigated. We have found that, for appropriate values of the system parameters, two compensation temperatures may occur in the present system.     

Monte Carlo simulation of Ising models by multispin coding on a vector computer

Rebbi's efficient multispin coding algorithm for Ising models is combined with the use of the vector computer CDC Cyber 205. A speed of 21.2 million updates per second is reached. This is comparable to that obtained by special- purpose computers.     

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.     

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 study of dynamic phase transition in Ising metamagnet driven by oscillating magnetic field

The dynamical responses of Ising metamagnet (layered antiferromagnet) in the presence of a sinusoidally oscillating magnetic field are studied by Monte Carlo simulation. The time average staggered magnetisation plays the role of dynamic order parameter. A dynamical phase transition was observed and a phase diagram was plotted in the plane formed by field amplitude and temperature. The dynamical phase boundary is observed to shrink inward as the relative antiferromagnetic strength decreases. The results are compared with that obtained from pure ferromagnetic system. The shape of dynamic phase boundary observed to be qualitatively similar to that obtained from previous meanfield calculations.     

Study of the three-dimensional Ising model on film geometry with the cluster Monte Carlo method

Topological properties of clusters are used to extract critical parameters. This method is tested for the bulk properties ofd=2 percolation and thed=2, 3 Ising model. For the latter we obtain an accurate value of the critical temperatureJ/k _B T _c=0.221617(18). In the case of thed=3 Ising model with film geometry the critical value of the surface coupling at the special transitions is determined as J_1c/J=1.5004(20) together with the critical exponents ₁ ^m =0.237(5) and=0.461(15).     

One-dimensional inhomogeneous Ising model: A new approach

We present a new method for the study of a one-dimensional inhomogeneous Ising chain with nonconstant nearest neighbor interactions. The external field required to produce a given magnetization profile is derived exactly. Some properties of the pair direct correlation function are derived. Our findings generalize previous results of Percus.     

Interface motion in a two-dimensional Ising model with a field

We determine by Monte Carlo simulations the width of an interface between the stable phase and the metastable phase in a two-dimensional Ising model with a magnetic field, in the case of nonconversed order parameter (Glauber dynamics). At zero temperature, the width increases ast with–1/3, as predicted by earlier theories. As temperature increases, the value of the effective exponent that we measure decreases toward the value 1/4, which is the value in the absence of magnetic field.     

Statistical and systematic errors in Monte Carlo sampling

We have studied the statistical and systematic errors which arise in Monte Carlo simulations and how the magnitude of these errors depends on the size of the system being examined when a fixed amount of computer time is used. We find that, depending on the degree of self-averaging exhibited by the quantities measured, the statistical errors can increase, decrease, or stay the same as the system size is increased. The systematic underestimation of response functions due to the finite number of measurements made is also studied. We develop a scaling formalism to describe the size dependence of these errors, as well as their dependence on the bin length (size of the statistical sample), both at and away from a phase transition. The formalism is tested using simulations of thed=3 Ising model at the infinite-lattice transition temperature. We show that for a 96×96×96 system noticeable systematic errors (systematic underestimation of response functions) are still present for total run lengths of 10⁶ Monte Carlo steps/site (MCS) with measurements taken at regular intervals of 10 MCS.This paper is dedicated to Jerry Percus on the occasion of his 65th birthday.     

Hysteretic behavior of Fe0.9−qMn0.1Alq-disordered alloys: a Monte Carlo study

In this work we have simulated the Fe_0.9−qMn_0.1Al_q alloy series with Al contents ranging from 10 up to 50 at%, and for several system sizes. In the simulation, the atoms are randomly distributed on a body-centered cubic according to the atomic disorder achieved through quenching techniques for the experimental samples. In computing the thermodynamic quantities such as the magnetization per site as a function of an external applied magnetic field, we have employed a Monte Carlo algorithm based on a Metropolis dynamics implemented with a random site-diluted Ising model. In this model, we have taken into account nearest-neighbor interactions for which both the ferromagnetic (Fe–Fe) and antiferromagnetic (Fe–Mn, Mn–Mn) interactions are present. From the simulation of the hysteresis loops at room temperature, the remanence and the coercive force as a function of the Al concentration have been obtained. Finally, a comparison with the previous experimental data on coercivity obtained by means of vibrating sample magnetometry is also carried out and discussed.