吸附与催化是当前热门研究领域,常见的计算结果包括催化剂结构,磁性,电子性质等,这些性质对催化剂有着决定性作用。本章华算科技朱老师将正式介绍催化领域的表面与低维结构模型构建,磁性体系处理,关键电子性质计算。
反铁磁(Antiferromagnetic, AFM)结构是磁性材料中最基本和最常见的磁有序状态之一。在反铁磁材料中,相邻原子的磁矩方向相反,相互抵消,因此宏观上不显示净磁矩。
反铁磁材料在自旋电子学、信息存储、磁传感器等领域具有重要应用价值。使用VASP(Vienna Ab initio Simulation Package)进行反铁磁结构计算是凝聚态物理和材料科学研究的重要手段。本文系统介绍VASP中反铁磁结构计算的理论基础、关键参数设置、计算流程和常见问题。
铁磁(Ferromagnetic, FM):所有原子磁矩平行排列,体系具有净磁矩,宏观上表现出自发磁性。反铁磁(Antiferromagnetic, AFM):相邻原子磁矩反向平行排列,正负磁矩相互抵消,净磁矩为零。典型的反铁磁结构包括:AFM: ↑ ↓ ↑ ↓ ↑ ↓ (相邻磁矩反向) 反铁磁与铁磁的能量差反映了材料的磁交换作用强度。交换作用能可通过下式计算:若ΔE,则反铁磁态更稳定;若ΔE>0,则铁磁态更稳定。
ISPIN是控制自旋极化的核心参数,反铁磁计算必须设置:ISPIN = 2 # 自旋极化计算(反铁磁/铁磁必需)
MAGMOM是设置反铁磁结构的核心参数,用于指定每个原子的初始磁矩:MAGMOM = 3*2.0 -3*-2.0 # 前3个原子磁矩为+2.0,后3个为-2.0 MAGMOM = N*1.0 -N*-1.0 # N个原子正向,N个原子反向 原子顺序:MAGMOM的数值顺序与POSCAR中的原子顺序严格对应初始值:建议设置为实验磁矩的1.2-1.5倍,如Fe通常设为2.0-3.0 μBISYM = 0 # 关闭对称性(反铁磁计算推荐) ISYM = 1 # 使用晶体对称性(可能破坏反铁磁有序) 建议:反铁磁计算通常使用ISYM=0,因为反铁磁序会破坏晶体的平移对称性,使用ISYM=1可能导致自旋极化设置被自动对称化而失效。LNONCOLLINEAR = .TRUE. # 开启非共线磁性计算 LSORBIT = .TRUE. # 开启自旋轨道耦合 对于某些反铁磁材料,特别是重元素体系,自旋轨道耦合可能影响磁各向异性。EDIFFG = -0.02 # 离子步收敛标准(力)
对于反铁磁计算,通常需要构造超胞(supercell)以容纳反平行磁矩排列:磁矩反向的原子必须位于超胞内,不能仅靠周期性边界实现SYSTEM = Structure Optimization (Non-magnetic) 目的:获得合理的几何结构,避免因初始磁性导致的结构畸变。SYSTEM = AFM Structure Optimization MAGMOM = 4*2.0 -4*-2.0 # 根据原子数调整 SYSTEM = AFM Static Calculation LWAVE = .TRUE. # 保存WAVECAR LCHARG = .TRUE. # 保存CHGCAR
反铁磁结构的基本原理:介绍了反铁磁与铁磁的区别、磁交换作用能、INCAR关键参数设置:ISPIN参数、MAGMOM参数、其他重要参数反铁磁计算完整流程:构建反铁磁超胞、结构优化流程、静态计算流程
下一章将正式引入本次教程的核心—VASP亚铁磁结构计算。我们将从亚铁磁结构的基本原理、INCAR关键参数设置、亚铁磁计算完整流程方面详细介绍。
声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!