On the mean-field Ising model in a random external field

We use a method developed by van Hemmen to obtain the free energy of the mean-field Ising model in a random external magnetic field. Some results of previous mean-field calculations are confirmed and generalized. The tricritical point in the global phase diagram is discussed in detail. We also consider different probability distributions of the random fields and provide some proofs regarding the conditions for the existence of a tricritical point.     

Fluctuations in the Curie-Weiss version of the random field Ising model

The fluctuations of the order parameter in the Curie-Weiss version of the Ising model with random magnetic field are computed. Away from criticality or at first-order critical points they have a Gaussian distribution with random (i. e.,sample-dependent) mean, thermal fluctuations contributing in same order as the fluctuations of the field; at second- or higher-order critical points, non-Gaussian sample-dependent distributions appear, and the fluctuations of the fields are enhanced, dominating over the thermal ones.     

Absence of tricritical behavior of the random field Ising model in a honyecomb lattice

In this letter, we study the behavior of the random field Ising model on a honeycomb lattice by means of the effective field theory. We obtain the phase diagram in the T$T$–H$H$ plane for clusters with one spin in a finite size cluster scheme and it is observed the absence of a tricritical point.     

Mixed spin Ising model with four-spin interaction and random crystal field

The effects of fluctuations of the crystal field on the phase diagram of the mixed spin-1/2 and spin-1 Ising model with four-spin interactions are investigated within the finite cluster approximation based on a single-site cluster theory. The state equations are derived for the two-dimensional square lattice. It has been found that the system exhibits a variety of interesting features resulting from the fluctuation of the crystal field interactions. In particular, for low mean value D of the crystal field, the critical temperature is not very sensitive to fluctuations and all transitions are of second order for any value of the four-spin interactions. But for relatively high D, the transition temperature depends on the fluctuation of the crystal field, and the system undergoes tricritical behaviour for any strength of the four-spin interactions. We have also found that the model may exhibit reentrance for appropriate values of the system parameters.     

Wang-Landau study of the critical behavior of the bimodal 3D random field Ising model

We apply the Wang-Landau method to the study of the critical behavior of the three-dimensional random field Ising model with a bimodal probability distribution. For high values of the random field intensity we find that the energy probability distribution at the transition temperature is double peaked, suggesting that the phase transition is of first order. On the other hand, the transition looks continuous for low values of the field intensity. In spite of the large sample to sample fluctuations observed, the double peak in the probability distribution is always present for high fields.     

One-dimensional Ising model in arbitrary external field

The external potential needed to produce an arbitrary equilibrium density profile for a one-dimensional lattice gas with nearest neighbor interactions is solved exactly. The resulting sequence of direct correlation functions is shown to be of short range, and in the ferromagnetic case the even members alternate in sign at zero spin. The even Ursell distributions in this case likewise alternate in sign.Supported in part by U.S. ERDA under contract E(11-1)-3077.     

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.     

横向随机晶场Ising模型的相图和磁化行为研究

在有效场理论和切断近似框架内，选择自旋S=1的二维方格子，研究横向随机晶场Ising模型的相图和磁化行为，重点是横向随机晶场浓度和晶场比率对相图和磁化的影响.给出了i>T-D_x空间的相图和m-T空间的磁化图.在晶场稀疏情况下，负晶场方向存在临界温度的峰值，正方向可出现重入现象.晶场比率取+0.5和-0.5时，磁有序相范围缩小，特别是晶场比率取-0.5时，随晶场浓度的降低，临界温度峰值从横向晶场负方向渡越到正方向.固定某一负晶场值，不同晶场比率的磁化行为有明显差异.同时与纵向稀疏晶场Ising模型结果进行有意义的比较. 关键词： 横向随机晶场Ising模型 相图 磁化行为     

Zero-temperature dynamics for the ferromagnetic Ising model on random graphs

We consider Glauber dynamics at zero temperature for the ferromagnetic Ising model on the usual random graph model on N vertices, with on average γ edges incident to each vertex, in the limit as N→∞. Based on numerical simulations, Svenson (Phys. Rev. E 64 (2001) 036122) reported that the dynamics fails to reach a global energy minimum for a range of values of γ. The present paper provides a mathematically rigorous proof that this failure to find the global minimum in fact happens for all γ>0. A lower bound on the residual energy is also given.     

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.     

Rigorous bounds of the Lyapunov exponents of the one-dimensional random Ising model

We find analytic upper and lower bounds of the Lyapunov exponents of the product of random matrices related to the one-dimensional disordered Ising model, using a deterministic map which transforms the original system into a new one with smaller average couplings and magnetic fields. The iteration of the map gives bounds which estimate the Lyapunov exponents with increasing accuracy. We prove, in fact, that both the upper and the lower bounds converge to the Lyapunov exponents in the limit of infinite iterations of the map. A formal expression of the Lyapunov exponents is thus obtained in terms of the limit of a sequence. Our results allow us to introduce a new numerical procedure for the computation of the Lyapunov exponents which has a precision higher than Monte Carlo simulations.     

The random field Ising model with an asymmetric and anisotropic trimodal probability distribution

The Ising model in the presence of a random field, drawn from the asymmetric and anisotropic trimodal probability distribution P(_hi)=pδ(_hi−_h0)+qδ(_hi+λ∗_h0)+rδ(_hi)$P\left(h_i\right)=p\delta (h_i-h_0)+q\delta (h_i+\lambda \ast h_0)+r\delta \left(h_i\right)$, is investigated. The partial probabilities p,q,r$p,q,r$ take on values within the interval [0,1]$[0,1]$ consistent with the constraint p+q+r=1$p+q+r=1$; asymmetric distribution, _hi$h_i$ is the random field variable with basic absolute value _h0$h_0$ (strength); λ$\lambda$ is the competition parameter, which is the ratio between the respective strength of the random magnetic field in the two principal directions (+z)$(+z)$ and (−z)$(-z)$ and is positive so that the random fields are competing, anisotropic distribution. This probability distribution is an extension of the bimodal one allowing for the existence in the lattice of non magnetic particles or vacant sites. The current random field Ising system displays mainly second order phase transitions, which, for some values of p,q$p,q$ and _h0$h_0$, are followed by first order phase transitions joined smoothly by a tricritical point; occasionally, two tricritical points appear implying another second order phase transition. In addition to these points, re-entrant phenomena can be seen for appropriate ranges of the temperature and random field for specific values of λ$\lambda$, p$p$ and q$q$. Using the variational principle, we write down the equilibrium equation for the magnetization and solve it for both phase transitions and at the tricritical point in order to determine the magnetization profile with respect to _h0$h_0$, considered as an independent variable in addition to the temperature.     

The random field Ising model with an asymmetric trimodal probability distribution

The Ising model in the presence of a random field is investigated within the mean field approximation based on Landau expansion. The random field is drawn from the trimodal probability distribution P(h_i)=pδ(h_i−h₀)+qδ(h_i+h₀)+rδ(h_i), where the probabilities p,q,r take on values within the interval [0,1] consistent with the constraint p+q+r=1 (asymmetric distribution), h_i is the random field variable and h₀ the respective strength. This probability distribution is an extension of the bimodal one allowing for the existence in the lattice of non magnetic particles or vacant sites. The current random field Ising system displays second order phase transitions, which, for some values of p,q and h₀, are followed by first order phase transitions, thus confirming the existence of a tricritical point and in some cases two tricritical points. Also, reentrance can be seen for appropriate ranges of the aforementioned variables. Using the variational principle, we determine the equilibrium equation for magnetization, solve it for both transitions and at the tricritical point in order to determine the magnetization profile with respect to h₀.     

Random infinite-volume Gibbs states for the Curie-Weiss random field Ising model

An approach to the definition of infinite-volume Gibbs states for the (quenched) random-field Ising model is considered in the case of a Curie-Weiss ferromagnet. It turns out that these states are random quasi-free measures. They are random convex linear combinations of the free product-measures shifted by the corresponding effective mean fields. The conditional self-averaging property of the magnetization related to this randomness is also discussed.This paper is dedicated to Robert A. Minlos on the occasion of his 60th birthday.     

The hierarchical random field Ising model

We introduce and solve explicitly a hierarchical approximation to the random field Ising model. This approximation is defined in terms of Peierls' contours. It exhibits a spontaneous magnetization ind>2 and illustrates some of the ideas used in the proof of that result for the real RFIM. Ind=2, there is no spontaneous magnetization, but a very slow decay of correlations. However, we argue that this latter property is an artifact of the approximation. For the real RFIM, we expect exponential decay of correlations for any value of the disorder.     

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.     

Ising metamagnet under staggered field

We report a Monte Carlo simulation of the layered Ising antiferromagnet under an external magnetic field. We show that under a staggered field, there occurs a phase transition from a metastable state which follows a Vogel-Fulcher law. For a staggered intrasublattice interaction a similar situation occurs.     

Diagrammatic expansion and metastability in the random-field Ising model

Within the perturbation diagrammatic expansion we discuss the origin of differences in determinations of the lower critical dimension of the random-field Ising model and show that below four dimensions metastability and hysteresis occur. We also explain the occurrence of a quasicritical d=2 behavior at weak random fields, which is responsible for local stability of the ordered state above two dimensions.