Three-dimensional universality class of the Ising model with power-law correlated critical disorder

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 a benchmark system for both 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 nu(d) = 1.13(5), eta(d) = 0.48(3) differ significantly from theoretical predictions. We find that depending on the 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 in halves of the system genuinely has multiple phase transitions, and thermodynamic properties can be flexibly tuned by adjusting the model parameters.