使用CI-NEB方法搜索Li3N中Li离子迁移的过渡态，初末态能量发生了改变

limpa · 发表于 Post on 2025-4-17 19:32:44

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

我的初末态优化后xyz文件、几何优化与过渡态搜索inp与ener文件随附在下。

geo_opt.inp (1.2 KB, 下载次数 Times of downloads: 1)

band_1.inp (2.11 KB, 下载次数 Times of downloads: 4)

li3n.neb-1.ener (1.56 KB, 下载次数 Times of downloads: 1)

opt_end_1.xyz.txt (4.52 KB, 下载次数 Times of downloads: 4)

opt_start_1.xyz.txt (4.52 KB, 下载次数 Times of downloads: 3)

maoxinxina · 发表于 Post on 2025-4-18 11:27:33

找到问题所在了吗

sobereva · 发表于 Post on 2025-4-23 15:32:21

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

limpa · 发表于 Post on 2025-4-23 17:25:52

感谢老师们。我将ALIGN_FRAMES与ROTATE_FRAMES设为F后，末态能量和结构优化后的能量一致了，跑出的过渡态也比较正常。

limpa · 发表于 Post on 2025-4-23 17:31:35

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

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

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

sobereva · 发表于 Post on 2025-4-23 21:29:49

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

参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）幻灯片

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

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