五边形MnX2单层中交变磁性与铁电性的共存与耦合Phys. Rev. B 113, 144405 (2026)五边形MnX2单层中交变磁性与铁电性的共存与耦合Interplay of Altermagnetism and Ferroelectricity in Pentagonal MnX2 Monolayers导读 导读交变磁性Altermagnetism是近年发现的新型磁性相兼具反铁磁的零净磁矩和铁磁的自旋劈裂能带但电场调控AM序一直是挑战。本文通过DFT计算和对称性分析首次在五边形MnX2XS, Se单层中揭示了g波交变磁性与II型铁电性的共存和耦合。配体二聚体超交换机制赋予室温Neel温度MnS2: 370K, MnSe2: 285KII型铁电性将Neel矢量锁定在面外极化方向实现电场稳定AM序。Janus MnSSe进一步增强了AM劈裂和铁电性。【一、前言背景】交变磁性Altermagnetism磁学领域的第三类磁性交变磁性Altermagnetism, AM是近年来发现的一种新型磁性相不同于传统的铁磁FM和反铁磁AFM。AM材料具有补偿的磁矩净磁化强度为零类似AFM但同时在动量空间中具有自旋劈裂的能带类似FM。这种无净磁矩的自旋劈裂是AM的核心特征源于磁空间群中连接相反自旋子晶格的对称操作不是平移或反演而是旋转因此称为alter-magnetism。AM材料兼具FM和AFM的优势无净磁矩抗外磁场干扰高密度存储和自旋劈裂能带可实现反常Hall效应、自旋流生成。AM已被预测存在于RuO2、MnTe、V2Se2O、CrSb等多种材料中但电场调控AM序仍然是一个挑战。本文的核心创新首次在五边形MnX2XS, Se单层中揭示了交变磁性与铁电性的共存和耦合。II型铁电性由d-p轨道杂化产生与Neel矢量的方向耦合使电场可以调控AM序为AM材料的电控提供了新途径。五边形二维材料从石墨烯到MnX2的新范式核心问题如何在二维极限下实现稳健的AM序和铁电性共存五边形结构Cairo tessellation提供了一种独特的原子构型。方法体系VASPDFT, GGAU, Ueff3.0 eV for Mn- PHONOPYDFPT动力学稳定性- Berry Phase铁电极化- MCSOLVERMC模拟TN- VAMPIRELLG自旋动力学模拟。关键发现(1) MnS2和MnSe2均为g波AM室温Neel温度TN~370K和285K(2) 通过配体二聚体介导的超交换机制稳定AM序(3) II型铁电性将Neel矢量锁定在面外极化方向实现电场稳定AM序(4) Janus MnSSe具有增强的AM劈裂和铁电性。五边形MnX2单层中交变磁性与铁电性的共存与耦合研究流程。晶体结构p-421m, 五边形- VASP DFTGGAU, Ueff3.0eV- 交换参数提取J1, J2, J3- Phonopy动力学稳定性- Berry Phase铁电极化- g波交变磁性补偿磁矩k空间自旋劈裂- 配体二聚体超交换 - II型铁电性Mn-d/S-p杂化- MC模拟TN: MnS2370K, MnSe2285K- LLG自旋动力学电场稳定AM序- 应变Janus工程。【二、研究方法】DFT计算设置VASP GGAU Phonopy Berry PhaseVASP计算PAW赝势GGA-PBE泛函ENCUT600 eVGamma-centered 19x19x1 Monkhorst-Pack k点网格。Mn 3d轨道的在位库仑相互作用Ueff3.0 eV基于文献中同类Mn化合物的经验值。Phonopy使用DFPT方法计算声子谱确认五边形MnX2单层的动力学稳定性。对于二维材料声子谱中无虚频是结构可合成的必要条件。Berry Phase方法用于计算铁电极化。对于二维材料单位体积使用原子层厚度去除真空层厚度。注意二维材料的铁电极化计算需要特别注意偶极修正和真空层厚度的影响。多尺度模拟从DFT到MC到LLG交换参数提取从DFT总能量计算中提取最近邻J1、次近邻J2和第三近邻J3交换参数。通过比较FM和不同AFM构型的总能量用Heisenberg模型拟合得到J值。MC模拟使用MCSOLVER在50x50超胞上进行Monte Carlo模拟确定Neel温度TN。MC模拟使用Metropolis算法在有限温度下采样自旋构型空间。LLG自旋动力学使用VAMPIRE求解Landau-Lifshitz-Gilbert方程验证电场对AM序的稳定作用。在LLG方程的Hamiltonian中加入-E*P项电场-极化耦合模拟电场对磁序的调控。Heisenberg模型用于描述MnX2中的磁相互作用。J1, J2, J3分别为最近邻、次近邻和第三近邻交换参数K为磁各向异性能MAE常数。LLG方程描述自旋在有效场Heff作用下的进动和阻尼。在VAMPIRE中Heff包括交换场、各向异性场和外加电场-极化耦合项。Berry Phase方法铁电极化的量子力学计算公式。对于二维材料需要对Berry联络在BZ中积分。【三、核心结果】图 1五边形MnS2和MnSe2单层的晶体结构。(a-b) MnS2和MnSe2的五边形晶格结构p-421m层群Mn和S/Se原子分别占据2b和4e Wyckoff位置。五边形结构由交替的Mn和S/Se二聚体组成形成Cairo tessellation图案。图 2MnX2的自旋劈裂能带结构和Brillouin区分析。(a-b) MnS2和MnSe2的自旋分辨能带结构展示g波AM的k空间自旋劈裂特征。(c) Brillouin区中的自旋劈裂分布展示AM的d波/g波对称性。g波交变磁性自旋劈裂的对称性分析MnX2中的AM属于g波类型l4的角动量对称性这意味着自旋劈裂在k空间中的角分布具有四重旋转对称性加上符号翻转。g波AM的特征是自旋劈裂在Gamma-M方向最大在Gamma-X方向为零节点。自旋劈裂的物理起源Mn的3d轨道通过配体S/Se二聚体介导的超交换相互作用产生AFM耦合。但由于Mn子晶格之间的旋转对称性连接而非平移或反演自旋向上和向下能带在k空间中发生劈裂。室温TN是AM材料实用化的关键MnS2的TN370K和MnSe2的TN285K均高于室温表明这些材料在实际器件工作温度下可保持AM序。这一高TN源于强的Mn-S/Se-Mn超交换耦合。图 3Neel矢量方向与面外极化强度的关系。(a) Neel矢量方向的定义。(b) 极化强度随Neel矢量方向的正弦变化展示II型铁电性的特征耦合。(c) 电场对Neel矢量方向的稳定作用验证电控AM序的可行性。II型铁电性与AM序的电控机制II型铁电性由d-p轨道杂化产生与I型铁电性由结构畸变产生有本质区别。在MnX2中Mn-d和S/Se-p轨道的杂化导致非对称的电荷分布产生面外方向的铁电极化。关键耦合机制铁电极化强度P依赖于Neel矢量L的方向遵循正弦关系P(L) ~ sin(theta)。这意味着通过翻转Neel矢量如通过外加磁场可以改变铁电极化方向反之通过外加电场改变极化方向可以稳定Neel矢量在特定方向。LLG模拟验证在LLG方程中加入-E*P项后外加电场可以有效地将Neel矢量锁定在面外方向即使在有限温度下也能防止热涨落导致的磁序翻转。这为AM材料的电控提供了可行的理论方案。图 4压缩应变对MnX2自旋劈裂的影响。(a) 不同应变下的自旋劈裂能带。(b) 自旋劈裂值和交换参数随应变的变化。(c) 应变对TN的影响。压缩应变增强磁性耦合提高AM序的稳定性。应变工程增强AM耦合和自旋劈裂压缩应变compressive strain显著增强了MnX2中的磁交换耦合。物理机制压缩应变减小了Mn-Mn原子间距增强了通过配体二聚体的超交换相互作用从而增大J1和J2值。自旋劈裂随压缩应变增大在-4%双轴压缩应变下MnS2和MnSe2的自旋劈裂值分别增加约30%和40%。这为通过衬底选择晶格失配应变来调控AM性能提供了实验可行的途径。TN也随压缩应变升高在-4%应变下MnS2的TN从370K升至约420K进一步确保了AM序在器件工作温度下的稳定性。图 5Janus MnSSe的晶体结构和自旋劈裂。(a) Janus MnSSe的晶体结构一侧为S原子另一侧为Se原子打破面外对称性。(b) 自旋劈裂能带结构展示比MnS2和MnSe2更大的自旋劈裂。(c) 增强的铁电性和AM耦合。Janus MnSSe对称性破缺增强AM和铁电性Janus结构一面S、一面Se打破了面外镜面对称性产生两个效果(1) 增强面外铁电极化两面不同的电负性产生内建电场(2) 增强AM自旋劈裂不对称的配体环境改变了d-p杂化。Janus MnSSe的自旋劈裂值显著大于MnS2和MnSe2表明对称性破缺是增强AM性能的有效策略。这为设计更高性能的AM材料提供了新思路通过打破特定的晶体对称性可以增强自旋劈裂而不改变磁基态。实验可行性Janus TMDs如MoSSe已通过CVD方法成功合成。MnSSe的合成可能通过类似的硒化/硫化策略实现。【DFT Tips】【DFT Tip 1】GGAU计算中Ueff值的选取Mn体系的经验本文对Mn 3d轨道使用Ueff3.0 eV。Mn的U值选择范围通常在2-5 eV取决于氧化态Mn2 vs Mn3 vs Mn4、配位环境和泛函。U值对Mn体系的影响(1) U增大 - d带更局域 - 带隙增大 - 磁矩增大(2) U过大5 eV可能导致虚假的Mott绝缘态(3) U过小2 eV可能无法正确描述关联效应。建议(1) 优先使用线性响应法计算U值Cococcioni de Gironcoli, PRB 71, 035105(2) 对U值进行敏感性分析如U2,3,4 eV报告关键结果对U的依赖(3) 对于Mn体系Ueff3-4 eV通常是合理的起点。【DFT Tip 2】二维材料的铁电极化计算Berry Phase方法的注意事项二维材料的铁电极化计算比三维材料更复杂因为(1) 真空层方向没有周期性Berry Phase的积分路径需要特别注意(2) 极化量纲C/m^2需要使用原子层厚度而非超胞c轴长度。VASP设置开启LBERRY.TRUE.设置IGPAR指向极化方向。对于面外极化z方向需要在INCAR中设置正确的真空层厚度作为有效层厚。常见错误(1) 使用超胞c轴长度而非原子层厚度计算极化导致极化值被严重低估(2) 未进行偶极修正LDIPOL.TRUE.导致极化值受真空层厚度影响。【DFT Tip 3】交换参数J的提取从DFT总能量到Heisenberg模型从DFT提取交换参数的标准方法(1) 构建多种磁构型FM, AFM1, AFM2, AFM3...(2) 计算每种构型的DFT总能量(3) 用Heisenberg模型E E0 - sum(J_ij * S_i * S_j)进行最小二乘拟合。对于MnX2每个Mn有S5/2需要至少4种磁构型来唯一确定J1, J2, J3。建议使用更多构型6-8种进行过定拟合评估拟合误差。常见陷阱(1) 磁构型不够多导致J值不唯一(2) 未考虑周期性边界条件导致的磁相互作用计数错误(3) 对于金属体系Heisenberg模型可能不适用需要RKKY模型。【DFT Tip 4】Monte Carlo模拟中的超胞尺寸选择本文使用50x50超胞2500个自旋进行MC模拟。超胞尺寸的选择需要平衡统计精度和计算成本。一般原则(1) 对于二维体系超胞尺寸至少为40x40以确保有限尺寸效应可忽略(2) 对于三维体系至少20x20x20(3) MC步数至少10^5用于平衡10^5-10^6用于采样。验证方法使用不同超胞尺寸如30x30, 50x50, 80x80计算TN确认TN不再随超胞尺寸变化。如果TN随尺寸变化5%需要增大超胞。注意MC模拟的TN通常略高于实验值因为MC忽略了量子涨落和自旋-轨道耦合的细节效应。【DFT Tip 5】交变磁性AM的DFT识别如何判断一个材料是AMAM的DFT识别需要三个步骤(1) 确认磁基态为补偿AFM净磁矩0(2) 计算自旋分辨能带确认k空间中存在自旋劈裂不同k点自旋向上和向下能带不简并(3) 确认连接相反自旋子晶格的对称操作是旋转而非平移或反演。VASP中的自旋分辨能带在非共线计算LNONCOLLINEAR.TRUE.中设置LORBIT11可以输出每个能带的自旋投影Sx, Sy, Sz。对于AMSx, Sy, Sz在k空间中的分布应展示特定的对称性如d波、g波。常见错误将SOC诱导的自旋劈裂Rashba/Dresselhaus误认为AM自旋劈裂。AM的自旋劈裂远大于SOC劈裂通常100 meV vs 10 meV且对称性不同。【DFT Tip 6】LLG模拟中的温度效应和电场耦合LLG方程描述的是零温自旋动力学。要模拟有限温度效应需要在LLG方程中加入随机热场stochastic LLG将温度以Langevin热噪声的形式引入。VAMPIRE中的设置temperature参数控制热噪声的幅度。电场耦合通过添加-E*P能量项实现其中P可以从DFT的Berry Phase计算获得。常见陷阱(1) 热噪声的幅度与温度和时间步长有关步长过大导致数值不稳定(2) 电场-极化耦合的强度需要从DFT计算中提取不能随意设定(3) LLG模拟的阻尼系数alpha通常取0.01-0.1对结果有显著影响。【DFT Tip 7】二维材料声子谱计算虚频的判断与处理对于二维材料声子谱中Gamma点附近可能出现小的虚频 -10 cm^-1这通常是由于数值误差而非真正的结构不稳定性。判断标准(1) 虚频在Gamma点附近且|omega| 10-20 cm^-1可能是数值误差(2) 虚频远离Gamma点或|omega| 50 cm^-1通常说明结构不稳定(3) 增大超胞尺寸可以消除Gamma点附近的数值虚频。解决方案(1) 增大超胞尺寸至少3x3x1(2) 提高力收敛标准 0.001 eV/A(3) 对于极性二维材料开启LO-TO劈裂修正LPHON_POLAR.TRUE.。【DFT Tip 8】Janus材料的DFT计算对称性破缺的处理Janus MnSSe打破了面外镜面对称性这需要在DFT计算中正确处理(1) 空间群从p-421m的层群降低到p-4无反演、无镜面(2) 偶极修正LDIPOL.TRUE., IDIPOL3是必需的因为不对称表面产生净偶极矩(3) 真空层厚度需要足够大20 A以避免相邻slab之间的偶极-偶极相互作用。常见错误忽略偶极修正导致虚假的层间相互作用表现为能带色散和功函数随真空层厚度变化。验证方法改变真空层厚度如从15A到25A确认能带结构和功函数收敛。【知识扩展】【知识扩展 1】交变磁性Altermagnetism从理论预测到实验验证【理论解释】AM是2019-2022年由Smejkal、Sinova和Jungwirth等人提出的新型磁性分类。在传统的FM/AFM二分法中AM被归类为AFM补偿磁矩但其自旋劈裂的能带结构与传统AFM有本质区别。AM的自旋劈裂源于磁空间群中连接相反自旋子晶格的旋转对称性而非平移或反演。【方法比较】AM vs FM两者都有自旋劈裂能带但AM无净磁矩抗干扰、高密度。AM vs AFM两者都有补偿磁矩但AM有自旋劈裂能带可实现反常Hall效应、自旋流。AM 无净磁矩 自旋劈裂能带。【经典参考】Smejkal et al., PRX 12, 011028 (2022)--AM理论框架Smejkal et al., PRX 12, 040501 (2022)--AM中的反常Hall效应Krempasky et al., Nature 626, 517 (2024)--MnTe中AM的ARPES实验验证。【迁移能力】AM的电控策略通过II型铁电性耦合可推广到其他AM材料特别是具有d-p杂化和结构不对称性的体系。【知识扩展 2】多铁性磁电耦合的物理机制【理论解释】多铁性材料同时具有铁电序和磁序。磁电耦合ME coupling是两者之间的相互作用允许电场控制磁性或磁场控制电极化。多铁性分为两类(1) I型多铁铁电性和磁性起源不同耦合弱(2) II型多铁铁电性由磁序驱动如自旋螺旋、交换伸缩耦合强。【MnX2中的耦合机制】II型铁电性源于Mn-d和S/Se-p的杂化极化方向与Neel矢量方向耦合。这种耦合是锁定的而非可切换的但提供了通过电场稳定特定磁序的途径。【经典参考】Cheong Mostovoy, Nat. Mater. 6, 13 (2007)--多铁性综述Tokura et al., RPP 77, 076501 (2014)--多铁性中的磁电耦合Spaldin Ramesh, Nat. Mater. 18, 203 (2019)--多铁性进展。【迁移能力】五边形结构中的AM-铁电耦合机制可推广到其他五边形二维材料如PdS2, PdSe2和Janus结构。【科研经验】【科研经验 1】新型磁性材料论文的证据链要求问题AM是新兴领域审稿人通常要求提供充分的证据来确认AM而非传统AFM。原因AM和AFM在磁矩补偿上是相同的区分它们需要更精细的电子结构证据。解决方案本文的策略值得借鉴(1) 自旋分辨能带展示k空间自旋劈裂AM的核心特征(2) 对称性分析说明为什么自旋劈裂是允许的磁空间群中连接子晶格的是旋转对称性(3) 自旋劈裂的角分布g波、d波等确认AM类型(4) 与已知AM材料的比较。建议AM论文的核心不是发现了新的磁性材料而是证明了这种磁性是AM而非传统AFM。如果自旋劈裂的证据不充分审稿人可能认为这只是普通的AFM。【科研经验 2】DFTMCLLG多尺度模拟的整合策略问题DFT、MC和LLG是三个不同尺度的模拟方法如何将它们整合成一个连贯的计算工作流原因DFT电子尺度- 交换参数J - MC统计力学尺度- TN, 磁化强度 - LLG连续介质尺度- 自旋动力学。每个步骤的误差都会传播到下一步。解决方案(1) DFT提取J值时应进行交叉验证不同磁构型、不同泛函(2) MC模拟的TN应与实验如果可用或DFT的磁各向异性分析比较(3) LLG模拟的阻尼参数alpha应从实验或第一性原理提取。建议多尺度模拟论文中应明确报告每个步骤的误差来源和传播。如果DFT的J值误差为10%MC的TN误差可能为15-20%这对LLG模拟的定量预测有显著影响。【如果是我我还会继续算】【继续算 1】Berry曲率与反常Hall效应AHE为什么值得算AM的自旋劈裂能带应产生非零的Berry曲率从而可能产生反常Hall效应即使净磁矩为零。这是AM区别于传统AFM的关键实验特征。能回答的问题MnX2的AHE有多大AHE的角度依赖性与AM的g波对称性是否一致适合体系所有AM材料。输入Wannier90 Berry曲率计算。【继续算 2】自旋Hall效应与自旋-电荷转换为什么值得算AM材料的自旋劈裂能带可能产生巨大的自旋Hall效应这在自旋电子学中具有重要应用价值。能回答的问题MnX2的自旋Hall电导率是多少是否与AM的对称性相关适合体系所有AM材料。输入Wannier90 Kubo公式计算SHC。【继续算 3】磁振子Magnon能带与AM特征为什么值得算AM的磁振子能带应展示手性劈裂不同于传统AFM这是AM的另一个可实验验证的特征。能回答的问题MnX2的磁振子能带是否展示手性劈裂劈裂的幅度和对称性与AM的g波特征是否一致适合体系所有AM材料。输入自旋波理论计算或DFT磁振子计算。【继续算 4】AM的磁各向异性MAE与外场响应为什么值得算AM的磁各向异性决定了磁序的稳定性方向和外场响应。对于二维材料MAE通常很小 1 meV/atom可能不足以抵抗热涨落。能回答的问题MnX2的MAE是多少面内还是面外易轴MAE是否足以在室温下稳定AM序适合体系所有二维磁性材料。输入非共线SOC DFT计算到中等。【继续算 5】其他五边形AM材料的高通量筛选为什么值得算本文证明五边形结构是AM和铁电性共存的优秀平台。是否存在其他五边形过渡金属硫族化合物也具备AM性质能回答的问题五边形CrX2、FeX2、CoX2是否也是AM它们的TN和自旋劈裂如何适合体系所有五边形TMDs。输入高通量DFT计算但可并行。【继续算 6】Janus材料的合成可行性与稳定性为什么值得算Janus MnSSe是理论预测的材料但其合成可行性和热力学稳定性尚未验证。能回答的问题Janus MnSSe是否在热力学上稳定是否存在相分离为MnS2和MnSe2的趋势合成温度窗口是多少适合体系所有Janus二维材料。输入AIMD 声子计算 凸包图分析。Jiang, Cao, Guo, Zeng, Zhang Chen, Sichuan Univ. | Phys. Rev. B 113, 144405 (2026) | 交变磁性 铁电性 MnX2 五边形结构 VASP MC LLG