交错磁体Cr2Se2O单层的超快自旋动力学

交错磁体Cr2Se2O单层的超快自旋动力学

NPJ COMPUT. MATER. 12, 241 (2026)

交错磁体Cr2Se2O单层的超快自旋动力学

Ultrafast Spin Dynamics in Altermagnetic Cr2Se2O Monolayer

导读 导读：交错磁体（Altermagnet）兼具铁磁体的自旋劈裂和反铁磁体的零净磁矩，是下一代自旋电子学的理想平台。本文通过rt-TDDFT+NAMD联合方法，首次模拟了单层Cr2Se2O在飞秒激光激发下的自旋动力学。核心发现：交错磁体相的自旋弛豫时间（718 fs）显著长于传统反铁磁相（202 fs和498 fs），根本原因在于SOC矩阵元受Neel型磁序对称性约束而减弱，抑制了自旋翻转过程。

一、前言背景

交错磁体：二维自旋电子学的新范式

交错磁体（Altermagnet, AM）是2022年以来凝聚态物理最受关注的发现之一。它同时兼具铁磁体的自旋劈裂能带和反铁磁体的零净磁矩，在THz级自旋动力学、无杂散场器件集成等方面具有天然优势。

将AM推向二维极限是当前的重要方向。单层Cr2Se2O（ML-Cr2Se2O）是近期提出的二维AM候选材料，具有P4/mmm空间群对称性，C4z和T分别破缺但C4zT联合对称性保留，导致d波形式的非相对论自旋劈裂。

本文通过rt-TDDFT + NAMD联合方法，首次从第一性原理模拟了ML-Cr2Se2O在超快飞秒激光激发下的自旋动力学，揭示了AM相与常规AFM相自旋弛豫时间的本质差异。

超快自旋动力学：从激光激发到自旋弛豫

核心问题：飞秒激光脉冲如何操控二维交错磁体的自旋？AM相的自旋弛豫机制与常规反铁磁体有何不同？

方法体系：VASP（DFT+U，基态电子结构） -> ELK（rt-TDDFT，超快激光激发） -> Hefei-NAMD（非绝热分子动力学，自旋弛豫过程） -> EspinS（Monte Carlo，磁转变温度）。

关键发现：(1) AM相自旋弛豫时间718 fs显著长于AFM2相（202 fs）和AFM3相（498 fs）；(2) AM相中SOC效应减弱，抑制了自旋翻转过程，延长了弛豫时间。

Cr2Se2O超快自旋动力学研究流程。VASP DFT+U（基态电子结构、MAE、交换耦合、Monte Carlo） -> ELK rt-TDDFT（超快激光激发，飞秒退磁化） -> Hefei-NAMD（非绝热分子动力学，自旋弛豫过程）。核心发现：AM相自旋弛豫时间718 fs显著长于AFM2相（202 fs）和AFM3相（498 fs）。

二、研究方法

VASP DFT+U：基态电子结构与磁性

计算设置：PBE泛函 + PAW赝势, ENCUT=600 eV, 15x15x1 k点（弛豫）, 25x25x1（电子结构）, 真空层15 A。Cr原子Hubbard Ueff=4 eV（与文献一致）。MAE通过非自洽SOC计算总能量差得到。

磁基态确定：在2x2x1超胞中对比FM、AFM1、AFM2、AFM3四种磁构型。AFM1为基态（能量最低），每Cr原子磁矩3.27 muB，与Cr3+（S=3/2）一致。DMI被Moriya对称性规则禁止，不考虑非共线磁态。

交换耦合：Heisenberg模型拟合，J1=-12.8 meV（最近邻Cr-Cr直接交换，强AFM），J2=3.9 meV（Cr-Se-Cr超交换，~90度，FM），J2'=-6.4 meV（Cr-O-Cr超交换，180度，AFM）。竞争耦合导致一定程度的磁阻挫。

Heisenberg自旋哈密顿量：J_{ij}为交换耦合参数，A为单离子各向异性

rt-TDDFT + NAMD：超快自旋动力学模拟

ELK rt-TDDFT：在LSDA+U水平上求解含时Kohn-Sham方程，引入线偏振激光脉冲（沿x轴或沿z轴），FWHM=6 fs, 时间步长0.0048 fs, k点5x5x1。模拟飞秒激光激发下的磁矩演化。

Hefei-NAMD：300K 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分量期望值

三、核心结果

图 1：ML-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 eV（spin-up和spin-down通道相同）。沿Gamma-M方向能带简并，但其他高对称点出现明显自旋劈裂。劈裂幅度约607 meV，与MnTe（1100 meV）、CrSb（705 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 muB（11 fs处），随后稳定。在~34 fs处出现二次退磁化（减至~1.17 muB），伴随Se/O磁矩的显著增强振荡。

物理机制：激光极化方向决定了对称性破缺的类型，不同极化方向下轨道-激光耦合的矩阵元不同，导致退磁化程度和自旋注入效率的显著差异。

图 3：NAMD模拟中的自旋动力学。(a,b) AM相中EPC和SOC矩阵元，数值越大耦合越强。(c) AM相中平均电子能量演化，插图为激发电子的自旋占据。(d,e) AM相和AFM2相的光激发自旋载流子动力学过程示意图。

自旋弛豫机制：AM vs AFM的对比分析

AM相的三条弛豫通道：Process 1（Cr-d^到Cr-d^, d-d跃迁，EPC=4.20 meV + SOC=8.4 meV）, Process 2（Cr-d^到Se-pv, 自旋翻转+电荷转移，SOC=6.59 meV主导）, Process 3（Se-pv到Se-pv, p-p跃迁，EPC=53.28 meV主导）。

AFM2相的四条弛豫通道：Process 1-3与AM类似，但新增Process 4（Cr-d^直接到终态Se-pv，SOC=29.27+27.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 Qb@ll

ELK是全势LAPW方法的rt-TDDFT实现，精度高但。对于2D材料，5x5x1的k点采样在rt-TDDFT中已经足够，因为含时演化对k点采样敏感度低于静态计算。

激光参数设置：FWHM=6 fs对应约0.3 eV的频谱宽度，频率0.2 a.u.（约5.44 eV）需匹配材料带隙。时间步长Deltat=0.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) 比较不同磁构型的总能量。

对于Cr2Se2O：2x2x1超胞中测试4种磁构型（FM, AFM1, AFM2, AFM3），AFM1（AM态）能量最低。DMI被Moriya规则禁止，排除非共线磁态。

常见错误：直接用FM初始磁矩计算AM材料，导致收敛到错误的磁态。务必手动设置MAGMOM标签。

【DFT Tip 4】Hubbard U对AM自旋劈裂的影响

Cr2Se2O使用Ueff=4 eV，这是Cr3+的典型值。但U值的选择会显著影响自旋劈裂的大小和带隙。建议：(1) 先用PBE（不加U）计算参考能带；(2) 用线性响应方法确定U；(3) 比较PBE和PBE+U的自旋劈裂模式。

注意：U值过大可能导致虚假的Mott绝缘态，掩盖AM金属性行为。对于Cr2Se2O（半导体），U=4 eV是合理的，但需要验证U值+/-1 eV范围内的结果是否定性一致。

不同赝势、不同版本VASP的U值不可直接迁移。即使同一元素，在不同化学环境（配位数、氧化态）中U值也可能不同。

【DFT Tip 5】MAE计算：SOC收敛与数值精度

MAE = E_theta - E_001，是总能量之差，通常为meV量级，对数值精度要求极高。建议：(1) 使用高精度收敛标准（EDIFF=1e-7或1e-8）；(2) 增加k点密度（至少25x25x1）；(3) 在非自洽SOC计算中读取自洽的CHGCAR。

Cr2Se2O的MAE=0.43 meV/atom，属于典型的3d过渡金属磁各向异性能量级。在SOC计算中，磁矩方向通过SAXIS标签控制。

常见陷阱：SOC自洽计算的能量比非SOC更高，这是正常的--SOC打开了额外的耦合通道。不要因为SOC能量更高就认为计算错误。

【DFT Tip 6】Monte Carlo模拟磁转变温度：EspinS

Cr2Se2O使用EspinS包计算Neel温度T_AM=160 K。EspinS基于经典Heisenberg模型，输入为DFT计算的交换耦合参数J_{ij}。

关键设置：20x20x1超胞（400个自旋），自旋在空间中随机旋转，Metropolis算法。T_AM的准确性取决于J_{ij}的精度和截断半径。

注意：经典Monte Carlo通常高估T_c，因为忽略了量子涨落。对于S=3/2的Cr3+，量子修正约为10-20%。实际T_AM可能低于160 K。与实验比较时需考虑这一系统偏差。

【DFT Tip 7】超交换耦合的GKA规则与DFT验证

Cr2Se2O中J1=-12.8 meV（Cr-Cr直接交换，AFM），J2=3.9 meV（Cr-Se-Cr，~90度，FM），J2'=-6.4 meV（Cr-O-Cr，180度，AFM）。这符合Goodenough-Kanamori-Anderson（GKA）规则：~90度超交换倾向于FM，180度超交换倾向于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-TDDFT+NAMD联合方案是理想选择：rt-TDDFT模拟激光激发（~10 fs），NAMD模拟后续自旋弛豫（~1 ps）。两者互补，覆盖了从激发到弛豫的完整时间尺度。

【经典参考】Dewhurst et al., Nano Lett. 18, 1842 (2018)--ELK rt-TDDFT；Zheng 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) 脉宽：FWHM=6-10 fs是典型值，过短则频谱太宽，过长则无法区分激发和弛豫过程。

建议：在论文中明确说明激光参数是"为捕获定性特征而选择"，而非直接对应实验条件。同时进行参数敏感性分析，验证主要结论对激光参数不敏感。

【科研经验 2】AM材料自旋弛豫时间的DFT预测：从定性到定量

问题：DFT+NAMD预测的自旋弛豫时间（718 fs）是否可靠？与实验比较时应注意什么？

原因：(1) NAMD基于经典核近似，忽略了核量子效应（如零点运动、隧穿）；(2) 表面跳跃算法忽略了电子-电子散射和退相干效应；(3) DFT带隙通常被低估（PBE），影响激发态能量和弛豫通道。

解决方案：(1) 使用HSE06杂化泛函校核带隙，评估带隙误差对弛豫时间的影响；(2) 比较不同温度（100K、300K、500K）下的弛豫时间，评估温度依赖性；(3) 从定性角度解读结果：AM > AFM3 > AFM2的弛豫时间排序比绝对值更有意义。

建议：将DFT+NAMD的弛豫时间视为定性趋势而非定量预测。与实验比较时，关注排序而非绝对值。如果实验能确认AM相弛豫时间确实长于AFM相，就是对理论的最大支持。

Guo, Gao, Zhao, Jiang | npj Comput. Mater. 12, 241 (2026) | 交错磁体 超快自旋动力学 rt-TDDFT NAMD