VASP结构优化后再做单点磁矩及能量不一致的问题

LAA

想请教一下大家，有没有遇到过类似的问题，一共109个原子的体系。
具体操作如下：
我在结构优化后，想再做个单点能（静态计算），就把结构优化后的CONTCAR文件复制为POSCAR，INCAR的NSW设置为1，然后其他文件保持一致。算出来发现电子步没有收敛，能量和磁矩差别也很大（设置的60步），遂又把NSW改为原来的100（也就是和结构优化一致了），然后结果发现总磁矩也是差了2-4。进一步看力的情况，力收敛设置的是-0.02（EDIFFG = -0.02），结构优化最后一步最大原子受力为0.019,但是把这个CONTCAR拿去再算一次，第一步的最大原子受力变为0.571了。按理来说应该也是0.019才对吧？感觉好像是初猜的不一样，导致磁矩不一样，想问下这个是什么情况，如何解决呢。谢谢大家！

参数说明：
ISPIN = 2，打开了自旋，但是没设置初始磁矩，用默认的。
NELECT ，因为体系少一个电子，所以这边总电子数-1
IBRION = 2      
ISIF = 2

LREAL = Auto

K点221，Gamma

白皛白

做单点的时候最好读取结构优化时的WAVECAR和CHGCAR。具体请看这个 https://www.bilibili.com/read/cv10064345/

ahxb

两次计算收敛到的波函数不相同，自然波函数的导数也不相同。想保持一致需要读结构优化完成后输出的WAVECAR。

LAA

白皛白 发表于 2023-10-16 10:18
做单点的时候最好读取结构优化时的WAVECAR和CHGCAR。具体请看这个 https://www.bilibili.com/read/cv100643 ...

谢谢~我有试过，这样就可以保证一致，实际情况是跑两个电子步就收敛了，但是每次都这样的话，那文件实在太大了

LAA

ahxb 发表于 2023-10-16 14:18
两次计算收敛到的波函数不相同，自然波函数的导数也不相同。想保持一致需要读结构优化完成后输出的WAVECAR ...

感谢~两次计算出的波函数出来不一样这个确实，但是主要是力场也不一样就感觉很离谱。因为第一次结构优化力已经收敛，再继续拿着个优化后的CONTCAR来再跑结构优化，居然第一步的力大了很多，理论上应该和上一次结构优化的最后一步的力一致才对吧

hebrewsnabla

LAA 发表于 2023-10-17 16:39
感谢~两次计算出的波函数出来不一样这个确实，但是主要是力场也不一样就感觉很离谱。因为第一次结构优化 ...

力是由波函数决定的。