Alma 发表于 2022-1-7 22:38:52

什么时候设置偶极修正?

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




vaspkit 发表于 2022-1-8 14:42:46

一般如果上面表面不对称,需加修正。

Alma 发表于 2022-1-8 16:49:38

那么就是官网的例子错了?
页: [1]
查看完整版本: 什么时候设置偶极修正?