搜索

x

留言板

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

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

具有大磁晶各向异性能的单层BaPb的室温量子反常霍尔效应

张蔚曦 李勇 田昌海 佘彦超

引用本文:
Citation:

具有大磁晶各向异性能的单层BaPb的室温量子反常霍尔效应

张蔚曦, 李勇, 田昌海, 佘彦超

Room-temperature quantum anomalous Hall effect in monolayer BaPb with large magnetocrystalline anisotropy energies

Zhang Wei-Xi, Li Yong, Tian Chang-Hai, She Yan-Chao
PDF
HTML
导出引用
  • 具有巨大应用潜力的二维材料在纳米技术领域引起了人们极大的研究兴趣. 基于第一性原理计算, 本文预测了二维六角晶格BaPb体系具有室温量子反常霍尔效应. 体系磁基态是铁磁半金属态, 并且自旋极化的Pb-p轨道导致体系的时间反演对称性破缺. 具有非零Chern数(C = 1)的单层BaPb中非平庸拓扑性来源于全自旋极化的px,y轨道形成的二次型非狄拉克能带. 不同于之前报道的pz轨道形成的狄拉克拓扑态很容易被衬底破坏, px,y轨道形成的σ键非常稳定. 当考虑自旋-轨道耦合作用时, 二次型的非狄拉克点打开了接近177.39 meV的非平庸带隙. 通过反常霍尔电导、陈数、贝里曲率和边缘态的计算, 证实了BaPb非平庸的拓扑性. 此外, 体系还表现出每个单元52.01 meV的大磁晶各向异性能.
    The quantum anomalous Hall effect is an intriguing quantum state that exhibits chiral edge states in the absence of a magnetic field. The chiral edge states are topologically protected and robust against electron scattering, which possesses great potential applications in designing low energy consumption and dissipation less spintronic devices. The experimental conditions are required to be very high, such as extremely low temperature (< 100 mK) due to the small band gap and the greatly accurate control of the extrinsic impurities. These greatly hinder their devices from being put into applications further. Hence, it would be meaningful to search for a new Chern insulator with a large band gap and high Curie temperature. According to the first-principles calculations, we predict the room temperature quantum anomalous Hall effect in the monolayer BaPb. The nontrivial topology of this new type of ferroelectric semi-metal material derives from fully spin-polarized quadratic non-Dirac bands. The quantum anomalous Hall effect can be realized in the monolayer BaPb with fully spin-polarized quadratic px,y non-Dirac bands with the nonzero Chern number (C = 1). Because of the trigonal symmetry of monolayer BaPb material, these bands composed of px,y orbitals are at the $ \varGamma $ point, which is different from the Dirac state formed by the pz orbital reported previously. In addition, it can still retain its original topological properties even if strongly hybridized with the substrate. The calculated phonon spectrum shows no imaginary frequency in the entire Brillouin zone, indicating that the monolayer BaPb system is dynamically stable. By using Monte Carlo simulation, we determine the Curie temperature of BaPb monolayer toreach up to 378 K. We also calculate the magnetic anisotropy energy of the BaPb cell, defined as $ \Delta E={E_{100}}-{E_{001}} $. Here, we consider two magnetization easy-axis directions, [100] and [001]. To our surprise, the MAE of monolayer BaPb is as high as 52.01 meV/cell by considering the spin-orbit coupling effect. Furthermore, the nontrivial band gap is opened with a magnitude of 177.39 meV when the spin-orbit coupling effect is included. The calculations of Berry curvature and edge states further prove that the monolayer BaPb system can realize the quantum anomalous Hall state. This discovery indicates that the monolayer BaPb materials can be used as a candidate for quantum anomalous Hall effect materials, thereby promoting the development of spintronics.
      通信作者: 佘彦超, ycshe@xtu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 12064038)、贵州省科研项目(批准号: KY[2019]179, KY[2017]315, ZK[2021]一般–034)、贵州省优秀青年科技人才项目(批准号: [2019]5673)和铜仁市科技局科学技术基金(批准号[2020]77)资助的课题
      Corresponding author: She Yan-Chao, ycshe@xtu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 12064038), the Scientific Research Program of Guizhou Province, China (Grant Nos. KY[2019]179, KY[2017]315, TK[2021] General projects -034), the Outstanding Young Science and Technology Talents of Guizhou Province, China (Grants No. [2019]5673), and the Science and Technology Foundation of Tongren Science and Technology Bureau, China (Grant No. [2020]77)
    [1]

    Lu J M, Zheliuk O, Leermakers I, Yuan N F Q, Zeitler U, Lwa K T, Ye J T 2015 Science 350 1353Google Scholar

    [2]

    Gong C, Zhang X 2019 Science 363 eaav4450

    [3]

    Deng Y, Yu Y, Shi M Z, Guo Z, Xu Z, Wang J, Chen X H, Zhang Y 2020 Science 367 895

    [4]

    高艺璇, 张礼智, 张余洋, 杜世萱 2018 物理学报 67 238101Google Scholar

    Gao Y X, Zhang L Z, Zhang Y Y, Du S X 2018 Acta Phys. Sin. 67 238101Google Scholar

    [5]

    徐依全, 王聪 2020 物理学报 69 184216Google Scholar

    Xu Y Q, Wang C 2020 Acta Phys. Sin. 69 184216Google Scholar

    [6]

    孙真昊, 管鸿明, 付雷, 沈波, 唐宁 2021 物理学报 70 027302

    Sun Z H, Guan H M, Fu L, Shen B, Tang N 2021 Acta Phys. Sin. 70 027302

    [7]

    Haldane F D M 1988 Phys. Rev. Lett. 61 2015Google Scholar

    [8]

    Liu C X, Qi X L, Dai X, Fang Z, Zhang S C 2008 Phys. Rev. Lett. 101 146802Google Scholar

    [9]

    Zhang H, Lazo C, Blügel S, Heinze S, Mokrousov Y 2012 Phys. Rev. Lett. 108 056802Google Scholar

    [10]

    Hu J, Zhu Z, Wu R 2015 Nano Lett. 15 2074Google Scholar

    [11]

    Li P, Li X, Zhao W, Chen H, Chen M X, Guo Z X, Feng J, Gong X G, MacDonald A H 2017 Nano Lett. 17 6195Google Scholar

    [12]

    Li P 2019 Phys. Chem. Chem. Phys. 21 6712Google Scholar

    [13]

    Onoda M, Nagaosa N 2003 Phys. Rev. Lett. 90 206601Google Scholar

    [14]

    Hasan M Z, Kane C L 2010 Rev. Mod. Phys. 82 3045Google Scholar

    [15]

    Luo W, Xiang H 2015 Nano Lett. 15 3230Google Scholar

    [16]

    Checkelsky J G, Yoshimi R, Tsukazaki A, Takahashi K S, Kozuka Y, Falson J, Kawaski M, Tokura Y 2014 Nat. Phys. 10 731Google Scholar

    [17]

    Bestwick A J, Fox E J, Kou X, Pan L, Wang K L, Goldhaber-Gordon D 2015 Phys. Rev. Lett. 114 187201Google Scholar

    [18]

    Chang C Z, Zhao W, Kim D Y, Zhang H, Assaf B A, Heiman D, Zhang S C, Liu C, Chan M H W, Moodera J S 2015 Nat. Mater. 14 473Google Scholar

    [19]

    Ou Y, Liu C, Jiang G, Feng Y, Zhao D, Wu W, Wang X X, Li W, Song C, Wang L L, Wang W, Wu W, Wang Y, He K, Ma X C, Xue Q K 2018 Adv. Mater. 30 1703062Google Scholar

    [20]

    Otrokov M M, Rusinov I P, Blanco-Rey M, Hoffmann M, Vyazovskaya A Y, Eremeev S V, Ernst A, Echenique P M, Arnau A, Chulkov E V 2019 Phys. Rev. Lett. 122 107202Google Scholar

    [21]

    Zhang D, Shi M, Zhu T, Xing D, Zhang H, Wang J 2019 Phys. Rev. Lett. 122 206401Google Scholar

    [22]

    Otrokov M M, Klimovskikh I I, Chulkov E V 2019 Nature 576 416Google Scholar

    [23]

    Wang S, Yu Z, Liu Y, Jiao Y, Guan S, Sheng X, Yang S A 2019 Phys. Rev. Materials 3 084201Google Scholar

    [24]

    Kong X, Li L, Leenaerts O, Wang W, Liu X, Peeters F M 2018 Nanoscale 10 8153Google Scholar

    [25]

    Gong C, Li L, Li Z, Ji H, Stern A, Xia Y, Cao T, Bao W, Wang C, Wang Y, Qiu Z Q, Cava R J, Louie S G, Xia J, Zhang X 2017 Nature 546 265Google Scholar

    [26]

    Huang B, Clark G, Navarro-Moratalla E, Klein D R, Cheng R, Seyler K L, Zhong D, Schmidgall E, McGuire M A, Cobden D H, Yao W, Xiao D, Jarillo-Herrero P, Xu X 2017 Nature 546 270Google Scholar

    [27]

    Zhang Z, Shang J, Jiang C, Rasmita A, Gao W, Yu T 2019 Nano Lett. 19 3138Google Scholar

    [28]

    Kresse G, Hafner J 1993 Phys. Rev. B 47 558Google Scholar

    [29]

    Blochl P E 1994 Phys. Rev. B 50 17953Google Scholar

    [30]

    Dudarev S L, Botton G A, Savrasov S Y, Humphreys C J, Sutton A P 1998 Phys. Rev. B 57 1505Google Scholar

    [31]

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

    [32]

    Mostofi A A, Yates J R, Lee Y S, Souza I, Vanderbilt D, Marzari N 2008 Comput. Phys. Commun. 178 685Google Scholar

    [33]

    Wu Q S, Zhang S N, Song H F, Troyer M, Suluyanov A A 2018 Comput. Phys. Commun. 224 405Google Scholar

    [34]

    Zhuang H L, Henning R G 2013 Appl. Phys. Lett. 103 212102Google Scholar

    [35]

    Evans R F, Fan W J, Chureemart P, Ostler T A, Ellis M O, Chantrell R W 2014 J. Phys. Condens. Matter. 26 103202Google Scholar

    [36]

    Thouless D J, Kohmoto M, Nightingale M P, den Nijs M 1982 Phys. Rev. Lett. 49 405Google Scholar

    [37]

    Zhang R W, Zhang C W, Ji W X, Yan S S, Yao Y G 2017 Nanoscale 9 8207Google Scholar

    [38]

    Reis F, Li G, Dudy L, Bauernfeind M, Glass S, Hanke W, Thomale R, Schafer J, Claessen R 2017 Science 357 287Google Scholar

  • 图 1  (a) 单层BaPb的晶格结构; (b) 单层BaPb声子谱. 计算中使用一个4 × 4的超胞

    Fig. 1.  (a) Lattice structure for monolayer BaPb; (b) the phonon dispersion curves monolayer BaPb. A 4 × 4 supercell is used in the calculation.

    图 2  单分子层BaPb的铁磁构型(a)和反铁磁构型(b); (c)利用蒙特卡罗模拟得到的单层BaPb的磁矩与温度的关系

    Fig. 2.  Two magnetic configurations for monolayer BaPb: (a) Ferromagnetic; (b) antiferromagnetic. (c) The average magnetic moment per unit cell with respect to temperature calculated for monolayer BaPb obtained by using Monte Carlo simulations.

    图 3  单分子层BaPb的自旋极化能带结构

    Fig. 3.  Spin-polarized band structure for monolayer BaPb.

    图 4  考虑SOC时单层BaPb的能带结构 (a) 磁矩平面内沿x方向; (b) 磁矩平面外沿z方向

    Fig. 4.  Band structures (including SOC) for monolayer BaPb: (a) With magnetic moments oriented in-plane along x; (b) with magnetic moments oriented out of the plane (along z directions).

    图 5  (a) 反常霍尔电导率随费米能级移动的变化; (b) 动量空间中带有SOC的Berry曲率; (c) 计算得到的单层BaPb半无限片的边缘状态

    Fig. 5.  (a) Anomalous Hall conductivity when the Fermi level is shifted around the original Fermi level; (b) the Berry curvature with SOC in the momentum space; (c) calculated edge state of a semi-infinite sheet for monolayer BaPb.

    图 6  (a), (b) h-BN/BaPb异质结(h-BN/BaPb)体系中单分子层BaPb的结构晶体结构的俯视图和侧视图; (c) h-BN/BaPb异质结自旋极化能带结构; (d) 考虑SOC效应时h-BN/BaPb异质结体系中单分子层BaPb的带结构(考虑SOC效应)

    Fig. 6.  Crystal structure of monolayer BaPb in the h-BN/BaPb heterostructure system: (a) Top and (b) side views. Here, the black solid lines denote the unit cell. (c) Spin-polarized band structure of h-BN/BaPb heterostructure; (d) the band structures (including SOC) of monolayer BaPb in the h-BN/BaPb heterostructure system.

    图 7  利用蒙特卡罗模拟得到的h-BN/BaPb异质结的磁矩与温度的关系

    Fig. 7.  Average magnetic moment per unit cell with respect to temperature calculated for the h-BN/BaPb heterojunction obtained by using Monte Carlo simulation.

  • [1]

    Lu J M, Zheliuk O, Leermakers I, Yuan N F Q, Zeitler U, Lwa K T, Ye J T 2015 Science 350 1353Google Scholar

    [2]

    Gong C, Zhang X 2019 Science 363 eaav4450

    [3]

    Deng Y, Yu Y, Shi M Z, Guo Z, Xu Z, Wang J, Chen X H, Zhang Y 2020 Science 367 895

    [4]

    高艺璇, 张礼智, 张余洋, 杜世萱 2018 物理学报 67 238101Google Scholar

    Gao Y X, Zhang L Z, Zhang Y Y, Du S X 2018 Acta Phys. Sin. 67 238101Google Scholar

    [5]

    徐依全, 王聪 2020 物理学报 69 184216Google Scholar

    Xu Y Q, Wang C 2020 Acta Phys. Sin. 69 184216Google Scholar

    [6]

    孙真昊, 管鸿明, 付雷, 沈波, 唐宁 2021 物理学报 70 027302

    Sun Z H, Guan H M, Fu L, Shen B, Tang N 2021 Acta Phys. Sin. 70 027302

    [7]

    Haldane F D M 1988 Phys. Rev. Lett. 61 2015Google Scholar

    [8]

    Liu C X, Qi X L, Dai X, Fang Z, Zhang S C 2008 Phys. Rev. Lett. 101 146802Google Scholar

    [9]

    Zhang H, Lazo C, Blügel S, Heinze S, Mokrousov Y 2012 Phys. Rev. Lett. 108 056802Google Scholar

    [10]

    Hu J, Zhu Z, Wu R 2015 Nano Lett. 15 2074Google Scholar

    [11]

    Li P, Li X, Zhao W, Chen H, Chen M X, Guo Z X, Feng J, Gong X G, MacDonald A H 2017 Nano Lett. 17 6195Google Scholar

    [12]

    Li P 2019 Phys. Chem. Chem. Phys. 21 6712Google Scholar

    [13]

    Onoda M, Nagaosa N 2003 Phys. Rev. Lett. 90 206601Google Scholar

    [14]

    Hasan M Z, Kane C L 2010 Rev. Mod. Phys. 82 3045Google Scholar

    [15]

    Luo W, Xiang H 2015 Nano Lett. 15 3230Google Scholar

    [16]

    Checkelsky J G, Yoshimi R, Tsukazaki A, Takahashi K S, Kozuka Y, Falson J, Kawaski M, Tokura Y 2014 Nat. Phys. 10 731Google Scholar

    [17]

    Bestwick A J, Fox E J, Kou X, Pan L, Wang K L, Goldhaber-Gordon D 2015 Phys. Rev. Lett. 114 187201Google Scholar

    [18]

    Chang C Z, Zhao W, Kim D Y, Zhang H, Assaf B A, Heiman D, Zhang S C, Liu C, Chan M H W, Moodera J S 2015 Nat. Mater. 14 473Google Scholar

    [19]

    Ou Y, Liu C, Jiang G, Feng Y, Zhao D, Wu W, Wang X X, Li W, Song C, Wang L L, Wang W, Wu W, Wang Y, He K, Ma X C, Xue Q K 2018 Adv. Mater. 30 1703062Google Scholar

    [20]

    Otrokov M M, Rusinov I P, Blanco-Rey M, Hoffmann M, Vyazovskaya A Y, Eremeev S V, Ernst A, Echenique P M, Arnau A, Chulkov E V 2019 Phys. Rev. Lett. 122 107202Google Scholar

    [21]

    Zhang D, Shi M, Zhu T, Xing D, Zhang H, Wang J 2019 Phys. Rev. Lett. 122 206401Google Scholar

    [22]

    Otrokov M M, Klimovskikh I I, Chulkov E V 2019 Nature 576 416Google Scholar

    [23]

    Wang S, Yu Z, Liu Y, Jiao Y, Guan S, Sheng X, Yang S A 2019 Phys. Rev. Materials 3 084201Google Scholar

    [24]

    Kong X, Li L, Leenaerts O, Wang W, Liu X, Peeters F M 2018 Nanoscale 10 8153Google Scholar

    [25]

    Gong C, Li L, Li Z, Ji H, Stern A, Xia Y, Cao T, Bao W, Wang C, Wang Y, Qiu Z Q, Cava R J, Louie S G, Xia J, Zhang X 2017 Nature 546 265Google Scholar

    [26]

    Huang B, Clark G, Navarro-Moratalla E, Klein D R, Cheng R, Seyler K L, Zhong D, Schmidgall E, McGuire M A, Cobden D H, Yao W, Xiao D, Jarillo-Herrero P, Xu X 2017 Nature 546 270Google Scholar

    [27]

    Zhang Z, Shang J, Jiang C, Rasmita A, Gao W, Yu T 2019 Nano Lett. 19 3138Google Scholar

    [28]

    Kresse G, Hafner J 1993 Phys. Rev. B 47 558Google Scholar

    [29]

    Blochl P E 1994 Phys. Rev. B 50 17953Google Scholar

    [30]

    Dudarev S L, Botton G A, Savrasov S Y, Humphreys C J, Sutton A P 1998 Phys. Rev. B 57 1505Google Scholar

    [31]

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

    [32]

    Mostofi A A, Yates J R, Lee Y S, Souza I, Vanderbilt D, Marzari N 2008 Comput. Phys. Commun. 178 685Google Scholar

    [33]

    Wu Q S, Zhang S N, Song H F, Troyer M, Suluyanov A A 2018 Comput. Phys. Commun. 224 405Google Scholar

    [34]

    Zhuang H L, Henning R G 2013 Appl. Phys. Lett. 103 212102Google Scholar

    [35]

    Evans R F, Fan W J, Chureemart P, Ostler T A, Ellis M O, Chantrell R W 2014 J. Phys. Condens. Matter. 26 103202Google Scholar

    [36]

    Thouless D J, Kohmoto M, Nightingale M P, den Nijs M 1982 Phys. Rev. Lett. 49 405Google Scholar

    [37]

    Zhang R W, Zhang C W, Ji W X, Yan S S, Yao Y G 2017 Nanoscale 9 8207Google Scholar

    [38]

    Reis F, Li G, Dudy L, Bauernfeind M, Glass S, Hanke W, Thomale R, Schafer J, Claessen R 2017 Science 357 287Google Scholar

  • [1] 秦文静, 徐波, 孙宝珍, 刘刚. 原子吸附的二维CrI3铁磁半导体电学和磁学性质的第一性原理研究. 物理学报, 2021, 70(11): 117101. doi: 10.7498/aps.70.20210090
    [2] 罗娅, 张耘, 梁金铃, 刘林凤. 铜铁镁三掺铌酸锂晶体的第一性原理研究. 物理学报, 2020, 69(5): 054205. doi: 10.7498/aps.69.20191799
    [3] 王靖. 手征马约拉纳费米子. 物理学报, 2020, 69(11): 117302. doi: 10.7498/aps.69.20200534
    [4] 梁金铃, 张耘, 邱晓燕, 吴圣钰, 罗娅. 铁镁共掺钽酸锂晶体的第一性原理研究. 物理学报, 2019, 68(20): 204205. doi: 10.7498/aps.68.20190575
    [5] 陈国祥, 樊晓波, 李思琦, 张建民. 碱金属和碱土金属掺杂二维GaN材料电磁特性的第一性原理计算. 物理学报, 2019, 68(23): 237303. doi: 10.7498/aps.68.20191246
    [6] 杨彪, 王丽阁, 易勇, 王恩泽, 彭丽霞. C, N, O原子在金属V中扩散行为的第一性原理计算. 物理学报, 2015, 64(2): 026602. doi: 10.7498/aps.64.026602
    [7] 黄艳平, 袁健美, 郭刚, 毛宇亮. 硅烯饱和吸附碱金属原子的第一性原理研究. 物理学报, 2015, 64(1): 013101. doi: 10.7498/aps.64.013101
    [8] 胡洁琼, 谢明, 张吉明, 刘满门, 杨有才, 陈永泰. Au-Sn金属间化合物的第一性原理研究. 物理学报, 2013, 62(24): 247102. doi: 10.7498/aps.62.247102
    [9] 王爱玲, 毋志民, 王聪, 胡爱元, 赵若禺. 新型稀磁半导体Mn掺杂LiZnAs的第一性原理研究. 物理学报, 2013, 62(13): 137101. doi: 10.7498/aps.62.137101
    [10] 赵立凯, 赵二俊, 武志坚. 5d过渡金属二硼化物的结构和热、力学性质的第一性原理计算. 物理学报, 2013, 62(4): 046201. doi: 10.7498/aps.62.046201
    [11] 杨理践, 刘斌, 高松巍, 陈立佳. 金属磁记忆效应的第一性原理计算与实验研究. 物理学报, 2013, 62(8): 086201. doi: 10.7498/aps.62.086201
    [12] 吴木生, 徐波, 刘刚, 欧阳楚英. 应变对单层二硫化钼能带影响的第一性原理研究 . 物理学报, 2012, 61(22): 227102. doi: 10.7498/aps.61.227102
    [13] 袁娣, 罗华锋, 黄多辉, 王藩侯. Zn,O共掺杂实现p型AlN的第一性原理研究. 物理学报, 2011, 60(7): 077101. doi: 10.7498/aps.60.077101
    [14] 罗礼进, 仲崇贵, 全宏瑞, 谭志中, 蒋青, 江学范. Heusler合金Mn2NiGe磁性形状记忆效应的第一性原理预测. 物理学报, 2010, 59(11): 8037-8041. doi: 10.7498/aps.59.8037
    [15] 孙源, 明星, 孟醒, 孙正昊, 向鹏, 兰民, 陈岗. 多铁材料BaCoF4电子结构的第一性原理研究. 物理学报, 2009, 58(8): 5653-5660. doi: 10.7498/aps.58.5653
    [16] 卢志鹏, 祝文军, 刘绍军, 卢铁城, 陈向荣. 非静水压条件下铁从α到ε结构相变的第一性原理计算. 物理学报, 2009, 58(3): 2083-2089. doi: 10.7498/aps.58.2083
    [17] 杨银堂, 武 军, 蔡玉荣, 丁瑞雪, 宋久旭, 石立春. p型K:ZnO导电机理的第一性原理研究. 物理学报, 2008, 57(11): 7151-7156. doi: 10.7498/aps.57.7151
    [18] 张金奎, 邓胜华, 金 慧, 刘悦林. ZnO电子结构和p型传导特性的第一性原理研究. 物理学报, 2007, 56(9): 5371-5375. doi: 10.7498/aps.56.5371
    [19] 丁少锋, 范广涵, 李述体, 肖 冰. 氮化铟p型掺杂的第一性原理研究. 物理学报, 2007, 56(7): 4062-4067. doi: 10.7498/aps.56.4062
    [20] 姚红英, 顾 晓, 季 敏, 张笛儿, 龚新高. SiO2-羟基表面上金属原子的第一性原理研究. 物理学报, 2006, 55(11): 6042-6046. doi: 10.7498/aps.55.6042
计量
  • 文章访问数:  951
  • PDF下载量:  54
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-01-04
  • 修回日期:  2021-03-17
  • 上网日期:  2021-06-07
  • 刊出日期:  2021-08-05

具有大磁晶各向异性能的单层BaPb的室温量子反常霍尔效应

  • 铜仁学院物理与电子工程系, 铜仁 554300
  • 通信作者: 佘彦超, ycshe@xtu.edu.cn
    基金项目: 国家自然科学基金(批准号: 12064038)、贵州省科研项目(批准号: KY[2019]179, KY[2017]315, ZK[2021]一般–034)、贵州省优秀青年科技人才项目(批准号: [2019]5673)和铜仁市科技局科学技术基金(批准号[2020]77)资助的课题

摘要: 具有巨大应用潜力的二维材料在纳米技术领域引起了人们极大的研究兴趣. 基于第一性原理计算, 本文预测了二维六角晶格BaPb体系具有室温量子反常霍尔效应. 体系磁基态是铁磁半金属态, 并且自旋极化的Pb-p轨道导致体系的时间反演对称性破缺. 具有非零Chern数(C = 1)的单层BaPb中非平庸拓扑性来源于全自旋极化的px,y轨道形成的二次型非狄拉克能带. 不同于之前报道的pz轨道形成的狄拉克拓扑态很容易被衬底破坏, px,y轨道形成的σ键非常稳定. 当考虑自旋-轨道耦合作用时, 二次型的非狄拉克点打开了接近177.39 meV的非平庸带隙. 通过反常霍尔电导、陈数、贝里曲率和边缘态的计算, 证实了BaPb非平庸的拓扑性. 此外, 体系还表现出每个单元52.01 meV的大磁晶各向异性能.

English Abstract

    • 在过去的十几年里, 二维材料由于其独特的物理化学性质, 在许多领域均有着广泛的应用, 例如纳米电子和光器件、能量存储、催化和化学传感器等[1-6]. Haldane[7]基于二维蜂窝状晶格提出了量子反常霍尔效应(quantum anomalous Hall effect, QAHE). 这一理论预测激发了研究者们对非耗散电子器件应用的极大兴趣[8-11]. 然而, 实现QAHE的难点在于既要满足时间反演对称性的破坏, 又要满足自旋-轨道耦合(spin-orbit coupling, SOC) [12-15]. 因此, 建立长程铁磁有序的体系成为最常规的设计思路. 目前, 该方法只能使QAHE在极低温条件下(通常低于100 mK)观察到, 因为长程磁序会被磁畴破坏[16-18]. 随后, V-I或Cr-V共掺杂和三掺杂的拓扑绝缘体系的掺杂方案被提出, 以提高QAHE观测温度, 但改善效果不明显[19-22]. 极端温度条件严重限制了自旋电子器件的实际应用. 利用QAHE在室温下寻找一些非常规磁性材料是非常重要的[23,24]. 幸运的是, 研究人员发现, 碱土金属体系作为一种独特的铁磁半金属, 可能为解决这一问题提供了一种潜在的方法. 这种材料的磁性来源于阴离子p轨道的自旋极化. 这与传统的磁性材料有很大的不同. 最近发现的二维铁磁绝缘体也提供了前所未有的机会来创造能谷劈裂, 从而在高温下更薄的体系中实现QAHE[25-27].

      本文基于密度泛函理论(density functional theory, DFT)预测室温量子反常霍尔效应可以在单层BaPb材料中实现. 这种新型铁磁半金属材料的非平庸拓扑性来自于全自旋极化的二次非狄拉克带. 由于单层BaPb材料的C3v对称性, 因此由px,y轨道组成的位于Γ点的非线性能带, 这与之前报道的pz轨道形成的线性狄拉克态不同, 即使与衬底杂化, 它仍能保持原有的拓扑性质. 计算得到的声子谱在整个布里渊区内没有虚频, 说明单分子层BaPb体系是动力学稳定的. 通过蒙特卡罗模拟, 确定了单层BaPb的居里温度可达378 K. 通过定义式$\Delta E={E_{100}}-{E_{001}}$, 考虑了两种磁化方向: 平面内的[100]方向和垂直于平面的[001]方向, 计算了BaPb晶胞的磁晶各向异性能(magnetic anisotropy energy, MAE). 令人惊讶的是, 考虑SOC效应, 单分子层BaPb的MAE高达52.01 meV/cell. 同时, 二次型的非狄拉克点打开177.39 meV的带隙. 贝利曲率和边缘态的计算进一步证明了单层BaPb系统可以实现量子反常霍尔态. 这一发现表明单层BaPb材料可以作为QAHE材料的候选材料, 从而促进了自旋电子学的发展.

    • 为了系统地研究单层BaPb的磁性、电子性质和拓扑性质, 采用了基于DFT的VASP软件包(the Vienna ab initio simulation package)[28,29]. 选取了广义梯度近似的Perdew-Burke-Ernzerh (PBE) 势作为交换关联势. 计算中使用了投影缀加平面波(projected augmented wave, PAW)方法处理电子与原子核之间的相互作用. 在结构优化中, 为了找到最稳定的结构, 离子驰豫的收敛标准设置为10–6 eV, 平面波截断能采用了500 eV, 并使用了Γ-centered方法产生K点, K点网格为12 × 12 × 1. 此外, 为了防止薄膜与重复单元薄膜直接的相互作用, 在垂直于平面向外的z方向, 设置了20 Å的真空层. 本文使用PHONOPY软件包的有限位移法来计算声子谱, 以保证结构的稳定性[30]. 利用半无限边界条件下格林函数的方法和Wannier Tools软件包对单层BaPb的边缘态、贝利曲率和反常霍尔电导进行了计算[31]. 对于贝利曲率积分, 采用了K点网格为120 × 120 × 1[32,33].

    • 优化后的单层BaPb结构如图1(a)所示. 原胞包含一个Ba原子和一个Pb原子, 形成类似石墨烯的平整的六角形蜂窝状结构. 其空间群为$P\bar 6 m2$, 点群为${D_{3 h}}$. 体系优化后晶格常数为6.01 Å, Ba—Pb键长为3.47 Å. 本文从形成能和结构的动力学稳定性两个方面分析了单层BaPb的稳定性. 首先通过定义${E_{\rm{f}}} = {E_{{\rm{2 D}}}}/{N_{{\rm{2 D}}}} - {E_{{\rm{3 D}}}}/{N_{{\rm{3 D}}}}$计算了形成能[34]. 式中${E_{2{\rm{D}}}}$ (${E_{{\rm{3 D}}}}$) 和${N_{{\rm{2 D}}}}$ (${N_{{\rm{3 D}}}}$)分别是二维(三维)结构的总能和各自原胞内所包含的原子数. 计算得到该结构的形成能为1.03 eV/atom, 大于众所周知的硅烯的0.76 eV/atom[35,36]. 从单分子层BaPb的声子谱可以推断其动力学稳定性. 随后采用4 × 4的超胞计算了其声子谱, 如图1(b)所示. 可以看到, 单层BaPb结构声子谱的所有振动模式在整个布里渊区均为正值, 不存在虚频, 表明其动力学是稳定的.

      图  1  (a) 单层BaPb的晶格结构; (b) 单层BaPb声子谱. 计算中使用一个4 × 4的超胞

      Figure 1.  (a) Lattice structure for monolayer BaPb; (b) the phonon dispersion curves monolayer BaPb. A 4 × 4 supercell is used in the calculation.

    • 通过DFT计算发现, 无自旋极化的单层BaPb的总能量远高于自旋极化的单层BaPb. 这意味着尽管单分子层BaPb系统不包含磁性原子, 但它仍然具有磁性. 为了解系统的磁性来源, 如图2(a)图2(b)所示, 计算了$2 \times 2$超胞的三种磁构型: 铁磁(ferromagnetic, FM)态、反铁磁(antiferromagnetic, AFM)态和非铁磁(non-magnetic, NM)态的能量, 以确定单层BaPb的磁基态.

      图  2  单分子层BaPb的铁磁构型(a)和反铁磁构型(b); (c)利用蒙特卡罗模拟得到的单层BaPb的磁矩与温度的关系

      Figure 2.  Two magnetic configurations for monolayer BaPb: (a) Ferromagnetic; (b) antiferromagnetic. (c) The average magnetic moment per unit cell with respect to temperature calculated for monolayer BaPb obtained by using Monte Carlo simulations.

      计算表明, 对于单层BaPb体系, FM态能量比NM态和AFM态分别低981.26和198.95 meV. 这表明, FM态具有最低能量, 体系基态为铁磁态, 磁矩主要来源于Pb原子的p轨道上. 此外, 考虑了两种磁化方向: 平面内的[100]方向和垂直于平面的[001]方向计算了BaPb原胞的MAE. MAE定义为$\Delta E={E_{100}}-{E_{001}}$, 其中${E_{100}}$${E_{001}}$分别表示磁化沿[100]和[001]方向的体系的总能. MAE为正值说明易磁化方向沿着z轴方向, 反之则是沿着x轴方向. 计算结果约为52.01 meV, 表明磁矩沿着[001]轴时能量最低. 因此, 将重点研究BaPb单层磁化沿z轴的拓扑性质.

      体系的铁磁居里温度对于QAHE的实验观测极为重要. 因此, 还利用蒙特卡罗(Monte Carlo, MC)模拟, 估计了BaPb单层的居里温度(${T_{\rm{c}}}$)[35]. 本文用海森伯模型来描述自旋哈密顿量, 即$H = $$ - \dfrac{1}{2}\displaystyle\sum\nolimits_{i, j} {{J_{ij}}{{\boldsymbol{S}}_i}} \cdot {{\boldsymbol{S}}_j} - K\sum\nolimits_i {{{(S_i^z)}^2}}$. 第一项和第二项分别表示到最近邻的短程交换作用(J )和体系的MAE (K ). 通过DFT计算可以得出交换作用参数J值为16.58 meV. 在这里, 本文的MC模拟是考虑周期边界条件在一个100 × 100的超胞上进行的, 从而得出单层BaPb的居里温度为378 K, 如图2(c)所示.

    • 单层BaPb最有趣的性质仍然反映在其电子结构上. 首先研究不考虑SOC作用的单层BaPb的能带结构. 如图3所示, 单分子层BaPb的能带结构是完全自旋极化的, 即, 自旋向上的能带展现出绝缘性, 而自旋向下的能带展现出金属性. 同时, 可以清楚地观测到自旋向下的能带在费米能级附近的Γ点上存在二次型的非狄拉克能带.

      图  3  单分子层BaPb的自旋极化能带结构

      Figure 3.  Spin-polarized band structure for monolayer BaPb.

      随后, 研究考虑SOC作用的磁矩分别沿x方向平面内和沿z方向平面外的单层BaPb的能带结构, 如图4所示. 结果表明, 磁矩沿z方向, 二次型的非狄拉克能带打开了全局域带隙177.39 meV, 体系呈现出拓扑绝缘态, 而在平面内沿x方向有磁矩的带隙结构则显示出金属态. 这意味着通过施加电场或构造异质结改变磁化方向, 可以实现体系从金属态到拓扑绝缘态的相变.

      图  4  考虑SOC时单层BaPb的能带结构 (a) 磁矩平面内沿x方向; (b) 磁矩平面外沿z方向

      Figure 4.  Band structures (including SOC) for monolayer BaPb: (a) With magnetic moments oriented in-plane along x; (b) with magnetic moments oriented out of the plane (along z directions).

    • 在这一部分, 通过计算贝利曲率和Chern数来研究单层BaPb的拓扑性质. 首先, 利用公式${\sigma _{xy}} = C{e^2}/h$计算出反常霍尔电导, 如图5(a)所示. Chern数C [36]可以通过对第一布里渊区的倒空间的贝利曲率积分$C \!= \!\dfrac{1}{{2{\text{π}}}}\!\displaystyle\int_{BZ} \!{{{\rm{d}}^2}k\varOmega (k)}$得到, 其中

      图  5  (a) 反常霍尔电导率随费米能级移动的变化; (b) 动量空间中带有SOC的Berry曲率; (c) 计算得到的单层BaPb半无限片的边缘状态

      Figure 5.  (a) Anomalous Hall conductivity when the Fermi level is shifted around the original Fermi level; (b) the Berry curvature with SOC in the momentum space; (c) calculated edge state of a semi-infinite sheet for monolayer BaPb.

      $ \varOmega ({\boldsymbol{k}}) \!=\! - \displaystyle\sum\limits_n {{f_n}} \sum\limits_{n' \ne n} \frac{2 {\rm Im} \langle \psi _{nk} | \nu _x |\psi _{n'k}\rangle \langle \psi _{nk} |\nu _y |\psi _{n'k}\rangle}{{({E_{n'}} - {E_n})}^2} $

      为倒空间的贝利曲率; ${\nu _x}$${\nu _y}$是电子沿x, y方向的速度算符. 比较有意思的是, 单层BaPb的反常霍尔电导曲线中可以观察到一个量子化的平台, 这是陈绝缘体的特征. 在图5(a)中的结果表明, 系统的Chern数为1, 这意味着在导带和价带之间只有一条边缘态连接. 同时, 量子反常霍尔态沿高对称路径分布的贝利曲率如图5(b)所示. 可以看到在倒空间中贝利曲率在Γ点附近表现出高对称性, 从另一个方面证明了系统可以实现QAHE.

      此外, Chern绝缘体最显著的特征是在纳米带边界处存在手性边缘态. 换句话说, 非平庸量子反常霍尔态将确保量子化手性边缘态的存在. 因此, 本文基于最局域Wannier函数构造了BaPb半无限纳米片的格林函数, 并计算了边缘处的局域态密度. 体系Chern数的绝对值等于其边缘态数目. 可以看到, 在费米能级附近的导带和价带之间存在一个拓扑保护的手性边缘态, 这与Chern数绝对值为1时的结果相一致, 如图5(c)所示. 此前计算出的在Γ点打开的局域间隙为177.39 meV, 这意味着单层BaPb中实现的QAHE有望超过室温, 大大高于最近在掺杂Cr的Bi2Se3薄膜中观察到QAHE的30 mK的温度.

    • 考虑到二维材料的电子能带结构和磁性质的稳定性对于量子器件的实验制备和发展非常重要[37,38]. 本部分以六方氮化硼(hexagonal boron nitride, h-BN)为衬底研究了单层BaPb的能带结构. 在h-BN衬底上生长的单层BaPb可能的晶格结构的俯视图和侧视图如图6(a)图6(b)所示. 定义晶格失配率为$\left| {{a_{{\rm{h}} \text- {\rm{BN}}}} - {a_{{\rm{BaPb}}}}} \right|/{a_{\text{h-}{\rm{BN}}}}$, 其中${a_{{\rm{h}} \text- {\rm{BN}}}}$${a_{{\rm{BaPb}}}}$分别为h-BN衬底和单层BaPb的晶格常数. 在本文的模型中, 使用了1 × 1的BaPb和$\sqrt 7 \times \sqrt 7 $的h-BN衬底, 计算得到两者的晶格失配率非常小, 为5.71%.

      图  6  (a), (b) h-BN/BaPb异质结(h-BN/BaPb)体系中单分子层BaPb的结构晶体结构的俯视图和侧视图; (c) h-BN/BaPb异质结自旋极化能带结构; (d) 考虑SOC效应时h-BN/BaPb异质结体系中单分子层BaPb的带结构(考虑SOC效应)

      Figure 6.  Crystal structure of monolayer BaPb in the h-BN/BaPb heterostructure system: (a) Top and (b) side views. Here, the black solid lines denote the unit cell. (c) Spin-polarized band structure of h-BN/BaPb heterostructure; (d) the band structures (including SOC) of monolayer BaPb in the h-BN/BaPb heterostructure system.

      图6(c)所示, 在h-BN/BaPb异质结构中, 费米能级附近的二次型非狄拉克带色散仍然存在. 考虑SOC效应时, 二次型非狄拉克能带打开了全局带隙82.11 meV, 如图6(d)所示. 这些结果表明, h-BN衬底对单分子层BaPb的能带结构几乎没有影响. 因此, 在狄拉克锥上观察到的SOC诱导的非平庸体带隙可能表明在该异质结构中可以实现量子反常霍尔绝缘态, 用于潜在的器件应用.

      最后, 还用蒙特卡罗模拟了h-BN/BaPb异质结中磁矩与温度的关系, 如图7所示. 结果表明, 由于基底导致的近邻效应可能影响异质结中电子跃迁和磁交换机制, 所以交换作用参数J值变为2.39 meV, 相应的Tc变为54 K. 但由于非磁性的基底h-BN与BaPb是通过范德瓦耳斯作用结合的, 因此体系的磁各向异性能几乎不变, 为51.32 meV.

      图  7  利用蒙特卡罗模拟得到的h-BN/BaPb异质结的磁矩与温度的关系

      Figure 7.  Average magnetic moment per unit cell with respect to temperature calculated for the h-BN/BaPb heterojunction obtained by using Monte Carlo simulation.

    • 本文利用第一性原理计算研究了单层BaPb作为新型Chern绝缘体候选的电子、磁性和拓扑性质. 声子谱计算验证了单层BaPb的动力学稳定性. 这种二维材料的负的形成能进一步证实了实验实现的可能性. 在估算居里温度为378 K左右的情况下, 它具有铁磁基态, 可在垂直于平面的坐标轴上得到巨大的MAE, 磁矩沿z方向时, 体系呈现出拓扑绝缘态, 而沿x方向则显示出金属态. 这意味着通过施加电场或构造异质结改变磁化方向, 可实现体系从金属态到拓扑绝缘态的相变. 最有趣的是, 在一个非零的Chern数(C = 1)下, 二次型非狄拉克点通过SOC效应打开了177.39 meV的非平庸拓扑间隙. 这种既具有全自旋极化的能带结构, 又有超大带隙, 而且还能实现金属态到拓扑绝缘态的相变结果对于将来纳米自旋电子器件的应用是非常有利的.

参考文献 (38)

目录

    /

    返回文章
    返回