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