# Spin-Lattice Coupling and Superconductivity in Fe Pnictides

Authors: T. Egami; B. V. Fine; D. Parshall; A. Subedi; D. J. Singh
Journal: Advances in Condensed Matter Physics
 (2010)

---

## Abstract

We consider strong spin-lattice and spin-phonon coupling in iron pnictides and discuss its implications on superconductivity. Strong magneto-volume effect in iron compounds has long been known as the Invar effect. Fe pnictides also exhibit this effect, reflected in particular on the dependence of the magnetic moment on the atomic volume of Fe defined by the positions of the nearest neighbor atoms. Through the phenomenological Landau theory, developed on the basis of the calculations by the density functional theory (DFT) and the experimental results, we quantify the strength of the spin-lattice interaction as it relates to the Stoner criterion for the onset of magnetism. We suggest that the coupling between electrons and phonons through the spin channel may be sufficiently strong to be an important part of the superconductivity mechanism in Fe pnictides.


---

## Body

## 1. Introduction

Recent discovery of superconductivity in Fe pnictides [1] with the critical temperature TC up to 55 K [2] caused enormous excitement in the field, for various reasons. First, this is the first noncuprate family of superconductors with TC above 40 K. Second, superconductivity appears when the antiferromagnetic (AFM) order is suppressed by doping [3], just as in the cuprates. Third, unlike the cuprates, strong electron correlations are not observed by spectroscopy [4], suggesting that the Mott physics may not be a necessary ingredient for the mechanism of high-temperature superconductivity. Finally, there is a large family of similar compounds that show superconductivity, making experimental research less restricted by chemical or materials issues. The field is making surprisingly fast development, partly because of the accumulated experience of working on the cuprates. It is possible that the origin of thesuperconductivity in this family of compounds may be easier to identify than for the cuprates, and the success of solving this problem hopefully will facilitate understanding of the cuprate problem.Even though the AFM order is suppressed by doping, spins are active in the doped Fe pnictide superconductors. Again, just as in the cuprates strong magnetic excitations, including the so-called resonance peak, are observed by inelastic neutron scattering [5–8]. Core level spectroscopy is consistent with the local Fe moment of about 1μB [4]. Interestingly, the density functional theory (DFT) calculations always predict the magnetic ground state (AFM or incommensurate order) for the experimental lattice constants [9–12]. The phonon dispersions observed by inelastic X-ray scattering are consistent with the DFT calculations only when the magnetic ground state is assumed [13]. Only in the collapsed phase of CaFe2As2, in which the c-axis lattice constant is reduced by as much as 10% compared to the magnetic state, does the material become truly spin-degenerate [14]. All these observations strongly suggest that the superconducting samples are locally and dynamically spin-polarized, and show strong dynamic Fe spin fluctuations. Although we do not have precise knowledge of their spin dynamics, judged from the absence of strong quasielastic scattering in neutron scattering with energy resolution of 1 meV, the time-scale of fluctuation must not be slower than 1 ps. This result supports the view that spins are involved in the mechanism of superconductivity, for instance though the spin-fluctuation mechanism [15–17].However, there are many puzzling, important questions that need to be answered before addressing the question of the mechanism: the first puzzle is the effect of doping. In the cuprates, doping is necessary for introducing mobile charge carriers, since the parent compounds are Mott-Hubbard insulators. In Fe pnictides, on the other hand, the parent compounds are already metallic, and doping does not appear to change the charge carrier density very much [9]. Rather, the main effect of doping is to suppress the AFM ground state. In the pnictide parent compounds strongly two-dimensional spin fluctuations are observed above TN [18], just as in the superconducting Fe pnictides [6, 8]. However, whereas the LaFeAsO (1111) type compounds are strongly two-dimensional [10], BaFe2As2(122) type compounds are much more three-dimensional [19]. The second curious behavior is that the observed magnetic moment on the AFM phase varies from a compound to compound, but is always significantly smaller than predicted by the DFT calculations [10–12]. Third, there is an interesting interplay between both the lattice and magnetism [20], and the lattice and superconductivity [21, 22]. In this article, we focus on the third point, that of the lattice effect. For the purpose of highlighting the essence of the effect, we use simple approximations, namely the Landau theory and the Stoner theory, using the results of the LDA calculations as a guide. We argue that through the spin-lattice coupling effect the lattice may play a much larger role than generally acknowledged in determining the properties of Fe pnictides, possibly including even the superconductivity.


## 2. Magneto-Volume Effect in Fe Pnictides

### 2.1. Dependence of Fe Moment on the Structure

It has long been known that the magnetic moment of transition metals depends on volume [23]. Because of the Pauli exclusion principle, the electron kinetic energy of the spin-polarized state is higher for parallel spins if the volume is the same, and volume expansion relaxes the kinetic energy. In some iron alloys, thermal volume expansion due to lattice anharmonicity cancels the decrease in volume associated in the decrease in spin-splitting, resulting in zero thermal expansion, widely known as the Invar behavior. The negative or zero thermal expansion is indeed observed for PrFeAsO [24]. The collapsed phase of CaFe2As2 is a dramatic case of such a magneto-volume effect. This compound shows AFM order below 140 K, but with the pressure of 0.4 GPa it undergoes the first-order phase transition into a nonmagnetic phase with the reduction in volume of 5% [14].In Fe pnictides, layers of Fe atoms are sandwiched by layers of pnictide such as As or P [1]. Thus if the layer-layer distance of pnictide is changed, the magnetic moment of Fe is strongly affected. This coupling of the Fe moment to the pnictide position in the lattice was recognized early by the DFT calculations [9, 20]. Figure 1 shows the calculated dependence of the Fe moment in Ba(Fe0.92Co0.08)2As2 on the separation between the As layer and the Fe layer (z). The calculations were done within virtual crystal approximation (VCA) and local density approximation (LDA) with general potential linearized augmented plane-wave (LAPW) method [25], including local orbitals [26]. LAPW sphere radii of 2.2a0,2.0a0,and2.0a0, where a0 is Bohr radius, were used for Ba, Fe, and As, respectively. To account for Co doping, an electron number of 26.08 was used for Fe. We used the experimentally reported tetragonal lattice parameters (a = 3.9625 Å, c = 13.0168 Å) [27]. In the calculated result, clearly there is a quantum critical point (QCP) for magnetism near zc=1.20Å as shown in Figure 1. The local exchange interaction is strong enough to spin-split the band by overcoming the kinetic energy cost only for z>zc. We obtain similar results from the calculation on undoped BaFe2As2. Compared also with other data [12, 20], the relation between z and Fe moment M appears to be rather insensitive to compositions, and the relation shown in Figure 1 appears to be a nearly universal property of the FeAs triple layer. This must be because the in-plane lattice constant, thus the Fe-Fe distance, is very similar within ±1% among many Fe pnictide compounds. Thus the parameter z, the As-Fe layer separation, is a good common measure of the magneto-volume effect. For instance, in the collapsed phase of CaFe2As2 the value of z is 1.23Å [14], close to the value of zc in Figure 1.Figure 1
The magnetic moment,M, as a function of the Fe-As layer separation, z, calculated for Ba(Fe0.92Co0.08)2As2. The solid curve is a fit by (3).


### 2.2. Landau Theory

Let us develop a Landau-type theory to describe the dependence of the local magnetic moment,M, on the As-Fe layer separation, z. We may write the magnetic free energy as(1)FM=AM2+BM4+Fs-l,
where Fs-l is the spin-lattice interaction energy expanded by z-zc,
(2)Fs-l=[α(z-zc)+β(z-zc)2]M2,
where α<0.Fs-l is negative only for z>zc. We retain only the terms with even powers of M because of the symmetry. We set A=0 so that zc is the critical point. By minimizing FM with respect to M we obtain,
(3)M2=|α|2B[(z-zc)+βα(z-zc)2].
The dashed line in Figure 1 gives the fit of this equation to the calculated moment. From this fit we obtain zc=1.20Å, α/2B=19.16μB2/Å, and α/β=-1.40Å. In CaFe2As2, the QCP is hidden because the nature of the transition is strongly first order. In the vicinity of zc, the moment is proportional to (z-zc)1/2.Recently, the magnetic moment was experimentally determined for the series of compounds,CeFe(As1-xPx)O [28]. Because As and P are isovalent, the replacement of As by P does not dope carriers to the system, but merely change the Fe(As/P) layer separation. Indeed (3) fits the results nicely as shown in Figure 2 [29]. Since the range of z is relatively, small β was neglected. The fit to the data gave zc=1.278Å,α/2B=11.67(μB/Å). Theoretical calculations for the CeFe(As,P)O system are in progress. Note that in the calculation presented in Figure 1, the lattice is assumed to be static. In reality, however, the zero-point phonons are present, with the amplitude of the order of 0.1 Å. This should shift the critical point upward, in better agreement with the data in Figure 2.Figure 2
Dependence of the experimentally determined magnetic moment on the Fe-As layer separation for CeFe(As,P)O (open circles) [28], and the fit by (3) (solid curve).


### 2.3. Phonon Softening

We now add the lattice elastic energy in order to consider the phonon softening due to the magneto-volume effect. The phonon to be considered here is the As Raman mode, in which As layers move against each other along thec-axis, either toward or away from the Fe layer. The magneto-elastic free energy is
(4)FS(z)=[α(z-zc)+β(z-zc)2]M2+K(z-z1)2=-α2(z-zc)22B[1+βα(z-zc)]2+K(z-z1)2.
By minimizing FS with respect to z we obtain the As position, zM,
(5)∂FS∂z|z=zM=-α2(zM-zc)B[1+2βα(zM-zc)](1+βα(zM-zc))+2K(zM-z1)=0.
Then by taking the second derivative we obtain the elastic stiffness renormalized by the spin-lattice interaction,
(6)K′=K{1-α22KB[1+6βα(zM-zc)+6(βα)2(zM-zc)2]}.
As shown in Figure 3, allowing spin polarization softens the As Raman phonon frequency by 15% at zM=1.36Å. This effect was noted earlier [12], and agrees with the experimental observations [13]. Note that the energy minimum of the DFT calculation systematically underestimates the lattice constant. Thus we obtain, α=-0.193eV/ÅμB2, β=0.137eV/Å2μB2, and z1=1.32Å. The value of z1 is in good agreement with the LDA calculation (1.32Å after correcting for the systematic underestimation), proving the internal consistency. Thus this theory elucidates how the magneto-volume interaction, (2) can induce magnetization as a function of the As-Fe layer separation, and softening of the As Raman phonon mode.Figure 3
Electron energy as a function of the Fe-As layer separation is calculated by LDA forBa(Fe0.92Co0.08)2As2, with spin polarization; the elastic constant is reduced by 29% compared to the state without spin polarization, accounting for the phonon softening by 15% due to spin polarization.


## 2.1. Dependence of Fe Moment on the Structure

It has long been known that the magnetic moment of transition metals depends on volume [23]. Because of the Pauli exclusion principle, the electron kinetic energy of the spin-polarized state is higher for parallel spins if the volume is the same, and volume expansion relaxes the kinetic energy. In some iron alloys, thermal volume expansion due to lattice anharmonicity cancels the decrease in volume associated in the decrease in spin-splitting, resulting in zero thermal expansion, widely known as the Invar behavior. The negative or zero thermal expansion is indeed observed for PrFeAsO [24]. The collapsed phase of CaFe2As2 is a dramatic case of such a magneto-volume effect. This compound shows AFM order below 140 K, but with the pressure of 0.4 GPa it undergoes the first-order phase transition into a nonmagnetic phase with the reduction in volume of 5% [14].In Fe pnictides, layers of Fe atoms are sandwiched by layers of pnictide such as As or P [1]. Thus if the layer-layer distance of pnictide is changed, the magnetic moment of Fe is strongly affected. This coupling of the Fe moment to the pnictide position in the lattice was recognized early by the DFT calculations [9, 20]. Figure 1 shows the calculated dependence of the Fe moment in Ba(Fe0.92Co0.08)2As2 on the separation between the As layer and the Fe layer (z). The calculations were done within virtual crystal approximation (VCA) and local density approximation (LDA) with general potential linearized augmented plane-wave (LAPW) method [25], including local orbitals [26]. LAPW sphere radii of 2.2a0,2.0a0,and2.0a0, where a0 is Bohr radius, were used for Ba, Fe, and As, respectively. To account for Co doping, an electron number of 26.08 was used for Fe. We used the experimentally reported tetragonal lattice parameters (a = 3.9625 Å, c = 13.0168 Å) [27]. In the calculated result, clearly there is a quantum critical point (QCP) for magnetism near zc=1.20Å as shown in Figure 1. The local exchange interaction is strong enough to spin-split the band by overcoming the kinetic energy cost only for z>zc. We obtain similar results from the calculation on undoped BaFe2As2. Compared also with other data [12, 20], the relation between z and Fe moment M appears to be rather insensitive to compositions, and the relation shown in Figure 1 appears to be a nearly universal property of the FeAs triple layer. This must be because the in-plane lattice constant, thus the Fe-Fe distance, is very similar within ±1% among many Fe pnictide compounds. Thus the parameter z, the As-Fe layer separation, is a good common measure of the magneto-volume effect. For instance, in the collapsed phase of CaFe2As2 the value of z is 1.23Å [14], close to the value of zc in Figure 1.Figure 1
The magnetic moment,M, as a function of the Fe-As layer separation, z, calculated for Ba(Fe0.92Co0.08)2As2. The solid curve is a fit by (3).


## 2.2. Landau Theory

Let us develop a Landau-type theory to describe the dependence of the local magnetic moment,M, on the As-Fe layer separation, z. We may write the magnetic free energy as(1)FM=AM2+BM4+Fs-l,
where Fs-l is the spin-lattice interaction energy expanded by z-zc,
(2)Fs-l=[α(z-zc)+β(z-zc)2]M2,
where α<0.Fs-l is negative only for z>zc. We retain only the terms with even powers of M because of the symmetry. We set A=0 so that zc is the critical point. By minimizing FM with respect to M we obtain,
(3)M2=|α|2B[(z-zc)+βα(z-zc)2].
The dashed line in Figure 1 gives the fit of this equation to the calculated moment. From this fit we obtain zc=1.20Å, α/2B=19.16μB2/Å, and α/β=-1.40Å. In CaFe2As2, the QCP is hidden because the nature of the transition is strongly first order. In the vicinity of zc, the moment is proportional to (z-zc)1/2.Recently, the magnetic moment was experimentally determined for the series of compounds,CeFe(As1-xPx)O [28]. Because As and P are isovalent, the replacement of As by P does not dope carriers to the system, but merely change the Fe(As/P) layer separation. Indeed (3) fits the results nicely as shown in Figure 2 [29]. Since the range of z is relatively, small β was neglected. The fit to the data gave zc=1.278Å,α/2B=11.67(μB/Å). Theoretical calculations for the CeFe(As,P)O system are in progress. Note that in the calculation presented in Figure 1, the lattice is assumed to be static. In reality, however, the zero-point phonons are present, with the amplitude of the order of 0.1 Å. This should shift the critical point upward, in better agreement with the data in Figure 2.Figure 2
Dependence of the experimentally determined magnetic moment on the Fe-As layer separation for CeFe(As,P)O (open circles) [28], and the fit by (3) (solid curve).


## 2.3. Phonon Softening

We now add the lattice elastic energy in order to consider the phonon softening due to the magneto-volume effect. The phonon to be considered here is the As Raman mode, in which As layers move against each other along thec-axis, either toward or away from the Fe layer. The magneto-elastic free energy is
(4)FS(z)=[α(z-zc)+β(z-zc)2]M2+K(z-z1)2=-α2(z-zc)22B[1+βα(z-zc)]2+K(z-z1)2.
By minimizing FS with respect to z we obtain the As position, zM,
(5)∂FS∂z|z=zM=-α2(zM-zc)B[1+2βα(zM-zc)](1+βα(zM-zc))+2K(zM-z1)=0.
Then by taking the second derivative we obtain the elastic stiffness renormalized by the spin-lattice interaction,
(6)K′=K{1-α22KB[1+6βα(zM-zc)+6(βα)2(zM-zc)2]}.
As shown in Figure 3, allowing spin polarization softens the As Raman phonon frequency by 15% at zM=1.36Å. This effect was noted earlier [12], and agrees with the experimental observations [13]. Note that the energy minimum of the DFT calculation systematically underestimates the lattice constant. Thus we obtain, α=-0.193eV/ÅμB2, β=0.137eV/Å2μB2, and z1=1.32Å. The value of z1 is in good agreement with the LDA calculation (1.32Å after correcting for the systematic underestimation), proving the internal consistency. Thus this theory elucidates how the magneto-volume interaction, (2) can induce magnetization as a function of the As-Fe layer separation, and softening of the As Raman phonon mode.Figure 3
Electron energy as a function of the Fe-As layer separation is calculated by LDA forBa(Fe0.92Co0.08)2As2, with spin polarization; the elastic constant is reduced by 29% compared to the state without spin polarization, accounting for the phonon softening by 15% due to spin polarization.


## 3. Generalized Stoner Condition

We now turn our attention to the Fe band splitting. We start with the classical Stoner theory of itinerant ferromagnetism [30], and generalize it to an antiferromagnet [31–33]. Even though more complex and accurate DFT calculations can be made, this is a useful exercise for simplicity and clarity of logic. In the Stoner theory the electron energy is given by(7)E=∫0εF↑εN↑(ε↑)dε↑+∫0εF↓εN↓(ε↓)dε↓+In↑n↓,
where N↑(ε) is the electron density of states for up spins, n↑ is the density of electrons with up spin, and ε↑ is the energy of an electron with up spin. Magnetization is given by M=n↑-n↓. If we start with the nonmagnetic state and introduce a small spin splitting by shifting the up spin by dε↑, the energy changes(8)dE=2N(εF)[1-IN(εF)]dε↑2,
which gives the Stoner criterion,(9)IN(εF)≥1.
Near the Stoner QCP,(10)M=dn↑-dn↓=2dn↑=2N(εF)dε↑,
thus(11)dE=1-IN(εF)2N(εF)M2,
and in (1),(12)A=1-IN(εF)2N(εF).
Thus A=0 at the QCP.In the case of an antiferromagnet, the same argument can be constructed for the staggered spin-split bands, except that the exchange energy termI includes the kinetic energy cost for modulating the spin split band by K=π/a [32, 33]. Therefore, we can interpret zc as the generalized Stoner QCP. In general for a spin density wave (SDW) ordering with Q, within the random phase approximation (RPA), the magnetic moment above the SDW ordering temperature is given by(13)M(Q)=χ0(Q)1-Iχ0(Q)hQ,
where M(Q) is the SDW amplitude, χ0(Q) is the generalized bare susceptibility, and hQ is the SDW mean-field [31]. Thus the generalized Stoner condition is given by [31–33](14)Iχ0(Q)≥1.
Note that χ0(Q=0)=N(εF), which recovers the Stoner condition for ferromagnetism. Now the effect of lattice strain is to modify the generalized Stoner condition. For instance, reducing the As-Fe separation increases the bandwidth, thus decreases the generalized susceptibility. This brings the system closer to the Stoner QCP, and reduces the magnetic moment. Because the Fe pnictides are close to the Stoner QCP, even a small lattice strain, including phonons, can have very significant effects. Similarly, the in-plane Fe phonon mode that is relevant to the tetragonal-orthorhombic structural phase transition could also be important to the magnetism of Fe pnictides. The strong anisotropy of the spin exchange in the Fe plane [34] is indicative of the strong spin-lattice coupling in the Fe plane.


## 4. Relation to Superconductivity

One of the most intriguing lattice effects on the superconductivity of Fe pnictides is the dependence of the critical temperature,TC, on the geometry of the FeAs4 tetrahedron [21, 22]. Data show that TC is strongly related to the As-Fe-As angle. Because the As-Fe-As angle is directly related to the Fe-As layer separation, z, in Figure 4, we plotted TC as a function of z, using the published results of crystallographic analysis. The results are shown also in Table 1, with references for the data. Clearly the behavior above za=1.4Å is different from that below za. Above zaTC is not much dependent on z, whereas there is strong dependence below. It is possible that the systems with z>za are regular BCS superconductors, and the enhancement is present only for z<za. Below za, there is recognizable correlation between TC and z, except for a few outliers. This correlation supports the idea that z is a universal parameter for the properties of FeAs compounds, including the superconductivity, regardless of the composition.Table 1
The values ofz, TC for various Fe pnictides. Na-111 means NaFeAs, K-Ba122 means K1-xBaxFe2As2, and so forth.
Compoundz(Å)TC (K)Ref.Na1111.4115.5[35]K-Ba1221.3738[36]F-La11111.3426[37]F-Nd11111.3746[38]F-Ce11111.3435[39]F-Pr11111.3247[40]F-La11111.3220[41]Pr11111.3247[40]O-Sm11111.3734[42]V-Nd11111.3851[43]Co-La11111.3414.3[44]Fe(Se0.416Te0.584)1.6514[45]Fe(Se0.493Te0.507)1.6214[45]Li1111.5110[46]Li1111.6818[47]F-Sm11111.3646[48]F-Tb11111.3844[49]V-Nd11111.3751[43]Co-SrFeAsF1.354[50]K-Ba1221.3838[36]Figure 4
Superconducting critical temprtureTC, as a function of the As-Fe layer separation, z.It is interesting to note that, from (2) and (3), Fs-p=-2BM4. Thus we plotted TC against M4 calculated by (3), not the experimental values of M, in Figure 5. Again a linear correlation appears to be present, although the correlation is far from perfect. On the other hand, the relation, TC∝M2, is observed if we plot TC against the moment experimentally determined for the CeFe(As, P)O system shown in Figure 2 [29]. An obvious implication of these correlations is that indeed magnetism is deeply involved in superconductivity, even though there is no static magnetic order in most of the superconducting samples. However, the involvement of spins in the mechanism may not be limited to the spin-fluctuation mechanism. The spin-lattice coupling could be involved in the superconductivity mechanism through the electron-phonon (e-p) coupling in the spin-channel.Figure 5
Superconducting critical temperature,TC, against M4 calculated by the LDA.For instance, the phonon can modify the electron hopping integraltij [51]. In the presence of antiferromagnetic correlation, the phonon-induced charge transfer from one spin sublattice to the other can be given by
(15)H=∑q,k,sΔtq,k[(aq+fk,s+gk+q,s+aqfk,s+gk-q,s)+H.c.],
where s refers to spin (↓or↑), Δtq,k is the derivative of t with respect to phonon displacement, aq is a phonon annihilation operator, and fk and gk are electron annihilation operators for the spin ↑ sublattice and the spin ↓ sublattice. The ↑ spin is a minority spin for the spin ↓ sublattice, so (15) results in the spin transfer between the majority band and the minority band. Note that this coupling is different from the usual spin-phonon coupling through the modification of the exchange integral J by phonons, and ultimately we will have to consider both [52]. Even though the conventional charge-channel e-p coupling (Fröhlich coupling) is weak for the Fe pnictides [53], the e-p coupling through the spin-channel, such as the one in (15), may be strong enough to make a difference. Other possibility includes the mechanism proposed by Kulic and Haghighirad [54] in which the on-site Hubbard repulsion, U, is modified by the dipolar moment of pnictide ions.Possible importance of the spin-phonon coupling in the cuprates was discussed earlier [55–57]. Whereas many people speculate that the spin-fluctuation is the common mechanism for the cuprates as well as for Fe pnictides, it is possible that the spin-channel e-p coupling is the common thread for high-temperature superconductors, although much remain to be researched before we reach any definite conclusion.


## 5. Conclusions

The conventional electron-phonon (e-p) coupling through the charge channel is quite small for the Fe pnictide compounds [53]. This led many to conclude that lattice and phonons are irrelevant to the superconductivity of the Fe pnictides. Consequently, the spin-fluctuation mechanism is regarded to be the leading mechanism to explain their high TC [15]. However, the lattice is intimately involved in the magnetism of this compound through the magneto-volume effect, as shown in this paper. This coupling has been known for a long time as Invar effect. The lattice controls the Stoner condition, and thus the onset of spin-splitting. Because the lattice is so intimately involved in magnetism, it is furthermore possible that the e-p coupling through the spin-channel is relevant, for instance involving the As Raman phonon mode or the in-plane Fe mode. The lattice effect may be much more important than generally assumed to the properties of Fe pnictides, including superconductivity.



---

*Source: 164916-2010-03-08.xml*