第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

亚铁磁(Ferrimagnetic, FiM)是磁性材料中一种特殊的磁有序状态,介于铁磁反铁磁之间。在亚铁磁材料中,不同子晶格的原子磁矩方向相反,但大小不相等,因此宏观上仍具有净磁矩。
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

引言

第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
亚铁磁材料在自旋电子学、微波器件、磁记录等领域具有重要应用价值。典型的亚铁磁材料包括铁氧体(如Fe₃O₄)、钇铁石榴石(YIG)等。本文系统介绍VASP计算亚铁磁结构的方法、关键参数设置和常见问题。
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

亚铁磁结构的基本原理

第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

亚铁磁与铁磁、反铁磁的区别

铁磁(Ferromagnetic, FM:所有原子磁矩平行同向排列,净磁矩最大
反铁磁(Antiferromagnetic, AFM:相邻原子磁矩反向平行,大小相等,净磁矩为零
亚铁磁(Ferrimagnetic, FiM:相邻原子磁矩反向排列,但大小不等,净磁矩不为零
FM:  ↑ ↑ ↑ ↑ ↑ ↑      (所有磁矩同向,净磁矩最大) 
AFM: ↑ ↓ ↑ ↓ ↑ ↓      (磁矩大小相等,方向相反,净磁矩为零) 
FiM: ↑↓↑ ↑↓↑ ↑↓↑      (磁矩大小不等,方向相反,净磁矩不为零) 
亚铁磁的典型例子是磁铁矿(Fe₃O₄,其中:
B位Fe³⁺离子磁矩:+5 μB(向下)
A位Fe³⁺离子磁矩:+5 μB(向上)
B位Fe²⁺离子磁矩:+4 μB(向下)
净磁矩:约4 μB/分子式
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

亚铁磁的微观机制

亚铁磁源于不同原子位点的磁矩大小差异。常见的产生机制包括:
不同价态离子:如Fe²⁺和Fe³⁺磁矩不同
不同元素:如Fe-Co、Fe-Ni合金中不同元素磁矩不同
晶体场效应:不同晶格位置的d轨道分裂不同
亚铁磁材料通常由两个或多个磁性子晶格组成,各子晶格磁矩反向排列但大小不等。
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

INCAR关键参数设置

第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

ISPIN参数

ISPIN是自旋极化的核心控制参数:
ISPIN = 2  # 自旋极化计算(铁磁/反铁磁/亚铁磁必需) 
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

MAGMOM参数

MAGMOM是亚铁磁计算的核心参数,需按原子顺序指定每个原子的初始磁矩:
# 示例:Fe3O4亚铁磁结构 
# A位Fe³⁺(向上),B位Fe³⁺(向下),B位Fe²⁺(向下) 
MAGMOM = 8*5.0 -16*-5.0 -8*-4.0  # Fe3O4单位晶胞 
# 通用格式 
MAGMOM = m1 -m2 -m3 …  # 第一个子晶格向上,其余向下 
MAGMOM = m1 -m2           # 双子晶格亚铁磁 
# 使用重复符 
MAGMOM = 8*5.0 -24*-5.0   # 8个原子磁矩为5,24个为-5 
设置要点
原子顺序严格对应:MAGMOM顺序与POSCAR原子顺序一致
子晶格区分:不同子晶格磁矩方向相反,但大小不等
磁矩符号:正值自旋向上,负值自旋向下
初始值估计:根据实验磁矩或原子磁矩理论值设置
单位:波尔磁矩(μB)
Fe₃O₄示例
# Fe3O4磁性离子磁矩 
# Fe³⁺(3d⁵): ~5 μB 
# Fe²⁺(3d⁶): ~4 μB 
# 在Fe3O4反尖晶石结构中 
# A位:8个Fe³⁺,向上 
# B位:16个Fe³⁺ + 8个Fe²⁺,向下 
MAGMOM = 8*5.0 -16*-5.0 -8*-4.0 
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

NUPDOWN参数

NUPDOWN用于控制自旋电子数差,对亚铁磁计算有一定帮助:
NUPDOWN = 差值  # 自旋向上与向下电子数之差 
说明
NUPDOWN设置总磁矩,MAGMOM设置各原子磁矩
两者应保持一致以保证收敛
VASP文档建议:MAGMOM总和应与NUPDOWN一致
示例
# Fe3O4每个分子式约4 μB净磁矩 
# 1个Fe3O4含28个电子 
NUPDOWN = 4  # 总磁矩约4 μB 
# 但MAGMOM需要分别指定各原子 
MAGMOM = 8*5.0 -16*-5.0 -8*-4.0 
# 总磁矩:8×5 + 16×(-5) + 8×(-4) = 40 – 80 – 32 = -72 μB 
# 需要调整比例 
MAGMOM = 8*5.0 -24*-5.0  # 净磁矩:40 – 120 = -80 μB 
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

其他重要参数

对称性处理
ISYM = 0    # 关闭对称性(亚铁磁计算推荐) 
ISYM = 1    # 使用对称性(可能被磁性破坏) 
建议:亚铁磁体系因磁矩差异会破坏晶体对称性,建议设置ISYM=0。
GGA+U参数
对于过渡金属氧化物等强关联体系:
ISPIN = 2 
LDAU = .TRUE. 
LDAUTYPE = 2 
LDAUL = 2 2 -1      # Fe的d轨道 
LDAUU = 4.0 4.0 0.0 # U值,可区分不同位点 
LDAUJ = 0.0 0.0 0.0 # J值 
MAGMOM = … 
混合参数
AMIX = 0.2      # 电荷混合 
BMIX = 0.0001   # 电荷二次混合 
AMIX_MAG = 0.4  # 磁矩混合 
BMIX_MAG = 0.001 # 磁矩二次混合 
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

亚铁磁计算完整流程

第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

准备工作

构建POSCAR文件
以Fe₃O₄反尖晶石结构为例:
Fe3O4 Ferrimagnetic 
1.0
8.389600  0.000000  0.000000 
0.000000  8.389600  0.000000 
0.000000  0.000000  8.389600 
Fe O 
12 16 
Direct
# A位Fe³⁺ (8个) 
0.25000000 0.25000000 0.25000000 
0.75000000 0.75000000 0.25000000 
# … 8个A位Fe 
# B位Fe³⁺ (16个) 
# … 16个B位Fe 
# O离子 (16个) 
# … 16个O 
要点
确保原子坐标正确
磁性原子数量与MAGMOM对应
亚铁磁结构的超胞需包含完整的反平行磁矩对
准备POTCAR
# 下载PBE泛函的PAW赝势 
# 对于Fe和O元素,选择对应的POTCAR 
# 将POTCAR合并到工作目录 
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

结构优化流程

第一步:非磁结构优化
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 
第二步:亚铁磁结构优化
SYSTEM = FiM 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 = 8*5.0 -16*-5.0 -8*-4.0 
# 可选:总磁矩约束 
NUPDOWN = -128    # 总磁矩:40 – 80 – 32 = -72 μB 
# 读取CHGCAR 
ICHARG = 11 
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

静态计算流程

SYSTEM = FiM Static Calculation 
PREC = Accurate 
ENCUT = 520 
EDIFF = 1E-6 
ISPIN = 2 
ISYM = 0 
ISMEAR = 0 
SIGMA = 0.01 
ICHARG = 11 
IBRION = -1 
NSW = 0 
LWAVE = .TRUE. 
LCHARG = .TRUE. 
LORBIT = 11       # 投影态密度 
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

能带计算流程

SYSTEM = FiM Band Structure 
PREC = Accurate 
ENCUT = 520 
EDIFF = 1E-6 
ISPIN = 2 
ISYM = 0 
ICHARG = 11 
LWAVE = .FALSE. 
LORBIT = 10 
IBRION = -1 
NSW = 0 
MAGMOM = 8*5.0 -16*-5.0 -8*-4.0 
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

本章要点总结

第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

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

亚铁磁结构的基本原理:介绍了亚铁磁与铁磁、反铁磁的区别亚铁磁的微观机制
INCAR关键参数设置ISPIN参数MAGMOM参数其他重要参数
亚铁磁计算完整流程准备工作结构优化流程能带计算流程
第四章:VASP吸附与催化教程(2026新版):亚铁磁结构计算讲解!

下一步学习建议

下一章将正式引入本次教程的核心—结构优化与原子磁矩分析。我们将从结构优化基本原理INCAR关键参数设置原子磁矩分析方法方面详细介绍

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

相关推荐