A travelling cluster approximation for lattice fermions strongly coupled to classical degrees of freedom

We suggest and implement a new Monte Carlo strategy for correlated models involving fermions strongly coupled to classical degrees of freedom, with accurate handling of quenched disorder as well. Current methods iteratively diagonalise the full Hamiltonian for a system of N sites with computation time τ<Subscript>N</Subscript> ∼N<Superscript>4</Superscript>. This limits achievable sizes to N ∼100. In our method the energy cost of a Monte Carlo update is computed from the Hamiltonian of a cluster, of size N<Subscript>c</Subscript>, constructed around the reference site, and embedded in the larger system. As MC steps sweep over the system, the cluster Hamiltonian also moves, being reconstructed at each site where an update is attempted. In this method τ<Subscript>N,Nc</Subscript> ∼NN<Subscript>c</Subscript> <Superscript>3</Superscript>. Our results are obviously exact when N<Subscript>c</Subscript>=N, and converge quickly to this asymptote with increasing N<Subscript>c</Subscript>, particularly in the presence of disorder. We provide detailed benchmarks on the Holstein model and the double exchange model. The `locality' of the energy cost, as evidenced by our results, suggests that several important but inaccessible problems can now be handled with control. This method forms the basis of our studies in Europhys. Lett. <Emphasis Type="Bold">68, 564 (2004), Phys. Rev. Lett. <Emphasis Type="Bold">94, 136601 (2005), and Phys. Rev. Lett. <Emphasis Type="Bold">96, 016602 (2006). Copyright EDP Sciences/Società Italiana di Fisica/Springer-Verlag 2006