Monte Carlo simulations of the site-diluted three-dimensional XY model

We present extensive Monte Carlo simulations of the classical XY model in three dimensions on a simple cubic lattice with periodic boundary conditions subject to quenched, random site dilution. By using a hybrid algorithm and finite-size scaling techniques we estimate the transition temperatures and the critical exponents associated with the diluted system. Our results for the critical exponents and universal cumulants are, within statistical errors, the same as for the zero dilution case and are independent of the amount of dilution, in agreement with the Harris criterion. The initial reduction of the critical temperature with dilution is also evaluated and compared to recent experimental results.