-
Heusler合金是一类具有丰富物理性质的功能材料, 其最早报道可以追述到1903年F. Heusler关于Cu2MnAl和Cu2MnSn系列磁性高化学有序度合金的研究[1], 之后的100多年里涌现出了大量的Heusler合金的相关实验和理论研究报道, 研究显示其具有磁电阻效应, 含稀土的Heusler合金具有超导电性、大磁感生应变效应以及半金属(half-metal)性等很多物理效应[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合金的居里温度达到1100 K, 磁矩为6.0
${{\text{μ}}_{\rm{B}}}$ , 是良好的自旋电子学器件候选材料. 除此之外, 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].
-
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.6922 nm. 在四方畸变计算过程中, 保持体系的体积不变, 沿c轴方向进行拉伸或压缩.
图 1 (a) Co2FeAl1-xSix (x = 0.25)的L21结构; (b) Co2FeAl0.75Si0.25的薄膜结构
Figure 1. (a) L21 structure of Co2FeAl1-xSix (x = 0.25); (b)thin film structure of Co2FeAl0.75Si0.25.
本文的第一性原理计算采用基于赝势平面波方法的vasp软件包[28], 交换关联能处理采用PBE-GGA(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.0 eV时与实验数据吻合较好. 因此, 本文计算中U值取为2 eV. 为了保证计算具有良好的收敛性, 平面波截断能选取为500 eV. 结构优化时k空间取点采用5 × 5 × 5网格, 力的收敛标准取为0.001 eV/Å. 总能和态密度计算采用15 × 15 × 15的k点网格, 自洽迭代收敛精度选取为1 × 10–5 eV/原子. 本文中的热电性能计算采用弛豫时间近似的玻尔兹曼输运理论, 计算采用BoltzTrap软件包[32], 计算过程中提高k点至10000, 基本理论如下. 电导张量定义为
${{{\sigma }}_{\alpha \beta }}\left( \varepsilon \right) = \frac{{\rm{1}}}{N}\sum\limits_{i,k} {{{{\sigma }}_{\alpha \beta }}\left( {i,k} \right)} \frac{{\sigma \left( {\varepsilon - {\varepsilon _{i,k}}} \right)}}{{d\varepsilon }},$ 其中N为采样点的数目,
${\varepsilon _{i,k}}$ 为电子能带能量.输运张量包含: 电导张量
${{{\varepsilon }}_{\alpha \beta }}$ , 群速度${\upsilon _{\alpha \beta }}$ , 电子热导$\kappa _{\alpha \beta }^{\rm{0}}$ 和Seebeck系数${S_{i,j}}$ . 它们均能通过电导率分别被计算:${{\sigma }_{\alpha \beta }}\left( {T;\mu } \right) = \frac{{\rm{1}}}{\varOmega }\int {{{\sigma }_{\alpha \beta }}\left( \varepsilon \right)} \left[ { - \frac{{\partial {f_\mu }\left( {T;\varepsilon } \right)}}{{{\partial _\varepsilon }}}} \right] d\varepsilon ,$ ${\nu _{\alpha \beta }}\left( {T;\mu } \right) = \frac{1}{{eT\varOmega }}\int {{s_{\alpha \beta }}\left( \varepsilon \right)\left( {\varepsilon - \mu } \right)\left[ { - \frac{{\partial {f_\mu }\left( {T;\varepsilon } \right)}}{{\partial \varepsilon }}} \right]} d\varepsilon ,$ $\kappa _{\alpha \beta }^{\rm{0}}\left( {T;\mu } \right) = \frac{1}{{eT\varOmega }}\!\int \!\!{{{{\sigma }}_{\alpha \beta }}\left( \varepsilon \right)\!{{\left( {\varepsilon - \mu } \right)}^2}\!\left[ { - \frac{{\partial {f_\mu }\left( {T;\varepsilon } \right)}}{{\partial \varepsilon }}} \right]}d\varepsilon ,$ 其中,
$\mu $ 是电子化学势,$\varOmega $ 是晶胞的体积, e是电子电荷, T是绝热温度,${f_\mu }$ 是分布函数. -
为了获得Co2FeAl1–xSix系列合金的平衡晶格常数, 分别对铁磁和反铁磁序下的Cu2MnAl型结构进行了优化, 能量最低点均出现在铁磁态, 优化结果如图2所示. 表1为平衡晶格常数、各原子分磁矩以及合金原胞总磁矩的计算结果. 可以看出三种化合物中Fe原子的总磁矩为分别为2.998, 3.037和3.094
${{\text{μ}}_{\rm{B}}}$ , 对合金总磁矩的贡献最大; Co原子总磁矩分别为1.262, 1.337和1.400${{\text{μ}}_{\rm{B}}}$ , 对合金总磁矩也有较大的贡献; Al和Si原子只贡献很小的磁矩, 磁矩几乎为零. Heusler合金的磁矩由Slater-Pauling(SP)定律Mt = Zt-24描述[33], 其中Mt为原胞总磁矩, Zt为原胞内的总价电子数. 其中各原子价电子布局分别为Co 3d74s2, Fe 3d64s2, Al 3s23p1, Si 3s23p2, 当Si原子的掺杂浓度为0.25, 0.5, 0.75时, 三种化合物按SP定律原胞总磁矩分别为5.25, 5.5, 5.75${{\text{μ}}_{\rm{B}}}$ , 与我们计算出的磁矩非常接近, 仅有约0.2${{\text{μ}}_{\rm{B}}}$ 的偏离. Heusler合金的磁性与其晶体结构密切相关, 如图1所示, 在Co2FeAl晶体中, Co原子占据(A, C)位, Fe原子占据B位, Al原子占据D位, 当Si原子替代Al原子后, 由于原子半径大小不同以及原子间相互作用的改变, 使得原子之间距离发生变化, 因此对3d电子杂化产生影响, 从而磁矩相应产生变化. 由表1可以看到, 随着Si原子浓度的增加, Fe和Co原子总磁矩也在逐渐增大, 因此总磁矩也在逐渐增加; Si和Al原子磁距变化很小. 说明合金的总磁矩主要受Fe和Co原子影响.a/Å mAl/${{\text{μ}}_{\rm{B}}}$ mSi/${{\text{μ}}_{\rm{B}}}$ mFe/${{\text{μ}}_{\rm{B}}}$ mCo/${{\text{μ}}_{\rm{B}}}$ Mt/${{\text{μ}}_{\rm{B}}}$ Co2FeAl0.75Si0.25 5.6520 –0.053 –0.039 2.998 1.262 5.473 Co2FeAl0.5Si0.5 5.6607 –0.046 –0.028 3.037 1.337 5.688 Co2FeAl0.25Si0.75 5.6406 –0.038 –0.012 3.094 1.400 5.891 表 1 Co2FeAl1-xSix合金在x = 0.25, 0.5, 0.75时的晶格参数及磁矩
Table 1. Lattice parameters and magnetic moments of Co2FeAl1-xSix alloys at x = 0.25, 0.5 and 0.75.
-
我们运用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—3 eV. 图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.65 eV), Co2FeAl0.5Si0.5(0.79 eV)和Co2FeAl0.25Si0.75(0.61 eV).
图 3 (a) Co2FeAl0.75Si0.25, (b) Co2FeAl0.5Si0.5和(c) Co2FeAl0.25Si0.75的能带结构
Figure 3. Energy band structure of (a) Co2FeAl0.75Si0.25, (b) Co2FeAl0.5Si0.5 and (c) Co2FeAl0.25Si0.75.
我们使用GGA + U方法对Co2FeAl1–xSix系列合金的电子结构进行了计算, 计算的态密度如图4(a)—(c)所示. 三种替代浓度下, 即Co2FeAl0.75Si0.25, Co2FeAl0.5Si0.5和Co2FeAl0.25Si0.75的向下自旋态在费米能级附近有明显的带隙, 向上自旋态是金属特性, 自旋极化率公式为
$P = ( {N_ \uparrow }{\rm{ - }}{N_ \downarrow }/{N_ \uparrow } $ $+ {N_ \downarrow } ) \times {\rm{100}}$ %, 其中N↑和N↓分别代表自旋向上和向下态在费米面处的态密度[34], 计算出三种化合物的自旋极化率P均为100%, 呈现很好的半金属特性. 由图中的各原子的轨道态密度可以看出, 各个合金的Co原子的向下自旋态在–0.5和1.2 eV处(x = 0.25), –1.25和1.02 eV处(x = 0. 5), –0.81 和0.82 eV处(x = 0.75)分别有两个峰, 这是由于在立方晶体场下3d电子的eg与t2g轨道产生交换劈裂所导致的, Fe原子的态密度同样具有此特点. 因此, 过渡金属Co原子和Fe原子都具有显著的3d电子劈裂特性. 由于强烈的3d轨道杂化, 使得总态密度在费米能级上下分别形成了明显的杂化峰, 位于费米能级下方的峰值出现在–0.51, –1.25 eV和–0.81 eV处, 费米能级以上峰值出现在1.18, 0.81和0.82 eV处. 可见, 合金带隙两侧的峰恰恰和Co原子以及Fe原子态密度中的那两个峰完全对应, 表明Co与Fe原子的3d电子劈裂以及它们之间的杂化是形成带隙的主要原因, 并决定了带隙的宽度. -
热电材料种类很多, 如Bi2Te3, PbTe, SiGe, GeTe及half-Heusler合金等. 迄今为止, 发现的已有100多种[35], 其中half-Heusler合金是一类重要的中温区热电材料(500—800 K), 由于Heusler合金在结构上与half-Heusler合金非常相近, 因此近年来其热电性能也引起了广泛关注[36,37].
我们采用弛豫时间近似的玻尔兹曼输运理论分别计算了Co2FeAl1–xSix系列合金在100, 300, 500, 900和1200 K温度下的热电输运特性, 由于本系列合金没有相关热电数据报道, 因此, 弛豫时间采取拟合与本体系接近的Co2FeSi的电导率实验数据得到[38], 拟合的弛豫时间结果为1.3 × 10–14 s. 计算的输运特性如图5(a)—(i)所示. 图5(a), (e)和(i)分别为Co2FeAl1–xSix在掺杂浓度x为0.25, 0.5和0.75下的向下自旋态(半导体特性)的Seebeck系数随化学势的变化. 如图所示, 分别在化学势正值区域(电子掺杂)和负值区域(空穴掺杂)出现峰值, 且峰值均随温度的增大而减小, 当T = 100 K时, Seebeck系数均可取得最大值. 当x = 0.25时, 位于负值区域的峰值为235.98
${{\text{μ}}{\rm{V}}} \cdot {{\rm{K}}^{ - {\rm{1}}}}$ , 位于正值区域的峰值的绝对值为197.70${{\text{μ}}{\rm{V}}} \cdot {{\rm{K}}^{ - {\rm{1}}}}$ , 当x = 0.5和x = 0.75时, 位于负值区域的峰值分别为329.86和376.61${{\text{μ}}{\rm{V}}} \cdot {{\rm{K}}^{ - 1}}$ , 位于正值区域的峰值的绝对值分别为278.18和328.93${{\text{μ}}{\rm{V}}} \cdot {{\rm{K}}^{ - 1}}$ , 因此, 三种替代浓度空穴掺杂下易获得较高的Seebeck系数.图 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)功率因子随化学势的变化
Figure 5. The transport properties with variation of chemical potential
$\mu $ for Co2FeAl1-xSix(x = 0.25, 0.5 and 0.75). The case of x = 0.25 corresponds to (a), (b), (c) and (d), and the case of x = 0.5 corresponds 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$\sigma $ , electronic thermal conductivity${\kappa _{\rm{e}}}$ and PF (${S^2}\sigma $ ), respectively.图5(b), (f)和(j)分别为x = 0.25, 0.5和0.75的电导率随化学势的变化. 可以看出在化学势–0.35—0.4 eV范围内, 电导均随着化学势的增加而增加, 随着温度的升高而升高, 显示出明显的半导体特征. 图5(c), (g)和(k)分别为x = 0.25, 0.5和0.75时的热导率随化学势的变化. 可以看出热导均随化学势的增加略有增加, 随温度升高而显著增加. 图5(d), (h)和(l)分别为x = 0.25, 0.5和0.75的功率因子随化学势的变化, 三种化合物的功率因子在化学势正值区域和负值区域各有一个峰值, 总趋势都是随着温度的升高而增大. 在温度T = 1200 K时达到最大值, x = 0.25, 0.5和0.75时负值区域的峰值分别为163.60, 173.09和176.29
${{\text{μ}}{\rm{W}}} \cdot {\rm{c}}{{\rm{m}}^{ - {\rm{1}}}} \cdot {{\rm{K}}^{ - 2}}$ , 正值区域的峰值分别为52.47, 53.00和54.63${{\text{μ}}{\rm{W}}} \cdot {\rm{c}}{{\rm{m}}^{ - {\rm{1}}}} \cdot {{\rm{K}}^{ - 2}}$ , 位于正值区域的峰值均大于位于负值区域的峰值, 因此空穴掺杂比电子掺杂更能获得高的功率因子.在室温300 K时, x = 0.25, 0.5和0.75时的化合物分别对应的功率因子值为28.70, 30.07, 30.74
${{\text{μ}}{\rm{W}}} \cdot {\rm{c}}{{\rm{m}}^{ - {\rm{1}}}} \cdot {{\rm{K}}^{ - 2}}$ , 这和广泛商用的热电材料PbTe在室温下的功率因子21.0${{\text{μ}}{\rm{W}}} \cdot {\rm{c}}{{\rm{m}}^{ - {\rm{1}}}} \cdot {{\rm{K}}^{ - 2}}$ 相接近[39], 所以这三种化合物有望成为有价值的候选热电材料. 总体来看, 随着Si原子替代浓度的增加, Seebeck系数峰值不断增大, 电导率、热导与功率因子的值均在逐渐增加, 由此可见Co2FeAl1–xSix在Si的掺杂浓度为0.75下的热电性能最好. -
为了验证Co2FeAl1–xSix系列合金的动力学稳定性, 我们计算了它们的声子谱, 如图6(a)—(c)所示. 对于L21相的Co2FeAl1–xSix系列合金, 每个立方晶胞具有16个原子, 因此有48个声子振动模式(3n个自由度), 其中3个为声学支, 其余是光学支[40]. 由图可以看出在X-
$\Gamma $ -L高对称方向之间, 色散曲线处于声子频率正值范围内, 没有出现虚频, 因此它们满足动力学稳定性. 图6(d)为计算的三种化合物的热容量随温度的变化曲线, 三种化合物的声子谱色散曲线极为相似, 因此图6(d)中三种化合物的定压热容Cv随温度的变化曲线也很相近, 热容在温度从零到大约350 K之间显着增加, 而在大于约350 K之后缓慢上升, 且基本趋于一个定值, 这与经典的Dulong-Petit定律相一致[41]. 在0—350 K范围, 比热随着 Si原子替代浓度增加而逐渐增加.图 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)的比热容随温度的变化
Figure 6. Full phonon spectra of Co2FeAl1-xSix (x = 0.25, 0.5 and 0.75) alloys in (a), (b) and (c). The temperature dependent heat capacity Cv with an inset graph showing the temperaturefrom 180 K to 250 K in (d).
-
对于立方相, 只有三个弹性常数是独立的, 即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由以下算式给出:
$\begin{split} &{G_{\rm{V}}} \!=\! \frac{1}{5}\!\left[ {\left( {{C_{11}} \!-\! {C_{12}}} \right) \!+\! 3{C_{44}}} \right]\!,{B_{\rm{V}}} \!=\! \frac{1}{3}\!\left( {{C_{11}}\!\!+\! 2{C_{12}}} \right)\!,\\ &{G_{\rm{R}}} \!=\! \frac{5}{{4({S_{11}} - {S_{12}}) + 3{S_{44}}}};\;{B_{\rm{R}}} = \frac{1}{{3{S_{11}} + 6{S_{12}}}}; \end{split}$ $\begin{split} {C_{44}} =\; & S_{44}^{ - 1},\;\;{{C_{11}}} - {C_{12}} = {\left( {{S_{11}} - {S_{12}}} \right)^{ - 1}}, \\ & {C_{11}} + 2{C_{12}} = {\left( {{S_{11}} + 2{S_{12}}} \right)^{ - 1}}; \end{split} $ $\begin{split} {G_{\rm{R}}} & = {\left[ {\frac{4}{5}{{\left( {{C_{11}} - {C_{12}}} \right)}^{ - 1}} + \frac{3}{5}C_{44}^{ - 1}} \right]^{ - 1}}, \\ {B_{\rm{R}}} & = \frac{1}{3}\left[ {\left( {{C_{11}} + 2{C_{12}}} \right)} \right] = {B_{\rm{V}}}. \end{split} $ 所以,
$\begin{split} G =& \frac{1}{2}\left( {{G_{\rm{V}}}\! +\! {G_{\rm{R}}}} \right) = \frac{{\rm{1}}}{{\rm{2}}}\left\{ {\frac{1}{5}\left[ {\left( {{C_{11}} \!-\! {C_{12}}} \right) \!+\! 3{C_{44}}} \right]} \right.\\ & + \left. {{{\left[ {\frac{4}{5}{{\left( {{C_{11}} - {C_{12}}} \right)}^{ - 1}} + \frac{3}{5}C_{44}^{ - 1}} \right]}^{ - 1}}} \right\} \\ B = & {B_{\rm{V}}} = {B_{\rm{R}}} = \left( {{C_{{\rm{11}}}} + {\rm{2}}{{\rm{C}}_{{\rm{12}}}}} \right)/3. \end{split}$ 利用第一性原理计算的弹性常数、体模量及剪切模量如表2所列. 三种化合物的弹性常数均满足上述机械稳定性条件: C11 > 0, C44 > 0, C11 > C12, C11 + 2C12 > 0, 说明它们在立方相下均是稳定的. 三种化合物在替代浓度x = 0.25, 0.5和0.75时, 分别对应的体模量B的值为193.77, 184.43和92.85 GPa, 可见随着Si原子的增加体模量在减小, 当x = 0.25时合金的体模量值最大, 即合金的抗压缩能力最好. 在x = 0.75时, 合金的体模量值最小, 抗压缩能力相对较弱. 三种化合物在替代浓度x = 0.25, 0.5和0.75时, 分别对应的剪切模量GH值为86.04, 101.25和100.42 GPa, 可以看到当x = 0.5时合金的剪切模量值最大, 即抵抗切应变的能力最大, 合金的刚性最强. 当x = 0.25时合金的剪切模量值最小, 合金的刚性最弱.
C11/GPa C12/GPa C44/GPa B/GPa GV/GPa GR/GPa GH/GPa Co2FeAl0.75Si0.25 247.38 166.97 142.33 193.77 101.48 70.60 86.04 Co2FeAl0.5Si0.5 266.15 143.57 141.73 184.43 109.55 92.94 101.25 Co2FeAl0.25Si0.75 176.46 51.042 137.67 92.85 107.69 93.14 100.42 表 2 计算的Co2FeAl1-xSix (x = 0.25, x = 0.5, x = 0.75)合金的弹性常数、体模量及剪切模量
Table 2. The calculated cubic elastic constant C11, C12, C44, shear modulus Gv, GR and GH in GPa.
-
近些年, 半金属薄膜材料越来越受到研究者的关注. 例如: Okamura等[42]将Co2FeAl作为铁磁层制备出Co2FeAl/Al-Ox/Co75Fe25磁隧道结, 得到47%的室温磁电阻. 文献[43-48]在MgO基底001面上制备了Co2FeAl0.5Si0.5 (CFAS)/Ag/CFAS的结构, 膜厚度为20 nm, 证明了使用CFAS Heusler合金作为铁磁电极在Ag作为底部缓冲层时可产生相对大的巨磁阻(MR), 并且通过增加两个CFAS层的有序度, 预期MR值的进一步增强.
本节对Co2FeAl1-xSix系列合金薄膜的电子结构进行了模拟计算, 建模采用如图1(a)所示的晶胞, 沿z方向扩充3倍再加真空层, 结构图如图1(b)所示. 当x = 0.25, x = 0.5, x = 0.75时, 膜的厚度分别为1.6956, 1.6982和1.6922 nm. 图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.1 eV处出现峰值, 化合物总的态密度主要来自过渡金属元素Fe和Co原子的3d电子之间的强烈的d-d杂化. 而且, 总态密度与Co原子的分态密度趋势一致, 这也说明薄膜材料的电子结构主要受Co原子的影响.
-
近年来, 铁磁形记忆合金(FSMAS)在磁场的作用下表现出很好的物理性质而备受关注, 目前已经报道过的Heusler型FSMAs主要有Ni-基, Mn-基, Fe-基和Co-基体系. 基于前人的结论可知, Heusler合金体系要发生马氏体相变需要满足一定的条件: 首先要有充足的相变驱动力即四方相与立方相之间的能量差值
$\Delta E$ ; 其次, 畸变度, 即晶格常数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时存在马氏体相变.图 8 (a) x = 0.25,(b) x = 0.5和(c) x = 0.75替代浓度下Co2FeAl1-xSix合金体相的总能量差
$\Delta E$ 与畸变度c/a的关系; (d) x = 0.25, (e) x = 0.5和(f) x = 0.75替代浓度下Co2FeAl1-xSix薄膜的驱动力$\Delta E$ 与畸变度c/a的关系Figure 8. Calculated total energies as a function of the c/a ratio for Co2FeAl1-xSix (x = 0.25, 0.5 and 0.75) Heusler alloys in (a), (b) and (c) andfilm materials in (d), (e) and (f).
-
我们知道, 在立方相到四方相的结构畸变的过程中, 由于原子间距发生改变, 原子轨道之间的杂化强度会发生变化, 从而影响电子态密度, 进而影响磁特性. 图9(a)—(c)分别为Co2FeAl1–xSix系列合金薄膜在掺杂浓度x = 0.25,0.5和0.75下总磁矩及各原子磁矩随畸变度的变化. 由图可知, 当x = 0.25时, 当c/a < 1即对c轴进行压缩时, 总磁矩随着c/a的值的增大而增大, 由图可知Fe原子的磁矩变化决定了总磁矩的变化, 在1.0 < c/a <1.05的范围内, 总磁矩先减小后增大, Fe和Co原子共同决定了总磁矩的变化, 在c/a> 1.05时, 由图可知总磁矩的变化趋势主要决定于Co原子的磁矩变化. 当x = 0.5时, 由图可知磁矩的主要变化在畸变度c/a在1.0到1.1之间, 在c/a < 1.0时, 总磁矩主要趋势由Fe和Co两种原子磁矩变化决定, 当1.0 < c/a < 1.1时, 总磁矩先减小后增加, 主要变化趋势由Fe和Co两种原子变化趋势决定, c/a > 1.1时, 总磁矩略有减小, 变化趋势主要由Co原子变化趋势决定. 当x = 0.75时, 在c/a < 1.06范围内磁矩未产生明显变化, 在1.08 < c/a < 1.15之间, 磁矩先增加后减小, 磁矩的主要变化趋势决定于Co原子的主要变化趋势, c/a > 1.15时, 随着畸变度的增加, 总磁矩略有减小, 而且总磁矩变化趋势主要由Co原子磁矩变化决定.
-
本文运用基于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两种过渡金属原子决定.
-
运用基于密度泛函理论的第一性原理方法, 对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两种过渡金属原子的磁矩变化所决定.In the recent decades, the half-metallic materials have become a research hotspot because of their unique electronic structure. The 100% spin polarization at the Fermi level makes them widely used in spintronic devices. The Co-based Heusler alloys belong to an important class of magnetic material, and Co2FeAl and Co2FeSi have been experimentally confirmed to be half-metallic materials with 100% spin polarization at the Fermi level, and the Co2FeSi has a high Curie temperature of 1100 K and a large magnetic moment of 6.0
${{\text{μ}}{\rm{B}}}$ , which is a good candidate for spintronic devices. We here choose and substitute Al atoms in Co2FeAl with Si atoms, and then carry out the theoretical predictions of Co2FeAl1–xSix (x = 0.25, 0.5, 0.75) for both bulk and film . In this paper, using the first principles calculations based on the density functional theory (DFT) we study the electronic structure, tetragonal distortion, elastic constants, phonon spectrum and thermoelectric properties of Co2FeAl1–xSix (x = 0.25, 0.5, 0.75) series alloys. The calculation results show that the electronic structure of Co2FeAl1–xSix (x = 0.25, 0.5, 0.75) series alloys are all half-metallic with 100% spin polarization, and the down spin states (semiconducting character) all exhibit good thermoelectric properties, and the power factor increases with the substitution concentration of Si atoms increasing. The calculated phonon spectrum does not have virtual frequency, indicating its dynamic stability, and all cubic phases fulfill the mechanical stability criteria, i.e. Born criteria: C11 > 0, C44 > 0, C11–C12 > 0, C11 + 2C12 > 0, and C12 < B < C11. With the variation of lattice constant ratio c/a, the lowest energy point of the structure for Co2FeAl1–xSix (x = 0.25, 0.5, 0.75) series alloys are all at c/a = 1, showing that the stability of the structure does not change with the variation of distortion c/a, and further the martensitic transformation cannot occur. For the Co2FeAl1–xSix (x = 0.25, 0.5, 0.75) series alloy thin films, the calculated electronic structures all show a high spin polarization, and it reaches 100% at x = 0.75, and for x = 0.75, the lowest energy point of the structure is at c/a = 1.2, suggesting the martensitic transformation in this structure. With the variation of the tetragonal distortion, the total magnetic moment also changes and it is mainly determined by the changes of atomic magnetic moment of transition-metals Fe and Co.-
Keywords:
- half-metallic /
- first principles /
- electronic structure /
- magnetism
[1] Heusler F 1903 Deut. Phys. Ges. 5 219
[2] Murray S J, Marioni M, Allen S M, O’Handley R C 2000 Appl. Phys. Lett. 77 886
Google Scholar
[3] Donni A, Fischer P, Fauth F, Convert P, Aoki Y, Sugawara H, Sato H 1999 Physica B 259 705
[4] Wu G H, Yu C H, Meng L Q, Chen J L, Yang F M, Qi S R, Zhan W S 1999 Appl. Phys. Lett. 75 2990
Google Scholar
[5] Saha B, Shakouri A, Sands T D 2018 Appl. Phys. Rev. 5 021101
Google Scholar
[6] Webster P J 1971 J. Phys. Chem. Solids 32 1221
Google Scholar
[7] Kübler J, William A R, Sommers C B 1983 Phys. Rev. B 28 1745
Google Scholar
[8] de Groot R A, Müller F M, van Engen P G, Buschow K H J 1983 Phys. Rev. Lett. 50 2024
Google Scholar
[9] Comtesse D, Geisler B, Entel P, Kratzer P, Szunyogh L 2014 Phys. Rev. B 89 094410
Google Scholar
[10] Fecher G H, Felser C 2007 J. Phys. D: Appl. Phys. 40 1582
Google Scholar
[11] Li X M, Li T, Chen Z F, Hui F, Li X S, Wang X R, Xu J B, Zhu H W 2017 Appl. Phys. Rev. 4 021306
Google Scholar
[12] Balli M, Jandl S, Fournier P, Kedous-Lebouc A 2017 Appl. Phys. Rev. 4 021305
Google Scholar
[13] Kainuma R, Imano Y, Ito W, Sutou Y, Morito H, Okamoto S, Kitakami O, Oikawa K, Fujita A, Kanomata T, Ishida K 2006 Nature 439957
[14] Yu S Y, Liu Z H, Liu G D, Chen J L, Cao Z X, Wu G H, Zhang B, Zhang X X 2006 Appl. Phys. Lett. 89 162503
Google Scholar
[15] Dubenko I, Pathak A K, Stadler S, Ali N, Kovarskii Y, Prudnikov V N, Perov N S, Granovsky A B 2009 Phys. Rev. B 80 092408
Google Scholar
[16] Karaca H E, Karaman I, Basaran B, Ren Y, Chumlyakov Y I, Maier H J 2009 Adv. Funct. Mater. 19 983
Google Scholar
[17] Chmielus M, Zhang X X, Witherspoon C, Dunand D C, Mullner P 2009 Nat. Mater. 8 863
Google Scholar
[18] Sarawate N, Dapino M 2006 Appl. Phys. Lett. 88 121923
Google Scholar
[19] Mañosa L, González-Alonso D, Planes A, Bonnot E, Barrio M, Tamarit J L, Aksoy S, Acet M 2010 Nat. Mater. 9 478
Google Scholar
[20] Barman S R, Chakrabarti A, Singh S, Banik S, Bhardwaj S, Paulose P L, Chalke B A, Panda A K, Mitra A, Awasthi A M 2008 Phys. Rev. B 78 134406
Google Scholar
[21] Zayak A T, Entel P, Rabe K M, Adeagbo W A, Acet M 2005 Phys. Rev. B 72 054113
Google Scholar
[22] 罗礼进, 仲崇贵, 董正超, 方靖淮, 周朋霞, 江学范 2010 物理学报 59 8037
Google Scholar
Luo L J, Zhong C G, Dong Z C, Fang J H, Zhou P X, Jiang X F 2010 Acta Phys. Sin. 59 8037
Google Scholar
[23] 罗礼进, 仲崇贵, 江学范, 方靖淮, 蒋青 2010 物理学报 59 521
Google Scholar
Luo L J, Zhong C G, Jiang X F, Fang J H, Jiang Q 2010 Acta Phys. Sin. 59 521
Google Scholar
[24] 罗礼进, 仲崇贵, 赵永林, 方靖淮, 周朋霞, 江学范 2011 物理学报 60 127502
Google Scholar
Luo L J, Zhong C G, Zhao Y L, Fang J H, Zhou P X, Jiang X F 2011 Acta Phys. Sin. 60 127502
Google Scholar
[25] 罗礼进, 仲崇贵, 董正超, 方靖淮, 周朋霞, 江学范 2012 物理学报 61 207503
Google Scholar
Luo L J, Zhong C G, Dong Z C, Fang J H, Zhou P X, Jiang X F 2012 Acta Phys. Sin. 61 207503
Google Scholar
[26] Luo H Z, Jia P Z, Liu G D, Meng F B, Liu H Y, Liu E K, Wang W H, Wu G H, 2013 Solid State Commun. 17044
[27] Luo H Z, Meng F B, Liu G D, Liu H Y, Jia P Z, Liu E K, Wang W H, Wu G H 2013 Intermetallics 38 139
Google Scholar
[28] Kress G, Hafner J 1993 Phys. Rev. B 47 558
Google Scholar
[29] 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
[30] Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865
Google Scholar
[31] Kandpal H C, Fecher G H, Felser C 2007 J. Phys. D: Appl. Phys. 40 1507
Google Scholar
[32] Madsen G K H, Singh D J 2006 Comput. Phys. Commum. 175 67
Google Scholar
[33] Galanakis I, Mavropoulos P, Dederichs P H 2006 J. Phys. D: Appl. Phys. 39 765
Google Scholar
[34] Sargolzaei M, Richter M, Koepernik K, Opahle I, Eschrig H, Chaplygin I 2006 Phys. Rev. B 74 224410
Google Scholar
[35] Jansen H J F, Freeman A J 1984 Phys. Rev. B 30 561
Google Scholar
[36] Li J, Li J, Zhang Q, Zhang Z D, Yang G, Ma H R, Lu Z M, Fang W, Xie H X, Liang C Y, Yin F X 2016 Comp. Mater. Sci. 125 183
Google Scholar
[37] Li J, Yang G, Yang Y M, Ma H R, Zhang Q, Zhang Z D, Fang W, Yin F X, Li J 2017 J. Magn. Magn. Mater. 442 371
Google Scholar
[38] Kourov N I, Marchenkov V V, Perevozchikova Y A, Weber H W 2017 Phys. Solid State 59 898
Google Scholar
[39] Bilc D I, Mahanti S D, Kanatzidis M G 2006 Phys. Rev. B 74 125202
Google Scholar
[40] Al S, Arikan N, Demir S, Iyig€or A 2018 Physica B 531 16
Google Scholar
[41] Li J, Zhang Z D, Sun Y B, Zhang J, Zhou G X, Luo H Z, Liu G D 2013 Physica B 409 35
Google Scholar
[42] Okamura S, Miyazaki A, Sugimoto S 2005 Appl. Phys. Lett. 86 232503
Google Scholar
[43] Zhu W H, Wu D, Zhao B C, Zhu Z D, Yang X D, Zhang Z Z, Jin Q Y 2017 Phys. Rev. Appl. 8 034012
[44] Hazra B K, Raja M M, Srinath S 2016 J. Phys. D: Appl. Phys. 49 065007
Google Scholar
[45] Xu Z, Zhang Z, Hu F, Liu E, Xu F 2016 Mater. Res. Express 3 116103
Google Scholar
[46] Yadav A, Chaudhary S 2015 J. Appl. Phys. 118 193902
Google Scholar
[47] Chen J, Sakuraba Y, Masuda K, Miura Y, Li S, Kasai S, Furubayashi T, Hono K 2017 Appl. Phys. Lett. 110 242401
Google Scholar
[48] Huang X F, Dai Z W, Huang L, Lu G D, Liu M, Piao H G, Kim D H, Yu S C, Pan L Q 2016 J. Phys.: Condens. Matter 284 76006
-
图 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
$\mu $ for Co2FeAl1-xSix(x = 0.25, 0.5 and 0.75). The case of x = 0.25 corresponds to (a), (b), (c) and (d), and the case of x = 0.5 corresponds 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$\sigma $ , electronic thermal conductivity${\kappa _{\rm{e}}}$ and PF (${S^2}\sigma $ ), respectively.图 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.5 and 0.75) alloys in (a), (b) and (c). The temperature dependent heat capacity Cv with an inset graph showing the temperaturefrom 180 K to 250 K in (d).
图 8 (a) x = 0.25,(b) x = 0.5和(c) x = 0.75替代浓度下Co2FeAl1-xSix合金体相的总能量差
$\Delta E$ 与畸变度c/a的关系; (d) x = 0.25, (e) x = 0.5和(f) x = 0.75替代浓度下Co2FeAl1-xSix薄膜的驱动力$\Delta E$ 与畸变度c/a的关系Fig. 8. Calculated total energies as a function of the c/a ratio for Co2FeAl1-xSix (x = 0.25, 0.5 and 0.75) Heusler alloys in (a), (b) and (c) andfilm materials in (d), (e) and (f).
表 1 Co2FeAl1-xSix合金在x = 0.25, 0.5, 0.75时的晶格参数及磁矩
Table 1. Lattice parameters and magnetic moments of Co2FeAl1-xSix alloys at x = 0.25, 0.5 and 0.75.
a/Å mAl/${{\text{μ}}_{\rm{B}}}$ mSi/${{\text{μ}}_{\rm{B}}}$ mFe/${{\text{μ}}_{\rm{B}}}$ mCo/${{\text{μ}}_{\rm{B}}}$ Mt/${{\text{μ}}_{\rm{B}}}$ Co2FeAl0.75Si0.25 5.6520 –0.053 –0.039 2.998 1.262 5.473 Co2FeAl0.5Si0.5 5.6607 –0.046 –0.028 3.037 1.337 5.688 Co2FeAl0.25Si0.75 5.6406 –0.038 –0.012 3.094 1.400 5.891 表 2 计算的Co2FeAl1-xSix (x = 0.25, x = 0.5, x = 0.75)合金的弹性常数、体模量及剪切模量
Table 2. The calculated cubic elastic constant C11, C12, C44, shear modulus Gv, GR and GH in GPa.
C11/GPa C12/GPa C44/GPa B/GPa GV/GPa GR/GPa GH/GPa Co2FeAl0.75Si0.25 247.38 166.97 142.33 193.77 101.48 70.60 86.04 Co2FeAl0.5Si0.5 266.15 143.57 141.73 184.43 109.55 92.94 101.25 Co2FeAl0.25Si0.75 176.46 51.042 137.67 92.85 107.69 93.14 100.42 -
[1] Heusler F 1903 Deut. Phys. Ges. 5 219
[2] Murray S J, Marioni M, Allen S M, O’Handley R C 2000 Appl. Phys. Lett. 77 886
Google Scholar
[3] Donni A, Fischer P, Fauth F, Convert P, Aoki Y, Sugawara H, Sato H 1999 Physica B 259 705
[4] Wu G H, Yu C H, Meng L Q, Chen J L, Yang F M, Qi S R, Zhan W S 1999 Appl. Phys. Lett. 75 2990
Google Scholar
[5] Saha B, Shakouri A, Sands T D 2018 Appl. Phys. Rev. 5 021101
Google Scholar
[6] Webster P J 1971 J. Phys. Chem. Solids 32 1221
Google Scholar
[7] Kübler J, William A R, Sommers C B 1983 Phys. Rev. B 28 1745
Google Scholar
[8] de Groot R A, Müller F M, van Engen P G, Buschow K H J 1983 Phys. Rev. Lett. 50 2024
Google Scholar
[9] Comtesse D, Geisler B, Entel P, Kratzer P, Szunyogh L 2014 Phys. Rev. B 89 094410
Google Scholar
[10] Fecher G H, Felser C 2007 J. Phys. D: Appl. Phys. 40 1582
Google Scholar
[11] Li X M, Li T, Chen Z F, Hui F, Li X S, Wang X R, Xu J B, Zhu H W 2017 Appl. Phys. Rev. 4 021306
Google Scholar
[12] Balli M, Jandl S, Fournier P, Kedous-Lebouc A 2017 Appl. Phys. Rev. 4 021305
Google Scholar
[13] Kainuma R, Imano Y, Ito W, Sutou Y, Morito H, Okamoto S, Kitakami O, Oikawa K, Fujita A, Kanomata T, Ishida K 2006 Nature 439957
[14] Yu S Y, Liu Z H, Liu G D, Chen J L, Cao Z X, Wu G H, Zhang B, Zhang X X 2006 Appl. Phys. Lett. 89 162503
Google Scholar
[15] Dubenko I, Pathak A K, Stadler S, Ali N, Kovarskii Y, Prudnikov V N, Perov N S, Granovsky A B 2009 Phys. Rev. B 80 092408
Google Scholar
[16] Karaca H E, Karaman I, Basaran B, Ren Y, Chumlyakov Y I, Maier H J 2009 Adv. Funct. Mater. 19 983
Google Scholar
[17] Chmielus M, Zhang X X, Witherspoon C, Dunand D C, Mullner P 2009 Nat. Mater. 8 863
Google Scholar
[18] Sarawate N, Dapino M 2006 Appl. Phys. Lett. 88 121923
Google Scholar
[19] Mañosa L, González-Alonso D, Planes A, Bonnot E, Barrio M, Tamarit J L, Aksoy S, Acet M 2010 Nat. Mater. 9 478
Google Scholar
[20] Barman S R, Chakrabarti A, Singh S, Banik S, Bhardwaj S, Paulose P L, Chalke B A, Panda A K, Mitra A, Awasthi A M 2008 Phys. Rev. B 78 134406
Google Scholar
[21] Zayak A T, Entel P, Rabe K M, Adeagbo W A, Acet M 2005 Phys. Rev. B 72 054113
Google Scholar
[22] 罗礼进, 仲崇贵, 董正超, 方靖淮, 周朋霞, 江学范 2010 物理学报 59 8037
Google Scholar
Luo L J, Zhong C G, Dong Z C, Fang J H, Zhou P X, Jiang X F 2010 Acta Phys. Sin. 59 8037
Google Scholar
[23] 罗礼进, 仲崇贵, 江学范, 方靖淮, 蒋青 2010 物理学报 59 521
Google Scholar
Luo L J, Zhong C G, Jiang X F, Fang J H, Jiang Q 2010 Acta Phys. Sin. 59 521
Google Scholar
[24] 罗礼进, 仲崇贵, 赵永林, 方靖淮, 周朋霞, 江学范 2011 物理学报 60 127502
Google Scholar
Luo L J, Zhong C G, Zhao Y L, Fang J H, Zhou P X, Jiang X F 2011 Acta Phys. Sin. 60 127502
Google Scholar
[25] 罗礼进, 仲崇贵, 董正超, 方靖淮, 周朋霞, 江学范 2012 物理学报 61 207503
Google Scholar
Luo L J, Zhong C G, Dong Z C, Fang J H, Zhou P X, Jiang X F 2012 Acta Phys. Sin. 61 207503
Google Scholar
[26] Luo H Z, Jia P Z, Liu G D, Meng F B, Liu H Y, Liu E K, Wang W H, Wu G H, 2013 Solid State Commun. 17044
[27] Luo H Z, Meng F B, Liu G D, Liu H Y, Jia P Z, Liu E K, Wang W H, Wu G H 2013 Intermetallics 38 139
Google Scholar
[28] Kress G, Hafner J 1993 Phys. Rev. B 47 558
Google Scholar
[29] 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
[30] Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865
Google Scholar
[31] Kandpal H C, Fecher G H, Felser C 2007 J. Phys. D: Appl. Phys. 40 1507
Google Scholar
[32] Madsen G K H, Singh D J 2006 Comput. Phys. Commum. 175 67
Google Scholar
[33] Galanakis I, Mavropoulos P, Dederichs P H 2006 J. Phys. D: Appl. Phys. 39 765
Google Scholar
[34] Sargolzaei M, Richter M, Koepernik K, Opahle I, Eschrig H, Chaplygin I 2006 Phys. Rev. B 74 224410
Google Scholar
[35] Jansen H J F, Freeman A J 1984 Phys. Rev. B 30 561
Google Scholar
[36] Li J, Li J, Zhang Q, Zhang Z D, Yang G, Ma H R, Lu Z M, Fang W, Xie H X, Liang C Y, Yin F X 2016 Comp. Mater. Sci. 125 183
Google Scholar
[37] Li J, Yang G, Yang Y M, Ma H R, Zhang Q, Zhang Z D, Fang W, Yin F X, Li J 2017 J. Magn. Magn. Mater. 442 371
Google Scholar
[38] Kourov N I, Marchenkov V V, Perevozchikova Y A, Weber H W 2017 Phys. Solid State 59 898
Google Scholar
[39] Bilc D I, Mahanti S D, Kanatzidis M G 2006 Phys. Rev. B 74 125202
Google Scholar
[40] Al S, Arikan N, Demir S, Iyig€or A 2018 Physica B 531 16
Google Scholar
[41] Li J, Zhang Z D, Sun Y B, Zhang J, Zhou G X, Luo H Z, Liu G D 2013 Physica B 409 35
Google Scholar
[42] Okamura S, Miyazaki A, Sugimoto S 2005 Appl. Phys. Lett. 86 232503
Google Scholar
[43] Zhu W H, Wu D, Zhao B C, Zhu Z D, Yang X D, Zhang Z Z, Jin Q Y 2017 Phys. Rev. Appl. 8 034012
[44] Hazra B K, Raja M M, Srinath S 2016 J. Phys. D: Appl. Phys. 49 065007
Google Scholar
[45] Xu Z, Zhang Z, Hu F, Liu E, Xu F 2016 Mater. Res. Express 3 116103
Google Scholar
[46] Yadav A, Chaudhary S 2015 J. Appl. Phys. 118 193902
Google Scholar
[47] Chen J, Sakuraba Y, Masuda K, Miura Y, Li S, Kasai S, Furubayashi T, Hono K 2017 Appl. Phys. Lett. 110 242401
Google Scholar
[48] Huang X F, Dai Z W, Huang L, Lu G D, Liu M, Piao H G, Kim D H, Yu S C, Pan L Q 2016 J. Phys.: Condens. Matter 284 76006
计量
- 文章访问数: 3806
- PDF下载量: 70
- 被引次数: 0