Efficient Numerical Self-Consistent Mean-Field Approach for Fermionic Many-Body Systems by Polynomial Expansion on Spectral Density

Abstract

We propose an efficient numerical algorithm to solve Bogoliubov de Gennes\nequations self-consistently for inhomogeneous superconducting systems with a\nreformulated polynomial expansion scheme. This proposed method is applied to\ntypical issues such as a vortex under randomly distributed impurities and a\nnormal conducting junction sandwiched between superconductors. With various\ntechnical remarks, we show that its efficiency becomes remarkable in\nlarge-scale parallel performance.\n