-
Among the spectra of rare-earth permanent magnetic materials, the Sm-Co-based alloys stand out with their excellent magnetic properties in high-temperature environments. However, the practical applications of these alloys in high-temperature settings face constraints due to their comparatively lower saturation magnetization and structural stability. In this study, Fe, Ni, Cu, and Zr are used as representative transition metal elements to investigate the effects of doping elements on the structural stability, magnetic properties, and electronic structure of SmCo3 alloy by first-principles calculations. The findings indicate that the doping of elements Ni, Cu, and Fe contributes positively to enhancing the structural stability of the SmCo3, while the introduction of Zr element has an adverse effect. Magnetic property calculations reveal that the incorporation of non-magnetic elements leads the total magnetic moment of the SmCo3 to decrease to a certain extent, whereas the introduction of magnetic elements can enhance the total magnetic moment. Notably, not all doped magnetic elements in the SmCo3 result in an increasing total magnetic moment. The underlying microscopic mechanisms are elucidated through electronic structure analysis. Finally, it is screened out that the transition element Fe is beneficial to improving the magnetic properties and structural stability of SmCo3, and the doping concentration (atomic percentage) in its unit cell ranges from 0 to 22.22%, the optimal doping concentration (atomic percentage) is predicted to be 18.52%.
1. 引 言
Sm-Co基稀土永磁合金具有优异的高温磁性能, 被广泛应用于航空航天、新能源和国防军工等领域, 尤其是在工作温度超过500 ℃的条件下有着无可替代的地位[1–6]. 因此, 其从1967年[7]发现至今得到了广泛的关注和研究. Sm-Co基稀土永磁合金具有多种构型, 如SmCo3, SmCo5和Sm2Co17等. SmCo3晶体结构可以由SmCo5晶体结构通过Sm原子替代一些Co原子位点得到[8], 其轴比c/a > 4.8[9], 具有较大的晶格各向异性, 有利于磁性能的提高. 但是由于其较低的饱和磁化强度[10]限制了最大磁能积, 从而阻碍了实际应用.
前人的研究表明, 掺杂合金元素是增强Sm-Co基稀土永磁合金的稳定性和磁性能的理想途径之一[11–14]. 实验方面, Guo等[15]采用Mn掺杂SmCo7纳米晶结构, 发现Mn元素掺杂既可以提高结构稳定性, 又可以提升体系综合磁性能, 并探究出最佳掺杂浓度. Wang等[16]以Fe, Ni和Nb为掺杂元素, 发现元素Fe和Ni可以有效提高Sm5Co19基纳米合金的综合磁性能. 理论计算方面, Das等[17]基于第一性原理计算, 在SmCo5体系中用Fe替代Co, 对不同浓度Fe掺杂SmCo5合金的晶体结构稳定性、磁化强度和磁晶各向异性能进行研究, 计算发现高浓度Fe的掺杂不利于晶体结构稳定性, 同时使磁晶各向异性能下降, 但是可以提高体系的饱和磁化强度, Liu和Altounian[18]也得出类似的结论. Landa等[19]通过第一性原理计算, 在SmCo5体系中用Fe替代大部分Co并使用Ni作为热力学稳定剂, 产生了具有优异磁性的SmCoNiFe3合金, 磁化强度与SmCo5合金相比从8.27 μB增至10 μB, 且具有高的居里温度和强的磁各向异性, 其磁各向异性能大小约为Nd2Fe14B合金的两倍. 无独有偶, Antoniou等[20]也同样发现在SmCo5体系中掺杂Ni可以提高掺杂体系的结构稳定性. 综上所述, 实验探索和理论计算两方面均表明掺杂过渡元素对于提升Sm-Co基稀土永磁合金的结构稳定性和磁性能有着非常重要的作用. 而第一性原理计算可以从体系的能量、晶体结构和电子结构等微观层面分析材料的稳定性和相关磁性能, 从而筛选出有利于提升Sm-Co基稀土永磁合金综合磁性能的掺杂元素, 有望从理论上为解决国内外稀土永磁新材料与新产品研发周期长和成本高的缺点提供有益指导.
目前, 国内外对于SmCo3体系的实验探索和理论研究还非常少, 尤其是掺杂过渡元素的理论研究, 且掺杂元素种类和掺杂浓度在SmCo3合金中的优先占位情况和对磁性能的影响机制尚不明确. 本文基于第一性原理计算, 在SmCo3体系中设计多种过渡元素的掺杂构型. 通过对掺杂体系的晶格常数、掺杂元素的优先占位、电子结构和磁矩等方面的计算, 得出了过渡元素掺杂对SmCo3体系结构稳定性和磁性能的作用关系, 并揭示了其微观机理. 此外, 基于掺杂结构稳定性、磁性能和掺杂元素的浓度三方面作用因素, 筛选出了最佳的掺杂元素和其掺杂浓度, 为新型高性能Sm-Co基永磁合金的研究和开发提供了有益的理论指导.
2. 模型构建和计算方法
2.1 模型构建
SmCo3的晶体结构属于PuNi3型三方(菱方)晶系, 空间群为Rˉ3m. 如图1(a)所示, 其单胞中有36个原子, 3个Sm原子占据3a晶格位置, 6个Sm原子占据6c晶格位置, 3个Co原子占据3b晶格位置, 6个Co原子占据6c晶格位置, 18个Co原子占据18h晶格位置. 在三维坐标系中, 5个不等价的晶格位置坐标分别为: Sm(3a)原子坐标为(0, 0, 0); Sm(6c)原子坐标为(2/3, 1/3, 0.193699); Co(3b)原子坐标为(1/3, 2/3, 1/6); Co(6c)原子坐标为(1/3, 2/3, 0.332163); Co(18h)原子坐标为(2/3, 5/6, 0.253176), 其实验晶格常数为a = b = 5.050 Å, c = 24.590 Å, c/a = 4.8693[10]. 掺杂过渡元素以Fe, Ni, Cu和Zr为代表, 建立未掺杂的SmCo3晶胞(1×1×1)和掺杂原子百分比为3.7%的Sm9Co26M (M = Fe, Ni, Cu, Zr)的模型, 而且在计算过程中考虑掺杂元素占据Co的3b, 6c和18h的晶格位置3种情况, 如图1所示.
图 1 (a) SmCo3的晶体结构和Sm9Co26M晶体结构中M占据(b) 3b, (c) 6c和(d) 18h晶格位置Fig. 1. Crystal structures of SmCo3 (a) and Sm9Co26M with M occupied (b) 3b site, (c) 6c and (d) 18h site.2.2 计算方法及相关参数
本研究的计算都是在VASP (Vienna ab-initio simulation package)[21,22]中的密度泛函理论(DFT)框架下, 通过平面波基投影增强波(PAW)方法[23,24]进行的. 使用广义梯度近似下的平面波(PBE)泛函[25,26]. 平面波截断能设置为500 eV, 原子的自洽总能量收敛标准和几何(弛豫)结构力的收敛标准分别设为10–5 eV和0.01 eV/Å, K点设置为9×9×3. 在结构弛豫过程中, 原子位置以及晶胞的形状和体积都被允许弛豫. 为了更好地描述d轨道和f轨道电子, 本文采用Dudarev提出的DFT+U方法[27], 对Sm, Fe, Co, Ni, Cu和Zr原子的有效U值分别设置为4.70 eV, 1.83 eV, 2.22 eV, 1.61 eV, 4.80 eV和1.50 eV, 这些U值的选取文献[28, 29]. 计算时各原子的电子组态分别为: Sm-5s25p65d16s2, Fe-3d64s2, Co-3d74s2, Ni-3d84s2, Cu-3d104s1和Zr-4d25s2.
Sm元素的赝势选择Sm_3进行计算, 也是VASP官方推荐使用的赝势[30,31]. 该赝势将4f电子作为原子核内电子, 不作为价电子, 本文研究的是掺杂过渡元素种类和浓度对SmCo3体系结构稳定性和磁矩的影响, 而不需要考虑4f-3d轨道电子相互作用的影响, 因此选择该赝势是合理的, 也是研究Sm-Co合金体系常用的处理方式.
3. 计算结果分析与讨论
3.1 掺杂元素占位与结构稳定性
首先, 对掺杂前后的模型进行结构优化, SmCo3晶胞的晶格常数a = b = 5.0123 Å, c = 24.6424 Å, 体积V = 536.17 Å3, 与实验晶格参数a = b = 5.050 Å, c = 24.590 Å, V = 543.09 Å3[10]相吻合. 如表1所列, 优化后的晶格常数与实验结果相差仅1%左右, 计算得到的SmCo3晶胞体积比实验体积约小2%, 说明计算结果是可靠的. 由于掺杂元素的原子半径不同, SmCo3掺杂体系的c/a略有不同, 且掺杂体系的晶胞体积随着掺杂元素的原子半径的增大而增大.
表 1 SmCo3和Sm9Co26M的晶格参数和晶胞体积Table 1. Lattice parameters and cell volumes of SmCo3 and Sm9Co26M.Model Lattice parameter/Å c/a V/Å3 a b c SmCo3 5.0123 5.0123 24.6424 4.9164 536.17 Experiments[10] 5.050 5.050 24.590 4.8693 543.09 Sm9Co26Ni 4.9852 4.9852 24.5962 4.9338 536.76 Sm9Co26Fe 5.0428 5.0428 24.6771 4.8935 537.62 Sm9Co26Cu 4.9793 4.9793 24.5967 4.9398 538.11 Sm9Co26Zr 5.0670 5.0670 24.8397 4.9022 548.71 为研究掺杂元素种类和占位对晶体结构稳定性的影响, 计算了掺杂Sm9Co26M体系的替代能, 具体公式为[32]
E(sub)=E(Sm9Co27−nMn)−E(SmCo3)+nE(Co)−nE(M), (1) 式中, E(Sm9Co27–nMn)表示掺杂体系的总能量, E(SmCo3)表示未掺杂体系的总能量, E(Co)和E(M)分别表示Co和过渡金属元素M处于单质状态时每个原子的基态能量, n表示原子个数. 替代能为负值表示该元素掺杂较容易, 且有利于提升掺杂体系的结构稳定性; 反之正值则表示该元素掺杂较难同时不利于体系结构稳定性. 元素在不同晶格位置取代Co原子时, 优先占据替代能较低的位点, 不同掺杂体系的替代能如图2所示. 从图2可以发现Fe, Ni, Cu和Zr都优先占据18h位点. 替代能越低, 晶体结构越稳定, Ni, Cu和Fe的掺杂提高了SmCo3结构的稳定性, 与之前的研究结果一致[17,20]; 而Zr的掺杂则不利于SmCo3结构的稳定性. 插图显示掺杂元素在优先占位的体系替代能, 由此可以得出掺杂元素对SmCo3结构的稳定能力由大到小依次为: Ni, Cu, Fe, Zr.
图 2 掺杂元素在3b, 6c和18h位置的替代能, 插图表示不同元素在优先占位的Sm9Co26M的替代能; 插图中红色虚线表示未掺杂SmCo3的基态能量, 其替代能为0作为参考值Fig. 2. Substitution energy of the doping element at 3b, 6c and 18h sites, the insert shows the substitution energy of Sm9Co26M with different elements at the preferential site. The red dashed line in the illustration represents the ground state energy of undoped SmCo3, and its substitution energy is 0 eV as a reference value.3.2 差分电荷密度分析
为阐明SmCo3掺杂M (M = Fe, Ni, Cu, Zr)后体系中各类原子间的电荷转移机制, 计算了不同M原子掺杂在18h晶位(001)晶面的差分电荷密度, 结果如图3所示. 图3(c), (d)中Ni和Cu原子替代Co原子后, 在Co-Ni和Co-Cu原子间聚集大量电子, 且倾向于向掺杂原子Ni和Cu原子转移. 这表明图3(c), (d)中Co-Ni和Co-Cu原子间金属键得到明显加强, 且大于图3(a)中本身的Co—Co键, 从而提升了掺杂体系稳定性; 图3(e)中Zr原子替代Co原子后, Zr原子周围表现为明显的电子缺失, 且在Co-Zr原子间的差分电荷密度没有重叠, 表明Co-Zr原子间的相互作用减弱, 从而降低了掺杂体系稳定性; 图3(b)中Fe原子替代Co原子与未掺杂图3(a)中Co原子之间电荷密度差异并不明显; 上述分析与3.1节中替代能计算结果定性一致. 造成原子间电子转移情况的不同主要归因于每个原子的电负性各不相同, 从大到小排列顺序依次为Ni, Cu, Co, Fe, Zr.
图 3 (001)面上(a) SmCo3和(b)—(e) Sm9Co26M (M = Fe, Ni, Cu, Zr)的差分电荷密度; 红色表示电荷的积累, 蓝色表示电荷的损耗Fig. 3. Difference charge density of (a) SmCo3 and (b)—(e) Sm9Co26M (M = Fe, Ni, Cu, Zr) on (001) plane; the red indicates enrichment of electrons , blue indicates loss of electrons.3.3 掺杂元素种类对基体总磁矩的影响
图4展现了各个掺杂体系总磁矩的计算结果, 图中虚线表示未掺杂SmCo3体系的总磁矩为38.74 μB, 与实验值37.35 μB接近[10]. 导致理论和实验磁矩的微小差别可能有以下两种原因: 一是在实验制备SmCo3合金的过程中, 并不能完美地控制Sm:Co的含量比为1∶3; 另一种是计算中只考虑了原子的自旋磁矩, 而忽略了原子的轨道磁矩. 从图4可知: Ni, Cu和Zr元素的掺杂使体系的总磁矩减弱, 主要原因是Cu和Zr元素是非磁性元素; 而Ni作为磁性元素并没有增大体系的总磁矩, 这与传统磁性元素[12,13]掺杂提升体系饱和磁化强度的认识不符; 磁性元素Fe掺杂使体系的总磁矩大幅增大. 此外, 当M原子在不同晶格位置取代Co原子时, 对总磁矩的影响也各不相同. 从图4可以看出, 掺杂元素Fe, Ni, Cu在不同位置总磁矩的差异很小. 但是掺杂元素Zr的体系在取代不同位置的Co时, 对总磁矩的影响很大, 其中掺杂在3b, 6c和18h晶格位置体系的总磁矩分别为38.73 μB, 37.43 μB和37.02 μB. 主要原因是Zr的原子半径远远大于Co原子和不同位置的晶格间隙不同, Zr原子掺杂后造成晶格内附近的Co原子位置发生不同程度的位移, 其中3b晶格位置间隙最小, 使周围Co原子之间距离减小, 增强了其周围Co原子的磁矩, 因此对体系总磁矩的减小程度产生了一定的补偿作用.
图 4 SmCo3和Sm9Co26M (M = Fe, Ni, Cu, Zr)的总磁矩; 红色虚线表示未掺杂SmCo3体系的总磁矩为38.74 μBFig. 4. Total magnetic moments of SmCo3 and Sm9Co26M (M = Fe, Ni, Cu, Zr); the red dotted line indicates that the total magnetic moment of the undoped SmCo3 system is 38.74 μB.3.4 偏态密度分析
为进一步揭示掺杂元素对SmCo3体系总磁矩影响的微观机制, 选择最稳定的掺杂结构(18h)进行元素电子的偏态密度(partial density of state)分析, 如图5所示. 灰色阴影代表上下自旋态密度对称部分, 而红色阴影代表不同部分. 体系的净磁矩源于上下自旋态密度的不对称性. 如图5(a), (b)所示, Fe在18h晶格位置的上下自旋态密度对称性远低于Co, 说明Fe原子磁矩大于Co原子磁矩, 所以Fe掺杂可以提升体系的总磁矩; 图5(c)表示Ni元素在18h晶格位置的偏态密度, 由于其费米能级附近上下自旋电子数较少, 且上下自旋态密度对称性较高, 说明其磁矩很小, 远远小于Co元素在18h晶格位置的磁矩, 因此Ni的掺杂降低了体系的总磁矩; 图5(d), (e)分别表示Cu元素和Zr元素在18h晶格位置的偏态密度, 两者上下自旋态密度都高度对称. 前者在费米能级附近上下自旋电子数稀少, 对应红色区域表示的自旋差值部分也非常少, 且在较深能级处出现了上下自旋对称灰色部分, 说明Cu原子本身具有的磁矩非常小. 而后者在费米能级附近上下自旋电子数相对较多, 且红色自旋差值部分也较大, 但是远远小于Co元素在该位置的红色区域, 说明Zr原子具有一定大小的磁矩但小于Co原子的磁矩. 原子Cu, Zr与原子Co的本征自旋磁矩方向相反, 造成体系总磁矩下降, 但是Zr原子磁矩大于Cu原子磁矩, 所以Zr原子掺杂对体系的总磁矩影响更大. 经过上述PDOS分析, 与上文中磁矩计算结果一致.
图 5 未掺杂体系中Co原子和掺杂体系中掺杂原子M在18h晶格位置的PDOSFig. 5. PDOS of Co atom in undoped system and doped atom M in different doped systems at 18h lattice position.3.5 元素的筛选
如图6所示, 给出不同掺杂体系最稳定掺杂结构(18h)的替代能和总磁矩. 蓝色虚线代表未掺杂SmCo3体系的总磁矩, 红色虚线代表未掺杂SmCo5体系的替代能, 其值为0. 从图6可以看出, 红线左侧掺杂元素的体系具有较低的替代能, 有利于PuNi3型晶体结构稳定性. 掺杂体系替代能从小到大排列依次是Ni , Cu, Fe以及不利于结构稳定的Zr. 所以从结构稳定性方面选择的最佳元素是Ni, Cu次之, 最后是Fe. 另一方面, 蓝色虚线以上的掺杂体系, 只有Sm9Co26Fe体系总磁矩比未掺杂的总磁矩高, 而蓝色虚线以下的体系总磁矩则比未掺杂的低. 因此, 在4种过渡金属元素中只有Fe元素掺杂既有利于提升结构稳定性, 又能提高SmCo3体系的总磁矩.
图 6 Sm9Co26M (M = Fe, Ni, Cu, Zr)合金的替代能和总磁矩Fig. 6. Substitution energy and total magnetic moment of the Sm9Co26M (M = Fe, Ni, Cu, Zr) alloys.3.6 掺杂元素浓度对结构稳定性和总磁矩的影响
通过上述对4种过渡元素掺杂SmCo3体系进行系统计算分析, 得到掺杂不同种类过渡元素对SmCo3体系结构稳定性、总磁矩和电子结构可以产生很大的影响. 在实际的Sm-Co稀土材料研发中, 仅仅了解掺杂元素的种类并不足以确保材料性能的最优化. 为实现最佳效果, 必须精确调控各组成的含量. 如在SmCo5体系中掺杂Fe时, 在合适的浓度范围内可以提高饱和磁化强度和居里温度, 但过量掺杂会导致饱和磁化强度、居里温度和结构稳定性下降, 使其综合磁性能受到影响[17,18,33]. 为了进一步分析4种过渡元素掺杂浓度对体系的影响, 并预测合适的掺杂浓度范围. 本文参照实验可行性掺杂浓度范围[34], 在SmCo3单原胞中, 掺杂原子百分比为0—22.22%范围内, 根据原子个数配比设置了计算原子百分比分别为x = 3.7%, 7.4%, 11.11%, 14.81%, 18.52%, 22.22%的掺杂体系.
如图7所示, 分别计算了SmCo3–xMx体系中替代能和体系总磁矩随着掺杂浓度的变化趋势(结构都是最稳定的), 红色(黑色)折线表示体系总磁矩(替代能)随着掺杂浓度的演化关系. 如图7(a)所示, 在SmCo3–xFex体系中, 可看出随着Fe掺杂浓度的增大, 掺杂体系的总磁矩也随之增大, 同时替代能随之降低, 表明掺杂体系结构也越来越稳定; 图7(b), (c)分别表示SmCo3–xNix和SmCo3–xCux体系的演化情况, 可以看出, 随着掺杂浓度的增大, 体系的替代能几乎呈线性方式下降, 表明体系的结构越来越稳定, 但总磁矩下降; 图7(d)表示SmCo3–xZrx体系, 伴随着掺杂浓度的增大, 体系的替代能几乎呈线性方式上升, 表明体系的结构越来越不稳定, 同时掺杂体系的总磁矩也在下降. 综合4种SmCo3基掺杂体系随着掺杂浓度变化的替代能和总磁矩特点, 可以预测在掺杂原子百分比为0—22.22%范围内, Fe元素为最佳的掺杂原子, 与3.5节叙述结果相对应, 且其最优的掺杂原子百分比为18.52%左右(考虑到实验的可行性掺杂浓度).
图 7 SmCo3–xMx (M = Fe (a), Ni (b), Cu (c), Zr (d))体系替代能和体系总磁矩随掺杂浓度的变化Fig. 7. Relationship between the substitution energy and total magnetic moment of SmCo3–xMx (M = Fe (a), Ni (b), Cu (c), Zr (d) system with doping concentration.4. 结 论
采用第一性原理计算方法, 研究过渡金属M = Fe, Ni, Cu, Zr掺杂的SmCo3型永磁合金的结构稳定性和磁性能. 主要结论如下:
1)通过计算各掺杂体系的替代能可得, 过渡元素Fe, Ni, Cu和Zr都优先占据18h晶格位置, 元素Ni, Cu和Fe的掺杂可以降低掺杂体系的替代能, 从而提升SmCo3体系结构稳定性; 而元素Zr掺杂不利于体系结构稳定性.
2)添加磁性元素Fe可以使体系总磁矩从38.74 μB增至40.16 μB, 而添加磁性元素Ni却使体系总磁矩从38.74 μB降低到37.67 μB. 这说明在SmCo3体系并不是所有的铁磁元素都可以增加体系的总磁矩. 掺杂Cu和Zr这种非磁性元素通常会在一定程度上减弱SmCo3体系的总磁矩.
3)在SmCo3单原胞中, 掺杂原子百分比为0—22.22%范围内, 通过计算获得了SmCo3不同掺杂体系结构稳定性和磁矩随着掺杂元素浓度变化规律. 最后筛选出既能提高掺杂体系结构稳定性, 又可以提升磁性能的元素: Fe, 并预测了其最合适的掺杂原子百分比为18.52%.
[1] Guo Z H, Pan W, Li W 2006 J. Magn. Magn. Mater. 303 396
Google Scholar
[2] 张杨, 宋晓艳, 徐文武, 张哲旭 2012 物理学报 61 016102
Google Scholar
Zhang Y, Song X Y, Xu W W, Zhang Z X 2012 Acta Phys. Sin. 61 016102
Google Scholar
[3] Gutfleisch O, Willard M A, Brück E, Chen C H, Sankar S G, Liu J P 2011 Adv. Mater. 23 821
Google Scholar
[4] 周相龙, 袁涛, 宋欣, 贾文涛, 马天宇 2020 中国材料进展 39 384
Google Scholar
Zhou X L, Yuan T, Song X, Jia W T, Ma T Y 2020 Mater. China 39 384
Google Scholar
[5] 张昌文, 李华, 董建敏, 王永娟, 潘凤春, 郭永权, 李卫 2005 物理学报 54 1814
Google Scholar
Zhang C W, Li H, Dong J M, Wang Y J, Pan F C, Guo Y Q, Li W 2005 Acta Phys. Sin. 54 1814
Google Scholar
[6] Xue Z Q, Liu L, Liu Z, Li M, Lee D, Chen R J, Guo Y Q, Yan A R 2016 Scr. Mater. 113 226
Google Scholar
[7] Strnat K, Hoffer G, Olson J, Ostertag W, Becker J J 1967 J. Appl. Phys. 38 1001
Google Scholar
[8] Streever R L 1979 Phys. Rev. B. 19 2704
Google Scholar
[9] Buschow K H J, Van der Goot A S 1968 J. Less Common Met. 14 323
Google Scholar
[10] Gaidukova I Y, Granovsky S A, Markosyan A S, Rodimin V E 2006 J. Magn. Magn. Mater. 301 124
Google Scholar
[11] Saito T, Nishio-Hamane D 2014 J. Alloys Compd. 585 423
Google Scholar
[12] Guo K, Lü H, Xu G J, Liu D, Wang H B, Liu X M, Song X Y 2022 Mater. Today Chem. 25 100983
Google Scholar
[13] Mao F, Lü H, Liu D, Guo K, Tang F W, Song X Y 2019 J. Alloys Compd. 810 151888
Google Scholar
[14] 毛斐, 吕皓, 唐法威, 郭凯, 刘东, 宋晓艳 2021 金属学报 57 948
Google Scholar
Mao F, Lü H, Tang F W, Guo K, Song X Y 2021 Acta. Metall. Sin. 57 948
Google Scholar
[15] Guo K, Lu H, Mao F, Liu D, Tang F W, Wang H B, Song X Y 2020 Nanoscale 12 5567
Google Scholar
[16] Wang D, Liu D, Hou C, Wang H B, Liu X, Song X Y 2017 J. Alloys Compd. 717 93
Google Scholar
[17] Das B, Choudhary R, Skomski R, Balasubramanian B, Pathak A K, Paudyal D, Sellmyer D J 2019 Phys. Rev. B 100 024419
Google Scholar
[18] Liu X B, Altounian Z 2011 Comput. Mater. Sci. 50 841
[19] Landa A, Söderlind P, Parker D, Åberg D, Lordi V, Perron A, Turchi P E A, Chouhan R K, Paudyal D, Lograsso T A 2018 J. Alloys Compd. 765 659
Google Scholar
[20] Antonioua E, Semprosa G, Gjokab M, Sarafidisa C, Polatogloua H M, Kioseogloua J 2021 J. Alloys Compd. 882 160699
Google Scholar
[21] Kresse G, Furthmüller J 1996 Comput. Mater. Sci. 6 15
Google Scholar
[22] Kresse G, Furthmüller J 1996 Phys. Rev. B 54 11169
Google Scholar
[23] Kresse G, Joubert D 1999 Phys. Rev. B 59 1758
[24] Blöchl P E 1994 Phys. Rev. B. 50 17953
Google Scholar
[25] Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865
Google Scholar
[26] Perdew J P, Chevary J A, Vosko S H, Jackson K A, Pederson M R, Singh D J, Fiolhais C 1992 Phys. Rev. B 46 6671
Google Scholar
[27] Dudarev S L, Botton G A, Savrasov S Y, Humphreys C J, Sutton A P 1998 Phys. Rev. B 57 1505
Google Scholar
[28] Van der Marel D, Sawatzky G A 1988 Phys. Rev. B 37 10674
Google Scholar
[29] Zhang Z F, Guo Y Z, Robertson J 2020 Phys. Rev. Mater. 4 054603
Google Scholar
[30] Brooks M S S, Eriksson O, Wills J M, Johansson B 1997 Phys. Rev. Lett. 79 2546
[31] Georg K, Martijn M, Jurgen F 2018 User Manual for VASP Version 5.4
[32] Liu B J, Wang H, Xu C, Liu X P, Zhang Q F, Zhan T L, Stamenov P, Coey J M D, Jiang C B 2020 Phys. Rev. Mater. 4 044406
Google Scholar
[33] Söderlind P, Landa A, Locht I L, Åberg D, Kvashnin Y, Pereiro M, Däne M, Turchi P E A, Antropov V P, Eriksson O 2017 Phys. Rev. B 96 100404
Google Scholar
[34] Raja A, Adhikary T, Al-Omari I A, Das G P, Ghosh S, Satapathy D K, Oraon A, Shield J E, Aich S 2020 J. Magn. Magn. Mater. 504 166645
Google Scholar
-
图 1 (a) SmCo3的晶体结构和Sm9Co26M晶体结构中M占据(b) 3b, (c) 6c和(d) 18h晶格位置
Figure 1. Crystal structures of SmCo3 (a) and Sm9Co26M with M occupied (b) 3b site, (c) 6c and (d) 18h site.
图 2 掺杂元素在3b, 6c和18h位置的替代能, 插图表示不同元素在优先占位的Sm9Co26M的替代能; 插图中红色虚线表示未掺杂SmCo3的基态能量, 其替代能为0作为参考值
Figure 2. Substitution energy of the doping element at 3b, 6c and 18h sites, the insert shows the substitution energy of Sm9Co26M with different elements at the preferential site. The red dashed line in the illustration represents the ground state energy of undoped SmCo3, and its substitution energy is 0 eV as a reference value.
图 3 (001)面上(a) SmCo3和(b)—(e) Sm9Co26M (M = Fe, Ni, Cu, Zr)的差分电荷密度; 红色表示电荷的积累, 蓝色表示电荷的损耗
Figure 3. Difference charge density of (a) SmCo3 and (b)—(e) Sm9Co26M (M = Fe, Ni, Cu, Zr) on (001) plane; the red indicates enrichment of electrons , blue indicates loss of electrons.
图 4 SmCo3和Sm9Co26M (M = Fe, Ni, Cu, Zr)的总磁矩; 红色虚线表示未掺杂SmCo3体系的总磁矩为38.74 μB
Figure 4. Total magnetic moments of SmCo3 and Sm9Co26M (M = Fe, Ni, Cu, Zr); the red dotted line indicates that the total magnetic moment of the undoped SmCo3 system is 38.74 μB.
图 5 未掺杂体系中Co原子和掺杂体系中掺杂原子M在18h晶格位置的PDOS
Figure 5. PDOS of Co atom in undoped system and doped atom M in different doped systems at 18h lattice position.
图 6 Sm9Co26M (M = Fe, Ni, Cu, Zr)合金的替代能和总磁矩
Figure 6. Substitution energy and total magnetic moment of the Sm9Co26M (M = Fe, Ni, Cu, Zr) alloys.
图 7 SmCo3–xMx (M = Fe (a), Ni (b), Cu (c), Zr (d))体系替代能和体系总磁矩随掺杂浓度的变化
Figure 7. Relationship between the substitution energy and total magnetic moment of SmCo3–xMx (M = Fe (a), Ni (b), Cu (c), Zr (d) system with doping concentration.
表 1 SmCo3和Sm9Co26M的晶格参数和晶胞体积
Table 1. Lattice parameters and cell volumes of SmCo3 and Sm9Co26M.
Model Lattice parameter/Å c/a V/Å3 a b c SmCo3 5.0123 5.0123 24.6424 4.9164 536.17 Experiments[10] 5.050 5.050 24.590 4.8693 543.09 Sm9Co26Ni 4.9852 4.9852 24.5962 4.9338 536.76 Sm9Co26Fe 5.0428 5.0428 24.6771 4.8935 537.62 Sm9Co26Cu 4.9793 4.9793 24.5967 4.9398 538.11 Sm9Co26Zr 5.0670 5.0670 24.8397 4.9022 548.71 
DownLoad: CSV
-
[1] Guo Z H, Pan W, Li W 2006 J. Magn. Magn. Mater. 303 396
Google Scholar
[2] 张杨, 宋晓艳, 徐文武, 张哲旭 2012 物理学报 61 016102
Google Scholar
Zhang Y, Song X Y, Xu W W, Zhang Z X 2012 Acta Phys. Sin. 61 016102
Google Scholar
[3] Gutfleisch O, Willard M A, Brück E, Chen C H, Sankar S G, Liu J P 2011 Adv. Mater. 23 821
Google Scholar
[4] 周相龙, 袁涛, 宋欣, 贾文涛, 马天宇 2020 中国材料进展 39 384
Google Scholar
Zhou X L, Yuan T, Song X, Jia W T, Ma T Y 2020 Mater. China 39 384
Google Scholar
[5] 张昌文, 李华, 董建敏, 王永娟, 潘凤春, 郭永权, 李卫 2005 物理学报 54 1814
Google Scholar
Zhang C W, Li H, Dong J M, Wang Y J, Pan F C, Guo Y Q, Li W 2005 Acta Phys. Sin. 54 1814
Google Scholar
[6] Xue Z Q, Liu L, Liu Z, Li M, Lee D, Chen R J, Guo Y Q, Yan A R 2016 Scr. Mater. 113 226
Google Scholar
[7] Strnat K, Hoffer G, Olson J, Ostertag W, Becker J J 1967 J. Appl. Phys. 38 1001
Google Scholar
[8] Streever R L 1979 Phys. Rev. B. 19 2704
Google Scholar
[9] Buschow K H J, Van der Goot A S 1968 J. Less Common Met. 14 323
Google Scholar
[10] Gaidukova I Y, Granovsky S A, Markosyan A S, Rodimin V E 2006 J. Magn. Magn. Mater. 301 124
Google Scholar
[11] Saito T, Nishio-Hamane D 2014 J. Alloys Compd. 585 423
Google Scholar
[12] Guo K, Lü H, Xu G J, Liu D, Wang H B, Liu X M, Song X Y 2022 Mater. Today Chem. 25 100983
Google Scholar
[13] Mao F, Lü H, Liu D, Guo K, Tang F W, Song X Y 2019 J. Alloys Compd. 810 151888
Google Scholar
[14] 毛斐, 吕皓, 唐法威, 郭凯, 刘东, 宋晓艳 2021 金属学报 57 948
Google Scholar
Mao F, Lü H, Tang F W, Guo K, Song X Y 2021 Acta. Metall. Sin. 57 948
Google Scholar
[15] Guo K, Lu H, Mao F, Liu D, Tang F W, Wang H B, Song X Y 2020 Nanoscale 12 5567
Google Scholar
[16] Wang D, Liu D, Hou C, Wang H B, Liu X, Song X Y 2017 J. Alloys Compd. 717 93
Google Scholar
[17] Das B, Choudhary R, Skomski R, Balasubramanian B, Pathak A K, Paudyal D, Sellmyer D J 2019 Phys. Rev. B 100 024419
Google Scholar
[18] Liu X B, Altounian Z 2011 Comput. Mater. Sci. 50 841
[19] Landa A, Söderlind P, Parker D, Åberg D, Lordi V, Perron A, Turchi P E A, Chouhan R K, Paudyal D, Lograsso T A 2018 J. Alloys Compd. 765 659
Google Scholar
[20] Antonioua E, Semprosa G, Gjokab M, Sarafidisa C, Polatogloua H M, Kioseogloua J 2021 J. Alloys Compd. 882 160699
Google Scholar
[21] Kresse G, Furthmüller J 1996 Comput. Mater. Sci. 6 15
Google Scholar
[22] Kresse G, Furthmüller J 1996 Phys. Rev. B 54 11169
Google Scholar
[23] Kresse G, Joubert D 1999 Phys. Rev. B 59 1758
[24] Blöchl P E 1994 Phys. Rev. B. 50 17953
Google Scholar
[25] Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865
Google Scholar
[26] Perdew J P, Chevary J A, Vosko S H, Jackson K A, Pederson M R, Singh D J, Fiolhais C 1992 Phys. Rev. B 46 6671
Google Scholar
[27] Dudarev S L, Botton G A, Savrasov S Y, Humphreys C J, Sutton A P 1998 Phys. Rev. B 57 1505
Google Scholar
[28] Van der Marel D, Sawatzky G A 1988 Phys. Rev. B 37 10674
Google Scholar
[29] Zhang Z F, Guo Y Z, Robertson J 2020 Phys. Rev. Mater. 4 054603
Google Scholar
[30] Brooks M S S, Eriksson O, Wills J M, Johansson B 1997 Phys. Rev. Lett. 79 2546
[31] Georg K, Martijn M, Jurgen F 2018 User Manual for VASP Version 5.4
[32] Liu B J, Wang H, Xu C, Liu X P, Zhang Q F, Zhan T L, Stamenov P, Coey J M D, Jiang C B 2020 Phys. Rev. Mater. 4 044406
Google Scholar
[33] Söderlind P, Landa A, Locht I L, Åberg D, Kvashnin Y, Pereiro M, Däne M, Turchi P E A, Antropov V P, Eriksson O 2017 Phys. Rev. B 96 100404
Google Scholar
[34] Raja A, Adhikary T, Al-Omari I A, Das G P, Ghosh S, Satapathy D K, Oraon A, Shield J E, Aich S 2020 J. Magn. Magn. Mater. 504 166645
Google Scholar
下载:
DownLoad:
Catalog
Metrics
- Abstract views: 3810
- PDF Downloads: 144
