We study the two-dimensional Ising model on networks with quenched topological (connectivity) disorder. In particular, we construct random lattices of constant coordination number and perform large-scale Monte Carlo simulations in order to obtain critical exponents using finite-size scaling relations. We find disorder-dependent effective critical exponents, similar to diluted models, showing thus no clear universal behavior. Considering the very recent results for the two-dimensional Ising model on proximity graphs and the coordination number correlation analysis suggested by Barghathi and Vojta [Phys. Rev. Lett. 113, 120602 (2014)PRLTAO0031-900710.1103/PhysRevLett.113.120602], our results indicate that the planarity and connectedness of the lattice play an important role on deciding whether the phase transition is stable against quenched topological disorder.
The Voronoi construction is ubiquitous across the natural sciences and engineering. In statistical mechanics, however, only its dual, the Delaunay triangulation, has been considered in the investigation of critical phenomena. In this paper we set to fill this gap by studying three prominent systems of classical statistical mechanics, the equilibrium spin-1/2 Ising model, the nonequilibrium contact process, and the conserved stochastic sandpile model on two-dimensional random Voronoi graphs. Particular motivation comes from the fact that these graphs have vertices of constant coordination number, making it possible to isolate topological effects of quenched disorder from node-intrinsic coordination number disorder. Using large-scale numerical simulations and finite-size scaling techniques, we are able to demonstrate that all three systems belong to their respective clean universality classes. Therefore, quenched disorder introduced by the randomness of the lattice is irrelevant and does not influence the character of the phase transitions. We report the critical points to considerable precision and, for the Ising model, also the first correction-to-scaling exponent.
We consider the two-dimensional randomly site diluted Ising model and the random-bond +/-J Ising model (also called the 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.
We study the nonequilibrium phase transition in the two-dimensional contact process on a randomly diluted lattice by means of large-scale Monte Carlo simulations for times up to 10;{10} and system sizes up to 8000x8000 sites. Our data provide strong evidence for the transition being controlled by an exotic infinite-randomness critical point with activated (exponential) dynamical scaling. We calculate the critical exponents of the transition and find them to be universal, i.e., independent of disorder strength. The Griffiths region between the clean and the dirty critical points exhibits power-law dynamical scaling with continuously varying exponents. We discuss the generality of our findings and relate them to a broader theory of rare region effects at phase transitions with quenched disorder. Our results are of importance beyond absorbing state transitions because, according to a strong-disorder renormalization group analysis, our transition belongs to the universality class of the two-dimensional random transverse-field Ising model.
Motivated by recent experiments on phase behavior of systems confined in porous media, we have studied the effect of randomness on the nature of the phase transition in the two-dimensional Potts model. To model the effects of the porous matrix we introduce a random distribution of couplings P(J(ij))=pdelta(J(ij)-J1)+(1-p)delta(J(ij)-J2) in the q state Potts Hamiltonian. An extensive Monte Carlo study is made on this system for q=5. We studied two different cases of disorder (a) J(1)/J(2)-->infinity and p=0.8 and (b) J(1)/J(2)=10 and p=0.5. We observed, in both cases, that the weak first order transition that appears in the pure case, changes to a second-order transition. A finite size scaling analysis shows that the correlation length exponent nu is close to 1 and the best fit to the dependence of the specific heat on system size is logarithmic. This suggests that both cases belong to the universality class of the Ising model. In contrast, the magnetic exponents point to a different universality class.
We report single-cluster Monte Carlo simulations of the Ising model on three-dimensional Poissonian random lattices with up to $128000\ensuremath{\approx}{50}^{3}$ sites which are linked together according to the Voronoi-Delaunay prescription. For each lattice size quenched averages are performed over 96 realizations. By using reweighting techniques and finite-size scaling analyses we investigate the critical properties of the model in the close vicinity of the phase transition point. Our random lattice data provide strong evidence that, for the available system sizes, the resulting effective critical exponents are indistinguishable from recent high-precision estimates obtained in Monte Carlo studies of the Ising model and ${\ensuremath{\varphi}}^{4}$ field theory on three-dimensional regular cubic lattices.
Using the two-dimensional Ashkin-Teller (AT) model, we compare criticality on three non-periodic lattices: the Smith-hat aperiodic tiling, Voronoi-Delaunay (VD) random triangulations, and uncorrelated diluted square lattices. The decay of the block-averaged coordination fluctuation
σQ with exponent
α is used to describe the connectivity disorder. The first two lattices share the same exponent
α, which differs from that of the third. We consider the regime where the correlation-length exponent
ν<1, where randomness is relevant according to the Harris criterion
dν≤2, but should be irrelevant in cases of the Smith-hat tiling and VD triangulations, where
αν>1, according to the Harris--Barghathi--Vojta (HBV) criterion. For the diluted lattice, we indeed find that the clean universality behavior breaks down along the entire critical line, indicating a crossover to a fixed line dominated by disorder, in line with both the Harris criterion and the HBV criterion. In contrast, both VD and Smith-hat lattices display critical exponents consistent with the clean AT universality class, as validated by a Coulomb-gas self-consistency check, violating the Harris criterion while conforming to the HBV criterion.