Application of polynomial-expansion Monte Carlo method to a spin-ice   Kondo lattice model
  We present the results of Monte Carlo simulation for a Kondo lattice model in which itinerant electrons interact with Ising spins with spin-ice type easy-axis anisotropy on a pyrochlore lattice. We demonstrate the efficiency of the truncated polynomial expansion algorithm, which enables a large scale simulation, in comparison with a conventional algorithm using the exact diagonalization. Computing the sublattice magnetization, we show the convergence of the data with increasing the number of polynomials and truncation distance. 
Hiroaki Ishizuka, Masafumi Udagawa, and Yukitoshi Motome
cond-mat.str-el