什么时候设置偶极修正?
VASP官网的例子中,在计算Ni(111)表面弛豫时没有用偶极修正,但是在计算CO在Ni表面的DOS和work function时使用了LDIPOL,另一个H2O在TiO2表面时计算表面弛豫时也没有用偶极修正。但是在大师兄科研网中提到:偶极矫正:LDIPOL = .TRUE. ; IDIPOL = 3 Slab模型一般都需要加上这个参数。网址如下:https://www.bigbrosci.com/2018/08/26/ex43/?highlight=dipol
看大师兄科研网的例子,无论是表面弛豫还是功函数的计算都加入了偶极修正。
那种方法是对的?
INCAR for Ni(111)
ISTART = 0ICHARG = 2 general:SYSTEM = CO adsorption on Ni(111)ENMAX= 400ISMEAR = 2; SIGMA = 0.2ALGO= FastEDIFF = 1E-6 dynamic:NSW=100POTIM = 0.2IBRION = 1
INCAR for CO on Ni(111)
general:ENMAX = 400SYSTEM = CO adsorption on Ni(111)ISMEAR = -5ALGO = Fast LDOS:LORBIT = 11 workfunction:IDIPOL = 3LDIPOL = .TRUE.LVHAR = .TRUE.#LVTOT = .TRUE.
INCAR for H2O on TiO2(110)
SYSTEM = H2O_TiO2
ENMAX = 400
ISMEAR =2
SIGMA = 0.05
EDIFF = 1e-6
EDIFFG = -0.05
IBRION = 2
POTIM = 0.5
NSW = 200
一般如果上面表面不对称,需加修正。 那么就是官网的例子错了?
页:
[1]