计算化学公社

标题: VASP结构优化后再做单点磁矩及能量不一致的问题 [打印本页]

作者
Author: LAA 时间: 2023-10-16 09:52
标题: VASP结构优化后再做单点磁矩及能量不一致的问题
想请教一下大家，有没有遇到过类似的问题，一共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

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

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

作者
Author: LAA 时间: 2023-10-17 16:36

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

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

作者
Author: LAA 时间: 2023-10-17 16:39

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

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

作者
Author: hebrewsnabla 时间: 2023-10-17 16:55

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

力是由波函数决定的。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)