Probing Band-center Anomaly with the Kernel Polynomial Method

We investigate the anomalous behavior of localization length of a non-interacting one-dimensional Anderson model at zero temperature. We report numerical calculations of the Thouless expression of localization length, based on the Kernel polynomial method (KPM), which has an O(N ) computational complexity, where N is the system size. The KPM results show an excellent agreement with perturbative result in large system size limit, confirming the validity of Thouless formula. Thus, contrary to the previous numerical results, the KPM approximation of the Thouless expression produce the correct localization length at the band center. The Thouless expression relates localization length in terms of density of states in a one-dimensional disordered system. By calculating the KPM estimates of density of states, we find a cusp-like behavior around the band center in the perturbative regime. This cusp-like singularity can not be obtained by an approximate analytical calculations with in the second-order approximations, reflects the band-center anomaly.