No Arabic abstract
We consider the three-dimensional $pm J$ model defined on a simple cubic lattice and study its behavior close to the multicritical Nishimori point where the paramagnetic-ferromagnetic, the paramagnetic-glassy, and the ferromagnetic-glassy transition lines meet in the T-p phase diagram (p characterizes the disorder distribution and gives the fraction of ferromagnetic bonds). For this purpose we perform Monte Carlo simulations on cubic lattices of size $Lle 32$ and a finite-size scaling analysis of the numerical results. The magnetic-glassy multicritical point is found at $p^*=0.76820(4)$, along the Nishimori line given by $2p-1={rm Tanh}(J/T)$. We determine the renormalization-group dimensions of the operators that control the renormalization-group flow close to the multicritical point, $y_1 = 1.02(5)$, $y_2 = 0.61(2)$, and the susceptibility exponent $eta = -0.114(3)$. The temperature and crossover exponents are $ u=1/y_2=1.64(5)$ and $phi=y_1/y_2 = 1.67(10)$, respectively. We also investigate the model-A dynamics, obtaining the dynamic critical exponent $z = 5.0(5)$.
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 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.
We report a high-precision finite-size scaling study of the critical behavior of the three-dimensional Ising Edwards-Anderson model (the Ising spin glass). We have thermalized lattices up to L=40 using the Janus dedicated computer. Our analysis takes into account leading-order corrections to scaling. We obtain Tc = 1.1019(29) for the critical temperature, u = 2.562(42) for the thermal exponent, eta = -0.3900(36) for the anomalous dimension and omega = 1.12(10) for the exponent of the leading corrections to scaling. Standard (hyper)scaling relations yield alpha = -5.69(13), beta = 0.782(10) and gamma = 6.13(11). We also compute several universal quantities at Tc.
The existence of an equilibrium glassy phase for charges in a disordered potential with long-range electrostatic interactions has remained controversial for many years. Here we conduct an extensive numerical study of the disorder-temperature phase diagram of the three-dimensional Coulomb glass model using population annealing Monte Carlo to thermalize the system down to extremely low temperatures. Our results strongly suggest that, in addition to a charge order phase, a transition to a glassy phase can be observed, consistent with previous analytical and experimental studies.
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.