Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Bond dilution in the 3D Ising model: a Monte Carlo study

Pierre-Emmanuel Berche, Christophe Chatelain, Wolfhard Janke

Published Jan 1, 2008• 44 citations
OpenAlex

Abstract

We study by Monte Carlo simulations the influence of bond dilution on the three-dimensional Ising model. This paradigmatic model in its pure version displays a second-order phase transition with a positive specific heat critical exponent $\alpha$ . According to the Harris criterion disorder should hence lead to a new fixed point characterized by new critical exponents. We have determined the phase diagram of the diluted model, starting from the pure model limit down to the neighbourhood of the percolation threshold. For the estimation of critical exponents, we have first performed a finite-size scaling study, where we concentrated on three different dilutions to check the stability of the disorder fixed point. We emphasize in this work the great influence of the cross-over phenomena between the pure, disorder and percolation fixed points which lead to effective critical exponents dependent on the concentration. In a second set of simulations, the temperature behaviour of physical quantities has been studied in order to characterize the disorder fixed point more accurately. In particular this allowed us to estimate ratios of some critical amplitudes. In accord with previous observations for other models this provides stronger evidence for the existence of the disorder fixed point since the amplitude ratios are more sensitive to the universality class than the critical exponents. Moreover, the question of non-self-averaging at the disorder fixed point is investigated and compared with recent results for the bond-diluted q = 4 Potts model. Overall our numerical results provide evidence that, as expected on theoretical grounds, the critical behaviour of the bond-diluted model is indeed governed by the same universality class as the site-diluted model.

Related papers

High-temperature series for the bond-diluted Ising model in 3, 4, and 5 dimensions

Meik Hellmund, Wolfhard Janke E

In order to study the influence of quenched disorder on second-order phase transitions, high-temperature series expansions of the susceptibility and the free energy are obtained for the quenched bond-diluted Ising model in $d=3--5$ dimensions. They are analyzed using different extrapolation methods tailored to the expected singularity behaviors. In $d=4$ and 5 dimensions we confirm that the critical behavior is governed by the pure fixed point up to dilutions near the geometric bond percolation threshold. The existence and form of logarithmic corrections for the pure Ising model in $d=4$ are confirmed, and our results for the critical behavior of the diluted system are in agreement with the type of singularity predicted by renormalization group considerations. In three dimensions we find large crossover effects between the pure Ising, percolation, and random fixed points. We estimate the critical exponent of the susceptibility to be $\ensuremath{\gamma}=1.305(5)$ at the random fixed point.

arXiv:cond-mat/0606320•Physical Review B•Oct 12, 2006•13 citations

The universality class of 3D site-diluted and bond-diluted Ising systems

Martin Hasenbusch B, Francesco Parisen Toldin, Andrea Pelissetto C +1 more

We present a finite-size scaling analysis of high-statistics Monte Carlo simulations of the three-dimensional randomly site-diluted and bond-diluted Ising model. The critical behaviour of these systems is affected by slowly decaying scaling corrections which make the accurate determination of their universal asymptotic behaviour quite hard, requiring an effective control of the scaling corrections. For this purpose we exploit improved Hamiltonians, for which the leading scaling corrections are suppressed for any thermodynamic quantity, and improved observables, for which the leading scaling corrections are suppressed for any model belonging to the same universality class. The results of the finite-size scaling analysis provide strong numerical evidence that phase transitions in three-dimensional randomly site-diluted and bond-diluted Ising models belong to the same randomly dilute Ising universality class. We obtain accurate estimates of the critical exponents, ν = 0.683(2), η = 0.036(1), α = −0.049(6), γ = 1.341(4), β = 0.354(1), δ = 4.792(6), and of the leading and next-to-leading correction-to-scaling exponents, ω = 0.33(3) and ω 2 = 0.82(8).

arXiv:cond-mat/0611707•Journal of Statistical Mechanics Theory and Experiment•Feb 21, 2007•68 citations

Relaxational dynamics in 3D randomly diluted Ising models

Martin Hasenbusch B, Andrea Pelissetto C, Ettore Vicari

We study the purely relaxational dynamics (model A) at criticality in three-dimensional disordered Ising systems whose static critical behaviour belongs to the randomly diluted Ising universality class. We consider the site-diluted and bond-diluted Ising models, and the ± J Ising model along the paramagnetic–ferromagnetic transition line. We perform Monte Carlo simulations at the critical point using the Metropolis algorithm and study the dynamic behaviour in equilibrium at various values of the disorder parameter. The results provide a robust evidence of the existence of a unique model-A dynamic universality class which describes the relaxational critical dynamics in all considered models. In particular, the analysis of the size dependence of suitably defined autocorrelation times at the critical point provides the estimate z = 2.35(2) for the universal dynamic critical exponent. We also study the off-equilibrium relaxational dynamics following a quench from to T c . In agreement with the field-theory scenario, the analysis of the off-equilibrium dynamic critical behaviour gives an estimate of z that is perfectly consistent with the equilibrium estimate z = 2.35(2).

arXiv:0709.4179•Journal of Statistical Mechanics Theory and Experiment•Nov 23, 2007•23 citations

Monte Carlo Studies of Three-Dimensional Bond-Diluted Ferromagnets

P. -E. Berche, Christophe Chatelain, Bertrand Berche +1 more

In this report we give an overview on recent results obtained from extensive Monte Carlo (MC) computer simulations of the 3D 2-state (Ising) and 4-state Potts models with bond-dilution. The motivation to study the 4-state Potts model derives from the fact that, in the pure case, this model is known to exhibit a fairly strong first-order transition, such that a disorder-induced softening to a second-order transition would give clear support of the theoretical expectations. Modeling the disorder by bond-dilution enables in the Ising case a test of the expected universality with respect to the type of disorder. Furthermore, for both models this choice facilitates a quantitative comparison with recent high-temperature series expansions for general random-bond q-state Potts models.

arXiv:cond-mat/0212504•Dec 19, 2002•1 citation

Self-averaging in the three-dimensional site diluted Heisenberg model at the critical point

A. Gordillo-Guerrero, Juan Jesus Ruiz-Lorenzo

We study the self-averaging properties of the three-dimensional site diluted Heisenberg model. The Harris criterion states that disorder is irrelevant since the specific heat critical exponent of the pure model is negative. According to some analytical approaches, this implies that the susceptibility should be self-averaging at the critical temperature ( R χ = 0). We have checked this theoretical prediction for a large range of dilution (including strong dilution) at criticality and we have found that the introduction of scaling corrections is crucial for obtaining ‘self-averageness’ in this model. Finally we have computed critical exponents and cumulants which compare very well with those of the pure model, supporting the universality predicted by the Harris criterion.

arXiv:cond-mat/0703820•Journal of Statistical Mechanics Theory and Experiment•Jun 20, 2007•15 citations

Three-dimensional randomly dilute Ising model: Monte Carlo results

Pasquale Calabrese, Víctor Martín-Mayor, Andrea Pelissetto C +1 more

We perform a high-statistics simulation of the three-dimensional randomly dilute Ising model on cubic lattices L3 with L< or =256. We choose a particular value of the density, x=0.8, for which the leading scaling corrections are suppressed. We determine the critical exponents, obtaining nu=0.683(3), eta=0.035(2), beta=0.3535(17), and alpha=-0.049(9), in agreement with previous numerical simulations. We also estimate numerically the fixed-point values of the four-point zero-momentum couplings that are used in field-theoretical fixed-dimension studies. Although these results somewhat differ from those obtained using perturbative field theory, the field-theoretical estimates of the critical exponents do not change significantly if the Monte Carlo result for the fixed point is used. Finally, we determine the six-point zero-momentum couplings, relevant for the small-magnetization expansion of the equation of state, and the invariant amplitude ratio R(+)(xi) that expresses the universality of the free-energy density per correlation volume. We find R(+)(xi)=0.2885(15).

arXiv:cond-mat/0306272•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Sep 30, 2003•53 citations
We study by Monte Carlo simulations the influence of bond dilution on the three-dimensional Ising model. This paradigmatic model in its pure version displays a second-order phase transition with a positive specific heat critical exponent α\alphaα . According to the Harris criterion disorder should hence lead to a new fixed point characterized by new critical exponents. We have determined the phase diagram of the diluted model, starting from the pure model limit down to the neighbourhood of the percolation threshold. For the estimation of critical exponents, we have first performed a finite-size scaling study, where we concentrated on three different dilutions to check the stability of the disorder fixed point. We emphasize in this work the great influence of the cross-over phenomena between the pure, disorder and percolation fixed points which lead to effective critical exponents dependent on the concentration. In a second set of simulations, the temperature behaviour of physical quantities has been studied in order to characterize the disorder fixed point more accurately. In particular this allowed us to estimate ratios of some critical amplitudes. In accord with previous observations for other models this provides stronger evidence for the existence of the disorder fixed point since the amplitude ratios are more sensitive to the universality class than the critical exponents. Moreover, the question of non-self-averaging at the disorder fixed point is investigated and compared with recent results for the bond-diluted q = 4 Potts model. Overall our numerical results provide evidence that, as expected on theoretical grounds, the critical behaviour of the bond-diluted model is indeed governed by the same universality class as the site-diluted model.
In order to study the influence of quenched disorder on second-order phase transitions, high-temperature series expansions of the susceptibility and the free energy are obtained for the quenched bond-diluted Ising model in d=3−−5d=3--5d=3−−5 dimensions. They are analyzed using different extrapolation methods tailored to the expected singularity behaviors. In d=4d=4d=4 and 5 dimensions we confirm that the critical behavior is governed by the pure fixed point up to dilutions near the geometric bond percolation threshold. The existence and form of logarithmic corrections for the pure Ising model in d=4d=4d=4 are confirmed, and our results for the critical behavior of the diluted system are in agreement with the type of singularity predicted by renormalization group considerations. In three dimensions we find large crossover effects between the pure Ising, percolation, and random fixed points. We estimate the critical exponent of the susceptibility to be $\ensuremath{\gamma}=1.305(5)$ at the random fixed point.
We present a finite-size scaling analysis of high-statistics Monte Carlo simulations of the three-dimensional randomly site-diluted and bond-diluted Ising model. The critical behaviour of these systems is affected by slowly decaying scaling corrections which make the accurate determination of their universal asymptotic behaviour quite hard, requiring an effective control of the scaling corrections. For this purpose we exploit improved Hamiltonians, for which the leading scaling corrections are suppressed for any thermodynamic quantity, and improved observables, for which the leading scaling corrections are suppressed for any model belonging to the same universality class. The results of the finite-size scaling analysis provide strong numerical evidence that phase transitions in three-dimensional randomly site-diluted and bond-diluted Ising models belong to the same randomly dilute Ising universality class. We obtain accurate estimates of the critical exponents, ν = 0.683(2), η = 0.036(1), α = −0.049(6), γ = 1.341(4), β = 0.354(1), δ = 4.792(6), and of the leading and next-to-leading correction-to-scaling exponents, ω = 0.33(3) and ω 2 = 0.82(8).
We study the purely relaxational dynamics (model A) at criticality in three-dimensional disordered Ising systems whose static critical behaviour belongs to the randomly diluted Ising universality class. We consider the site-diluted and bond-diluted Ising models, and the ± J Ising model along the paramagnetic–ferromagnetic transition line. We perform Monte Carlo simulations at the critical point using the Metropolis algorithm and study the dynamic behaviour in equilibrium at various values of the disorder parameter. The results provide a robust evidence of the existence of a unique model-A dynamic universality class which describes the relaxational critical dynamics in all considered models. In particular, the analysis of the size dependence of suitably defined autocorrelation times at the critical point provides the estimate z = 2.35(2) for the universal dynamic critical exponent. We also study the off-equilibrium relaxational dynamics following a quench from to T c . In agreement with the field-theory scenario, the analysis of the off-equilibrium dynamic critical behaviour gives an estimate of z that is perfectly consistent with the equilibrium estimate z = 2.35(2).
In this report we give an overview on recent results obtained from extensive Monte Carlo (MC) computer simulations of the 3D 2-state (Ising) and 4-state Potts models with bond-dilution. The motivation to study the 4-state Potts model derives from the fact that, in the pure case, this model is known to exhibit a fairly strong first-order transition, such that a disorder-induced softening to a second-order transition would give clear support of the theoretical expectations. Modeling the disorder by bond-dilution enables in the Ising case a test of the expected universality with respect to the type of disorder. Furthermore, for both models this choice facilitates a quantitative comparison with recent high-temperature series expansions for general random-bond q-state Potts models.
We study the self-averaging properties of the three-dimensional site diluted Heisenberg model. The Harris criterion states that disorder is irrelevant since the specific heat critical exponent of the pure model is negative. According to some analytical approaches, this implies that the susceptibility should be self-averaging at the critical temperature ( R χ = 0). We have checked this theoretical prediction for a large range of dilution (including strong dilution) at criticality and we have found that the introduction of scaling corrections is crucial for obtaining ‘self-averageness’ in this model. Finally we have computed critical exponents and cumulants which compare very well with those of the pure model, supporting the universality predicted by the Harris criterion.
We perform a high-statistics simulation of the three-dimensional randomly dilute Ising model on cubic lattices L3 with L< or =256. We choose a particular value of the density, x=0.8, for which the leading scaling corrections are suppressed. We determine the critical exponents, obtaining nu=0.683(3), eta=0.035(2), beta=0.3535(17), and alpha=-0.049(9), in agreement with previous numerical simulations. We also estimate numerically the fixed-point values of the four-point zero-momentum couplings that are used in field-theoretical fixed-dimension studies. Although these results somewhat differ from those obtained using perturbative field theory, the field-theoretical estimates of the critical exponents do not change significantly if the Monte Carlo result for the fixed point is used. Finally, we determine the six-point zero-momentum couplings, relevant for the small-magnetization expansion of the equation of state, and the invariant amplitude ratio R(+)(xi) that expresses the universality of the free-energy density per correlation volume. We find R(+)(xi)=0.2885(15).