Multiscaling in the 3D critical site-diluted Ising ferromagnet

We study numerically the appearance of multiscaling behavior in the 3D ferromagnetic Ising site-diluted model, in the form of a multifractal distribution of the decay exponents for the spatial correlation functions at the critical temperature. We have computed the exponents of the long-distance decay of higher moments of the correlation function, up to the 10th power, by studying three different quantities: global susceptibilities, local susceptibilities and correlation functions. We have found very clear evidence of multiscaling behavior.