total drift太大的问题
在做一个结构的弛豫的时候遇到了这个问题
total drift: 19.934551 34.617080 2.916653
一般这个都是接近0的,之前也一直以为它是上面力的总和。
在我这个例子里,无论力和压力收敛的再好,这个total drift始终很大,非常困扰。
从vaspwiki上还有冲浪得到的结果来看,他们提出的方法包括
增加截断能
增大FFT网格(NGXF)
关掉实空间投影 LREAL
三个都试了一圈无果。后来在组里师弟的帮助下发现是由于ISMEAR的SIGMA取的太大的缘故。
我之前取得都是0.2,改成0.02之后就没有这个问题了。
VASP 加SOC后受力异常的问题
最近在计算WSe2中的SOC,本来想的是这么流行的体系应该很好算,结果VASP狠狠给我上了一课。弛豫好的结构放到SOC里受力都是个位数的 eV/Angstrom。怎么改参数都不行,搜也搜不着,而且别人都用VASP算这个发了文章了。没办法只能一个参数一个参数的测试,最后发现竟然是vasp_ncl
版本的并行的问题。去掉NCORE
参数后受力就正常了,目测是vasp_ncl
的并行的问题。可以考虑使用KPAR
并行,这个经测试没问题
顺便搜了一下vasp并行策略的问题,根据官方的PPT里说的,KPAR是k点并行是最有效率的,NPAR是能带并行次之。
计算SOC的参数
#SOC related input
ISYM = -1
# ISPIN = 2
LSORBIT = .TRUE.
# LNONCOLLINEAR = .TRUE. 打开LSORBIT后这个会自动打开
GGA_COMPAT = .FALSE.
LMAXMIX = 6 ##我这里是WSe2,涉及f轨道需要设为6
VASP优化结构中断问题
我这里是因为没有设置EDIFF,VASP在优化结构的时候还没达到力收敛标准的时候就能量收敛了(因为程序取的默认值比我要的偏大)。所以报这个错误。加上一个EDIFF参数就好了。
ZBRENT: fatal error in bracketing please rerun with smaller EDIFF, or copy CONTCAR to POSCAR and continue
网友评论