Migdal-Eliashberg calculation considering the retardation effect in phonon-mediated superconductors

We formulate an efficient scheme to perform Migdal-Eliashberg calculation considering the retardation effect from first principles^[1]. While the conventional approach requires a huge number of Matsubara frequencies, we show that the intermediate representation (IR) of the Green's function^[2] and the sparse sampling method^[3] dramatically reduce the numerical cost to solve the linearized gap equation. Without introducing any empirical parameter, we demonstrate that we can reproduce the experimental superconducting transition temperature of elemental Nb (~ 10 K). The present result indicates that our approach has a superior performance for many superconductors for which T_c is lower than the order of 10 K. Our scheme based on IR reduces the computation time cost and memory footprint by 20 times, compared to the conventional method using 4000 Matsubara frequencies.

[1] T. Wang, T. Nomoto, Y. Nomura, H. Shinaoka, J. Otsuki, T. Koretsune and R. Arita, Phys. Rev. B 102, 134503 (2020)
[2] H. Shinaoka, J. Otsuki, M. Ohzeki and K. Yoshimi, Phys. Rev. B 96, 035147 (2017)
[3] J. Li, M. Wallerberger, N. Chikano, C.-N. Yeh, E. Gull, and H. Shinaoka, Phys. Rev. B 101, 035144 (2020)