vasp入門到精通[31]計算缺陷(或固定原子)

這一節講一下計算缺陷的例子,就是有的體系會出現空位。最常見的是氧空位。這一節以TiO2為例子介紹氧空位的計算。

主要是介紹小編寫的一個python腳本,怎麼去固定原子計算,固定原子計算,固定原子計算!!!(重要的事情說三遍)

工具/原料

vasp

方法/步驟

以TiO2為例子,單胞如下圖:

做養空位很簡單,就是去掉其中的一個O原子。但是實際上氧空位含量很少,也就是要做一個超胞,這個看你做多少含量的氧空位了,比如;

一個原胞里里有4個O原子,你拿掉一個,就是25%的O空位,這個數值太大了

那如果我想做個比例小的怎麼辦,很簡單,建立超胞就可以了,如果我建立一個2x2x2的超胞,拿掉一個O原子,那麼O空位就是3.125%,如下圖,把黃色的原子拿掉

然後進行優化計算,但是超胞的計算量太大,沒有必要全部優化,我們可以把下面幾層原子固定住:

在POSCAR中做如下操作:

如上圖在Direct上面加入Selective:

Selective

Direct

然後在遠離空位的層固定住就可以了,比如你去掉的原子坐標是0.1 0.1 0.1

那麼你超胞的下面可以固定住

不固定的在坐標後面加上T T T

固定的在坐標後面加上F F F

INCAR就是正常優化:

計算完就可以看一下能量。

你可以分別去掉不同的O原子,對比哪一個空位更穩定

小編寫了個小腳本,固定原子的,比如原子幾百個,你怎麼修改POSCAR呢,難道要一行行去加參數麼,太麻煩,所以小編寫了個python的小腳本,如果有需要就留言吧,我下一節就放上去,如果沒人留言我就當大家不感興趣了,就不再更新這個小腳本教程了

ok,下期見!!!

本文內容整理自網絡, 文中所有觀點看法不代表淘大白的立場