在催化、吸附体系的理论计算中,不少催化剂、活性位点及低维表面结构天然带有磁性特征,自旋电子效应会直接改变体系的电子分布、界面作用与反应活性,是影响催化性能的关键因素。因此在 VASP 计算中,磁性体系的合理设置、自旋态判定与结果分析,成为表面建模、电子性质计算不可或缺的一环。本章华算科技朱老师将正式介绍催化领域的表面与低维结构模型构建,磁性体系处理,关键电子性质计算。
铁磁(Ferromagnetic, FM)结构是磁性材料中最基本和最常见的磁有序状态之一。在铁磁材料中,所有原子的磁矩在同一方向上平行排列,产生自发磁化,宏观上表现出强磁性。
铁磁材料广泛应用于信息存储、磁传感器、电机、变压器等领域。使用VASP(Vienna Ab initio Simulation Package)进行铁磁结构计算是凝聚态物理和材料科学研究的常规手段。
铁磁(Ferromagnetic, FM):所有原子磁矩平行同向排列,体系具有净磁矩,宏观上表现出自发磁性。反铁磁(Antiferromagnetic, AFM):相邻原子磁矩反向平行排列,正负磁矩相互抵消,净磁矩为零。顺磁(Paramagnetic, PM):原子具有固有磁矩,但在无外场时磁矩随机取向,无净磁矩。FM: ↑ ↑ ↑ ↑ ↑ ↑ (所有磁矩同向,净磁矩最大) AFM: ↑ ↓ ↑ ↓ ↑ ↓ (相邻磁矩反向,净磁矩为零) 根据居里温度(Curie Temperature, Tc)和磁化强度,铁磁材料可分为:铁磁有序来源于交换相互作用(Exchange Interaction),这是量子力学效应导致的相邻原子磁矩之间的耦合。铁磁交换作用能可表示为:其中Jij>0表示铁磁耦合占优,相邻磁矩倾向于平行排列。
ISPIN是控制自旋极化的核心参数,铁磁计算必须设置:ISPIN = 2 # 自旋极化计算(铁磁/反铁磁必需) MAGMOM是设置铁磁结构的核心参数,用于指定每个原子的初始磁矩:MAGMOM = 8*2.0 # 8个原子,每个初始磁矩为+2.0μB MAGMOM = N*1.0 # N个原子,每个初始磁矩为+1.0μB 原子顺序:MAGMOM的数值顺序与POSCAR中的原子顺序严格对应初始值:建议设置为实验磁矩的1.2-1.5倍,如Fe通常设为2.0-3.0 μBISYM = 1 # 使用晶体对称性(铁磁计算推荐) ISYM = 0 # 关闭对称性(当对称性被磁性破坏时) 建议:铁磁计算通常使用ISYM=1,因为铁磁序不会破坏晶体的平移对称性,使用对称性可加快计算。但需注意,如果初始磁矩设置导致对称性降低,ISYM=1可能被自动调整为ISYM=0。EDIFFG = -0.02 # 离子步收敛标准(力) LREAL = Auto # 自动选择实空间投影(大体系推荐) LREAL = .TRUE. # 显式开启实空间投影 LREAL = .FALSE. # 使用倒易空间投影 BMIX_MAG = 0.001 # 磁矩二次混合参数
2.866500 0.000000 2.866500 0.000000 2.866500 2.866500 2.866500 2.866500 0.000000 0.0000000000 0.0000000000 0.0000000000 0.5000000000 0.5000000000 0.5000000000 确保晶格常数合理(如Fe的bcc结构晶格常数约为2.86 Å)SYSTEM = Structure Optimization (Non-magnetic) 目的:获得合理的几何结构,避免因初始磁性导致的结构畸变。SYSTEM = FM Structure Optimization MAGMOM = 2*2.5 # 2个Fe原子,初始磁矩各2.5μB SYSTEM = FM Static Calculation
铁磁结构的基本原理:介绍了铁磁与反铁磁的区别、铁磁材料的分类、磁交换相互作用INCAR关键参数设置:ISPIN参数、MAGMOM参数、其他重要参数铁磁计算完整流程:构建POSCAR文件、结构优化流程、静态计算流程
下一章将正式引入本次教程的核心—VASP反铁磁结构计算。我们将从反铁磁结构的基本原理、INCAR关键参数设置、反铁磁计算完整流程方面详细介绍VASP反铁磁结构计算。
声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!