Visualization of atomic-scale phenomena in superconductors: application   to FeSe
  We propose a simple method of calculating inhomogeneous, atomic-scale phenomena in superconductors which makes use of the wave function information traditionally discarded in the construction of tight-binding models used in the Bogoliubov-de Gennes equations. The method uses symmetry based first principles Wannier functions to visualize the effects of superconducting pairing on the distribution of electronic states over atoms within a crystal unit cell. Local symmetries lower than the global lattice symmetry can thus be exhibited as well, rendering theoretical comparisons with scanning tunneling spectroscopy data much more useful. As a simple example, we discuss the geometric dimer states observed near defects in superconducting FeSe. 
Peayush Choubey, T. Berlijn, A. Kreisel, C. Cao and P. J. Hirschfeld
cond-mat.supr-con