潘 妍 ，于海林 ，b，錢 斌 ，b，江學(xué)范 ，b，郭 磊 ，呂 群 ，谷全超

（a.常熟理工學(xué)院 物理與電子工程學(xué)院；b.江蘇省新型功能材料重點建設(shè)實驗室，江蘇 常熟 215500）

自從發(fā)現(xiàn)鐵磷族化合物和鐵硫族化合物的超導(dǎo)電性和反鐵磁序有關(guān)以來[1-8]，人們認(rèn)為超導(dǎo)配對機(jī)制與自旋漲落密切相關(guān)[9-13]，超導(dǎo)電性和自旋漲落有關(guān)在超導(dǎo)態(tài)的自旋共振方面得到了進(jìn)一步的證實[11-15]. 為了闡明自旋漲落在超導(dǎo)配對機(jī)制上的本質(zhì)作用，在未摻雜的母體化合物的磁性，以及這類材料的磁性和超導(dǎo)電性相互作用方面做了許多研究工作.雖然鐵磷族化合物和鐵硫族化合物有相似的層狀結(jié)構(gòu)和費米面的拓?fù)湫?，它們未摻雜母體化合物的反鐵磁序（AFM）是完全不同的.例如，LaOFeAs[16]和BaFe2As2[17]表現(xiàn)出“單條紋”（single-stripe）反鐵磁序，而鐵硫族元素化合物中的Fe1+yTe[18]顯示出“雙條紋”（double-stripe）形狀的反鐵磁序.

這兩類材料在反鐵磁序結(jié)構(gòu)上的顯著差別，引發(fā)了關(guān)于磁性起源方面的大量爭論.對于鐵磷族化合物超導(dǎo)配對機(jī)制的解釋并不適用于Fe1+yTe.有趣的是除了費米面驅(qū)動的磁性機(jī)制外，多種基于局域磁矩超交換相互作用的模型也相繼被提出來[19-24].最近，考慮局域磁矩和巡游電子耦合的模型也被提了出來[19,22,25-26].此外，Yin等提出了一個統(tǒng)一的微觀模型來解釋各種反鐵磁關(guān)聯(lián)現(xiàn)象[27].該模型認(rèn)為局域磁矩和巡游電子共存，局域自旋的反鐵磁超交換耦合與雙交換鐵磁相互作用競爭，“single-stripe”和“double-stripe”反鐵磁序都可以用陰離子距Fe平面的高度來解釋.

近來，Moon和Choi利用第一性原理計算，發(fā)現(xiàn)陰離子高度對這類材料基態(tài)磁結(jié)構(gòu)有重要影響[28].他們的工作顯示如果陰離子Te的高度ZTe小于臨界值（～1.71 ?），F(xiàn)e1+yTe的反鐵磁序?qū)摹皊ingle-stripe”向“dou?ble-stripe”轉(zhuǎn)變，這和實驗結(jié)果完全一致.同時還發(fā)現(xiàn)，從總能的角度看，隨著Zanion的增加，有向鐵磁態(tài)轉(zhuǎn)變的趨勢.最近，錢斌等在實驗上首次合成了具有與鐵磷族化合物類似結(jié)構(gòu)的新型鐵磷族化合物CuFeSb[29]，并發(fā)現(xiàn)其具有完全不同于目前研究中的鐵磷族化合物和鐵硫族化合物的金屬性反鐵磁序，而是表現(xiàn)出金屬性鐵磁序.然而，目前為止，還沒有針對CuFeSb的磁性和電子特性的理論研究.本文采用第一性原理對CuFeSb的磁性和電子特性進(jìn)行計算，闡明其金屬性鐵磁特性產(chǎn)生的物理機(jī)制.

1 計算方法

密度泛函（DFT）下的廣義梯度近似（GGA）是目前較為準(zhǔn)確的電子結(jié)構(gòu)計算的理論方法[30-31]，它是電子密度及其梯度的泛函[31]，并把密度分布的空間變化包括在其中，實現(xiàn)了可較大幅度減少局域密度近似（LDA）誤差的目的[32].離子和電子間的相互作用采用投影綴加平面波（projector augmented wave，PAW）[33]方法來處理.這種方法可以相當(dāng)程度地減少過渡金屬或第一行元素的原子所必需的平面波的數(shù)量，并且對于磁性化合物而言，PAW方法要優(yōu)于超軟贗勢方法.本文采用的計算軟件為Vienna ab-initio Simulation package（VASP），截斷能選取500 eV，K點采樣選擇Monkhorst-Pack方法.對鐵磁相的計算，采用的是CuFeSb的四方單包，其K點采樣使用8×8×5，而對反鐵磁相的計算采用的2×2×1的超包，其K點采樣使用4×4×5.

2 結(jié)果與討論

2.1 結(jié)構(gòu)特性

CuFeSb的四方結(jié)構(gòu)屬于P4/nmm空間群，如圖1（a）所示.在對原子結(jié)構(gòu)進(jìn)行弛豫優(yōu)化時，我們選擇的晶格常數(shù)為實驗值a=3.934 ?，c=6.251 ?[29]，在本工作中，我們只對原子的內(nèi)部坐標(biāo)進(jìn)行了弛豫，計算結(jié)果如表1所示.從表1中我們發(fā)現(xiàn)計算結(jié)果和實驗值[29]能夠很好地吻合.弛豫以后的Sb原子距Fe原子的距離為1.839 ?，遠(yuǎn)遠(yuǎn)高于Moon和Choi提出的臨界值[28].這是CuFeSb表現(xiàn)出鐵磁序的一個原因，這一點和Moon等通過計算得出的趨勢是一致的.

2.2 磁性和電子特性

對于CuFeSb的磁性和電子特性的計算，其鐵磁相采用的是如圖1（a）所示的四方相單包.而對反鐵磁序的計算，采用的2×2×1的超包，其Fe原子平面磁矩的排列，我們考慮了三種構(gòu)型，即C型、E型和G型，如圖1（b）-（d）所示. 我們對弛豫以后的總能量進(jìn)行比較發(fā)現(xiàn)，鐵磁相具有最低的總能，即CuFeSb的基態(tài)磁結(jié)構(gòu)為鐵磁相.對其磁性的計算我們得到其四方相單包的磁矩為 4.50 μB，這和前期的實驗結(jié)果[29]能夠很好地吻合.

我們對CuFeSb的能帶和態(tài)密度進(jìn)行了詳細(xì)計算，如圖2和圖3所示.從圖2中我們可以看出CuFeSb表現(xiàn)出典型的金屬性，這和實驗結(jié)果完全一致.自旋向上的能帶結(jié)構(gòu)和自旋向下的能帶結(jié)構(gòu)總體上相近，如圖2（右）和（左）所示.其主要差別在于在-4 eV～0.0 eV范圍內(nèi)，自旋向下的能帶和自旋向上的能帶相比，其色散性更明顯，這一點主要由Fe原子d軌道上不同自旋方向的電子的占據(jù)情況決定.如圖3所示，相比于Cu原子，F(xiàn)e原子的3d軌道其局域性減小，巡游性增強(qiáng).其自旋向下的電子態(tài)主要由反鍵軌道占據(jù)，這樣導(dǎo)致其自旋向下的能帶較自旋向上的能帶有更大的色散性.從Fe的平面內(nèi)軌道dx2和dxy來看，其自旋向下的電子態(tài)基本上完全由反鍵軌道占據(jù)，這表明平面內(nèi)Fe原子間的相互作用是比較強(qiáng)的.考慮到Sb原子距Fe原子平面的高度較大，導(dǎo)致Fe原子和Sb原子間的相互作用減弱，從而相應(yīng)導(dǎo)致平面內(nèi)Fe原子的相互作用增強(qiáng).

圖1 (a)CuFeSb四方相單包結(jié)構(gòu)示意圖,黑色表示Fe原子,深灰色表示銅原子,灰色表示Sb原子.(b)、(c)、(d)分別表示Fe原子平面原子磁矩的排列,白色和黑色分別表示磁極化方向相反的兩類Fe原子.(b)表示C型(collinear),(c)表示E型(bicollinear),(d)表示G型(checkerboard)

表1 弛豫以后的CuFeSbC的結(jié)構(gòu)參數(shù)

圖2 鐵磁相CuFeSb的能帶和相應(yīng)的態(tài)密度.(左)自旋向下的能帶,(中)態(tài)密度圖,紅色曲線表示自旋向下的態(tài)密度,黑色曲線表示自旋向上的態(tài)密度,(右)自旋向上的能帶

3 總結(jié)

利用第一性原理計算對CuFeSb的晶體結(jié)構(gòu)、磁性和電子特性進(jìn)行了研究.從能帶和態(tài)密度的計算結(jié)果來看，CuFeSb表現(xiàn)出典型的金屬性，這一點和實驗結(jié)果完全一致.通過對Fe原子的d軌道進(jìn)行深入分析，發(fā)現(xiàn)其d軌道的局域性相對減弱，巡游性增強(qiáng)，并發(fā)現(xiàn)其平面內(nèi)Fe原子的相互作用較強(qiáng)，其金屬性鐵磁序與平面內(nèi)Fe原子的相互作用增強(qiáng)有關(guān).

圖3 Fe(左)和Cu(右)原子的d軌道態(tài)密度圖

[1]Zhao J,Huang Q,Cruz C de la,et al.Structural and magnetic phase diagram of CeFeAsO1-xFxand its relation to high-temperature superconductivity[J].Nat Mater,2008,7:953.

[2]Luetkens H,Klauss H H,Kraken M,et al.The electronic phase diagram of the LaO1-xFxFeAs superconductor[J].Nat Mater,2009,8:305.

[3]Drew A J,Niedermayer C,Baker P J,et al.Coexistence of static magnetism and superconductivity in SmFeAsO1-xFx as revealed by muon spin rotation[J].Nat Mater,2009,8:310.

[4]Chen H,Ren Y,Qiu Y,et al.Coexistence of the spin-density wave and superconductivity in Ba1?xKxFe2As2[J].Europhys Lett,2009,85:17006.

[5]Chu J H,Analytis J G,Kucharczyk C,et al.Determination of the phase diagram of the electron-doped superconductor Ba(Fe1?xCox)2As[J].Phys Rev B,2009,79:014506.

[6]Nandi S,Kim M G,Kreyssig A,et al.Anomalous Suppression of the Orthorhombic Lattice Distortion in Superconducting Ba(Fe1-xCox)2As2Single Crystals[J].Phys Rev Lett,2010,104:057006.

[7]Liu T J,Hu J,Qian B,et al.From(π,0)magnetic order to superconductivity with(π,π)magnetic resonance in Fe1.02Te1?xSex[J].Nat Mater,2010,9:718.

[8]Katayama N,Ji S,Louca D,et al.Investigation of the Spin-Glass Regime between the Antiferromagnetic and Superconducting Phases in Fe1+ySexTe1-x[J].J Phys Soc Jpn,2010,79:113702.

[9]Ahilan K,Ning F L,Imai T,et al.19F NMR investigation of the iron pnictide superconductor LaFeAsO0.89F0.11[J].Phys Rev B,2008,78:100501.

[10]Nakai Y,Ishida K,Kamihara Y,et al.Evolution from Itinerant Antiferromagnet to Unconventional Superconductor with Fluorine Doping in LaFeAs(O1-xFx)Revealed by 75As and 139La Nuclear Magnetic Resonance[J].J Phys Soc Jpn,2008,77:073701.

[11]Christianson A D,Goremychkin E A,Osborn R,et al.Unconventional superconductivity in Ba0.6K0.4Fe2As2from inelastic neutron scattering[J].Nature,2008,456:930.

[12]Lumsden M D,Christianson A D,Parshall D,et al.Two-dimensional resonant magnetic excitation in BaFe1.84Co0.16As2[J].Phys Rev Lett,2009,102:107005.

[13]Chi S,Schneidewind A,Zhao J,et al.Inelastic Neutron-Scattering Measurements of a Three-Dimensional Spin Resonance in the FeAs-Based BaFe1.9Ni0.1As2Superconductor[J].Phys Rev Lett,2009,102:107006.

[14]Mook H A,Lumsden M D,Christianson A D,et al.Neutron scattering patterns show Superconductivity in FeTe0.5Se0.5likely results from itinerant electron fluctuations[J].e-print arXiv,2009,0904:2178.

[15]Qiu Y,Bao W,Zhao Y,et al.Spin Gap and Resonance at the Nesting Wave Vector in Superconducting FeSe0.4Te0.6[J].Phys Rev Lett,2009,103:067008.

[16]Cruz C de la,Huang Q,Lynn J W,et al.Magnetic order close to superconductivity in the iron-based layered LaO1-xFxFeAs systems[J].Nature,2008,453:899.

[17]Huang Q,Qiu Y,Bao W,et al.Neutron-Diffraction Measurements of Magnetic Order and a Structural Transition in the Parent BaFe2As2Compound of FeAs-Based High-Temperature Superconductor[J].Phys Rev Lett,2008,101:257003.

[18]Bao W,Qiu Y,Huang Q,et al.Tunable(δπ,δπ)-Type Antiferromagnetic Order in α-Fe(Te,Se)Superconductors[J].Phys Rev Lett,2009,102:247001.

[19]Ma F,Lu Z-Y,Xiang T.Arsenic-bridged antiferromagnetic superexchange interactions in LaFeAsO[J].Phys Rev B,2008,78:224517.

[20]Si Q,Abrahams E.Strong Correlations and Magnetic Frustration in the High Tc Iron Pnictides[J].Phys Rev Lett,2008,101:076401.

[21]Yildirim T.Origin of the 150-K Anomaly in LaFeAsO:Competing Antiferromagnetic Interactions,Frustration,and a Structural Phase Transition[J].Phys Rev Lett,2008,101:057010.

[22]Ma F,Ji W,Hu J,et al.First-Principles Calculations of the Electronic Structure of Tetragonal α-FeTe and α-FeSe Crystals:Evi-dence for a Bicollinear Antiferromagnetic Order[J].Phys Rev Lett,2009,102:177003.

[23]Hu J,Xu B,Liu W,et al.Unified minimum effective model of magnetic properties of iron-based superconductors[J].Phys Rev B,2012,85:144403.

[24]Hu J,Ding H.Local antiferromagnetic exchange and collaborative Fermi surface as key ingredients of high temperature superconductors[J].e-print arXiv,2011,1107:1334.

[25]Wu J,Phillips P,Castro Neto A H C.Theory of the Magnetic Moment in Iron Pnictides[J].Phys Rev Lett,2008,101:126401.

[26]Johannes M D,Mazin I I.Microscopic origin of magnetism and magnetic interactions in ferropnictides[J].Phys Rev B,2009,79:220510.

[27]Yin W-G,Lee C-C,Ku W.Unified Picture for Magnetic Correlations in Iron-Based Superconductors[J].Phys Rev Lett,2010,105:107004.

[28]Moon C-Y,Choi H J.Chalcogen-Height Dependent Magnetic Interactions and Magnetic Order Switching in FeSexTe1-x[J].Phys Rev Lett,2010,104:057003.

[29]Qian B,Lee J,Hu J,et al.Ferromagnetism in CuFeSb:Evidence of competing magnetic interactions in iron-based superconductors[J].Phys Rev B,2012,85:144427.

[30]Marlo M.Density-functional study of bulk and surface properties of titanium nitride using different exchange-correlation functionals[J].Phys Rev B,2000,62:2899.

[31]李震宇，賀偉，楊金龍.密度泛函理論及其數(shù)值方法新進(jìn)展[J].化學(xué)進(jìn)展，2005，17（2）：192-202.

[32]黃美純.密度泛函理論的若干進(jìn)展[J].物理學(xué)進(jìn)展，2000，20（3）：199-219.

[33]Bl?chl P E.Projector augmented-wave method[J].Phys Rev B,1994,50:17953.