Electronic and lattice dynamical properties of the iron-based superconductors LiFeAs and NaFeAs

The electronic structure and lattice vibrational frequencies of the newly discovered superconductors, LiFeAs and NaFeAs, are calculated within density functional theory. We show that in the vicinity of the Fermi energy, the density of states is dominated by contributions from Fe 3d states. We also calculate the electron-phonon coupling strength and show that it is too weak to account for the observed values of the superconducting transition temperatures. This seems to indicate that the iron-based superconductors are not of the conventional type.

leads to the reasonable belief that strong electronic correlations are important in these systems, and that superconductivity in these compounds is somehow connected to magnetic fluctuations. 31,32,33,34,35,36,37,38,39,40,41,42,43 Indeed, the electron-phonon coupling in LaOFeAs was estimated to be too small 44 to give rise to superconductivity within the conventional BCS formulation. 45 Recently, a third class of iron-based superconductors was discovered. LiFeAs and NaFeAs were found to superconduct below 18 K and 9 K, respectively. 46,47,48,49 It turns out that in these two compounds, no magnetic order is detected at all temperatures. In some sense, these two compounds are important with regards to understanding the mechanism of superconductivity in ironbased superconductors. the absence of spin density wave (SDW) transition, on the one hand, and the relatively low T c in comparison with the first two classes of iron-based superconductors, on the other hand, make these two compounds possible candidates for being conventional BCS superconductors.
Band structure calculations, using local density approximation (LDA) within density functional theory (DFT), were recently reported for LiFeAS. 50,51 It was found that LiFeAs is semi-metallic, and that the density of states (DOS) near the Fermi level is dominated by the Fe 3d states. Thus, the electronic structure of stoichiometric LiFeAs is similar to that of the parent compounds of the first class, with a hole cylinder at the Brillouin zone (BZ) center, electron cylinders at the BZ corners, and an electronic DOS that decreases strongly with increasing energy in the vicinity of the Fermi energy.
In this work we report DFT calculations of the electronic and lattice properties of LiFeAs and NaFeAs. In particular, we calculate the electron-phonon coupling strength and show that it is too weak to account for the superconducting transition temperatures observed in these compounds. Our calculations, together with previous calculations 44 of the electron-phonon coupling strength in LaOFeAs, seem to indicate that iron-based superconductors are not of the conventional type. arXiv:0812.1215v1 [cond-mat.supr-con] 5 Dec 2008

II. METHOD
The electronic structure calculations are carried out using the all-electron full-potential linear augmented plane wave (FP-LAPW) method as implemented in WIEN2K code. 52 The exchange-correlation potential was calculated using the generalized gradient approximation (GGA) as proposed by Pedrew, Burke, and Ernzerhof (PBE). 53 The radii of the muffin-tin spheres for the various atoms were chosen so that the neighboring spheres almost touch each other. We set the parameter R M T K max =7, where R M T is the smallest muffin-tin radius and K max is a cutoff wave vector. The valence wave functions inside the muffin-tin spheres are expanded in terms of spherical harmonics up to l max = 10, and in terms of plane waves with a wave vector cutoff K max in the interstitial region. The charge density is Fourier expanded up to G max =13a 0 −1 , where a 0 is the Bohr radius. Convergence of the self-consistent field calculations is attained with a total energy convergence tolerance of 0.01 mRy.
The calculation of the frequencies of the vibrational modes and the electron-phonon coupling parameter was performed using ultrasoft pseudopotentials and an expansion of the wave function of the valence electrons in terms of plane waves, with an energy cutoff of 30 Rydbergs. 54 In both the electronic and lattice calculations, the experimental values of the low-temperature lattice constants and atomic positions 46,49 are used. For both compounds, the crystal is tetragonal with space group P4/nmm. In LiFeAs, the lattice constants are a = 3.76982Å, c = 6.30693Å, whereas in NaFeAs, a = 3.94729Å, and c = 6.99112Å.

III. RESULTS AND DISCUSSION
Our results for the electronic structure calculations for LiFeAs and NaFeAs are summarized in figures 1 and 2, respectively, where the electronic density of states (DOS) is displayed. For LiFeAs, our calculated DOS is similar to that reported earlier. 50,51 In both LiFeAs and NaFeAs, the DOS plots show some generic features that are common to the parent compounds of the iron-based superconductors: A DOS that is dominated by the Fe 3d states in the vicinity of the Fermi energy, with only a small contribution from the As and alkali metal states, and that is strongly decreasing with energy near the Fermi energy. It should be noted that in a unit cell of LiFeAs, for example, there are two Li, two Fe, and two As atoms. Thus to get the total atomic DOS, the values of the atomic DOS shown in 1 should be multiplied by 2. The total DOS is the sum of the total atomic DOS and the DOS in the interstitial region.
Since the DOS at the Fermi energy, N(E F ), is ∼4 states/eV in both LiFeAs and NaFeAs, which is not very small, and because of the relatively lower T c compared with the other iron-based superconductors, one may won-der whether electron-phonon coupling may lie behind the mechanism for superconductivity in LiFeAs and NaFeAs. This notion may be given more credence by the observation that a sodium atom is about three times more massive than a lithium atom, so that if the attractive electron-electron interaction is mediated by the alkalimetal atomic vibrations, then this difference in the mass could explain the difference in the values of T c between the two compounds via the well-known isotope effect.
To test this idea, we carried out a calculation of the phonon dispersion curves and the electron-phonon coupling strength in these compounds. Since the crystallographic point group in LiFeAs and NaFeAs is D 4h , the vibrational modes at Γ, the BZ center, are decomposed according to the following irreducible representations The acoustic modes, with vanishing frequency at Γ, the BZ center, transform according to the A 2u and E u irreducible representations. Excluding the acoustic modes, we are left with 15 modes with nonzero frequencies; among these, the symmetric ones are Raman-active, while the antisymmetric modes are infrared-active. The calculated frequencies of the Raman-and infrared-active modes at the Γ point of the BZ are given in Table I. The phonon dispersion curves in LiFeAs, plotted along high symmetry directions in the BZ, are shown in Fig. 3, and the corresponding curves in NaFeAs are given in Fig. 4. Our results for the phonon frequencies at the BZ center may be checked by Raman scattering and infrared absorption experiments, while the phonon dispersion curves may be checked by neutron scattering measurements.
We calculated the electron-phonon coupling parameter λ and found it to be 0.29 and 0.27 for LiFeAs and NaFeAs, respectively. For conventional superconductors, where the attractive electron-electron interaction is mediated by phonons, the transition temperature is given by the Allen and Dynes 55 modified McMillan's formula 56 where ω log is the logarithmic average phonon frequency, expressed in degrees Kelvin, and µ * is the Coulomb pseudopotential parameter, usually taken to be ∼0.13. Since ω log in LiFeAs and NaFeAs is of the order of 100 K, the resulting value of T c is much less than 1 K. We conclude that the electron-phonon coupling is too weak to account for superconductivity in this class of iron-based superconductors.
There are some puzzling questions that beset this third class of iron-based superconductors.
In the parent compounds of the first two classes of ironbased superconductors, magnetic order is established at low temperatures, where the Fe magnetic moments adopt a collinear antiferromagnetic (c-AFM) order; this is unequivocally confirmed by both neutron  diffraction measurements 18,19,20,21 and DFT calculations. 32,33,34,35,36,37,38,39,40 In the first class, it is only upon electron doping through the replacement of a small percentage of oxygen atoms with fluorine atoms, or the removal of a small percentage of oxygen atoms, that the magnetic order is suppressed and superconductivity emerges. In the second class of iron-based superconductors, magnetic order is suppressed by hole doping through the replacement of some alkaline earth atoms with alkali atoms. To better understand the situation in the third class of iron-based superconductors, we carried out spinpolarized DFT calculations on stoichiometric LiFeAs and NaFeAs, in addition to the calculations reported above for the nonmagnetic phases of these compounds. We considered, within GGA and GGA+U, various spin arrangements on the Fe sites. Similarly to the case of the first two classes, we find that the c-AFM phase, with a spin-stripes pattern, is indeed the lowest energy phase. Within GGA, the energy of the c-AFM phase in LiFeAs is lower than the AFM phase by 0.081 eV per Fe atom (eV/Fe), lower than the ferromagnetic (FM) phase by 0.085 eV/Fe, and lower than the nonmagnetic phase by 0.123 eV/Fe. For NaFeAs, the energy of the c-AFM phase is lower than the AFM phase by 0.048 eV/Fe, lower than the FM phase by 0.205 eV/Fe, and lower than the nonmagnetic phase by 0.182 eV/Fe. The differences are even greater within GGA+U, where onsite Coulomb repulsion is taken into account. Thus, according to our DFT calculations, stoichiometric LiFeAs and NaFeAs should be similar to the parent compounds of the first two classes, and they should not superconduct; instead, at low temperature, the stoichiometric compounds should display magnetic order. Deviations from stoichiometry, on the other hand, may suppress the magnetic order, making way for superconductivity, just like doping does in the first two classes of iron-based compounds. In the case of LiFeAs, it is indeed the case that the synthesized superconducting compounds were not stoichiometric, the chemical formula being Li 1−x FeAs. 48 The situation is less clear in the case of NaFeAs, 49 but we believe, on the basis of our spinpolarized calculations, and the absence of any detectable magnetic order at low temperatures, that the synthesized NaFeAs samples must also be nonstoichiometric, and that further studies on sample characterization are necessary in this regards.

IV. CONCLUSIONS
In conclusion, we have presented the results of electronic structure calculations on LiFeAs and NaFeAs, members of a new class of superconducting compounds. In similarity to other iron-based superconductors, the density of states in the vicinity of the Fermi energy is found to be dominated by contributions from the Fe 3d states. We have also calculated the Raman and infrared phonon frequencies at the Brillouin zone center, as well the phonon dispersion curves along high symmetry directions in the Brillouin zone. We have evaluated the electron-phonon coupling parameter in LiFeAs and NaFeAs, and found its value to be too small to account for the observed superconducting transition temperatures in these compounds. Our results, taken together with previous estimates of the electron-phonon coupling strength in LaOFeAs, seem to suggest clearly that iron-based superconductors are not of the conventional type, where the attractive electron-electron interaction is mediated by phonons.