We study sample-to-sample fluctuations in a critical two-dimensional Ising model with quenched random ferromagnetic couplings. Using replica calculations in the renormalization group framework we derive explicit expressions for the probability distribution function of the critical internal energy and for the specific heat fluctuations. It is shown that the disorder distribution of internal energies is Gaussian, and the typical sample-to-sample fluctuations as well as the average value scale with the system size L like ∼Llnln(L). In contrast, the specific heat is shown to be self-averaging with a distribution function that tends to a δ peak in the thermodynamic limit L→∞. While previously a lack of self-averaging was found for the free energy, we here obtain results for quantities that are directly measurable in simulations, and implications for measurements in the actual lattice system are discussed.
Abstract In this paper, we review recent results on sample-to-sample fluctuations in a critical Ising model with quenched random ferromagnetic couplings. In particular, in terms of the renormalized replica Ginzburg–Landau Hamiltonian in dimensions D < 4 an explicit expression for the probability distribution function (PDF) of the critical free energy fluctuations is derived. Next, using known fixed-point values for the renormalized coupling parameters the universal curve for such PDF in the dimension D = 3 is obtained. For the specific case of the two-dimensional Ising model, using replica calculations in the renormalization group framework, we derive explicit expressions for the PDF of the critical internal energy and for the specific heat fluctuations. It is shown that the disorder distribution of internal energy is Gaussian, and its typical sample-to-sample fluctuations as well as its average value scale with the system size L like ~ Llnln(L). In contrast, the specific heat is shown to be self-averaging with a distribution function that tends to a δ-peak in the thermodynamic limit L → ∞.
We consider the sample to sample fluctuations that occur in the value of a thermodynamic quantity
P in an ensemble of finite systems with quenched disorder, at equilibrium. The variance of
P,
VP, which characterizes these fluctuations is calculated as a function of the systems' linear size
l, focusing on the behavior at the critical point. The specific model considered is the bond-disordered Ashkin-Teller model on a square lattice [Phys. Rev. 64, 178 (1943)]. Using extensive Monte Carlo simulations, several bond-disordered Ashkin-Teller models were examined, including the bond-disordered Ising model and the bond-disordered four-state Potts model. It was found that far from criticality all thermodynamic quantities which were examined (energy, magnetization, specific heat, susceptibility) are strongly self-averaging, that is ${V}_{P}\ensuremath{\sim}{l}^{\ensuremath{-}d}$ (where
d=2 is the dimension). At criticality though, the results indicate that the magnetization
M and the susceptibility $\ensuremath{\chi}$ are non-self-averaging, i.e., $\frac{{V}_{\ensuremath{\chi}}}{{\ensuremath{\chi}}^{2}}, \frac{{V}_{M}}{{M}^{2}}\ensuremath{\nrightarrow}0$. The energy
E at criticality is clearly weakly self-averaging, that is ${V}_{E}\ensuremath{\sim}{l}^{\ensuremath{-}{y}_{v}}$ with
0<yv<d. Less conclusively, and possibly only as a transient behavior, the specific heat too is found to be weakly self-averaging. A phenomenological theory of finite size scaling for disordered systems is developed, based on physical considerations similar to those leading to the Harris criterion. Its main prediction is that when the specific heat exponent $\ensuremath{\alpha}<0$ ($\ensuremath{\alpha}$ of the disordered model) then, for a quantity
P which scales as ${l}^{\ensuremath{\rho}}$ at criticality, its variance
VP will scale asymptotically as ${l}^{2\ensuremath{\rho}+\frac{\ensuremath{\alpha}}{\ensuremath{\nu}}}$. The theory is not applicable in the asymptotic limit ($l\ensuremath{\rightarrow}\ensuremath{\infty}$) to the bond-disordered Ashkin-Teller model where $\frac{\ensuremath{\alpha}}{\ensuremath{\nu}}={0}_{+}$. Nonetheless in the accessible range of lattice sizes we found very good agreement between the theory and the data for ${V}_{\ensuremath{\chi}}$ and
VE. The theory may also be compatible with the data for the variance of the magnetization
VM and the variance of the specific heat
VC, but evidence for this is less convincing.
We evaluate by Monte Carlo simulations various singular thermodynamic quantities
X for ensembles of quenched random Ising and Ashkin-Teller models. The measurements are taken at
Tc and we study how the distributions
P(X) (and, in particular, their relative squared width,
RX) over the ensemble depend on the system size
l. The Ashkin-Teller model was studied in the regime where bond randomness is irrelevant and we found weak self-averaging; ${R}_{X}\ensuremath{\sim}{l}^{\ensuremath{\alpha}/\ensuremath{\nu}}\ensuremath{\rightarrow}0,$ where $\ensuremath{\alpha}<0$ and $\ensuremath{\nu}$ are the exponents (of the pure model fixed point) governing the transition. For the site-dilute Ising model on a cubic lattice, known to be governed by a random fixed point, we find that
RX tends to a constant, as predicted by Aharony and Harris. We tested whether this constant is universal. However, this constant is different for canonical and grand canonical disorder. We identify the pseudocritical temperature of each sample
i, Tc(i,l), as the temperature at which the susceptibility reaches its maximal value ${\ensuremath{\chi}}^{\mathrm{max}}.$ The distribution of these sample dependent
Tc(i,l) was investigated; we found that its variance scales as $[\ensuremath{\delta}{T}_{c}(l){]}^{2}\ensuremath{\sim}{l}^{\ensuremath{-}2/\ensuremath{\nu}}.$ Our previously proposed finite size scaling ansatz for disordered systems was tested and found to hold. We did observe deviations from a single function, which imply that sample dependent scaling functions are needed. These deviations are, however, relatively small and hence to obtain a fixed statistical error it may be more computationally efficient to measure ${\ensuremath{\chi}}^{\mathrm{max}}$ than the commonly used $\ensuremath{\chi}{(T}_{c}^{\ensuremath{\infty}}).$
We discuss the non-self-averaging phenomena in the critical point of weakly disordered Ising ferromagnet. In terms of the renormalized replica Ginzburg-Landau Hamiltonian in dimensions D<4, we derive an explicit expression for the probability distribution function (PDF) of the critical free-energy fluctuations. In particular, using known fixed-point values for the renormalized coupling parameters, we obtain the universal curve for such PDF in the dimension D=3. It is demonstrated that this function is strongly asymmetric: its left tail is much slower than the right one.
The eight-vertex model solved by Baxter is shown to be equivalent to two Ising models with nearest-neighbor coupling interacting with one another via a four-spin coupling term. The critical properties of the model in the weak-coupling limit are in agreement with the scaling hypothesis. In this limit where $\ensuremath{\alpha}\ensuremath{\rightarrow}0$, the critical indices obey $\frac{\ensuremath{\gamma}}{{\ensuremath{\gamma}}_{0}}=\frac{\ensuremath{\beta}}{{\ensuremath{\beta}}_{0}}=\frac{\ensuremath{\nu}}{{\ensuremath{\nu}}_{0}}=1\ensuremath{-}\frac{1}{2}\ensuremath{\alpha}$, $\frac{\ensuremath{\delta}}{{\ensuremath{\delta}}_{0}}=\frac{\ensuremath{\eta}}{{\ensuremath{\eta}}_{0}}=1$, with the subscripts zero denoting the index values for the ordinary two-dimensional Ising model.