在上一章中，华算科技朱老师详细介绍了VASP示例任务。VASP计算一般从结构优化开始，然后对优化后的结构进行性质计算。本章将介绍VASP结构优化计算实例，让大家更好的了解计算流程，具体包括结构优化、半导体性质案例、表面与催化案例。
吸附体系主要分为气相吸附（如CO吸附在金属表面）和液相/电化学吸附（如水分子在氧化石墨烯上的吸附）。对于气相吸附，通常采用真空层模型；对于电化学体系，则需加入电极势能修正。
VASP默认采用投影增强平面波（PAW）方法。对于过渡金属表面，常用的交换关联泛函（XC）是PBE-GGA。如果涉及强关联系统（如FeO表面），可能需要加上U修正（DFT+U）。
表面计算通常需要较大的真空层（>15 Å）和较密的K点网格（如(4x4x1)），因此计算资源需求高于体相计算。收敛性测试（能量截断、K点密度）是保证结果可靠性的前提。
构建表面模型是吸附计算的第一步，直接决定了后续吸附位点的可行性。
使用Materials Studio、VESTA或ASE等工具，将体相晶体切割成特定的晶面（如(111)、(100)等）。对于金属表面，常见的高对称性面是(111)面。
为了模拟半无限的表面，需要在计算模型中加入真空层。一般建议：
层数：至少3-5层原子（取决于表面类型）。通常固定底层1-2层，放松表面层。
真空层：>15 Å，以消除周期性复制产生的相互作用。
对于研究缺陷位点吸附（如空位、取代），需要先在表面构建缺陷结构，再进行吸附。缺陷位点通常具有更高的活性。
表面通常包含多个高对称性位点。选择合适的吸附位点是确定吸附构型的关键。
以金属(111)面为例，主要有以下四种高对称性位点：
Top位
吸附原子直接位于表面原子上方（上方原子）。适用于原子吸附或小分子键长较短的情况。
Bridge位
吸附原子位于两个相邻表面原子之间的桥接位置。适用于双原子分子（如O2、N2）的平行吸附。
Hollow位（凹位）
位于三个或四个表面原子形成的空洞中心。分为：
Fcc Hollow
位于面心立方结构的凹位。
Hcp Hollow
位于六方密堆结构的凹位。
对于环状分子或大分子，Hollow位常是最稳定的吸附位点。
除了位点本身，分子的取向（Orientation）同样重要。例如，CO分子可以垂直吸附（C端或O端向下），也可以平行吸附。对于环状分子（如苯），可以平面平行于表面（π-π相互作用）或垂直吸附。
在构建吸附模型时，需要人为设置吸附物与表面原子之间的初始距离（Initial Height）。通常：
原子吸附：距离约为2.0-2.5 Å。
分子吸附：距离约为3.0-3.5 Å（取决于键长和范德华半径）。
如果距离设置过近，优化计算可能会直接导致重叠错误；距离过远，则计算收敛时间会增加。
1.打开VESTA：加载已优化好的表面POSCAR文件。
2.选择位点：使用“Add Atom”功能，在目标位点上添加吸附原子。
3.调整距离：通过编辑坐标手动调整z轴坐标，确保符合上述建议距离。
4.保存：将新文件保存为吸附体系的POSCAR。
构建好POSCAR后，需要准备INCAR、KPOINTS和POTCAR文件。
针对吸附优化，常用的INCAR参数如下：
SYSTEM = CO_adsorption_on_Au111
ISTART = 0
ICHARG = 2
ENCUT = 400
EDIFF = 1E-5
EDIFFG = -0.02 # 力收敛阈值 (eV/Å)
ISMEAR = 0
SIGMA = 0.05
ISIF = 2 # 固定晶格，仅优化原子位置
IBRION = 2 # 共轭梯度法
NSW = 100 # 最大离子步数
LCHARG = .TRUE. # 输出电荷密度
LWAVE = .FALSE. # 不输出波函数
关键点：ISIF=2非常重要，它确保表面晶格参数固定，仅允许原子（包括吸附物）移动。
对于大真空层模型，建议使用较稀疏的z方向网格：
Automatic mesh
0
Monkhorst-Pack
4 4 1
0 0 0
这里的1表示z方向只有一个k点。
确保包含了表面材料和吸附物的所有元素。例如，吸附CO在Au(111)上，需要包含Au、C、O三种元素的POTCAR。
吸附能是衡量吸附强度的关键指标，计算公式为：
E_{ads} = E_{slab+adsorbate} – E_{slab} – E_{adsorbate}
其中：
E_{slab+adsorbate}：吸附体系的总能量（优化后）。
E_{slab}：纯净表面的总能量（单独计算）。
E_{adsorbate}：吸附物（如CO）的总能量（在同样大小的真空盒中单独计算）。
注意：为了消除计算误差，E_{slab}和E_{adsorbate}的计算条件（KPOINTS、ENCUT等）必须与E_{slab+adsorbate}完全一致。
模型选择：介绍了如何选择计算方法、计算资源、计算模型
吸附位点选取与建模：介绍了吸附位点，分子构型，初始层间距
输入文件：介绍了INCAR、POSCAR、POTCAR、KPOINTS文件
下一章将正式引入本次教程的核心—吸附能计算。我们将从吸附能意义、计算流程、结果分析方面详细介绍VASP 表面与催化案例，敬请期待！
在学习和实操过程中遇到问题？欢迎关注DFT计算参谋长，获取更多独家学习资源、实用脚本工具包，并加入活跃的学习交流群，与众多同行一起交流进步！
#华算科技##第一性原理#DFT#VASP#MS#高斯#表面#吸附#催化#建模#吸附点位#VASP基础教程#分子构型
