计算化学公社

标题: 求助：过渡态能量计算与单点能结果不一致 [打印本页]

作者
Author: wrtgcn 时间: 2019-9-4 10:33
标题: 求助：过渡态能量计算与单点能结果不一致
粗算下得到了收敛的过渡态，跑精算时个别结构的电子步难以收敛，就想保存粗算结果的WAVECAR帮助计算。对已收敛的结构做了同一级别下的单点计算，却发现得到的能量趋势与过渡态的能量结果并不一致。单点能结果是能量一直向下，过渡态计算的结果则为能量先增再减。
想请问各位老师和同学，过渡态计算实际也是优化计算，为什么得到的能量趋势会与单点计算不同？非常感谢大家。

作者
Author: 卡开发发 时间: 2019-9-4 14:15
优化计算的过程中粗略得到的势能面和精确的势能面有差异啊，要不然就不用强调收敛性测试的问题了啊。

作者
Author: wrtgcn 时间: 2019-9-4 14:53

卡开发发发表于 2019-9-4 14:15
优化计算的过程中粗略得到的势能面和精确的势能面有差异啊，要不然就不用强调收敛性测试的问题了啊。

感谢老师的回答。我的ENCUT和收敛标准都是按文献中取的，跑粗算的时候只减小了K点数并固定了大部分原子。没想到算出来的能量趋势都不正确。还想请问您，这种情况是否应该提高精度再跑粗算？以及如果优化计算得到的是粗糙的势能面，得到能量的相对大小关系都不正确，是否总是有必要在优化后重新计算单点能？初学VASP，可能问题比较基础，再次感谢老师。

作者
Author: 卡开发发 时间: 2019-9-4 16:25

wrtgcn 发表于 2019-9-4 14:53
感谢老师的回答。我的ENCUT和收敛标准都是按文献中取的，跑粗算的时候只减小了K点数并固定了大部分原子。 ...

1、即便粗略计算精度也不要放得太宽，否则给出的势能面可能定性都不一定对（比如毛糙），另外我不知道如果文献没做固定你为啥要做固定，这么做也基本不会减小计算量。

2、势能面定性大小都不正确的情况得到的构型都不一定是你想要的，这时候算得的单点能也没啥意义，如果结构比较接近也建议提高精度优化一下再算单点。单点自洽计算在离子步出现难收敛以及做变胞优化的情况都非常有必要单独做一下。

3、我也很难保证文献给的结果一定正确。

作者
Author: wrtgcn 时间: 2019-9-4 16:53

卡开发发发表于 2019-9-4 16:25
1、即便粗略计算精度也不要放得太宽，否则给出的势能面可能定性都不一定对（比如毛糙），另外我不知道如 ...

感谢您的指导。固定原子主要是参考了Learn VASP The Hard Way教程，里面说粗算时取的K点数太少可能会出现原子乱跑的情况，所以固定了那些不是特别关注的部分。

可能我的体系比较特殊，里面有臭氧分子，在过渡态计算中还碰到了磁矩反复变化影响能量的情况。现在尝试固定了总磁矩，好不容易跑出了能量趋势，没想到势能面的定性正确都没达到。

关于您说的第二点，可以理解为在过渡态计算中（精度足够且结构没太大问题时），如果出现了离子步难以收敛的情况，可以算一个单点再读取WAVECAR作为初猜来帮助收敛吗？

作者
Author: 卡开发发 时间: 2019-9-4 19:08

wrtgcn 发表于 2019-9-4 16:53
感谢您的指导。固定原子主要是参考了Learn VASP The Hard Way教程，里面说粗算时取的K点数太少可能会出现 ...

1、

里面说粗算时取的K点数太少可能会出现原子乱跑的情况

我觉得精度尽量不要太低就可以避免这个问题，固定很多情况并不能解决什么实际问题。

2、60步电子步还没收敛的趋势应该考虑下结构或磁序方面是否设置正确，参数是否合理等。

作者
Author: wrtgcn 时间: 2019-9-4 19:29

卡开发发发表于 2019-9-4 19:08
1、
我觉得精度尽量不要太低就可以避免这个问题，固定很多情况并不能解决什么实际问题。

非常感谢老师的帮助。

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