In the supplemental materials we justify our choice of the number of Chebychev moments used within the kernel polynomial method, show some preliminary results for the large coupling behavior, discuss possible correlation effects in the local density
of states, estimate the spin relaxation length and introduce the goodness of fit probability that is used to assess the quality of the fits.
We study effects of classical magnetic impurities on the Anderson metal-insulator transition numerically. We find that a small concentration of Heisenberg impurities enhances the critical disorder amplitude $W_{rm c}$ with increasing exchange couplin
g strength $J$. The resulting scaling with $J$ is analyzed which supports an anomalous scaling prediction by Wegner due to the combined breaking of time-reversal and spin-rotational symmetry. Moreover, we find that the presence of magnetic impurities lowers the critical correlation length exponent $ u$ and enhances the multifractality parameter $alpha_0$. The new value of $ u$ improves the agreement with the value measured in experiments on the metal-insulator transition (MIT) in doped semiconductors like phosphor-doped silicon, where a finite density of magnetic moments is known to exist in the vicinity of the MIT. The results are obtained by a finite-size scaling analysis of the geometric mean of the local density of states which is calculated by means of the kernel polynomial method. We establish this combination of numerical techniques as a method to obtain critical properties of disordered systems quantitatively.
