反铁磁拓扑材料Cr2Se3的合成、物理性质与能带结构

反铁磁拓扑材料Cr2Se3的合成、物理性质与能带结构PHYS. REV. B 114, 024404 (2026)反铁磁拓扑材料Cr2Se3的合成、物理性质与能带结构Synthesis, Physical Properties, and Band Structure of Antiferromagnetic Cr2Se3导读 导读Cr2Se3是一种具有非中心对称三方结构的反铁磁拓扑材料。本文通过实验合成CVT单晶与第一性原理计算WIEN2KWannier90相结合系统研究了其结构、磁性、输运和拓扑性质。核心发现(1) Cr2Se3存在两个反铁磁相变温度TN1~38K和TN2~42K(2) 反常Hall效应得到实验和理论的双重支持(3) Z2(1;000)确认为强拓扑半金属拓扑性质由反演对称性保护。一、前言背景反铁磁拓扑材料从MnBi2Te4到Cr2Se3反铁磁拓扑材料由于自旋与拓扑序的独特耦合效应在低功耗自旋电子学器件和量子计算领域具有重要应用前景。MnBi2Te4体系的实验突破证实了AFM拓扑态结合的潜力但该体系存在机械稳定性差、化学计量比敏感等限制。Cr2Se3是几十年前就被确认的反铁磁材料但系统性的物理性质Hall电阻率、热容等和拓扑特征研究一直缺失。本文通过实验合成第一性原理计算填补了这一空白。核心发现(1) Cr2Se3存在两个反铁磁转变温度TN1~38K和TN2~42K(2) 可能存在的反常Hall效应得到理论计算支持(3) Z2(1;000)确认为强拓扑半金属。实验-理论联合研究从晶体生长到拓扑表征实验方面采用化学气相输运CVT法生长Cr2Se3单晶碘为输运剂。通过XRD、EDX、XPS、SQUID、PPMS进行全面的结构和物性表征。理论方面使用WIEN2K全势LAPW方法计算电子结构Wu-Cohen GGA泛函包含SOC。通过Wannier90构建最大局域化Wannier函数紧束缚模型WannierTools计算拓扑性质和反常Hall电导率。关键方法对比不同于多数VASP研究本文使用WIEN2K全势方法对于含重元素Se的体系全势方法在SOC描述上具有优势。Cr2Se3反铁磁拓扑材料研究流程。CVT单晶生长 - 结构/磁性/输运/热容全面表征 - WIEN2K DFT全势LAPW, Wu-Cohen GGASOC - Wannier90WannierTools拓扑分类、表面态、反常Hall电导率。核心发现两次AFM相变 Z2(1;000)强拓扑半金属 反常Hall效应。二、研究方法实验方法CVT单晶生长与多维度物性表征晶体生长化学气相输运CVT法Cr:Se2:3摩尔比碘为输运剂。双温区炉800/700度C2天升温10天恒温自然冷却。获得具有明显金属光泽的单晶。结构表征室温XRDRigaku, Cu Kalpha四圆单晶XRDBruker D8 Venture, 150K, Mo Kalpha空间群R-3No. 148ab6.2544(2) A, c17.2712(11) A。XPS确认Cr3和Se2-的化学态。物性测量SQUID磁强计MPMS3测量磁化率和磁化曲线PPMS测量电阻率、Hall效应和热容。DFT计算方法WIEN2K全势LAPW与Wannier函数WIEN2K设置全势线性化缀加平面波FP-LAPW方法Wu-Cohen GGA泛函非经验泛函在晶格常数和晶体结构描述上优于PBE包含SOC。AFM磁构型Cr1(Cr2)与Cr3自旋反平行。Wannier函数从WIEN2K能带构建Cr-d和Se-p轨道的最大局域化Wannier函数MLWFs生成紧束缚模型。WannierTools计算Z2拓扑不变量Wilson loop方法、表面态谱函数、反常Hall电导率。反常Hall电导率sigma_A^xy -e^2/hbar * sum_n integral dk Omega_n^z(k) f(E_n(k))其中Omega_n^z为Berry曲率z分量。两带载流子模型纵向电阻率和横向Hall电阻率n_h/n_e为载流子浓度mu_h/mu_e为迁移率Hall电导率从纵向和横向电阻率计算Debye模型热容的晶格贡献theta_D为Debye温度三、核心结果图 1Cr2Se3的晶体结构和化学表征。(a) 晶体结构示意图。(b) EDX谱确认Cr:Se40.6%:59.4%。(c) 室温XRD仅(00l)峰确认高结晶度。(d) 150K下四圆XRD的(hk0)倒空间反射展示三方对称性。(e,f) Cr 2p和Se 3d的高分辨XPS谱。Cr2Se3的晶体结构与化学态Cr2Se3为三方晶系非中心对称空间群R-3No. 148可视为刚玉型结构的畸变变体。Se原子形成近密堆积的阴离子框架Cr3占据八面体间隙位CrSe6八面体通过边共享和角共享模式连接。XPS Cr 2p谱展示Cr3在配位场中的特征多重态分裂P1-P5所有5个分量属于同一Cr3氧化态。Se 3d谱确认Se2-的阴离子态。150K到室温范围内无结构相变。图 2Cr2Se3的磁性。(a) ZFC磁化率随温度变化0.1T, ab面和c轴插图为低温放大。(b) 逆磁化率Curie-Weiss拟合。(c,d) 选定温度下的等温磁化曲线。两次反铁磁相变与磁各向异性零场冷ZFC磁化率在TN1~38K和TN2~42K处显示两次连续磁转变。TN1以下c轴磁化率显著大于ab面展示easy-axis沿c轴的单轴磁各向异性。Curie-Weiss拟合mu_eff6.6 muB (ab) / 6.7 muB (c)远大于Cr3纯自旋值3.87 muB表明轨道角动量未完全淬灭。Weiss温度theta_P~-185K (ab) / -175K (c)确认为主导反铁磁交换。等温磁化曲线呈线性7T下仍未饱和这是反铁磁基态的典型特征。加热/冷却曲线完全重合排除一级相变。图 3Cr2Se3的输运和热容。(a) 面内电阻率随温度变化插图为低温放大。(b) 零场低温热容插图为低温放大。(c) 70-250K区间的Debye模型拟合插图为Cp/T vs T^2线性拟合。(d) 磁热容Cmag和磁熵Smag。电阻率异常与热容分析电阻率随温度降低单调减小展示金属性导电行为。在TN1~38K和TN2~42K处出现两个明显异常局部极小值和极大值反映磁性与电子自由度的强耦合。热容低温拟合Cpgamma*Tbeta*T^3得到电子热容系数gamma7.12 mJ/mol K^2beta0.54 mJ/mol K^4。Debye模型拟合70-250K得到theta_D~240K。磁熵Smag在60K以上渐近饱和至约9 J/mol K接近R*ln(2S1)R*ln(4)~11.5 J/mol K的理论值确认磁自由度在顺磁态完全释放。图 4Cr2Se3的Hall效应和磁输运。(a) 纵向电阻率rho_xx和(b) 横向Hall电阻率rho_yx随磁场变化。(c,d) 载流子浓度n和迁移率mu的温度依赖性。(e) Hall电导率sigma_xy。(f) 9T下提取的Hall电导率温度依赖性。(g) 理论计算的sigma_A^xy随化学势变化。反常Hall效应的实验与理论证据两带载流子模型在2K下n_h1.21x10^24 m^-3, n_e4.59x10^23 m^-3, mu_h0.0844 m^2/Vs, mu_e0.0095 m^2/Vs。空穴载流子浓度和迁移率均高于电子。Hall电导率sigma_xy展示S形磁场依赖性偏离线性行为在高场倾向于饱和。这暗示反常Hall效应的可能存在。实验提取的AHC约72 S/cm2K, 9T与理论值约54 S/cm定性一致。理论计算确认有限的本征反常Hall电导率来源于Brillouin区中非平庸的Berry曲率分布。AHC可在三个晶轴方向探测到理论-实验的定性一致进一步支持反常Hall效应的存在。图 5Cr2Se3的电子结构和拓扑性质。(a) 无SOC的能带。(b) 含SOC的能带。(c) (111)表面沿K-Gamma-M的表面态。(d) (111)表面在费米能级处的二维谱函数A(k)。(e) R-3空间群的Brillouin区。强拓扑半金属的证据Z2(1;000)能带结构两条能带穿过费米能级导带和价带在Gamma点附近接触和重叠。引入SOC后两条能带在费米能级处分裂。表面态在(111)表面Gamma点附近出现费米面谱函数显示Gamma点周围明显的亮环展示拓扑非平庸特征。Z2拓扑不变量Wilson loop方法计算得到Z21kz0平面其他平面为零拓扑指数为(1;000)确认为强拓扑半金属。非磁相和AFM相具有相同的拓扑分类表明拓扑性质由反演对称性而非时间反演对称性保护。AFM序的额外作用AFM相中破缺的时间反演对称性使反常Hall效应和自旋依赖的表面散射成为可能磁性作为拓扑性质的功能调控参数而非前提条件。DFT Tips【DFT Tip 1】WIEN2K vs VASP全势方法的选择WIEN2K使用全势LAPW方法不引入赝势近似对含重元素体系如Se的SOC描述更精确。VASP使用PAW赝势计算效率更高但依赖于赝势质量。对于Cr2Se3这类含4p重元素Se的体系全势方法在SOC矩阵元计算上有优势。但如果需要大超胞如声子计算、缺陷计算VASP的平面波效率更高。建议如果计算重点是SOC相关的拓扑性质Berry曲率、AHC、Z2全势方法更可靠。如果计算重点是结构弛豫或大体系VASP更实用。【DFT Tip 2】Wu-Cohen GGA超越PBE的交换关联泛函Wu-Cohen GGA是非经验泛函在晶格常数和晶体结构描述上显著优于PBE。对于Cr2Se3准确的晶格常数对后续能带和拓扑性质计算至关重要。在VASP中可通过GGAWC设置。但注意Wu-Cohen在形成能、吸附能等热力学性质上的表现不一定优于PBE。对于拓扑性质计算WC和PBE通常给出定性一致的结果。建议如果材料的结构参数对拓扑性质敏感如应变诱导的拓扑相变优先使用WC或PBEsol。如果不敏感PBE即可。【DFT Tip 3】Wannier90拟合从DFT能带到紧束缚模型WIEN2K到Wannier90的接口WIEN2WANNIER需要手动设置投影轨道。Cr2Se3的投影选择Cr-d轨道5个和Se-p轨道3个共18个Wannier轨道。关键步骤(1) 确认能带解缠disentanglement窗口覆盖所有目标能带(2) 检查Wannier函数的空间局域性spread5 A^2为佳(3) 验证内插能带与DFT能带的一致性最大偏差10 meV。常见陷阱Wannier拟合的初始投影选择不当可能导致虚假的Wannier函数或错误的拓扑不变量。始终检查Wannier函数的实空间分布和对称性。【DFT Tip 4】Z2拓扑不变量的计算Wilson loop方法WannierTools使用Wilson loop方法计算Z2。关键步骤(1) 在kz0和kzpi平面各计算Wilson loop的演化(2) 数Wilson loop穿过参考线的奇偶次数(3) 得到Z2(nu0; nu1 nu2 nu3)。对于Cr2Se3Z2(1;000)表示强拓扑绝缘体/半金属nu01是关键。nu01意味着在Brillouin区中任意kz平面上的时间反演极化是非平庸的。注意Z2计算要求体系具有时间反演对称性或至少反演对称性。对于AFM体系需要确认磁构型是否保留了有效的反演对称性。【DFT Tip 5】反常Hall电导率的DFT计算AHC -e^2/hbar * sum_n integral dk Omega_n^z(k) f(E_n(k))。在Wannier90紧束缚模型基础上WannierTools通过密集k点积分计算AHC。关键参数k点密度至少200x200x200内插能量展宽smearing约0.01 eV。AHC对k点密度和展宽参数敏感需要收敛性测试。常见陷阱AHC在带隙内是量子化的整数倍e^2/h但在金属体系中是非量子化的。对于金属性Cr2Se3AHC的绝对值取决于费米能级附近的Berry曲率分布对带隙和费米能级位置敏感。【DFT Tip 6】AFM体系的磁构型设置超胞与对称性Cr2Se3的AFM磁构型Cr1(3a, 0,0,0)和Cr2(3b, 0,0,0.5)自旋平行Cr3(6c, 0,0,1/3)自旋反平行。在WIEN2K中需通过case.inorb和case.indm文件设置。注意AFM构型下Cr1和Cr2的磁矩大小可能不同由于不等价Wyckoff位置但应接近。总磁矩应为零精确补偿。如果总磁矩不为零检查磁构型设置和自洽收敛。建议先计算非磁态再添加AFM磁序。比较非磁态和AFM态的总能量确认AFM态能量更低本文中约低344 meV/atom。【DFT Tip 7】两带模型拟合Hall数据的数值技巧两带模型有4个自由参数n_h, n_e, mu_h, mu_e对rho_xx和rho_yx同时拟合。多参数模型的拟合质量可能由参数简并度而非物理精度决定。建议(1) 使用全局优化算法如模拟退火或遗传算法而非局部梯度下降(2) 进行参数不确定性分析如bootstrap或MCMC(3) 检查拟合参数的物理合理性载流子浓度应在10^23-10^25 m^-3量级迁移率应在0.001-0.1 m^2/Vs量级。注意两带模型的良好拟合不代表排除了其他模型如三带模型或包含AHE项的模型。本文作者也明确承认了这一局限性。【DFT Tip 8】热容分析电子晶格磁的分离低温热容Cpgamma*Tbeta*T^3。gamma给出电子贡献正比于DOS(E_F)beta给出晶格和反铁磁磁振子的联合贡献两者在低温下都遵循T^3。注意不能从beta单独分离晶格和磁振子贡献本文明确说明Cpgamma*Tbeta*T^3拟合仅用于可靠提取电子gamma*T项。晶格贡献通过Debye模型在更高温度70-250K独立拟合。常见错误很多初学者直接从beta计算Debye温度而不考虑磁振子贡献对于AFM材料这是错误的。知识扩展【知识扩展 1】反铁磁拓扑材料对称性保护与Z2分类【理论解释】拓扑绝缘体的Z2分类最初基于时间反演对称性TRS。在AFM体系中TRS可能被破缺但反演对称性IS和TRS的联合操作PT可以保护拓扑态。Cr2Se3的拓扑性质由IS保护而非TRS。【方法比较】TRS保护的拓扑绝缘体如Bi2Se3表面态为奇数个Dirac锥被TRS保护。IS保护的拓扑材料Z2拓扑分类在IS存在时仍有效但表面态可能被磁性破坏。AFM拓扑材料结合了两者优势体态拓扑保护磁输运自由度。【经典体系】MnBi2Te4AFM拓扑绝缘体QAHE在低温实现、Mn3Sn非共线AFM大AHE、CuMnAsAFM Dirac半金属。Cr2Se3加入这一家族但具有更简单的二元组成和更好的结构稳定性。【经典参考】Hasan Kane, RMP 82, 3045 (2010)--拓扑绝缘体综述Deng et al., Science 367, 895 (2020)--MnBi2Te4中QAHE实验Smejkal et al., PRX 12, 011028 (2022)--AFM拓扑物理。【知识扩展 2】化学气相输运CVT法生长单晶【原理】CVT利用输运剂如碘在高温区与源材料反应生成挥发性中间产物中间产物在低温区分解沉积形成单晶。输运驱动力是温度梯度导致的化学势差。【优势】相比助熔剂法和Bridgman法CVT可以在较低温度下生长高熔点材料单晶且晶体质量高、缺陷少。特别适合含硫族元素S、Se、Te的过渡金属化合物。【关键参数】输运剂种类和浓度碘是最常用的浓度约2-5 mg/cm^3、温度梯度通常50-100度C、生长时间7-14天。温度梯度过大导致多晶过小导致生长速率慢。【经典参考】Nitsche, J. Phys. Chem. Solids 17, 163 (1960)--CVT方法奠基性工作Schmidt Gruehn, Chem. Unserer Zeit 14, 15 (1980)--CVT原理综述。科研经验【科研经验 1】AFM拓扑材料的实验挑战如何区分正常Hall和反常Hall问题在AFM材料中净磁矩为零如何区分正常Hall效应NHE和反常Hall效应AHE原因AHE通常与净磁化强度成正比AFM中净磁矩为零AHE只能来自动量空间的Berry曲率。但NHE在低场下是线性的AHE在低场下可能非线性两者在实验数据中难以分离。解决方案(1) 使用两带模型或多带模型拟合场依赖的rho_xx和rho_yx同时约束NHE和AHE分量(2) 辅助DFT计算AHC提供独立的理论预测(3) 如果可能使用更高磁场14T使AHE趋于饱和便于分离。建议实验数据不足以唯一确定AHE时应明确说明不确定性。本文承认由于实验设备限制无法确认是否达到完全饱和这是诚实的科学态度值得学习。【科研经验 2】拓扑性质的完整证据链从DFT到实验问题仅凭DFT计算的Z2指数或表面态能否声称材料是拓扑材料原因DFT可能低估带隙PBE或给出错误的带序Z2不变量的计算依赖于能带拓扑对带隙闭合/打开的定性特征敏感。表面态可能来自平庸的表面悬挂键而非拓扑。解决方案(1) 使用HSE06或GW校核带隙和带序(2) 在多组Wannier投影下验证Z2的鲁棒性(3) ARPES实验验证表面态(4) 输运实验验证拓扑特征如SdH振荡、弱反局域化。建议对于Cr2Se3目前仅完成了DFT层面的拓扑分类。实验层面的ARPES和输运验证是下一步的关键。如果是我我还会继续算【继续算 1】HSE06杂化泛函验证能带和拓扑性质为什么值得算PBE/WC泛函可能低估带隙影响能带反转和拓扑分类。HSE06可以提供更准确的带隙和带序验证Z2(1;000)的鲁棒性。能回答的问题HSE06的带隙和带序是否与GGA一致Z2分类是否保持不变适合体系所有半导体/半金属拓扑材料。输入HSE06计算WIEN2K或VASP【继续算 2】ARPES可观测量的DFT模拟表面态谱函数为什么值得算WannierTools计算的表面态谱函数可以直接与ARPES实验比较。进一步计算不同光子能量和偏振下的矩阵元效应可以更真实地模拟ARPES信号。能回答的问题表面态的色散和自旋纹理是否可被ARPES分辨哪些光子能量最适合探测表面态适合体系所有拓扑材料。输入Wannier90WannierTools光电发射矩阵元计算。【继续算 3】VASP全电子计算Wannier90与WIEN2K的交叉验证为什么值得算用VASP重复DFT计算可以交叉验证结果对方法的依赖性。如果VASP和WIEN2K给出相同的拓扑分类结论更加可靠。能回答的问题PAW赝势和全势LAPW在Cr2Se3的能带拓扑上是否一致适合体系所有材料。输入VASP PBESOCWannier90。【继续算 4】声子谱计算验证结构稳定性为什么值得算本文未计算声子谱无法确认Cr2Se3的动力学稳定性。声子计算是材料预测论文的标准要求。能回答的问题Cr2Se3是否有虚频Gamma点的光学声子模式是什么适合体系所有材料。输入PhonopyVASPDFPT或有限位移法到高。【继续算 5】磁各向异性与自旋波激发磁振子谱为什么值得算Cr2Se3具有两个AFM相变和明显的磁各向异性磁振子谱可以揭示不同AFM相的磁激发特征。计算磁振子谱有助于理解磁转变的微观机制。能回答的问题AFM(L)和AFM(H)的磁振子谱有何不同磁振子能隙与磁各向异性有何关系适合体系所有磁性材料。输入VASPTB2JSpinW。【继续算 6】Berry曲率偶极矩与非互易输运为什么值得算Cr2Se3的反演对称性破缺可能导致Berry曲率偶极矩产生非线性的非互易输运效应如非线性Hall效应。这在AFM拓扑材料中是一个新兴方向。能回答的问题Cr2Se3的Berry曲率偶极矩有多大非线性Hall效应是否可观测适合体系反演对称性破缺的拓扑材料。输入Wannier90Berry曲率偶极矩计算。Qian, Shen, Zhou et al. | Phys. Rev. B 114, 024404 (2026) | 反铁磁拓扑材料 Cr2Se3 CVT WIEN2K Wannier90