A pronounced local in-gap zero-energy bound state (ZBS) has been observed by recent scanning tunneling microscopy (STM) experiments on the interstitial Fe impurity (IFI) and its nearest-neighboring (nn) sites in $mathrm{FeTe_{0.5}Se_{0.5}}$ supercond
ucting (SC) compound. By introducing a new impurity mechanism, the so-called tunneling impurity, and based on the Bogoliubove-de Gennes (BDG) equations we investigated the low-lying energy states of the IFI and the underlying Fe-plane. We found the peak of ZBS does not shift or split in a magnetic field as long as the tunneling parameter between IFI and the Fe-plane is sufficiently small and the Fe-plane is deep in the SC state. Our results are in good agreement with the experiments. We also predicted that modulation of spin density wave (SDW), or charge density wave (CDW) will suppress the intensity of the ZBS.
Effects of disorder on electron-doped iron pnictides are investigated systematically based on self-consistent Bogoliubov-de Gennes equations. Multiply impurities with same scattering potential (SP) are randomly distributed in a square lattice. Probab
ility distribution functions of normalized order parameters for different impurity concentrations $delta_{imp}$, different electron doping concentrations $delta$ are investigated for given SPs. Samples are found to be very robust against weak SP, in which order parameters do not have qualitative change even at very large $delta_{imp}$. While strong SP is able to easily break down the order parameters. For moderate SP, variations of order parameters on and around impurities strongly depend on $delta$, however the distribution functions of normalized order parameters have similar behavior as $delta_{imp}$ increases. Compared with superconducting (SC) order, the magnetic order is more sensitive to multi-impurity effect. The spatial spin density wave pattern has already been destroyed before the system loses its superconductivity. Dependence of SC order on temperature is similar to that of impurity-free case, with the critical temperature being remarkably suppressed for high $delta_{imp}$.
By making use of renormalized mean-field theory, we investigate possible superconducting symmetries in the ground states of t1-t2-J1-J2 model on square lattice. The superconducting symmetries of the ground states are determined by the frustration amp
litude t2/t1 and doping concentration. The phase diagram of this system in frustration-doping plane is given. The order of the phase transitions among these different superconducting symmetry states of the system is discussed.
