Microsoft Word - 4002-11881-1-LE-F 225 American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) ISSN (Print) 2313-4410, ISSN (Online) 2313-4402 © Global Society of Scientific Research and Researchers http://asrjetsjournal.org/ Electronic and Superconductivity Properties of FeSe Superconductor Sumal Chandra * Department of Physics, Mawlana Bhashani Science and Technology University Santosh, Tangail-1902, Bangladesh Email: sumalchandra@yahoo.com Abstract We have investigated theoretically electronic and superconducting properties of FeSe superconductor. The calculated specific heat effective mass is found m*/mband = 2.15. This result is comparable with angle resolved photoemission spectroscopy data with the other Fe-pnictdes compound. The electron-phonon coupling constant is found λ = 0.7 for this superconductor. This result is also comparable with the results of linearized augmented plane wave method. It is observed that the FeSe system is more correlated system compared to the Fe element. Keywords: FeSe Superconductors; Specific Heat; Effective Mass; Electron-phonon coupling constant. 1. Introduction The discovery of superconductivity in Fe-pnictdes [1] compounds has attracted tremendous interest the research activity in strongly correlated electron system. PbO type-FeSe phase has been reported to be a superconductor with transition temperature, Tc ~ 8 K [2]. The newly discovered superconductor FeSe has been subjected to theoretical investigation by Subedi and his colleagues [3]. On the other hand the compounds FeSe1-xTex (x = 0- 1), where Te substitution has an effect on superconductivity, have been investigated experimentally [4-5]. It was found that Tc increases with Te doping, reaching a maximum Tc∼15 K at about 50-70% substitution, and then decreases with more Te doping. Recently a systematic theoretical examination of the effect of Te doping on electronic properties in FeSe1-xTex system are discussed in [6]. The electronic properties namely the band effective mass m*/mband is not established well for Fe and Fe-pnictdes compound such as FeSe compound. ------------------------------------------------------------------------ * Corresponding author. American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2018) Volume 41, No 1, pp 225-232 226 The experimental mass enhancement m*/mband= 1.3−2.1 is observed in LaOFeP compound within the angle resolved photoemission spectroscopy (ARPES) data [7]. It is also comparable to the ARPES data [8] as well as the experimental data obtained by the low temperature specific heat data 1.4-2.1 [9] of Fe element. However, it is significant lower that the renormalization seen in most transition metal oxides including the cuprate high temperature superconductors [10-11] and is not clear for the FeSe system. Moreover the experimental estimation mass enhancement is found m*/mband= 6-20 for FeSe0.42Te0.58 compound from the ARPES study [12]. On the hand electron-phonon coupling constant play an important role in superconductivity according to the McMillan equation. At ambient conditions, the electron-phonon coupling constant is found for FeSe compound λ= 0.603 -0.65[13] by inverting the McMillan equation using the measured Debye temperature θD. Other computational papers [3, 14] found this value of approximately λ = 0.15 using linear response theory. It is not clear the source of discrepancy between these calculations. The effective mass enhancement and electron- phonon coupling constant for FeSe compound are not established well. It is therefore important to fully study the electronic structure of FeSe compound, and its implication to specific heat effective carrier mass as well as electron-phonon coupling constant. We already performed the electronic (band and density of state (DOS)) and elastic properties (transverse vt, longitudinal vl wave velocity, and θD) calculations of FeSe1-xTex [6,15] compounds using plane wave pesoduepotential method (PWPP) within the fame work of the density functional theory [DFT] [16] implemented in CASTEP code [17]. The exchange correlation term is considered by the local density approximation (LDA). We also performed tight binding liner muffin-tin-orbitals method (TB-LMTO) with atomic sphere approximation (ASA) [18] in order to understand the basic electronic structure of FeSe compound [19]. We used these DOS and Debye temperature θD for calculating the effective mass and electron phonon coupling constant λ. In this work we studied the Sommerfeld coefficient and specific heat effective mass for FeSe compound. Finally we calculated the electron-phonon coupling constant λ for FeSe superconductor by inverting the McMillan equation using Debye temperature θD and experimental critical temperature Tc. We discussed the results of effective mass and electron-phonon coupling constant in both the PWPP and LMTO methods. 2. Results and Discussions 2.1 Electronic properties of FeSe compound In order to explain the electronic properties, we have taken the calculated value of DOS from our calculations [6, 19]. We obtained the DOS values in the PWPP and LMTO methods 1.49 and 1.81 eV −1 on a per Fe both spins basis, respectively [6, 19]. This result is comparable to the 1.41 [20] using SPRKKR code. We can evaluate the Sommerfeld constant γ from density of state. The value of N(EF) is directly proportional to the Sommerfeld constant [21] American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2018) Volume 41, No 1, pp 225-232 227 � = 1 3 ��� �� ��. 1� Here N(EF) is the density of state in the Fermi level. � is the Boltzmann constant. We can calculate the Sommerfeld constant γ by using the formula (1). We obtained the γ values 3.51 [6] and 4.26 (mJ/molK 2 ) in the PWPP and LMTO methods, respectively. On the other hand the experimental Sommerfeld constant for 9.17mJ/mol K 2 obtained low temperature specific heat [2]. The experimental estimations of γ for FeSe compound is quite lager than the calculated results. This is due to the electron-phonon interactions. The electron-phonon interactions are not included in the band structure calculations. The calculated electronic specific heat ������can be compared with the available experimental data on the electronic specific heat co- efficient ���� . The differences between ������and ���� are usually attributed to the renormalization of one- electron effective masses due to the electron-phonon interaction. Therefore, we obtain the general expression of effective mass [21] ���� ������ = �∗ ����� . 2� We can calculate the specific heat effective mass from eq. (2), and obtain m*/mband= 2.61 in PWPP method [6]. On the other hand we obtain m*/mband= 2.15 in LMTO method. A sheet dependent mass enhancement m*/mband= 1.3−2.1 is observed in 1111 Fe-pnictides of the form LaOFeP within the ARPES. The specific heat effective mass of FeSe compound is larger than the ARPES value of LaOFeP compound in PWPP method [6]. The LMTO method value is comparable with the ARPES data. Table 1: The effective mass enhancement of different elements. LaOFeP FeSe Fe m*/m 1.3-2.1[7] 2.61[6] 2.15 a 1.4-2.1[9] a This with LMTO method 2.2 Superconductivity properties of FeSe compound The DOS calculations are used to calculate superconductivity parameters. We can calculate the electron-ion matrix element known as the Hopfield parameter, η�[22], using the following formula [23] for each atom η� = � �� < � >� . (3) Where N(EF) is the total DOS per spin at EF. < � >�is the electron-ion matrix element for each atom type, which is calculated by the Gaspari and Gyorffy theory [23]. It is given by the following equation American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2018) Volume 41, No 1, pp 225-232 228 < � >�= � ���� ��"2 # + 1�%&'� ()*+� − ()���)*+� �)� �)*+ +�,��) +�,� � )./ . 4� Where �)�are the per spin angular momenta (l) components of the DOS at EF for atom type i, �) +�,� are the so- called free-scatterer DOS for atom type i, and ()� are scattering phase shifts calculated at the muffin-tin radius and at EF for atom type i. Free-scatterer DOS are calculated by �) +� = 1 � � 2# + 1�2 3�4)� 3, � 56 / � 73, 5� and scattering phase shifts are calculated by tan () <=, �� = >) ?<=� − >) ?<=�@) ?, <=� ')A ?<=� − ')A ?<=�@) ?, <=�, 6� where RS is the muffin-tin radius, >) are spherical Bessel functions, η)are spherical Newnaun functions, @) = CDE CD is the logarithmic derivative of the radial wave function, 4), One can evaluate at RS for different energies, ? = 1 �. The 4) can compute by solving the radial wave equation at each k point in the Brillouin zone. The Hopfield parameter is then used to calculate the electron-phonon coupling constant, which is obtained using the following equation from McMillan’s strong-coupling theory [24] λ = " η� F� < G� >, 7� � �.+ where Mi is the atomic mass of atom type i, and <ω 2 > is the average of the squared phonon frequency taken from the calculated Debye temperature of Chandra and his colleagues [15]. The Debye temperature is related to the phonon frequency by < G� >= + � IJ� . (8) The Debye temperature is proportional to the mean sound velocity vm [25]: IJ = � KL M N� OPQR ST U , (9) where h is Planck’s constant, the volume of unit cell of FeSe system is V0 =71.766ÅN and n is the number of atoms in unit cell, vm is the average wave velocity and can be obtained from the transverse vt and longitudinal wave velocity vl, respectively [25] American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2018) Volume 41, No 1, pp 225-232 229 ST = W+N X � YZU + + YDU [\ + N] . (10) We have taken vt and vl values from [14], and obtained the values of vm as well as IJ. The estimated values of Debye temperature is 285 K [15]. We used this value for calculation of electron-phonon coupling constant. The McMillan equation for critical superconductivity transition temperature is given by the [24] _̂ = `a +.Ob exp fg+./O +*λ� λgµ∗ +*λ� h. (11) We have also taken the experimental Tc = 8K for FeSe from Ref. [3]. The µ∗is the Coulomb pseudopotential, given by the Bennemann-Garland equation [26] µ∗ = 0.13 j kl� +*j kl� . (12) From the above expressions (7) and (11) we can determine the coupling constant λ into two ways. One method is first find the value Hopfield parameter ηi, and phonon frequency by conversion of Debey temperatures. The last one is from direct conversion of superconducting critical temperature Tc and Debey temperature. In our calculation we have chosen the last one for determine the coupling constant. Our calculated Debye temperature is 285 K [15]. The critical transition temperature Tc is given 8 K for FeSe superconductor [3]. The calculated electron phonon- coupling constant is 0.7. This is comparable to the pressure dependent value 0.61-0.71 [27]. At ambient conditions, our electron-phonon coupling constant calculation λ is consistent with the value calculated by Ksenofontov and his colleagues [13], λ = 0.65, by inverting the McMillan equation using the measured θD. Other calculations based on the linear response theory [3, 14] find values of approximately λ = 0.15. It is not clear to the source of discrepancy between the band calculations and linear response theory. Other iron-based superconducting material, LaFeAsO, electron-phonon coupling constants of approximately 0.2 [28, 29] have been reported. The LaFeAsO element is known to be on the verge of magnetic instability. This suggests that spin fluctuations are important in this material. Table 2: The electron-phonon coupling constant of different elements. FeSe LaFeAsO λλλλ 0.7 a 0.61-0.71[27] 0.65[13] 0.15[3] 0.2 [28, 29] American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2018) Volume 41, No 1, pp 225-232 230 3. Conclusions In the present work we investigated electronic and superconductivity properties theoretically based on the first- principles DFT total energy calculations in PWPP and LMTO methods. The specific heat effective mass enhancement factor of FeSe compound is found to be 2.61 in FeSe compound in the PWPP and 2.15 in the LMTO methods, respectively. It is comparable to a sheet dependent mass enhancement m*/mband= 1.3−2.1 is observed in Fe- pnictides of the form LaOFeP. The effective mass of FeSe compound is larger than in the Fe compound. It indicates that FeSe compound is more correlated compare to the Fe element. The calculated electron-phonon coupling constant is 0.7. This is comparable to the pressure dependent value 0.61-0.71. At ambient conditions, our electron-phonon coupling constant calculation λ = 0.603 is consistent with the value calculated by Ksenofontov and his colleagues λ = 0.65, by inverting the McMillan equation using the measured θD. Other calculations based on the linear response theory find values of approximately λ = 0.15. It is not clear the source of discrepancy between the band calculations and the linear response theory. 4. Recommendations We have performed these calculations in the non-magnetic cases for FeSe compound. The calculations using the ferromagnetic, antiferromagnetic, and other magnetic orders should be considered for further study, but are beyond the scope of this paper. References [1] Y. Kamihara,T. Watanabe, M. Hirano, and H. Hosono (2008) Iron-Based Layered Superconductor La[O1-xFx]FeAs (x = 0.05−0.12) with Tc = 26 K, J. Amer. Chem. Soc. 130, 3296. [2] F. C. Hsu, J. Y. Luo, K.W. Yeh, T. K. Chen, T. W. Huang, P. M. Wu, Y.C. Lee, Y. L. Huang, Y.Y. Chu, D.C. Yan, and M. K. Wu (2008) Superconductivity in the PbO-type structure α-FeSe, Proc. Natl. Acad. Sci. U.S.A. 38, 14262. [3] A. Subedi, L. Zhang, D. J. Singh, and M. H. Du (2008) Density functional study of FeS, FeSe, and FeTe: Electronic structure, magnetism, phonons, and superconductivity, Phys. Rev. B 78, 134514. [4] Y. Mizuguchi, F. Tomioka, S. Tsuda, T. Yamaguchi, and Y. Takano (2009) Substitution Effects on FeSe Superconductor, J. Phys. Soc. Jpn. 78, 074712. [5] K. W. Yeh, H. C. Hsu, T. W. Huang, P. M. Wu, Y. L. Huang, T. K. Chen, J. Y. Luo, and M. K. Wu (2008) Se and Te Doping Study of the FeSe Superconductors, J. Phys. Soc. Jpn. Suppl. C 77, 19. [6] Sumal Chandra (2018) Electronic Properties of FeSe1-xTex Superconductors, J. Sci. Res. Phys. Math. Sci. Vol. 5(1) ISSN:2349-7149. [7] D. H. Lu, M. Yi, S. -K. Mo, S. A. Erickson, J. Analytis, J, -H. Chu, D. J. Sing, Z. Hussain, T. H. American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2018) Volume 41, No 1, pp 225-232 231 Geballe, I. R. Fisher, and Z.-X. Shen (2008) Electronic structure of the iron-based superconductor LaOFeP, Nature 455, 81. [8] J. Sànchez- Barriga, J. Fink, V. Boni, I. Di Marco, J. Braun, J. Minàr, A. Varykhalov, O. Rader, V. Bellini, F. Manghi, M.I. Katsnelson, A. I. Lichtenstein, O. Erikson, W. Eberhardt, and H. A. Dürr (2009) Strength of Correlation Effects in the Electronic Structure of Iron, Phys. Rev. Lett. 103, 267203. [9] C. H. Cheng, C. T. Wei, and P. A. Beck (1960) Low-temperature specific heat of body-centered cubic alloys of 3d transition elements, Phys. Rev. 120, 426. [10] M. M. Qazilbash, J. J. Hamlin, R. E. Baumbach, L. Zhang, D. J. Sing, M. B. Maple, and D. N. Basov (2009) Electron correlations in the Iron-pnictides, Nature Physics 5, 647. [11] X. Y. Zhou, T. Yoshida, A. Lanzara, P. V. Bogdanov, S. A. Kellar, K. M.shen, W. L. Yang, F. Ronning, T. Sasagawa, T. Kakeshita, T. Noda, H. Eisaki, S. Uchida, C. T. Lin, F. Zhou, J.W. Xiong, W. X. Ti, Z. X. Zhao, A. Fujimori, Z. Hussain, and Z. –X. Shen (2003) High-temperature superconductors: Universal nodal Fermi velocity, Nature 423, 398. [12] A. Tamai, A. Y. Ganin, E. Rozbicki, J. Bacsa, W. Meevasana, P. D. C. King, M. Caffio, R. Schaub, S. Margadonna, K. Prassides, M. J. Rosseinsky, and F. Baumberger, (2010) Strong electron correlations in the normal state of the iron based FeSe0.42Te0.58 superconductor observed by angle- resolved photomission spectroscopy, Phys. Rev. Lett. 104, 097002. [13] V. Ksenofontov, G. Wortmann, A. I. Chumakov, T. Gasi,S. Medvedev, T. M. McQueen, R. J. Cava, and C. Felser (2010) Density of phonon states in superconducting FeSe as a function of temperature and pressure, Phys. Rev. B 81, 184510. [14] TimurBazhirov and Marvin L. Cohen (2012) Spin resolved electron phonon coupling in FeSe and KFe2Se2, Phys. Rev. B 86, 134517. [15] S. Chandra, and A. K. M. A. Islam (2010) Elastic properties of mono- and polycrystalline PbO-type FeSe1-xTex (x=0-1.0): A first principles study, Phycica C 470, 2072 [16] W. Kohn and L. J. Sham (1965) Self-Consistent Equations Including Exchange and Correlation Effects,Phys. Rev. A 140, 1133-1138. [17] S.J. Clark, M.D. Segall, C.J. Pickard, P.J. Hasnip, M.J. Probert, K. Refson, M.C.Payne (2005) First principles methods using CASTEP, Z. Kristallogr. 220 (5-6) 567-570. [18] O. Jepsen and O. K. Andersen (2000) The Stuttgart TB-LMTO Program (Federal Republic of Germany). [19] Sumal Chandra, (2013) Construction of tight-binding real space Hamiltonian and its application to American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) (2018) Volume 41, No 1, pp 225-232 232 correlated electron system, Master’s thesis, University of the Ryukyus, Japan. [20] J. Kumar, S. Auluck, P. K. Ahluwalia, and V. P. S. Awana (2012) Chalcogen height dependence of magnetism and Fermiology in FeTexSe1−x, Supercond. Sci. Technol. 25, 095002. [21] C. P. Poole, H. A. Farach, R. J. Creswick, and R. Prozorov, Superconductivity 2nd ed., Elsevier, Netherlands, 2007, p. 12. [22] J. Hopfield (1969)Angular momentum and transition metal superconductivity, Phys. Rev. 186, 443. [23] G. Gaspari and B. Gyorffy (1972) Electron-Phonon interactions, d Resonances, and Superconductivity in transition metals, Phys. Rev. Lett. 28, 801. [24] W. McMillan (1968) Transition temperature and strong-couple superconductor, Phys. Rev. 167, 331. [25] O. L. Anderson (1963) A simplified method for calculating the Debye temperature from elastic constant, J. Phys. Chem. Solids. 24, 909. [26] K. H. Bennemann and J. W. Garland (1972) Theory for superconductivity in d band metal, AIP Conf. Proc. 4, 103. [27] A. P. Koufos, D. A. Papaconstantopoulos, and M. J. Mehl (2014) First-principles study of the electronic structure of iron-selenium: Implications for electron-phonon superconductivity Phys. Rev. B 89, 035150. [28] L. Boeri, O. V. Dolgov, and A. A. Golubov (2008) Is LaFeAsO1-xFx an electron-phonon Superconductor?, Phys. Rev. Lett.101, 026403. [29] I. I. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du (2008) Unconventional superconductivity with a sign reversal in the order parameter of LaFeAsO1-xFx , Phys. Rev. Lett. 101, 057003.