计算化学公社

标题: 使用CI-NEB方法搜索Li3N中Li离子迁移的过渡态，初末态能量发生了改变 [打印本页]

作者
Author: limpa 时间: 2025-4-17 19:32
标题: 使用CI-NEB方法搜索Li3N中Li离子迁移的过渡态，初末态能量发生了改变
老师们，我在搜索Li3N中Li离子迁移的过渡态，构建了（3*3*2）的超胞，挖去两个相邻离子作为初末态并进行优化，然后使用CI-NEB方法搜索过渡态。按理来说初末态能量应该一样，但是问题在于在ener文件中初态能量不变，末态能量发生了增大，甚至高于中间过渡态。通过调整精度与结构优化一致后也还是这种结果。想请教老师们问题在哪儿，望不吝赐教。

我的初末态优化后xyz文件、几何优化与过渡态搜索inp与ener文件随附在下。
(, 下载次数 Times of downloads: 1)
(, 下载次数 Times of downloads: 4)
(, 下载次数 Times of downloads: 1)
(, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 3)

作者
Author: maoxinxina 时间: 2025-4-18 11:27
找到问题所在了吗

作者
Author: sobereva 时间: 2025-4-23 15:32
末态能量一直是-575.87026没变。更后面位数轻微波动是并行运算带来的随机性所致。如果你指ener文件更后面的列，那是NEB点之间的间距。

作者
Author: limpa 时间: 2025-4-23 17:25
感谢老师们。我将ALIGN_FRAMES与ROTATE_FRAMES设为F后，末态能量和结构优化后的能量一致了，跑出的过渡态也比较正常。

作者
Author: limpa 时间: 2025-4-23 17:31

sobereva 发表于 2025-4-23 15:32
末态能量一直是-575.87026没变。更后面位数轻微波动是并行运算带来的随机性所致。如果你指ener文件更后面的 ...

老师，我这里表述有问题。我的意思是我对末态结构优化后能量为-575.87067，但在跑CI-NEB时能量变成了-575.87026，所以比较困惑。但把ALIGN_FRAMES与ROTATE_FRAMES设为F后，能量与结构优化一致了。

但我还是困惑什么时候应该将ALIGN_FRAMES与ROTATE_FRAMES关闭，还不太清楚这两项的意义，希望老师解答。

作者
Author: sobereva 时间: 2025-4-23 21:29

limpa 发表于 2025-4-23 17:31
老师，我这里表述有问题。我的意思是我对末态结构优化后能量为-575.87067，但在跑CI-NEB时能量变成了-575 ...

参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）幻灯片
(, 下载次数 Times of downloads: 2)

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