交错磁体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）。

计算设置：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。

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 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模拟中的激光参数选择

人工智能

交错磁体Cr2Se2O单层超快自旋动力学研究显示，AM相自旋弛豫时间（718fs）显著长于反铁磁相（202fs和498fs），源于Neel型磁序对称性约束SOC矩阵元，抑制自旋翻转。飞秒激光极化方向显著影响退磁化过程，为二维自旋电子学器件设计提供新机制。

文章

问答

视频

教程

学习中心

腾讯云实验室

直播

竞赛

腾讯云代码分析专区

腾讯iOA零信任安全管理系统专区

腾讯云架构师技术同盟交流圈

腾讯云数据库专区

腾讯云智能顾问专区

腾讯云原生专区

腾讯混元专区

腾讯云TCE专区

腾讯云Lighthouse专区

腾讯云HAI专区

腾讯云Edgeone专区

腾讯云存储专区

腾讯云智能专区

腾讯轻联专区 

腾讯云开发专区

TAPD专区

腾讯轻量云游戏服专区

EdgeOne AI 安全实战专区

腾讯云最具价值专家

腾讯云架构师技术同盟

腾讯云创作之星

腾讯云开发者先锋

腾讯云代码助手

云原生构建

TAPD 敏捷项目管理

Cloud Studio

SDK中心

API中心

命令行工具

涵盖代码开发、场景应用、自动测试全流程，助你从零构建专属AI助手

一站式MCP教程库，解锁AI应用新玩法

聚焦“写作效率、视觉美观与运行性能”三方面进行全面升级，为您提供更高效、稳定的创作环境

社区富文本&Markdown编辑器全新改版上线，欢迎大家体验!

诚挚邀请您参与本次调研，分享您的真实使用感受与建议。您的反馈至关重要，感谢您的支持与参与！

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

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

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐