Monte Carlo study of the change of critical temperature in a diluted Ising model due to configuration disorder

Monte Carlo simulations on a diluted Ising Hamiltonian were used to obtain the susceptibility of virtual binary samples (conformed by ferromagnetic atoms and non-magnetic atoms) on a bcc lattice. Samples size of L=10L=10 were constructed with different configuration order using random mixtures. The susceptibility curves illustrate that the critical temperature decreases when the disorder in the samples increase. From fittings of the exchange interaction versus the concentration of non-magnetic atoms it was possible to reasonably describe the magnetic phase diagram of the Fe1−qAlqFe1−qAlq alloys.