第46卷第2期 2017年3月 内蒙古师范大学学报(自然科学汉文版) Journal of Inner Mongolia Normal University(Natural Science Edition) Vo1.46 No.2 Mar.2017 Fe14 T2(T—Cr,Mn,Co,Ni)二元合金结构 与力学性质的第一性原理研究 马 爽 ,董 雪 ,武晓霞 ,管鹏飞。,那日苏 (1.内蒙古师范大学物理与电子信息学院,内蒙古呼和浩特010022; 2.内蒙古科技大学数理与生物工程学院,内蒙古包头014010;3.北京计算科学研究中心,北京100193) 摘要:运用以密度泛函理论为基础的投影缀加波方法,研究Fe Tz(T=Cr,Mn,Co,Ni)二元合金不同相 结构的总能、磁性和力学性质.结果表明,Fe, T。(T—Cr,Mn,Co,Ni)Z.元合金的体心立方相(bcc相)的铁磁态能 量最低,应为体系的基态.铁磁性的面心立方相(fcc相)具有高自旋和低自旋的特性,而反铁磁性的fce相具有更 低的能量,相对铁磁态更为稳定.弹性常数的计算表明,铁磁性fee相的HS态和铁磁Fe Ni fee相LS态在力学 上不稳定,其他相均在力学上稳定. 关键词:Fe T (T Cr,Mn,Co,Ni)二元合金;相稳定;力学性质;第一性原理计算 中图分类号:0 41 文献标志码:A 文章编号:1001-8735(2017)02-0182-05 doi:10.3969/j.issn.1001-8735.2017.02.007 二元铁合金具有复杂的相图和丰富的物理特性,主要原因在于铁合金的结构与磁有序性有强烈的关 联[1 ].不同结构的Fe及其合金可通过衬底约束和间隙原子掺杂等方法得到[5 ],不同的结构产生不同的 磁性及功能,其中包括早年的Invar合金[4 和永磁钢『5 等,以及可电场调控磁有序性的Fe金属纳米结构[6]. 因此,了解不同结构Fe合金的组分、原胞体积以及磁有序性和相稳定性之间关联特性,对功能合金的设计 及制备有一定的参考价值.到目前为止,就二元Fe合金的不同相结构及其磁有序性的研究较多[7 引,但大 部分仅关注某一种合金或者某一种磁性状态.Okatov等[7]研究了不同磁有序状态下金属Fe由bcc到fcc的 结构相变过程,发现fee结构的磁性为非共线的反铁磁(antKerromagnetic,AFM)态,而具有共线的反铁磁 序的Fe在c/a 1.1的面心四方(face centered tetragonal,let)结构中具有最低的能量.然而他们没有分析 元胞体积对磁有序性的影响,而fcc相的铁磁(ferromagnetic,FM)态具有高自旋(High Spin,HS)和低自旋 (Low Spin,LS)的特性[8].Wang等。。 对fee结构Fe7oTm2(Tm—Cr,Ni,Ta,Zr)合金体系的研究发现,体系 中产生Y.3<lll>晶界时体系能量较低,同时具有奇异的磁性.Gebhardt等[1o]的研究结果表明,高Mn含量 的fee相FeMn合金在室温下呈现AFM态.Mirzoev等口1]发现,在bcc相FM态的Fe中固溶Mn原子组分 >3 时,近邻原子的耦合由AFM转变为近似FM耦合.Filiberto等 发现,Co组分较小的FeCo合金的 fee结构为铁磁性,且保留了Fe的HS和LS特性,而fee相AFM态是否稳定尚无定论.Benjamin等 纠发 现,bcc相的FeNi合金是FM态,从二元相图Ira]可知,Ni在Fe的bcc相固溶度,此外对于低温时fcc相的 FeNi合金的结构和磁性了解较少.Shi等预测,fee相的Fe1-xPd 合金[14]和Fe1-xCu 合金ll5]可能以亚稳态 的形式存在.综上,Fe及Fe 一 T 合金的不同结构的磁性和稳定性有待系统研究.最近董雪等口 研究了不 同组分Fe 一 T 合金的能量和结构稳定性,但尚缺乏铁磁性fee相的结构和力学性能方面的结果.因此,本 文用密度泛函理论(DFT)研究Fe和Fe1 T。(T:Cr,Mn,Co,Ni)二元合金的fee相FM态和AFM态以及 bee相FM态的能量、磁性和力学性质. 收稿日期:2016-06—02 基金项目:国家自然科学基金资助项目(11464037);内蒙古自然科学基金资助项目(2016MS0113) 作者简介:马爽(1992一),女,内蒙古通辽市人,内蒙古师范大学硕士研究生 通信作者:那日苏(1981一),男,内蒙古通辽市人,内蒙古师范大学副教授,博士,主要从事凝聚态理论研究 第2期 马 爽等:Fe T2(T—Cr,Mn,Co,Ni)二元合金结构与力学性质的第一性原理研究 1 模型与计算 1.1计算模型和力学性质计算 对于bcc相的Fe T (T—Cr,Mn,Co,Ni)二元合金,采用2×2×2的16个原子的超胞模型,基矢大小 n一6一C,掺杂的T(T—Cr,Mn,Co,Ni)原子占据Fe原子的(0,0,0)和(0.5,0.5,0.5)位置.对于fee相 Fe T。(T—Cr,Mn,Co,Ni)二元合金,采用2×2×2的bct超胞模型,基矢a=b≠c,其中c/a约为1.414,根 据不同的磁性结构对c/a进行充分弛豫. 为了了解Fe T (T—Cr,Mn,Co,Ni)二元合金力学性质及其力学稳定性,对弹性常数进行了计算口 ]. 对于立方晶系的Fe T (T—Cr,Mn,Co,Ni)二元合金,有C 、C 和C 三个独立的弹性常数,其弹性自由能 能量在任意形变下的变化为 知E一寺cl1(P;+8l+P;)+C】2(8l e2+e2P3+P1 83)+ 1 C44(P:+P;+Pi)+O(e。), 其中C 、C 可组合成为体模量B一(C +2C )/3和四方剪切模量C 一(C ~C。 )/z.体模量可通各向同 性的体积膨胀并用物态方程拟合得到.四方剪切模量C 和C 可通过保持体积不变的正交形变D。和单斜形 变D 得到,分别为 D。+ —1 0 fH o 1一 0 ] 0 I, D + —l f 1 o ] 0 l, 【0 o( )和△E( )一2VC 1.2数值计算 0 [1~ ] J 【0 0 [1一 ] J 其中 和5 分别为正应变和剪切应变的大小.这两种形变下晶体的弹性能的变化分别为△E( )一2VC :+ +0( ),其中O( )为略去项.最后,求出的C 、C 。和C 需要满足立方晶系 的力学稳定性条件c >0,C >l c z I,c 1+2c 2>0. 使用密度泛函理论(Density Functional Theory)DT]框架下的投影缀加波方法(Projector Augmented wave method) ¨ 和VASP程序包口 ,计算合金体系的总能量以及弹性常数.选用PAW势描述价电子和离 子实之间的相互作用,采用广义梯度近似(Generalized Gradient Approximation)为基础的PBE(Perdew— Burke—Ernzerhof)交换关联势形式进行计算口 .为了确保计算过程中的收敛性,对于bee相的二元合金,平 面波基函数截断能设为500 eV,电子弛豫的自洽收敛标准设为10一eV/atom,原子之间的相互作用力收敛 判据设为0.1 eV/nm,布里渊区积分采用Monkhorst—Pack的特殊撒点方式[2¨,k点网格设为12×12×12. 对于fcc相的二元合金,截断能设为530 eV,k点设为12×12×8,其余均与bcc相的设置相同.考虑到计算过 程中FM态与AFM态自旋极化的相互作用,分别设置了不同的磁性状态,进而对不同体积下的Fe T。 (T—Cr,Mn,Co,Ni)二元合金充分弛豫,得到优化后的结构,再应用Murnaghan方程拟合物态方程 ,得 到体系的体模量B和平衡体积.弹性常数的形变应变量 均取0.00~土0.04,并通过E~ 关系拟合得到弹 性常数c 和C . 2 计算结果与讨论 2.1 Fe T (T===Cr,Mn,Co。Ni)二元合金能量与磁矩 图1是Fe和Fe T (T—Cr,Mn,Co,Ni)二元合金fcc相FM态、AFM态和bcc相FM态中每个原子的 能量、磁矩随体积变化的关系图.由图1可以看出,Fe和Fe。 T (T—Cr,Mn,Co,Ni)二元合金bcc相FM态 到fcc相AFM态体积的减小,同时由表1的数据可知,体积变化率在5 左右.体积缩小的原因主要在于fee 结构是密排结构,在其平衡态附近原子间距较小而产生反铁磁耦合,且该磁性态呈亚稳性. 在能量方面,对于Fe和Fe T (T:Cr,Mn,Co,Ni)二元合金而言,bce相FM态相对于fee相FM态和 AFM态稳定;对于fee相的FM态,都发现存在HS和LS的特性,即合金的磁矩随着体积的增加存在一个 突变,突变前后(LS和HS)在能量上均出现一个亚稳态.这与Wang等[9]发现fcc相的FM态Fe,以及 Filiberto等[8 发现fee相的FM态FeCo合金具有HS和LS特性相一致;对于fcc相的AFM态,Okatov[ ] ・184・ 内蒙古师范大学学报(自然科学汉文版) 第46卷 和Gebhardt等 。 发现,Fe、FeMn合金fcc相的磁性是AFM态,这也与图1中Fe、Fel4 Mn2相符合,即fcc 相AFM态要比FM态稳定,而对于fcc相AFM态的Fe T (T=Cr,Co,Ni)合金,从图1来看,Fe Cr 、 Fe14Mnz、Fe14Co2的fcc相AFM态要比FM态稳定.然而,fcc相FM态和AFM态,bcc相FM态的Fe和 Fe, Tz(T:==Cr,Mn,Co,Ni)二元合金的稳定性,还需要结合力学性质进行判别,即合金的力学稳定性.从图1 的右列给出的磁矩随体积变化的规律可以看出,fcc相AFM态Fe和Fe Cr 合金的磁矩基本为零,而 Fe T。(T—Mn,Co,Ni)合金的总磁矩不为零是因为部分局域磁矩未抵消,可能是这三种合金构成fcc结构 时AFM态并非严格的共线反铁磁态.fcc相FM态合金的磁矩均有一个跃变,跃变前后磁矩的差别大于 1 .值得注意的是,无论是bee还是fee结构,Cr的磁矩永远与Fe的磁矩相反,导致总磁矩小于其他合金. : -3.O 2.4 ∞7.9 嘉-刚8.0 .1.8 盖1.2 O.6 O.O 8.2 8-3 .: 一2.O i 6 ∞8.0 砉- 8.21 .1.2 量0.8 0.4 0.0 8.3 堪.4 > ∞ 毒 文 2.5 2.O 》 ∞1.5 裔 盖1.0 0.5 0.0 2.4 2.O > ∞1.6 宣1.2 O.8 O.4 V/10 3nm3 F/10 3nm3 图1不同磁有序性fec和bee相Fe和Fe T (T—Cr,Mn,Co,Ni)二元合金能量、磁矩随体积的变化 Fig.1 Volume dependence of energy and magnetic moment of Fe and Fe14Tz(T—Cr,Mn,Co,Ni)alloys with different magnetic state in fee structure and in bcc structure 2.2 Fel T2(T—Cr。Mn.Co,Ni)二元合金的力学稳定性 表1给出计算得到的fee相FM态、AFM态和bcc相FM态的Fe和Fe Tz(T=Cr,Mn。Co,Ni)--元合 金每个原子的各项参数,并与现有的不同方法计算的结果作了对比.可以看出,表1所列的计算结果与其他 方法和其他文献报道的结果较为接近,保证了计算的准确性,进而说明后续的弹性常数计算结果是可靠的. 对于Fe和Fe14T2(T—Cr,Mn,Co,Ni)二元合金,除了fee相FM HS态和Fe14 Ni2fee相FM LS态不满足立 第2期 马 爽等:Fe T。(T—Cr,Mn,Co,Ni)二元合金结构与力学性质的第一性原理研究 ・185・ 方晶系的稳定性条件外,其他相均符合稳定性条件,即在力学上是稳定的.另外,力学稳定的fcc相Fe和 Fe T。(T—Cr,Mn,Co,Ni)二元合金的弹性常数C 和C 均大于相应的bcc相,说明fcc相力学性能要比 bcc相更硬.这与含碳的钨钢、铬钢和钴钢等合金材料,从fcc顺磁相淬火时部分或全部转变为bct铁磁性的 马氏体相,且该相的力学和磁学性能都很硬等实验结果相一致_5J. 表1 fec相FM态、AFM态。以及bcc相FM态的Fe和Fe14T2(T—Cr。Mn。Co。Ni}二元合金的参数 Tab.1 Various parameters of the ferromagnetic and antiferromagnetic state of face centered cubic phase- and ferromagnetic state of body centered cubic phase of Fe and Fel4Tz(T—CrIMn-Co-Ni)alloys a/nm fee—FM—LS (PAW)[9] 0.35 一 VllO—snm3 10.53 lO.48 B/GPa 199.9 166.1 E/eV —8.08 —8.16 / B 1.O7 0.98 C /GPa 95.6 一 C44/GPa 97.9 一 Cl】/GPa 327.3 一 Clz/GPa 136.2 一 fcc-FM—HS Fe (PAW)[93 fcc-AFM fct~AFM (PAW)[9] fcc—FM—LS 0.36 — 0.35 0.34 一 0.35 12.OO 12.03 1O.6O 1O.78 1I.43 1O.8O 166.2 165.9 194.7 211.9 I76 147.2 —8.O8 —8.16 —8.11 —8.21 —8.32 —8.28 2.6O 2.57 0 0 一 0.78 一89.1 一 <0 l14.2 一 85.1 —102 一 113.8 267.8 一 173.8 47.35 一 一 364.2 一 260.6 225.6 一 一 135.8 一 9O.44 fee—FM—HS fcc-Fe7oCr2[9] Fel4Cr2 fcc—AFM bcc—FM 0.36 一 035 .11.49 1O.34 lO.85 l1.5O 103.8 224.5 236.5 172.4 —8.26 —8.11 —8.30 —8.38 1.26 0.74 —0.05 1.74 —39.7 一 84.81 55.82 51.77 一 241.3 124.4 5O.9 一 349.5 246.9 130.3 一 18O.0 135.2 0.28 10 ̄Cr)[2a] fee-FM—LS fee-FM—HS fee—AFM Fe14Mn2 O・29 0.35 0.36 0.35 一 10.34 12.O2 1O.96 一 225.4 一 181.9 一 —8.19 一8.1l —8.25 一 0.43 Z.13 0.11 791 103.6 —105 123.1 120.7221.8 —55.6 219.8 78 363.5 一 346.0 127. 9156.4 一 99.81 bcc-FM 0.28 l1.26 219.3 —8.29 2.03 59.17 113.8 298.2 179.8 1O%Mn)[23] 0・0. 28 一一 一一 一一 ~ 77.81 120.9 9 285. 5 129. 9fee-FM—LS fee-FM—HS fee-AFM Fel4Coz bcc—FM 0.35 0.36 0.34 0.28 10.59 l1.79 10.70 11.47 163.9 161.6 210.2 188.0 —7.94 —7.96 —7.97 —8.11 1.13 2.43 0.16 2.31 9O.58 —91.6 66.57 65.55 169.3 一 217.6 107.6 284.7 一 299.0 275.4 103.5 一 165.9 144.3 10 ̄Co)[zs]Icc-FM—LS fcc-FM—HS fcc-Fe70Ni2[93 (0.・290.35 0.36 — —一 10.93 11.81 1092 .179. 144.5 147.0 185.6 一 —7.77 —7.78 —7.99 一1.48 2.33 0.78 41 —26.5 —5O.2 一 108.7 133.5 15.1 一 227.1 109.2 8O.12 一 130 ̄3 162.1 180.4 一 Fel4Niz fee-AFM bcc-FM (EMTo. 035 .10.87 203.7 —7.77 0.36 66.80 199.2 292.7 159.1 0.29 11.66 161.2 —7.88 2.33 47.40 102.3 224.4 129.6 10%Ni)[zs] O・29 — 148.9 一 一 62.16 108.2 231.8 107.5 此外,关于Fe的fee相AFM态的计算中发现,其能量最低点的位置发生了微小的偏移,随着c/a的增 加能量减小,所以C 为负值,此状态不稳定,而c/a=1.07的fct结构是AFM态的稳态.这与Herper[。]认为 fee相AFM态在奈尔温度以上不存在及OkatovL7 认为c/a 1.1的fct相AFM态能量最低相一致,也就是 说,fee相AFM态是一个可能但非严格的亚稳态,不可能通过热处理等方法得到块体材料的反铁磁亚稳态, 而需用衬底加以约束 ].合金化对与Fe的fee相AFM态具有稳定化作用,从表1可以看出,所有合金的反 铁磁fee态都是力学稳定的.而得到反铁磁态的前提条件是能够使合金的晶胞体积缩小,实验上最好的办法 是将薄膜生长在合适的衬底之上. ・ 186・ 内蒙古师范大学学报(自然科学汉文版) 第46卷 3 小结 运用密度泛函理论(DFT)研究了不同结构、不同磁性状态的Fe T。(T—Cr,Mn,Co,Ni)二元合金的能 量、磁矩性质以及力学性质.(1)在能量和磁矩方面,Fe和Fe T。(T—Cr,Mn,Co,Ni)二元合金的bcc相 FM态比fee相FM态、AFM态更稳定,同时在fee相的FM态中都发现具有HS和LS的特性,合金的fee 相AFM态要比FM态稳定.(2)Fe的fee相AFM态是Fe的一个可能的亚稳态,但由于力学不稳定,需用 衬底约束Fe薄膜获得fee反铁磁态.(3)对Fe和Fe T (T—Cr,Mn,Co,Ni)二元合金力学性质的计算可 知,仅有fee相FM HS态和Fe Ni fee相FM LS态在力学上不稳定,其他相在力学上是稳定的.由于铁磁 fcc至铁磁bcc容易通过c/a弛豫的Bain路径(Bain Path)而产生相变,因此尽管合金的fee铁磁态力学稳 定,但不可能在低温条件下出现,可能的办法是通过衬底束缚或者用间隙原子掺杂膨胀晶胞体积得到. 参考文献: [1-I Grimvall G.Thermophysical Properties of Materials[M].The Netherlands:Elsevier Science B V,1999:27—45. [2] Vitos L.Computational Quantum Mechanics for Materials Engineers:The EMTO Method and Applications[M]. London:Springer-Verlag,2007;103—107. P.Ab initio full—potential study of the structural and magnetic phase stability of iron [3] Herper H C,Hoffmann E,Entel[J].Phys Rev B,1999,60:3839-3848. [4] van Schilfgaarde M,Abrikosov I A,Johansson B.Origin of the Invar effect in ion-rnickel alloys[J].Nature,1999,400:46—49. shima T,Makino N.Studies on industrialization of mt permanent magnet(Ⅲ):On the Magnetic Properties and the E53 MiIndustrial Applications[M].Japan:Iron and Steel,1965,43:557,647,726. [6] Gerhard L,Yamada T K,Balashov T,et a1.Magnetoelectric coupling at metal surfaces[J].Nature Nanotch,2010(5):e792. [73 Okatov S V,Kuznetsov A R,Gornostyrev Y N,et a1.Effect of magnetic state on the a transition in iron:First—principle calculation of the bain transformation path[J].Phys Rev B,2008,79:094111(1-4). liberto O C,Aguayo A,de Coss R.Effects of Co doping on the metamagnetic states of the ferromagnetic fee Fe—Co alloy [8] Fi[J].J Phys Condens Matter,2013,25:026001-026009. .Lattice distortion induced anomalous ferromagnetism and electronic structure in [9] Wang W Y,Shang S L,Wang Y,et a1fcc Fe and Fe-TM(TM—Cr,Ni,Ta,Zr)alloys[J].Mater Chem Phys,2015,162:748—756. [1O] El1] Gebhardt T,Music D,Ekholm M,et a1.Elastic properties of fee Fe-Mn-X(X—Al,Si)alloys studied by theory and experiment[J].Acta Mater,2011,59:3145—3155. Mirzoev A A,Mirzaev D A,Yalolov M M.Energy of mixing and magnetic state of components of Fe-Mn alloys:A first— principles calculation for the ground state[J].Phys Met Metallogr,2006,101:341—348. [12] Benjamin G,Stefan B,Brtissing F,et a1.Development of magnetic moments in Fel Ni alloys[J].J Phys Condens Matter,2011,23:254210—254218. [13] H.Ⅱ.梁基谢夫.金属二元系相图手册[M].郭青蔚,译.北京:化学工业出版社,2008:408,443,593,597. ,[143 施一生.Fe1-xPdx合金电子结构和磁性的理论研究[J].物理学报,2003,52:993—998. -Cu析出相结构稳定性的影响[J].中国有色金属学报,2013,23:2840—2845. [15] 温玉锋,孙坚.合金元素镍对FCC Fe武晓霞,等.7相Fe 一 L(T—Cr,Mn,Co,Ni)合金的第一性原理研究EJ].内蒙古师范大学学报:自然科学 [163 董雪,马爽,汉文版,2016,45(2):199-205. [173 Kohn W,Sham L J.Self-Consistent Equations Including Excangeh and orCrelation Effctes口].Phys.Rev,1965,140:1133-1138. [18] BlOchl P E.Projector augmented-wave method[J].Phys Rev B,1994,50(24):17953—17979. [19] Kresse G,Hanfner J.Ab initio molecular dynamics for liquid metals[J].Phys Rev B,1993,47(1):558—561. dew J P,Burke K,Ernzerhof M.Generalized Gradient Approximation Made Simple[J].Phys Rev Lett,1996, [2o3 Per77(18):3865-3868. st H J,Pack J D.On Special Points for Brillouin Zone Integrations[J].Phys Rev B,1976,13:5188—5192. [213 Monkhornaghan F D.The Compressibility of Media Under Extreme Pressures[M].USA:Proc Natl Acad Sci,1194:244—247. [z2] Mur elastic constants of ferromagnetic bcc Fe-based random [23] Zhang H L,Punkkinen M P】,Johansson B,et a1.Single-crystalalloys from first—principles theory[J].Phys Rev B,2010,81:184105(1—14). (下转第191页) 第2期 鲍 迪等:MnFe。…Z P。_7 Ge。. 。(Z—Co,Ni)系列化合物的结构与磁性 Structure and Magnetic Properties of MnFe1-xZ Po) 77Geo23(Z=Co,Ni..BAO Di,BAO Jian—hua,LI Ying-jie,0.Tegus (College of Physics and Electronic Information,InnerMongolia Normal University,Hohhot 010022,China) Abstract:The MnFe1-xZ P077Ge023(Z—Co,Ni)compounds were prepared by high energy ball milling ..and a solid—state reaction.The structure and magnetocaloric effect and electronic states of Fe atom investi— gated by using X-ray diffraction。magnetic measurement and XANES technology.X-ray diffraction shows that all compounds crystallize in the Fe2 P—type hexagonal structure with space group P62m symmetry.The Lattice parameter a decreases and c increases with increasing the Co/Ni content.Magnetic measurements show that the Curie temperature of MnFeP0.77 Geo.z3 compound is 360 K.The Curie temperature,thermal hysteresis and isothermal magnetic entropy change of MnFe1…Z P0.77 Ge0.23(Z:Co,Ni)compounds decrease with increasing the Co/Ni content.The analysis of Fe L3 edge XANES indicates,the Fe+2 valence state increase and Fe+3 valence state decrease when substituting Co/Ni for Fe. Key words:MnFe(P,Ge)compounds;Curie temperature;magnetic entropy change;Fe L edge XANES 【责任编辑陈汉忠】 (上接第186页) First Principles Study on Structur al and Mechanical Properties of Binary Fe14T2(T=Cr, Mn,Co,Ni)Alloys MA Shuang ,DONG Xue ,WU Xiao—xia ,GUAN Peng-fei。,B.Narsu (1.College of Physical and Electronic Information,Inner Mongolia Normal University,Hohhot 010022,China; 2.College of Mathematical and Biological Engineering。Inner Mongolia University of Science and Technology, Baotou 040017,Inner Mongolia,China;3.Beijing Computational Science Research Center,Beijing 100193,China) Abstract:The projector augmented wave method based on the Density Functional Theory is employed to study the total energy,magnetic and mechanical properties of Fe“T2(T—Cr,Mn,Co,Ni)alloys with different phase structures.The results indicate that the body centered cubic phase(bcc phase)with ferro— magnetism of Fe14 T2(T=Cr,Mn,Co,Ni)binary alloy having the lowest energy and thus should be the ground state.Ferromagnetic face centered cubic phase(fcc phase)of these binary alloys show high spin (HS)and low spin(LS)characteristics,while antiferromagnetic fcc phase has a relatively lower energy in comparison with its ferromagnetic counterpart,hence antiferromagnetic phase could be more energetically more stable.The calculated elastic constants indicate that ferromagnetic fcc Fe and the HS state of ferro— magnetic Fel4 Ni2 alloy are mechanically unstable,other phase with different magnetic order are stable mechanically. Key words:Fe“T2(T—Cr,Mn,Co,Ni)alloys;phase stability;mechanical property;first—principles calculation 【责任编辑陈汉忠】