计算化学公社

 找回密码 Forget password
 注册 Register
Views: 8506|回复 Reply: 17
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] 计算吸附过程发现物理吸附吸附能比化学吸附更大,这可能吗?

[复制链接 Copy URL]

30

帖子

0

威望

213

eV
积分
243

Level 3 能力者

本帖最后由 xilemiki 于 2022-3-17 22:51 编辑

各位老师好,我在计算糠醇在Ir2-Gra表面吸附的过程中出现了物理吸附比化学吸附吸附能更大的现象,按照常识来说,不应该是化学吸附吸附能更大吗(物理吸附dE=-0.42,化学吸附dE=0.03,两者DFT能量相差0.45)?
从计算完成后的arc图可以看出来图1C和Ir之间确实成键了,图2也确实为物理吸附,实在不明白为什么会出现这种现象,希望各位老师帮忙解答一下,谢谢!

1.png (50.28 KB, 下载次数 Times of downloads: 32)

图1

图1

2.png (31.45 KB, 下载次数 Times of downloads: 37)

2.png

INCAR

1.01 KB, 下载次数 Times of downloads: 9

254

帖子

0

威望

1143

eV
积分
1397

Level 4 (黑子)

2#
发表于 Post on 2022-3-18 01:09:58 | 只看该作者 Only view this author
IDIPOL = 3,打开了偶极之后,计算的能量就要十分小心(详细可以查看电子部收敛时的能量和该步最后总结能量的差别),偶极校正一般没有特别情况不需要打开,结果有时候会难以想象。
努力挖掘科学的本质

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

3#
发表于 Post on 2022-3-18 07:16:43 | 只看该作者 Only view this author
这种事明显说明自己的计算有问题,没有什么物理上的原因
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

30

帖子

0

威望

213

eV
积分
243

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2022-3-18 10:18:22 | 只看该作者 Only view this author
dizzy:差值为2.6eV 这合理吗

3.png (15.68 KB, 下载次数 Times of downloads: 34)

3.png

30

帖子

0

威望

213

eV
积分
243

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2022-3-18 10:21:21 | 只看该作者 Only view this author
sobereva 发表于 2022-3-18 07:16
这种事明显说明自己的计算有问题,没有什么物理上的原因

谢谢社长 我也觉得我有问题 我再找找原因

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

6#
发表于 Post on 2022-3-18 15:09:36 | 只看该作者 Only view this author
没有化学吸附能一定比物理吸附能更大这个规则,只是大部分情况更大,而且对于课本上的简单例子都是更大而已,课本不会教给你例外情况。
你可以这么想:有的体系有化学吸附,化学吸附能比物理吸附能大;有的体系没有化学吸附,如果强行摆成化学吸附的结构,吸附质会因为能量太高而自发远离表面,回到物理吸附的结构。那么如果有一个体系介于这两种体系之间,就有可能化学吸附结构虽然存在但是比物理吸附结构的能量高。具体到你的体系,因为化学吸附会破坏呋喃环的芳香性,所以导致化学吸附不稳定,但是没有不稳定到完全不存在的程度。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

69

帖子

1

威望

982

eV
积分
1071

Level 4 (黑子)

7#
发表于 Post on 2022-3-18 15:36:04 | 只看该作者 Only view this author
大概率是你吸附的基底结构变了,这部分能量在你计算化学吸附的时候被忽视了。
信息不灭,熵永增

30

帖子

0

威望

213

eV
积分
243

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2022-3-18 16:28:32 | 只看该作者 Only view this author
wzkchem5 发表于 2022-3-18 15:09
没有化学吸附能一定比物理吸附能更大这个规则,只是大部分情况更大,而且对于课本上的简单例子都是更大而已 ...

谢谢解答,但是为什么在其它同类型的表面上没有出现这种现象呢(掺杂Ru,Co等过渡金属)?在其他过渡金属的表面上化学吸附吸附能是比物理吸附吸附能大的

30

帖子

0

威望

213

eV
积分
243

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2022-3-18 16:32:09 | 只看该作者 Only view this author
EdwardLimit 发表于 2022-3-18 15:36
大概率是你吸附的基底结构变了,这部分能量在你计算化学吸附的时候被忽视了。

谢谢解答,吸附上糠醇后基底结构确实会发生变化,但是计算吸附能的时候不就是dE=Eads-Esur-Egas吗?其中的Esur指的难道不是结构优化后的表面能量吗?

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

10#
发表于 Post on 2022-3-18 18:16:04 | 只看该作者 Only view this author
xilemiki 发表于 2022-3-18 09:28
谢谢解答,但是为什么在其它同类型的表面上没有出现这种现象呢(掺杂Ru,Co等过渡金属)?在其他过渡金属 ...

首先,一个分子在一个基底上的化学吸附结构未必是唯一的(还是那句话,课本里的简单情况都是唯一的,课本不会给你教不唯一的情形),可能对于Ru、Co也存在能量高的化学吸附态,反过来对于Ir可能也存在能量低的化学吸附态,必须对每个金属都把所有可能的化学吸附构型找全了才能讨论。
其次,如果已经确定把所有化学吸附构型都找全了,结论仍然如此的话,那一个可能的因素是Ir的原子尺寸,或者说原子轨道尺寸的问题,可能这个化学吸附要求金属-金属距离比较近,或者要求金属-碳键比较短,否则不容易形成比较强的化学吸附
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

69

帖子

1

威望

982

eV
积分
1071

Level 4 (黑子)

11#
发表于 Post on 2022-3-18 22:19:07 | 只看该作者 Only view this author
xilemiki 发表于 2022-3-18 16:32
谢谢解答,吸附上糠醇后基底结构确实会发生变化,但是计算吸附能的时候不就是dE=Eads-Esur-Egas吗?其中 ...

你算基底单点能Esur的时候取得是最开始的基底还是变化后的基底,要是变化后的基底当我没说。但就看你那两张图,不可能前者是正的还比后者大0.45 ev.
信息不灭,熵永增

30

帖子

0

威望

213

eV
积分
243

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2022-3-19 23:27:27 | 只看该作者 Only view this author
EdwardLimit 发表于 2022-3-18 22:19
你算基底单点能Esur的时候取得是最开始的基底还是变化后的基底,要是变化后的基底当我没说。但就看你那两 ...

Esur是结构优化的基底的能量,吸附导致后会导致基底变形,基底变形的Esur不是基底的基态能量阿,应该不能这么算吧,而且后续吸附物变化基底也会发生不同程度的变形,那Esur的值岂不是一直在变..

30

帖子

0

威望

213

eV
积分
243

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2022-3-19 23:34:06 | 只看该作者 Only view this author
EdwardLimit 发表于 2022-3-18 22:19
你算基底单点能Esur的时候取得是最开始的基底还是变化后的基底,要是变化后的基底当我没说。但就看你那两 ...

我没有计算基底的单点能,使用的是结构优化的得到的基底的基态能量

147

帖子

0

威望

1592

eV
积分
1739

Level 5 (御坂)

14#
发表于 Post on 2022-3-21 09:46:29 | 只看该作者 Only view this author
xilemiki 发表于 2022-3-19 23:34
我没有计算基底的单点能,使用的是结构优化的得到的基底的基态能量

我也刚刚遇到了,CO2一个化学吸附在底物上,co2已经明显被活化了;但是存在另外一种物理吸附,CO2距离底物大概2.5A,没有被活化,但两者吸附能算下来只相差了0.01eV。而且跟你的情况一样,化学吸附时,底物几乎没有变化,但物理吸附的时候,底物产生了轻微形变,能量升了0.5eV。。。
一蓑烟雨任平生

30

帖子

0

威望

213

eV
积分
243

Level 3 能力者

15#
 楼主 Author| 发表于 Post on 2022-3-24 09:38:18 | 只看该作者 Only view this author
ymeng 发表于 2022-3-21 09:46
我也刚刚遇到了,CO2一个化学吸附在底物上,co2已经明显被活化了;但是存在另外一种物理吸附,CO2距离底 ...

我们老板告诉我化学能吸附太小是因为这个吸附分子是我放的比较近然后成键,这种有点类似于人为施力,能量高一点是正常的,他以前也遇到过这种

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 19:57 , Processed in 0.180944 second(s), 23 queries , Gzip On.

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