Computing microscopic structures of inhomogeneities in superconductors Original Research Article

We present an efficient computational method to solve selfconsistently the Bogoliubov–de Gennes equations of weak coupling superconductivity. As a function of system size, the scaling of the CPU time required by the scheme is shown to be preferable compared to the methods commonly used. Also, the scheme allows taking into account nonlocal pairing interactions without additional computational cost. The favourable scaling behavior enables computation of microscopic electronic structures for ranges of physical parameter values previously inaccessible.