计算化学公社

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

[GROMACS] 请教:固体表面吸附模拟成键判断

[复制链接 Copy URL]

29

帖子

0

威望

275

eV
积分
304

Level 3 能力者

本帖最后由 PHOENIXHMX 于 2023-5-26 10:29 编辑

想要请教各位老师,我在gromacs模拟固体表面水分子吸附的模拟,我该如何判断水分子是否与表面分子成键呢?用RDF或是其他方法?有没有软件能够显示成键情况?感谢各位老师解答~

62

帖子

0

威望

1663

eV
积分
1725

Level 5 (御坂)

2#
发表于 Post on 2023-5-26 11:50:50 | 只看该作者 Only view this author
Gromacs应该不能模拟成键吧,形成新的化学键说明有化学反应进行,这种情况如果做分子动力学的话,需要用到反应力场。
gromacs没有反应力场,而lammps有。如果体系不大的话,也可以尝试用CP2K做AIMD。

146

帖子

0

威望

1597

eV
积分
1743

Level 5 (御坂)

3#
发表于 Post on 2023-5-26 15:23:42 | 只看该作者 Only view this author
吸附水分子的话跟成键应该没啥关系吧,正常情况应该都是弱相互作用啊,难道你想说的是氢键?
由衷感谢每位帮助我的好心人

29

帖子

0

威望

275

eV
积分
304

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2023-5-26 16:13:19 | 只看该作者 Only view this author
本帖最后由 PHOENIXHMX 于 2023-5-26 16:17 编辑
dzdhp 发表于 2023-5-26 15:23
吸附水分子的话跟成键应该没啥关系吧,正常情况应该都是弱相互作用啊,难道你想说的是氢键?

因为我在看相关的表面吸附文献时,会产生Mg-O键和F-H键,就是图上这种,所以我就想在gromacs模拟中能否判断这种键的出现。

202305261613135439..png (102.24 KB, 下载次数 Times of downloads: 2)

202305261613135439..png

146

帖子

0

威望

1597

eV
积分
1743

Level 5 (御坂)

5#
发表于 Post on 2023-5-26 16:34:10 | 只看该作者 Only view this author
PHOENIXHMX 发表于 2023-5-26 16:13
因为我在看相关的表面吸附文献时,会产生Mg-O键和F-H键,就是图上这种,所以我就想在gromacs模拟中能否判 ...

文章发来看看,
由衷感谢每位帮助我的好心人

6万

帖子

99

威望

6万

eV
积分
125152

管理员

公社社长

6#
发表于 Post on 2023-5-27 05:38:08 | 只看该作者 Only view this author
PHOENIXHMX 发表于 2023-5-26 16:13
因为我在看相关的表面吸附文献时,会产生Mg-O键和F-H键,就是图上这种,所以我就想在gromacs模拟中能否判 ...

仅当表面吸附作用光靠原子电荷描述的静电相互作用结合范德华作用就能够合理描述时,靠GROMACS基于非反应性的经典力场的模拟才说得过去。否则用诸如CP2K跑DFT级别的第一性原理动力学研究才是妥当的。是哪种情况必须清楚具体交代当前研究的体系。

至于统计,rdf或者自己写VMD脚本都可以,对于展现特定结构的表面成键作用还可以用IGMH或IRI来考察
使用IRI方法图形化考察化学体系中的化学键和弱相互作用
http://sobereva.com/598http://bbs.keinsci.com/thread-23457-1-1.html
使用Multiwfn做IGMH分析非常清晰直观地展现化学体系中的相互作用
http://sobereva.com/621http://bbs.keinsci.com/thread-28147-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

29

帖子

0

威望

275

eV
积分
304

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2023-5-27 08:12:19 | 只看该作者 Only view this author
sobereva 发表于 2023-5-27 05:38
仅当表面吸附作用光靠原子电荷描述的静电相互作用结合范德华作用就能够合理描述时,靠GROMACS基于非反应 ...

明白了,谢谢sob老师

29

帖子

0

威望

275

eV
积分
304

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2023-5-27 08:58:48 | 只看该作者 Only view this author
sobereva 发表于 2023-5-27 05:38
仅当表面吸附作用光靠原子电荷描述的静电相互作用结合范德华作用就能够合理描述时,靠GROMACS基于非反应 ...

想要再请教老师一个问题,我现在模拟氟化钙(CaF2)表面水分子的运动分布情况,查看水分子的分层,用哪种力场模拟比较合适?

6万

帖子

99

威望

6万

eV
积分
125152

管理员

公社社长

9#
发表于 Post on 2023-5-28 01:57:39 | 只看该作者 Only view this author
PHOENIXHMX 发表于 2023-5-27 08:58
想要再请教老师一个问题,我现在模拟氟化钙(CaF2)表面水分子的运动分布情况,查看水分子的分层,用哪种 ...

最好是用CP2K做AIMD
用力场也不是不行(比如UFF的LJ参数结合REPEAT方法算的CaF2的slab体系的原子电荷,并结合OPC水模型),但最好有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

29

帖子

0

威望

275

eV
积分
304

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2023-5-29 11:06:23 | 只看该作者 Only view this author
sobereva 发表于 2023-5-28 01:57
最好是用CP2K做AIMD
用力场也不是不行(比如UFF的LJ参数结合REPEAT方法算的CaF2的slab体系的原子电荷, ...

sob老师好,因为我没做过AIMD,可以请教一下具体的流程要怎么做或者是有没有例子提供参考,谢谢了。

6万

帖子

99

威望

6万

eV
积分
125152

管理员

公社社长

11#
发表于 Post on 2023-5-29 12:28:28 | 只看该作者 Only view this author
PHOENIXHMX 发表于 2023-5-29 11:06
sob老师好,因为我没做过AIMD,可以请教一下具体的流程要怎么做或者是有没有例子提供参考,谢谢了。

几句话不可能说清楚,使用CP2K需要了解的知识相当多

北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/workshop/KFP_content.html)里面有非常详细的讲解,在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

255

帖子

1

威望

2682

eV
积分
2957

Level 5 (御坂)

12#
发表于 Post on 2023-5-29 12:40:41 | 只看该作者 Only view this author
PHOENIXHMX 发表于 2023-5-26 16:13
因为我在看相关的表面吸附文献时,会产生Mg-O键和F-H键,就是图上这种,所以我就想在gromacs模拟中能否判 ...

看贴的图片涉及到水分子断键了,这个gromacs实现不了。要么找水分子在材料表面裂解的过渡态,要么跑AIMD。如果不是大量水分子,前者或许更好一些。
我唯一知道的就是我一无所知,但我是化学小迷弟

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 13:53 , Processed in 0.177498 second(s), 23 queries , Gzip On.

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