吸附能可能低至-48eV吗？可能的原因是什么呢？

卡开发发

jiangning198511 发表于 2022-1-10 16:53
Co的POTCAR为啥不一致，POTCAR:    PAW_PBE Co 02Aug2007；POTCAR:   PAW_PBE Co 06Sep2000？

应该就是这个问题，要不然三个体系的EATOM差距不会那么大。

Daisy_zju

jiangning198511 发表于 2022-1-10 16:53
Co的POTCAR为啥不一致，POTCAR:    PAW_PBE Co 02Aug2007；POTCAR:   PAW_PBE Co 06Sep2000？

谢谢老师指出的问题，POTCAR的不一致应该是因为一个是我用Vaspkit特地生成的，一个是生成KPOINTS的时候产生的，计算的时候也没有检查两者的一致性（万分没想到）。

Daisy_zju

卡开发发 发表于 2022-1-10 17:10
应该是楼上说的赝势问题，其他影响应该是次要的。道理上说，ALGO一般而言影响很小，因为只是求解本征值的算 ...

应该是的，我之后重算一遍，谢谢老师的指导和分享！

jiangning198511

Daisy_zju 发表于 2022-1-10 21:00
谢谢老师指出的问题，POTCAR的不一致应该是因为一个是我用Vaspkit特地生成的，一个是生成KPOINTS ...

统一用VASPKIT生成 应该就一致了

Daisy_zju

jiangning198511 发表于 2022-1-11 08:57
统一用VASPKIT生成 应该就一致了

老师您好，不好意思再麻烦您一下，我昨晚统一赝势之后重新计算，现在吸附结构的计算还没完成但是几近收敛，发现虽然吸附能有所改善，但是数值还是在-34eV左右，附件里打包了我的吸附质（OUTCAR-1）、纯表面(OUTCAR-2)和吸附结构（OUTCAR-3）的结果以及吸附过程优化的参数文件。不知道老师能不能再帮我看一看，谢谢！

卡开发发

Daisy_zju 发表于 2022-1-11 17:20
老师您好，不好意思再麻烦您一下，我昨晚统一赝势之后重新计算，现在吸附结构的计算还没完成但是几近收敛 ...

这次的问题来源于只有OUTCAR-3的体系做了DFT-D3(zero)的修正（也就是IVDW=11），而其他的两个都没做，这部分的能量贡献为Edisp (eV)  -34.45662。

另一个问题是纯Co体系的ENAUG和其他两个不相同，可能是简单设置PREC导致的，不过这个不是导致你结果异常的最大因素。

jiangning198511

Daisy_zju 发表于 2022-1-11 17:20
老师您好，不好意思再麻烦您一下，我昨晚统一赝势之后重新计算，现在吸附结构的计算还没完成但是几近收敛 ...

你的问题是输入文件中关键参数要保持一致，最好保持一致(当然与系统相关的需要修改)，另外 需要极度仔细的检查自己的输出文件

Daisy_zju

卡开发发 发表于 2022-1-12 09:04
这次的问题来源于只有OUTCAR-3的体系做了DFT-D3(zero)的修正（也就是IVDW=11），而其他的两个都没做，这 ...

谢谢老师的指导，修改之后吸附能正常了，麻烦老师了！

Daisy_zju

jiangning198511 发表于 2022-1-12 13:09
你的问题是输入文件中关键参数要保持一致，最好保持一致(当然与系统相关的需要修改)，另外 需要极度仔细 ...

是的，经卡开发发老师提醒后发现小分子计算的INCAR出了点问题，谢谢老师的指导！