搜索

x

留言板

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

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

Be, Mg, Mn掺杂CuInO2形成能的第一性原理研究

莫曼 曾纪术 何浩 张喨 杜龙 方志杰

Be, Mg, Mn掺杂CuInO2形成能的第一性原理研究

莫曼, 曾纪术, 何浩, 张喨, 杜龙, 方志杰
PDF
HTML
导出引用
导出核心图
  • 研制开发新型的光电材料对促进社会经济发展具有重要的科学意义和实用价值. 利用宽禁带CuInO2铟基材料实现全透明光电材料是目前深入研究的热点. 通过基于密度泛函的第一性原理计算方法, 本文计算出掺杂元素Mg, Be, Mn在CuInO2的形成能. 计算结果表明, 施主类缺陷(如掺杂元素替代Cu原子或进入间隙位置)由于较高的形成能和较深的跃迁能级, 很难在CuInO2材料中出现N型导电; 而受主缺陷中, 在氧原子化学势极大的情况下, Mg原子替代In能成为CuInO2理想的受主缺陷. 计算结果可为制备性能优异的CuInO2材料提供指导.
      通信作者: 方志杰, nnfang@semi.ac.cn
    • 基金项目: 国家自然科学基金(批准号: 11464003, 11864005)、广西自然科学基金(批准号: 2017GXNSFAA198315)、柳州市科技计划项目(批准号: 2016B040202)和广西高校中青年教师基础能力提升项目(批准号: 2018KY0324)资助的课题.
    [1]

    Ginley D S, Bright C 2000 MRS Bull. 25 8

    [2]

    Wang L J, Wang W Z, Chen Y L, Yao L Z, Zhao X, Shi H L, Cao M S, Liang Y J 2018 ACS Appl. Mat. Interfaces 10 11652

    [3]

    Cao M S, Wang X X, Cao W Q, Fang X Y, Wen B, Yuan J 2018 Small 14 1

    [4]

    Chen Y L, Wang L Y, Wang W Z, Cao M S 2017 Appl. Catal. B 209 110

    [5]

    Kawazoe H, Yasukawa M, Hyodo H 1997 Nature 389 939

    [6]

    Yanagi H, Inoue S, Ueda K, Kawazoe H, Hosono H 2000 J. Appl. Phys. 88 4159

    [7]

    Nakanishi A, Katayama-Yoshida H, Ishikawa T, Shimizu K 2016 J. Phys. Soc. Jpn. 85 094711

    [8]

    Jedidi A, Rasul S, Masih D, Cavallo L, Takanabe K 2015 J. Mater. Chem. A 3 19085

    [9]

    Nie X, Wei S H, Zhang S B 2002 Phys. Rev. Lett. 88 066405

    [10]

    Hamada I, Katayama-Yoshida H 2006 Physsica B 377 808

    [11]

    Jiang H F, Zhu X B, Lei H C, Li G, Yang Z R 2011 Thin Solids Film 519 2559

    [12]

    Shimode M, Sasaki M, Mukaida K 2000 J. Soli. Stat. Chem. 151 16

    [13]

    Liu Q J, Liu Z T, Feng L P 2010 Physica B 405 2028

    [14]

    Singh M, Mehta B R, Varandani D, Singh V N 2009 J. Appl. Phys. 106 053709

    [15]

    Ye F, Cai X M, Dai F P, Zhang D P, Fan P, Liu L J 2011 Adv. Mater. Res. 239 242

    [16]

    Shin D, Foord J S, Payne D J 2009 Phys. Rev. B 80 233105

    [17]

    Varandani D, Singh B, Mehta B, Singh M, Singh V, Nand G D 2010 J. Appl. Phys. 107 103703

    [18]

    Roland G, John R 2011 Phys. Rev. B 84 035125

    [19]

    Godinho K G, Morgan B J, Allen J P, Scanlon D O, Watson G W 2011 J. Phys.: Cond. Matter 23 334201

    [20]

    Yao Y, Xie G, Song N, Yu X H, Li R X 2011 Adv. Mater. Res. 399 401

    [21]

    Liu L, Bai K W, Gong H, Wu P 2005 Phys. Rev. B 72 125204

    [22]

    Blouchl P E 1994 Phys. Rev. B 50 17953

    [23]

    Kresse G, Joubert J 1999 Phys. Rev. B 59 1758

    [24]

    Hohenberg P, Kohn W 1964 Phys. Rev. B 136 864

    [25]

    Kohn W, Sham L J 1965 Phys. Rev. A 140 1133

    [26]

    Kresse G, Furthmüller J 1996 Phys. Rev. B 54 11169

    [27]

    Kresse G, Furthmüller J 1996 Comp. Mater. Sci. 6 15

    [28]

    Monkhorst H J, Pack J D 1976 Phys. Rev. B 13 5188

    [29]

    Pack J D, Monkhorst H J 1977 Phys. Rev. B 16 1748

    [30]

    Murnaghan F D, Natl P 1944 Acad. Sci. USA 30 244

    [31]

    Chu S, Hollberg L, Bjorkholm J E, Cable A, Ashkin A 1985 Phys. Rev. Lett. 55 48

    [32]

    Wei S H 2004 Comput. Mat. Sci. 30 337

  • 图 1  CuInO2晶体结构图, 图中红色原子为O原子, 灰色原子为In原子, 棕色原子为Cu原子 (a)黄色掺杂原子替代Cu原子的情况; (b)绿色掺杂原子替代In原子的情况

    Fig. 1.  The crystal structure of the CuInO2, the red atoms are O atoms, the brown atoms are Cu atoms, the purple atoms are In atoms: (a) Substituting yellow dopant atom for Cu atom; (b) substituting green dopant atom for In atom.

    图 2  掺杂形成能在cation-poor, anion-rich的变化图

    Fig. 2.  The change of doping formation energies under cation-poor, anion-rich condition.

    图 3  掺杂形成能在cation-rich, anion-poor的变化图

    Fig. 3.  The change of doping formation energies under cation-rich, anion-poor condition.

    图 4  掺杂元素在CuInO2的缺陷跃迁能级

    Fig. 4.  The calculated transition energy levels of the extrinsic defects in CuInO2.

    表 1  CuInO2晶格常数理论计算值与实验值

    Table 1.  Theoretical values and experimental values of lattice constants in CuInO2.

    CuInO2 晶格常数理论计算值 晶格常数实验值
    a 3.297 Å a 3.292 Å
    c 17.192 Å c 17.388 Å
    u 0.1047 Å u 0.1061 Å
    c/a 5.214 c/a 5.282
    下载: 导出CSV

    表 2  掺杂元素周围的弛豫变化情况, 键距后括号为掺杂原子的近邻原子种类和个数

    Table 2.  The surrounding atoms of dopant, kinds and numbers of dopantxs nearest neighbor atoms in parentheses.

    掺杂种类 最近邻键距 掺杂种类 最近邻键距 掺杂种类 最近邻键距
    BeCu 1.49 Å(O, 2) MgCu 1.85 Å(O, 2) MnCu 1.77Å(O, 2)
    BeIn 1.96 Å(O, 6) MgIn 2.12 Å(O, 6) MnIn 1.93Å(O, 6)
    Bei 1.62 Å(O, 3) Mgi 2.04 Å(O, 3) Mni 1.99Å(O, 3)
    下载: 导出CSV

    表 3  CuInO2的掺杂元素形成能

    Table 3.  The calculated formation energies of dopants in CuInO2.

    掺杂类型 形成能 掺杂类型 总能/eV 形成能 掺杂类型 总能/eV 形成能
    BeCu –2.20 MgCu –650.87 –2.58 MnCu –654.63 1.76
    BeIn 0.79 MgIn –651.36 –1.59 MnIn –656.54 1.33
    Bei –1.05 Mgi –654.27 –1.27 Mni –657.39 3.71
    下载: 导出CSV
  • [1]

    Ginley D S, Bright C 2000 MRS Bull. 25 8

    [2]

    Wang L J, Wang W Z, Chen Y L, Yao L Z, Zhao X, Shi H L, Cao M S, Liang Y J 2018 ACS Appl. Mat. Interfaces 10 11652

    [3]

    Cao M S, Wang X X, Cao W Q, Fang X Y, Wen B, Yuan J 2018 Small 14 1

    [4]

    Chen Y L, Wang L Y, Wang W Z, Cao M S 2017 Appl. Catal. B 209 110

    [5]

    Kawazoe H, Yasukawa M, Hyodo H 1997 Nature 389 939

    [6]

    Yanagi H, Inoue S, Ueda K, Kawazoe H, Hosono H 2000 J. Appl. Phys. 88 4159

    [7]

    Nakanishi A, Katayama-Yoshida H, Ishikawa T, Shimizu K 2016 J. Phys. Soc. Jpn. 85 094711

    [8]

    Jedidi A, Rasul S, Masih D, Cavallo L, Takanabe K 2015 J. Mater. Chem. A 3 19085

    [9]

    Nie X, Wei S H, Zhang S B 2002 Phys. Rev. Lett. 88 066405

    [10]

    Hamada I, Katayama-Yoshida H 2006 Physsica B 377 808

    [11]

    Jiang H F, Zhu X B, Lei H C, Li G, Yang Z R 2011 Thin Solids Film 519 2559

    [12]

    Shimode M, Sasaki M, Mukaida K 2000 J. Soli. Stat. Chem. 151 16

    [13]

    Liu Q J, Liu Z T, Feng L P 2010 Physica B 405 2028

    [14]

    Singh M, Mehta B R, Varandani D, Singh V N 2009 J. Appl. Phys. 106 053709

    [15]

    Ye F, Cai X M, Dai F P, Zhang D P, Fan P, Liu L J 2011 Adv. Mater. Res. 239 242

    [16]

    Shin D, Foord J S, Payne D J 2009 Phys. Rev. B 80 233105

    [17]

    Varandani D, Singh B, Mehta B, Singh M, Singh V, Nand G D 2010 J. Appl. Phys. 107 103703

    [18]

    Roland G, John R 2011 Phys. Rev. B 84 035125

    [19]

    Godinho K G, Morgan B J, Allen J P, Scanlon D O, Watson G W 2011 J. Phys.: Cond. Matter 23 334201

    [20]

    Yao Y, Xie G, Song N, Yu X H, Li R X 2011 Adv. Mater. Res. 399 401

    [21]

    Liu L, Bai K W, Gong H, Wu P 2005 Phys. Rev. B 72 125204

    [22]

    Blouchl P E 1994 Phys. Rev. B 50 17953

    [23]

    Kresse G, Joubert J 1999 Phys. Rev. B 59 1758

    [24]

    Hohenberg P, Kohn W 1964 Phys. Rev. B 136 864

    [25]

    Kohn W, Sham L J 1965 Phys. Rev. A 140 1133

    [26]

    Kresse G, Furthmüller J 1996 Phys. Rev. B 54 11169

    [27]

    Kresse G, Furthmüller J 1996 Comp. Mater. Sci. 6 15

    [28]

    Monkhorst H J, Pack J D 1976 Phys. Rev. B 13 5188

    [29]

    Pack J D, Monkhorst H J 1977 Phys. Rev. B 16 1748

    [30]

    Murnaghan F D, Natl P 1944 Acad. Sci. USA 30 244

    [31]

    Chu S, Hollberg L, Bjorkholm J E, Cable A, Ashkin A 1985 Phys. Rev. Lett. 55 48

    [32]

    Wei S H 2004 Comput. Mat. Sci. 30 337

  • [1] 陈东运, 高明, 李拥华, 徐飞, 赵磊, 马忠权. MoO3/Si界面区钼掺杂非晶氧化硅层形成的第一性原理研究. 物理学报, 2019, 68(10): 103101. doi: 10.7498/aps.68.20190067
    [2] 李宗宝, 王霞, 樊帅伟. Cu/N表面沉积共掺杂TiO2光催化剂作用机理的理论研究. 物理学报, 2014, 63(15): 157102. doi: 10.7498/aps.63.157102
    [3] 李虹, 王绍青, 叶恒强. Nb掺杂对γ-TiAl抗氧化能力影响的第一性原理研究. 物理学报, 2009, 58(13): 224-S229. doi: 10.7498/aps.58.224
    [4] 王泽普, 付念, 于涵, 徐晶威, 何祺, 郑树凯, 丁帮福, 闫小兵. 铟掺杂钨位增强钨酸铋氧空位光催化效率. 物理学报, 2019, 68(21): 217102. doi: 10.7498/aps.68.20191010
    [5] 唐冬华, 薛林, 孙立忠, 钟建新. B在Hg0.75Cd0.25Te中掺杂效应的第一性原理研究. 物理学报, 2012, 61(2): 027102. doi: 10.7498/aps.61.027102
    [6] 张梅玲, 陈玉红, 张材荣, 李公平. 内在缺陷与Cu掺杂共存对ZnO电磁光学性质影响的第一性原理研究. 物理学报, 2019, 68(8): 087101. doi: 10.7498/aps.68.20182238
    [7] 郝红飞, 王静, 孙锋, 张澜庭. Dy在Nd2Fe14B晶格中的占位及其对Fe原子磁矩影响的第一性原理计算. 物理学报, 2013, 62(11): 117501. doi: 10.7498/aps.62.117501
    [8] 徐金荣, 王影, 朱兴凤, 李平, 张莉. N掺杂和N-V共掺杂锐钛矿相TiO2的第一性原理研究. 物理学报, 2012, 61(20): 207103. doi: 10.7498/aps.61.207103
    [9] 黄祖飞, 明 星, 王春忠, 孟 醒, 陈 岗, 金胜哲. 二价金属元素掺杂对LiCoO2体系电子输运性质的影响. 物理学报, 2007, 56(10): 6008-6012. doi: 10.7498/aps.56.6008
    [10] 邹林儿, 矶 守, 浜中广见, 杜丽萍, 陈抱雪, 孙 蓓, 陈 直. 掺杂As2S8非晶态薄膜波导的光阻断效应. 物理学报, 2008, 57(6): 3593-3599. doi: 10.7498/aps.57.3593
    [11] 周传仓, 刘发民, 丁芃, 钟文武, 蔡鲁刚, 曾乐贵. 钶铁矿型MnNb2O6的熔盐法合成、钒掺杂与磁性研究. 物理学报, 2011, 60(4): 048101. doi: 10.7498/aps.60.048101
    [12] 吴木生, 徐波, 刘刚, 欧阳楚英. Cr和W掺杂的单层MoS2电子结构的第一性原理研究. 物理学报, 2013, 62(3): 037103. doi: 10.7498/aps.62.037103
    [13] 曹娟, 崔磊, 潘靖. V,Cr,Mn掺杂MoS2磁性的第一性原理研究. 物理学报, 2013, 62(18): 187102. doi: 10.7498/aps.62.187102
    [14] 徐晶, 梁家青, 李红萍, 李长生, 刘孝娟, 孟健. Ti掺杂NbSe2电子结构的第一性原理研究. 物理学报, 2015, 64(20): 207101. doi: 10.7498/aps.64.207101
    [15] 袁国才, 陈曦, 黄雨阳, 毛俊西, 禹劲秋, 雷晓波, 张勤勇. Mg2Si0.3Sn0.7掺杂Ag和Li的热电性能对比. 物理学报, 2019, 68(11): 117201. doi: 10.7498/aps.68.20190247
    [16] 王先杰, 隋 郁, 千正男, 刘志国, 苗继鹏, 黄喜强, 吕 喆, 朱瑞滨, 程金光, 苏文辉. Fe位Al掺杂对Sr2FeMoO6磁结构和磁输运性质的影响. 物理学报, 2006, 55(2): 849-853. doi: 10.7498/aps.55.849
    [17] 王庆宝, 张仲, 徐锡金, 吕英波, 张芹. N, Fe, La三掺杂锐钛矿型TiO2能带调节的理论与实验研究. 物理学报, 2015, 64(1): 017101. doi: 10.7498/aps.64.017101
    [18] 陶鹏程, 黄燕, 周孝好, 陈效双, 陆卫. 掺杂对金属-MoS2界面性质调制的第一性原理研究. 物理学报, 2017, 66(11): 118201. doi: 10.7498/aps.66.118201
    [19] 朱玥, 李永成, 王福合. Li掺杂对MgH2(001)表面H2分子扩散释放影响的第一性原理研究. 物理学报, 2016, 65(5): 056801. doi: 10.7498/aps.65.056801
    [20] 王冠仕, 林彦明, 赵亚丽, 姜振益, 张晓东. (Cu,N)共掺杂TiO2/MoS2异质结的电子和光学性能:杂化泛函HSE06. 物理学报, 2018, 67(23): 233101. doi: 10.7498/aps.67.20181520
  • 引用本文:
    Citation:
计量
  • 文章访问数:  436
  • PDF下载量:  11
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-12-24
  • 修回日期:  2019-04-10
  • 上网日期:  2019-05-01
  • 刊出日期:  2019-05-20

Be, Mg, Mn掺杂CuInO2形成能的第一性原理研究

  • 1. 广西科技大学理学院, 柳州 545006
  • 2. 广西科技大学材料科学与工程研究中心, 柳州 545006
  • 通信作者: 方志杰, nnfang@semi.ac.cn
    基金项目: 国家自然科学基金(批准号: 11464003, 11864005)、广西自然科学基金(批准号: 2017GXNSFAA198315)、柳州市科技计划项目(批准号: 2016B040202)和广西高校中青年教师基础能力提升项目(批准号: 2018KY0324)资助的课题.

摘要: 研制开发新型的光电材料对促进社会经济发展具有重要的科学意义和实用价值. 利用宽禁带CuInO2铟基材料实现全透明光电材料是目前深入研究的热点. 通过基于密度泛函的第一性原理计算方法, 本文计算出掺杂元素Mg, Be, Mn在CuInO2的形成能. 计算结果表明, 施主类缺陷(如掺杂元素替代Cu原子或进入间隙位置)由于较高的形成能和较深的跃迁能级, 很难在CuInO2材料中出现N型导电; 而受主缺陷中, 在氧原子化学势极大的情况下, Mg原子替代In能成为CuInO2理想的受主缺陷. 计算结果可为制备性能优异的CuInO2材料提供指导.

English Abstract

    • 透明导电氧化物TCO (transparent conductive oxide)不仅具有良好的电学性能, 也具有良好的光学透射性能. 人们将透明导电氧化物广泛应用在太阳能电池、平面显示、特殊功能窗口涂层等光电器件领域[1]. 常见的透明导电氧化物TCO薄膜如ZnO, ln2O3, SnO2[2-4]作为新一类性能优越、用途广泛的TCO薄膜材料体系, 已经工业化并得到广泛使用. 然而, 目前广泛研究和应用的TCO材料大多数是N型导电材料. 为了能够制备出透明二极管以及更多透明导电器件, 需要研制性能优越、工艺稳定以及能够实用化的P型导电TCO材料. 如果将P型TCO材料与N型TCO材料进行复合得到透明导电器件, 将使得透明导电氧化物在应用方面得到很大的发展. 然而与数量众多的N型TCO材料相比, P型TCO材料非常少, 而且P型TCO的导电性能与N型TCO的导电性能相差很大, 通常情况下P型TCO的电导率要比N型TCO的电导率低3—4个数量级. 由于缺乏性能良好的P型TCO材料, 目前仍无法实现具有良好性能的全透明的P-N结, 因此, 研制性能优越、工艺稳定的P型TCO材料成为拓展透明导电氧化物应用所必须面对的课题.

      1997年, 日本的科学家首次设计出具有较高载流子浓度和高电导率的P型TCO材料—CuAlO2[5], 这一研究成果发表在《Nature》上. CuAlO2具有铜铁矿结构, 其室温电导率为0.95 S·cm–1, Seebeck系数为+138 (正值表明材料具有P型导电特性), 光学禁带宽度为3.5 eV, 在500 nm厚度下, CuAlO2的可见光透射率达到80%, 表明该体系可以满足透明电子学的要求. 虽然CuAlO2 的电导率与ZnO等N型TCO材料的电导率相比仍比较低, 但是这一具有里程碑意义的研究发现还是极大地吸引了人们对P型TCO材料的兴趣[6]. 在发现CuAlO2具有P型导电性能之后, 一系列P型铜铁矿类TCO材料也陆续被发现具有P型导电特性[7-14].

      在目前P型铜铁矿类TCO材料中, CuInO2作为具有极大带隙值的P型TCO材料, 而且具有双性掺杂的性能, 即可进行N型和P型掺杂. 因此, CuInO2受到许多人的重视. 但是众多实验结果[15-17]表明CuInO2材料高的电阻率和低的透光率仍旧是制约其发展的两大难题, 如何进一步提高CuInO2材料的导电性能并同时保证其较高的光学性能是以后努力研究的方向之一, 对CuInO2掺杂将成为解决这一问题的有效途径. 而选用何种类型的掺杂元素、采用何种掺杂技术则是制备高质量和高发光性能铟基光电材料的关键所在. 目前对CuInO2掺杂改性的研究已有报道[18-20], 其中吴平等[21]用Sn, Ca元素对CuInO2材料里的In元素进行替位掺杂研究, 发现Sn, Ca掺杂之后的杂质能级都是属于深能级, 因此对提高CuInO2材料的导电性帮助不大. 到目前为止, 对CuInO2材料的掺杂元素选择主要为Sn, Ca元素, 而对具有其他的掺杂金属则还没有涉及. 为了更好地了解对CuInO2材料的掺杂性能, 本文选择掺杂元素Be, Mg, Mn对CuInO2进行掺杂改性研究, 分别对替代Cu、替代In、进入间隙位三种掺杂模式进行基于密度泛函的第一性原理计算, 研究掺杂元素在CuInO2材料中的微观作用机理, 研究成果可为设计和开发高效率的新型光电材料提供崭新的研究思路.

    • CuInO2的晶体结构属于R3m空间群, 晶胞中具有3个特征结构单元: 1)平行C轴分布的O-Cu-O哑铃结构; 2)垂直C轴的六角Cu层; 3) InO6 八面体结构, 其中In位于八面体内(见图1). 总能计算中原子赝势则采用 投影缀加波(PAW)赝势[22,23], 交换关联势采用局域密度近似(LDA)[24,25], 采用软件包VASP进行计算[26,27]. 能带结构计算采用包含四个原子的原胞, 在计算的超晶胞中则包括108个原子, 计算过程中所有的原子都被弛豫到H-F受力最小. 截断能设置为400 eV. 布里渊区的K点设置为1 × 1 × 1 Monkhorst–Pack网格[28,29]. 对于带电荷的缺陷, 计算中引入一个均匀电荷背景以保持整个周期性超晶胞的电中性, 而且所有超晶胞计算都采用同样的K点设置. 单质的化学势也是在LDA和PAW势的总能计算中获得. CuInO2材料中所考虑的价电子为Cu的3d和4p电子, O的2s和2p电子, In的5d和6s电子; 对于掺杂元素Be, Mg, Mn, 价电子包括Be的2s电子, Mg的3s和3p电子, Mn的3d和4s电子.

      图  1  CuInO2晶体结构图, 图中红色原子为O原子, 灰色原子为In原子, 棕色原子为Cu原子 (a)黄色掺杂原子替代Cu原子的情况; (b)绿色掺杂原子替代In原子的情况

      Figure 1.  The crystal structure of the CuInO2, the red atoms are O atoms, the brown atoms are Cu atoms, the purple atoms are In atoms: (a) Substituting yellow dopant atom for Cu atom; (b) substituting green dopant atom for In atom.

      为了提高和调整CuInO2材料P型电导率, 我们计算了掺杂形成能和缺陷跃迁能级. 根据形成能的定义[30,31], 形成能由单质原子的化学势和电子的费米能决定:

      $\begin{split} \Delta {H_{\rm{f}}}(\alpha ,q) =\; & \Delta E(\alpha ,q) + {n_{{\rm{Cu}}}}{\mu _{{\rm{Cu}}}} + {n_{{\rm{In}}}}{\mu _{{\rm{In}}}}\\ & + {n_{\rm{O}}}{\mu _{\rm{O}}} + {n_A}{\mu _A} + q{E_{\rm{F}}}, \end{split} $

      $\begin{split} \Delta E(\alpha ,q) =\;&E(\alpha ,q) - {E_0}({\rm{CuIn}}{{\rm{O}}_{\rm{2}}}) + {n_{{\rm{Cu}}}}{\mu ^0_{{\rm{Cu}}}}\\ &+ {n_{{\rm{In}}}}{\mu ^0_{{\rm{In}}}} + {n_{\rm{O}}}{\mu ^0_{\rm{O}}} + {n_A}{\mu ^0_A} + q{E_{{\rm{VBM}}}}, \end{split} $

      此处E (α, q)为有缺陷α并且带电荷q的超晶胞总能; E0 (CuInO2)为没有任何缺陷的超晶胞总能; nCu, nO, nIn, nA分别表示在形成缺陷时组分Cu, O, In和掺杂元素A从超原胞中转移到原子池中的原子数; q为所带有的电荷数; $\mu _{{\rm{Cu}}}^{\rm{0}},\mu _{{\rm{In}}}^0,\mu _{\rm{O}}^0$为标准化学势, 可分别视为Cu (面心体)、In (面心体)、O (氧气)和掺杂元素的体系总能; EF为体系费米能, 从价带顶VBM开始, 取值范围为价带顶到导带底之间. 在化学势的计算中, 因为单质化学势是与实验中材料生长过程有关的, 因此具有相应的范围.

      对于CuInO2材料, 组分中的化学势范围必须满足以下条件:

      1) 每种组分的化学势不能大于该元素单质构成的体材料化学势, 即

      $ {\mu _{{\rm{Cu}}}} \leqslant 0,{\mu _{{\rm{In}}}} \leqslant 0,{\mu _{\rm{O}}} \leqslant 0; $

      2) 各组分的化学势相加为CuInO2,

      $ {\mu _{{\rm{Cu}}}} + {\mu _{{\rm{In}}}} + 2{\mu _{\rm{O}}} = \Delta {H_{\rm{f}}}{\rm{(CuIn}}{{\rm{O}}_{\rm{2}}}{\rm{),}} $

      CuInO2总能E0(CuInO2) = –651.22 eV, 标准化学势μCu = –4.71 eV, μIn = –3.23 eV, μO = 5.24 eV, 计算所得的CuInO2材料的形成能ΔHf (CuInO2) = –5.00 eV;

      3) 不能够产生组分之间的第二相,

      $ \begin{array}{l} 2{\mu _{{\rm{In}}}} + 3{\mu _{\rm{O}}} \leqslant \Delta {H_{\rm{f}}}{\rm{(I}}{{\rm{n}}_{\rm{2}}}{{\rm{O}}_{\rm{3}}}{\rm{),}}\\ 2{\mu _{{\rm{Cu}}}} + {\mu _{\rm{O}}} \leqslant \Delta {H_{\rm{f}}}{\rm{(C}}{{\rm{u}}_{\rm{2}}}{\rm{O);}} \end{array} $

      4)掺杂的元素不能与原材料的元素组成第二相, 即

      $ {\mu _A} + {\mu _{\rm{O}}} \leqslant \Delta {H_{\rm{f}}}{\rm{(}}A{\rm{O)}}{\rm{.}} $

      计算结果In2O3形成能ΔHf (In2O3) = –9.98 eV, Cu2O形成能ΔHf(Cu2O) = –1.67 eV. 因此, CuInO2材料中的三种成分元素的化学势取值范围分别为: –0. 84 ≤ μCu ≤ 0, –4.99 ≤ µIn ≤ –2.49, –1. 67 ≤ μO ≤ 0.

    • 计算中, 我们对CuInO2的晶格参数(a, c)和原子内部参数u进行了优化计算, CuInO2的晶格参数ac是根据Murnaghan状态方程来优化拟合计算总能量, 而原子内部参数u是根据原子的受力状况来进行优化, CuInO2晶格常数理论计算值与实验值[12]表1.

      CuInO2 晶格常数理论计算值 晶格常数实验值
      a 3.297 Å a 3.292 Å
      c 17.192 Å c 17.388 Å
      u 0.1047 Å u 0.1061 Å
      c/a 5.214 c/a 5.282

      表 1  CuInO2晶格常数理论计算值与实验值

      Table 1.  Theoretical values and experimental values of lattice constants in CuInO2.

      掺杂元素的进入使得CuInO2晶体内部平衡状态受到了破坏, 在掺杂元素周围的原子将发生不同程度的弛豫. 表2列出了CuInO2中掺杂元素周围的弛豫变化情况, 表中MgCu表示Mg原子以替位掺杂形式进入Cu原子位置, MgIn表示Mg原子以替位掺杂形式进入In原子位置, Mgi则表示Mg原子进入CuInO2间隙位置. 键距后括号为掺杂原子的近邻原子个数和种类. 从表2可以看出, 与CuInO2的CuO键长1.85 Å(O, 2), InO键2.18 Å(O, 6)相比, Mg元素的替换并没有引起周围氧原子大的弛豫变化, 而掺杂元素Be, Mn元素则引起了邻域原子较大的弛豫变化.

      掺杂种类 最近邻键距 掺杂种类 最近邻键距 掺杂种类 最近邻键距
      BeCu 1.49 Å(O, 2) MgCu 1.85 Å(O, 2) MnCu 1.77Å(O, 2)
      BeIn 1.96 Å(O, 6) MgIn 2.12 Å(O, 6) MnIn 1.93Å(O, 6)
      Bei 1.62 Å(O, 3) Mgi 2.04 Å(O, 3) Mni 1.99Å(O, 3)

      表 2  掺杂元素周围的弛豫变化情况, 键距后括号为掺杂原子的近邻原子种类和个数

      Table 2.  The surrounding atoms of dopant, kinds and numbers of dopantxs nearest neighbor atoms in parentheses.

      表3是掺杂下CuInO2在电中性状态下当q = 0和${\mu _i} = 0$时的外在掺杂形成能$\Delta E(\alpha ,q)$, 由此得出在不同的化学势条件下缺陷形成能可以通过(1)式推导. 由于掺杂形成能与化学势、费米能都有关系, 可分为几种极值情况进行分析. 在CuInO2材料中, 化学势条件分两种情况, 一种为有阴离子化学势极大(阳离子化学势极小)的情况, 另外一种则为阴离子化学势极小(而阳离子化学势极大)的情况.

      掺杂类型 形成能 掺杂类型 总能/eV 形成能 掺杂类型 总能/eV 形成能
      BeCu –2.20 MgCu –650.87 –2.58 MnCu –654.63 1.76
      BeIn 0.79 MgIn –651.36 –1.59 MnIn –656.54 1.33
      Bei –1.05 Mgi –654.27 –1.27 Mni –657.39 3.71

      表 3  CuInO2的掺杂元素形成能

      Table 3.  The calculated formation energies of dopants in CuInO2.

      图2图3分别为氧化学势极大和极小情况下, Mg, Mn, Be掺杂元素形成能随费米能的变化图, 纵坐标表示形成能(形成能在零以下表示易于形成), 实线表示掺杂的形成能随着费米能变化的过程, 数字表示带电情况, 小实心点表示跃迁能级, VBM为价带顶, CBM为导带底, 费米能取值在VBM与CBM之间. 当费米能EF移动接近导带底时, 根据(1)式, 受主类(Mg, Mn, Be替换In原子)的掺杂形成能就会降低, 而施主类(Mg, Mn, Be替换Cu原子或Mg, Mn, Be进入间隙位)的掺杂形成能将会增高. 由图2 (氧化学势极大)可以看出, 氧化学势极大情况下, 从掺杂元素的形成能分析, Mg, Mn原子都容易替换CuInO2里面的In原子, 随着费米能级向导带底移动, Be原子也有替换In原子的可能; 对Mg, Mn, Be替换Cu原子的情况, 计算结果显示Mg, Mn, Be替换Cu原子的行为具有较高的形成能, 表明这些掺杂元素很难替换Cu原子; 同时也因为极高的形成能, Mg, Mn, Be掺杂元素进入CuInO2间隙位置的情况不会发生. 而从图3可以看出, 在氧化学势极小的情况, 虽然Be, Mg都易于进入CuInO2, 但是会同时形成受主类和施主类缺陷, 这极大地影响了材料的导电性能和效率; 随着费米能EF向导带底移动, MnIn这类形式的掺杂也能发生, 而MnCu, Mni的形成能较大, 难以形成在CuInO2材料中.

      图  2  掺杂形成能在cation-poor, anion-rich的变化图

      Figure 2.  The change of doping formation energies under cation-poor, anion-rich condition.

      图  3  掺杂形成能在cation-rich, anion-poor的变化图

      Figure 3.  The change of doping formation energies under cation-rich, anion-poor condition.

      图4为各种类型掺杂元素进入CuInO2后所形成的缺陷能级, 在确定掺杂元素能级位置时, 我们采用芯能级对齐方式[32]进行能级位置确定. 从图4可看出, Ⅱ族元素Be, Mg替代In元素之后所形成的缺陷能级都很接近价带顶, MgIn和BeIn的位置分别为价带顶上0.05 eV和0.06 eV, 均属于缺陷浅能级, 而MnIn的位置位于价带顶上0.4 eV, 属于深能级. 这表明MgIn和BeIn很适合作为受主在CuInO2中产生空穴载流子, 再结合图2图3的形成能计算结果, Mg原子替换In原子在氧化学势极大的情况下能大量生成且不存在自补偿效应(MgCu和Mgi在该情况不会形成), 因此在氧化学势极大的情况下, MgIn是CuInO2材料中最为理想的受主, 并能较好地提高材料的P型导电性能. 为了增加CuInO2中的P型导电性, 通过原子层外延生长的控制, 在MBE或MOCVD晶体生长过程中实现减少各类金属的偏压以及增加氧偏压的方法以增加Mg原子替代In原子的概率, 这种选择性位置的掺杂方法能够提高材料的P型导电性能. 而CuInO2的外在掺杂N型缺陷都具有较高的形成能或较高的跃迁能级, 因此, 对于Mg, Mn, Be元素, 无论何种形式的掺杂均很难在CuInO2中实现N型导电性能.

      图  4  掺杂元素在CuInO2的缺陷跃迁能级

      Figure 4.  The calculated transition energy levels of the extrinsic defects in CuInO2.

    • 利用第一性原理计算方法, 研究掺杂元素Mg, Be, Mn对CuInO2材料的作用机理, 分别对替代Cu、替代In、进入间隙位三种掺杂模式进行计算. 计算结果发现, Mg原子的替位掺杂不会引起CuInO2体材料内部原子弛豫, 在氧化学势极大的情况下, MgIn能大量生成且缺陷位置在价带顶上方0.05 eV, 同时也避免了Mgi, MgCu的形成, MgIn因此能成为CuInO2的理想受主缺陷, 能较好地提高材料的P型导电性能. 另外, 因为较高的形成能和深的跃迁能级, Mg, Be, Mn三种元素替代Cu和进入间隙位的掺杂模式都不会给材料带来N型导电性能.

参考文献 (32)

目录

    /

    返回文章
    返回