计算化学公社

标题: vasp，结构优化和自洽计算结果不一致问题 [打印本页]

作者
Author: Shana 时间: 2022-5-15 11:59
标题: vasp，结构优化和自洽计算结果不一致问题
各位好，有时候做完结构优化后，用弛豫完的结构做自洽计算，得到的结果和结构优化最后的离子步结果不一致，差距在零点几~几eV，而有时又十分吻合，请问，这是什么原因造成的？无论是结构优化还是自洽计算，电子步都是收敛的。

作者
Author: jiangning198511 时间: 2022-5-15 12:50
优化晶胞参数了吗？如果晶胞体积有变化的话两次计算的平面波数量会不完全一样导致误差；当然这个问题不大

作者
Author: 吴双 时间: 2022-5-15 14:29
优化还有自洽的时候电子和离子步的收敛线是一致的么

作者
Author: Miles 时间: 2022-5-15 14:34
你所指的结果不一致是能量不一致吗？如果是的话，结构优化是离子步中嵌套电子步，当最后一个离子步收敛时，计算就停止了，认为达到收敛，但此时电子步并未收敛，所以需要静态自洽（单点能计算），单点计算时结构不变（NSW=0/1或 IBRION=-1），只是电子步优化，收敛后即最准确的能量。两者一般来说相差很小。

作者
Author: Shana 时间: 2022-5-15 15:31

jiangning198511 发表于 2022-5-15 12:50
优化晶胞参数了吗？如果晶胞体积有变化的话两次计算的平面波数量会不完全一样导致误差；当然这个问题不 ...

优化了，明白您的意思了。这样的话，在结构弛豫完后，最后再做一次自洽计算比较好吧

作者
Author: Shana 时间: 2022-5-15 15:31

吴双发表于 2022-5-15 14:29
优化还有自洽的时候电子和离子步的收敛线是一致的么

基本一致

作者
Author: Shana 时间: 2022-5-15 15:34

Miles 发表于 2022-5-15 14:34
你所指的结果不一致是能量不一致吗？如果是的话，结构优化是离子步中嵌套电子步，当最后一个离子步收敛时， ...

明白您的意思了，我隐约有这种感觉，您这一解释我就清晰了。看来，每一次结构优化完毕后，再做一次自洽计算比较好

作者
Author: Miles 时间: 2022-5-15 16:41

Shana 发表于 2022-5-15 15:34
明白您的意思了，我隐约有这种感觉，您这一解释我就清晰了。看来，每一次结构优化完毕后，再做一次自洽计 ...

本着严谨的态度，计算形成能、吸附能、分离能等能量时，都应该做单点计算得到电子步收敛时的能量。当然有的偷懒只用结构优化的能量来表示，差异也不明显

作者
Author: Shana 时间: 2022-5-15 17:27

Miles 发表于 2022-5-15 16:41
本着严谨的态度，计算形成能、吸附能、分离能等能量时，都应该做单点计算得到电子步收敛时的能量。当然有 ...

好的，明白了

作者
Author: Enrus 时间: 2024-11-19 18:06
转个计算关于该问题的帖子：https://www.bilibili.com/opus/497026122446620883

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