这一节讲一下计算缺陷的例子，就是有的体系会出现空位。最常见的是氧空位。这一节以TiO2为例子介绍氧空位的计算。主要是介绍小编写的一个python脚本，怎么去固定原子计算，固定原子计算，固定原子计算！！！（重要的事情说三遍）
工具/原料
vasp
方法/步骤
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，下期见！！！
上一篇：白内障做了一只眼睛怎样配镜
下一篇：烈焰单人副本诀窍分享