Lack of self-averaging in critical disordered systems

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$, $V_{P}$, 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. Using 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 the energy, magnetization, specific heat and susceptibility are strongly self averaging, that is $V_{P}\sim l^{-d}$ (where $d=2$ is the dimension). At criticality though, the results indicate that the magnetization $M$ and the susceptibility $\chi$ are non self averaging, i.e. $\frac{V_{\chi}}{\chi^{2}}, \frac{V_{M}}{M^{2}}\not \rightarrow 0$. The energy $E$ at criticality is weakly self averaging, that is $V_{E}\sim l^{-y_{v}}$ with $0
Comment: 33 pages, RevTex, 16 figures in tar compressed form included, Submitted to Phys. Rev. E The figures which were missing are now included, in a uuencoded tar compressed form