计算化学公社

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

[CP2K] 胺分子负载在石墨烯表面的吸附剂,用AIMD可以描述化学反应吗?

[复制链接 Copy URL]

69

帖子

0

威望

270

eV
积分
339

Level 3 能力者

本帖最后由 qwe1832518773 于 2020-10-20 13:47 编辑

做了一个胺分子负载在石墨烯表面的吸附剂,然后在上方放了一些CO2分子,用cp2k做AIMD,跑了快1.8ps了,丝毫看不见CO2与胺有什么想要反应的预兆,甚至CO2分子都不往胺分子之间的空隙中扩散,这是为什么呢?体系总共242个原子,基组用的BASIS_MOLOPT,势函数用的GTH_POTENTIALS,温度600K,步长2fs。

0.2ps 1.7ps



1242

帖子

1

威望

3672

eV
积分
4934

Level 6 (一方通行)

2#
发表于 Post on 2020-10-20 11:54:10 | 只看该作者 Only view this author
本帖最后由 biogon 于 2020-10-20 14:55 编辑

能的,但是这个摆法显然不容易观察到有反应发生


236

帖子

1

威望

3805

eV
积分
4061

Level 6 (一方通行)

3#
发表于 Post on 2020-10-20 12:20:31 | 只看该作者 Only view this author
AIMD能做化学反应。可以尝试调节温度范围。或者增加反应时间

71

帖子

0

威望

1782

eV
积分
1853

Level 5 (御坂)

4#
发表于 Post on 2020-10-20 12:25:43 | 只看该作者 Only view this author
一种,CO2初始就放在胺分子周围;再一种,上面侧加一段时间压力差,使CO2有向下走的趋势,但是不知是否符合反应条件,是不是能轻松实现。

69

帖子

0

威望

270

eV
积分
339

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2020-10-20 12:38:20 | 只看该作者 Only view this author
随心恣意 发表于 2020-10-20 11:22
请问楼主解决这个问题了吗?

准备试试楼下两位的回复

69

帖子

0

威望

270

eV
积分
339

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2020-10-20 12:40:20 | 只看该作者 Only view this author
wolfli369 发表于 2020-10-20 12:20
AIMD能做化学反应。可以尝试调节温度范围。或者增加反应时间

感谢回复,我用600K是想让CO2分子扩散快一些,因为想看到CO2扩散进入胺分子之间空隙。但是目前看不出这个趋势,是因为我的胺分子间隔太小了的缘故吗?

69

帖子

0

威望

270

eV
积分
339

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2020-10-20 12:41:18 | 只看该作者 Only view this author
朙天儿 发表于 2020-10-20 12:25
一种,CO2初始就放在胺分子周围;再一种,上面侧加一段时间压力差,使CO2有向下走的趋势,但是不知是否符合 ...

感谢回复!打算先降一下温度试试看

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

8#
发表于 Post on 2020-10-20 13:19:22 | 只看该作者 Only view this author
biogon 发表于 2020-10-20 11:54
能的,但是这个摆法显然不容易观察到有反应发生

贴图方式不对,其他人看不到。请重新编辑帖子并阅读置顶的新社员必读贴了解怎么正确贴图,此问题在这里特意强调了:http://bbs.keinsci.com/thread-18961-1-1.html
北京科音自然科学研究中心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

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

9#
发表于 Post on 2020-10-20 13:25:49 | 只看该作者 Only view this author
VMD的图形窗口里按住鼠标左键拖动,松手的时候如果仍然有速度,体系就会一直旋转

AIMD最大价值之一就是能描述化学反应的动态过程,当然能描述化学反应

能否反应和原子间相对运动速度、相对朝向关系极大,不是很容易发生的反应(诸如自由基复合等)显然不是随便一摆就能自发跑出来的。多看看相关文章研究了解基本知识、模拟条件的设置。
北京科音自然科学研究中心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

553

帖子

0

威望

3324

eV
积分
3877

Level 5 (御坂)

10#
发表于 Post on 2020-10-20 13:27:00 | 只看该作者 Only view this author
常用的时间尺度下,brute-force MD只能跨越kbT level的能垒,所以靠裸跑MD看见化学反应本身就很难。

考虑用bias sampling方法吧。

71

帖子

0

威望

1782

eV
积分
1853

Level 5 (御坂)

11#
发表于 Post on 2020-10-20 13:36:03 | 只看该作者 Only view this author
qwe1832518773 发表于 2020-10-20 12:41
感谢回复!打算先降一下温度试试看

只降温恐怕没效果。

69

帖子

0

威望

270

eV
积分
339

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2020-10-20 13:37:42 | 只看该作者 Only view this author
sobereva 发表于 2020-10-20 13:25
VMD的图形窗口里按住鼠标左键拖动,松手的时候如果仍然有速度,体系就会一直旋转

AIMD最大价值之一就是 ...

感谢sob老师回复,基础知识确实有所欠缺,多请教一个问题,那降低温度(比如降到300K)是否会有利于观察到化学反应呢?
之前想的是:我既然不知道怎么摆放更好的话,我就索性给一个高温(600K),这样分子运动会更快,会更容易看到化学反应。
以上两种想法是孰优孰劣呢,存在什么问题呢

69

帖子

0

威望

270

eV
积分
339

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2020-10-20 13:38:19 | 只看该作者 Only view this author
k64_cc 发表于 2020-10-20 13:27
常用的时间尺度下,brute-force MD只能跨越kbT level的能垒,所以靠裸跑MD看见化学反应本身就很难。

考 ...

感谢回复,cp2k小白先记下了,暂时还不太了解您说的这俩要怎么实现

69

帖子

0

威望

270

eV
积分
339

Level 3 能力者

14#
 楼主 Author| 发表于 Post on 2020-10-20 13:44:52 | 只看该作者 Only view this author
朙天儿 发表于 2020-10-20 13:36
只降温恐怕没效果。

啊这样,但是光给CO2施加一个向下的力或者初速度要怎么实现呢

71

帖子

0

威望

1782

eV
积分
1853

Level 5 (御坂)

15#
发表于 Post on 2020-10-20 14:04:19 | 只看该作者 Only view this author
qwe1832518773 发表于 2020-10-20 13:44
啊这样,但是光给CO2施加一个向下的力或者初速度要怎么实现呢

大概知道反应温度压力吧,设置就尽量按照实际情况,而且初始摆放尽量靠近反应位点,不知道的位置话可以提出胺分子预测一下反应位点。不过这结构不太复杂,你的初始CO2向下平移,O贴近N多一些。
预测反应位点也可以成为工作的一部分。

本版积分规则 Credits rule

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

GMT+8, 2025-8-18 02:29 , Processed in 0.239856 second(s), 23 queries , Gzip On.

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