Co2-基 Heusler 合金 Co2FeAl1–xSix(x = 0.25, x =0.5, x = 0.75）的结构、电子结构及热电特性的第一性原理研究*
2019-03-16杨艳敏李佳马洪然杨广毛秀娟李聪聪
杨艳敏 李佳 马洪然 杨广 毛秀娟 李聪聪
(河北工业大学理学院,天津 300401）
运用基于密度泛函理论的第一性原理方法,对 Co2FeAl1–xSix(x=0.25,0.5,0.75）系列 Heusler合金的电子结构、四方畸变、弹性常数,声子谱以及热电特性进行了计算研究.结果显示,Co2FeAl1–xSix系列合金的电子结构均为半金属特性,向下自旋态(半导体性)均呈现良好的热电特性,并且随着硅原子浓度的增加功率因子随之增加.计算的声子谱不存在虚频,均满足动力学稳定性条件,弹性常数均满足玻恩稳定性条件,机械稳定性均良好.随着晶格常数c/a的比值变化,体系的能量最低点均出现在c/a=1处,即结构稳定性不随畸变度c/a的变化而变化,说明不存在马氏体相变.此系列合金薄膜的电子结构呈现较高的自旋极化率,在替代浓度 x=0.75 时自旋极化率达到 100%,且当 x=0.75 时薄膜在畸变度 c/a=1.2 时存在马氏体相变.随着晶格畸变度的改变,总磁矩也发生变化,且主要由Fe和Co两种过渡金属原子的磁矩变化所决定.
1 引 言
Heusler合金是一类具有丰富物理性质的功能材料,其最早报道可以追述到1903年F.Heusler关于Cu2MnAl和Cu2MnSn系列磁性高化学有序度合金的研究[1],之后的100多年里涌现出了大量的Heusler合金的相关实验和理论研究报道,研究显示其具有磁电阻效应,含稀土的Heusler合金具有超导电性、大磁感生应变效应以及半金属(halfmetal)性等很多物理效应[2−4].半金属特性是近些年发现的,因其同时具有较大的磁矩和较高的居里温度已经成为研究热点,半金属材料自旋极化的两条子能带分别呈现金属性和半导体特性,以至于载流子的自旋极化率为100%,这为操纵电子自旋自由度开辟了新道路,因此被广泛关注.半金属材料、半金属和半导体异质结构[5]等都可以应用于隧道结、自旋阀和磁随机存储器等自旋相关设备中,因此半金属材料是一类具有广泛应用前景的自旋电子学器件材料.
Co-基Heusler合金是一类重要的磁性材料,如 Co2MnAl,Co2MnSi和 Co2MnSn等在 20世纪70年代就已经被合成并深入研究[6].其中Co2MnAl和Co2MnSn最先被发现在费米能级处少数自旋态消失,导致其具有特殊的载流子输运性质[7],同时de Groot等[8]提出了所谓的半金属铁磁体(HMF)的概念,即这种材料的电子在费米能级处表现出100%的自旋极化.Co-基Heusler合金由于具有强的磁性和高的居里温度被高度关注,其中Co2Fe-基Heusler合金是自旋电子学领域研究最多的化合物,例如Co2FeAl[9]和Co2FeSi[10]已被实验证实属于半金属材料,具有100%自旋极化率,且Co2FeSi合金的居里温度达到 1100K,磁矩为 6.0,是良好的自旋电子学器件候选材料.除此之外,Si在集成电路和器件中有着广泛的应用,在Si上直接生长的超级氧化物具有优异的机械和热稳定性以及电绝缘性[11],并且近来发现含 Si元素的LaFe13－xSixHy氢化物作为磁性制冷剂具有出色的性能[12].所以在本文中,我们选用Si原子对Co2FeAl进行掺杂.另外,铁磁形状记忆合金(FSMAs)在磁场的作用下表现出丰富的物理性质,如磁场驱动形状记忆效应、磁熵变、磁场诱发应变、磁电阻和霍尔效应等,并且这些特性在磁驱动器、传感器、磁制冷等方面具有潜在应用价值[13−19],因此,寻找新的FSMA具有重要的意义.近年来,Barman等[20]的计算和实验研究表明,Ga2MnNi具有马氏体相变而成为新型的Heusler型FSMA,这为我们寻找新的Heusler型FSMA提供了重要的启示.
本文运用基于密度泛函理论(density functional theory,DFT)的第一性原理方法,研究了Co2FeAl1－xSix系列Heusler合金的电子结构、热电特性、四方畸变,声子谱及弹性特性.并且模拟了Co2FeAl1－xSix系列合金的薄膜性质,计算了其电子结构、四方相变及磁特性等.计算结果表明,Co2FeAl1－xSix系列合金均为半金属材料,声子谱的计算结果表明其动力学是稳定的,且满足力学稳定性条件.Co2FeAl1－xSix系列合金薄膜的电子结构显示其具有较高的自旋极化率,且Co2FeAl0.25Si.75满足马氏体相变的条件[21−27].
2 结构与计算方法
Heusler合金属于立方 L21结构,空间群为Fm3m,化学式可表示为X2YZ,其中X和Y为过渡金属元素,Z为主族元素(也称sp元素),各原子在一个单胞结构中的Wyckoff位置为 A(0,0,0),B(1/4,1/4,1/4）,C(1/2,1/2,1/2） 和 D(3/4,3/4,3/4）.一般主族元素占据D位,过渡金属元素占据A,B和C三个位置.Heusler合金原子占位规律往往是含较多3d电子的过渡金属元素原子占据(A,C)位,称为 Cu2MnAl型,含较少 3d 电子的过渡金属元素原子占据 (A,B),称为 Hg2CuTi型.在Co2FeAl1－xSix系列合金中,Co 原子比 Fe原子多一个 3d 电子,因此 Co 原子占据 (A,C)位,Fe原子占据B位,Al和Si原子占据D位,形成高有序的 L21面心结构,例如替代浓度 x=0.25 的结构如图1所示.我们设置铁磁和反铁磁两种磁序进行能量优化,来确定最稳定的结构以及磁序.此外,我们模拟了相应的薄膜特性,薄膜结构采用在图1基础上,在[001]方向周期性扩展三层,然后加15Å的真空层,薄膜的厚度在掺杂浓度 x=0.25,x=0.5 和 x=0.75 下分别为 1.6956,1.6982 和1.6922nm.在四方畸变计算过程中,保持体系的体积不变,沿c轴方向进行拉伸或压缩.
图1 (a)Co2FeAl1－xSix(x=0.25）的 L21 结构;(b)Co2FeAl0.75Si0.25 的薄膜结构Fig.1.(a)L21structure of Co2FeAl1－xSix(x=0.25）;(b)thin film structure of Co2FeAl0.75Si0.25.
本文的第一性原理计算采用基于赝势平面波方法的vasp软件包[28],交换关联能处理采用PBEGGA(Perdew-Burke-Ernzerh,generalized gradient approximations)方案[29−31],原子核和电子之间的相互作用采用 PAW(projector augmented wave)势方法,选用GGA近似能显著提高计算的基态性质,能给出金属正确的基态,且计算的晶格常数更贴近实验值.考虑到过渡金属元素d电子的库仑排斥,计算时加U值进行修正,即采用GGA+U方法[28],这样计算得到的过渡金属d电子符合局域特征,且带隙增大,计算结果更为准确.Gercsi等[10]研究了Co2FeSi的电子结构与磁特性,测试结果显示U值范围取1.8—2.0eV时与实验数据吻合较好.因此,本文计算中 U 值取为 2eV.为了保证计算具有良好的收敛性,平面波截断能选取为500eV.结构优化时 k 空间取点采用 5×5×5 网格,力的收敛标准取为 0.001eV/Å.总能和态密度计算采用 15×15×15 的 k 点网格,自洽迭代收敛精度选取为 1×10–5eV/原子.本文中的热电性能计算采用弛豫时间近似的玻尔兹曼输运理论,计算采用BoltzTrap软件包[32],计算过程中提高k点至10000,基本理论如下.电导张量定义为
3 结果与讨论
3.1 Co2FeAl1–xSix系列合金的结构特性与磁特性
为了获得Co2FeAl1–xSix系列合金的平衡晶格常数,分别对铁磁和反铁磁序下的Cu2MnAl型结构进行了优化,能量最低点均出现在铁磁态,优化结果如图2所示.表1为平衡晶格常数、各原子分磁矩以及合金原胞总磁矩的计算结果.可以看出三种化合物中Fe原子的总磁矩为分别为2.998,3.037 和 3.094,对合金总磁矩的贡献最大;Co原子总磁矩分别为 1.262,1.337 和 1.400,对合金总磁矩也有较大的贡献;Al和Si原子只贡献很小的磁矩,磁矩几乎为零.Heusler合金的磁矩由Slater-Pauling(SP)定律 Mt=Zt－24 描述[33],其中Mt为原胞总磁矩,Zt为原胞内的总价电子数.其中各原子价电子布局分别为 Co3d74s2,Fe 3d64s2,Al3s23p1,Si3s23p2,当 Si原子的掺杂浓度为 0.25,0.5,0.75 时,三种化合物按 SP 定律原胞总磁矩分别为 5.25,5.5,5.75,与我们计算出的磁矩非常接近,仅有约 0.2的偏离.Heusler合金的磁性与其晶体结构密切相关,如图1所示,在Co2FeAl晶体中,Co 原子占据 (A,C)位,Fe 原子占据B位,Al原子占据 D位,当Si原子替代 Al原子后,由于原子半径大小不同以及原子间相互作用的改变,使得原子之间距离发生变化,因此对3d电子杂化产生影响,从而磁矩相应产生变化.由表 1可以看到,随着 Si原子浓度的增加,Fe和Co原子总磁矩也在逐渐增大,因此总磁矩也在逐渐增加;Si和Al原子磁距变化很小.说明合金的总磁矩主要受Fe和Co原子影响.
图2 Co2FeAl1－xSix合 金 在 铁 磁 态 (FM)和 反 铁 磁 态(AFM) 下的晶格常 数优 化曲线(a)x=0.25;(b)x=0.5;(c)x=0.75Fig.2.Optimization curves of lattice constant for Co2FeAl1－ xSixalloy under ferromagnetic and antiferromagnetic magnetic order.
表1 Co2FeAl1－xSix 合金在 x=0.25,0.5,0.75 时的晶格参数及磁矩Table1. Lattice parameters and magnetic moments of Co2FeAl1－xSixalloys at x=0.25,0.5and0.75.
3.2 Co2FeAl1－xSix系列合金的电子结构
我们运用 GGA+U的方法对 Co2FeAl1–xSix(x=0.25,0.5,0.75）的能带结构进行了计算,Co2FeAl1–xSix(x=0.25,0.5,0.75）系列合金体相在布里渊区沿着特殊K点高对称方向的电子能带结构如图3(a)—(c)所示,图中的虚线表示费米能级,将费米能级处选为能量零点,能带计算选取的能量范围是–3—3eV.图 3(a)—(c)中,红色曲线为自旋向下方向的能带结构,蓝色线为自旋向上方向的能带结构,显示出在费米能级附近三种合金自旋向下方向均可以看到清晰的带隙,这表明Co2FeAl1–xSix(x=0.25,0.5,0.75）系列合金体相在自旋向上方向均具有金属性质,自旋向下方向均有半导体特性,即材料为半金属特性.从图3(a)—(c)可以看出,Co2FeAl1–xSix(x=0.25,0.5,0.75） 系列合金的电子能带的色散是非常的相似,这是因为它们的sp原子都具有相同的原子壳层.电子能带结构的带隙宽度分别为Co2FeAl0.75Si0.25(0.65eV),Co2FeAl0.5Si0.5(0.79eV)和 Co2FeAl0.25Si0.75(0.61 eV).
我们使用 GGA+U 方法对 Co2FeAl1–xSix系列合金的电子结构进行了计算,计算的态密度如图 4(a)—(c)所示.三种替代浓度下,即 Co2FeAl0.75Si0.25,Co2FeAl0.5Si0.5和 Co2FeAl0.25Si0.75的 向 下自旋态在费米能级附近有明显的带隙,向上自旋态是金属特性,自旋极化率公式为%,其中N↑和N↓分别代表自旋向上和向下态在费米面处的态密度[34],计算出三种化合物的自旋极化率P均为100%,呈现很好的半金属特性.由图中的各原子的轨道态密度可以看出,各个合金的Co原子的向下自旋态在–0.5和1.2eV处 (x=0.25）,–1.25 和 1.02eV 处 (x=0.5）,–0.81和 0.82eV 处 (x=0.75）分别有两个峰,这是由于在立方晶体场下3d电子的eg与t2g轨道产生交换劈裂所导致的,Fe原子的态密度同样具有此特点.因此,过渡金属Co原子和Fe原子都具有显著的3d电子劈裂特性.由于强烈的3d轨道杂化,使得总态密度在费米能级上下分别形成了明显的杂化峰,位于费米能级下方的峰值出现在–0.51,–1.25eV 和–0.81eV 处,费米能级以上峰值出现在 1.18,0.81 和 0.82eV 处.可见,合金带隙两侧的峰恰恰和Co原子以及Fe原子态密度中的那两个峰完全对应,表明Co与Fe原子的3d电子劈裂以及它们之间的杂化是形成带隙的主要原因,并决定了带隙的宽度.
图3 (a) Co2FeAl0.75Si0.25, (b) Co2FeAl0.5Si0.5和 (c)Co2FeAl0.25Si0.75的能带结构Fig.3.Energy band structure of(a) Co2FeAl0.75Si0.25,(b)Co2FeAl0.5Si0.5and(c)Co2FeAl0.25Si0.75.
图4 (a) Co2FeAl0.75Si0.25,(b) Co2FeAl0.5Si0.5 和 (c)Co2FeAl0.25Si0.75的总态密度和分态密度Fig.4.Thetotaland atom-projected density of states for Heusler alloys Co2FeAl1－xSix(x=0.25,0.5,0.75）film in(a),(b)and(c).
3.3 Co2FeAl1–xSix系列合金的热电输运特性
热电材料种类很多,如 Bi2Te3,PbTe,SiGe,GeTe 及 half-Heusler合金等.迄今为止,发现的已有100多种[35],其中half-Heusler合金是一类重要的中温区热电材料(500—800K),由于Heusler合金在结构上与half-Heusler合金非常相近,因此近年来其热电性能也引起了广泛关注[36,37].
我们采用弛豫时间近似的玻尔兹曼输运理论分别计算了 Co2FeAl1–xSix系列合金在 100,300,500,900 和 1200K 温度下的热电输运特性,由于本系列合金没有相关热电数据报道,因此,弛豫时间采取拟合与本体系接近的Co2FeSi的电导率实验数据得到[38],拟合的弛豫时间结果为 1.3×10–14s.计算的输运特性如图 5(a)—(i)所示.图 5(a),(e)和 (i)分别为 Co2FeAl1–xSix在掺杂浓度 x为0.25,0.5和0.75下的向下自旋态(半导体特性)的Seebeck系数随化学势的变化.如图所示,分别在化学势正值区域(电子掺杂)和负值区域(空穴掺杂)出现峰值,且峰值均随温度的增大而减小,当T=100K 时,Seebeck 系数均可取得最大值.当 x=0.25时,位于负值区域的峰值为235.98,位于正值区域的峰值的绝对值为197.70,当 x=0.5 和 x=0.75 时,位于负值区域的峰值分别为329.86和376.61,位于正值区域的峰值的绝对值分别为278.18和328.93,因此,三种替代浓度空穴掺杂下易获得较高的Seebeck系数.
图 5(b),(f)和 (j)分 别 为 x= 0.25,0.5 和0.75的电导率随化学势的变化.可以看出在化学势–0.35—0.4eV范围内,电导均随着化学势的增加而增加,随着温度的升高而升高,显示出明显的半导体特征.图 5(c),(g)和 (k)分别为 x=0.25,0.5和0.75时的热导率随化学势的变化.可以看出热导均随化学势的增加略有增加,随温度升高而显著增加.图 5(d),(h)和 (l)分别为 x=0.25,0.5和0.75的功率因子随化学势的变化,三种化合物的功率因子在化学势正值区域和负值区域各有一个峰值,总趋势都是随着温度的升高而增大.在温度 T=1200K 时达到最大值,x=0.25,0.5 和0.75时负值区域的峰值分别为163.60,173.09和176.29,正值区域的峰值分别为52.47,53.00 和 54.63,位于正值区域的峰值均大于位于负值区域的峰值,因此空穴掺杂比电子掺杂更能获得高的功率因子.
图5 Co2FeAl0.75Si0.25 向下自旋态的 (a)Seebeck 系数,(b)电导,(c)热导和 (d) 功率因子随化学势的变化;Co2FeAl0.5Si0.5 向下自旋态的 (e)Seebeck 系数,(f) 电导,(g) 热导和 (h)功率因子随化学势的变化;Co2FeAl0.25Si0.75 向下自旋态的 (i)Seebeck 系数,(j)电导,(k)热导和(l)功率因子随化学势的变化Fig.5.The transport properties with variation of chemical potential for Co2FeAl1－xSix(x=0.25,0.5and0.75）.The case of x=0.25corresponds to(a),(b),(c)and(d),and the case of x=0.5corresponds to(e),(f),(g)and(h),and the case of x=0.75 corresponds to (i),(j),(k) and (l).The four columns from left to right correspond to the Seebeck coefficients S,electrical conductivity ,electronic thermal conductivity and PF( ),respectively.
在室温 300K 时,x=0.25,0.5 和 0.75 时的化合物分别对应的功率因子值为28.70,30.07,30.74,这和广泛商用的热电材料PbTe在室温下的功率因子21.0相接近[39],所以这三种化合物有望成为有价值的候选热电材料.总体来看,随着Si原子替代浓度的增加,Seebeck系数峰值不断增大,电导率、热导与功率因子的值均在逐渐增加,由此可见Co2FeAl1–xSix在Si的掺杂浓度为0.75下的热电性能最好.
3.4 Co2FeAl1–xSix系列合金的声子振动谱及比热容
为了验证Co2FeAl1–xSix系列合金的动力学稳定性,我们计算了它们的声子谱,如图6(a)—(c)所示.对于 L21相的 Co2FeAl1–xSix系列合金,每个立方晶胞具有16个原子,因此有48个声子振动模式(3n个自由度),其中3个为声学支,其余是光学支[40].由图可以看出在X--L高对称方向之间,色散曲线处于声子频率正值范围内,没有出现虚频,因此它们满足动力学稳定性.图6(d)为计算的三种化合物的热容量随温度的变化曲线,三种化合物的声子谱色散曲线极为相似,因此图6(d)中三种化合物的定压热容Cv随温度的变化曲线也很相近,热容在温度从零到大约350K之间显着增加,而在大于约350K之后缓慢上升,且基本趋于一个定值,这与经典的Dulong-Petit定律相一致[41].在 0—350K 范围,比热随着Si原子替代浓度增加而逐渐增加.
3.5 Co2FeAl1–xSix系列合金的弹性性质及力学稳定性
对于立方相,只有三个弹性常数是独立的,即C11,C12和C44.一个稳定的立方相应该满足的机械稳定性标准,即 Born 标准为:C11>0,C44>0,C11>C12,C11+ 2C12>0.已 知 Voigt模 量 和Reuss模量是实际有效模量的上限和下限,Voigt-Reus-Hill模量是Voigt模量和Reuss模量的算术平均[41].对于立方相,Voigt剪切模量 GV,Reuss剪切模量GR,Voigt-Reuss-Hill剪切模量GH以及体模量B由以下算式给出:
图6 Co2FeAl1－xSix 合 金 在 x=0.25,0.5,0.75 时 的 声 子 谱 及 比 热 容 (a)Co2FeAl1－xSix(x=0.25）,(b)Co2FeAl1－xSi x(x=0.5） 和 (c)Co2FeAl1－xSi x(x=0.75）的声子谱;(d)Co2FeAl1－ xSix(x=0.25,0.5,0.75） 的比热容随温度的变化Fig.6.Full phonon spectra of Co2FeAl1－xSix(x=0.25,0.5and0.75）alloys in(a),(b)and(c).The temperature dependent heat capacity Cvwith an inset graph showing the temperaturefrom180K to250K in(d).
表2 计算的 Co2FeAl1－xSix(x=0.25,x=0.5,x=0.75）合金的弹性常数、体模量及剪切模量Table2. The calculated cubic elastic constant C11,C12,C44,shear modulus Gv,GRand GHin GPa.
利用第一性原理计算的弹性常数、体模量及剪切模量如表2所列.三种化合物的弹性常数均满足上述机械稳定性条件:C11>0,C44>0,C11>C12,C11+2C12>0,说明它们在立方相下均是稳定的.三种化合物在替代浓度 x=0.25,0.5 和0.75时,分别对应的体模量 B的值为 193.77,184.43和 92.85GPa,可见随着 Si原子的增加体模量在减小,当x=0.25时合金的体模量值最大,即合金的抗压缩能力最好.在 x=0.75 时,合金的体模量值最小,抗压缩能力相对较弱.三种化合物在替代浓度 x=0.25,0.5 和 0.75 时,分别对应的剪切模量 GH值为 86.04,101.25 和 100.42GPa,可以看到当x=0.5时合金的剪切模量值最大,即抵抗切应变的能力最大,合金的刚性最强.当x=0.25时合金的剪切模量值最小,合金的刚性最弱.
3.6 Co2FeAl1–xSix系列合金薄膜的电子结构
近些年,半金属薄膜材料越来越受到研究者的关注.例如:Okamura 等[42]将 Co2FeAl作为铁磁层制备出Co2FeAl/Al-Ox/Co75Fe25磁隧道结,得到 47%的室温磁电阻.文献[43−48]在 MgO基底001面上制备了Co2FeAl0.5Si0.5(CFAS)/Ag/CFAS的 结 构,膜 厚 度 为 20nm,证 明 了 使 用 CFAS Heusler合金作为铁磁电极在Ag作为底部缓冲层时可产生相对大的巨磁阻(MR),并且通过增加两个CFAS层的有序度,预期MR值的进一步增强.
图7 (a) Co2FeAl0.75Si0.25, (b) Co2FeAl0.5Si0.5和 (c)Co2FeAl0.25Si0.75薄膜的总态密度和原子分态密度Fig.7.Thetotaland atom-projected density of states for Co2FeAl1－ xSix(x= 0.25,0.5 and 0.75） film in (a),(b)and(c).
本节对Co2FeAl1－xSix系列合金薄膜的电子结构进行了模拟计算,建模采用如图1(a)所示的晶胞,沿z方向扩充3倍再加真空层,结构图如图 1(b)所示.当 x=0.25,x=0.5,x=0.75 时,膜的厚度分别为 1.6956,1.6982和 1.6922nm.图 7(a)—(c) 分别为 x=0.25,x=0.5,x=0.75时的态密度,费米能级处的电子自旋极化率分别为 43.09%,59.07%,100%,由此看出随着 Si原子的增加,自旋极化率在不断升高,且当 x=0.75 时,合金的自旋极化率达到100%,属于半金属,可见Co2FeAl1–xSix系列合金的薄膜也可以呈现出很好的半金属特性.图 7(a)—(c)显示,掺杂浓度 x=0.25,0.5,0.75 时,自旋向下方向的态密度分别在费米能级之上 0.35,0.31,0.1eV 处出现峰值,化合物总的态密度主要来自过渡金属元素Fe和Co原子的3d电子之间的强烈的d-d杂化.而且,总态密度与Co原子的分态密度趋势一致,这也说明薄膜材料的电子结构主要受Co原子的影响.
3.7 Co2FeAl1－xSix系列合金体相及其薄膜的四方畸变
图8 (a)x=0.25,(b)x=0.5 和 (c)x=0.75 替代浓度下 Co2FeAl1－xSix合金体相的总能量差 与畸变度 c/a 的关系;(d)x=0.25,(e)x=0.5 和 (f)x=0.75 替代浓度下 Co2FeAl1－xSix薄膜的驱动力 与畸变度 c/a 的关系Fig.8.Calculated total energies as a function of the c/a ratio for Co2FeAl1－xSix(x=0.25,0.5and0.75）Heusler alloys in(a),(b)and(c)andfilm materials in(d),(e)and(f).
近年来,铁磁形记忆合金(FSMAS)在磁场的作用下表现出很好的物理性质而备受关注,目前已经报道过的Heusler型FSMAs主要有Ni-基,Mn-基,Fe-基和 Co-基体系.基于前人的结论可知,Heusler合金体系要发生马氏体相变需要满足一定的条件:首先要有充足的相变驱动力即四方相与立方相之间的能量差值;其次,畸变度,即晶格常数c与a的比值也是关键,满足马氏体相变的适当的畸变度 c/a通常在1.2—1.3之间[18].图8(a)—(f)分别为Co2FeAl1–xSix系列合金体相及其薄膜的相变驱动力随畸变度的变化.可以看到,体相的Co2FeAl1–xSix系列合金能量最低点均出现在c/a=1 处,没有发生马氏体相变.Co2FeAl1–xSix系列合金薄膜的能量最低点均未出现在畸变度c/a=1 处,当掺杂浓度 x=0.25,0.5,0.75 时,能量最低点分别出现在畸变度 c/a=1.05,c/a=1.02,c/a=1.2 处,即在 x=0.25,x=0.5 时不存在马氏体相变,而在x=0.75时存在马氏体相变.
3.8 Co2FeAl1–xSix系列合金薄膜的磁性在四方畸变下的响应
我们知道,在立方相到四方相的结构畸变的过程中,由于原子间距发生改变,原子轨道之间的杂化强度会发生变化,从而影响电子态密度,进而影响磁特性.图 9(a)—(c)分别为 Co2FeAl1–xSix系列合金薄膜在掺杂浓度x=0.25,0.5和0.75下总磁矩及各原子磁矩随畸变度的变化.由图可知,当x=0.25 时,当 c/a<1 即对 c 轴进行压缩时,总磁矩随着c/a的值的增大而增大,由图可知Fe原子的磁矩变化决定了总磁矩的变化,在1.0
图9 (a)x=0.25,(b)x=0.5 和 (c)x=0.75 替代浓度下Co2FeAl1－xSix合金薄膜的总磁矩及各原子总磁矩随畸变度的变化Fig.9.The total magnetic moment and the magnetic moment of each atom of Co2FeAl1－xSixfilm change with distortion at x=0.25,x=0.5and x=0.75in(a),(b)and(c).
4 结 论
本文运用基于DFT的第一性原理方法,对Co2FeAl1–xSix(x=0.25,0.5,0.75） 系列合金体相及其薄膜的结构特性、电子结构,磁性及热电输运特性进行了计算研究.通过分析计算结果得出以下结论:Co2FeAl1–xSix(x=0.25,0.5,0.75）立方相均为半金属材料,且均机械稳定性及动力学稳定性良好,当 x=0.25 时材料抗压缩能力最好,当 x=0.5时材料的刚性最强,三种材料在随着Si原子替代浓度的增加,自旋向下 (半导体特性)态的Seebeck系数峰值不断增大,电导、热导与功率因子的值均在逐渐增加,是有价值的候选热电材料.合金的总磁矩主要由Fe和Co两种过渡金属原子磁矩所决定.在畸变度 c/a=0.3—1.4 范围内,三种替代浓度下的体相均未出现马氏体相变.Co2FeAl1－xSix(x=0.25,0.5,0.75） 系列合金薄膜的自旋极化率随着Si原子的增加而增加,当x=0.75 时自旋极化率为 100%,且当 x=0.75 时结构的能量最低点在c/a=1.2处,存在马氏体相变.薄膜材料的总磁矩也主要由Fe和Co两种过渡金属原子决定.
