计算化学公社

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

[VASP] 请问AIMD模拟电催化CO2RR,CO2总是游离在溶剂中怎么处理?

[复制链接 Copy URL]

8

帖子

0

威望

199

eV
积分
207

Level 3 能力者

各位老师好,我最近在使用AIMD模拟电催化CO2RR时,CO2一直游离在溶剂中,不吸附在表面。文章中的机理如图所示,
我尝试过很多初始构型,分别控制CO2与界面、金属阳离子的距离,但保证CO2最初都是吸附在slab上的,可AIMD模拟后期CO2均游离在溶剂中,不与金属阳离子接触,也不与H2O分子反应。按理说我的初始构型是机理图中的第二步,那么接下来应该生成O-C-O-H,可一直游离在溶剂中。
我的理解是:机理图中,CO2从游离状态到吸附状态是得到电子的,我没有给CO2电子,单纯的将CO2放置在slab上是不可行的,所以不能吸附在表面?如果需要给予CO2电子的话,具体应当如何操作呢?
请各位老师指导~




8

帖子

0

威望

199

eV
积分
207

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2022-11-9 14:13:04 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-9 05:30
AIMD只能模拟ps时间尺度的反应,靠硬堆计算资源可以做到ns级,但也就仅此而已了。
你再想想CO2RR的TOF一 ...

十分感谢老师~
CO2反应所需的时间的确远远大于AIMD所能模拟的时间。
我现在重复的是Koper教授在Nat Catal发表的一篇文章,https://doi.org/10.1038/s41929-021-00655-5 文章中用的就是2ps,我猜想他们是不是利用了一些特殊的方法让反应在2ps内完成。后期还没办法我就换个方向。
谢谢老师~

8

帖子

0

威望

199

eV
积分
207

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2022-11-9 13:44:15 | 只看该作者 Only view this author
asd9998876 发表于 2022-11-8 22:27
你加个金属阳离子体系不是有净电荷了吗,有没有加抗衡离子

谢谢老师~
我是加的抗衡离子是OH

1万

帖子

0

威望

9859

eV
积分
22096

Level 6 (一方通行)

5#
发表于 Post on 2022-11-9 05:30:54 | 只看该作者 Only view this author
新建文件夹 发表于 2022-11-8 08:32
老师您好,非常感谢您的回复~
请问有什么办法可以增强CO2在催化剂表面的吸附吗?无论我怎么调节初始构型 ...

AIMD只能模拟ps时间尺度的反应,靠硬堆计算资源可以做到ns级,但也就仅此而已了。
你再想想CO2RR的TOF一般是多少,我没有具体概念,但肯定远不到10^9 s^-1对吧?
所以你跑好几个月的AIMD,连一个CO2分子都没反应,是极其正常的,因为现实中反应就是这么慢。这种情况下说明AIMD不适合研究这个体系,需要改用过渡态理论来研究。
做AIMD的一个常见误区就是不预先估计反应的时间尺度,明明需要至少毫秒级的时间才能发生的反应,却用只能研究纳秒以内时间尺度的方法来计算,那当然不反应。

评分 Rate

参与人数
Participants 4
eV +7 收起 理由
Reason
ylbt1989 + 2 我很赞同
pikapikachu + 2 我很赞同
银河牛仔 + 1 我很赞同
新建文件夹 + 2 谢谢

查看全部评分 View all ratings

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
本团队长期招收研究生,有意者可私信联系

42

帖子

0

威望

928

eV
积分
970

Level 4 (黑子)

4#
发表于 Post on 2022-11-8 22:27:06 | 只看该作者 Only view this author
你加个金属阳离子体系不是有净电荷了吗,有没有加抗衡离子

8

帖子

0

威望

199

eV
积分
207

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2022-11-8 15:32:58 | 只看该作者 Only view this author
Eudaimonia 发表于 2022-11-8 09:10
吸附这个过程应该是不需要电荷参与的

你这个过程很像是CO2在催化剂表面吸附太弱了/竞争不过溶剂,导致一 ...

老师您好,非常感谢您的回复~
请问有什么办法可以增强CO2在催化剂表面的吸附吗?无论我怎么调节初始构型最后都是游离在溶剂中,困扰接近一个月了
添加金属阳离子是为了探究不同阳离子对CO2RR的影响;
您最后一句我不太理解,通过SCF计算电荷的位置和电子落在CO2上有什么联系吗?计算出最低能量轨道,vasp可以设置电子位置吗?
麻烦您了~

210

帖子

3

威望

2110

eV
积分
2380

Level 5 (御坂)

2#
发表于 Post on 2022-11-8 09:10:18 | 只看该作者 Only view this author
吸附这个过程应该是不需要电荷参与的

你这个过程很像是CO2在催化剂表面吸附太弱了/竞争不过溶剂,导致一定时间AIMD后脱离

另外你的金属阳离子加进去的原因是什么呢?有没有设置恰当的电荷,单纯想控制电子落在CO2上比较困难,一般来说都是靠SCF计算来确定电荷在哪个能量最低的轨道上

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
新建文件夹 + 2 谢谢

查看全部评分 View all ratings

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

GMT+8, 2026-2-21 01:32 , Processed in 0.177088 second(s), 25 queries , Gzip On.

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