计算金属表面Li的扩散时，加了DFTD3色散校正使得能垒为负，不知为何？

1138711019

第一次计算NEB这样，反复计算了很多次，最后检查色散校正的时候发现了，不知道各位大佬有什么解决方吗？

卡开发发

你的两次计算的两个始末态的结构是一样的，还是各自在各自的方法（也就是分别考虑DFTD3）下进行的优化？
如果没有各自用各自的方法分别优化，那么应该分别优化后分别重新做路径的插值。
如果是各自用了各自的方法优化，很显然第二张图的中还存在其他极小点，建议把这些结构挑出来重新优化一下，可能能拿到能量更低的两个结构，以此作为始末态建立路径插值后再来做NEB试试。

1138711019

卡开发发 发表于 2022-9-29 08:05
你的两次计算的两个始末态的结构是一样的，还是各自在各自的方法（也就是分别考虑DFTD3）下进行的优化？
...

谢谢老师，采用的各自用各自的方法优化的，加了DFTD3进行几何优化一直优化不到能量最低那个态，然后计算NEB的时候就发现后面一个点能量比之前几何优化的还低（这两个点差距很小，只有0.244埃，而且提高几何优化的K点，降低一点NEB的K点也会出现这个现象）。
另外，我发现把NEB过程中DFTD3最后一步校正能量手动去掉（可以得到能垒0.15eV），几乎也能得到 不加DFTD3校正的能垒（0.14eV）。而加了DFTD3校正的能量波动很小（约0.03~0.08eV左右），是否能说这个表面扩散受色散作用为主要呢？

卡开发发

1138711019 发表于 2022-9-29 09:46
谢谢老师，采用的各自用各自的方法优化的，加了DFTD3进行几何优化一直优化不到能量最低那个态，然后计算N ...

你如果把NEB当中能量低点抽出来优化一下试试看呢？而不要用原来的初始结构。如果还不行，建议做一下收敛性测试看看能量和力的计算在你的计算精度下是否充足。

1138711019

卡开发发 发表于 2022-9-29 13:04
你如果把NEB当中能量低点抽出来优化一下试试看呢？而不要用原来的初始结构。如果还不行，建议做一下收敛 ...

好的，我先再算一次NEB。
谢谢老师~