搜索

x

留言板

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

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

应变对(Ga, Mo)Sb磁学和光学性质影响的理论研究

潘凤春 林雪玲 王旭明

引用本文:
Citation:

应变对(Ga, Mo)Sb磁学和光学性质影响的理论研究

潘凤春, 林雪玲, 王旭明

First-principles study of strain effect on magnetic and optical properties in (Ga, Mo)Sb

Pan Feng-Chun, Lin Xue-Ling, Wang Xu-Ming
PDF
HTML
导出引用
  • 近年来, 作为一种自旋电子学领域的关键材料, 具有高温本征铁磁性的稀磁半导体受到了广泛的关注. 为探索能够提高本征铁磁性居里温度(Curie temperature, TC)的方法, 本文运用第一性原理LDA+U方法研究了应变对Mo掺杂GaSb的电子结构、磁学及光学性质的影响. 研究结果表明: –6%—2.5%应变范围下GaSb半导体材料具有稳定的力学性能, 压应变下GaSb材料的可塑性、韧性增强, 有利于GaSb半导体材料力学性能的提升; 应变对Mo替代Ga缺陷(MoGa)的电子结构有重要的影响, –3%至–1.2%应变范围下MoGa处于低自旋态(low spin state, LSS), 具有1${\mu _{{\rm{B}}}}$的局域磁矩, –1.1%—2%应变范围下MoGa处于高自旋态(high spin state, HSS), 具有3${\mu _{{\rm{B}}}}$的磁矩; 不管是LSS还是HSS, MoGa产生局域磁矩之间的耦合都是铁磁耦合, 但铁磁耦合的强度和物理机制不同, 适当的压应变可有效提高铁磁耦合强度, 这有利于实现高TC的GaSb基磁性半导体; Mo可极大提高GaSb半导体材料的电极化能力, 这有利于光生电子-空穴对的形成和分离, 提高掺杂体系对长波光子的光电转化效率; Mo引入的杂质能级使电子的带间跃迁对所需要吸收光子的能量变小, 掺杂体系光学吸收谱的吸收边发生了红移, 拉应变可进一步提升(Ga,Mo)Sb体系在红外光区的光学性能.
    In recent decades, as a kind of critical material in spintronics field, the diluted magnetic semiconductor with high temperature and intrinsic ferromagnetism has attracted extensive attention. In order to explore the approach to enhancing Curie temperature (TC), the LDA+U method of the first-principles calculation is adopted to study the effect of strain on electronic structure, magnetic and optical properties in Mo doped GaSb system. The results indicate that the structure of GaSb is stable with strain in a range of –6%—2.5%. Plasticity and toughness of GaSb increase under compressive strains, which is conducive to the improvement of the mechanical properties. The strain affects the electronic structure of MoGa greatly. In a range from –3% to –1.2%, MoGa is in the low spin state (LSS) with a 1${\mu _{{\rm{B}}}}$local magnetic moment, while in a range of –1.1%—2%, MoGa is in high spin state (HSS) with a 3${\mu _{{\rm{B}}}}$moment. The magnetic interactions between MoGa and MoGa are all ferromagnetic for LSS and so is the case for HSS, although they are different in coupling intensity and mechanism. In particular, appropriate compressive strains can improve the strength of ferromagnetic coupling effectively and are favorable for the preparation of the GaSb-based diluted magnetic semiconductors with high Curie temperatures and inherent ferromagnetism. Moreover, we find that Mo can greatly improve the polarization capability of GaSb and play a vital role in forming and separating the electron-hole pairs, and thus further improving the photoelectric conversion capability for long wave photons. The energy required to absorb photons for inter-band transition of electrons decreases because of the impurity levels induced by Mo, which leads the absorption edge to be red-shifted. The optical properties of (Ga,Mo)Sb in infrared region are further enhanced by the tensile strain.
      通信作者: 王旭明, wang_xm@126.com
    • 基金项目: 国家自然科学基金(批准号: 11764032, 11665018)资助的课题.
      Corresponding author: Wang Xu-Ming, wang_xm@126.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 11764032, 11665018).
    [1]

    Ohno H 1998 Science 281 951Google Scholar

    [2]

    Dietl T, Ohno H, Matsukura F, Cibert J, Ferrand D 2000 Science 287 1019Google Scholar

    [3]

    Dietl T 2010 Nat. Mater. 9 965Google Scholar

    [4]

    Sato K, Bergqvist L, Kudrnovsky J, Dederichs P H, Eriksson O, Turek I, Sanyal B, Bouzerar G, Katayama-Yoshida H, Dinh V A, Fukushima T, Kizaki H, Zeller R 2010 Rev. Modern Phys. 82 1633Google Scholar

    [5]

    Dietl T, Ohno H 2014 Rev. Mod. Phys. 86 187Google Scholar

    [6]

    王少霞, 赵旭才, 潘多桥, 庞国旺, 刘晨曦, 史蕾倩, 刘桂安, 雷博程, 黄以能, 张丽丽 2020 物理学报 69 197101Google Scholar

    Wang S X, Zhao X C, Pan D Q, Pang G W, Liu C X, Shi L Q, Liu G A, Lei B C, Huang Y N, Zhang L L 2020 Acta Phys. Sin. 69 197101Google Scholar

    [7]

    姚仲瑜, 孙丽, 潘孟美, 孙书娟, 刘汉军 2018 物理学报 67 217501Google Scholar

    Yao Z Y, Sun L, Pan M M, Sun S J, Liu H J 2018 Acta Phys. Sin. 67 217501Google Scholar

    [8]

    You J Y, Gu B, Maekawa S, Su Gang 2020 Phys. Rev. B 102 094432Google Scholar

    [9]

    Takeda T, Suzuki M, Anh L D, Tu N, Schmitt T, Yoshida S, Sakano M, Ishizaka K, Takeda Y, Fijimori S 2020 Phys. Rev. B 101 155142Google Scholar

    [10]

    Takeda T, Sakamoto S, Araki K, Fujisawa Y, Anh L, Tu N, Takeda Y, Fujimori S, Fujimori A, Tanaka M 2020 Phys. Rev. B 102 245203Google Scholar

    [11]

    Sena N, Dussan A, Mesa F, Castao E, Gonzalez-Hernandez R 2016 J. Appl. Phys. 120 051704Google Scholar

    [12]

    Kudrnovsky J, Turek I, Drchal V, Maca F, Weinberger P, Bruno P 2004 Phys. Rev. B 69 115208Google Scholar

    [13]

    Da Pieve F, Di Matteo S, Rangel T, Giantomassi M, Lamoen D, Rignanese G M, Gonze X 2013 Phys. Rev. Lett. 110 136402Google Scholar

    [14]

    Chou H, Lin C P, Huang J C A, Hsu H S 2008 Phys. Rev. B 77 245210Google Scholar

    [15]

    Sato K, Dederichs P H, Katayama-Yoshida H, Kudrnovsky J 2004 J. Phys.: Condens. Matter 16 S5491Google Scholar

    [16]

    Tu N T, Hai P N, Anh L D, Tanaka M 2015 Phys. Rev. B 92 144403Google Scholar

    [17]

    Tu N T, Hai P N, Le D A, Tanaka M 2014 Appl. Phys. Lett. 105 132402Google Scholar

    [18]

    Tu N T, Hai P N, Le D A, Tanaka M 2016 Appl. Phys. Lett. 108 192401Google Scholar

    [19]

    Kondrin M V, Gizatullin V R, Popova S V, Lyapin A G, Brazhkin V V, Ivanov V Y, Pronin A A, Lebed Y B, Sadykov R A 2011 J. Phys. Condens. Matter 23 446001Google Scholar

    [20]

    Ganesan K, Pendyala N B, Koteswara Rao K S R, Venkataraman V, Bhat H L 2010 Semicond. Sci. Technol. 25 105003Google Scholar

    [21]

    Boishin G I, Sullivan J M, Whitman L J 2005 Phys. Rev. B 71 193307Google Scholar

    [22]

    Abe E, Matsukura F, Yasuda H, Ohno Y, Ohno H 2000 Physica E 7 981Google Scholar

    [23]

    Medvedeva J E 2006 Phys. Rev. Lett. 97 086401Google Scholar

    [24]

    Park C Y, Yoon S G, Jo Y H, Shin S C 2009 Appl. Phys. Lett. 95 122502Google Scholar

    [25]

    Park C Y, You C Y, Jeon K R, Shin S C 2012 Appl. Phys. Lett. 100 222409Google Scholar

    [26]

    Egbo K O, Adesina A E, Ezeh C V, Liu C P, Yu K M 2021 Phys. Rev. Mater. 5 094603Google Scholar

    [27]

    Lu W J, Sun Y P, Zhao B C, Zhu X B, Song W H 2006 Phys. Rev. B 73 174425Google Scholar

    [28]

    Dwivedi G D, Tseng K F, Chan C L, Shahi P, Louremban J, Chatterjee B, Ghosh A K, Yang H D, Chatterjee S 2010 Phys. Rev. B 82 134428Google Scholar

    [29]

    Linpeng X, Karin T, Durnev M V, Glazov M M, Schott R, Wieck A D, Ludwig A, Fu K M C 2021 Phys. Rev. B 103 115412Google Scholar

    [30]

    Patel K, Prosandeev S, Xu B, Xu C S, Bellaiche L 2021 Phys. Rev. B 103 094103Google Scholar

    [31]

    Zhong H X, Xiong W Q, Lv P F, Yu J, Yuan S J 2021 Phys. Rev. B 103 085124Google Scholar

    [32]

    Goel S, Anh L D, Ohya S, Tanaka M 2019 Phys. Rev. B 99 014431Google Scholar

    [33]

    Rawat K, Fong D D, Aidhy D S 2021 J. Appl. Phys. 129 095301Google Scholar

    [34]

    Breev I D, Poshakinskiy A V, Yakovleva V V, et al. 2021 Appl. Phys. Lett. 118 084003Google Scholar

    [35]

    Hohenberg P, Kohn W 1964 Phys. Rev. 136 B864Google Scholar

    [36]

    Perdew J P, Wang Y 1992 Phys. Rev. B 45 13244Google Scholar

    [37]

    Segall M D, Lindan P J D, Probert M J, Pickard C J, Hasnip P J, Clark S J, Payne M C 2002 J. Phys.: Condens. Matter 14 2717Google Scholar

    [38]

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

    [39]

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

    [40]

    Straumanis M E, Kim C D 1965 J. Electrochem. Soc. 112 112Google Scholar

    [41]

    Straumanis M E, Kim C D 1965 J. Appl. Phys. 36 3822Google Scholar

    [42]

    Jakowetz W, Ruhle W, Breuninger, Pilkuhn M 1972 Phys. Status Solidi (A) 12 169Google Scholar

    [43]

    Van de Walle C G, Neugebauer J 2004 J. Appl. Phys. 95 3851Google Scholar

    [44]

    Zota C B, Kim S H, Yokoyama M, Takenaka M, Takagi S 2012 Appl. Phys. Express 5 071201Google Scholar

    [45]

    Yokoyama M, Nishi K, Kim S, Yokoyama H, Takenaka M, Takagi S 2014 Appl. Phys. Lett. 104 093509Google Scholar

    [46]

    Watt J P 1980 J. Appl. Phys. 51 1520Google Scholar

    [47]

    Reuss A 1929 Z. Ang. Math. Mech. 9 49Google Scholar

    [48]

    Wu Z J, Zhao E J, Xiang H P, Hao X F, Liu X J, Meng J 2007 Phys. Rev. B 76 054115Google Scholar

    [49]

    Hill B R 1952 Proc. Phys. Soc. A65 349

    [50]

    Li K Y, Ding Z S, Xue D F 2011 Phys. Status Solidi (B) 248 1227Google Scholar

    [51]

    Hao Y L, Li S J, Sun S Y, Zheng C Y, Yang R 2007 Acta Biomater. 3 277Google Scholar

    [52]

    Pugh S F 2009 Philos. Mag. 45 823Google Scholar

    [53]

    Ling T, Yan D Y, Wang H, Jiao Y, Hu Z P, Zheng Y, Zheng L R, Mao J, Liu H, Du X W, Jaroniec M, Qiao S Z 2017 Nat. Commun. 8 1509Google Scholar

    [54]

    Petrie J R, Jeen H, Barron S C, Meyer T L, Lee H N 2016 J. Am. Chem. Soc. 138 7252Google Scholar

    [55]

    Wolfschmidt H, Bussar R, Stimming U 2008 J. Phys. Condens. Matter 20 374127Google Scholar

    [56]

    Tan Y, Liu P, Chen L, Cong W, Ito Y, Han J, Guo X, Tang Z, Fujita T, Hirata A, Chen M 2014 Adv. Mater. 26 8023Google Scholar

    [57]

    王志伟, 郝永芹, 李洋, 谢检来, 王霞, 晏长岭, 魏志鹏, 马晓辉 2018 发光学报 39 935Google Scholar

    Wang Z W, Hao Y Q, Li Y, Xie J L, Wang X, Yan C L, Wei Z P, Ma X H 2018 Chin. J. Lumi. 39 935Google Scholar

    [58]

    Zhang P, Kim Y H, Wei S H 2019 Phys. Rev. Appl. 11 054058Google Scholar

    [59]

    Pan F C, Lin X L, Wang X M 2021 Chin. Phys. B 30 096105Google Scholar

    [60]

    Dietl T, Ohno H 2001 Physica E Low Dimens. Syst. Nanostruct. 9 185Google Scholar

    [61]

    Matsukura F, Ohno H, Dietl T 2002 Handbook of Magnetic Materials 14 1Google Scholar

    [62]

    Yakunin A M, Silov A Y, Koenraad P M, Wolter J H, Roy W V, Boeck J D 2003 Superlattices Microst. 34 539Google Scholar

    [63]

    Khvostikov V P, Khvostikova O A, Gazaryan P Y, et al. 2007 J. Sol. Energy Engineer. 129 291Google Scholar

    [64]

    Kim J M, Dutta P S, Brown E, Borrego J M, Greiff P 2013 Semicond. Sci. Technol. 28 065002Google Scholar

    [65]

    Liu Z, Qiu K 2017 Energy 141 892Google Scholar

    [66]

    Lou Y Y, Zhang X L, Huang A B, Wang Y 2017 Sol. Energy Mater. Sol. Cells 172 124Google Scholar

    [67]

    Li M Z, Chen X L, Li H L, et al. 2018 Chin. Phys. B 27 078101Google Scholar

    [68]

    黄昆 1988 固体物理学 (北京: 高等教育出版社) 第437−452页

    Huang K 1988 Solid State Physics (Beijing: Higher Education Press) pp437−452 (in Chinese)

    [69]

    潘凤春, 林雪玲, 曹志杰, 李小伏 2019 物理学报 68 184202Google Scholar

    Pan F C, Lin X L, Cao Z J, Li X F 2019 Acta Phys. Sin. 68 184202Google Scholar

    [70]

    Uspenskii Y A, Kulatov E T 2009 J. Magn. Magn. Mater. 321 931Google Scholar

    [71]

    杨卫霞, 张贺翔, 潘凤春, 林雪玲 2021 人工晶体学报 50 2019Google Scholar

    Yang W X, Zhang H X, Pan F C, Lin X L 2021 J. Synth. Cryst. 50 2019Google Scholar

  • 图 1  GaSb超晶胞结构, 其中大的绿色球代表Ga原子, 小的紫色球代表Sb原子

    Fig. 1.  Structure of GaSb supercell, where the big green balls and small purple balls denote Ga and Sb atoms, respectively.

    图 2  力学稳定性判据

    Fig. 2.  Criteria of mechanical stability.

    图 3  (Ga, Mo)Sb体系在HSS和LSS下的态密度图 (a) 零应变下的TDOS; (b)零应变下Mo-4d电子的PDOS; (c) –1.2%应变下的TDOS; (d) –1.2%应变下Mo-4d电子的PDOS

    Fig. 3.  Density of states (DOS) of (Ga, Mo)Sb in HSS and LSS: (a) TDOS and (b) PDOS of Mo-4d without strain; (c) the TDOS and (d) the PDOS of Mo-4d under –1.2% strain.

    图 4  (a) 应变下, (Ga, Mo)Sb体系中掺杂Mo原子和近邻Sb原子的磁矩; (b) 应变下, S(12)和S(23)结构的ΔE

    Fig. 4.  (a) Magnetic moments contributed by Mo and neighbor Sb (N Sb) under strains and (b) ΔE of S(12) and S(23) under strains, respectively.

    图 5  等值面为0.01 e/Å3的空间自旋密度图 (a) HSS下MoGa; (b) LSS下MoGa; (c) HSS下S(12)结构FM耦合; (d) LSS下S(23)结构FM耦合

    Fig. 5.  The spatial spin density for (a) MoGa under HSS, (b) MoGa under LSS, (c) FM coupling of S(12) under HSS, (d) FM coupling of S(23) under LSS, respectively. The isovalue is set to 0.01 e/Å3.

    图 6  d电子轨道四面体晶场分裂示意图

    Fig. 6.  Diagram of Mo-4d orbital splitting under tetrahedral crystal field.

    图 7  (Ga, Mo)Sb体系的能带图 (a)零应变下HSS; (b) –1.2%应变下LSS. 黑色实线表示自旋向上的能带, 红色虚线表示自旋向下的能带. 在能量为0处的水平蓝色虚线表示费米能级

    Fig. 7.  Band structures of (Ga, Mo)Sb in (a) HSS without strain and (b) LSS under –1.2% strain. Black solid lines and red dotted lines denote spin-up and spin-down bands, respectively. The horizontal blue dotted lines located at 0 energy are Fermi levels.

    图 8  (a) 应变下, GaSb和(Ga, Mo)Sb的光学吸收谱; (b) 应变下, GaSb和(Ga, Mo)Sb的复介电函数, 实线代表实部, 虚线代表虚部

    Fig. 8.  (a) Absorption spectra and (b) complex dielectric function of GaSb and (Ga, Mo)Sb under strains. In panel (b) solid and dotted lines denote real and imaginary parts, respectively.

    表 1  应变下GaSb的体积模量B、剪切模量G、泊松比$\gamma $B/G

    Table 1.  Bulk modulus B, shear modulus G, Poisson ratio $\gamma $ and B/G under strains.

    Strain/%B/GPaG/GPa$\gamma $B/G
    –6148.2137.670.383.93
    –5129.6643.080.353.01
    –3.7109.9343.850.322.51
    –290.0343.300.292.08
    068.7747.210.221.46
    0.568.7347.170.221.47
    164.2645.470.211.48
    1.561.3540.380.201.52
    255.4938.020.201.53
    下载: 导出CSV
  • [1]

    Ohno H 1998 Science 281 951Google Scholar

    [2]

    Dietl T, Ohno H, Matsukura F, Cibert J, Ferrand D 2000 Science 287 1019Google Scholar

    [3]

    Dietl T 2010 Nat. Mater. 9 965Google Scholar

    [4]

    Sato K, Bergqvist L, Kudrnovsky J, Dederichs P H, Eriksson O, Turek I, Sanyal B, Bouzerar G, Katayama-Yoshida H, Dinh V A, Fukushima T, Kizaki H, Zeller R 2010 Rev. Modern Phys. 82 1633Google Scholar

    [5]

    Dietl T, Ohno H 2014 Rev. Mod. Phys. 86 187Google Scholar

    [6]

    王少霞, 赵旭才, 潘多桥, 庞国旺, 刘晨曦, 史蕾倩, 刘桂安, 雷博程, 黄以能, 张丽丽 2020 物理学报 69 197101Google Scholar

    Wang S X, Zhao X C, Pan D Q, Pang G W, Liu C X, Shi L Q, Liu G A, Lei B C, Huang Y N, Zhang L L 2020 Acta Phys. Sin. 69 197101Google Scholar

    [7]

    姚仲瑜, 孙丽, 潘孟美, 孙书娟, 刘汉军 2018 物理学报 67 217501Google Scholar

    Yao Z Y, Sun L, Pan M M, Sun S J, Liu H J 2018 Acta Phys. Sin. 67 217501Google Scholar

    [8]

    You J Y, Gu B, Maekawa S, Su Gang 2020 Phys. Rev. B 102 094432Google Scholar

    [9]

    Takeda T, Suzuki M, Anh L D, Tu N, Schmitt T, Yoshida S, Sakano M, Ishizaka K, Takeda Y, Fijimori S 2020 Phys. Rev. B 101 155142Google Scholar

    [10]

    Takeda T, Sakamoto S, Araki K, Fujisawa Y, Anh L, Tu N, Takeda Y, Fujimori S, Fujimori A, Tanaka M 2020 Phys. Rev. B 102 245203Google Scholar

    [11]

    Sena N, Dussan A, Mesa F, Castao E, Gonzalez-Hernandez R 2016 J. Appl. Phys. 120 051704Google Scholar

    [12]

    Kudrnovsky J, Turek I, Drchal V, Maca F, Weinberger P, Bruno P 2004 Phys. Rev. B 69 115208Google Scholar

    [13]

    Da Pieve F, Di Matteo S, Rangel T, Giantomassi M, Lamoen D, Rignanese G M, Gonze X 2013 Phys. Rev. Lett. 110 136402Google Scholar

    [14]

    Chou H, Lin C P, Huang J C A, Hsu H S 2008 Phys. Rev. B 77 245210Google Scholar

    [15]

    Sato K, Dederichs P H, Katayama-Yoshida H, Kudrnovsky J 2004 J. Phys.: Condens. Matter 16 S5491Google Scholar

    [16]

    Tu N T, Hai P N, Anh L D, Tanaka M 2015 Phys. Rev. B 92 144403Google Scholar

    [17]

    Tu N T, Hai P N, Le D A, Tanaka M 2014 Appl. Phys. Lett. 105 132402Google Scholar

    [18]

    Tu N T, Hai P N, Le D A, Tanaka M 2016 Appl. Phys. Lett. 108 192401Google Scholar

    [19]

    Kondrin M V, Gizatullin V R, Popova S V, Lyapin A G, Brazhkin V V, Ivanov V Y, Pronin A A, Lebed Y B, Sadykov R A 2011 J. Phys. Condens. Matter 23 446001Google Scholar

    [20]

    Ganesan K, Pendyala N B, Koteswara Rao K S R, Venkataraman V, Bhat H L 2010 Semicond. Sci. Technol. 25 105003Google Scholar

    [21]

    Boishin G I, Sullivan J M, Whitman L J 2005 Phys. Rev. B 71 193307Google Scholar

    [22]

    Abe E, Matsukura F, Yasuda H, Ohno Y, Ohno H 2000 Physica E 7 981Google Scholar

    [23]

    Medvedeva J E 2006 Phys. Rev. Lett. 97 086401Google Scholar

    [24]

    Park C Y, Yoon S G, Jo Y H, Shin S C 2009 Appl. Phys. Lett. 95 122502Google Scholar

    [25]

    Park C Y, You C Y, Jeon K R, Shin S C 2012 Appl. Phys. Lett. 100 222409Google Scholar

    [26]

    Egbo K O, Adesina A E, Ezeh C V, Liu C P, Yu K M 2021 Phys. Rev. Mater. 5 094603Google Scholar

    [27]

    Lu W J, Sun Y P, Zhao B C, Zhu X B, Song W H 2006 Phys. Rev. B 73 174425Google Scholar

    [28]

    Dwivedi G D, Tseng K F, Chan C L, Shahi P, Louremban J, Chatterjee B, Ghosh A K, Yang H D, Chatterjee S 2010 Phys. Rev. B 82 134428Google Scholar

    [29]

    Linpeng X, Karin T, Durnev M V, Glazov M M, Schott R, Wieck A D, Ludwig A, Fu K M C 2021 Phys. Rev. B 103 115412Google Scholar

    [30]

    Patel K, Prosandeev S, Xu B, Xu C S, Bellaiche L 2021 Phys. Rev. B 103 094103Google Scholar

    [31]

    Zhong H X, Xiong W Q, Lv P F, Yu J, Yuan S J 2021 Phys. Rev. B 103 085124Google Scholar

    [32]

    Goel S, Anh L D, Ohya S, Tanaka M 2019 Phys. Rev. B 99 014431Google Scholar

    [33]

    Rawat K, Fong D D, Aidhy D S 2021 J. Appl. Phys. 129 095301Google Scholar

    [34]

    Breev I D, Poshakinskiy A V, Yakovleva V V, et al. 2021 Appl. Phys. Lett. 118 084003Google Scholar

    [35]

    Hohenberg P, Kohn W 1964 Phys. Rev. 136 B864Google Scholar

    [36]

    Perdew J P, Wang Y 1992 Phys. Rev. B 45 13244Google Scholar

    [37]

    Segall M D, Lindan P J D, Probert M J, Pickard C J, Hasnip P J, Clark S J, Payne M C 2002 J. Phys.: Condens. Matter 14 2717Google Scholar

    [38]

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

    [39]

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

    [40]

    Straumanis M E, Kim C D 1965 J. Electrochem. Soc. 112 112Google Scholar

    [41]

    Straumanis M E, Kim C D 1965 J. Appl. Phys. 36 3822Google Scholar

    [42]

    Jakowetz W, Ruhle W, Breuninger, Pilkuhn M 1972 Phys. Status Solidi (A) 12 169Google Scholar

    [43]

    Van de Walle C G, Neugebauer J 2004 J. Appl. Phys. 95 3851Google Scholar

    [44]

    Zota C B, Kim S H, Yokoyama M, Takenaka M, Takagi S 2012 Appl. Phys. Express 5 071201Google Scholar

    [45]

    Yokoyama M, Nishi K, Kim S, Yokoyama H, Takenaka M, Takagi S 2014 Appl. Phys. Lett. 104 093509Google Scholar

    [46]

    Watt J P 1980 J. Appl. Phys. 51 1520Google Scholar

    [47]

    Reuss A 1929 Z. Ang. Math. Mech. 9 49Google Scholar

    [48]

    Wu Z J, Zhao E J, Xiang H P, Hao X F, Liu X J, Meng J 2007 Phys. Rev. B 76 054115Google Scholar

    [49]

    Hill B R 1952 Proc. Phys. Soc. A65 349

    [50]

    Li K Y, Ding Z S, Xue D F 2011 Phys. Status Solidi (B) 248 1227Google Scholar

    [51]

    Hao Y L, Li S J, Sun S Y, Zheng C Y, Yang R 2007 Acta Biomater. 3 277Google Scholar

    [52]

    Pugh S F 2009 Philos. Mag. 45 823Google Scholar

    [53]

    Ling T, Yan D Y, Wang H, Jiao Y, Hu Z P, Zheng Y, Zheng L R, Mao J, Liu H, Du X W, Jaroniec M, Qiao S Z 2017 Nat. Commun. 8 1509Google Scholar

    [54]

    Petrie J R, Jeen H, Barron S C, Meyer T L, Lee H N 2016 J. Am. Chem. Soc. 138 7252Google Scholar

    [55]

    Wolfschmidt H, Bussar R, Stimming U 2008 J. Phys. Condens. Matter 20 374127Google Scholar

    [56]

    Tan Y, Liu P, Chen L, Cong W, Ito Y, Han J, Guo X, Tang Z, Fujita T, Hirata A, Chen M 2014 Adv. Mater. 26 8023Google Scholar

    [57]

    王志伟, 郝永芹, 李洋, 谢检来, 王霞, 晏长岭, 魏志鹏, 马晓辉 2018 发光学报 39 935Google Scholar

    Wang Z W, Hao Y Q, Li Y, Xie J L, Wang X, Yan C L, Wei Z P, Ma X H 2018 Chin. J. Lumi. 39 935Google Scholar

    [58]

    Zhang P, Kim Y H, Wei S H 2019 Phys. Rev. Appl. 11 054058Google Scholar

    [59]

    Pan F C, Lin X L, Wang X M 2021 Chin. Phys. B 30 096105Google Scholar

    [60]

    Dietl T, Ohno H 2001 Physica E Low Dimens. Syst. Nanostruct. 9 185Google Scholar

    [61]

    Matsukura F, Ohno H, Dietl T 2002 Handbook of Magnetic Materials 14 1Google Scholar

    [62]

    Yakunin A M, Silov A Y, Koenraad P M, Wolter J H, Roy W V, Boeck J D 2003 Superlattices Microst. 34 539Google Scholar

    [63]

    Khvostikov V P, Khvostikova O A, Gazaryan P Y, et al. 2007 J. Sol. Energy Engineer. 129 291Google Scholar

    [64]

    Kim J M, Dutta P S, Brown E, Borrego J M, Greiff P 2013 Semicond. Sci. Technol. 28 065002Google Scholar

    [65]

    Liu Z, Qiu K 2017 Energy 141 892Google Scholar

    [66]

    Lou Y Y, Zhang X L, Huang A B, Wang Y 2017 Sol. Energy Mater. Sol. Cells 172 124Google Scholar

    [67]

    Li M Z, Chen X L, Li H L, et al. 2018 Chin. Phys. B 27 078101Google Scholar

    [68]

    黄昆 1988 固体物理学 (北京: 高等教育出版社) 第437−452页

    Huang K 1988 Solid State Physics (Beijing: Higher Education Press) pp437−452 (in Chinese)

    [69]

    潘凤春, 林雪玲, 曹志杰, 李小伏 2019 物理学报 68 184202Google Scholar

    Pan F C, Lin X L, Cao Z J, Li X F 2019 Acta Phys. Sin. 68 184202Google Scholar

    [70]

    Uspenskii Y A, Kulatov E T 2009 J. Magn. Magn. Mater. 321 931Google Scholar

    [71]

    杨卫霞, 张贺翔, 潘凤春, 林雪玲 2021 人工晶体学报 50 2019Google Scholar

    Yang W X, Zhang H X, Pan F C, Lin X L 2021 J. Synth. Cryst. 50 2019Google Scholar

  • [1] 李发云, 杨志雄, 程雪, 甄丽营, 欧阳方平. 单层缺陷碲烯电子结构与光学性质的第一性原理研究. 物理学报, 2021, 70(16): 166301. doi: 10.7498/aps.70.20210271
    [2] 叶建峰, 秦铭哲, 肖清泉, 王傲霜, 何安娜, 谢泉. Ti, V, Co, Ni掺杂二维CrSi2材料的电学、磁学及光学性质的第一性原理研究. 物理学报, 2021, 70(22): 227301. doi: 10.7498/aps.70.20211023
    [3] 王闯, 赵永红, 刘永. Ga1–xCrxSb (x = 0.25, 0.50, 0.75) 磁学和光学性质的第一性原理研究. 物理学报, 2019, 68(17): 176301. doi: 10.7498/aps.68.20182305
    [4] 潘凤春, 林雪玲, 曹志杰, 李小伏. Fe, Co, Ni掺杂GaSb的电子结构和光学性质. 物理学报, 2019, 68(18): 184202. doi: 10.7498/aps.68.20190290
    [5] 赵佰强, 张耘, 邱晓燕, 王学维. Cu,Fe掺杂LiNbO3晶体电子结构和光学性质的第一性原理研究. 物理学报, 2016, 65(1): 014212. doi: 10.7498/aps.65.014212
    [6] 潘凤春, 林雪玲, 陈焕铭. C掺杂金红石相TiO2的电子结构和光学性质的第一性原理研究. 物理学报, 2015, 64(22): 224218. doi: 10.7498/aps.64.224218
    [7] 骆最芬, 岑伟富, 范梦慧, 汤家俊, 赵宇军. BiTiO3电子结构及光学性质的第一性原理研究. 物理学报, 2015, 64(14): 147102. doi: 10.7498/aps.64.147102
    [8] 何静芳, 郑树凯, 周鹏力, 史茹倩, 闫小兵. Cu-Co共掺杂ZnO光电性质的第一性原理计算. 物理学报, 2014, 63(4): 046301. doi: 10.7498/aps.63.046301
    [9] 谢知, 程文旦. TiO2纳米管电子结构和光学性质的第一性原理研究. 物理学报, 2014, 63(24): 243102. doi: 10.7498/aps.63.243102
    [10] 程旭东, 吴海信, 唐小路, 王振友, 肖瑞春, 黄昌保, 倪友保. Na2Ge2Se5电子结构和光学性质的第一性原理研究. 物理学报, 2014, 63(18): 184208. doi: 10.7498/aps.63.184208
    [11] 程和平, 但加坤, 黄智蒙, 彭辉, 陈光华. 黑索金电子结构和光学性质的第一性原理研究. 物理学报, 2013, 62(16): 163102. doi: 10.7498/aps.62.163102
    [12] 王寅, 冯庆, 王渭华, 岳远霞. 碳-锌共掺杂锐钛矿相TiO2 电子结构与光学性质的第一性原理研究 . 物理学报, 2012, 61(19): 193102. doi: 10.7498/aps.61.193102
    [13] 杨春燕, 张蓉, 张利民, 可祥伟. 0.5NdAlO3-0.5CaTiO3电子结构及光学性质的第一性原理计算. 物理学报, 2012, 61(7): 077702. doi: 10.7498/aps.61.077702
    [14] 宋庆功, 刘立伟, 赵辉, 严慧羽, 杜全国. YFeO3的电子结构和光学性质的第一性原理研究. 物理学报, 2012, 61(10): 107102. doi: 10.7498/aps.61.107102
    [15] 李晓凤, 姬广富, 彭卫民, 申筱濛, 赵峰. 高压下固态Kr弹性性质、电子结构和光学性质的第一性原理计算. 物理学报, 2009, 58(4): 2660-2666. doi: 10.7498/aps.58.2660
    [16] 关丽, 李强, 赵庆勋, 郭建新, 周阳, 金利涛, 耿波, 刘保亭. Al和Ni共掺ZnO光学性质的第一性原理研究. 物理学报, 2009, 58(8): 5624-5631. doi: 10.7498/aps.58.5624
    [17] 林竹, 郭志友, 毕艳军, 董玉成. Cu掺杂的AlN铁磁性和光学性质的第一性原理研究. 物理学报, 2009, 58(3): 1917-1923. doi: 10.7498/aps.58.1917
    [18] 邢海英, 范广涵, 周天明. p,n型掺杂剂与Mn共掺杂GaN的电磁性质. 物理学报, 2009, 58(5): 3324-3330. doi: 10.7498/aps.58.3324
    [19] 毕艳军, 郭志友, 孙慧卿, 林 竹, 董玉成. Co和Mn共掺杂ZnO电子结构和光学性质的第一性原理研究. 物理学报, 2008, 57(12): 7800-7805. doi: 10.7498/aps.57.7800
    [20] 段满益, 徐 明, 周海平, 沈益斌, 陈青云, 丁迎春, 祝文军. 过渡金属与氮共掺杂ZnO电子结构和光学性质的第一性原理研究. 物理学报, 2007, 56(9): 5359-5365. doi: 10.7498/aps.56.5359
计量
  • 文章访问数:  329
  • PDF下载量:  29
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-12-15
  • 修回日期:  2022-01-13
  • 上网日期:  2022-01-28
  • 刊出日期:  2022-05-05

应变对(Ga, Mo)Sb磁学和光学性质影响的理论研究

  • 宁夏大学物理与电子电气工程学院, 银川 750021
  • 通信作者: 王旭明, wang_xm@126.com
    基金项目: 国家自然科学基金(批准号: 11764032, 11665018)资助的课题.

摘要: 近年来, 作为一种自旋电子学领域的关键材料, 具有高温本征铁磁性的稀磁半导体受到了广泛的关注. 为探索能够提高本征铁磁性居里温度(Curie temperature, TC)的方法, 本文运用第一性原理LDA+U方法研究了应变对Mo掺杂GaSb的电子结构、磁学及光学性质的影响. 研究结果表明: –6%—2.5%应变范围下GaSb半导体材料具有稳定的力学性能, 压应变下GaSb材料的可塑性、韧性增强, 有利于GaSb半导体材料力学性能的提升; 应变对Mo替代Ga缺陷(MoGa)的电子结构有重要的影响, –3%至–1.2%应变范围下MoGa处于低自旋态(low spin state, LSS), 具有1${\mu _{{\rm{B}}}}$的局域磁矩, –1.1%—2%应变范围下MoGa处于高自旋态(high spin state, HSS), 具有3${\mu _{{\rm{B}}}}$的磁矩; 不管是LSS还是HSS, MoGa产生局域磁矩之间的耦合都是铁磁耦合, 但铁磁耦合的强度和物理机制不同, 适当的压应变可有效提高铁磁耦合强度, 这有利于实现高TC的GaSb基磁性半导体; Mo可极大提高GaSb半导体材料的电极化能力, 这有利于光生电子-空穴对的形成和分离, 提高掺杂体系对长波光子的光电转化效率; Mo引入的杂质能级使电子的带间跃迁对所需要吸收光子的能量变小, 掺杂体系光学吸收谱的吸收边发生了红移, 拉应变可进一步提升(Ga,Mo)Sb体系在红外光区的光学性能.

English Abstract

    • 作为自旋电子学器件的关键材料, 磁性半导体同时利用了电子的电荷属性和自旋属性而具有本征的半导体性质和铁磁性, 已受到了广泛而持续的关注[1-7]. 其中III-V族磁性半导体很容易与III-V族非磁性半导体GaAs, AlAs, GaP, InP和InSb等结合形成异质结构, 呈现出自旋相关的散射、层间相互耦合作用、隧穿磁阻等现象. 有研究表明III-V族GaSb基磁性半导体具有可能高的TC、优异的电、磁、光等物理特性, 从而成为磁性半导体中研究的热门材料[8-11]. 第一过渡族金属(V, Cr, Mn, Fe, Co, Ni等)通常被选为III-V族和氧化物等磁性半导体的掺杂剂[12-14], 一方面是其3d电子壳层是部分占据的, 可提供局域磁矩, 另一方面这些过渡族金属元素的原子半径和金属阳离子的半径相匹配, 很容易形成替代掺杂. Dietl等[2]利用平均场理论预测当Mn掺杂含量和空穴浓度达到一定水平时, GaAs基磁性半导体的TC可以提高到室温以上. Sato等[15]从理论上预测Mn掺杂GaSb基磁性半导体的TC与Mn的掺杂浓度成正比. Tu等[16-18]在实验上证实Fe掺杂GaSb基磁性半导体的TC随着Fe掺杂浓度的增大而升高, 当Fe掺杂浓度达到25%时, (Ga, Fe)Sb的TC达到了340 K. 由此可见, 增大掺杂浓度是提高GaSb基磁性半导体TC的有效手段, 为了进一步提高过渡族金属的掺杂浓度, 在样品的制备过程中需要利用高温高压等极端条件对样品的生长氛围进行控制[19,20]. 然而, 过高的掺杂浓度容易在样品中形成过渡族金属团簇及过渡族金属化合物等磁性第二相[21,22], 因此有必要寻找另外的方法实现具有高TC和本征磁性的GaSb基磁性半导体. 从原子半径匹配的角度考虑, 第二过渡族金属如Zr, Nb, Mo等, 更容易替代半导体中的金属阳离子形成替代掺杂, 并且这些过渡族金属相是非磁性的, 排除了磁性半导体中磁性来源于金属团簇的可能性. Medvedeva[23]利用第一性原理计算的方法, 研究了Mo掺杂In2O3的电子结构和磁光特性, 揭示了磁相互作用在导电性以及控制光吸收的Burstein-Moss位移方面的重要作用. Park等[24,25]在Mo掺杂的In2O3半导体中观察到了室温铁磁性. Egbo等[26]研究发现Mo掺杂的In2O3体系位于导带底部电子的有效质量变小, 这提高了载流子的迁移率. Lu等[27]在空气中烧结的LaMn0.96Mo0.04O3中观察到的TC为238 K, 认为样品的铁磁性来源于由Mo6+替代Mn3+后诱发的Mn3+和Mn2+之间的双交换作用. Dwivedi等[28]第一次在Mo掺杂的CoFe2O4中发现了磁有序和铁电性共存的现象, 掺杂体系的磁化强度随着Mo掺杂浓度的增大而增强, 并把掺杂体系的巨介电常数归因于Maxwell-Wagner弛豫机制. 可以看出, Mo作为一种掺杂元素可以实现半导体的铁磁性, 并改善掺杂体系介电常数和光学吸收谱等光学性能. 最近有研究表明[29-34], 通过对半导体材料施加应变来调节其物理特性是一种可行的方法. Linpeng等[29]指出应变对受主束缚空穴自旋弛豫和退相具有重要的作用. 在双层膜中, 研究者们发现应变可以诱导材料从半导体态到金属态的转变[31]. (Ga0.8, Fe0.2)Sb材料中的磁晶各向异性表现出了对外延应变的依赖性, 并且当应变由拉应变变为压应变时其正负性发生了改变[32]. 通过选择不同衬底生长的半导体异质结构, 在其内部往往具有拉应变或压应变, 研究应变对半导体物理特性的影响具有重要的意义. 作为III-V族半导体材料, GaSb具有禁带宽度小、电子和空穴迁移率高等优点, 使其在半导体自旋电子器件、热光伏电池、红外激光器、红外探测器等领域有重要的应用前景. 在GaSb基磁性半导体领域, 存在掺杂元素种类少(主要为Mn, Fe掺杂), 磁矩起源及磁相互作用机制并没有得到充分研究, 研究结论比较浅显等显著问题. 基于以上事实, 本文研究了应变对Mo掺杂GaSb半导体的电子结构、磁学性质及光学性质的影响, 探讨应变及Mo掺杂对GaSb半导体物理特性影响的微观物理机制, 为实验上制备集磁、光功能于一体的新型半导体器件提供一定参考.

    • 本文采用基于密度泛函理论的CASTEP[35-37]软件进行研究. 计算体系为64个原子GaSb超晶胞, 包含32个Ga原子和32个Sb原子. 电子体系波函数采用平面波波函数展开, 平面波波函数截止能量为380 eV, 计算体系电荷密度和总能量在布里渊区积分进行, 基态能量采用Pulay密度混合法, 自洽精度为$5.0 \times {10^{ - 7}}~{{\rm{eV/atom}}}$, 交换关联泛函采用LDA-CA-PZ泛函, K空间网格点采用Monkhorst-Pack[38,39]方案, 数值取为$3 \times 3 \times 3$. 采用LDA+U方案来修正LDA/GGA交换关联泛函对GaSb半导体禁带宽度低估的这一问题, 对Ga-3d电子和Sb-5p电子的库仑位能进行修正, 当库仑位能的数值为$U_{{{\text{Ga-}}}3{{\rm{d}}}}^{} = 2.5\;{{\rm{eV}}}$${U_{{{\text{Sb-}}}5{{\rm{p}}}}} = 3.0\;{{\rm{eV}}}$时计算得到的结果同实验值符合最好, 此时得到的晶格常数${a_0} = 5.912$Å (实验值为6.096 Å, 1 Å = 0.1 nm), 禁带宽度Eg = 0.809 eV (实验值为0.810 eV)[40-42].

      Mo掺杂GaSb有两种基本的缺陷类型, 分别是MoGa和Mo替代Sb缺陷(MoSb). 通常用缺陷的形成能来衡量晶体中缺陷的浓度含量, 文中两种替代缺陷的形成能可由下式计算[43]:

      $ {E}_{{\rm{formation}}} = {E}_{{\rm{defect}}} - {E}_{{\rm{perfect}}} - {\mu }_{{\rm{Mo}}} + {\mu }_{X}, $

      式中${E_{{{\rm{defect}}}}}$表示包含一个MoGa(MoSb)的超晶胞的总能量, ${E_{{{\rm{perfect}}}}}$表示完好的GaSb超晶胞的总能量, ${\mu _{{{\rm{Mo}}}}}$表示Mo原子的化学势, ${\mu _X}$表示Ga原子或Sb原子的化学势. 缺陷浓度C与形成能${E_{{{\rm{formation}}}}}$的关系表达式为[43]

      $ C = {N}_{{\rm{sites}}}{N}_{{\rm{config}}}{{\rm{exp}}}\left[ - {E}_{{\rm{formation}}}/\left({k}_{{\rm{B}}}T\right)\right], $

      式中${N_{{{\rm{sites}}}}}$为单位体积晶格中缺陷位置的数目, ${N_{{{\rm{config}}}}}$为产生缺陷的等效数量, ${k_{{\rm{B}}}}$是玻尔兹曼常数, T是热力学温度. 由上式可知, 当T一定时, 缺陷浓度与形成能的大小有关, 某一缺陷的形成能越小, 相应的缺陷浓度就越大. 计算结果表明, MoGa的形成能是1.576 eV, MoSb的形成能是3.606 eV, 显然在热平衡状态下, MoGa的浓度高于MoSb. 实验上和理论上已经证实[17,18,44,45], MoGa为主的掺杂GaSb体系可以稳定存在, 因此本文主要对MoGa的电子结构、磁学和光学特性进行理论计算研究. 图1给出了包含64个原子的GaSb超晶胞结构示意图, 图中阿拉伯数字表示MoGa的位置, 包含MoGa的GaSb体系记为(Ga, Mo)Sb. 本文应变的定义为

      $ \varepsilon = \frac{{a - {a_0}}}{{{a_0}}} \times 100{\text{%}}, $

      式中$a$${a_0}$分别表示有应变时和无应变时的晶格常数, $\varepsilon > 0$表示拉应变, $\varepsilon < 0$表示压应变.

      图  1  GaSb超晶胞结构, 其中大的绿色球代表Ga原子, 小的紫色球代表Sb原子

      Figure 1.  Structure of GaSb supercell, where the big green balls and small purple balls denote Ga and Sb atoms, respectively.

    • 材料结构的稳定性、刚性、韧性等力学特性取决于材料对外力的响应, 因此研究材料的弹性性质是研究材料本构关系的重要内容. 材料的体积模量(bulk modulus, B)和剪切模量(shear modulus, G)通常有两种不同的计算方法, 分别是Voigt提出的晶体边界上的应力连续模型[46] 和Reuss提出的晶体边界上的应变连续模型[47,48]. Hill[49]证明了Voigt模型和Reuss模型的计算结果分别对应晶体弹性系数的上下限, 并提出了晶体弹性系数取Voigt模型和Reuss模型计算结果算术平均值的Voigt-Reuss-Hill模型. 在弹性变化范围内, 晶体的应力和应变成正比, 弹性模量矩阵是对称矩阵, 满足${C_i}_j = {C_{ji}}$. GaSb半导体属于对称性最高的立方晶系, 独立变化的张量元素个数为3个, 分别为${C_{11}}$, ${C_{44}}$${C_{12}}$. Voigt-Reuss-Hill模型下体积模量和剪切模量的计算公式及力学稳定性判据(criteria, CRIT)如下:

      $ B = ({C_{11}} + 2{C_{12}})/3, $

      $ \begin{split} G = \dfrac{{C}_{11} - {C}_{12} + 3{C}_{44}}{10} + \dfrac{5({C}_{11} - {C}_{12}){C}_{44}}{8{C}_{44} + 6({C}_{11} - {C}_{12})}, \end{split} $

      $ {{\text{CRIT 1:}}}\begin{array}{*{20}{c}} {} \end{array}~ {C_{11}} > 0, $

      $ {{\text{CRIT 2:}}}\begin{array}{*{20}{c}} {} \end{array}~ {C_{44}} > 0, $

      $ {{\text{CRIT 3:}}}\begin{array}{*{20}{c}} {} \end{array}~ {C_{11}} - \left| {{C_{12}}} \right| > 0. $

      图2给出了–14%—5%应变区间力学稳定性判据的变化趋势, 可以看出当应变范围为–6%—2.5%时, GaSb半导体具有稳定的力学性能. 表1列出了–6%—2%应变范围下GaSb半导体的体积模量、剪切模量、泊松比和B/G的变化数据.

      图  2  力学稳定性判据

      Figure 2.  Criteria of mechanical stability.

      Strain/%B/GPaG/GPa$\gamma $B/G
      –6148.2137.670.383.93
      –5129.6643.080.353.01
      –3.7109.9343.850.322.51
      –290.0343.300.292.08
      068.7747.210.221.46
      0.568.7347.170.221.47
      164.2645.470.211.48
      1.561.3540.380.201.52
      255.4938.020.201.53

      表 1  应变下GaSb的体积模量B、剪切模量G、泊松比$\gamma $B/G

      Table 1.  Bulk modulus B, shear modulus G, Poisson ratio $\gamma $ and B/G under strains.

      体积模量反映了材料在弹性系下对外界均一性压缩的抵抗能力. 由表1可以看出, GaSb材料的体积模量随着压应变的增大而增大, 随着拉应变的增大而减小. 晶体的体积模量与组成晶体化学键的伸缩能力有关, 微观水平下晶体的体积模量由化学键的强度和压缩性来决定[50]. 计算结果表明, 随着压应变的增大, Ga原子的失电子能力逐渐减弱, Ga—Sb化学键的键布居逐渐增大, 这说明Ga—Sb化学键的离子性减弱, 共价性增强. 从另一个角度看, 原子抓电子能力越强, 价电子云在外力作用下越难发生移动, 因此一定程度上可以认为化学键的共价性越强, 材料对外界压缩的抵抗力也就越强. 同时从表1可以看出, 剪切模量在零应变状态下的数值最大为47.21, 且随着拉应变或压应变的增大而变小, 表示材料在应变下抗剪切形变的能力逐渐减弱, 可塑性增强.

      泊松比表示纵向应力所引起的横向应变与相应纵向应变之比的绝对值, 是反映材料横向变形的弹性常数. 从表1可以看出, GaSb材料的泊松比$\gamma $数值随着压应变的增大而逐渐增大. 按照断裂行为的判据[51], 泊松比高的材料具有高的韧性, 因此可认为随着压应变的增大, GaSb半导体材料的韧性增强, 发生脆性断裂的可能性小. 根据Pugh判据[52], B/G可以表示材料的塑性变形能力, 从表1可以看出, B/G随着压应变或拉应变的增大而逐步升高, 说明应变下GaSb半导体材料的延展性变好.

      综上所述, –6%—2.5%应变范围内, GaSb半导体材料具有稳定的力学性能. 在此应变范围内, 随着压应变的增大, GaSb材料的抗压性增强, 抵御塑性形变的能力降低, 材料具有较好的可塑性. 泊松比的数值表明, 压应变下GaSb半导体材料的韧性增强, 发生脆性断裂的可能性较小. 可以看出, 适当的压应变有利于GaSb半导体材料力学性能的提升, 这对GaSb半导体在极端环境下保持物理性能的稳定性是有利的.

      半导体材料在制备的过程中, 一般可通过材料设计, 包括引入缺陷[53,54], 构建失配的异质结构[55,56]等方法在材料的内部引入压应变或拉应变. 王志伟等[57]利用等离子体增强化学气相沉积法在晶格失配较大的GaSb衬底上沉积SiO2薄膜, 通过改变薄膜沉积时的工艺条件, 如反应温度、射频功率、反应压强等获得了具有不同内部应变的SiO2薄膜. Goel等[32]利用分子束外延法在不同缓冲层上制备了具有不同应变的(Ga, Fe)Sb磁性半导体材料($\varepsilon = - 1.7$%, 缓冲层AlSb; $\varepsilon = 0.23$%, 缓冲层In0.5Ga0.5As; $\varepsilon = 3.84$%, 缓冲层GaAs), 并认为通过不同的缓冲层可生长出具有不同应变的(Ga, Fe)Sb磁性半导体材料. Kondrin 等[19]在实验上研究了(Ga, Cr)Sb体系在高压(6—8 GPa)下的电子输运特性和磁特性, 第一性原理计算结果表明, 当施加静水压力为8 GPa时, GaSb体系的应变可以达到$\varepsilon = - 3.5$%, 可见对样品施加高压也是使其内部产生压应变的有效手段. 基于半导体材料晶格常数的考虑, 可认为通过以AlSb, CdSe, InSb等(InAs, InP, AlAs等)作为缓冲层可制备出具有不同拉应变(压应变)的(Ga, Mo)Sb半导体材料. 考虑到太大的拉应变或压应变在材料制备的过程难以实现, 本文主要研究–3%—2%应变范围内(Ga, Mo)Sb体系的磁学及光学性能.

    • 研究了不同应变下包含一个MoGa的(Ga, Mo)Sb体系的电子结构和磁特性. 计算结果表明, 当应变在–1.1%—2%范围内变化时, (Ga, Mo)Sb体系的态密度图(density of states, DOS)和磁矩的分布具有相似性. 图3(a)图3(b)给出了零应变下(Ga, Mo)Sb体系的总态密度图(total DOS, TDOS)和Mo-4d电子的投影态密度图(projected DOS, PDOS). 从图3(a)可以看出, 零应变下(Ga, Mo)Sb体系的自旋向上电子和自旋向下电子的态密度在零点费米能级处不对称, 发生了自旋劈裂, 说明(Ga, Mo)Sb体系中产生了净磁矩. 图3(b)给出了Mo-4d电子的PDOS, Mo-4d电子占据自旋向上态的电子数目多于占据自旋向下态的电子数目, 表明(Ga, Mo)Sb体系中的Mo原子产生了自旋向上(正)的净磁矩. 图4(a)给出了不同应变下(Ga, Mo)Sb体系中磁矩的分布来源, 计算结果显示, 零应变下Mo原子产生3.04${\mu _{{\rm{B}}}}$的磁矩, 近邻Sb原子产生–0.13${\mu _{{\rm{B}}}}$的磁矩, 其余磁矩由次近邻Ga原子所贡献, 一个MoGa产生总磁矩的大小为3${\mu _{{\rm{B}}}}$, 可以看出磁矩主要由掺杂Mo原子所贡献. 为了形象展示磁矩的这一分布, 图5(a)给出了–1%应变下包含一个MoGa的(Ga, Mo)Sb体系的空间自旋密度分布图, 图中的蓝色部分表示由Mo原子产生向上的净自旋, 黄色部分是Mo近邻4个Sb原子产生向下的净自旋. –1.1%—2%应变范围内, 一个MoGa产生总的磁矩都为3${\mu _{{\rm{B}}}}$, 只是磁矩的分布有所不同, 见图4(a), 此时MoGa处于高自旋态(high spin state, HSS).

      图  3  (Ga, Mo)Sb体系在HSS和LSS下的态密度图 (a) 零应变下的TDOS; (b)零应变下Mo-4d电子的PDOS; (c) –1.2%应变下的TDOS; (d) –1.2%应变下Mo-4d电子的PDOS

      Figure 3.  Density of states (DOS) of (Ga, Mo)Sb in HSS and LSS: (a) TDOS and (b) PDOS of Mo-4d without strain; (c) the TDOS and (d) the PDOS of Mo-4d under –1.2% strain.

      图  4  (a) 应变下, (Ga, Mo)Sb体系中掺杂Mo原子和近邻Sb原子的磁矩; (b) 应变下, S(12)和S(23)结构的ΔE

      Figure 4.  (a) Magnetic moments contributed by Mo and neighbor Sb (N Sb) under strains and (b) ΔE of S(12) and S(23) under strains, respectively.

      图  5  等值面为0.01 e/Å3的空间自旋密度图 (a) HSS下MoGa; (b) LSS下MoGa; (c) HSS下S(12)结构FM耦合; (d) LSS下S(23)结构FM耦合

      Figure 5.  The spatial spin density for (a) MoGa under HSS, (b) MoGa under LSS, (c) FM coupling of S(12) under HSS, (d) FM coupling of S(23) under LSS, respectively. The isovalue is set to 0.01 e/Å3.

      当应变在–3%至–1.2%范围变化时, MoGa转变为低自旋态(low spin state, LSS), 此时(Ga, Mo)Sb体系产生的总磁矩为1${\mu _{{\rm{B}}}}$. 从图3(c)可以看出(Ga, Mo)Sb体系的TDOS自旋向上电子占据数目和自旋向下电子占据数目不相等, 表明体系中产生了净磁矩. 图3(d)给出了Mo-4d电子的PDOS, 同样表明掺杂体系中Mo原子产生了正的净磁矩. 比较图3(b)图3(d)可以看出, HSS下Mo-4d电子在零点费米能级处的一个自旋向上的态密度峰是部分占据的, 而LSS下Mo-4d电子在零点费米能级处有一个电子部分占据的自旋向下的态密度峰. 从图4(a)可以看出, –3%至–1.2%应变范围下(Ga, Mo)Sb体系的磁矩主要由Mo原子所贡献, Mo原子近邻Sb原子产生自旋向下的净磁矩, 且数值很小. LSS下MoGa的空间自旋密度分布图, 见图5(b), 可以看出, Mo原子周围4个近邻Sb原子产生的磁矩很小.

      显然, 掺杂的Mo原子在一定程度上诱导了GaSb半导体的自旋极化, 应变对MoGa产生的局域磁矩大小和分布都有影响. MoGa中, Mo原子周围最近邻的4个Sb原子构成正四面体, Mo原子位于正四面体的中心. 以正立方体为参照系, 坐标原点位于正立方体的中心, 坐标沿着与正立方体的棱平行的方向伸展, 见图6. 4个最近邻Sb原子位于正立方体互不相邻的4个顶点上, dxy, dxz和dyz这3条d电子轨道的波瓣指向正立方体各棱的中心, 距离4个Sb原子较近, 受到较大的排斥作用使轨道能量升高较多, 分裂后这3条d电子轨道的能量简并, 用群论的符号记为${{{\rm{t}}}_2}$轨道. 另外两个${{{\rm{d}}}_{{z^2}}}$${{{\rm{d}}}_{{x^2} - {y^2}}}$轨道的波瓣指向正立方体的面心, 距离4个Sb原子较远, 受到的排斥力较小, 轨道能量升高较少, 分裂后这两条轨道的能量简并, 记为e轨道. 四面体晶场作用下d电子轨道的分裂能即为${{{\rm{t}}}_2}$轨道和e轨道之间的能量差, 记为${{{\varDelta}} _{{\rm{t}}}}$.

      图  6  d电子轨道四面体晶场分裂示意图

      Figure 6.  Diagram of Mo-4d orbital splitting under tetrahedral crystal field.

      Mo原子的电子组态是[Kr]4d55s1, 当Mo原子替代Ga原子后, Mo原子失去3个价电子与周围近邻的4个Sb原子成键, 此时Mo离子的电子组态变为[Kr]4d3. Mo-4d电子轨道在四面体晶场的作用下分裂成能量较低的e轨道和能量较高的${{{\rm{t}}}_2}$轨道, Mo剩余的3个4d电子的排布有两种可能: 2个电子排在自旋向上的e轨道, 1个电子排在自旋向下的e轨道, 电子结构可以表示为${{{\rm{e}}}^{ \uparrow \uparrow \downarrow }}$($ \uparrow $表示自旋向上的电子, $ \downarrow $表示自旋向下的电子), 这种方式称为LSS, 此时电子成对要克服成对能${E_{{\rm{P}}}}$; 为了避免成对能, 2个电子在自旋向上的e轨道, 1个电子在自旋向上的${{{\rm{t}}}_2}$轨道, 电子结构可以表示为${{{\rm{e}}}^{ \uparrow \uparrow }}{{\rm{t}}}_2^ \uparrow $, 这种方式称为HSS. 显然, 当${{{\varDelta}} _{{\rm{t}}}} > {E_{{\rm{P}}}}$时电子排布采取LSS, 反之采取HSS.

      –1.1%—2%应变范围下, Mo-4d电子的轨道分裂能${{{\varDelta}} _{{\rm{t}}}} < {E_{{\rm{P}}}}$, Mo-4d电子处于HSS, 此时(Ga, Mo)Sb体系的能带结构图见图7(a). 从图7(a)可以看出, ${{\rm{t}}}_{2 + }^{}$(“+”表示自旋向上的轨道)轨道跨过零点费米能级(蓝色虚线), 对应图3(b)中零点费米能级处被电子部分占据的自旋向上态. –3%—–1.2%应变范围下, Mo-4d电子的轨道分裂能${{{\varDelta}} _{{\rm{t}}}} > $$ {E_{{\rm{P}}}}$, Mo-4d电子处于LSS, 此时(Ga, Mo)Sb体系的能带结构图见图7(b). 从图7(b)可以看出, ${{{\rm{e}}}_ - }$(“–“表示自旋向下的轨道)轨道跨过零点费米能级, 对应图3(d)中零点费米能级处被电子部分占据的自旋向下态. Mo-4d电子从HSS到LSS的转变, 是由于不同应变下Mo原子和近邻Sb原子的紧密程度不同所导致. 应变由2%至–3%变化的过程中, Mo原子和4个近邻的Sb原子越来越紧密, 受到四面体晶场的影响越来越大, d电子轨道的分裂能也就越来越大, 当应变达到–1.2%时, 分裂能${{{\varDelta}} _{{\rm{t}}}} > {E_{{\rm{P}}}}$, 这是Mo-4d电子由HSS向LSS转变的临界应变.

      图  7  (Ga, Mo)Sb体系的能带图 (a)零应变下HSS; (b) –1.2%应变下LSS. 黑色实线表示自旋向上的能带, 红色虚线表示自旋向下的能带. 在能量为0处的水平蓝色虚线表示费米能级

      Figure 7.  Band structures of (Ga, Mo)Sb in (a) HSS without strain and (b) LSS under –1.2% strain. Black solid lines and red dotted lines denote spin-up and spin-down bands, respectively. The horizontal blue dotted lines located at 0 energy are Fermi levels.

      为了研究MoGa产生局域磁矩之间的耦合特性, 在包含64原子GaSb超晶胞结构中引入两个MoGa, 为了避免第一性原理计算中镜像缺陷之间不合理的相互作用对计算结果的影响, 两个MoGa之间的距离不应该超过掺杂体系晶格常数的一半. 设计了两种结构: 用两个Mo原子替代图1中的1号和2号Ga原子, 记为结构S(12); 用两个Mo原子替代2号和3号Ga原子形成结构S(23). 零应变下, S(12)和S(23)结构中两个Mo原子之间的距离分别为4.180 Å和5.912 Å, 其中S(23)结构中两个MoGa都在(010)晶面内. 设定两种结构下Mo原子的初始自旋, 就可以计算铁磁(ferromagnetism, FM)耦合和反铁磁(anti ferromagnetism, AFM)耦合状态下体系的总能量, MoGa的磁耦合特性可以根据FM状态下体系总能量(${E_{{{\rm{FM}}}}}$)和AFM状态下体系总能量(${E_{{{\rm{AFM}}}}}$)的能量差来判断, $\Delta E = {E_{{{\rm{AFM}}}}} - {E_{{{\rm{FM}}}}}$, 若$\Delta E > 0$, 表明MoGa产生的局域磁矩倾向于FM耦合, 反之为AFM耦合.

      图4(b)给出了$\Delta E$随应变的变化趋势, $\Delta E$可以表征局域磁矩之间的耦合强度和相应磁序的稳定性. 从图4(b)可以看出, S(12)和S(23)两种结构在不同应变下, $\Delta E > 0$, $\Delta E$越大表明体系FM耦合相互作用越强, TC越高[58]. 图4(b)的结果表明, 无论是LSS还是HSS, 掺杂体系局域磁矩之间的耦合都倾向于FM耦合, 其FM耦合机制可以用d-d交换模型和p-d交换模型来解释[59]. –3%至–1.2%应变范围下MoGa处于LSS, 其电子结构为${{{\rm{e}}}^{ \uparrow \uparrow \downarrow }}$, ${{{\rm{e}}}_{{\rm{ - }}}}$轨道跨过零点费米能级, 见图3(d)图7(b), 不同MoGa${{{\rm{e}}}_{{\rm{ - }}}}$轨道上电子之间的交换并不需要改变电子的自旋方向, 因此LSS下MoGa产生局域磁矩之间的耦合是FM耦合, 这种由Mo-4d电子直接交换产生的相互作用称为d-d交换作用. –1.1%—2%应变范围下MoGa处于HSS, 其电子结构为${{{\rm{e}}}^{ \uparrow \uparrow }}{{\rm{t}}}_2^ \uparrow $, ${{\rm{t}}}_{2 + }^{}$轨道跨过零点费米能级, 见图3(b)图7(a), ${{\rm{t}}}_{2 + }^{}$轨道上电子之间的交换作用(d-d交换作用)导致MoGa产生的局域磁矩之间的耦合同样是FM耦合, 虽然HSS下e轨道是半占据的, 但由于${{{\rm{e}}}_{{\rm{ + }}}}$轨道位于零点费米能级下–1.25 eV处, 见图7(a), 因此${{{\rm{e}}}_{{\rm{ + }}}}$轨道上的电子并不会产生直接的交换作用.

      图4(b)中的两条$\Delta E$曲线展现了左低相差较大和右高相差较小的特点. 从图4(b)可以看出, –3%至–1.2%应变范围下, S(12)结构的$\Delta E$(黑色方块线)较小, –1.2%应变下的数值为122.3 meV, 当应变增大到–1%时, $\Delta E$跃变到212.6 meV. 对比图4(a)可以看出, MoGa也由LSS转变到HSS. 对于S(23)结构, –3%至–1.2%应变范围内的$\Delta E$(红色圆点线)低于同应变范围内S(12)结构的$\Delta E$, –1.2%应变下的数值只有75.8 meV, 当应变增大到–1%时, $\Delta E$跃变到179.6 meV, 此时S(23)结构中MoGa由LSS转变到HSS. S(12)和S(23)两种结构在LSS下(–3%至–1.2%)的$\Delta E$均明显小于其在HSS下(–1%—2%)的$\Delta E$, 这是因为LSS下局域磁矩之间相互作用只能通过d-d交换作用来实现, 但对于HSS状态, 除了d-d交换作用, 还有Mo-4d电子通过近邻Sb-5p电子产生的间接交换作用, 这种交换作用称为p-d交换作用. 图4(a)的计算结果表明, LSS下掺杂Mo原子近邻Sb-p电子产生的磁矩(–0.05${\mu _{{\rm{B}}}}$)远小于HSS下近邻Sb-p电子产生的磁矩(–0.12${\mu _{{\rm{B}}}}$), 因此HSS下局域磁矩之间的FM耦合是d-d交换作用和p-d交换作用共同支配的结果, HSS下FM耦合相互作用更强一些, 这导致了较大的$\Delta E$. 对于LSS, 相同应变下S(12)和S(23)两种结构体系的$\Delta E$相差较大, 这是由d-d交换作用的短程性的特点所导致[4], S(23)结构下两个MoGa的距离较远, d-d交换作用的强度迅速减弱. 对于HSS, 相同应变下S(12)和S(23)两种结构体系的$\Delta E$变化相对较小, 这是因为HSS下除了d-d交换作用, 还有p-d交换作用, 而p-d交换作用具有长程性[4], 因此HSS下MoGa的间距对其产生局域磁矩之间耦合强度的影响没有LSS下那样明显. 值得关注的是, 在–1.1%—0%应变范围内, Mo处于HSS, 此时局域磁矩之间的FM耦合强度明显大于零应变时的情况, 说明在此应变范围内对实现高TC本征铁磁性(Ga, Mo)Sb半导体是极其有利的.

      对于III-V族磁性半导体各种磁现象的解释, Dietl等[2,60]建立了空穴(hole, h)载流子导致局域化自旋磁性相互作用的平均场模型. 从图4(a)可以看出, –1.2%应变下Mo产生的磁矩为1.03${\mu _{{\rm{B}}}}$, –1%应变下Mo产生的磁矩为3.01${\mu _{{\rm{B}}}}$, 因此(Ga, Mo)Sb体系中Mo原子价态为+3价, 此时Mo-4d电子数目为3个. 根据Hund定则, 半满的d壳层(5个d电子)是比较稳定的, 因此MoGa中心将束缚2个电子在4d壳层, 再加上2个弱束缚的空穴构成的状态比较稳定, 此时MoGa的电子态为A0(d5+2h), 这里A0表示电中性中心, 这样Mo原子替代3价Ga原子, 既提供了局域化的自旋, 又提供了空穴载流子, 由于p-d杂化, Mo杂质带附近的价带空穴自旋与Mo局域自旋方向相反. 由于价带波函数比较扩展, 一个空穴可以与一定数目的Mo局域自旋结合, 当空穴浓度足够高时, 就可以产生以空穴载流子为媒介的自旋与自旋之间的相互作用. 在(Ga, Mn)As磁性半导体中, Mn的价态为+3价, Mn替代Ga缺陷(MnGa)中心的电子态为A0(d5+h)[61,62], 可以看出, 相比于Mn原子, Mo原子在增大空穴载流子浓度方面更具有优势, 这对于提高GaSb基磁性半导体的TC是极为有利的.

      图5(c)给出了HSS下S(12)结构两个MoGa产生局域磁矩FM耦合的自旋密度空间分布图. 可以看出磁矩之间的相互作用是通过Mo-4d电子波函数的交叠(d-d交换作用)和Sb-5p电子波函数的交叠(p-d交换作用)来实现. 比较图5(a)图5(c)可知, 图5(c)中的两个MoGa均处于HSS. 图5(d)给出了LSS下S(23)结构两个MoGa产生局域磁矩FM耦合的自旋密度空间分布图, LSS下Sb-5p电子对磁矩的贡献很小, MoGa产生磁矩之间的相互作用主要通过d-d交换作用而非p-d交换作用来实现. 比较图5(b)图5(d)可知, 图5(d)中两个MoGa均处于LSS.

    • GaSb由于其较低的禁带宽度, 使其可以与多种辐射体的光谱相匹配, 因此以GaSb材料为基础的热光伏电池得到了广泛的研究[63,64]. 为了使GaSb基热光伏电池更好地与温度较低的辐射器相匹配, 基于GaSb相关的三元、四元III-V族化合物被发现并用来制作热光伏电池. 虽然多元GaSb相关化合物的禁带宽度最低达到了0.53 eV, 有效拓展了红外光谱的响应范围, 但由于较大的表面复合速度而限制了热电转换效率的进一步提升[65-67]. 目前用于热光伏电池的GaSb基半导体材料仍然存在禁带宽度偏大, 不能很好地和热辐射器相匹配以及实际转换效率不高等问题, 因此有必要研究(Ga, Mo)Sb体系的光学性质, 拓展GaSb基半导体材料在红外热光伏电池、红外半导体激光器、探测器等领域的应用.

      半导体在线性响应范围的光学性质通常用复介电函数$ \varepsilon (\omega ) = {\varepsilon _1}(\omega ) + {\rm{i}} {\varepsilon _2}(\omega ) $来描述, 其中${\varepsilon _2}(\omega )$表示复介电函数的虚部, 由价电子在占据轨道和非占据轨道之间的跃迁来计算; ${\varepsilon _1}(\omega )$表示复介电函数的实部, 由${\varepsilon _2}(\omega )$所满足的Kramers-Krönig色散关系确定[68]. 具体函数表达式为

      $ {\varepsilon }_{1}(\omega )=1+\frac{2}{{\rm{π}}}p{\displaystyle \int _{0}^{\infty }\frac{s{\varepsilon }_{2}(s)}{{s}^{2}-{\omega }^{2}}}{\rm{d}}s, $

      $ {\varepsilon }_{2}(\omega )=-\frac{2\omega }{{\rm{π}}}p{\displaystyle \int _{0}^{\infty }\frac{{\varepsilon }_{1}(s)}{{s}^{2}-{\omega }^{2}}}{\rm{d}}s, $

      式中$p$为积分主值, $\omega $为入射光子的圆频率, s为积分变量. 积分主值p的表达式为

      $ p\int_0^\infty {} = \mathop {\lim }\limits_{\delta \to 0} \left(\int_0^{\omega - \delta } {} + \int_{\omega + \delta }^\infty {} \right). $

      利用复介电函数与复折射率之间的函数关系(${(n(\omega ) + {{\rm{i}}}k(\omega ))^2} = {\varepsilon _1}(\omega ) +{\rm{i}} {\varepsilon _2}(\omega )$)和Kramers-Krönig色散关系可以得到吸收系数$\alpha (\omega )$[68], 函数表达式如下:

      $ \alpha (\omega ) = 2\omega k(\omega )/c $

      式中c为光速, 普通折射率$n(\omega )$为复折射率的实部, 决定光衰减的$k(\omega )$为复折射率的虚部.

      图8(a)给出了应变下GaSb和(Ga, Mo)Sb的光学吸收谱. 可以看出GaSb体系(黑色实线)对远红外光区光子没有吸收, 其光学吸收谱的吸收边落在了1667 nm处, 此波长光子对应的能量是电子从价带顶跃迁到导带底所需要的最小能量, 即0.809 eV. –2%—2%应变范围内, (Ga, Mo)Sb体系在大于300 nm范围内对入射光子的吸收幅度相比于未掺杂GaSb体系均有很大的提升, 且其光学吸收谱的吸收边落在了远红外光区, 即发生了红移现象. 从图7可以看出, Mo掺杂引入的杂质能级位于价带和导带之间, 且跨过零点费米能级, 电子由价带向导带跃迁时以杂质能级为桥梁, 跃迁所需要的能量变小, 因此提升了(Ga, Mo)Sb体系对长波光子的吸收能力, 这是(Ga, Mo)Sb体系光学吸收谱的吸收边发生红移的原因. 从图8(a)还可以看出, (Ga, Mo)Sb体系对可见光区和红外光区光子的吸收幅度随着拉应变的增大而增大, 随着压应变的增大而减小. 这说明应变对(Ga, Mo)Sb体系的光吸收也有影响, 适当的拉应变有利于提升(Ga, Mo)Sb体系对可见光区和红外光区光子的吸收幅度.

      图  8  (a) 应变下, GaSb和(Ga, Mo)Sb的光学吸收谱; (b) 应变下, GaSb和(Ga, Mo)Sb的复介电函数, 实线代表实部, 虚线代表虚部

      Figure 8.  (a) Absorption spectra and (b) complex dielectric function of GaSb and (Ga, Mo)Sb under strains. In panel (b) solid and dotted lines denote real and imaginary parts, respectively.

      复介电函数的实部${\varepsilon _1}(\omega )$表征了半导体材料在外电场作用下的极化程度, 实部数值越大说明体系的电极化能力越强, 在入射光子能量为零时对应的数值为静态介电常数. 从图8(b)可以看出, 零应变下GaSb复介电函数实部(黑色实线)在入射光子能量为0时的数值为19.43, 零应变下(Ga, Mo)Sb体系(绿色实线)所对应的数值为32.84. 这说明MoGa可以有效提高(Ga, Mo)Sb体系的介电性能, 掺杂体系电极化能力的提升和光生电场强度的变大有利于光生电子-空穴对的形成和分离, 提高了掺杂体系对长波光子的光电转换能力. –2%—2%应变范围内, (Ga, Mo)Sb体系在0—1.75 eV区间的实部数值均大于零应变GaSb体系所对应实部数值, 且(Ga, Mo)Sb体系的静态介电常数随着拉应变的增大而增大, 随着压应变的增大而减小. 这同样说明应变对(Ga, Mo)Sb体系的光电转换效率有重要影响, 适当的拉应变有利于(Ga, Mo)Sb体系光学性能的提升.

      复介电函数的虚部${\varepsilon _2}(\omega )$表征的是半导体内部形成电偶极子所消耗的能量, 虚部数值越大说明价电子处于激发态的数目越多, 反映了半导体材料中电子的受激跃迁程度. 从图8(b)可以看出, –2%—2%应变范围内, (Ga, Mo)Sb体系在0—1.51 eV能量区间的虚部数值均大于零应变下GaSb体系(黑色虚线)对应能量区间的虚部数值, 表明(Ga, Mo)Sb体系对红外、远红外光区光子的吸收能力增强. 零应变(Ga, Mo)Sb体系的虚部(绿色虚线)在光子能量为0的数值为0.87, 而GaSb体系对应虚部的数值为0. 这说明(Ga, Mo)Sb体系对能量为“0”的光子就有响应, 其原因在于Mo引入的杂质能级是部分占据的且跨过零点费米能级, 电子在简并的杂质能级之间跃迁并不需要吸收光子即可发生, 此即电子的带内响应, 虽然Mo掺杂引入杂质能级上电子的带内响应对掺杂体系的光学性质有一定的影响, 起主要作用的是电子在不同能级之间的直接跃迁[69]. –2%—2%应变范围内, (Ga, Mo)Sb体系在低能区的数值随着拉应变的增大而增大, 随着压应变的增大而降低, 说明拉应变可以提升(Ga, Mo)Sb体系电子的受激跃迁程度, 有利于材料光催化性能的提升.

      (Ga, Mo)Sb体系光学性质部分的计算模型为包含一个MoGa的64个原子GaSb超晶胞结构体系, 此时掺杂Mo原子的摩尔浓度为3.125%. 若考虑(Ga, Mo)Sb体系磁性半导体具有较高的TC, Mo原子的掺杂浓度不应低于6.25%, 这也正是S(23)结构中Mo原子的掺杂浓度. 虽然掺杂浓度高, 费米能级处杂质能级上电子的数目多, 有利于(Ga, Mo)Sb体系对红外光区、远红外光区光子的吸收; 但也有其不利的一面, 由掺杂元素引入杂质能级的附加势是作用距离仅为一两个原子的短程势[68], 高浓度掺杂下密集的Mo原子可以成为电子-空穴对的有效复合中心, 从而使载流子的寿命大大降低, 影响材料的光学性能. 采用低温分子束外延生长方法, 浓度高达25%的(Ga, Fe)Sb样品实验上已制备成功[16,17], 理论上关于铁磁性计算的掺杂浓度高达50%[11,70], 若综合考虑(Ga, Mo)Sb体系的磁学和光学性质, 可认为在均匀掺杂前提下, 提高掺杂浓度既可以增大体系中的空穴载流子浓度, 又可以避免电子-空穴对复合中心的形成[69,71].

      以上分析结果表明, MoGa的掺杂有利于提高GaSb半导体材料中载流子的受激跃迁程度, 同时掺杂体系的电极化能力和光生电场强度的增强有利于光生电子-空穴对的产生和分离, 改善了(Ga, Mo)Sb体系的光催化特性, 提升了(Ga, Mo)Sb体系对可见光区、红外光区光子的吸收幅度, 拉应变可进一步提高(Ga, Mo)Sb体系的光学性能.

    • 运用第一性原理LDA+U的方法研究了不同应变下Mo掺杂GaSb的电子结构、磁学和光学性质. 结论如下: 1) –6%—2.5%应变范围内GaSb半导体材料具有稳定的力学性能. 随着压应变的增大, GaSb材料的抗压性和可塑性增强, 材料发生脆性断裂的可能性较小, 适当的压应变有利于GaSb半导体材料力学性能的提高, 这对GaSb半导体材料在极端环境下保持稳定的物理性能是有利的. 2) 应变对MoGa的电子结构有重要的影响. –3%至–1.2%应变范围下MoGa处于LSS并产生1${\mu _{{\rm{B}}}}$的局域磁矩, 局域磁矩之间的耦合是FM耦合, 强度较弱, FM耦合的微观物理机制可用d-d交换作用模型来解释. –1.1%—2%应变范围下MoGa处于HSS并产生3${\mu _{{\rm{B}}}}$的磁矩, 可用d-d交换作用模型和p-d交换作用模型来解释局域磁矩之间的FM耦合, HSS下磁矩之间的FM耦合具有很高的稳定性. HSS下的FM耦合强度远大于LSS下的FM耦合强度, 对于HSS, $\Delta E$随着压应变增大而变大, –1%应变下的$\Delta E$高达212.6 meV, 表明对于(Ga, Mo)Sb体系适当的压应变可以使MoGa处于HSS, 且MoGa产生局域磁矩之间的FM耦合强度也会明显增强. 3) MoGa引入的杂质能级可以作为电子从价带向导带跃迁的桥梁, 电子的带间跃迁对所需要吸收光子的能量变小, 从而提升了(Ga, Mo)Sb体系对长波光子的吸收能力, 光学吸收谱的吸收边落在了远红外光区, 有效拓展了(Ga, Mo)Sb体系对光子的吸收范围; 零应变下, MoGa使GaSb半导体的静态介电常数由19.43提升至32.84, 表明MoGa的引入可以大大提高(Ga, Mo)Sb体系的电极化能力, 有利于(Ga, Mo)Sb体系光生电子-空穴对的产生和分离; (Ga, Mo)Sb体系复介电函数虚部数值的增大表明其价电子具有很高的受激跃迁程度. MoGa提高了(Ga, Mo)Sb体系对长波光子的吸收及其光电转换效率, 提高了材料的光催化特性. 此外, 拉应变可进一步提升(Ga, Mo)Sb在红外光区的光学性能.

      通过对应变下(Ga, Mo)Sb体系磁学和光学性质影响的研究分析, 我们认为在样品的生长过程中, 施加一定的应变可提升(Ga, Mo)Sb体系的磁学及光学性能, 期望我们的研究为实验上制备新型的GaSb基半导体材料提供一定的理论参考.

参考文献 (71)

目录

    /

    返回文章
    返回