VAPS结构优化（3）——Slab模型的固定

二维金属材料在计算的时候往往需要切面，默认的结构优化是允许所有的原子在三个方向自由弛豫，但是有的时候就需要固定部分原子。比如Slab模型，一般都是固定2~3层，弛豫的原子层数大于或者等于固定的原子层数。这里给出了固定slab模型POSCAR的部分方法（这里默认大家都已经切好Slab模型了）。

1. 手动方法（最耗时，不建议）

• 在POSCAR第7行和第8行之间插入一行， Selective Dynamics 前面讲过VASP只认第一个字母，也就是S是必须的。

• 加入Selective之后，需要在每一行原子坐标的后面加上T或者F表示允许和禁止移动。这里我们需要加三个T或者F，表示在x y z 三个方向上选择性固定原子的移动，可以使用notepad编辑器加入，也可以使用vi编辑器，这里就不多加赘述

2. 根据原子层数固定，使用vaspkit——402（注意此方法只适用于固定z方向的slab模型）

流程：vaspkit因为是POSCAR,所以选1 Threshold是vaspkit的总层数判断标准，在不同的标准下层数会判断的不太一样，选择Threshold的数值,这里选择1，判断出来5层我们固定三层，所以输入3检查一下生成的 POSCAR_fix 文件，如果没有问题：cp POSCAR_fix POSCAR

3. 根据原子z坐标高度固定原子层数，vaspkit——403

首先要弄清楚我们需要固定哪些原子，和这些原子的z坐标，比如 Au(100)固定下两层原子。找到我们要固定的原子中z坐标最大及最小的原子的坐标

4. 使用脚本固定, 知乎作者在这篇文章中，给出了一个脚本，https://zhuanlan.zhihu.com/p/397858779，有兴趣可以调用。（建议使用）