No Arabic abstract
For the 2D Ising model, we analyzed dependences of thermodynamic characteristics on number of spins by means of computer simulations. We compared experimental data obtained using the Fisher-Kasteleyn algorithm on a square lattice with $N=l{times}l$ spins and the asymptotic Onsager solution ($Ntoinfty$). We derived empirical expressions for critical parameters as functions of $N$ and generalized the Onsager solution on the case of a finite-size lattice. Our analytical expressions for the free energy and its derivatives (the internal energy, the energy dispersion and the heat capacity) describe accurately the results of computer simulations. We showed that when $N$ increased the heat capacity in the critical point increased as $lnN$. We specified restrictions on the accuracy of the critical temperature due to finite size of our system. Also in the finite-dimensional case, we obtained expressions describing temperature dependences of the magnetization and the correlation length. They are in a good qualitative agreement with the results of computer simulations by means of the dynamic Metropolis Monte Carlo method.
We study critical behavior of the diluted 2D Ising model in the presence of disorder correlations which decay algebraically with distance as $sim r^{-a}$. Mapping the problem onto 2D Dirac fermions with correlated disorder we calculate the critical properties using renormalization group up to two-loop order. We show that beside the Gaussian fixed point the flow equations have a non trivial fixed point which is stable for $0.995<a<2$ and is characterized by the correlation length exponent $ u= 2/a + O((2-a)^3)$. Using bosonization, we also calculate the averaged square of the spin-spin correlation function and find the corresponding critical exponent $eta_2=1/2-(2-a)/4+O((2-a)^2)$.
We consider the two-dimensional randomly site diluted Ising model and the random-bond +-J Ising model (also called Edwards-Anderson model), and study their critical behavior at the paramagnetic-ferromagnetic transition. The critical behavior of thermodynamic quantities can be derived from a set of renormalization-group equations, in which disorder is a marginally irrelevant perturbation at the two-dimensional Ising fixed point. We discuss their solutions, focusing in particular on the universality of the logarithmic corrections arising from the presence of disorder. Then, we present a finite-size scaling analysis of high-statistics Monte Carlo simulations. The numerical results confirm the renormalization-group predictions, and in particular the universality of the logarithmic corrections to the Ising behavior due to quenched dilution.
The phase diagram of the random field Ising model on the Bethe lattice with a symmetric dichotomous random field is closely investigated with respect to the transition between the ferromagnetic and paramagnetic regime. Refining arguments of Bleher, Ruiz and Zagrebnov [J. Stat. Phys. 93, 33 (1998)] an exact upper bound for the existence of a unique paramagnetic phase is found which considerably improves the earlier results. Several numerical estimates of transition lines between a ferromagnetic and a paramagnetic regime are presented. The obtained results do not coincide with a lower bound for the onset of ferromagnetism proposed by Bruinsma [Phys. Rev. B 30, 289 (1984)]. If the latter one proves correct this would hint to a region of coexistence of stable ferromagnetic phases and a stable paramagnetic phase.
We use large-scale Monte Carlo simulations to test the Weinrib-Halperin criterion that predicts new universality classes in the presence of sufficiently slowly decaying power-law-correlated quenched disorder. While new universality classes are reasonably well established, the predicted exponents are controversial. We propose a method of growing such correlated disorder using the three-dimensional Ising model as benchmark systems both for generating disorder and studying the resulting phase transition. Critical equilibrium configurations of a disorder-free system are used to define the two-value distributed random bonds with a small power-law exponent given by the pure Ising exponent. Finite-size scaling analysis shows a new universality class with a single phase transition, but the critical exponents $ u_d=1.13(5), eta_d=0.48(3)$ differ significantly from theoretical predictions. We find that depending on details of the disorder generation, disorder-averaged quantities can develop peaks at two temperatures for finite sizes. Finally, a layer model with the two values of bonds spatially separated to halves of the system genuinely has multiple phase transitions and thermodynamic properties can be flexibly tuned by adjusting the model parameters.
The site-diluted transverse field Ising model in two dimensions is studied with Quantum-Monte-Carlo simulations. Its phase diagram is determined in the transverse field (Gamma) and temperature (T) plane for various (fixed) concentrations (p). The nature of the quantum Griffiths phase at zero temperature is investigated by calculating the distribution of the local zero-frequency susceptibility. It is pointed out that the nature of the Griffiths phase is different for small and large Gamma.