The nearest-neighbor-interaction spin-1 Ising model is investigated within the damage-spreading approach. Exact relations involving quantities computable through damage-spreading simulations and thermodynamic properties are derived for such a model, defined in terms of a very general Hamiltonian that covers several spin-1 models of interest in the literature. Such relations presuppose translational invariance and hold for any ergodic dynamical procedure, leading to an efficient tool for obtaining thermodynamic properties. The implementation of the method is illustrated through damage-spreading simulations for the ferromagnetic spin-1 Ising model on a square lattice. The two-spin correlation function and the magnetization are obtained, with precise estimates of their associated critical exponents and of the critical temperature of the model, in spite of the small lattice sizes considered. These results are in good agreement with the universality hypothesis, with critical exponents in the same universality class of the spin- 12 Ising model. The advantage of the present method is shown through a significant reduction of finite-size effects by comparing its results with those obtained from standard Monte Carlo simulations.
The nearest-neighbor-interaction spin-1 Ising model is investigated within the damage-spreading approach. Exact relations involving quantities computable through damage-spreading simulations and thermodynamic properties are derived for such a model, defined in terms of a very general Hamiltonian that covers several spin-1 models of interest in the literature. Such relations presuppose translational invariance and hold for any ergodic dynamical procedure, leading to an efficient tool for obtaining thermodynamic properties. The implementation of the method is illustrated through damage-spreading simulations for the ferromagnetic spin-1 Ising model on a square lattice. The two-spin correlation function and the magnetization are obtained, with precise estimates of their associated critical exponents and of the critical temperature of the model, in spite of the small lattice sizes considered. These results are in good agreement with the universality hypothesis, with critical exponents in the same universality class of the spin- 12 Ising model. The advantage of the present method is shown through a significant reduction of finite-size effects by comparing its results with those obtained from standard Monte Carlo simulations.
The nearest-neighbor-interaction ferromagnetic Ashkin-Teller model is investigated on a square lattice through a powerful computational method for dealing with correlation functions in magnetic systems. This technique, which is based on damage-spreading numerical simulations, makes use of exact relations involving special kinds of damage and correlation functions, as well as the corresponding order parameters of the model. The computation of correlation functions, which represents usually a hard task in standard Monte Carlo simulations, due to large fluctuations, turns out to be much simpler within the present approach. We concentrate our analysis along the Baxter line, well known for its continuously varying critical exponents; seven different points along this line are investigated. The critical exponents associated with correlation functions along the Baxter line are successfully evaluated, by means of numerical methods, within damage-spreading simulations. The efficiency of this method is confirmed through precise estimates of the critical exponents associated with the order parameters (magnetization and polarization), as well as with their corresponding correlation functions, in spite of the small lattice sizes considered.
We investigate the damage or Hamming distance between two configurations of Ising spins. We find an exact relation between the difference of the two possible types of damage and the spin-spin correlation function, which is generally valid. For the specific case of ferromagnetic interactions, heat bath dynamics and same sequence of random numbers, this relation involves only one type of damage. The susceptibility and the magnetization can also be expressed in terms of the damage. Numerical determination of the damage for the 2 d Ising model is not only an efficient way to calculate correlation functions but also gives access to spin fluctuations visualized as clusters of damaged sites which have a fractal dimension d −β/ν at T c and whose size distribution is also related to static exponents.
A powerful computational method for dealing with correlation functions in magnetic systems, based on damage-spreading simulations, is reviewed and tested, by investigating the q-state Potts ferromagnet, on a square lattice, at criticality. Exact relations involving special kinds of damage and the spin-spin correlation function, as well as the magnetization, are used. The efficiency of the method arises with a significant reduction of the finite-size effects, with respect to conventional Monte Carlo simulations. Correlation functions, which represent usually a hard task within this latter procedure, appear to be much more easily estimated through the present damage-spreading simulations. The effectiveness of the technique is illustrated by an accurate estimate of the exponent eta, of the spin-spin correlation function, for q=2, 3, and 4, with rather small lattice sizes. In the cases q > or = 5, an analysis of the magnetization is consistent with the well-known first-order phase transition.
Some exact results are derived relating quantities computable by the so-called damage spreading method and thermodynamic functions for the N -color Ashkin–Teller model. The results are valid for any ergodic dynamics. Since we restrict our analysis to the ferromagnetic case, the results are also valid for any translationally invariant lattice. The relations derived can be used in order to determine the N -color Ashkin–Teller critical exponents numerically, providing better accuracy and requiring less computational effort than standard Monte Carlo simulations.
For various Ising models two approaches are discussed, one is that of simulating lattices, also called gauging on exact equations, the other is that of calculating analytical expressions for the boundary free energy of Ising lattices. The first approach allows to conjecture a solution for some Ising models, that have sofar not been solved, once some exact partial result for the problem is known. The second approach aims at furnishing such a partial result in the form of a condition for the critical temperature. An example of such a result was recently given for the 2D Ising square lattice with nearest and next-nearest-neighbor interactions. The critical line that separates the ordered (ferromagnetic) phase from the disordered (paramagnetic) phase showed good agreement in the moderate and strong nearest neighbor coupling limit with several results obtained by Monte Carlo, transfer matrix and series expansion results. We extend the discussion of the critical line, finding an excellent fit, now also in other points, like the Padé point, as well as cusp behavior at the Onsager point where the lattice decouples into two 2D square lattices with only nearest-neighbor interaction. Combination of this result with a geometrical argument in the simulation approach leads to a critical exponent 2−2≈0.5858, comparable to the exponent 4/7≈0.5714 found from renormalization arguments.
The short-time critical dynamics of propagation of damage in the Ising ferromagnet in two dimensions is studied by means of Monte Carlo simulations. Starting with equilibrium configurations at $T=\ensuremath{\infty}$ and magnetization M=0, an initial damage is created by flipping a small amount of spins in one of the two replicas studied. In this way, the initial damage is proportional to the initial magnetization M0 in one of the configurations upon quenching the system at TC, the Onsager critical temperature of the ferromagnetic-paramagnetic transition. It is found that, at short times, the damage increases with an exponent ${\ensuremath{\theta}}_{D}=1.915(3)$, which is much larger than the exponent $\ensuremath{\theta}=0.197$ characteristic of the initial increase of the magnetization M(t). Also, an epidemic study was performed. It is found that the average distance from the origin of the epidemic (⟨R2(t)⟩) grows with an exponent ${z}^{\ensuremath{\ast}}\ensuremath{\approx}\ensuremath{\eta}\ensuremath{\approx}1.9$, which is the same, within error bars, as the exponent ${\ensuremath{\theta}}_{D}$. However, the survival probability of the epidemics reaches a plateau so that $\ensuremath{\delta}=0$. On the other hand, by quenching the system to lower temperatures one observes the critical spreading of the damage at ${T}_{D}\ensuremath{\simeq}0.51{T}_{C}$, where all the measured observables exhibit power laws with exponents ${\ensuremath{\theta}}_{D}=1.026(3)$, $\ensuremath{\delta}=0.133(1)$, and ${z}^{\ensuremath{\ast}}=1.74(3)$.