# Magneto-orbital coupling in iron pnictides.

###### Abstract

A magneto-orbital coupling mechanism is proposed to account for the weak energy gap at the Fermi energy in the ordered SDW state of a realistic three band model for iron pnictides involving , , and Fe orbitals. The orbital mixing terms between the and orbitals, which are important in reproducing the orbital composition of the elliptical electron pockets at and , are shown to play a key role in the energy gap formation in the SDW state.

## 1 Introduction

The iron pnictides exhibit a rich phase diagram [1, 2] including magnetic, structural (tetragonal-to-orthorhombic), and superconducting phase transitions [3]. The magnetic state exhibits ordering of Fe moments in the - plane, with a concomitant structural distortion ¿ possibly correlated with the ferro orbital order as seen in ARPES studies [4]. Inelastic neutron scattering experiments [5, 6] have shown well defined spin wave excitations with energy scale 200 meV, persisting even above the Néel temperature [7]. Evidently, short range antiferromagnetic (AF) and ferromagnetic (F) order remain in the and directions respectively even above the disordering temperature for long-range magnetic order, which may account for the narrow nematic phase [8, 9] above the Néel temperature where the ferro orbital order [4] and structural distortion survive, as well as the temperature dependence of the measured anisotropies in and directions of magnetic excitations and resistivity [4, 7, 10].

Evidently, the complex multi-orbital character of the underlying microscopic description of iron pnictides plays an important role in understanding their macroscopic behavior as seen from the interplay of magnetic, orbital, structural, and transport properties. Another important macroscopic behavior of iron pnictides is associated with the energy gap in the ordered magnetic state, as seen in optical conductivity, transport and scanning tunneling microscopy studies; investigating the role of the multi-orbital character on the magnetic-state electronic structure is therefore of interest.

For the non-magnetic state, intensive investigations of the electronic structure of these materials by first-principle calculations [11, 12, 13, 14, 15] and angle resolved photoemission spectroscopy (ARPES) experiments [16, 17, 4, 18, 19] have revealed that the Fermi surface consists of two nearly circular hole pockets around the center and elliptical electron pockets around the corner of the BZ which are primarily contributed by , , and Fe orbitals.

A variety of experimental techniques have confirmed the opening of a small energy gap in the excitation spectrum associated with the onset of the spin-densityâwave (SDW) order in iron pnictides [20]. The Fermi surface (FS) undergoes a complex multi-orbital reconstruction through the non-magnetic to antiferromagnetic transition. ARPES experiments have revealed that the electronic states are strongly modified across the SDW transition [21, 22, 23]. Bands get folded and hybridized, resulting in band splitting and opening of SDW gaps. Consequently, the Fermi surface breaks up into small droplet-like structures, supporting only a fraction of the original FS areas. These results are consistent with Quantum Oscillation measurements [24, 25] showing the reconstruction and drastic reduction of Fermi surface area, which can be as small as of the original BZ area for .

In addition, investigations of optical and transport properties by optical spectroscopy [26, 27] and broadband spectroscopic ellipsometry [28] measurements on single-crystalline iron pnictides have found an energy gap for electronic excitations in the SDW ordered state. The carrier density is reduced substantially in the ordered state with consequent decrease in optical conductivity. Also, as inferred from the results of infrared studies, short-range AF correlations have been suggested as being responsible for the pseudogap in iron superconductors [29, 30]. Moreover, scanning tunneling microscopy (STM) studies have also seen an energy gap in the density of states near the Fermi energy [31, 32].

Close proximity of the superconducting phase and magnetically ordered parent phase in iron pnictides have naturally led to comparisons with cuprates [29]. However, the magnitude of the experimentally measured energy gap at the Fermi energy (up to meV for ) is significantly smaller than the usual SDW gap associated with AF ordering with staggered magnetization and local Coulomb interaction , as in the Hubbard model representation for cuprates. Moreover, these SDW gaps open at energies typically away from the Fermi energy in realistic multi-band tight binding models which yield the correct FS structure in agreement with DFT calculations and ARPES studies in the non-magnetic state.

In this paper, a magneto-orbital coupling mechanism is proposed for a realistic three-band model at half filling for iron pnictides involving , , and orbitals of Fe, which is shown to account for the weak non-conventional energy gap at the Fermi energy due to a composite effect of the orbital mixing (hybridization) between the and orbitals and the SDW magnetic order. For realistic three-band-model parameters, the energy gap is found to be in quantitative agreement with experiments. The gapped SDW state in the realistic three band model at half filling therefore provides, in analogy with cuprates, a suitable reference state for theoretically investigating effects of electron/hole doping in iron pnictides.

## 2 Three-orbital model and () SDW state

We consider a minimal three-orbital model [33] involving , and Fe orbitals. The tight binding Hamiltonian in the plane-wave basis is defined as:

(1) |

where

(2) |

are the tight-binding matrix elements in the unfolded BZ (). Here, and are the intra-orbital hoppings for () along () and () directions, respectively, and are the intra and inter-orbital hoppings along diagonal direction for and , and are intra-orbital NN and NNN hoppings for , while and the NN and NNN hybridization between and . Finally, is the energy difference between the and degenerate orbitals.

The Fermi surface for the three-orbital Hamiltonian (1) for the values of hopping parameters given in Table 1 is shown in Figure 1(a). Here, the Fermi energy is kept at eV corresponding to half filling. There are two circular hole pockets around the center and elliptical electron pockets around () and () in the unfolded BZ. The two hole pockets involve primarily the and orbitals, while the electron pockets centered at () [()] arise mainly from the hybridization of the and [] orbitals. All of these features are in good agreement with results from DFT calculations and ARPES experiments.

We now consider the () ordered magnetic (SDW) state of this model. The various electron-electron interaction terms included are:

(3) | |||||

where () refer to the local spin (charge) density operators for orbital . The first and second terms are the intra-orbital and inter-orbital Coulomb interactions respectively, the third term is the Hund’s coupling and the fourth term the “pair-hopping” term.

Extending the two-sublattice basis approach for the SDW state in a single-band model [34] to a composite three-orbital, two-sublattice basis, the Hartree-Fock (HF) level Hamiltonian matrix in this composite basis (A A A B B B) is obtained as:

for spin , where

are the band energies corresponding to the hopping terms along different directions, and the self-consistent exchange fields are defined as in terms of sublattice magnetization for orbital .

The calculated spin wave spectral function [33] in the SDW state for the three-orbital model is shown in Fig. 1(b). Evidently, spin wave excitations are highly dispersive, and do not decay into the particle-hole continuum. The energy scale of spin excitations is 200 meV with a well-defined maximum at the ferromagnetic zone boundary []. These features of spin wave excitations are in excellent agreement with results from inelastic neutron scattering measurements, confirming that our three-orbital model is realistic.

## 3 Electronic structure in the SDW state

Figure 2 shows the the evolution of the orbital resolved density of states (DOS) near the Fermi energy in the SDW state with increasing interaction strength . The corresponding band structures along symmetry directions of the BZ are shown in Fig. 3. Figures 2(a) and 3(a) correspond to the PM state for which the model yields correct FS structure as shown in Fig. 1(a).

As seen in Fig. 2, a large gap opens for the band in the SDW state. This band gap is a conventional SDW gap, which opens at wave vector in the AF direction, as seen from the band dispersion plot shown in Fig. 3. Similar SDW gaps open up for the and bands at . However, these band gaps appear at energies far away from the Fermi energy. Besides these three conventional SDW band gaps, there is a weak band gap involving the and orbitals which opens up near the Fermi energy, as seen in Fig. 3, which is found to be robust with respect to small variations in the hopping parameters.

It is important to note that while partial density of states for and orbitals are symmetric in the PM state [Fig. 2(a)], they become anisotropic in the SDW state [2(b)]. The DOS for the orbital has a peak further away from Fermi energy () than that for and has a relatively smaller contribution at . This is in excellent agreement with results of DMFT studies [35] and provides a good understanding of the ferro orbital ordering between the two orbitals [36, 37]. Moreover, the band dispersion, gap formation, as well as the largest gap for [Fig. 2(c) and 2(d)] agree well with the self-consistent calculations for multi-orbital models [38, 39].

Physically, this weak energy gap near the Fermi energy involving the and orbitals is the most important. One of the key contributing factors for this energy gap formation is the orbital mixing between the and orbitals due to the hopping terms and . Only will be retained in the discussion below as is negligible in comparison (see Table 1). Fig. 4 shows the SDW state band structure with and without , confirming the role of the orbital mixing term in the gap formation. It is important to note that the opening of this gap for the and orbitals due to the orbital mixing term critically requires the presence of the SDW magnetic ordering, as is evident from the absence of any gap in the non-magnetic state (), even with the mixing term included [Fig. 4 (a)]. For simplicity, we have taken equal exchange fields ( eV, eV) for all three orbitals.

To summarize, both the SDW order as well as the orbital mixing term between the and orbitals are critically important in the formation of the weak energy gap near the Fermi energy. Only the term (non-magnetic state) or only the SDW order () are independently not sufficient to open the gap, both must be present together.

## 4 Magneto-orbital coupling

A magneto-orbital coupling mechanism is proposed below to account for this non-conventional weak energy gap near the Fermi energy involving the and orbitals. In the () ordered SDW state with AF ordering in the direction, the orbital mixing term connects and orbitals on opposite sublattices, highlighting the coupling between the orbital and magnetic sectors.

The composite effect of the orbital mixing between the and orbitals and the () SDW order on opening the weak band gap near the Fermi energy can be clearly seen by going over to a reduced basis involving the relevant and bands and contracting over the sublattice basis to obtain the band mixing terms in the magnetic state.

We focus on the direction of the Brillouin zone which is the AF ordering direction, and set . As the orbital completely decouples in this case, we start with the part of the Hamiltonian matrix given in Eq. (5) involving the (A B A B) part of the orbital-sublattice basis for the and orbitals:

where is the mixing term between and orbitals on opposite sublattices, and the bare band energy terms are and .

In the absence of the above mixing term , the and orbital sectors are decoupled and in the () ordered SDW state, the SDW band energies for the upper (+) and lower (-) bands are given by:

(7) |

The exchange fields and open band gaps at where . However, these gaps open at energies away from the Fermi energy, as discussed earlier [see Fig. 4]. For close to zero, which is far away from where the band gaps opens, the band energies are nearly unchanged for small [see Fig. 4(c)], and are given by:

(8) |

Out of these four bands, we focus on the two relevant SDW bands (upper band) and (lower band) which lie near the Fermi energy. We now consider an effective Hamiltonian matrix in the reduced basis consisting of these two eigenstates and . The band mixing term in this reduced basis is obtained by evaluating the matrix elements using the corresponding eigenvectors in the two-sublattice basis:

(9) |

for the and band sectors of the Hamiltonian matrix in Eq. (8). This yields

where

(11) |

is a magneto-orbital coupling term involving the orbital mixing term and the SDW state matrix element of representing the opposite sublattice operator. Here and are the and sublattice components of the eigenvectors for the two orbitals:

(12) |

We note here that the majority/minority amplitudes on the two sublattices get exchanged for the upper band and lower band.

In the absence of any SDW order (), the magneto-orbital coupling term as the A and B sublattice amplitudes are identical in magnitude (), so the negative sign in Eq. (12) yields an exact cancellation. However, for finite and , the sublattice asymmetry in the magnetic state results in , which leads to strong mixing at -points where the two band energies and are nearly degenerate, resulting in the opening of the non-conventional weak energy gap near the Fermi energy. For and , we obtain:

(13) |

explicitly showing the composite dependence of the magneto-orbital coupling term on both the orbital mixing and SDW order, and thus accounting for the absence of the band gap near the Fermi energy if either of these two terms vanishes, as seen in Fig. 4.

## 5 Conclusions

A magneto-orbital coupling mechanism was proposed and shown to account for the opening of a non-conventional weak energy gap near the Fermi energy in the SDW state of a realistic three band model for iron pnictides at half filling. The coupling term involves a composite dependence on both the orbital mixing (hybridization) terms between the and orbitals as well as the SDW order parameter, thus accounting for the absence of the energy gap if either of these two terms is absent. As the orbital mixing terms also play an important role in the mixed composition ( and ) of the elliptical electron pockets in the non-magnetic-state FS structure, the magneto-orbital coupling provides a subtle link between the non-magnetic-state FS features and the SDW state weak energy gap. Furthermore, the non-conventional nature of the weak energy gap for quasiparticle excitations with orbital mixing provides a sharp contrast from the conventional strong SDW gap in cuprate antiferromagnets.

## Acknowledgements

SG and NR acknowledge financial support from Council of Scientific and Industrial Research, India.

## References

## References

- [1] J. Zhao, Q. Huang, C. de la Cruz, S. Li, J. W. Lynn, Y. Chen, M. A. Green, G. F. Chen, G. Li, Z. Li, J. L. Luo, N. L. Wang, and P. Dai, Nature Mat. 7, 953 (2008).
- [2] S. Nandi, M. G. Kim, A. Kreyssig, R. M. Fernandes, D. K. Pratt, A. Thaler, N. Ni, S. L. Bud’ko, P. C. Canfield, J. Schmalian, R. J. McQueeney, and A. I. Goldman, Phys. Rev. Lett. 104, 057006 (2010).
- [3] R. M. Fernandes, A. V. Chubukov, J. Schmalian, Nature Phys. 10, 97 (2014).
- [4] M. Yi, D. Lu, J.-H. Chu, J. G. Analytis, A. P. Sorini, A. F. Kemper, B. Moritz, S.-K. Mo, R. G.Moore, M. Hashimoto, W.-S. Lee, Z. Hussain , T. P. Devereaux, I. R. Fisher, and Z-X Shen, Proc. Nat. Acad. Science 108, 6878 (2011).
- [5] R. A. Ewings, T. G. Perring, J. Gillett, S. D. Das, S. E. Sebastian, A. E. Taylor, T. Guidi, and A. T. Boothroyd, Phys. Rev. B 83, 214519 (2011).
- [6] L. W. Harriger, H. Q. Luo, M. S. Liu, C. Frost, J. P. Hu, M. R. Norman, and P. Dai, Phys. Rev. B 84, 054544 (2011).
- [7] X. Lu, J. T. Park, R. Zhang, H. Luo, A. H. Nevidomskyy, Q. Si, and P. Dai, Science 345, 657 (2014).
- [8] S. Kasahara, H. J. Shi, K. Hashimoto, S. Tonegawa, Y. Mizukami, T. Shibauchi, K. Sugimoto, T. Fukuda, T. Terashima, A. H. Nevidomskyy, and Y. Matsuda, Nature 486, 382 (2012).
- [9] R. Zhou, Z. Li, J. Yang, D. L. Sun, C. T. Lin, and Guo-qing Zheng, Nature Comm.4, 2265 (2013).
- [10] J.-H. Chu, J. G. Analytis, K. D. Greve, P. L. McMahon, Z. Islam, Y. Yamamoto, and I. R. Fisher, Science 329, 824 (2010).
- [11] I. I. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, Phys. Rev. Lett. 101, 057003 (2008).
- [12] K. Haule, J. H. Shim, and G. Kotliar, Phys. Rev. Lett. 100, 226402 (2008).
- [13] D. J. Singh, and M.-H. Du, Phys. Rev. Lett. 100, 237003 (2008).
- [14] I. A. Nekrasov, Z. V. Pchelkina, and M. V. Sadovskii, JETP letters 88, 144 (2008).
- [15] H.-J. Zhang, G. Xu, X. Dai, and Z. Fang Z, Chin. Phys. Lett. 26, 017401 (2009).
- [16] M. Yi, D. H. Lu, J. G. Analytis, J.-H. Chu, S.-K. Mo, R.-H. He, M. Hashimoto, R. G. Moore, I. I. Mazin, D. J. Singh, Z. Hussain, I. R. Fisher, and Z.-X. Shen, Phys. Rev. B 80, 174510 (2009).
- [17] T. Kondo, R. M. Fernandes, R. Khasanov, Chang Liu, A. D. Palczewski, Ni Ni, M. Shi, A. Bostwick, E. Rotenberg, J. Schmalian, S. L. Bud’ko, P. C. Canfield, and A. Kaminski, Phys. Rev. B 81, 060507(R) (2010).
- [18] V. Brouet, M. F. Jensen, P.-H. Lin, A. Taleb-Ibrahimi, P. Le Fèvre, F. Bertran, C.-H. Lin, Wei Ku, A. Forget, and D. Colson, Phys. Rev. B 86, 075123 (2012).
- [19] A. A. Kordyuk, V. B. Zabolotnyy, D. V. Evtushinsky, A. N. Yaresko, B. Büchner, and S. V. Borisenko, J. Supercond. Nov. Magn. 26, 2837 (2013).
- [20] D. C. Johnston, Adv. Phys. 59, 803 (2010)
- [21] S. de Jong et al, Euro Phys. Lett. 89, 27007 (2010).
- [22] B. Zhou et al, Phys. Rev. B 81, 155124 (2010).
- [23] L. X. Yang et al, Phys. Rev. Lett. 102, 107002 (2009).
- [24] S. E. Sebastian, J. Gillett, N. Harrison, P. H. C. Lau, D. J. Singh, C. H. Mielke, and G. G. Lonzarich, J. Physics: Condens. Matter 20, 42203 (2008).
- [25] J. G. Analytis, R. D. McDonald, J.-H. Chu, S. C. Riggs, A. F. Bangura, C. Kucharczyk, M. Johannes, and I. R. Fisher, Phys. Rev. B 80, 064507 (2009).
- [26] W. Z. Hu, J. Dong, G. Li, Z. Li, P. Zheng, G. F. Chen, J. L. Luo, and N. L. Wang, Phys. Rev. Lett. 101, 257005 (2008).
- [27] D. Wu et al, Phys. Rev. B 79, 155103 (2009).
- [28] A. Charnukha, D. Pröpper, T. I. Larkin, D. L. Sun, Z. W. Li, C. T. Lin, T. Wolf, B. Keimer, and A. V. Boris, Phys. Rev. B 88, 184511 (2013).
- [29] S. J. Moon et al, Phys. Rev. B 90, 014503 (2014)
- [30] S. J. Moon et al, Phys. Rev. Lett. 109, 027006 (2012).
- [31] H. Zhang et al, Phys. Rev. B 81, 104520 (2010).
- [32] A. Dutta, Anupam, Z. Hossain, and A. K. Gupta, J. Physics: Condens. Matter 25, 375602 (2013).
- [33] S. Ghosh and A. Singh, New J. Phys. 17, 063009 (2015).
- [34] N. Raghuvanshi, and A. Singh, J. Physics: Condens. Matter 22, 422202 (2010).
- [35] Z. P. Yin, K. Haule, and G. Kotliar, Nature Phys. 7, 294 (2011).
- [36] Y. K. Kim et al, Phys. Rev. Lett. 111, 217001 (2013).
- [37] S. Ghosh and A. Singh, J. Appl. Phys. 115, 103907 (2014).
- [38] N. Plonka, A. F. Kemper, S. Graser, A. P. Kampf, and T. P. Devereaux, Phys. Rev. B 88, 174518 (2013).
- [39] M. Yi et al, Nature Comm. 5, 3711 (2014).