应变与关联效应调控单层CeBr2谷极化行为

上一篇

下一篇

王一丁, 陈相, 李勇, 佘彦超, 张蔚曦. 应变与关联效应调控单层CeBr2谷极化行为[J]. 物理学报, 2026, 75(15): 150706. doi: 10.7498/aps.75.20260276
引用本文: 王一丁, 陈相, 李勇, 佘彦超, 张蔚曦. 应变与关联效应调控单层CeBr2谷极化行为[J]. 物理学报, 2026, 75(15): 150706. doi: 10.7498/aps.75.20260276
Yiding WANG, Xiang CHEN, Yong LI, Yanchao SHE, Weixi ZHANG. Strain and correlation effects modulate the valley polarization behavior in monolayer CeBr2[J]. Acta Physica Sinica, 2026, 75(15): 150706. doi: 10.7498/aps.75.20260276
Citation: Yiding WANG, Xiang CHEN, Yong LI, Yanchao SHE, Weixi ZHANG. Strain and correlation effects modulate the valley polarization behavior in monolayer CeBr2[J]. Acta Physica Sinica, 2026, 75(15): 150706. doi: 10.7498/aps.75.20260276

应变与关联效应调控单层CeBr2谷极化行为

    作者简介: 王一丁: 15281928253@163.com .
    通讯作者: E-mail: sheyanchao316@163.com.;  E-mail: zhangwwxx@sina.com.

Strain and correlation effects modulate the valley polarization behavior in monolayer CeBr2

    Corresponding authors: Corresponding author E-mail: sheyanchao316@163.com.;  Corresponding author E-mail: Zhangwwxx@sina.com.
  • 摘要: 谷自由度作为电子内禀自由度之一, 在低功耗信息器件中具有重要应用前景. 实现稳定的谷极化是谷电子学发展的关键, 其中铁谷材料因其本征磁有序可自发产生谷极化而备受关注. 本文基于密度泛函理论, 系统研究了二维CeBr2单层的结构稳定性、磁学性质、电子结构及其谷极化行为. 计算表明, 单层CeBr2具有良好的热力学和动力学稳定性, 其磁基态为铁磁序, 居里温度高于室温, 且表现出明显的面内磁各向异性($ \text{MAE}\approx 111\;\text{meV} $). 自旋-轨道耦合作用引起KK' 能谷退简并, 在价带顶产生约38.7 meV的自发能谷劈裂, 显著高于室温热涨落能量, 具备室温下稳定的谷极化能力. 此外, 计算揭示了该材料在能谷附近具有非零的反常霍尔电导与自旋霍尔电导, 支持反常谷霍尔效应的实现. 通过调控Hubbard U参数与施加双轴应变, 研究发现电子关联强度有利于增强磁各向异性与谷极化; 而在±5%的应变范围内, 谷极化与磁各向异性均可被连续调控, 且铁磁基态保持稳定. 本工作不仅为二维铁谷材料体系提供了一个具有大谷极化、高稳定性的新成员, 也为多场调控的谷自旋电子器件设计提供了理论依据.
  • 加载中
  • 图 1  (a)单层CeBr2晶体结构的俯视图和侧视图; (b)单层CeBr2的电子局域函数(ELF); (c)沿着第一布里渊区高对称路径计算的单层CeBr2的声子色散曲线; (d)具有倒易晶格矢量$ {\boldsymbol{b}}_{1} $$ {\boldsymbol{b}}_{2} $的蜂窝状晶格的第一布里渊区; (e)单层CeBr2的Ce原子轨道能级分布示意图

    Figure 1.  (a) Top and side views of the crystal structure of monolayer CeBr2; (b) electron localization function (ELF) of monolayer CeBr2; (c) phonon dispersion curves of monolayer CeBr2 calculated along the high-symmetry paths in the first Brillouin zone; (d) first Brillouin zone of the honeycomb lattice with reciprocal lattice vectors $ {\boldsymbol{b}}_{1} $ and $ {\boldsymbol{b}}_{2} $; (e) schematic illustration of the orbital energy-level distribution of Ce atoms in monolayer CeBr2.

    图 2  (a)单层CeBr2的FM, AFM自旋构型; (b)磁化方向位于xz平面上的单层CeBr2的MAE的角度依赖性; (c)单层CeBr2 的MAE不同f轨道贡献3D 柱状图

    Figure 2.  (a) FM and AFM spin configurations of monolayer CeBr2; (b) angular dependence of the MAE of monolayer CeBr2 with the magnetization direction confined to the xz plane; (c) three-dimensional bar chart showing the f-orbital-resolved contributions to the MAE of monolayer CeBr2.

    图 3  (a)单层CeBr2在无自旋-轨道耦合(SOC)下的能带结构, 绿色实线与蓝色虚线分别代表自旋向上和自旋向下的能带; (b)有SOC效应的单层CeBr2的能带图像; (c), (d)单层CeBr2的反常霍尔电导率和自旋霍尔电导率图像

    Figure 3.  (a) Band structure of monolayer CeBr2 without spin-orbit coupling (SOC), where the green solid lines and blue dashed lines represent the spin-up and spin-down energy bands, respectively; (b) projected band structure of monolayer CeBr2 with SOC; (c) anomalous Hall conductivity of monolayer CeBr2; (d) spin Hall conductivity of monolayer CeBr2.

    图 4  (a)单层CeBr2 AFM与FM磁态的能量差$ \Delta {E}_{{{}_{\text{AFM-FM}}}} $随Hubbard U的变化; (b)单层CeBr2的MAE随Hubbard U的变化; (c)单层CeBr2价带谷劈裂随Hubbard U的变化; (d)单层CeBr2带隙随Hubbard U的变化

    Figure 4.  (a) Energy difference between AFM and FM states of monolayer CeBr2 as a function of Hubbard U; (b) MAE of monolayer CeBr2 versus Hubbard U; (c) valley splitting at the valence band maximum and conduction band minimum of monolayer CeBr2 as a function of Hubbard U; (d) band gap of monolayer CeBr2 versus Hubbard U.

    图 5  (a)单层CeBr2 AFM与FM磁态的能量差$ \Delta {E}_{{{}_{\text{AFM-FM}}}} $随双轴应变的变化; (b)单层CeBr2的MAE随双轴应变的变化; (c)单层CeBr2价带谷劈裂随双轴应变的变化; (d)单层CeBr2的Ce原子轨道磁矩随双轴应变的变化

    Figure 5.  (a) Energy difference between AFM and FM states of monolayer CeBr2 as a function of strain; (b) MAE of monolayer CeBr2 as a function of strain; (c) valley splitting of monolayer CeBr2 as a function of strain; (d) orbital magnetic moment of the Ce atom of monolayer CeBr2 as a function of strain.

  • [1] Schaibley J R, Yu H Y, Clark G, Rivera P, Ross J S, Seyler K L, Wang Y, Xu X D 2016 Nat. Rev. Mater. 1 16055 doi: 10.1038/natrevmats.2016.55
    [2] Xiao D, Chang M C, Niu Q 2010 Rev. Mod. Phys. 82 1959 doi: 10.1103/RevModPhys.82.1959
    [3] Tong W Y, Gong S J, Wan X G, Duan C G 2016 Nat. Commun. 7 13612 doi: 10.1038/ncomms13612
    [4] Huang B, Clark G, Navarro Moratalla E, et al. 2017 Nature 546 270 doi: 10.1038/nature22391
    [5] Gong C, Li Z L, Li Z, Ji H W, Stern A, Xia Y, Cao T, Bao W, Wang C Z, Wang Y, Qiu Z Q, Cava R J, Louie S G, Xia J, Zhang X 2017 Nature 546 265 doi: 10.1038/nature22060
    [6] Sheng K, Chen Q, Yuan H K, Wang Z Y 2022 Phys. Rev. B 105 075304 doi: 10.1103/PhysRevB.105.075304
    [7] Sun H B, Ren Y W, Wu C, Dong P Q, Zhang W X, Wu Y Z, Li P 2025 Phys. Rev. Appl. 23 034032 doi: 10.1103/PhysRevApplied.23.034032
    [8] Wang Y D, Sun H B, Wu C, Zhang W X, Guo S D, She Y C, Li P 2025 Phys. Rev. B 111 085432 doi: 10.1103/PhysRevB.111.085432
    [9] She Y C, Wang Y D, Sun H B, Wu C, Zhang W X, Li P 2026 Phys. Rev. B 113 035420 doi: 10.1103/8c5h-31j6
    [10] 孙真昊, 管鸿明, 付雷, 沈波, 汤宁 2021 物理学报 70 027302 doi: 10.7498/aps.70.20201415 Sun Z H, Guan H M, Fu L, Shen B, Tang N 2021 Acta Phys. Sin. 70 027302 doi: 10.7498/aps.70.20201415
    [11] Mermin N D, Wagner H 1966 Phys. Rev. Lett. 17 1133 doi: 10.1103/PhysRevLett.17.1133
    [12] Jungwirth T, Sinova J, Manchon A, Marti X, Wunderlich J, Felser C 2018 Nat. Phys. 14 200 doi: 10.1038/s41567-018-0063-6
    [13] Guo S D, Zhang L G, Zhang Y W, Li P, Wang G Z 2024 Phys. Rev. B 110 024416 doi: 10.1103/PhysRevB.110.024416
    [14] Hedjar H, Meskine S, Boukortt A, Bennacer H, Benzidane M R 2022 Comput. Condens. Matter. 30 e00632 doi: 10.1016/j.cocom.2021.e00632
    [15] Cheng H X, Zhou J, Ji W, Zhang Y N, Feng Y P 2021 Phys. Rev. B 103 125121 doi: 10.1103/PhysRevB.103.125121
    [16] Sheng K, Yuan H K, Wang Z Y 2022 Phys. Chem. Chem. Phys. 24 3865 doi: 10.1039/D1CP05097D
    [17] Wu H F, Yan R F, Chen Y K, An Y K 2024 J. Magn. Magn. Mater. 611 172600 doi: 10.1016/j.jmmm.2024.172600
    [18] Chen Y K, Zhao X S, An Y K 2024 Phys. Rev. B 109 125421 doi: 10.1103/PhysRevB.109.125421
    [19] Anderson P W 1950 Phys. Rev. 79 350 doi: 10.1103/PhysRev.79.350
    [20] Dong L, Lou J, Shenoy V B 2017 ACS nano 11 8242 doi: 10.1021/acsnano.7b03313
    [21] Lang H F, Zhang S Q, Liu Z R 2016 Phys. Rev. B 94 235306 doi: 10.1103/PhysRevB.94.235306
    [22] Zhao P, Ma Y D, Lei C A, Wang H, Huang B B, Dai Y 2019 Appl. Phys. Lett. 115 261605 doi: 10.1063/1.5129311
    [23] 谢伟锋, 王利波, 许雄, 岳云亮, 夏华艳, 何龙辉, 王辉 2025 物理学报 74 227502 doi: 10.7498/aps.74.20251066 Xie W F, Wang L B, Xu X, Yue Y L, Xia H Y, He L H, Wang H 2025 Acta Phys. Sin. 74 227502 doi: 10.7498/aps.74.20251066
    [24] 朱靖欣, 郭三栋 2025 原子与分子物理学报 42 026005 doi: 10.19855/j.1000-0364.2025.026005 Zhu J X, Guo S D 2025 J. At. Mol. Phys. 42 026005 doi: 10.19855/j.1000-0364.2025.026005
    [25] Blöchl P E 1994 Phys. Rev. B 50 17953 doi: 10.1103/PhysRevB.50.17953
    [26] Kresse G, Furthmuller J 1996 Phys. Rev. B 54 11169 doi: 10.1103/PhysRevB.54.11169
    [27] Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865 doi: 10.1103/PhysRevLett.77.3865
    [28] Li S J, Hou Y F, Zhou M, Zheng F W, Shao X H, Zhang P 2023 J. Appl. Phys. 134 213901 doi: 10.1063/5.0178739
    [29] Togo A, Tanaka I 2015 Scr. Mater. 108 1 doi: 10.1016/j.scriptamat.2015.07.021
    [30] Marzari N, Mostofi A A, Yates J R, Souza I, Vanderbilt D 2012 Rev. Mod. Phys. 84 1419 doi: 10.1103/RevModPhys.84.1419
    [31] Li P, Wu C, Peng C, Yang M T, Xun W 2023 Phys. Rev. B 108 195424 doi: 10.1103/PhysRevB.108.195424
    [32] Thouless D J, Kohmoto M, Nightingale M P, den Nijs M 1982 Phys. Rev. Lett. 49 405 doi: 10.1103/PhysRevLett.49.405
  • 加载中
图( 6)
计量
  • 文章访问数:  109
  • HTML全文浏览数:  109
  • PDF下载数:  3
  • 施引文献:  0
出版历程
  • 收稿日期:  2026-02-25
  • 刊出日期:  2026-08-05

应变与关联效应调控单层CeBr2谷极化行为

    通讯作者: E-mail: sheyanchao316@163.com.; 
    通讯作者: E-mail: zhangwwxx@sina.com.
    作者简介: 王一丁: 15281928253@163.com
  • 1. 铜仁学院物理与电子工程系, 铜仁 554300
  • 2. 吉首大学物理与机电工程学院, 吉首 416000

摘要: 谷自由度作为电子内禀自由度之一, 在低功耗信息器件中具有重要应用前景. 实现稳定的谷极化是谷电子学发展的关键, 其中铁谷材料因其本征磁有序可自发产生谷极化而备受关注. 本文基于密度泛函理论, 系统研究了二维CeBr2单层的结构稳定性、磁学性质、电子结构及其谷极化行为. 计算表明, 单层CeBr2具有良好的热力学和动力学稳定性, 其磁基态为铁磁序, 居里温度高于室温, 且表现出明显的面内磁各向异性($ \text{MAE}\approx 111\;\text{meV} $). 自旋-轨道耦合作用引起KK' 能谷退简并, 在价带顶产生约38.7 meV的自发能谷劈裂, 显著高于室温热涨落能量, 具备室温下稳定的谷极化能力. 此外, 计算揭示了该材料在能谷附近具有非零的反常霍尔电导与自旋霍尔电导, 支持反常谷霍尔效应的实现. 通过调控Hubbard U参数与施加双轴应变, 研究发现电子关联强度有利于增强磁各向异性与谷极化; 而在±5%的应变范围内, 谷极化与磁各向异性均可被连续调控, 且铁磁基态保持稳定. 本工作不仅为二维铁谷材料体系提供了一个具有大谷极化、高稳定性的新成员, 也为多场调控的谷自旋电子器件设计提供了理论依据.

English Abstract

    • 谷自由度作为电子电荷与自旋之外的重要内禀自由度, 因其在动量空间具有高鲁棒性和可选择调控的特点, 已成为发展超快、低功耗非易失性信息器件的关键研究方向之一, 推动了谷电子学在凝聚态物理和材料科学领域的研究热潮[13]. 实现稳定的谷极化是谷电子学走向应用的核心前提, 通常依赖于对称性破缺的调控手段. 其中, 铁谷材料凭借其本征磁性打破时间反演对称性, 无需外场即可产生自发谷极化, 被认为是实现谷电子学与自旋电子学融合的理想载体[46].

      二维磁性材料以其原子级厚度、高柔韧性和易集成等优势, 为铁谷体系的构筑提供了理想平台[710]. 然而, 传统二维铁磁材料普遍面临居里温度偏低、磁各向异性弱以及谷极化调控方式有限等问题; 反铁磁材料虽具有净磁矩为零、抗干扰能力强等优点, 却因保持空间反演与时间反演联合对称性而难以实现自发谷极化, 制约了其实际应用[1113].

      近年来, 稀土基卤化物因局域f电子所带来的强磁各向异性和显著谷极化, 成为铁谷材料研究的重要体系[1416]. 其中Ce基卤化物表现尤为突出: 例如CeI2单层为室温铁谷半导体, 谷极化可达208 meV, 居里温度达374 K[6]; Janus结构的2H-CeXY单层不仅居里温度高达511—540 K, 同时具有29.1—78.7 meV的谷极化, 且其磁各向异性可通过应变进行调控[17,18]. 这些材料的铁磁性主要源于Ce原子与卤素原子之间的$ \text{d-p-d} $间接交换作用, 谷极化则由自旋-轨道耦合与磁有序共同诱导, 并可通过外场进一步调控[1921].

      尽管CeBr2与上述已报道的铁谷材料具有相似的层状六角晶格, 同为稀土基卤化物家族的重要成员, 但其结构稳定性、磁基态、谷极化行为以及外场调控规律等方面仍缺乏系统性研究, 相关物理机制尚不明确[22]. 与交错磁体中依赖元素替换或异质结实现物理性质调控的策略[23]、Janus LuClI中依赖压电效应实现谷极化调控的策略[24]不同, CeBr2作为本征稀土基二维材料, 其强磁各向异性、室温铁磁序及谷极化行为有望通过应变、电子关联等多种外场直接调控, 无需复杂材料工程. 这一特性使其成为谷电子学材料体系扩充的理想候选. 为此, 本研究基于密度泛函理论计算, 系统研究CeBr2的结构稳定性、磁学性质、电子结构及谷极化特性, 通过声子谱确认其稳定性, 揭示其磁基态、居里温度、磁各向异性及轨道贡献机制, 阐明自旋-轨道耦合诱导的谷极化与能谷劈裂行为, 并进一步探究Hubbard U、应变等条件对磁性与谷极化的调控规律, 验证该体系中可能存在的反常谷霍尔效应与自旋霍尔效应. 本研究旨在扩充二维铁谷材料体系, 深入揭示CeBr2中谷、自旋与轨道之间的耦合机制, 为未来多场调控的低功耗谷自旋电子器件提供理论依据与材料选择.

    • 本文基于密度泛函理论(density functional theory, DFT)开展第一性原理计算, 所有计算均采用维也纳从头算模拟软件包VASP(Vienna ab-initio simulation package)[25,26]完成, 并使用投影增强波(projector augmented-wave, PAW)赝势描述离子实与价电子之间的相互作用. 电子交换关联作用采用广义梯度近似(generalized gradient approximation, GGA)下的 Perdew-Burke-Ernzerhof(PBE)泛函进行处理[27]. 平面波基组的截断能设置为600 eV. 结构优化过程中, 体系的几何构型通过完全弛豫获得, 其收敛判据设定为每个原子所受残余力小于0.005 eV/Å, 总能量变化小于10–6 eV. 为消除周期性边界条件在面外方向引入的人为层间相互作用, 在垂直于二维材料平面的c轴方向设置了厚度为20 Å的真空层. 考虑到Ce原子4f电子存在显著的强关联效应, 计算中引入GGA+U方法进行修正, 其中有效哈伯德参数$ {U}_{\text{eff}}=U-J $取值为7.0 eV[28]. 布里渊区积分采用 Monkhorst-Pack方法生成的k点网格, 并通过收敛性测试以确保计算结果的可靠性. 为进一步验证所研究结构的动力学稳定性, 基于有限位移法并借助PHONOPY软件包[29], 在$ 3\times 3\times 1 $超胞模型下计算了体系的声子色散关系. 霍尔电导率均通过WANNIER90软件包中编码的最大局域化Wannier函数计算[30].

    • 单层CeBr2的晶体结构如图1(a)所示, 其空间群为$ P\overline{6}m2 $, 属于典型的六方晶系结构. 该体系呈现出“Br-Ce-Br”三明治式层状构型, 其中Ce原子层位于中间, 并在面内形成蜂窝状排列; 上下两层Br原子分别位于Ce原子层的两侧, 与Ce原子共同构成三角棱柱型配位环境. 这种层状配位结构为体系中磁交换相互作用及能谷相关物理行为的产生提供了必要的几何基础. 经过结构优化后, 单层CeBr2的平衡晶格常数为4.08 Å, Ce—Br键长为3.03 Å, 对应的Ce—Br—Ce键角约为84.6°. 上述结构参数与稀土卤化物体系中常见的局域配位特征一致, 表明该结构在几何构型上具有良好的合理性与稳定性. 为进一步分析体系的成键特性, 计算了其电子局域函数(electron localization function, ELF), 结果如图1(b)所示. 可以看到, 电子主要局域分布在Br原子周围, 而Ce与Br原子之间的ELF值相对较低, 表明电子从Ce原子向Br原子发生明显转移, Ce—Br键呈现出以离子键为主的成键特征. 这一结果与稀土基卤化物中常见的电荷转移行为相一致, 也反映了Ce原子4f电子的高度局域化特性. 为验证单层CeBr2的动力学稳定性, 进一步计算了其声子色散关系, 如图1(c)所示. 在整个布里渊区范围内未观察到虚频模式, 表明该体系在动力学上是稳定的. 声子谱中高频振动模式主要来源于质量较轻的Br原子, 而低频区域则对应于Ce原子的面内振动模式, 这与体系中不同原子质量及其配位环境密切相关.

      此外, 图1(e)给出了CeBr2体系中Ce原子的轨道能级分布情况. 可以看出, 在晶体场作用及强电子关联效应的共同影响下, Ce2+离子的4f与5d轨道发生了明显的能级分裂. 其中, 4f轨道(如$ {\text{f}}_{{{x}^{2}}-3{{y}^{2}}} $, $ {\text{f}}_{3{{x}^{2}}y-{{y}^{3}}} $, $ {\text{f}}_{xyz} $, $ {\text{f}}_{z{{x}^{2}}-z{{y}^{2}}} $等)表现出复杂的多重分裂特征, 反映了其在低对称配位环境中的强关联本征性质; 5d轨道(包括$ {\text{d}}_{xy} $, $ {\text{d}}_{xz} $, $ {\text{d}}_{yz} $, $ {\text{d}}_{{{x}^{2}}-{{y}^{2}}} $轨道对以及$ {\text{d}}_{{{z}^{2}}} $单轨道)同样在晶体场作用下产生能级差异. 图中箭头标注了各轨道电子占据数的变化情况, 整体揭示了Ce2+不同轨道在局域配位场中的分裂规律及其电子填充特征, 为后续磁性与谷相关物理性质的分析提供了微观能级基础.

    • 单层CeBr2的磁性主要起源于Ce原子的电子组态及其4f电子的强关联特性. 孤立Ce原子的价电子组态为$ {\text{4f}}^{1}{\text{5d}}^{1}{\text{6s}}^{2} $. 在形成单层CeBr2的过程中, 每个Ce原子向周围6个Br原子转移2个电子, 从而形成Ce2+离子, 其对应的电子组态为$ {\text{4f}}^{1}{\text{5d}}^{1}{\text{6s}}^{0} $. 在三角棱柱型晶体场作用下, Ce原子的4f与5d轨道发生能级分裂, 其中半填充的$ 4{\text{f}}_{{{x}^{2}}-3{{y}^{2}}} $$ 5{\text{d}}_{{{z}^{2}}} $轨道共同贡献了体系的局域磁矩, 使得每个Ce原子产生约$ 2{\mu }_{\text{B}} $的理论磁矩. 计算结果进一步显示Ce原子的局域磁矩仅为$ 1.47{\mu }_{\text{B}} $, 即Ce-4f壳层的磁矩约为$ 1{\mu }_{\text{B}} $, Ce-5d壳层的磁矩约为$ 0.4{\mu }_{\text{B}} $. 相比之下, Br原子的磁矩贡献可以忽略不计. 这表明体系的净磁矩主要来源于Ce原子的f轨道电子, 其通过间接交换机制在晶格中建立长程磁有序.

      为确定单层CeBr2的磁基态, 构建了$ 2\times 2\times 1 $超胞, 并分别考虑铁磁(ferromagnetic, FM)和反铁磁(antiferromagnetic, AFM)两种磁构型, 如图2(a) 所示. 参考已有稀土基卤化物体系(如CeI2)中关于强关联参数的取值经验, 采用 GGA+U方法($ {U}_{\text{eff}}=7.0\, \, \, \, \text{eV} $)计算了两种磁态的总能. 结果表明, FM态相较于AFM态能量低约285.54 meV, 明确指示单层CeBr2的磁基态为铁磁序.

      对于二维磁性材料而言, 易磁化方向在决定其磁矩稳定取向及潜在自旋电子学应用中起着关键作用. 为此, 进一步计算了体系的磁各向异性能(magnetic anisotropy energy, MAE), 其定义为$ E_\text{MAE}={E}_{001}-{E}_{100} $, 其中$ {{E}}_{100} $$ {{E}}_{001} $分别表示磁矩沿面外[001]方向和面内[100]方向取向时的总能. 当EMAE > 0时, 表明易磁化轴位于面内; 反之, EMAE < 0则对应易磁化轴在面外. 计算得到的EMAE为111 meV, 表明体系的易磁化方向位于面内, 即[100]方向. 结合Ce原子所处的三角棱柱晶体场对称性, 随磁矩取向角的变化磁各向异性能$ E_{\text{MAE}\left(\theta \right)}$可用(1)式进行拟合[31]:

      其中$ \theta $为磁矩方向与面外z轴之间的夹角, $ {K}_{1} $$ {K}_{2} $为磁各向异性常数. 如图2(b) 所示, 曲线最小值对应 θ = 90°即磁矩沿面内方向, 最大值对应 θ = 0°即磁矩沿面外方向. MAE在xz平面内呈现明显的各向异性特征, 其数值对磁矩在平面内的具体取向表现出显著依赖性.

      由于f轨道的自旋-轨道耦合(spin-orbit coupling, SOC)远远强于d轨道, 所以针对f轨道, 探究了其在MAE调控中的核心作用. 我们聚焦f轨道间的相互作用矩阵, 提取主导MAE的非对角矩阵元并构建轨道分辨三维柱状图(图2(c)). 正方向表示倾向于面内(x)贡献, 负方向表示倾向于面外(z)贡献. 从柱状图分布可见, 不同f轨道组合的MAE贡献呈现明显分化: 高幅值正贡献柱体集中分布于特定轨道组合(如$ {\text{f}}_{z{{x}^{2}}}{\text{-f}}_{{{z}^{3}}} $, $ {\text{f}}_{x{{z}^{2}}}{\text{-f}}_{y{{z}^{2}}} $), 对应非对角矩阵元绝对值最大, 是决定体系MAE大小的主导通道; 而部分轨道组合呈现负贡献, 柱体幅值显著低于正贡献通道. 这一分布特征直接反映出Ce原子f轨道间自旋-轨道耦合的非均匀性, f轨道相较于d轨道更复杂的能级分裂与轨道重叠模式, 进一步塑造了各轨道组合贡献的差异化, 最终共同决定了体系的磁各向异性能特征.

      为进一步定量评估单层CeBr2铁磁基态的热稳定性, 基于海森伯模型与平均场理论估算了其居里温度 $ {T}_{\text{c}} $. 考虑到体系的磁基态为铁磁序, 且最近邻Ce原子间的磁交换作用占主导, 可采用仅包含最近邻交换相互作用的哈密顿量:

      其中 $ J $为最近邻磁交换耦合常数, $ {S}_{i} $$ {S}_{j} $分别为格点$ i $和格点$ j $上Ce原子的自旋算符. 由 DFT 计算得到的铁磁与反铁磁构型能量差$ \Delta E={E}_{\text{AFM}}- {E}_{\text{FM}}=284.54\;\text{meV} $可推导出$ J $的表达式. 对于六角蜂窝晶格, 铁磁与反铁磁态的能量差满足:

      其中$ z=3 $为最近邻配位数, $ S $为Ce离子的自旋量子数. 由Ce原子的局域磁矩$ m=1.47{\mu }_{\text{B}} $, 结合 $ m=gS{\mu }_{\text{B}} $(取朗德因子$ g\approx 2 $), 可得$ S\approx 0.735 $. 代入(3)式解得

      在平均场近似下, 二维六角晶格的居里温度可表示为

      计算得到 $ T_{\text{C}}^{\text{MF}}\approx 390\;\text{K} $, 该值显著高于室温, 与实验可实现的室温应用需求相吻合.

    • 为深入揭示单层CeBr2的铁谷特性及其输运行为, 系统研究了其电子能带结构、自发谷劈裂以及反常谷霍尔效应. 图3(a)展示了未考虑SOC时的自旋分辨能带结构. 单层CeBr2表现为双极磁性半导体, 具有间接带隙(~0.44 eV), 价带顶(valence band maximum, VBM)位于KK' 点且保持自旋简并, 而导带底(conduction band minimum, CBM)位于M点. 费米能级附近的自旋向上与自旋向下能带呈现显著劈裂, 自旋极化率接近100%, 这源于Ce原子局域磁矩与电子态的强耦合, 为谷与自旋自由度的协同调控奠定了基础.

      引入SOC后(图3(b)), 时间反演对称性与空间反演对称性被同时打破, 导致K$K' $点的能谷发生退简并. 价带顶处出现约38.7 meV的能谷劈裂, 该值超室温热扰动能(~26 meV), 表明CeBr2具备室温下稳定的自发能谷极化能力. 轨道投影分析表明, 费米能级附近的能带主要由Ce-5d轨道($ {\text{d}}_{xy} $, $ {\text{d}}_{xz} $, $ {\text{d}}_{yz} $, $ {\text{d}}_{{{x}^{2}}-{{y}^{2}}} $$ {\text{d}}_{{{z}^{2}}} $)贡献, Br-p轨道贡献微弱, 体现了Ce-5d轨道在低能电子结构中的主导作用.

      由于单层CeBr2的铁磁性打破了自旋向上能带和自旋向下能带之间的自旋简并, SOC哈密顿量仅涉及相同自旋态的相互作用. 为描述能谷劈裂的物理起源, 构建了有效SOC哈密顿量模型[6]:

      其中λ为自旋-轨道耦合相关常数, 决定耦合强度; $ {\hat{S}}_{Z'} $$ {\hat{L}}_{Z} $分别是自旋和轨道角动量的$ {Z}^{\prime}/Z $分量; θ∈[0, π)和Ф∈[0, 2π)两者定义自旋取向, 其分别对应于极角和方位角; $ {\hat{L}}_{+}={\hat{L}}_{x}+\text{i}{\hat{L}}_{y} $, $ {\hat{L}}_{-}={\hat{L}}_{x}- \text{i}{\hat{L}}_{y} $. 考虑到Ce原子d轨道在费米能级附近的主导作用, 且$ {\text{d}}_{{{z}^{2}}} $轨道的轨道角动量为零、对能谷劈裂贡献可忽略, 模型聚焦于$ {\text{d}}_{xy} $, $ {\text{d}}_{{{x}^{2}}-{{y}^{2}}} $两类轨道的耦合效应, SOC为微扰项. 因此, 价带能谷劈裂由下式给出[6]:

      该劈裂源于SOC诱导下KK' 点轨道波函数的不对称耦合: 在K点, $ {\text{d}}_{xy} $$ {\text{d}}_{{{x}^{2}}-{{y}^{2}}} $轨道的杂化强度显著高于K' 点, 而$ {\text{d}}_{xz} $$ {\text{d}}_{yz} $轨道的耦合则呈现相反趋势, 这种轨道杂化的空间不对称性最终导致能谷能级分裂, 与图1(e)中Ce原子轨道能级分裂规律形成自洽关联.

      为探究CeBr2在谷电子学器件中的应用潜力, 本文通过Wannier函数进一步计算了其AHC与SHC. 对于二维铁磁系统, 第n个能带的贝利曲率可以写为[32]

      从Kubo公式导出的总Berry曲率是所有占据态的总和[32]:

      其中$ {f}_{nk} $为费米-狄拉克分布, $ {v}_{x/y} $是速度算符. 通过在2D布里渊区(BZ)上积分Berry曲率计算二维本征反常霍尔电导$ {\sigma }_{xy} $[32]:

      图3(c)为采用3.87 Å材料厚度3D化处理后的等效反常霍尔电导, AHC随能量变化呈现明显的谷依赖特征: 在费米能级附近–0.5 eV处出现约8 S/cm的极大值, 随后在0.5 eV附近降至约–13 S/cm的极小值, 之后逐渐回升. 该结果表明单层CeBr2在费米能级附近具有非零的反常霍尔电导, 预示着可观测的反常谷霍尔效应, 为谷极化电荷的输运调控提供了物理依据. 图3(d)为计算得到的自旋霍尔电导(spin Hall conductivity, SHC)随能量变化的曲线, 其能量轴与图3(b)的能带结构完全对应. 在能隙区间(费米能级附近, 能量约为0 eV处), 自旋霍尔电导呈现明显的平台状特征, 且在整个带隙范围内保持稳定、无明显波动. 这一结果表明单层 CeBr2在能隙区间具有显著的自旋霍尔响应, 为其在自旋输运器件中的应用提供了物理依据.

      综上, 单层CeBr2具有由SOC诱导的价带谷劈裂(38.7 meV), 并同时展现出明显的反常霍尔与自旋霍尔响应, 实现了谷、自旋、轨道与电荷自由度的多重耦合, 为其在低功耗、多场调控的谷自旋电子器件中的应用提供了理论基础.

      为系统考察电子关联作用对CeBr2物理性质的调控规律, 开展了Hubbard U参数(U = 4—8 eV)的依赖性计算. 结果表明, 电子关联强度对体系的磁学性质与电子结构均具有显著影响.

      在磁学性质方面, FM与反铁磁AFM态的能量差$ \Delta {E}_{{{}_{\text{AFM-FM}}}} $U值的变化呈现明显的双向极值特征, 如图4(a)所示: U = 5 eV 时能量差降至约–0.72 eV的极小值(表明此时AFM磁态相对FM的稳定性最强), 随后在U = 6 eV时跃升至约 0.27 eV的极大值(表明此时FM磁态相对AFM的稳定性最强), 且在U = 6—8 eV区间保持正值, 说明该关联强度范围内两种磁态的稳定性差异显著, 电子关联强度的变化可实现磁基态在AFM与FM之间的翻转. 磁各向异性能(MAE)则随着U的增大而增大, 如图4(b)所示, 反映出关联强度有利于增强磁各向异性; 当 U 增大至7—8 eV时, MAE的增长逐渐趋于饱和, 说明过强的电子关联对磁各向异性的进一步增强作用有限.

      在电子结构方面, 如图4(c), (d)所示, 价带的谷劈裂随U增大单调上升. 这一非对称调控行为主要源于Ce-4f与 Br-p轨道在不同关联强度下的杂化演化: 增强的电子关联使4f轨道更局域, 进而调制其与Br-p轨道的杂化强度, 从而影响不同能谷的轨道组成与能量顺序. 与此同时, 体系的带隙随U 增大而持续减小(见图4(d)). 这主要是因为, 在强关联电子体系中, Hubbard U参数的增大通常对应局域电子-电子库仑排斥作用的增强, 这将对体系的能带结构产生显著影响. 一方面, 关联修正会引起占据态(价带)与未占据态(导带)之间的能量重排, 表现为占据态能级下移、未占据态能级上移, 从而在静态平均场层面上倾向于扩大单粒子能隙. 另一方面, 对于以窄带f/d电子为主导的强关联体系, 增强的电子关联效应往往伴随着能带宽度的压缩以及有效质量的增大, 从而削弱能级分裂对能隙扩大的贡献. 在这种情况下, 价带顶与导带底在能量空间中的相对位置可能出现重新靠近, 导致计算得到的有效带隙随U的增大反而减小. 这种由关联效应引起的能隙重整化行为在多种稀土及过渡金属化合物中均已有报道. 在本文的计算中, 单CeBr2的带隙随Hubbard U参数的增大呈现持续减小的趋势, 表明在该体系中, 关联增强所引发的能带窄化效应占据主导地位. 该结果反映了强关联区域内价带与导带在能量上的相对接近, 预示着体系电子跃迁特性及低能电子结构将发生显著重构.

      综上所述, Hubbard U作为调控CeBr2电子关联强度的关键参数, 可有效调节其磁基态稳定性、磁各向异性、能谷劈裂及带隙大小, 当U = 7—8 eV时体系展现出最优的磁性与谷极化特性. 该结果为后续基于关联强度的材料设计与器件调控提供了清晰的理论依据.

    • 应变作为调控二维材料物理性质的有效外场手段, 在CeBr2中展现出对磁序、磁各向异性、谷极化及轨道磁矩的连续调控能力. 本文系统研究了双轴应变($ \varepsilon =(a-{a}_{0})/{a}_{0}\times 100{\text{%}}$, $ {a}_{0} $为平衡晶格常数)在$ \pm 5\% $范围内对单层CeBr2多物理性质的调控行为. 在未施加双轴应变时, FM态为体系的磁基态. 在ε = –5%—+5%的应变范围内(图5(a)), FM 态始终保持能量最低, 表明其磁基态在该应变区间内具有良好稳定性, 为应变可调自旋器件的磁态保持提供了基础.

      MAE随双轴应变呈现单调变化(图5(b)): 从ε = –5%时的169 meV减小至ε = +5%时的153 meV, 说明压缩应变可增强体系的磁各向异性. 这一变化源于应变引起的晶格畸变对Ce原子f轨道晶体场及自旋-轨道耦合强度的调制, 为通过应变实现磁化方向切换提供了可能.

      价带能谷劈裂值始终为正, 并且在无应变时达到最大值, 施加双轴应变(无论是压缩还是拉伸)均导致谷劈裂减小. 这一行为源于应变对Ce-5d轨道晶体场环境的扰动: 无应变时, Ce原子处于最优配位构型, 5d轨道在KK' 点的杂化不对称性最强, 自旋-轨道耦合诱导的谷劈裂最大; 施加应变后, 晶格畸变破坏这一最优构型, 削弱轨道杂化不对称性, 导致谷劈裂降低. 该结果表明, 双轴应变工程可在一定范围内连续调控谷极化强度, 但无应变状态对应最优谷劈裂性能.

      此外, Ce原子的轨道磁矩随双轴应变发生细微但系统的变化(图5(d)), 从$ -0.974{\mu }_{\text{B}} $(ε = –5%)逐渐变化至$ -0.982{\mu }_{\text{B}} $(ε = +5%). 这一变化反映了双轴应变对Ce局域电子结构的微扰, 其趋势与MAE的双轴应变响应相呼应, 共同揭示了双轴应变通过调制轨道自由度来调控磁性与谷相关的物理微观机制. 综上, 双轴应变可在一定范围内连续调节CeBr2的磁各向异性、谷极化及轨道磁矩, 且不改变其铁磁基态, 展现出该体系在多功能应变调控器件中的应用潜力.

    • 本文基于第一性原理计算, 系统研究了单层CeBr2的结构稳定性、磁学性质、电子结构及其谷极化与输运特性. 研究表明, 单层CeBr2为一种新型的本征二维铁谷材料, 具有良好的结构、热力学和动力学稳定性. 其铁磁基态由Ce-4f/5d与Br-p轨道的超交换作用所主导, 并表现出明显的面内磁各向异性(MAE = 111 meV). 在自旋-轨道耦合作用下, 体系在价带顶KK' 点处产生约38.7 meV的自发能谷劈裂, 该值超室温热涨落能量, 保障了谷极化在室温下的稳定性. 计算同时揭示了该材料在能谷附近具有非零的反常霍尔电导与自旋霍尔电导, 预示着其可观测的电荷与自旋输运效应, 支持反常谷霍尔效应的实现. 此外, 本研究通过调节Hubbard U参数与施加双轴应变, 系统考察了电子关联与外部应变对物性的调控规律. 结果表明, 适当的电子关联强度(U = 7—8 eV)能实现磁各向异性与谷极化的最优化; 而在±5%的双轴应变范围内, 谷极化与磁各向异性均可被连续调控, 且铁磁基态保持稳定, 展现出多物理场协同调控的潜力. 本工作不仅为二维铁谷材料家族提供了一个具有大谷极化、高稳定性的新成员, 而且从轨道耦合、电子关联及外场调控等多个维度, 深化了对稀土基卤化物中谷、自旋、轨道与电荷自由度之间耦合机制的理解. 这些结果为设计基于CeBr2的多场调控、低功耗谷自旋电子器件提供了可靠的理论依据与材料平台.

    参考文献 (32)

目录

/

返回文章
返回