Fast Fourier-Chebyshev approach to real-space simulations of the Kubo formula

Santiago Giménez de Castro,João M. Viana Parente Lopes,Aires Ferreira,D. A. Bahamon
DOI: https://doi.org/10.1103/PhysRevLett.132.076302
2024-02-19
Abstract:The Kubo formula is a cornerstone in our understanding of near-equilibrium transport phenomena. While conceptually elegant, the application of Kubo's linear-response theory to interesting problems is hindered by the need for algorithms that are accurate and scalable to large lattice sizes beyond one spatial dimension. Here, we propose a general framework to numerically study large systems, which combines the spectral accuracy of Chebyshev expansions with the efficiency of divide-and-conquer methods. We use the hybrid algorithm to calculate the two-terminal conductance and the bulk conductivity tensor of 2D lattice models with over $10^7$ sites. By efficiently sampling the microscopic information contained in billions of Chebyshev moments, the algorithm is able to accurately resolve the linear-response properties of complex systems in the presence of quenched disorder. Our results lay the groundwork for future studies of transport phenomena in previously inaccessible regimes.
Mesoscale and Nanoscale Physics
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop an efficient and scalable numerical method for studying linear - response transport phenomena in two - dimensional lattice models. Specifically, the authors propose an algorithm, the Fast Fourier - Chebyshev method, which combines Chebyshev expansion and the Fast Fourier Transform (FFT), to calculate the conductivity tensor described by the Kubo formula. ### Main problems 1. **Computational complexity**: - Traditional methods require \(O(DM^2)\) number of operations when calculating double Chebyshev expansion moments, which is very time - consuming when dealing with large - scale systems. - The new method proposed in the paper reduces the computational complexity to \(O(DM \log M)\), making the calculation of large - scale systems possible. 2. **Accuracy and resolution**: - The new method can accurately analyze the linear - response characteristics of complex systems while maintaining high - energy resolution, especially in the presence of disorder. - Through efficient sampling, the algorithm can accurately capture microscopic information, thus better understanding the transport properties of the system. 3. **Scope of application**: - This method can be applied to various two - dimensional lattice models, including materials such as graphene and twisted bilayer graphene, to study their transport properties such as conductivity and Hall conductivity. - In particular, it provides a new tool for the study of quantum critical phenomena such as the quantum Hall effect and Anderson localization. ### Key contributions - **Algorithm innovation**: The introduction of the Fast Fourier Transform (FFT) and the divide - and - conquer strategy significantly improves computational efficiency. - **Wide applicability**: It is suitable for linear - response transport calculations of various two - dimensional materials and complex systems. - **High precision**: It can perform full - spectrum calculations with high - energy resolution and capture subtle physical features. In summary, by proposing the Fast Fourier - Chebyshev method, this paper solves the bottlenecks of traditional methods in terms of computational complexity and accuracy, and provides a new and effective tool for studying linear - response transport phenomena in complex systems.