计算化学公社

标题: 小白求助，vasp计算LaNi5（100）面吸附CO吸附能时，计算收敛后OSZICAR的E0数值和E.... [打印本页]

作者
Author: CCYY 时间: 2024-10-24 21:20
标题: 小白求助，vasp计算LaNi5（100）面吸附CO吸附能时，计算收敛后OSZICAR的E0数值和E....
本帖最后由 CCYY 于 2024-10-24 21:32 编辑

本人初次发帖，是初学vasp计算的一名小白，希望大家可以看一个这样的问题。
在 VASP 计算中，当系统达到收敛时，E0 和 E 应该是非常接近的，但不一定完全相等。

E0 通常指的是在离子优化或电子优化中计算得到的能量，表示基态能量。
E 是当前计算步骤的能量

然而我在复现别人文章支持信息中的吸附能数据时，发现OSZICAR中的能量出现了问题，收敛后的能量过大；在收敛过程中E和E0差别很大，如果用最后的能量的话，吸附能达到了-8.9eV，如果用过程中的E能量计算，吸附能是合理的。
下面是我附上的INCAR、POSCAR、OUTCAR、OSZICAR
(, 下载次数 Times of downloads: 1)

[attach]100269[/attach]
[attach]100270[/attach]
[attach]100271[/attach]
[attach]100272[/attach]

作者
Author: 白皛白 时间: 2024-10-25 17:25
E是当前自由能，F是总自由能（除去参考原子的能量），E0是σ—>0的能量我的理解是消除了展宽影响的能量，然后OUTCAR里还有个energy without entropy这个不含熵的能量等于F-TS，这个TS在OUTCAR里搜索EENTRO可以找到，S实际上是电子熵。这里面E0和1/2*(F+不含熵的能量）很接近(如果INCAR里ISMEAR=-5的话F=不含熵的能量=E0)。
所以感觉E和E0之间的差距实际上取决于当前自由能和总自由能之间的差距，关于E和F之间具体有什么联系我也查了很久也没有搞明白，希望有老师路过这个帖子可以解惑，从我个人经验上来看体系的原子数越多这两个差距就越大，金属原子多的体系差距也会更大，我认为E和F有差距甚至很大的差距应该是正常现象，不是导致吸附能不合理的原因。
感觉一个方面是这个盒子本身太小了，b方向只有4埃的长度，这很可能导致相邻盒子的CO之间相互作用，还有一个可能是催化剂本身在吸附CO的过程中重构了，这个要看CONTCAR才能知道。

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