vasp入门到精通[31]计算缺陷（或固定原子）

1、以TiO2为例子，单胞如下图：

2、做养空位很简单，就是去掉其中的一个O原子。但是实际上氧空位含量很少，也就是要做一个超胞，这个看你做多少含量的氧空位了，比如;一个原胞里里有4个O原子，你拿掉一个，就是25%的O空位，这个数值太大了

3、那如果我想做个比例小的怎么办，很简单，建立超胞就可以了，如果我建立一个2x2x2的超胞，拿掉一个O原子，那么O空位就是3.125%,如下图，把黄色的原子拿掉

4、然后进行优化计算，但是超胞的计算量太大，没有必要全部优化，我们可以把下面几层原子固定住：在POSCAR中做如下操作：

5、如上图在Direct上面加入Selective：SelectiveDirect然后在远离空位的层固定住就可以了，比如你去掉的原子坐标是0.1 0.1 0.1那么你超胞的下面可以固定住不固定的在坐标后面加上T T T固定的在坐标后面加上F F F

6、INCAR就是正常优化：

7、计算完就可以看一下能量。你可以分别去掉不同的O原子，对比哪一个空位更稳定

8、小编写了个小脚本，固定原子的，比如原子几百个，你怎么修改POSCAR呢，难道要一行行去加参数么，涂尉郑堑太麻烦，所以小编写了个python的小脚本，如果有需要统潇瘵侃就留言吧，我下一节就放上去，如果没人留言我就当大家不感兴趣了，就不再更新这个小脚本教程了

9、ok，下期见！！！