计算化学公社

标题: NEB过渡态搜寻插点的结构能量比初态和末态都要低 [打印本页]

作者
Author: lxx-QUST 时间: 2023-12-4 08:57
标题: NEB过渡态搜寻插点的结构能量比初态和末态都要低
本帖最后由 lxx-QUST 于 2023-12-4 09:12 编辑

NEB过渡态搜寻插点的结构能量比初态和末态都低，插点之后计算出来的能量与路径的关系图如图一所示，我把01和07又进行了一遍结构优化，发现得到的能量与路径的关系图，没有找到问题具体所在，图二是我的INCAR设置，希望各位老师帮忙看一下，如何解决，谢谢！

作者
Author: pal 时间: 2023-12-4 09:12
说明不是从初始结构直接变成终态结构的呗，得优化那个能量更低的结构

作者
Author: lxx-QUST 时间: 2023-12-4 09:14

pal 发表于 2023-12-4 09:12
说明不是从初始结构直接变成终态结构的呗，得优化那个能量更低的结构

您好，那有没有什么好的方法来寻找那个能量更低的结构呢？

作者
Author: pal 时间: 2023-12-4 09:55

lxx-QUST 发表于 2023-12-4 09:14
您好，那有没有什么好的方法来寻找那个能量更低的结构呢？

直接把得到路径里那个能量更低的结构拿出来优化呗，只要那个结构不是太离谱

作者
Author: lxx-QUST 时间: 2023-12-4 10:08

pal 发表于 2023-12-4 09:55
直接把得到路径里那个能量更低的结构拿出来优化呗，只要那个结构不是太离谱

我尝试了几次您的方法，但是最后计算完NEB还是得到了一样的能量与路径关系图

作者
Author: pal 时间: 2023-12-4 11:14

lxx-QUST 发表于 2023-12-4 10:08
我尝试了几次您的方法，但是最后计算完NEB还是得到了一样的能量与路径关系图

所以你没有把01作为初态，07作为终态吧

作者
Author: zmjsce 时间: 2023-12-4 11:34
如果体系带有磁性的话，您可以看看IS TS FS三个结构优化后的mag参数，如果mag有显著差别的话能量可能不太准的

作者
Author: lxx-QUST 时间: 2023-12-4 14:15
本帖最后由 lxx-QUST 于 2023-12-4 14:21 编辑

pal 发表于 2023-12-4 11:14
所以你没有把01作为初态，07作为终态吧

我是把01作为初态，07作为终态，这样做是不对的吗？

作者
Author: lxx-QUST 时间: 2023-12-4 14:25

zmjsce 发表于 2023-12-4 11:34
如果体系带有磁性的话，您可以看看IS TS FS三个结构优化后的mag参数，如果mag有显著差别的话能量可能不太准 ...

谢谢您啦，我刚看了一下我优化之后的mag值他们三个是一样的

作者
Author: 哈库那摩塔塔 时间: 2023-12-4 19:43
你是00初态，08末态，01能量比00低，考虑01做初态，另外你初末态是差一点多埃吗，怎么插了这么多点，一般不是初末坐标差值÷0.8得插点数吗

作者
Author: lxx-QUST 时间: 2023-12-4 20:08

哈库那摩塔塔发表于 2023-12-4 19:43
你是00初态，08末态，01能量比00低，考虑01做初态，另外你初末态是差一点多埃吗，怎么插了这么多点，一般不 ...

对，初末态差1.05埃。
插点个数是因为老师要求我5个或者7个，然后我选择了7个。

作者
Author: 卡开发发 时间: 2023-12-4 21:53

lxx-QUST 发表于 2023-12-4 20:08
对，初末态差1.05埃。
插点个数是因为老师要求我5个或者7个，然后我选择了7个。

两个能量低的抓出来重新优化下，看看能量会不会更低，如果是这样，用这两个来插点找过渡态看看。

作者
Author: lxx-QUST 时间: 2023-12-5 07:59

卡开发发发表于 2023-12-4 21:53
两个能量低的抓出来重新优化下，看看能量会不会更低，如果是这样，用这两个来插点找过渡态看看。

您说的方法我重复了四次，上面的图是我第五次按照您的方法算出来的结果，所以我一直不知道问题出在哪

作者
Author: 123wo 时间: 2023-12-5 08:36
用能量较低的重新作为基态做一下试试

作者
Author: lxx-QUST 时间: 2023-12-5 09:12

123wo 发表于 2023-12-5 08:36
用能量较低的重新作为基态做一下试试

好的，谢谢啦

作者
Author: 卡开发发 时间: 2023-12-5 10:17

lxx-QUST 发表于 2023-12-5 07:59
您说的方法我重复了四次，上面的图是我第五次按照您的方法算出来的结果，所以我一直不知道问题出在哪

你可能需要做一下收敛性测试。

作者
Author: lxx-QUST 时间: 2023-12-7 09:09

123wo 发表于 2023-12-5 08:36
用能量较低的重新作为基态做一下试试

好的，谢谢啦

作者
Author: Observer-Aya 时间: 2024-4-16 10:49

lxx-QUST 发表于 2023-12-5 07:59
您说的方法我重复了四次，上面的图是我第五次按照您的方法算出来的结果，所以我一直不知道问题出在哪

您好，请问这个问题解决了吗，我最近也遇到了类似的问题，明明末态和初态几何优化后的能量相同，过渡态搜索时得到的末态却比初态和过渡态都高

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