VASP（Vienna Ab initio Simulation Package）是凝聚态物理和材料计算领域最常用的第一性原理计算软件之一。结构优化与原子磁矩分析是材料计算研究中最基础且最重要的两个环节。
结构优化能够找到体系能量最低的稳定构型，为后续性质计算提供可靠的几何基础；原子磁矩分析则是理解材料磁性质的关键。本文华算科技朱老师将系统介绍VASP结构优化的完整流程、关键参数设置、原子磁矩提取方法及常见问题解决方案。
结构优化的目标是通过调整原子位置和晶胞参数，使体系总能量最小化，找到势能面上的能量极小值点。在密度泛函理论（DFT）框架下，结构优化涉及电子步（SCF迭代）和离子步的交替优化。
VASP提供了多种离子步优化算法，各有特点：
|
方法 |
IBRION参数 |
算法 |
适用场景 |
|
最速下降法 |
0 |
RMM-DIIS |
初始结构较远时 |
|
共轭梯度法 |
1 |
CG |
大体系、内存有限 |
|
准牛顿法 |
2 |
标准优化 |
常规结构优化 |
|
自适应方法 |
3 |
自动步长 |
复杂势能面 |
准牛顿法（IBRION=2）是最常用的方法，它通过Hessian矩阵更新来加速收敛，平衡了速度与精度，适用于大多数体系。
ISIF参数控制优化内容：
|
ISIF值 |
优化内容 |
|
2 |
仅优化原子位置，晶胞固定 |
|
3 |
优化原子位置+晶胞形状+晶胞体积 |
|
4 |
优化原子位置+晶胞形状，固定体积 |
|
6 |
优化原子位置+晶胞体积，固定形状 |
|
7 |
优化晶胞形状+体积，原子位置固定 |
常规结构优化推荐使用ISIF=3，可全面优化几何结构。
# 结构优化基本设置
IBRION = 2 # 离子步优化算法（2=准牛顿法）
NSW = 200 # 最大离子步数
EDIFFG = -0.02 # 力收敛标准（负值=力收敛）
ISIF = 3 # 优化类型（3=位置+体积+形状）
POTIM = 0.5 # 步长参数
PREC = Accurate # 精度设置
ENCUT = 500 # 平面波截断能（需测试收敛）
PREC = Accurate # 高精度计算
EDIFF = 1E-6 # 电子步收敛标准
EDIFFG = -0.02 # 离子步收敛标准（负值表示力收敛）
ENCUT = 520 # 截断能（建议测试收敛）
ISMEAR = 0 # 展宽方法（0=四面体，1/2=Gaussian）
SIGMA = 0.05 # 展宽宽度
ALGO = Normal # 电子优化算法
ISPIN = 2 # 开启自旋极化计算
MAGMOM = N*2.5 # 初始磁矩设置
ISYM = 1 # 对称性处理（磁性体系建议ISYM=1或0）
LORBIT = 11 # 投影态密度和磁矩输出
LWAVE = .TRUE. # 保存波函数
LCHARG = .TRUE. # 保存电荷密度
POSCAR文件构建
Fe BCC Ferromagnetic Optimization
1.0
2.8665 0.0000 2.8665
0.0000 2.8665 2.8665
2.8665 2.8665 0.0000
Fe
2
Direct
0.0000000000 0.0000000000 0.0000000000
0.5000000000 0.5000000000 0.5000000000
POSCAR构建要点：
1.确保原子位置合理，避免原子重叠距离小于1.5 Å
2.晶格常数接近实验值（Fe的bcc结构约为2.86 Å）
3.磁性原子位置与MAGMOM设置一致
4.使用Direct坐标便于后续观察和修改
POTCAR准备
# 下载PBE泛函的PAW赝势
# 合并多个元素的POTCAR
cat Fe/POTCAR O/POTCAR > POTCAR
第一步：非磁结构优化
# 非磁结构优化INCAR
SYSTEM = Non-magnetic Structure Optimization
PREC = Accurate
ENCUT = 520
EDIFF = 1E-6
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 = Magnetic Structure Optimization
PREC = Accurate
ENCUT = 520
EDIFF = 1E-6
EDIFFG = -0.02
ISPIN = 2 # 开启自旋极化
ISYM = 1
ISMEAR = 0
SIGMA = 0.05
IBRION = 2
NSW = 200
ISIF = 3
POTIM = 0.5
ALGO = Normal
MAGMOM = 2*2.5 # 初始磁矩设置
ICHARG = 11 # 读取CHGCAR加速收敛
LWAVE = .TRUE.
LCHARG = .TRUE.
结构优化完成后，进行静态计算获取精确磁矩：
# 静态计算INCAR
SYSTEM = Static Calculation After Optimization
PREC = Accurate
ENCUT = 520
EDIFF = 1E-6
ISPIN = 2
ISYM = 1
ISMEAR = 0
SIGMA = 0.01
ICHARG = 11 # 读取优化后的CHGCAR
IBRION = -1 # 不进行离子步
NSW = 0
LWAVE = .TRUE.
LCHARG = .TRUE.
LORBIT = 11 # 输出投影磁矩
结构优化基本原理：介绍了结构优化基本原理、INCAR关键参数设置
INCAR关键参数设置：核心优化参数、精度与收敛参数、自旋极化参数（磁性体系）
原子磁矩分析方法：磁矩数据来源、OUTCAR磁矩提取、OSZICAR磁矩分析
下一章将正式引入本次教程的核心—。我们将从VASP功函数计算功函数的物理定义、计算执行步骤、后处理与数据提取方面详细介绍VASP功函数计算。
