搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Co2-基Heusler合金Co2FeAl1–xSix(x = 0.25, x = 0.5, x = 0.75)的结构、电子结构及热电特性的第一性原理研究

杨艳敏 李佳 马洪然 杨广 毛秀娟 李聪聪

Co2-基Heusler合金Co2FeAl1–xSix(x = 0.25, x = 0.5, x = 0.75)的结构、电子结构及热电特性的第一性原理研究

杨艳敏, 李佳, 马洪然, 杨广, 毛秀娟, 李聪聪
PDF
HTML
导出引用
  • 运用基于密度泛函理论的第一性原理方法, 对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两种过渡金属原子的磁矩变化所决定.
      通信作者: 李佳, jiali@hebut.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 61671199)、中国博士后基金(批准号: 61671199)、河北省博士后特别资助(批准号: 2016M601243)和国家春晖计划(批准号: Z2017024)资助的课题.
    [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

    [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

    [5]

    Saha B, Shakouri A, Sands T D 2018 Appl. Phys. Rev. 5 021101

    [6]

    Webster P J 1971 J. Phys. Chem. Solids 32 1221

    [7]

    Kübler J, William A R, Sommers C B 1983 Phys. Rev. B 28 1745

    [8]

    de Groot R A, Müller F M, van Engen P G, Buschow K H J 1983 Phys. Rev. Lett. 50 2024

    [9]

    Comtesse D, Geisler B, Entel P, Kratzer P, Szunyogh L 2014 Phys. Rev. B 89 094410

    [10]

    Fecher G H, Felser C 2007 J. Phys. D: Appl. Phys. 40 1582

    [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

    [12]

    Balli M, Jandl S, Fournier P, Kedous-Lebouc A 2017 Appl. Phys. Rev. 4 021305

    [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

    [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

    [16]

    Karaca H E, Karaman I, Basaran B, Ren Y, Chumlyakov Y I, Maier H J 2009 Adv. Funct. Mater. 19 983

    [17]

    Chmielus M, Zhang X X, Witherspoon C, Dunand D C, Mullner P 2009 Nat. Mater. 8 863

    [18]

    Sarawate N, Dapino M 2006 Appl. Phys. Lett. 88 121923

    [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

    [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

    [21]

    Zayak A T, Entel P, Rabe K M, Adeagbo W A, Acet M 2005 Phys. Rev. B 72 054113

    [22]

    罗礼进, 仲崇贵, 董正超, 方靖淮, 周朋霞, 江学范 2010 物理学报 59 8037

    Luo L J, Zhong C G, Dong Z C, Fang J H, Zhou P X, Jiang X F 2010 Acta Phys. Sin. 59 8037

    [23]

    罗礼进, 仲崇贵, 江学范, 方靖淮, 蒋青 2010 物理学报 59 521

    Luo L J, Zhong C G, Jiang X F, Fang J H, Jiang Q 2010 Acta Phys. Sin. 59 521

    [24]

    罗礼进, 仲崇贵, 赵永林, 方靖淮, 周朋霞, 江学范 2011 物理学报 60 127502

    Luo L J, Zhong C G, Zhao Y L, Fang J H, Zhou P X, Jiang X F 2011 Acta Phys. Sin. 60 127502

    [25]

    罗礼进, 仲崇贵, 董正超, 方靖淮, 周朋霞, 江学范 2012 物理学报 61 207503

    Luo L J, Zhong C G, Dong Z C, Fang J H, Zhou P X, Jiang X F 2012 Acta Phys. Sin. 61 207503

    [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

    [28]

    Kress G, Hafner J 1993 Phys. Rev. B 47 558

    [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

    [30]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865

    [31]

    Kandpal H C, Fecher G H, Felser C 2007 J. Phys. D: Appl. Phys. 40 1507

    [32]

    Madsen G K H, Singh D J 2006 Comput. Phys. Commum. 175 67

    [33]

    Galanakis I, Mavropoulos P, Dederichs P H 2006 J. Phys. D: Appl. Phys. 39 765

    [34]

    Sargolzaei M, Richter M, Koepernik K, Opahle I, Eschrig H, Chaplygin I 2006 Phys. Rev. B 74 224410

    [35]

    Jansen H J F, Freeman A J 1984 Phys. Rev. B 30 561

    [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

    [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

    [38]

    Kourov N I, Marchenkov V V, Perevozchikova Y A, Weber H W 2017 Phys. Solid State 59 898

    [39]

    Bilc D I, Mahanti S D, Kanatzidis M G 2006 Phys. Rev. B 74 125202

    [40]

    Al S, Arikan N, Demir S, Iyig€or A 2018 Physica B 531 16

    [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

    [42]

    Okamura S, Miyazaki A, Sugimoto S 2005 Appl. Phys. Lett. 86 232503

    [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

    [45]

    Xu Z, Zhang Z, Hu F, Liu E, Xu F 2016 Mater. Res. Express 3 116103

    [46]

    Yadav A, Chaudhary S 2015 J. Appl. Phys. 118 193902

    [47]

    Chen J, Sakuraba Y, Masuda K, Miura Y, Li S, Kasai S, Furubayashi T, Hono K 2017 Appl. Phys. Lett. 110 242401

    [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

  • 图 1  (a) Co2FeAl1xSix (x = 0.25)的L21结构; (b) Co2FeAl0.75Si0.25的薄膜结构

    Fig. 1.  (a) L21 structure of Co2FeAl1-xSix (x = 0.25); (b)thin film structure of Co2FeAl0.75Si0.25.

    图 2  Co2FeAl1-xSix合金在铁磁态(FM)和反铁磁态(AFM)下的晶格常数优化曲线 (a) x = 0.25; (b) x = 0.5; (c) x = 0.75

    Fig. 2.  Optimization curves of lattice constant for Co2FeAl1-xSix alloy under ferromagnetic and antiferromagnetic magnetic order.

    图 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.5 and (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).

    图 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  Co2FeAl1xSix合金在x = 0.25, 0.5, 0.75时的声子谱及比热容 (a) Co2FeAl1xSix (x = 0.25), (b) Co2FeAl1xSi x (x = 0.5)和(c) Co2FeAl1xSi x (x = 0.75)的声子谱; (d) Co2FeAl1 xSix (x = 0.25, 0.5, 0.75)的比热容随温度的变化

    Fig. 6.  Full phonon spectra of Co2FeAl1xSix (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).

    图 7  (a) Co2FeAl0.75Si0.25, (b) Co2FeAl0.5Si0.5和(c) Co2FeAl0.25Si0.75薄膜的总态密度和原子分态密度

    Fig. 7.  Thetotaland atom-projected density of states for Co2FeAl1xSix (x = 0.25, 0.5 and 0.75) film in (a), (b) and (c).

    图 8  (a) x = 0.25,(b) x = 0.5和(c) x = 0.75替代浓度下Co2FeAl1xSix合金体相的总能量差$\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 Co2FeAl1xSix (x = 0.25, 0.5 and 0.75) Heusler alloys in (a), (b) and (c) andfilm materials in (d), (e) and (f).

    图 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-xSix film change with distortion at x = 0.25, x = 0.5 and x = 0.75 in (a), (b) and (c).

    表 1  Co2FeAl1xSix合金在x = 0.25, 0.5, 0.75时的晶格参数及磁矩

    Table 1.  Lattice parameters and magnetic moments of Co2FeAl1xSix alloys at x = 0.25, 0.5 and 0.75.

    amAl/${{\text{μ}}_{\rm{B}}}$mSi/${{\text{μ}}_{\rm{B}}}$mFe/${{\text{μ}}_{\rm{B}}}$mCo/${{\text{μ}}_{\rm{B}}}$Mt/${{\text{μ}}_{\rm{B}}}$
    Co2FeAl0.75Si0.255.6520–0.053–0.0392.9981.2625.473
    Co2FeAl0.5Si0.55.6607–0.046–0.0283.0371.3375.688
    Co2FeAl0.25Si0.755.6406–0.038–0.0123.0941.4005.891
    下载: 导出CSV

    表 2  计算的Co2FeAl1xSix (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/GPaC12/GPaC44/GPaB/GPaGV/GPaGR/GPaGH/GPa
    Co2FeAl0.75Si0.25247.38166.97142.33193.77101.4870.6086.04
    Co2FeAl0.5Si0.5266.15143.57141.73184.43109.5592.94101.25
    Co2FeAl0.25Si0.75176.4651.042137.6792.85107.6993.14100.42
    下载: 导出CSV
  • [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

    [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

    [5]

    Saha B, Shakouri A, Sands T D 2018 Appl. Phys. Rev. 5 021101

    [6]

    Webster P J 1971 J. Phys. Chem. Solids 32 1221

    [7]

    Kübler J, William A R, Sommers C B 1983 Phys. Rev. B 28 1745

    [8]

    de Groot R A, Müller F M, van Engen P G, Buschow K H J 1983 Phys. Rev. Lett. 50 2024

    [9]

    Comtesse D, Geisler B, Entel P, Kratzer P, Szunyogh L 2014 Phys. Rev. B 89 094410

    [10]

    Fecher G H, Felser C 2007 J. Phys. D: Appl. Phys. 40 1582

    [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

    [12]

    Balli M, Jandl S, Fournier P, Kedous-Lebouc A 2017 Appl. Phys. Rev. 4 021305

    [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

    [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

    [16]

    Karaca H E, Karaman I, Basaran B, Ren Y, Chumlyakov Y I, Maier H J 2009 Adv. Funct. Mater. 19 983

    [17]

    Chmielus M, Zhang X X, Witherspoon C, Dunand D C, Mullner P 2009 Nat. Mater. 8 863

    [18]

    Sarawate N, Dapino M 2006 Appl. Phys. Lett. 88 121923

    [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

    [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

    [21]

    Zayak A T, Entel P, Rabe K M, Adeagbo W A, Acet M 2005 Phys. Rev. B 72 054113

    [22]

    罗礼进, 仲崇贵, 董正超, 方靖淮, 周朋霞, 江学范 2010 物理学报 59 8037

    Luo L J, Zhong C G, Dong Z C, Fang J H, Zhou P X, Jiang X F 2010 Acta Phys. Sin. 59 8037

    [23]

    罗礼进, 仲崇贵, 江学范, 方靖淮, 蒋青 2010 物理学报 59 521

    Luo L J, Zhong C G, Jiang X F, Fang J H, Jiang Q 2010 Acta Phys. Sin. 59 521

    [24]

    罗礼进, 仲崇贵, 赵永林, 方靖淮, 周朋霞, 江学范 2011 物理学报 60 127502

    Luo L J, Zhong C G, Zhao Y L, Fang J H, Zhou P X, Jiang X F 2011 Acta Phys. Sin. 60 127502

    [25]

    罗礼进, 仲崇贵, 董正超, 方靖淮, 周朋霞, 江学范 2012 物理学报 61 207503

    Luo L J, Zhong C G, Dong Z C, Fang J H, Zhou P X, Jiang X F 2012 Acta Phys. Sin. 61 207503

    [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

    [28]

    Kress G, Hafner J 1993 Phys. Rev. B 47 558

    [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

    [30]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865

    [31]

    Kandpal H C, Fecher G H, Felser C 2007 J. Phys. D: Appl. Phys. 40 1507

    [32]

    Madsen G K H, Singh D J 2006 Comput. Phys. Commum. 175 67

    [33]

    Galanakis I, Mavropoulos P, Dederichs P H 2006 J. Phys. D: Appl. Phys. 39 765

    [34]

    Sargolzaei M, Richter M, Koepernik K, Opahle I, Eschrig H, Chaplygin I 2006 Phys. Rev. B 74 224410

    [35]

    Jansen H J F, Freeman A J 1984 Phys. Rev. B 30 561

    [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

    [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

    [38]

    Kourov N I, Marchenkov V V, Perevozchikova Y A, Weber H W 2017 Phys. Solid State 59 898

    [39]

    Bilc D I, Mahanti S D, Kanatzidis M G 2006 Phys. Rev. B 74 125202

    [40]

    Al S, Arikan N, Demir S, Iyig€or A 2018 Physica B 531 16

    [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

    [42]

    Okamura S, Miyazaki A, Sugimoto S 2005 Appl. Phys. Lett. 86 232503

    [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

    [45]

    Xu Z, Zhang Z, Hu F, Liu E, Xu F 2016 Mater. Res. Express 3 116103

    [46]

    Yadav A, Chaudhary S 2015 J. Appl. Phys. 118 193902

    [47]

    Chen J, Sakuraba Y, Masuda K, Miura Y, Li S, Kasai S, Furubayashi T, Hono K 2017 Appl. Phys. Lett. 110 242401

    [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

  • [1] 姚仲瑜, 孙丽, 潘孟美, 孙书娟. 第一性原理研究semi-Heusler合金CoCrTe和CoCrSb的半金属性和磁性. 物理学报, 2016, 65(12): 127501. doi: 10.7498/aps.65.127501
    [2] 姚仲瑜, 孙丽, 潘孟美, 孙书娟, 刘汉军. 第一性原理研究half-Heusler合金VLiBi和CrLiBi的半金属铁磁性. 物理学报, 2018, 67(21): 217501. doi: 10.7498/aps.67.20181129
    [3] 许佳玲, 贾利云, 靳晓庆, 郝兴楠, 马丽, 侯登录. 系列CoMnZnZ四元Heusler化合物的结构和半金属铁磁性. 物理学报, 2019, 68(15): 157501. doi: 10.7498/aps.68.20190207
    [4] 文黎巍, 王玉梅, 裴慧霞, 丁俊. Sb系half-Heusler合金磁性及电子结构的第一性原理研究. 物理学报, 2011, 60(4): 047110. doi: 10.7498/aps.60.047110
    [5] 赵晶晶, 祁欣, 刘恩克, 朱伟, 钱金凤, 李贵江, 王文洪, 吴光恒. Co50Fe25-xMnxSi25系列合金的结构、磁性和半金属性研究. 物理学报, 2011, 60(4): 047108. doi: 10.7498/aps.60.047108
    [6] 程志梅, 王新强, 王风, 鲁丽娅, 刘高斌, 段壮芬, 聂招秀. 三元化合物ZnCrS2电子结构和半金属铁磁性的第一性原理研究. 物理学报, 2011, 60(9): 096301. doi: 10.7498/aps.60.096301
    [7] 王少霞, 赵旭才, 潘多桥, 庞国旺, 刘晨曦, 史蕾倩, 刘桂安, 雷博程, 黄以能, 张丽丽. 过渡金属(Cr, Mn, Fe, Co)掺杂对TiO2磁性影响的第一性原理研究. 物理学报, 2020, 69(19): 197101. doi: 10.7498/aps.69.20200644
    [8] 罗礼进, 仲崇贵, 方靖淮, 赵永林, 周朋霞, 江学范. Heusler合金Mn2 NiAl的电子结构和磁性对四方畸变的响应及其压力响应. 物理学报, 2011, 60(12): 127502. doi: 10.7498/aps.60.127502
    [9] 杜音, 王文洪, 张小明, 刘恩克, 吴光恒. 铁基Heusler合金Fe2Co1-xCrxSi的结构、磁性和输运性质的研究. 物理学报, 2012, 61(14): 147304. doi: 10.7498/aps.61.147304
    [10] 黄有林, 侯育花, 赵宇军, 刘仲武, 曾德长, 马胜灿. 应变对钴铁氧体电子结构和磁性能影响的第一性原理研究. 物理学报, 2013, 62(16): 167502. doi: 10.7498/aps.62.167502
    [11] 姚仲瑜, 傅军, 龚少华, 张月胜, 姚凯伦. 晶格各向同性应变对闪锌矿结构CrS和CrSe的半金属性和磁性的影响. 物理学报, 2011, 60(12): 127103. doi: 10.7498/aps.60.127103
    [12] 段满益, 徐 明, 周海平, 沈益斌, 陈青云, 丁迎春, 祝文军. 过渡金属与氮共掺杂ZnO电子结构和光学性质的第一性原理研究. 物理学报, 2007, 56(9): 5359-5365. doi: 10.7498/aps.56.5359
    [13] 罗礼进, 仲崇贵, 江学范, 方靖淮, 蒋青. Heusler合金Ni2MnSi的电子结构、磁性、压力响应及四方变形的第一性原理研究. 物理学报, 2010, 59(1): 521-526. doi: 10.7498/aps.59.521
    [14] 颜送灵, 唐黎明, 赵宇清. 不同组分厚度比的LaMnO3/SrTiO3异质界面电子结构和磁性的第一性原理研究. 物理学报, 2016, 65(7): 077301. doi: 10.7498/aps.65.077301
    [15] 马振宁, 蒋敏, 王磊. Mg-Y-Zn合金三元金属间化合物的电子结构及其相稳定性的第一性原理研究. 物理学报, 2015, 64(18): 187102. doi: 10.7498/aps.64.187102
    [16] 张加宏, 马 荣, 刘 甦, 刘 楣. 掺杂MgCNi3超导电性和磁性的第一性原理研究. 物理学报, 2006, 55(9): 4816-4821. doi: 10.7498/aps.55.4816
    [17] 张富春, 张威虎, 董军堂, 张志勇. Cr掺杂ZnO纳米线的电子结构和磁性. 物理学报, 2011, 60(12): 127503. doi: 10.7498/aps.60.127503
    [18] 胡洁琼, 谢明, 张吉明, 刘满门, 杨有才, 陈永泰. Au-Sn金属间化合物的第一性原理研究. 物理学报, 2013, 62(24): 247102. doi: 10.7498/aps.62.247102
    [19] 杨彪, 王丽阁, 易勇, 王恩泽, 彭丽霞. C, N, O原子在金属V中扩散行为的第一性原理计算. 物理学报, 2015, 64(2): 026602. doi: 10.7498/aps.64.026602
    [20] 陈国祥, 樊晓波, 李思琦, 张建民. 碱金属和碱土金属掺杂二维GaN材料电磁特性的第一性原理计算. 物理学报, 2019, 68(23): 237303. doi: 10.7498/aps.68.20191246
  • 引用本文:
    Citation:
计量
  • 文章访问数:  1182
  • PDF下载量:  38
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-09-03
  • 修回日期:  2018-12-24
  • 上网日期:  2019-02-01
  • 刊出日期:  2019-02-20

Co2-基Heusler合金Co2FeAl1–xSix(x = 0.25, x = 0.5, x = 0.75)的结构、电子结构及热电特性的第一性原理研究

    基金项目: 国家自然科学基金(批准号: 61671199)、中国博士后基金(批准号: 61671199)、河北省博士后特别资助(批准号: 2016M601243)和国家春晖计划(批准号: Z2017024)资助的课题.

摘要: 运用基于密度泛函理论的第一性原理方法, 对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两种过渡金属原子的磁矩变化所决定.

English Abstract

    • 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, 其中XY为过渡金属元素, 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) Co2FeAl1xSix (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原子影响.

      amAl/${{\text{μ}}_{\rm{B}}}$mSi/${{\text{μ}}_{\rm{B}}}$mFe/${{\text{μ}}_{\rm{B}}}$mCo/${{\text{μ}}_{\rm{B}}}$Mt/${{\text{μ}}_{\rm{B}}}$
      Co2FeAl0.75Si0.255.6520–0.053–0.0392.9981.2625.473
      Co2FeAl0.5Si0.55.6607–0.046–0.0283.0371.3375.688
      Co2FeAl0.25Si0.755.6406–0.038–0.0123.0941.4005.891

      表 1  Co2FeAl1xSix合金在x = 0.25, 0.5, 0.75时的晶格参数及磁矩

      Table 1.  Lattice parameters and magnetic moments of Co2FeAl1xSix alloys at x = 0.25, 0.5 and 0.75.

      图  2  Co2FeAl1-xSix合金在铁磁态(FM)和反铁磁态(AFM)下的晶格常数优化曲线 (a) x = 0.25; (b) x = 0.5; (c) x = 0.75

      Figure 2.  Optimization curves of lattice constant for Co2FeAl1-xSix alloy under ferromagnetic and antiferromagnetic magnetic order.

    • 我们运用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}}$%, 其中NN分别代表自旋向上和向下态在费米面处的态密度[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电子劈裂以及它们之间的杂化是形成带隙的主要原因, 并决定了带隙的宽度.

      图  4  (a) Co2FeAl0.75Si0.25,(b) Co2FeAl0.5Si0.5和(c) Co2FeAl0.25Si0.75的总态密度和分态密度

      Figure 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).

    • 热电材料种类很多, 如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  Co2FeAl1xSix合金在x = 0.25, 0.5, 0.75时的声子谱及比热容 (a) Co2FeAl1xSix (x = 0.25), (b) Co2FeAl1xSi x (x = 0.5)和(c) Co2FeAl1xSi x (x = 0.75)的声子谱; (d) Co2FeAl1 xSix (x = 0.25, 0.5, 0.75)的比热容随温度的变化

      Figure 6.  Full phonon spectra of Co2FeAl1xSix (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, C12C44. 一个稳定的立方相应该满足的机械稳定性标准, 即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/GPaC12/GPaC44/GPaB/GPaGV/GPaGR/GPaGH/GPa
      Co2FeAl0.75Si0.25247.38166.97142.33193.77101.4870.6086.04
      Co2FeAl0.5Si0.5266.15143.57141.73184.43109.5592.94101.25
      Co2FeAl0.25Si0.75176.4651.042137.6792.85107.6993.14100.42

      表 2  计算的Co2FeAl1xSix (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原子的影响.

      图  7  (a) Co2FeAl0.75Si0.25, (b) Co2FeAl0.5Si0.5和(c) Co2FeAl0.25Si0.75薄膜的总态密度和原子分态密度

      Figure 7.  Thetotaland atom-projected density of states for Co2FeAl1xSix (x = 0.25, 0.5 and 0.75) film in (a), (b) and (c).

    • 近年来, 铁磁形记忆合金(FSMAS)在磁场的作用下表现出很好的物理性质而备受关注, 目前已经报道过的Heusler型FSMAs主要有Ni-基, Mn-基, Fe-基和Co-基体系. 基于前人的结论可知, Heusler合金体系要发生马氏体相变需要满足一定的条件: 首先要有充足的相变驱动力即四方相与立方相之间的能量差值$\Delta E$; 其次, 畸变度, 即晶格常数ca的比值也是关键, 满足马氏体相变的适当的畸变度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替代浓度下Co2FeAl1xSix合金体相的总能量差$\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 Co2FeAl1xSix (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原子磁矩变化决定.

      图  9  (a) x = 0.25,(b) x = 0.5和(c) x = 0.75替代浓度下Co2FeAl1-xSix合金薄膜的总磁矩及各原子总磁矩随畸变度的变化

      Figure 9.  The total magnetic moment and the magnetic moment of each atom of Co2FeAl1-xSix film change with distortion at x = 0.25, x = 0.5 and x = 0.75 in (a), (b) and (c).

    • 本文运用基于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两种过渡金属原子决定.

参考文献 (48)

目录

    /

    返回文章
    返回