计算化学公社

 找回密码 Forget password
 注册 Register

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

查看数: 9217 | 评论数: 23 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-1-10 10:54

正文摘要:

各位老师好我又来求助了,本人在计算苯胺分子在钴金属表面的吸附能,计算选取slab模型为3层,真空层为20A,K点为2*2*1,优化前后的初末态结构都很合理,但是计算出的吸附能为-48.4185 eV,实在想不通为什么会出现这 ...

回复 Reply

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

是的,经卡开发发老师提醒后发现小分子计算的INCAR出了点问题,谢谢老师的指导!
Daisy_zju 发表于 Post on 2022-1-12 13:42:25
卡开发发 发表于 2022-1-12 09:04
这次的问题来源于只有OUTCAR-3的体系做了DFT-D3(zero)的修正(也就是IVDW=11),而其他的两个都没做,这 ...

谢谢老师的指导,修改之后吸附能正常了,麻烦老师了!
jiangning198511 发表于 Post on 2022-1-12 13:09:15
Daisy_zju 发表于 2022-1-11 17:20
老师您好,不好意思再麻烦您一下,我昨晚统一赝势之后重新计算,现在吸附结构的计算还没完成但是几近收敛 ...

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

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
Daisy_zju + 4 谢谢老师指导

查看全部评分 View all ratings

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

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

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

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
Daisy_zju + 4 谢谢老师指导

查看全部评分 View all ratings

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

统一用VASPKIT生成 应该就一致了
Daisy_zju 发表于 Post on 2022-1-10 21:02:04
卡开发发 发表于 2022-1-10 17:10
应该是楼上说的赝势问题,其他影响应该是次要的。道理上说,ALGO一般而言影响很小,因为只是求解本征值的算 ...

应该是的,我之后重算一遍,谢谢老师的指导和分享!
Daisy_zju 发表于 Post on 2022-1-10 21:00:16
jiangning198511 发表于 2022-1-10 16:53
Co的POTCAR为啥不一致,POTCAR:    PAW_PBE Co 02Aug2007;POTCAR:   PAW_PBE Co 06Sep2000?

谢谢老师指出的问题,POTCAR的不一致应该是因为一个是我用Vaspkit特地生成的,一个是生成KPOINTS的时候产生的,计算的时候也没有检查两者的一致性(万分没想到)。
卡开发发 发表于 Post on 2022-1-10 17:21:13
本帖最后由 卡开发发 于 2022-1-10 17:30 编辑
jiangning198511 发表于 2022-1-10 16:53
Co的POTCAR为啥不一致,POTCAR:    PAW_PBE Co 02Aug2007;POTCAR:   PAW_PBE Co 06Sep2000?

应该就是这个问题,要不然三个体系的EATOM差距不会那么大
卡开发发 发表于 Post on 2022-1-10 17:10:34
本帖最后由 卡开发发 于 2022-1-10 17:26 编辑
Daisy_zju 发表于 2022-1-10 16:20
十分抱歉老师,经您提醒才发现OUTCAR-1文件放错了(之前的那个是苯胺的OUTCAR),现在这个是正确的NIONS= ...
应该是楼上说的赝势问题,其他影响应该是次要的。道理上说,ALGO一般而言影响很小,因为只是求解本征值的算法不一样,不应该有特别显著的差异。SIGMA的差异则可以通过SCF最后一步的“entropy T*S    EENTRO”这个字段来得到SIGMA的影响,看起来应该是可以忽略的。

评分 Rate

参与人数
Participants 3
eV +7 收起 理由
Reason
978142355 + 2 正解
Daisy_zju + 2 谢谢
壹零壹室掃地僧 + 3 谢谢分享

查看全部评分 View all ratings

jiangning198511 发表于 Post on 2022-1-10 16:53:23
Daisy_zju 发表于 2022-1-10 16:20
十分抱歉老师,经您提醒才发现OUTCAR-1文件放错了(之前的那个是苯胺的OUTCAR),现在这个是正确的NIONS= ...

Co的POTCAR为啥不一致,POTCAR:    PAW_PBE Co 02Aug2007;POTCAR:   PAW_PBE Co 06Sep2000?

评分 Rate

参与人数
Participants 2
eV +4 收起 理由
Reason
Daisy_zju + 2 谢谢
卡开发发 + 2 我很赞同

查看全部评分 View all ratings

Daisy_zju 发表于 Post on 2022-1-10 16:20:00
卡开发发 发表于 2022-1-10 15:43
不是上面说的问题,三个OUTCAR里面给出的原子数(NIONS字段):
OUTCAR1:14
OUTCAR2:168

十分抱歉老师,经您提醒才发现OUTCAR-1文件放错了(之前的那个是苯胺的OUTCAR),现在这个是正确的NIONS=23。也存在同样的问题,即E(adsorption)=E(slab+mol)-E(slab)-E(mol)=-1275.1269
-(-1098.8605)-(
-129.0427 )= -47.2237 eV。另外,老师您的意思是说不是上面的ALGO和sigma的设置问题,但是这么低的吸附能也实在是太匪夷所思了

outcar.tar.gz

1.24 MB, 下载次数 Times of downloads: 4

卡开发发 发表于 Post on 2022-1-10 15:43:44
Daisy_zju 发表于 2022-1-10 15:12
好的老师,非常不好意思!

不是上面说的问题,三个OUTCAR里面给出的原子数(NIONS字段):
OUTCAR1:14
OUTCAR2:168
OUTCAR3:191
你可能需要仔细检查一下你的模型。
Daisy_zju 发表于 Post on 2022-1-10 15:12:14
sobereva 发表于 2022-1-10 15:04
此帖四处违规:

VASP是第一性原理程序,相关帖子不要发到量子化学板块,给你移动了,以后注意

好的老师,非常不好意思!

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-17 21:11 , Processed in 0.162211 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list