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

 

PHYS. REV. B 114, 024404 (2026)

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

Synthesis, Physical Properties, and Band Structure of Antiferromagnetic Cr2Se3

导读 导读：Cr2Se3是一种具有非中心对称三方结构的反铁磁拓扑材料。本文通过实验合成（CVT单晶）与第一性原理计算（WIEN2K+Wannier90）相结合，系统研究了其结构、磁性、输运和拓扑性质。核心发现：(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描述上具有优势。

figure

Cr2Se3反铁磁拓扑材料研究流程。CVT单晶生长 -> 结构/磁性/输运/热容全面表征 -> WIEN2K DFT（全势LAPW, Wu-Cohen GGA+SOC） -> Wannier90+WannierTools（拓扑分类、表面态、反常Hall电导率）。核心发现：两次AFM相变 + Z2=(1;000)强拓扑半金属 + 反常Hall效应。

二、研究方法

实验方法：CVT单晶生长与多维度物性表征

晶体生长：化学气相输运（CVT）法，Cr:Se=2:3摩尔比，碘为输运剂。双温区炉：800/700度C，2天升温+10天恒温+自然冷却。获得具有明显金属光泽的单晶。

结构表征：室温XRD（Rigaku, Cu Kalpha），四圆单晶XRD（Bruker D8 Venture, 150K, Mo Kalpha），空间群R-3（No. 148），a=b=6.2544(2) A, c=17.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分量。

formula

两带载流子模型：纵向电阻率和横向Hall电阻率，n_h/n_e为载流子浓度，mu_h/mu_e为迁移率

formulaformula

Hall电导率从纵向和横向电阻率计算

formula

Debye模型：热容的晶格贡献，theta_D为Debye温度

三、核心结果

figure

图 1：Cr2Se3的晶体结构和化学表征。(a) 晶体结构示意图。(b) EDX谱确认Cr:Se=40.6%:59.4%。(c) 室温XRD，仅(00l)峰，确认高结晶度。(d) 150K下四圆XRD的(hk0)倒空间反射，展示三方对称性。(e,f) Cr 2p和Se 3d的高分辨XPS谱。

Cr2Se3的晶体结构与化学态

Cr2Se3为三方晶系，非中心对称空间群R-3（No. 148），可视为刚玉型结构的畸变变体。Se原子形成近密堆积的阴离子框架，Cr3+占据八面体间隙位，CrSe6八面体通过边共享和角共享模式连接。

XPS Cr 2p谱展示Cr3+在配位场中的特征多重态分裂（P1-P5），所有5个分量属于同一Cr3+氧化态。Se 3d谱确认Se2-的阴离子态。150K到室温范围内无结构相变。

figure

图 2：Cr2Se3的磁性。(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_eff=6.6 muB (ab) / 6.7 muB (c)，远大于Cr3+纯自旋值3.87 muB，表明轨道角动量未完全淬灭。Weiss温度theta_P~-185K (ab) / -175K (c)，确认为主导反铁磁交换。

等温磁化曲线呈线性，7T下仍未饱和，这是反铁磁基态的典型特征。加热/冷却曲线完全重合，排除一级相变。

figure

图 3：Cr2Se3的输运和热容。(a) 面内电阻率随温度变化，插图为低温放大。(b) 零场低温热容，插图为低温放大。(c) 70-250K区间的Debye模型拟合，插图为Cp/T vs T^2线性拟合。(d) 磁热容Cmag和磁熵Smag。

电阻率异常与热容分析

电阻率随温度降低单调减小，展示金属性导电行为。在TN1~38K和TN2~42K处出现两个明显异常（局部极小值和极大值），反映磁性与电子自由度的强耦合。

热容低温拟合：Cp=gamma*T+beta*T^3，得到电子热容系数gamma=7.12 mJ/mol K^2，beta=0.54 mJ/mol K^4。Debye模型拟合（70-250K）得到theta_D~240K。

磁熵Smag在60K以上渐近饱和至约9 J/mol K，接近R*ln(2S+1)=R*ln(4)~11.5 J/mol K的理论值，确认磁自由度在顺磁态完全释放。

figure

图 4：Cr2Se3的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_h=1.21x10^24 m^-3, n_e=4.59x10^23 m^-3, mu_h=0.0844 m^2/Vs, mu_e=0.0095 m^2/Vs。空穴载流子浓度和迁移率均高于电子。

Hall电导率sigma_xy展示"S形"磁场依赖性，偏离线性行为，在高场倾向于饱和。这暗示反常Hall效应的可能存在。实验提取的AHC约72 S/cm（2K, 9T），与理论值约54 S/cm定性一致。

理论计算确认有限的本征反常Hall电导率，来源于Brillouin区中非平庸的Berry曲率分布。AHC可在三个晶轴方向探测到，理论-实验的定性一致进一步支持反常Hall效应的存在。

figure

图 5：Cr2Se3的电子结构和拓扑性质。(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方法计算得到Z2=1（kz=0平面），其他平面为零，拓扑指数为(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中可通过GGA=WC设置。但注意：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函数的空间局域性（spread<5 A^2为佳）；(3) 验证内插能带与DFT能带的一致性（最大偏差<10 meV）。

常见陷阱：Wannier拟合的初始投影选择不当可能导致虚假的Wannier函数或错误的拓扑不变量。始终检查Wannier函数的实空间分布和对称性。

【DFT Tip 4】Z2拓扑不变量的计算：Wilson loop方法

WannierTools使用Wilson loop方法计算Z2。关键步骤：(1) 在kz=0和kz=pi平面各计算Wilson loop的演化；(2) 数Wilson loop穿过参考线的奇偶次数；(3) 得到Z2=(nu0; nu1 nu2 nu3)。

对于Cr2Se3，Z2=(1;000)表示强拓扑绝缘体/半金属，nu0=1是关键。nu0=1意味着在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），但在金属体系中是非量子化的。对于金属性Cr2Se3，AHC的绝对值取决于费米能级附近的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】热容分析：电子+晶格+磁的分离

低温热容Cp=gamma*T+beta*T^3。gamma给出电子贡献（正比于DOS(E_F)），beta给出晶格和反铁磁磁振子的联合贡献（两者在低温下都遵循T^3）。

注意：不能从beta单独分离晶格和磁振子贡献！本文明确说明Cp=gamma*T+beta*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拓扑材料结合了两者优势：体态拓扑保护+磁输运自由度。

【经典体系】MnBi2Te4（AFM拓扑绝缘体，QAHE在低温实现）、Mn3Sn（非共线AFM，大AHE）、CuMnAs（AFM 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分辨？哪些光子能量最适合探测表面态？

适合体系：所有拓扑材料。输入：Wannier90+WannierTools+光电发射矩阵元计算。

【继续算 3】VASP全电子计算+Wannier90：与WIEN2K的交叉验证

为什么值得算：用VASP重复DFT计算可以交叉验证结果对方法的依赖性。如果VASP和WIEN2K给出相同的拓扑分类，结论更加可靠。

能回答的问题：PAW赝势和全势LAPW在Cr2Se3的能带拓扑上是否一致？

适合体系：所有材料。输入：VASP PBE+SOC+Wannier90。

【继续算 4】声子谱计算：验证结构稳定性

为什么值得算：本文未计算声子谱，无法确认Cr2Se3的动力学稳定性。声子计算是材料预测论文的标准要求。

能回答的问题：Cr2Se3是否有虚频？Gamma点的光学声子模式是什么？

适合体系：所有材料。输入：Phonopy+VASP（DFPT或有限位移法）到高。

【继续算 5】磁各向异性与自旋波激发：磁振子谱

为什么值得算：Cr2Se3具有两个AFM相变和明显的磁各向异性，磁振子谱可以揭示不同AFM相的磁激发特征。计算磁振子谱有助于理解磁转变的微观机制。

能回答的问题：AFM(L)和AFM(H)的磁振子谱有何不同？磁振子能隙与磁各向异性有何关系？

适合体系：所有磁性材料。输入：VASP+TB2J+SpinW。

【继续算 6】Berry曲率偶极矩与非互易输运

为什么值得算：Cr2Se3的反演对称性破缺可能导致Berry曲率偶极矩，产生非线性的非互易输运效应（如非线性Hall效应）。这在AFM拓扑材料中是一个新兴方向。

能回答的问题：Cr2Se3的Berry曲率偶极矩有多大？非线性Hall效应是否可观测？

适合体系：反演对称性破缺的拓扑材料。输入：Wannier90+Berry曲率偶极矩计算。

Qian, Shen, Zhou et al. | Phys. Rev. B 114, 024404 (2026) | 反铁磁拓扑材料 Cr2Se3 CVT WIEN2K Wannier90