2003/08/14 by Nobuo Furukawa, Yukitoshi Motome
Physics and Astronomy · #cond-mat.str-el #cond-mat.stat-mech
paper · pdf · doi:10.1143/jpsj.73.1482
published as J. Phys. Soc. Jpn. 73, 1482 (2004) · 9 pages with 4 figs
arxiv created 2003/08/14 · arxiv updated 2009/12/01
An improved algorithm is proposed for Monte Carlo methods to study fermion systems interacting with adiabatical fields. To obtain a weight for each Monte Carlo sample with a fixed configuration of adiabatical fields, a series expansion using Chebyshev polynomials is applied. By introducing truncations of matrix operations in a systematic and controlled way, it is shown that the cpu time is reduced from O(N3) to O(N) where N is the system size. Benchmark results show that the implementation of the algorithm makes it possible to perform systematic investigations of critical phenomena using system-size scalings even for an electronic model in three dimensions, within a realistic cpu timescale.