交错磁体Cr2Se2O单层的超快自旋动力学NPJ COMPUT. MATER. 12, 241 (2026)交错磁体Cr2Se2O单层的超快自旋动力学Ultrafast Spin Dynamics in Altermagnetic Cr2Se2O Monolayer导读 导读交错磁体Altermagnet兼具铁磁体的自旋劈裂和反铁磁体的零净磁矩是下一代自旋电子学的理想平台。本文通过rt-TDDFTNAMD联合方法首次模拟了单层Cr2Se2O在飞秒激光激发下的自旋动力学。核心发现交错磁体相的自旋弛豫时间718 fs显著长于传统反铁磁相202 fs和498 fs根本原因在于SOC矩阵元受Neel型磁序对称性约束而减弱抑制了自旋翻转过程。一、前言背景交错磁体二维自旋电子学的新范式交错磁体Altermagnet, AM是2022年以来凝聚态物理最受关注的发现之一。它同时兼具铁磁体的自旋劈裂能带和反铁磁体的零净磁矩在THz级自旋动力学、无杂散场器件集成等方面具有天然优势。将AM推向二维极限是当前的重要方向。单层Cr2Se2OML-Cr2Se2O是近期提出的二维AM候选材料具有P4/mmm空间群对称性C4z和T分别破缺但C4zT联合对称性保留导致d波形式的非相对论自旋劈裂。本文通过rt-TDDFT NAMD联合方法首次从第一性原理模拟了ML-Cr2Se2O在超快飞秒激光激发下的自旋动力学揭示了AM相与常规AFM相自旋弛豫时间的本质差异。超快自旋动力学从激光激发到自旋弛豫核心问题飞秒激光脉冲如何操控二维交错磁体的自旋AM相的自旋弛豫机制与常规反铁磁体有何不同方法体系VASPDFTU基态电子结构 - ELKrt-TDDFT超快激光激发 - Hefei-NAMD非绝热分子动力学自旋弛豫过程 - EspinSMonte Carlo磁转变温度。关键发现(1) AM相自旋弛豫时间718 fs显著长于AFM2相202 fs和AFM3相498 fs(2) AM相中SOC效应减弱抑制了自旋翻转过程延长了弛豫时间。Cr2Se2O超快自旋动力学研究流程。VASP DFTU基态电子结构、MAE、交换耦合、Monte Carlo - ELK rt-TDDFT超快激光激发飞秒退磁化 - Hefei-NAMD非绝热分子动力学自旋弛豫过程。核心发现AM相自旋弛豫时间718 fs显著长于AFM2相202 fs和AFM3相498 fs。二、研究方法VASP DFTU基态电子结构与磁性计算设置PBE泛函 PAW赝势, ENCUT600 eV, 15x15x1 k点弛豫, 25x25x1电子结构, 真空层15 A。Cr原子Hubbard Ueff4 eV与文献一致。MAE通过非自洽SOC计算总能量差得到。磁基态确定在2x2x1超胞中对比FM、AFM1、AFM2、AFM3四种磁构型。AFM1为基态能量最低每Cr原子磁矩3.27 muB与Cr3S3/2一致。DMI被Moriya对称性规则禁止不考虑非共线磁态。交换耦合Heisenberg模型拟合J1-12.8 meV最近邻Cr-Cr直接交换强AFMJ23.9 meVCr-Se-Cr超交换~90度FMJ2-6.4 meVCr-O-Cr超交换180度AFM。竞争耦合导致一定程度的磁阻挫。Heisenberg自旋哈密顿量J_{ij}为交换耦合参数A为单离子各向异性rt-TDDFT NAMD超快自旋动力学模拟ELK rt-TDDFT在LSDAU水平上求解含时Kohn-Sham方程引入线偏振激光脉冲沿x轴或沿z轴FWHM6 fs, 时间步长0.0048 fs, k点5x5x1。模拟飞秒激光激发下的磁矩演化。Hefei-NAMD300K AIMD轨迹5 ps, 1 fs步长取最后2 ps计算非绝热耦合NAC。在自旋非绝热spin-diabatic表象中展开Kohn-Sham轨道同时求解EPC和SOC驱动的电子跃迁。自旋动力学物理图像激光激发后光生载流子从Cr-d(CBM)弛豫到Se-p(VBM)过程中SOC驱动自旋翻转。AM相中SOC矩阵元受Neel型磁序对称性约束而减弱导致弛豫时间延长。含时Kohn-Sham方程A_ext(t)为激光矢势B_s为有效磁场外场XC场sigma为Pauli矩阵NAMD含时薛定谔方程H_0为无自旋哈密顿量H_SOC为自旋轨道耦合项自旋布居与平均电子能量rho_j为能带占据数sigma_{j,z}为自旋z分量期望值三、核心结果图 1ML-Cr2Se2O的晶体结构、磁性和电子性质。(a) 俯视图和侧视图蓝色/粉色/黄色球分别表示Cr/Se/O原子alpha和beta标记两个相反自旋的Cr。(b) MAE的角度依赖性展示easy-axis在[001]方向。(c) 自旋极化能带结构展示d波形式的自旋劈裂。(d-f) Cr-alpha, Cr-beta, Se, O原子的轨道投影DOS。能带结构与d波自旋劈裂ML-Cr2Se2O为间接带隙半导体带隙1.01 eVspin-up和spin-down通道相同。沿Gamma-M方向能带简并但其他高对称点出现明显自旋劈裂。劈裂幅度约607 meV与MnTe1100 meV、CrSb705 meV等已知AM材料相当。自旋劈裂的对称性分析C4z和T分别破缺C4zT保留。能带劈裂呈现d波特征4个节点特定k方向的spin-up和spin-down能带交换。PDOS分析揭示Cr-3d_{yz}alpha子晶格与Se-p_{xy}/O-p_z耦合Cr-3d_{xz}beta子晶格与Se-p_{xy}/O-p_z耦合轨道间耦合的不对称性导致了AM序。图 2超快激光诱导磁化动力学的时间演化。(a,b) 垂直极化激光极化方向垂直于自旋下Cr、Se、O原子的磁矩演化。(c,d) 平行极化激光下的磁矩演化。展示激光极化方向对退磁化行为的显著影响。激光诱导的退磁化与自旋转移垂直极化激光极化垂直于自旋Cr原子在激光脉冲持续时间内经历完全退磁化磁矩接近零。Se和O原子通过超交换和电荷转移从Cr原子获得瞬态磁矩~0.11 muB激光关闭后迅速衰减。平行极化激光极化平行于自旋Cr原子部分退磁化至~1.44 muB11 fs处随后稳定。在~34 fs处出现二次退磁化减至~1.17 muB伴随Se/O磁矩的显著增强振荡。物理机制激光极化方向决定了对称性破缺的类型不同极化方向下轨道-激光耦合的矩阵元不同导致退磁化程度和自旋注入效率的显著差异。图 3NAMD模拟中的自旋动力学。(a,b) AM相中EPC和SOC矩阵元数值越大耦合越强。(c) AM相中平均电子能量演化插图为激发电子的自旋占据。(d,e) AM相和AFM2相的光激发自旋载流子动力学过程示意图。自旋弛豫机制AM vs AFM的对比分析AM相的三条弛豫通道Process 1Cr-d^到Cr-d^, d-d跃迁EPC4.20 meV SOC8.4 meV, Process 2Cr-d^到Se-pv, 自旋翻转电荷转移SOC6.59 meV主导, Process 3Se-pv到Se-pv, p-p跃迁EPC53.28 meV主导。AFM2相的四条弛豫通道Process 1-3与AM类似但新增Process 4Cr-d^直接到终态Se-pvSOC29.2727.90 meV加速弛豫且EPC在Process 1和3中更强18.10 meV和14.36 meV导致弛豫时间202 fs显著短于AM相718 fs。AM相弛豫时间延长的根本原因Neel型磁序保留了更多对称性约束了SOC矩阵元降低了自旋翻转效率。AFM3相498 fs介于两者之间进一步验证了磁序-对称性-SOC-弛豫时间的关联。DFT Tips【DFT Tip 1】rt-TDDFT计算ELK vs Octopus vs QbllELK是全势LAPW方法的rt-TDDFT实现精度高但。对于2D材料5x5x1的k点采样在rt-TDDFT中已经足够因为含时演化对k点采样敏感度低于静态计算。激光参数设置FWHM6 fs对应约0.3 eV的频谱宽度频率0.2 a.u.约5.44 eV需匹配材料带隙。时间步长Deltat0.0048 fs约0.2 a.u.是rt-TDDFT的典型值满足数值稳定性。常见陷阱rt-TDDFT中激光场强不能直接与实验比较。模拟中故意使用强场以确保足够的电子激发目的是捕获定性特征而非定量预测。不要在论文中声称模拟的激光强度与实验一致。【DFT Tip 2】NAMD中的自旋非绝热表象与SOC处理Hefei-NAMD是少数支持自旋非绝热表象的NAMD代码。在spin-diabatic表象中波函数用自旋极化的Kohn-Sham轨道展开SOC矩阵元直接在自旋空间中计算从VASP的PROCAR读取sigma_{x,y,z}。关键参数AIMD轨迹长度至少2 ps用于NAC计算时间步长1 fs温度300 K。NAC从波函数重叠计算SOC矩阵元从PROCAR提取。自旋布居定义为 rho_j * sigma_{j,z}。注意NAMD基于经典核近似表面跳跃未考虑核量子效应和电子-电子散射。在百飞秒时间尺度上这些效应可能重要但当前方法难以精确捕获。【DFT Tip 3】交错磁体磁构型的DFT预处理在计算AM材料前必须确认磁构型满足自旋空间群对称性。步骤(1) 在非磁态弛豫结构(2) 根据对称性确定可能的磁构型(3) 在INCAR中设置MAGMOM每个原子指定自旋方向和大小(4) 比较不同磁构型的总能量。对于Cr2Se2O2x2x1超胞中测试4种磁构型FM, AFM1, AFM2, AFM3AFM1AM态能量最低。DMI被Moriya规则禁止排除非共线磁态。常见错误直接用FM初始磁矩计算AM材料导致收敛到错误的磁态。务必手动设置MAGMOM标签。【DFT Tip 4】Hubbard U对AM自旋劈裂的影响Cr2Se2O使用Ueff4 eV这是Cr3的典型值。但U值的选择会显著影响自旋劈裂的大小和带隙。建议(1) 先用PBE不加U计算参考能带(2) 用线性响应方法确定U(3) 比较PBE和PBEU的自旋劈裂模式。注意U值过大可能导致虚假的Mott绝缘态掩盖AM金属性行为。对于Cr2Se2O半导体U4 eV是合理的但需要验证U值/-1 eV范围内的结果是否定性一致。不同赝势、不同版本VASP的U值不可直接迁移。即使同一元素在不同化学环境配位数、氧化态中U值也可能不同。【DFT Tip 5】MAE计算SOC收敛与数值精度MAE E_theta - E_001是总能量之差通常为meV量级对数值精度要求极高。建议(1) 使用高精度收敛标准EDIFF1e-7或1e-8(2) 增加k点密度至少25x25x1(3) 在非自洽SOC计算中读取自洽的CHGCAR。Cr2Se2O的MAE0.43 meV/atom属于典型的3d过渡金属磁各向异性能量级。在SOC计算中磁矩方向通过SAXIS标签控制。常见陷阱SOC自洽计算的能量比非SOC更高这是正常的--SOC打开了额外的耦合通道。不要因为SOC能量更高就认为计算错误。【DFT Tip 6】Monte Carlo模拟磁转变温度EspinSCr2Se2O使用EspinS包计算Neel温度T_AM160 K。EspinS基于经典Heisenberg模型输入为DFT计算的交换耦合参数J_{ij}。关键设置20x20x1超胞400个自旋自旋在空间中随机旋转Metropolis算法。T_AM的准确性取决于J_{ij}的精度和截断半径。注意经典Monte Carlo通常高估T_c因为忽略了量子涨落。对于S3/2的Cr3量子修正约为10-20%。实际T_AM可能低于160 K。与实验比较时需考虑这一系统偏差。【DFT Tip 7】超交换耦合的GKA规则与DFT验证Cr2Se2O中J1-12.8 meVCr-Cr直接交换AFMJ23.9 meVCr-Se-Cr~90度FMJ2-6.4 meVCr-O-Cr180度AFM。这符合Goodenough-Kanamori-AndersonGKA规则~90度超交换倾向于FM180度超交换倾向于AFM。DFT中提取交换耦合参数的方法(1) 构建不同磁构型的超胞(2) 计算总能量差(3) 映射到Heisenberg模型解线性方程组得到J_{ij}。常见陷阱超胞不够大J_{ij}受周期性镜像相互作用影响。建议至少使用包含最近邻和次近邻交换路径的超胞大小。【DFT Tip 8】SOC对AM自旋弛豫的对称性约束AM相中SOC矩阵元受Neel型磁序对称性约束这是AM自旋弛豫时间延长的根本原因。在自旋群理论中AM相保留了更多幺正对称性操作这些操作约束了SOC在自旋空间中的矩阵元。验证方法在DFT中计算SOC矩阵元比较AM相和AFM相中自旋翻转通道的SOC强度。VASP中通过LSORBIT.TRUE.和LORBMOM.TRUE.输出轨道磁矩和SOC矩阵元。注意SOC矩阵元不仅是能量依赖的也是k依赖的。自旋弛豫过程涉及能带边缘的特定k点而非整个Brillouin区的平均。知识扩展【知识扩展 1】交错磁性Altermagnetism的理论基础【理论解释】传统磁性分类基于自旋排列的实空间对称性FM平行排列净磁矩不为零AFM反平行排列净磁矩为零。交错磁体突破这一范式实空间自旋反平行净磁矩为零但倒空间能带出现非相对论自旋劈裂。【对称性根源】AM序来源于自旋空间群spin group而非磁空间群magnetic group。在自旋群描述中自旋空间和实空间的操作是独立的允许非平庸的联合操作如C4T产生动量依赖的自旋劈裂。自旋劈裂的形式d波、g波、i波由磁对称性决定。【发现历程】2021年Smejkal、Sinova、Jungwirth在PRB提出altermagnetism概念2022年PRX发表AM综述2024年Nature发表MnTe和RuO2的实验验证。至今已有数百种AM候选材料被理论预测。【经典参考】Smejkal et al., PRX 12, 040501 (2022)Krempasky et al., Nature 626, 517 (2024)Zhu et al., Nature 626, 523 (2024)Liu et al., PRX 12, 021016 (2022)。【迁移能力】AM概念广泛适用于2D材料Cr2Se2O、V2Se2O、Kagome体系KV2Se2O、3D材料MnTe、RuO2、CrSb和f电子体系。【知识扩展 2】rt-TDDFT与NAMD超快动力学的两大计算范式【理论解释】rt-TDDFT直接求解含时Kohn-Sham方程描述电子密度在激光场中的实时演化可捕获亚飞秒尺度的电子动力学但通常冻结核运动。NAMD通过表面跳跃surface hopping算法同时处理电子和核的动力学可以描述电子-声子耦合驱动的弛豫过程。【方法比较】rt-TDDFT优势在于完全第一性原理、直接捕获激光-物质相互作用局限在于需飞秒级时间步长、核运动通常冻结。NAMD优势在于同时处理电子核动力学、可模拟皮秒级弛豫过程局限在于经典核近似、表面跳跃算法的近似性。【联合使用】本文的rt-TDDFTNAMD联合方案是理想选择rt-TDDFT模拟激光激发~10 fsNAMD模拟后续自旋弛豫~1 ps。两者互补覆盖了从激发到弛豫的完整时间尺度。【经典参考】Dewhurst et al., Nano Lett. 18, 1842 (2018)--ELK rt-TDDFTZheng et al., WIREs Comput. Mol. Sci. 9, e1411 (2019)--Hefei-NAMD综述Tully, J. Chem. Phys. 93, 1061 (1990)--表面跳跃算法。科研经验【科研经验 1】rt-TDDFT模拟中的激光参数选择问题rt-TDDFT中激光场强、频率、脉宽如何选择模拟结果与实验的可比性如何原因rt-TDDFT中的激光场强通常远大于实验值因为模拟需要足够的电子激发以观察定性效应。频率需匹配材料带隙。脉宽决定了频谱宽度和激发的时间分辨率。解决方案(1) 频率选择略大于带隙的光子能量确保价带电子能激发到导带(2) 场强从低到高扫描选择能产生足够激发但不会导致数值不稳定的场强(3) 脉宽FWHM6-10 fs是典型值过短则频谱太宽过长则无法区分激发和弛豫过程。建议在论文中明确说明激光参数是为捕获定性特征而选择而非直接对应实验条件。同时进行参数敏感性分析验证主要结论对激光参数不敏感。【科研经验 2】AM材料自旋弛豫时间的DFT预测从定性到定量问题DFTNAMD预测的自旋弛豫时间718 fs是否可靠与实验比较时应注意什么原因(1) NAMD基于经典核近似忽略了核量子效应如零点运动、隧穿(2) 表面跳跃算法忽略了电子-电子散射和退相干效应(3) DFT带隙通常被低估PBE影响激发态能量和弛豫通道。解决方案(1) 使用HSE06杂化泛函校核带隙评估带隙误差对弛豫时间的影响(2) 比较不同温度100K、300K、500K下的弛豫时间评估温度依赖性(3) 从定性角度解读结果AM AFM3 AFM2的弛豫时间排序比绝对值更有意义。建议将DFTNAMD的弛豫时间视为定性趋势而非定量预测。与实验比较时关注排序而非绝对值。如果实验能确认AM相弛豫时间确实长于AFM相就是对理论的最大支持。Guo, Gao, Zhao, Jiang | npj Comput. Mater. 12, 241 (2026) | 交错磁体 超快自旋动力学 rt-TDDFT NAMD