An On-Line Library of Extended High-Temperature Expansions of Basic Observables for the Spin-S Ising Models on Two- and Three-Dimensional Lattices

We present an on-line library of unprecedented extension for high-temperature expansions of basic observables in the Ising models of general spin S, with nearest-neighbor interactions. We have tabulated through order ²⁵ the series for the nearest-neighbor correlation function, the susceptibility and the second correlation moment in two dimensions on the square lattice, and, in three dimensions, on the simple-cubic and the body-centered cubic lattices. The expansion of the second field derivative of the susceptibility is also tabulated through ²³ for the same lattices. We have thus added several terms (from four up to thirteen) to the series already published for spin S = 1/2, 1, 3/2, 2, 5/2, 3, 7/2, 4, 5, .     

A Solvable Decorated Ising Lattice Model

A decorated lattice is suggested and the Ising model on it with three kinds of interactions K₁, K₂, and K₃ is studied. Using an equivalent transformation, the square decorated Ising lattice is transformed into a regular square Ising lattice with nearest-neighbor, next-nearest-neighbor, and four-spin interactions, and the critical fixed point is found at K₁=0.5769, K₂=-0.0671, and K₃=0.3428, which determines the critical temperature of the system. It is also found that this system and the regular square Ising lattice, and the eight-vertex model belong to the same universality class.     

Potts models and random-cluster processes with many-body interactions

Known differential inequalities for certain ferromagnetic Potts models with pair interactions may be extended to Potts models with many-body interactions. As a major application of such differential inequalities, we obtain necessary and sufficient conditions on the set of interactions of such a Potts model in order that its critical point be astrictly monotonic function of the strengths of interactions. The method yields some ancillary information concerning the equality of certain critical exponents for Potts models; this amounts to a small amount of rigorous universality. These results are achieved in the context of a Fortuin-Kasteleyn representation of Potts models with many-body interactions. For such a Potts model, the corresponding random-cluster process is a (random) hypergraph.     

Computation of critical exponents for two-dimensional ising model on a cellular automaton

Static critical exponents for the two-dimensional Ising model are computed on a cellular automaton. The analysis of the data within the framework of the finite-size scaling theory reproduces their well-established values.     

The Percolation Transition for the Zero-Temperature Stochastic Ising Model on the Hexagonal Lattice

On the planar hexagonal lattice , we analyze the Markov process whose state (t), in , updates each site v asynchronously in continuous time t0, so that _v (t) agrees with a majority of its (three) neighbors. The initial _v (0)'s are i.i.d. with P[ _v (0)=+1]=p[0,1]. We study, both rigorously and by Monte Carlo simulation, the existence and nature of the percolation transition as t and p1/2. Denoting by ⁺(t,p) the expected size of the plus cluster containing the origin, we (1) prove that ⁺(,1/2)= and (2) study numerically critical exponents associated with the divergence of ⁺(,p) as p1/2. A detailed finite-size scaling analysis suggests that the exponents and of this t= (dependent) percolation model have the same values, 4/3 and 43/18, as standard two-dimensional independent percolation. We also present numerical evidence that the rate at which (t)() as t is exponential.     

On the critical behavior of the magnetization in high-dimensional Ising models

We derive rigorously general results on the critical behavior of the magnetization in Ising models, as a function of the temperature and the external field. For the nearest-neighbor models it is shown that ind4 dimensions the magnetization is continuous atT _c and its critical exponents take the classical values=3 and=1/2, with possible logarithmic corrections atd=4. The continuity, and other explicit bounds, formally extend tod>3 1/2. Other systems to which the results apply include long-range models ind=1 dimension, with 1/|x–y| couplings, for which 2/(–1) replacesd in the above summary. The results are obtained by means of differential inequalities derived here using the random current representation, which is discussed in detail for the case of a nonvanishing magnetic field.Research supported in part by NSF grant PHY-8301493 A02, and by a John S. Guggenheim Foundation fellowship (M.A.).     

A numerical study of the hierarchical Ising model: High-temperature versus epsilon expansion

We study numerically the magnetic susceptibility of the hierarchical model with Ising spins (=±1) above the critical temperature and for two values of the epsilon parameter. The integrations are performed exactly using recursive methods which exploit the symmetries of the model. Lattices with up to 2¹⁸ sites have been used. Surprisingly, the numerical data can be fitted very well with a simple power law of the form (1-/ ₀)^g for thewhole temperature range considered. This approximate law implies a simple approximate formula for the coefficients of the high-temperature expansion, and, more importantly, approximate relations among the coefficients themselves. We found that some of these approximate relations hold with errors less then 2%. On the other hand,g differs significantly from the critical exponent calculated with the epsilon expansion, even when the fit is restricted to intervals closer to ^c. We discuss this discrepancy in the context the renormalization group analysis of the hierarchical model.     

Fast algorithm for the simulation of Ising models

Using a new microcanonical algorithm efficiently vectorized on a Cray XMP, we reach a simulation speed of 1.5 nsec per update of one spin, three times faster than the best previous method known to us. Data for the nonlinear relaxation with conserved energy are presented for the two-dimensional Ising model.     

High-temperature series analysis of the p-state Potts glass model on d-dimensional hypercubic lattices

We analyze recently extended high-temperature series expansions for the “Edwards-Anderson” spin-glass susceptibility of the p-state Potts glass model on d-dimensional hypercubic lattices for the case of a symmetric bimodal distribution of ferro- and antiferromagnetic nearest-neighbor couplings . In these star-graph expansions up to order 22 in the inverse temperature , the number of Potts states p and the dimension d are kept as free parameters which can take any value. By applying several series analysis techniques to the new series expansions, this enabled us to determine the critical coupling K_c and the critical exponent of the spin-glass susceptibility in a large region of the two-dimensional (p,d)-parameter space. We discuss the thus obtained information with emphasis on the lower and upper critical dimensions of the model and present a careful comparison with previous estimates for special values of p and d. Received: 25 May 1998 / Revised and Accepted: 11 August 1998     

镶嵌正方晶格上具有次近邻耦合作用Ising模型的临界性质

采用部分格点自旋消约变换,将镶嵌正方晶格上具有最近邻耦合作用K1和次近邻耦合作用K2的Ising模型变换成等效的具有最近邻、次近邻和四体耦合作用的正方Ising晶格.发现系统的临界点在(K1C,K2C)=(0.5125,0.2134),由此决定系统的临界温度,幷讨论了系统的普适性.     

The one-dimensional kinetic Ising model: A series expansion study

Exact power series expansions (through eight terms) in the time are derived for relaxation in the one-dimensional Ising model with nearest-neighbor interactions for a general rate parameter where the activation energy is a variable fraction of the energy required to break nearest-neighbor bonds. It is found that the qualitative nature of the relaxation is very dependent on this parameter, varying from nearly simple exponential decay (as with Glauber dynamics) for an intermediate value of this parameter, to an initial rate of change that is either much slower or faster than a simple exponential at the extremes of the range of variation of the parameter. The rate equations for the limit of rapid internal diffusion (internal equilibration) are integrated for several special values of the rate parameter. In general the internal equilibration approximation is not a good representation of the relaxation except when the relaxation is similar to Glauber dynamics.     

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).     

On the equivalence of different order parameters for Ising ferromagnets

In a recent note Barber showed, for a spin-1/2 Ising system with ferromagnetic pair interactions, that some critical exponents of the triplet order parameter _{i j k} are the same as those of the magnetization _i . Here we prove such results for all odd correlations and dispense with the requirement of pair interactions. We also prove that the critical temperatureT _c , defined as the temperature below which there is a spontaneous magnetization, is for fixed even spin interactionsJ _e independent of the way in which the odd interactionsJ _o approach zero from above. This is achieved by using only the simplest, Griffiths-Kelley-Sherman (GKS), inequalities, which apply to the most general many-spin, ferromagnetic interactions.Research supported in part by NSF Grant #MPS 75-20638.     

On the corrections to scaling in three-dimensional Ising models

The leading correction-to-scaling amplitudes for the spin-1/2, nearest-neighbor sc, bcc, and fee Ising models are considered with the particular aim of determining their signs. On the basis of previous two-variable series analyses by Chen, Fisher, and Nickel and renormalization group=4–d expansions, it is concluded that the correction amplitudes for the susceptibility, correlation length, specific heat, and spontaneous magnetization arenegative for all three lattices. Thus, for example, the effective exponent _eff(T) asymptotically approaches the true susceptibility exponent fromabove. Other earlier and more recent high-temperature series and field-theoretic analyses are seen to be consistent with this result. However, the usual nonasymptotic, perturbative field-theoretic approaches are essentially committed to positive correction amplitudes. The question of the signs therefore relates directly to the applicability of these non-asymptotic field-theoretic calculations to three-dimensional Ising models as well as to different experimental systems.     

Ising Model on an Infinite Ladder Lattice

In this paper we propose an Ising model on an infinite ladder lattice, which is made of two infinite Ising spin chains with interactions. It is essentially a quasi-one-dimessional Ising model because the length of the ladder lattice is infinite, while its width is finite. We investigate the phase transition and dynamic behavior of Ising model on this quasi-one-dimessional system. We use the generalized transfer matrix method to investigate the phase transition of the system. It is found that there is no nonzero temperature phase transition in this system. At the same time, we are interested in Glauber dynamics. Based on that, we obtain the time evolution of the local spin magnetization by exactly solving a set of master equations.     

A Mixed-Transfer-Matrix Method for Simulating Normal Conductor/Perfect Insulator/Perfect Conductor Random Networks

We develop a mixed-transfer-matrix approach for computing the macroscopic conductivity of a three-constituent normal conductor/perfect insulator/perfect conductor random network. This is applied to two-dimensional and three-dimensional samples at a percolation threshold. Such networks are simulated in order to test whether a diluted percolating network of normal conducting bonds remains in the same universality class of critical behavior when a finite fraction of those bonds have been replaced by perfectly conducting bonds. Also tested by such simulations is whether a percolating mixture of normal and perfectly conducting bonds remains in the same universality class of critical behavior when a finite fraction of the normal bonds are replaced by perfectly insulating bonds. These questions are crucial for some recently published exact results which connect the macroscopic electrical and elastic responses of percolating networks.     

非平衡相变的临界标度理论及普适性

综述了作者近年来在非平衡相变临界（ Ｎ Ｐ Ｃ） 标度理论及普适性研究的进展。主要包括一般 Ｎ Ｐ Ｃ 系统规格化模型,局域序参量的概率分布,广义势的临界渐近形式,空时有关函数及其临界奇异行为。论证了 Ｎ Ｐ Ｃ 系统的临界可标度性,导出了一组普适的 Ｎ Ｐ Ｃ 标度关系,由之计算出的４ 种 Ｎ Ｐ Ｃ 普适类的临界指数与目前已知的实验及理论结果吻合得非常好。此外,还讨论了非平衡相变临界标度理论的普适性,将平衡相变临界标度理论作为一种特殊极限情况含于同一理论体系中。     

Ising Model on an Infinite Ladder Lattice

In this paper we propose an Ising model on an infinite ladder lattice, which is made of two infinite Ising spin chains with interactions. It is essentially a quasi-one-dimessional Ising model because the length of the ladder lattice is infinite, while its width is finite. We investigate the phase transition and dynamic behavior of Ising model on this quasi-one-dimessional system. We use the generalized transfer matrix method to investigate the phase transition of the system. It is found that there is no nonzero temperature phase transition in this system. At the same time, we are interested in Glauber dynamics. Based on that, we obtain the time evolution of the local spin magnetization by exactly solving a set of master equations.