VASP外加电场下SCF难以收敛的解决方法
`LDIPOL=.TRUE.
IDIPOL=1/2/3/4
EFIELD= ?
DIPOL=0.5 0.5 0.5`
外加电场下需要以上代码,其中LDIPOL负责偶极矫正,是需要开的
对于IDIPOL=1-3,偶极矩将分别平行于第一、第二或第三晶格矢量的方向计算。总能量的修正计算为当前超单元中单极子/偶极子和四极子与放置在超单元中的同一偶极子之间的能量差,对应的晶格矢量接近无穷大。此标志应用于板计算。
对于IDIPOL=4,将计算所有方向上的完整偶极矩,并将总能量的修正计算为当前超级单体中的单极/偶极/四极和放置在真空中的同一单极/偶极/四极之间的能量差,使用该标志计算孤立分子。
EFIELD 的设置需要注意单位的换算,单位是eV/A
DIPOL对应的值是你材料的位置,质心或者几何中心,可以通过MS给出,或者自己简单估计一个
VASP里会设置虚拟电板在中间,可以尝试调整一下结构,避免材料直接放在晶胞正中
SCF难以收敛有如下原因和解决方法:
1、EFIELD设置过大,建议从0.1或者更小开始,后面的WAVECAR和CHGCAR使用前面成功的
2、DIPOL设置不对,应该设置到材料上,不能再真空层里
3、ALGO改成Normal或Fast,EDIFF降低精度,提高SCF次数
4、CG方法(IBRION=2)也可能是原因,可以试试IBRION=1

浙公网安备 33010602011771号