## In-gap Bound States Induced by Interstitial Fe Impurities in Iron-based Superconductors

Degang Zhang 1, 2

1 College of Physics and Electronic Engineering, Sichuan Normal University, Chengdu 610101, China 2 Institute of Solid State Physics, Sichuan Normal University, Chengdu 610101, China

Based on a two-orbit four-band tight binding model, we investigate the low-lying electronic states around the interstitial excess Fe ions in the iron-based superconductors by using T-matrix approach. It is shown that the local density of states at the interstitial Fe impurity (IFI) possesses a strong resonance inside the gap, which seems to be insensitive to the doping and the pairing symmetry in the Fe-Fe plane, while a single or two resonances appear at the nearest neighboring (NN) Fe sites. The location and height of the resonance peaks only depend on the hopping t and the pairing parameter ∆ I between the IFI and the NN Fe sites. These in-gap resonances are originated in the Andreev's bound states due to the quasiparticle tunneling through the IFI, leading to the change of the magnitude of the superconducting order parameter. When both t and ∆ I are small, this robust zero-energy bound state near the IFI is consistent with recent scanning tunneling microscopy observations.

PACS numbers: 71.10.Fd, 71.18.+y, 71.20.-b, 74.20.-z

The discovery of different families of the iron-based superconductors provided us a new platform to study high temperature superconductivity [1-6]. It has been reported that the superconducting transition temperature T c can reach as high as 55K [2]. Similar to the cuprate superconductors, the iron-based superconductors also have a layer structure and the superconductivity comes from the Cooper pairs in the Fe-Fe plane by doping electrons or holes. However, different from the oxygen ions at the Cu-Cu bonds in the cuprates, two As/Se/Te ions in each unit cell locate just above (below) the center of one face of the Fe square lattice and below (above) the center of the neighboring face, respectively (see Fig. 1). Such a combination of the Fe ions and the As/Se/Te ions produces a more complex energy band structure. Angle resolved photoemission spectroscopy (ARPES) experiments have revealed that the iron-based superconductors usually possess two hole Fermi surfaces around (0,0), i.e. Γ point, and two electron Fermi surfaces around ( π, π ), i.e. M point [7-16]. We note that due to different dopings, one or even two hole Fermi surfaces disappear in some families of the iron-based superconductors. In order to explain these Fermi surface characteristics, twoorbital, three-orbital, and five-orbital tight binding models have been proposed [17-20].

It is known that the impurity effects play a crucial role in exploring superconductivity [21]. The strong zeroenergy bound state (ZBS) at the Zn ion and the fourfold symmetric electronic states around the single impurity, observed by scanning tunneling microscopy (STM) [22], demonstrate that the superconducting order parameter in the cuprates has a d-wave symmetry. The in-gap resonances induced by a nonmagnetic impurity was also used to confirm the s + -pairing symmetry in the iron-based superconductors [18]. It is clear that these resonance peaks are originated in the Andreev's bound states due

FIG. 1: (Color online) Schematic lattice structure of the FeAs/Se/Te layers with each unit cell containing two Fe (A and B) and two As/Se/Te (A and B) ions. The As/Se/Te ions A and B are located just below and above the center of each face of the Fe square lattice, respectively. Here, t 1 is the nearest neighboring (NN) hopping between the same orbitals d xz or d yz . t 2 and t 3 are the next nearest neighboring hoppings between the same orbitals mediated by the As/Se/Te ions B and A, respectively. t 4 is the next nearest neighboring hopping between the different orbitals. The IFI is situated at the symmetric position of As/Se/Te ion about the Fe-Fe plane. t and ∆ I are the hopping and the pairing parameters between the same orbitals of the IFI and the NN Fe sites, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0567v1-figure-e661f9a7-6e2e-43dd-88d6-5bb47b6b4835.jpg)

to the quasiparticle scattering between the Fermi surfaces or parts of the Fermi surface with opposite phases of the superconducting order parameter.

Very recently, Yin et al. also used STM to study the electronic states near the interstitial excess Fe ions in the iron-based superconductor Fe(Te,Se) [23]. The as-grown Fe(Te,Se) single crystals usually contain a large amount of excess Fe ions situated randomly at the interstitial

sites, which are the symmetric points of the Te/Se ions about the Fe-Fe plane (see Fig. 1). They found that the differential conductance at the interstitial Fe impurity (IFI) also has a strong ZBS. The height of this resonance decays rapidly with the distance from the IFI and vanishes at about 10 ˚ . A It is especially surprising that the ZBS at the IFI is not affected by a magnetic field up to 8 Tesla. Such a ZBS induced by out of the superconducting plane impurities was never reported before, and its origin is unclear up to now.

In this work, we provide an explanation for the interesting observation of the robust ZBS. Because the big family of the iron-based superconductors has a similar energy band structure, here we employ the two-orbit fourband tight binding model presented in Ref. [18] to explore the origin of the ZBS near the IFI. We remember that this energy band structure was constructed by starting two degenerate orbitals d xz and d yz of per Fe ion and two Fe ions in each unit cell, and fitted well the ARPES observations. Because the asymmetry of the above and below As/Se/Te ions in the Fe-As/Se/Te surface layer was taken into account, the model naturally explained several important STM observations, e.g. in-gap impurity resonances [18,24], the bound state at negative energy in the vortex core [25,26], the domain walls [27-30], etc., and especially repeated the phase diagram measured by nuclear magnetic resonance and neutron scattering experiments [31-33].

It is obvious that the IFI is different from the nonmagnetic, magnetic, and Kondo impurities in the superconducting plane. Because Fe ion is assumed to possess two electron channels, the electrons in the Fe-Fe plane can tunnel through the IFI. Now we know that superconductivity in the iron-based superconductors originates from electron pairing between the next nearest neighboring Fe sites. So only one electron of a Cooper pair can arrive at the IFI in the superconducting state. Therefore, the Hamiltonian describing the IFI in the iron-based superconductors has the form

$$H = H _ { 0 } + H _ { \text{BCS} }$$

$$- t \sum _ { \alpha \sigma } \{ c _ { \alpha \sigma } ^ { + } [ c _ { A \alpha, 0 0 \sigma } + c _ { B \alpha, 0 0 \sigma } + c _ { A \alpha, 0 - 1 \sigma } + c _ { B \alpha, 1 0 \sigma } ] + \text{h.c.} \} \quad \text{$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdots$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cd.$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdcol$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdor$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdod$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdock$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdok$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdo$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cd<$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdt$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdet$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdore$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdol$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdon$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdode$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cd0$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdlet$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdos$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$\cdot$$$

$$+ \Delta _ { I } \sum _ { \alpha } \{ c _ { \alpha \uparrow } ^ { + } [ c _ { A \alpha, 0 0 \downarrow } ^ { + } + c _ { A \alpha, 0 - 1 \downarrow } ^ { + } + c _ { B \alpha, 0 0 \downarrow } ^ { + } + c _ { B \alpha, 1 0 \downarrow } ^ { + } ] \quad \text{the} \quad \text{ing}$$

$$+ c _ { \alpha, \downarrow } ^ { + } [ c _ { A \alpha, 0 0 \uparrow } ^ { + } + c _ { A \alpha, 0 - 1 \uparrow } ^ { + } + c _ { B \alpha, 0 0 \uparrow } ^ { + } + c _ { B \alpha, 1 0 \uparrow } ^ { + } ] + h. c. \}, \ ( 1 ) \quad \text{ferent} \\ \text{energy}.$$

where H 0 is the two-orbit four-band tight binding model proposed in Ref. [18], H BCS is the mean field BCS pairing Hamiltonian in the Fe-Fe plane, c + A B α,ijσ ( ) ( c A B α,ijσ ( ) ) creates (destroys) an α electron with spin σ (= ↑ or ↓ )

in the unit cell { i, j } of the sublattice A (B), α = 0(1) represents the degenerate orbital d xz ( d yz ), and t and ∆ are the hopping and pairing parameters between the I same orbitals of the IFI and the nearest neighboring (NN) Fe sites, respectively.

We note that the sub-Hamiltonian H 0 + H BCS can be diagonalized by suitable transformations. Introducing first the Fourier transformation c A B α,ijσ ( ) = 1 √ N ∑ k c A B α, ( ) k σ e i ( k x x i + k y y j ) with N the number of unit cells and the canonical transformations c Aα, k σ = ∑ uv ( -1) αv a u, k Γ u, k φ uv, k σ and c Bα, k σ = ∑ uv ( -1) αv /epsilon1 ∗ T, k Γ u, k φ uv, k σ with u, v = 0 and 1, a u, k = 1 2 ( /epsilon1 A, k -/epsilon1 B, k ) + ( -1) u √ 1 4 ( /epsilon1 A, k -/epsilon1 B, k ) 2 + /epsilon1 T, k /epsilon1 ∗ T, k , /epsilon1 A, k = -2( t 2 cos k x + t 3 cos k y ) , /epsilon1 B, k = -2( t 2 cos k y + t 3 cos k x ), /epsilon1 T, k = -t 1 [1 + e ik x + e ik y + e i ( k x + k y ) ], and Γ u, k = √ 2( a 2 u, k + /epsilon1 T, k /epsilon1 ∗ T, k ), and then taking the Bogoliubov transformations φ uv k ↑ = ∑ ν =0 1 , ( -1) ν ξ uv k ν ψ uv k ν and φ + uv - ↓ k = ∑ ν =0 1 , ξ uv k 1 -ν ψ uv k ν with constants ξ uv k ν to be determined below, the total Hamiltonian H finally becomes

$$H = \sum _ { u v { \mathbf k } \nu } ( - 1 ) ^ { \nu } \Omega _ { u v { \mathbf k } } \psi _ { u v { \mathbf k } \nu } ^ { + } \psi _ { u v { \mathbf k } \nu }$$

$$\ n g \exp _ { N } + \frac { 1 } { \sqrt { N } } \sum _ { u v \alpha k \nu } ( \mathcal { M } _ { u v k \nu, \alpha } c ^ { + } _ { \alpha \uparrow } \psi _ { u v k \nu } + \mathcal { N } _ { u v k \nu, \alpha } c ^ { + } _ { \alpha \downarrow } \psi ^ { + } _ { u v k \nu } + \text{h.c.} ),$$

(2) where Ω uv k = √ ( E uv k -µ ) 2 +∆ 2 uv k , E uv k = 1 2 ( /epsilon1 A, k + /epsilon1 B, k ) -2( -1) v t 4 (cos k x + cos k y ) + ( -1) u √ 1 4 ( /epsilon1 A, k -/epsilon1 B, k ) 2 + /epsilon1 T, k /epsilon1 ∗ T, k , M uv k ν,α = ( -1) αv τ ∗ u, k [ - -( 1) ν tξ uv k ν + ∆ I ξ uv k 1 -ν ], N uv k ν,α = ( -1) αv τ u, k [ -tξ uv k 1 -ν + ( -1) ν ∆ I ξ uv k ν ], τ u, k = [(1 + e ik y ) a u, k + (1 + e -ik x ) /epsilon1 T, k ] / Γ u, k , ξ uv k 0 ξ uv k 1 = ∆ uv k / (2Ω uv k ), and ξ 2 uv k ν = 1 2 [1+( -1) ν ( E uv k -µ / ) Ω uv k ]. Here, µ is the chemical potential to be determined by doping and ∆ uv k is the superconducting order parameter in the Fe-Fe plane. We can see from the Hamiltonian (2) that the hopping t may form new electron-fermion pairs while the pairing parameter ∆ I can also induce a hopping term. In other words, new quasiparticles are excited if electrons tunnel through the IFI. It is expected that such a novel mechanism produces the ZBS observed by the STM [23].

The Hamiltonian H in Eq. (2) has a quadratic form of the fermionic operators and can be exactly solved by using the T-matrix approach [18,21]. The analytic expression of the LDOS near the IFI has been obtained at different bias voltages. In our calculations, we have used the energy band parameters: t 1 = 0 5 . , t 2 = 0 2 . , t 3 = -1 0, . and t 4 = 0 02. . We note that the iron-based superconductor Fe(Te,Se) in the STM experiment [23] has a small gap of about 1.5 meV. In order to clearly observe the change of the impurity resonance peaks with the bias

FIG. 2: (Color online) The LDOS at the IFI as a function of the bias voltage eV under different t and ∆ I for the s + -pairing symmetry ∆ uv k = 1 2 ∆ (cos 0 k x +cos k y ) with ∆ 0 = 5 8 . meV at optimal electron doping ( µ = -0 54), where there are . two Fermi surface sheets around Γ point.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0567v1-figure-779ca478-99f8-46e8-9ad4-0b66f60dfeb9.jpg)

voltage, here we employ a large gap, i.e. 5.8 meV, in the electron-doped BaFe 2 -x Co As x 2 [34,35], where the excess Fe or Co ion is also regarded as an IFI. In Fig. 2, we plot the LDOS at the IFI as a function of the bias voltage eV under different t and ∆ I for the s + -pairing symmetry at the chemical potential µ = -0 54, where there are two . Fermi surface sheets around the center of the Brillouin zone, i.e. Γ point. Obviously, when t or ∆ I is small, i.e. &lt; ∼ 3 meV, the LDOS has a strong resonance peak at zero energy, and the superconducting coherence peaks do not show up. With increasing t or ∆ , the sharp zero I energy resonance peak moves toward positive or negative energy, respectively, and the superconducting coherence peaks appear at the symmetric bias voltages in Fig. 2(a) and (b). The height of the coherence peak at the positive energy is higher (lower) than that at the negative energy for the parameter t (∆ ). We note that when I t &lt; ∼ 4 meV and ∆ I = 0 meV, where the coherence peak at the negative energy disappears, the LDOS at the IFI has a singleor double-peak structure at the positive energies, very similar to the STM observations in the black regions of the superconducting BaFe 2 -x Co As x 2 sample [35]. When t = ∆ I =3 and 5 meV, the LDOS possesses the strong ZBS and the symmetric coherence peaks, shown in Fig. 2(c). However, the superconducting energy gap at the IFI is smaller that in the bulk, i.e. ∆ = 5 8 meV, and be0 . comes larger if t or ∆ I increases. When t is large enough to neglect ∆ I , the superconducting energy gap at the IFI is equal to ∆ 0 , and there is no in-gap bound state in Fig. 2(d).

Fig. 3 shows the LDOS at the NN Fe sites of the IFI for the s + -pairing symmetry under µ = -0 54. . When t (∆ ) I ∼ 2 meV and ∆ ( ) = 0 meV, there is a ZBS I t appeared in the LDOS (see Fig. 3(a) and (b)). With

FIG. 3: (Color online) The LDOS at the nearest neighboring Fe sites as a function of the bias voltage eV under different t and ∆ I for the s + -pairing symmetry ∆ uv k = 1 2 ∆ (cos 0 k x + cos k y ) with ∆ 0 = 5 8 meV at optimal electron doping ( . µ = -0 54), . where there are two Fermi surface sheets around Γ point.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0567v1-figure-539235f9-a278-4fa8-9fde-f8debf2fa300.jpg)

increasing t or ∆ , the single resonance splits into two I peaks at symmetric energies. The resonance peak at positive energy has higher height than that at negative energy. It is interesting that the curves of the LDOS in Fig. 3(a) are same with those in Fig. 3(b) if t and ∆ I are replaced each other. In other words, t and ∆ I play the same role for the LDOS at the NN Fe sites of the IFI. This is the result that the Cooper pairs in the Fe-Fe plane tunnel through the IFI. We note that when t = ∆ , there is I always a ZBS. If t = ∆ I is large, a strong resonance peak shows up at positive bias voltages in Fig. 3(c). However, for t &gt;&gt; ∆ , the in-gap resonances move outside of the I superconducting energy gap (see Fig. 3(d)).

In order to understand the variation of the in-gap bound states with the doping, we plot the LDOS curves near the IFI for the s + -pairing symmetry under µ =-0.3 and -0.45, where there is no Fermi surface sheet and one Fermi surface sheet around Γ point, respectively, in Fig. 4. Obviously, in both cases, the LDOS at the IFI in Fig. 4(a) and (b) has a sharp zero energy resonance, and the superconducting coherence peaks are completely suppressed when t = ∆ . Meanwhile, there is only a ZBS I at the NN Fe sites of the IFI in Fig. 4(c) and (d). The larger the parameter t = ∆ , the higher the zero energy I resonance peak.

Now we investigate the low-lying electronic states induced by the IFI under different pairing symmetry in the bulk. In Fig. 5, we present the LDOS near the IFI at different physical parameters for the s ++ pairing symmetry. After comparing Fig. 5 with Fig. 2(c), Fig. 3(c), and Fig. 4, we find that the LODS curves for the s + -pairing symmetry are identical with those for the s ++ pairing symmetry in the presence of the same physical

FIG. 4: (Color online) The LDOS at the IFI in (a) and (b) and the nearest neighboring Fe sites in (c) and (d) as a function of the bias voltage eV under µ = -0 3, where there is no Fermi . surface sheet around Γ point, and µ = -0 45, where there is . one Fermi surface sheet around Γ point, for the s + -pairing symmetry ∆ uv k = 1 2 ∆ (cos 0 k x +cos k y ) with ∆ 0 = 5 8 meV. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0567v1-figure-b5277ee9-e27b-4dfe-9b3c-b5b2f2a42b3a.jpg)

parameters. This means that the in-gap resonance peaks induced by the IFI is irrelative to the symmetry of the superconducting order parameter. Therefore, this out of the superconducting plane IFI is apparently different from the Zn impurity in cuprates [21,22] and the nonmagnetic impurity in the iron-based superconductors [18,24].

In summary, we have explored the impact of the IFI on the electronic states in the iron-based superconductors. This IFI is dealt as two electron channels with the orbital index α , which can also break the Cooper pairs when electrons tunnel through it. It is undoubted that re-forming new Cooper pairs with different magnitudes of the superconducting order parameter produces the in-gap resonances, which are determined by t and ∆ I and are independent of the pairing symmetry in the bulk. This implies that the resonance peaks induced by the IFI are not affected by a magnetic field. Obviously, the origin of these in-gap bound states is due to the Andreev reflection with the change of the magnitude rather than the phase of the superconducting order parameter, which is never reported in the superconductivity literature. We note that when both t and ∆ I are small, the LDOS near the IFI always has a single resonance peak at zero bias voltage, which is also irrelative to the doping. Such a robust ZBS is consistent with recent STM observations in the iron-based superconductor Fe(Te,Se) [23].

The author would like to thank Jiaxin Yin and Ang Li for useful discussions. This work was supported by the Sichuan Normal University.

FIG. 5: (Color online) The LDOS at the IFI in (a)-(c) and the nearest neighboring Fe sites in (d)-(f) as a function of the bias voltage eV under different t and ∆ I for the s ++ pairing symmetry ∆ uv k = 1 2 ∆ (cos 0 | k x +cos k y ) | with ∆ 0 = 5 8 meV . at different electron dopings.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0567v1-figure-7b51aa40-a31c-4db1-9496-01a2d1348de3.jpg)

- [1] Y. Kamihara et al. , J. Am. Chem. Soc. 130 , 3296 (2008).
- [2] Z. A. Ren et al. , Chin. Phys. Lett. 25 , 2215 (2008).
- [3] X. H. Chen et al. , Nature (London) 453 , 761 (2008).
- [4] C. de la Cruz et al. , Nature (London) 453 , 899 (2008).
- [5] G. F. Chen et al. , Phys. Rev. Lett. 100 , 247002 (2008).
- [6] H.-H Wen et al. , Europhys. Lett. 82 , 17009 (2008).
- [7] H. Ding et al. , Europhys. Lett. 83 , 47001 (2008).
- [8] D. H. Lu et al. , Nature (London) 455 , 81 (2008).
- [9] C. Liu et al. , Phys. Rev. Lett. 101 , 177005 (2008).
- [10] T. Kondo et al. , Phys. Rev. Lett. 101 , 147003 (2008).
- [11] D. V. Evtushinsky et al. , Phys. Rev. B 79 , 054517 (2009).
- [12] K. Nakayama et al. , Europhys. Lett. 85 , 67002 (2009) .
- [13] V. Zabolotnyy et al. , Nature (London) 457 , 569 (2009).
- [14] K. Terashima et al. , PNAS 106 , 7330 (2009) .
- [15] Y. Sekiba et al. , New J. Phys. 11 , 025020 (2009).
- [16] H. Miao et al. , Phys. Rev. B 85 , 094506 (2012).
- [17] S. Raghu et al. , Phys. Rev. B 77 , 220503 (2008).
- [18] Degang Zhang, Phys. Rev. Lett. 103 , 186402 (2009); ibid. 104, 089702 (2010).
- [19] P. A. Lee and X.-G. Wen, Phys. Rev. B 78 , 144517 (2008).
- [20] K. Kuroki et al. , Phys. Rev. Lett. 101 , 087004 (2008).
- [21] A. V. Balatsky, I. Vekhter, and J.-X. Zhu, Rev. Mod. Phys. 78 , 373 (2006).
- [22] S.H. Pan et al. , Nature (London) 403 , 746 (2000).
- [23] J.-X. Yin et al. , arXiv:1403.1027.
- [24] S. Grothe et al. , Phys. Rev. B 86 , 174503 (2012).
- [25] L. Shan et al. , Nature Phys. 7 , 325 (2011).

- [26] Yi Gao et al. , Phys. Rev. Lett. 106 , 027004 (2011).
- [27] T.-M Chuang et al. , Science 327 , 181 (2010).
- [28] Guorong Li et al. , arXiv:1006.5907.
- [29] Huaixiang Huang et al. , Phys. Rev. B 83 , 134517 (2011).
- [30] Bo Li et al. , New J. Phys. 15 , 103018 (2013).
- [31] Y. Laplace et al. , Phys. Rev. B 80 , 140501 (2009).
- [32] M.-H. Julien et al. , Europhys. Lett. 87 , 37001 (2009).
- [33] Tao Zhou, Degang Zhang, and C. S. Ting, Phys. Rev. B 81 , 052506 (2010).
- [34] Y. Yin et al. , Phys. Rev. Lett. 102 , 097002 (2009).
- [35] Ang Li et al. , (to be published).
Degang Zhang
In-gap Bound States Induced by Interstitial Fe Impurities in Iron-based Superconductors
Based on a two-orbit four-band tight binding model, we investigate the
low-lying electronic states around the interstitial excess Fe ions in the
iron-based superconductors by using T-matrix approach. It is shown that the
local density of states at the interstitial Fe impurity (IFI) possesses a
strong resonance inside the gap, which seems to be insensitive to the doping
and the pairing symmetry in the Fe-Fe plane, while a single or two resonances
appear at the nearest neighboring (NN) Fe sites. The location and height of the
resonance peaks only depend on the hopping $t$ and the pairing parameter
$\Delta_I$ between the IFI and the NN Fe sites. These in-gap resonances are
originated in the Andreev's bound states due to the quasiparticle tunneling
through the IFI, leading to the change of the magnitude of the superconducting
order parameter. When both $t$ and $\Delta_I$ are small, this robust
zero-energy bound state near the IFI is consistent with recent scanning
tunneling microscopy observations.