第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

吸附与催化是当前热门研究领域,常见的计算结果包括催化剂结构,磁性,电子性质等,这些性质对催化剂有着决定性作用。本章华算科技朱老师将正式介绍催化领域的表面与低维结构模型构建,磁性体系处理,关键电子性质计算。
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

引言

第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
反铁磁(Antiferromagnetic, AFM)结构是磁性材料中最基本和最常见的磁有序状态之一。在反铁磁材料中,相邻原子的磁矩方向相反,相互抵消,因此宏观上不显示净磁矩。
反铁磁材料在自旋电子学、信息存储、磁传感器等领域具有重要应用价值。使用VASP(Vienna Ab initio Simulation Package)进行反铁磁结构计算是凝聚态物理和材料科学研究的重要手段。本文系统介绍VASP中反铁磁结构计算的理论基础、关键参数设置、计算流程和常见问题。
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

反铁磁结构的基本原理

第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

反铁磁与铁磁的区别

铁磁(Ferromagnetic, FM:所有原子磁矩平行排列,体系具有净磁矩,宏观上表现出自发磁性。
反铁磁(Antiferromagnetic, AFM:相邻原子磁矩反向平行排列,正负磁矩相互抵消,净磁矩为零。典型的反铁磁结构包括:
FM:  ↑ ↑ ↑ ↑ ↑ ↑   (所有磁矩同向) 
AFM: ↑ ↓ ↑ ↓ ↑ ↓   (相邻磁矩反向) 
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

磁交换作用能

反铁磁与铁磁的能量差反映了材料的磁交换作用强度。交换作用能可通过下式计算:
ΔE=EAFM-EFM
ΔE,则反铁磁态更稳定;若ΔE>0,则铁磁态更稳定。
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

INCAR关键参数设置

第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

 ISPIN参数

ISPIN是控制自旋极化的核心参数,反铁磁计算必须设置:
ISPIN = 2  # 自旋极化计算(反铁磁/铁磁必需) 
说明
默认值为1(非自旋极化),仅适用于无磁体系
计算磁性体系(铁磁、反铁磁)时必须设置为2
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

MAGMOM参数

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 μB
单位:磁矩单位为波尔磁矩(μB)
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

其他重要参数

对称性处理
ISYM = 0    # 关闭对称性(反铁磁计算推荐) 
ISYM = 1    # 使用晶体对称性(可能破坏反铁磁有序) 
建议:反铁磁计算通常使用ISYM=0,因为反铁磁序会破坏晶体的平移对称性,使用ISYM=1可能导致自旋极化设置被自动对称化而失效。
自旋轨道耦合
LNONCOLLINEAR = .TRUE.  # 开启非共线磁性计算 
LSORBIT = .TRUE.        # 开启自旋轨道耦合 
SAXIS = 0 0 1           # 自旋量化轴方向 
对于某些反铁磁材料,特别是重元素体系,自旋轨道耦合可能影响磁各向异性。
结构优化参数
ISTART = 0          # 从头开始计算 
PREC = Accurate     # 高精度 
ENCUT = 450-520     # 截断能 
EDIFF = 1E-5        # 电子步收敛标准 
EDIFFG = -0.02      # 离子步收敛标准(力) 
ISIF = 2/3          # 优化原子位置和晶胞 
IBRION = 2          # 准牛顿法优化 
NSW = 100-300       # 离子步最大数 
POTIM = 0.1-0.5     # 时间步长 
ISMEAR = 0          # 四面体方法(绝缘体) 
SIGMA = 0.05        # 展宽宽度 
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

反铁磁计算完整流程

第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

准备工作

构建反铁磁超胞
对于反铁磁计算,通常需要构造超胞(supercell)以容纳反平行磁矩排列:
# 示例:构建2×1×1超胞用于反铁磁计算 
# 原胞:1个磁性原子 
# 超胞:2个磁性原子,可实现AFM排列 
超胞构建原则
磁矩反向的原子必须位于超胞内,不能仅靠周期性边界实现
超胞大小应确保反铁磁相邻原子在超胞内成对出现
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

结构优化流程

第一步:非磁结构优化
# INCAR设置 
SYSTEM = Structure Optimization (Non-magnetic) 
PREC = Accurate 
ENCUT = 520 
EDIFF = 1E-5 
EDIFFG = -0.02 
ISPIN = 1         # 非自旋极化 
ISYM = 1 
ISMEAR = 0 
SIGMA = 0.05 
IBRION = 2 
NSW = 200 
ISIF = 3 
POTIM = 0.5 
ALGO = Normal 
目的:获得合理的几何结构,避免因初始磁性导致的结构畸变。
第二步:反铁磁结构优化
# INCAR设置 
SYSTEM = AFM Structure Optimization 
PREC = Accurate 
ENCUT = 520 
EDIFF = 1E-5 
EDIFFG = -0.02 
ISPIN = 2         # 自旋极化 
ISYM = 0          # 关闭对称性 
ISMEAR = 0 
SIGMA = 0.05 
IBRION = 2 
NSW = 200 
ISIF = 3 
POTIM = 0.5 
ALGO = Normal 
# 反铁磁磁矩设置 
MAGMOM = 4*2.0 -4*-2.0  # 根据原子数调整 
# 读取第一步的CHGCAR 
ICHARG = 11 
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

静态计算流程

结构优化完成后,进行静态能量计算:
# INCAR设置 
SYSTEM = AFM Static Calculation 
PREC = Accurate 
ENCUT = 520 
EDIFF = 1E-6 
ISPIN = 2 
ISYM = 0 
ISMEAR = 0 
SIGMA = 0.01 
ICHARG = 11       # 读取CHGCAR 
IBRION = -1       # 不进行离子步 
NSW = 0 
LWAVE = .TRUE.    # 保存WAVECAR 
LCHARG = .TRUE.   # 保存CHGCAR 
LORBIT = 11       # 输出投影态密度 
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

本章要点总结

第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

VASP反铁磁结构计算知识点梳理

反铁磁结构的基本原理:介绍了反铁磁与铁磁的区别磁交换作用能
INCAR关键参数设置ISPIN参数MAGMOM参数其他重要参数
反铁磁计算完整流程构建反铁磁超胞结构优化流程静态计算流程
第三章:VASP吸附与催化教程(2026新版):反铁磁结构计算全解析!

下一步学习建议

下一章将正式引入本次教程的核心—VASP亚铁磁结构计算。我们将从亚铁磁结构的基本原理INCAR关键参数设置亚铁磁计算完整流程方面详细介绍

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
(0)
上一篇 2026年5月29日 下午2:47
下一篇 2026年6月1日 上午11:53

相关推荐