No Arabic abstract
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.
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 analyze a controversial question about the universality class of the three-dimensional Ising model with long-range-correlated disorder. Whereas both analytical and numerical studies performed so far support an extended Harris criterion (A. Weinrib, B. I. Halperin, Phys. Rev. B 27 (1983) 413) and bring about the new universality class, the numerical values of the critical exponents found so far differ essentially. To resolve this discrepancy we perform extensive Monte Carlo simulations of a 3d Ising magnet with non-magnetic impurities arranged as lines with random orientation. We apply Wolff cluster algorithm accompanied by a histogram reweighting technique and make use of the finite-size scaling to extract the values of critical exponents governing the magnetic phase transition. Our estimates for the exponents differ from the results of the two numerical simulations performed so far and are in favour of a non-trivial dependency of the critical exponents on the peculiarities of long-range correlations decay.
We discuss universal and non-universal critical exponents of a three dimensional Ising system in the presence of weak quenched disorder. Both experimental, computational, and theoretical results are reviewed. Special attention is paid to the results obtained by the field theoretical renormalization group approach. Different renormalization schemes are considered putting emphasis on analysis of divergent series obtained.
The random-field Ising model (RFIM), one of the basic models for quenched disorder, can be studied numerically with the help of efficient ground-state algorithms. In this study, we extend these algorithm by various methods in order to analyze low-energy excitations for the three-dimensional RFIM with Gaussian distributed disorder that appear in the form of clusters of connected spins. We analyze several properties of these clusters. Our results support the validity of the droplet-model description for the RFIM.
We use a non-equilibrium simulation method to study the spin glass transition in three-dimensional Ising spin glasses. The transition point is repeatedly approached at finite velocity $v$ (temperature change versus time) in Monte Carlo simulations starting at a high temperature. The normally problematic critical slowing-down is not hampering this kind of approach, since the system equilibrates quickly at the initial temperature and the slowing-down is merely reflected in the dynamic scaling of the non-equilibrium order parameter with $v$ and the system size. The equilibrium limit does not have to be reached. For the dynamic exponent we obtain $z = 5.85(9)$ for bimodal couplings distribution and $z=6.00(10)$ for the Gaussian case, thus supporting universal dynamic scaling (in contrast to recent claims of non-universal behavior).