吸附与催化是当前热门研究领域，常见的计算结果包括催化剂结构，磁性，电子性质等，这些性质对催化剂有着决定性作用。本章华算科技朱老师将正式介绍催化领域的表面与低维结构模型构建，磁性体系处理，关键电子性质计算。
VASP（Vienna Ab initio Simulation Package）是一种广泛应用的基于密度泛函理论（DFT）的第一性原理计算软件。在材料科学研究中，无磁性结构计算是基础且重要的一环，尤其适用于研究非磁性材料、金属、半导体等体系。
在DFT计算框架中，电子结构计算可分为自旋极化（spin-polarized）和非自旋极化（non-spin-polarized）两种模式。无磁结构计算通常采用非自旋极化模式，即不考虑电子自旋自由度对体系能量的影响。
非自旋极化计算（ISPIN=1）：
不区分自旋向上和自旋向下电子
电子密度对所有自旋方向相同
计算效率较高，适用于非磁性体系
自旋极化计算（ISPIN=2）：
分别处理自旋向上和自旋向下电子
可计算体系的磁矩和磁性性质
适用于磁性体系或需要研究磁性的情况
值得注意的是，VASP中的”无磁”计算实际上对应于”顺磁”（paramagnetic）状态的近似描述。顺磁性指原子或分子具有固有磁矩，但在无外磁场时磁矩随机取向，宏观上不显示磁性。在DFT计算中，ISPIN=1计算不考虑磁矩分布，而实际顺磁状态可能需要通过特殊方法处理，如使用无序局部矩（DLM）方法或超胞平均。
ISPIN是控制自旋极化的核心参数，其取值决定了计算类型：
ISPIN = 1 # 非自旋极化计算（无磁）
ISPIN = 2 # 自旋极化计算（磁）
设置要点：
默认值为1，若不使用MAGMOM参数，通常保持默认即可
对于无磁结构优化和静态计算，必须设置ISPIN=1
若计算磁性体系，建议先用ISPIN=1优化结构，再用ISPIN=2进行磁性计算
结构优化参数
对于无磁结构的几何优化，合理的参数设置至关重要：
ISTART = 0 # 从头开始计算
ICHARG = 2 # 读取原子电荷密度
PREC = Accurate # 高精度
ENCUT = 450-520 # 截断能（根据POTCAR调整）
EDIFF = 1E-4 # 电子步收敛标准
EDIFFG = -0.02~0.05 # 离子步收敛标准（力）
ISIF = 2/3 # 优化原子位置和晶胞
IBRION = 2/3 # 优化算法（准牛顿法或共轭梯度法）
NSW = 100-300 # 离子步最大数
POTIM = 0.1-0.5 # 时间步长
ISMEAR = 0/1 # 展宽方法
SIGMA = 0.1-0.2 # 展宽宽度
对称性处理
ISYM参数控制对称性的使用，对无磁计算有重要影响：
ISYM = 0 # 不使用对称性（适用于低对称性或磁性结构）
ISYM = 1 # 使用晶体对称性
ISYM = 2 # 使用晶体对称性但不处理原子对称性
建议：无磁结构优化时，通常使用ISYM=1以加快计算；若体系存在畸变或需要精确处理时，可使用ISYM=0。
实空间投影
对于大体系或金属体系，实空间投影可提高效率：
LREAL = Auto # 自动选择实空间投影
LREAL = .TRUE. # 显式开启实空间投影
LREAL = .FALSE.# 使用倒易空间投影
在某些情况下，需要在无磁计算基础上进行非共线磁性计算，相关参数包括：
LNONCOLLINEAR = .TRUE. # 开启非共线磁性
LSORBIT = .TRUE. # 开启自旋轨道耦合
SAXIS = 0 0 1 # 自旋量化轴方向
MAGMOM = N*0 0 0 # 每个原子三个磁矩分量（x, y, z）
结构文件构建
首先需要构建正确的POSCAR文件，包含：
晶格矢量（3行）
原子种类及数量（1行）
坐标类型选择（1行，如Direct或Selective Dynamics）
原子坐标（N行）
要点：确保结构合理性，避免原子重叠或异常大的晶格参数。
POTCAR准备
根据元素种类选择正确的赝势文件。VASP官网提供多种赝势库：
PAW_PBE：推荐使用的标准赝势
PAW：高精度赝势
GGA：梯度校正赝势
初步优化
# INCAR设置
SYSTEM = Structure Optimization
PREC = Normal
ENCUT = 400
EDIFF = 1E-4
ISPIN = 1
ISYM = 1
ISMEAR = 0
SIGMA = 0.1
IBRION = 2
NSW = 100
ISIF = 3
EDIFFG = -0.05
POTIM = 0.5
ALGO = Normal
高精度优化
初步优化后，进行更高精度的结构优化：
# INCAR设置
SYSTEM = High Precision Optimization
PREC = Accurate
ENCUT = 520
EDIFF = 1E-5
EDIFFG = -0.02
ISPIN = 1
ISYM = 2
ISMEAR = 0
SIGMA = 0.05
IBRION = 3
NSW = 200
ISIF = 3
ALGO = Fast
LREAL = .FALSE.
结构优化完成后，进行静态能量和性质计算：
# INCAR设置
SYSTEM = Static Calculation
PREC = Accurate
ENCUT = 520
EDIFF = 1E-6
ISPIN = 1
ISYM = 2
ISMEAR = 0
SIGMA = 0.01
IBRION = -1
NSW = 0
LWAVE = .TRUE. # 保存WAVECAR
LCHARG = .TRUE. # 保存CHGCAR
LORBIT = 11 # 输出投影态密度
无磁结构计算的基本原理：介绍了自旋极化与非自旋极化、无磁与顺磁的区别
INCAR关键参数设置：ISPIN参数、结构优化参数、对称性处理
计算流程与步骤：准备工作、结构优化流程、静态计算流程
下一章将正式引入本次教程的核心—VASP铁磁结构计算。我们将从铁磁结构的基本原理、INCAR关键参数设置、铁磁计算完整流程方面详细介绍VASP铁磁结构计算。
