Electron-phonon interaction in Fe-based superconductors: Coupling of magnetic moments with phonons in LaFeAsOF
The coupling of Fe magnetic moments in LaFeAsOF with the As phonon is calculated. We present first principles calculations of the atomic and electronic structure of LaFeAsO as a function of electron doping. We perform calculations using the virtual crystal approximation as well as supercell calculations with F substitutional impurity atoms. The results validate the virtual crystal approximation for the electronic structure near the Fermi level. Its is found that the electronic density of states at the Fermi level is maximum for x=0.125, enhancing the electron-phonon interaction. An additional increase of the electron-phonon parameter is obtained if the coupling between the phonon and the Fe magnetic moment is included. It is found that the electron-phonon interaction can be one order of magnitude larger than its value if no spin resolution is included in the calculation. The implications of these results on the superconducting transition are discussed
pacs:71.20.-b, 71.38.-k, 74.20.-z, 74.70.-b
Since the discovery of the Fe-based superconductors a few years ago Kamihara et al. (2006, 2008), several properties common to most of the compounds have been well established and general agreement has been reached. However, there are very important aspects that have not been settled Norman (2008). The most important one, is the mechanism responsible for their superconducting behavior. Early theoretical work Mazin et al. (2008); Boeri et al. (2008) ruled out any phonon’s role. More recent calculations Boeri et al. (2010) reveal that, in the non-magnetic phase, phonons coupled with magnetism contribute to superconductivity although not enough to explain the high critical temperature. Another aspect to be clarified involves the Fe magnetic moment. The low temperature stripped antiferromagnetic order in the parent compound is generally accepted, although the magnitude of the calculated magnetic moment is much larger than the observed one Korotin et al. (2008). Fluctuations mechanisms as well as orbital magnetic ordering Bascones et al. (2010); Lee et al. (2010) have been proposed to explain this discrepancy. The coupling of the Fe magnetic moments with phonons has been clearly established at least in CaFeAs, in the undoped orthorombic phase Hahn et al. (2009), and in SmFeAsOF Tacon et al. (2009) . Also the possibility of a non-zero Fe magnetic moment in the superconducting phase has not been ruled out Reznik et al. (2009).
Here we propose that the magnetic moment in Fe atoms does not disappear in the superconducting phase. Instead, what disappears is the long-range collinear antiferromagnetic order while keeping the As-mediated antiferromagnetic local coupling between the Fe atoms. To this end, we have calculated in the manner described in Yndurain and Soler (2009) how the energy barrier between two equivalent collinear antiferromagnetic configurations varies with the extra electron concentration in LaFeAsOF. We assume that a non-collinear configuration takes place between the two minima as discussed in Yndurain and Soler (2009). The results of the calculation are shown in Figure 1. The results indicate that the barrier between two equivalent antiferromagnetic arrangement is very small beyond x=0.1 up to the disappearance of the magnetic moment around x=0.3. This low barrier can be responsible for the disappearance of the antiferromagnetic long-range order since the system may fluctuate between the two equivalent magnetic arrangements.
To study how the electronic structure, and eventually, the electron-phonon interaction varies with doping we restrict ourself to the antiferromagnetic collinear arrangement since the non-collinear one has a similar electronic structure (see Yndurain and Soler (2009)) and its calculation is much more computational demanding. To obtain the electronic structure we have performed density functional Hohenberg and Kohn (1964); Kohn and Sham (1965) calculations using the SIESTA code Ordejon et al. (1996); Soler et al. (2002) which uses localized orbitals as basis functions Sankey and Niklewski (1989). In our calculation we use a double polarized basis set, non-local norm conserving pseudopotentials and a local density approximation (LDA) for exchange and correlation. The calculations are performed with stringent criteria in the electronic structure convergence (down to in the density matrix), Brillouin zone sampling (up to 18000 -points), real space grid (energy cut-off of 500 Ryd) and equilibrium geometry (residual forces lower than eV/Å). Due to the rapid variation of the density of states at the Fermi level, we used a polynomial smearing method Methfessel and Paxton (1989). To simulate the effect of doping we use the virtual crystal approximation (VCA) Bellaiche and Vanderbilt (2000).
The calculated total densities of states for various fluorine contents are drawn in Figure 2. We first observe the antiferromagnetic pseudo gap at the Fermi level for x=0. As the excess of electrons increases with the fluorine concentration the size of the gap decreases and eventually disappears between x=0.25 and x=0.30. It is interesting to notice that the antiferromagnetic gap lies fully bellow the Fermi level for . Moreover, the peak in the density of states at the antiferromagnetic gap upper edge crosses the Fermi level between x=0.10 and x=0.15.
Before we go into the discussion of the electron-phonon interaction parameter we analyze to what extent the virtual crystal approximation is a good starting point to study the electronic structure of these materials. To this end, we have calculated the equilibrium atomic structure and the electronic density of states for various super-cells size substituting an oxygen atom by a fluorine one. In this manner we have performed ”exact” calculations for x=0.25, 0.125 and 0.0625. The calculated density of states after full relaxation of the atomic positions are shown in Figure 3 which shows an excellent agreement with the VCA results. The VCA density of states reproduces the main features of the super-cell calculations, in particular, the peak in the density of sates crossing the Fermi level at around x=0.125. From these results it is clear that the details of how impurities are incorporated in the oxygen layer does not affect much the electronic structure around the Fermi level which is mainly dominated by the Fe electrons. The main effect of doping impurities at the oxygen layer is to pump electrons into the As-Fe-As layer, as anticipated Takahashi et al. (2008) .
We have also calculated the influence of an F impurity on the electronic charges and magnetic moments of the Fe atoms. In Table 1 we show the results for a supercell FeAsLaOF which corresponds to x=0.0625. We notice that the extra charge is distributed rather evenly through the Fe atoms and does not distinguish between those closer to the impurity and the other Fe atoms. On the other hand, although the magnetic moment of the Fe atoms closer to the F impurity display a magnetic moment larger than the others, the difference is rather small and we can rule out any relevant polaron formation. The antiferromagnetic coupling between successive Fe layer agrees with the experimental magnetic order of this material. The above results stress the appropriateness of the VCA approximation. Also, there is no significant influence of the presence of F impurities in the geometry of the Fe-As layers next to it. We find that the Fe-As distance depends on the extra charge at the As-Fe-As layers, rather than on the details concerning the oxygen layer and the substitutional impurities.
|* Fe atoms closest to the F impurity|
|Q: 8.301; : 1.247|
To address the electron-phonon interaction in this system and the superconducting parameter, we have considered the symmetric out-of-plane As mode at , in which the Fe atoms remain fixed whereas the As atoms move perpendicularly to the FeAs layers, expanding and compressing the Fe-As bonds. For this phonon mode the diagonal electron-phonon matrix element can be written Khan and Allen (1984):
where is the polarization vector, is the frequency of the vibrational mode involved and is the deformation potential for the state
The parameter associated to a vibrational mode can be written as a sum of the contributions of the electronic bands crossing the Fermi level in the form:
where is the nth-band density of states at the Fermi level, the phonon frequency, the mass of the arsenic atoms and is the variation of the energy band at the Fermi energy the phonon with the phonon amplitude:
Before we calculate the parameter it is worth looking at the deformation potential associated to the phonon. Although a supercell approximation could be more precise Noffsinger et al. (2009), the calculations are performed within the VCA. In Figure 4 we show, for the fluorine concentration of 10%, how the electronic band structure is perturbed by the presence of the phonon. In the case of a non magnetic calculation (Figure 4 (a)) the effect of the phonon in the bands at the Fermi level is small in all the bands but one (). On the contrary, in the case of the calculation for the antiferromagnetic configuration (Figure 4 (b)), all the bands are perturbed by the presence of the phonon, and, in addition, there is a flattening of the bands due to the antiferromagnetic pseudo-gap. The effect of the phonon is to modulate the size of the Fe magnetic moment with the vibration amplitude and therefore inducing an additional shift to the bands. A similar behavior was previously discussed Yndurain and Soler (2009).
We have calculated the parameter associated to this phonon in the approximation discussed above and within the VCA. The results, along with the variation of the density of states at the Fermi level, are shown in Figure 5, where a phonon energy of 25 meV is considered. We observe that the antiferromagnetic pseudo-gap induces an increase of the density of states at the Fermi energy with respect to the non-magnetic calculation (see Figure 2). This increase would enhance the parameter by a factor of two at the most. However, if the deformation potential is included, the parameter can be one order of magnitude larger than its value for the non-magnetic calculation (Figure 5 (b)). In addition, if we consider the renormalization of the parameter due to the renormalization of the phonon mode then which for the maximum value at x=0.125 becomes . This is indeed a substantial increase with respect to the non-magnetic calculation but not large enough to explain by itself the superconducting critical temperature. All the vibrational modes should be included in the calculation to obtain the total parameter. The connection between the Fe-As bond length and the magnetic moments in Fe atoms is important (as discussed by Yildirim Yildirim (2009)) and is relevant in the electron-phonon interaction.
In summary, we propose that in the superconducting phase of LaFeAsOF, the Fe atoms have a finite magnetic moment that fluctuates between two equivalent collinear antiferromagnetic configurations. The presence of this moment varies substantially the density of states near the Fermi level and enhances dramatically the electron-phonon interaction at least fort the As phonon mode. The maximum of the calculated parameter takes place at x=0.125. These results suggest that electron-phonon interaction, coupled with the Fe magnetic moments, must be carefully revised, before ruling out its connection with Fe-based compounds superconductivity.
I am indebted to J.M. Soler, for very enlightening discussions and for a critical reading of the manuscript. I also thank E. Bascones, M. J. Calderon, M.L.Cohen, G. Gomez-Santos, S.G. Louie, D. Sanchez-Portal and B. Valenzuela for very helpful comments about this work. This work was supported by the Spanish Ministry of Science and Innovation through grants FIS2009-12712 and CSD2007-00050.
- Y. Kamihara, H. Hirmatsu, M. Hirano, R. Kawamura, H. Yanagi, T. Kamiya, and H. Hosono, J. Am. Chem. Soc. 128, 10012 (2006).
- Y. Kamihara, T. Watanable, M. Hirano, and H. Hosono, J. Am. Chem. Soc. 130, 3296 (2008).
- M. R. Norman, Physics 1, 21 (2008).
- I. I.Mazin, M. D.Johannes, L. Boeri, K. Koepernik, and D. J. Singh, Phys. Rev. B 78, 085104 (2008).
- L. Boeri, O. V. Dolgov, and A. A. Golubov, Phys. Rev. Lett. 101, 026403 (2008).
- L. Boeri, M. Calandra, I. I.Mazin, O. V.Dolgov, and F. Mauri, Phys. Rev. B 82, 020506 (2010).
- M. A. Korotin, S. V. Streltsov, A. O. Shorikov, and V. I. Anisimov, Journal of Experimental and Theoretical Physics 107, 649 (2008).
- E. Bascones, M. J. Calderon, and B. Valenzuela, Phys. Rev. Lett. 104, 227201 (2010).
- C.-C. Lee, W.-G. Yin, and W. Ku, Phys. Rev. Lett. 103, 267001 (2009).
- S. Hahn, Y. Lee, N. Ni, P. Canfield, A. I. Goldman, R. McQueeney, B. Harmon, A. A, B.M-Leu, E. Alp, et al., Phys. Rev. B 79, 220511(R) (2009).
- M. Le Tacon, T. Forrest, C. Ruegg, A. Bosak, A. Walters, R. Mittal, H. M. Ronnow, N. Zhigadlo, S. Katrych, J. Karpinski, et al., Phys. Rev. B 80, 220504(R) (2009).
- D. Reznik, K. Lokshin, D. C. Mitchell, D. Parshall, W. Dmowski, D. Lamago, R. Heid, K.-P. Bohnen, A. S. Sefat, M. A. McGuire, et al., Phys. Rev. B 80, 214534 (2009).
- F. Yndurain and J. M. Soler, Phy. Rev. B 79, 134506 (2009).
- P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964).
- W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965).
- P. Ordejon, E. Artacho, and J. M. Soler, Phys. Rev. B 53, R10441 (1996).
- J. Soler, E. Artacho, J. Gale, A. Garcia, J. Junquera, P. Ordejon, and D. Sanchez-Portal, J. Phys.: Condens. Matter 14, 2745 (2002).
- O. F. Sankey and D. J. Niklewski, Phys. Rev. B 40, 3979 (1989).
- M. Methfessel and A. T. Paxton, Phys. Rev. B 40, 3616 (1989).
- L. Bellaiche and D. Vanderbilt, Phys. Rev. B 61, 7877 (2000).
- H. Takahashi, K. Igawa, K. Arii, Y. Kamihara, M. Hirano, and H. Hosono, Nature 453, 376 (2008).
- F. S. Khan and P. B. Allen, Phys. Rev. B 29, 3341 (1984).
- J. Noffsinger, F. Giustino, S. G. Louie, and M. L. Cohen, Phys. Rev. Lett. 102, 147003 (2009).
- T. Yildirim, Phys. Rev. Lett. 102, 037003 (2009).