计算化学公社

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

[CP2K] AIMD是否平衡的分析及氢键数目计算

[复制链接 Copy URL]

30

帖子

0

威望

213

eV
积分
243

Level 3 能力者

本帖最后由 fly1992 于 2022-5-10 16:52 编辑

想请问大家几个关于aimd的问题,1.我用社长的Multiwfn生成的cp2k的inp文件,仿照文献中模拟Cu-H2O界面,模拟了20000fs,用VMD分析轨迹的RMSD,得到如下的结果:

我看了社长之前的帖子,对于判断不同的体系是否达到平衡,RMSD分析应该用不同的标准,像这种固液界面的体系应该用什么样的标准?另外,这里应该选择all还是选择每一种元素单独查看?
2.判断这种固液界面模拟是否达到平衡,除了社长所提到的温度,能量,轨迹,我看到文献中这种固液界面的模拟还有用费米能级和静电势来判断模拟是否达到平衡,这里所说的费米能级和静电势应该采用什么标准进行判断呢?
3.vmd中有一个功能是Hydrogen Bond,这个里面section1 和section2 指的是两个氢键的两个部分吗?如果我想得到水分子之前的氢键数目,这里是不是应该分别填我想查看区域的O?
这里是一篇文献(Nature material, vol 18, 2019, 697)中的图,这里所提到的Ndonor应该跟氢键数目是相同的意思吗?如果是一个意思的话,section1和section2分别应该填什么?


6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

2#
发表于 Post on 2022-5-10 17:51:30 | 只看该作者 Only view this author
1 体系密度、水之间的氢键数,Cu特定距离内O或H数目随时间的变化等,都可以用来判断平衡

2 费米能级和静电势跟判断平衡没明显关系。如果你看到什么文章这么做,看文章里具体怎么说的

3 都填water
北京科音自然科学研究中心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 能力者

3#
 楼主 Author| 发表于 Post on 2022-5-10 19:52:15 | 只看该作者 Only view this author
sobereva 发表于 2022-5-10 17:51
1 体系密度、水之间的氢键数,Cu特定距离内O或H数目随时间的变化等,都可以用来判断平衡

2 费米能级和静 ...

谢谢社长回复,那这种体系,RMSD大概在什么范围也可以当作体系平衡的判据?

30

帖子

0

威望

213

eV
积分
243

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2022-5-10 19:54:12 | 只看该作者 Only view this author
sobereva 发表于 2022-5-10 17:51
1 体系密度、水之间的氢键数,Cu特定距离内O或H数目随时间的变化等,都可以用来判断平衡

2 费米能级和静 ...

还有,社长,我刚试了您说的两个section都填water,没有氢键存在,那是不是填element H O也是可以的?

240

帖子

1

威望

4087

eV
积分
4347

Level 6 (一方通行)

5#
发表于 Post on 2022-5-11 15:19:17 | 只看该作者 Only view this author
fly1992 发表于 2022-5-10 19:54
还有,社长,我刚试了您说的两个section都填water,没有氢键存在,那是不是填element H O也是可以的?

请问田院士这篇Nature Materials模拟过程中,怎么调节体系的Potential?能给解答一下吗?

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

6#
发表于 Post on 2022-5-11 17:48:20 | 只看该作者 Only view this author
fly1992 发表于 2022-5-10 19:52
谢谢社长回复,那这种体系,RMSD大概在什么范围也可以当作体系平衡的判据?

我不知道你说的RMSD是指什么
看曲线是否在水平线上下均匀波动就完了
北京科音自然科学研究中心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

威望

6万

eV
积分
125127

管理员

公社社长

7#
发表于 Post on 2022-5-11 17:49:12 | 只看该作者 Only view this author
fly1992 发表于 2022-5-10 19:54
还有,社长,我刚试了您说的两个section都填water,没有氢键存在,那是不是填element H O也是可以的?

你先在representation界面里把drawing method设成H-bond,恰当设置氢键判断阈值,看有没有显示出氢键来
北京科音自然科学研究中心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 能力者

8#
 楼主 Author| 发表于 Post on 2022-5-16 13:33:29 | 只看该作者 Only view this author
wolfli369 发表于 2022-5-11 15:19
请问田院士这篇Nature Materials模拟过程中,怎么调节体系的Potential?能给解答一下吗?

在界面处放置不同数目的Na进行处理

240

帖子

1

威望

4087

eV
积分
4347

Level 6 (一方通行)

9#
发表于 Post on 2022-5-16 14:32:42 | 只看该作者 Only view this author
fly1992 发表于 2022-5-16 13:33
在界面处放置不同数目的Na进行处理

麻烦问一下,根据钠离子数目怎么计算电势呢

30

帖子

0

威望

213

eV
积分
243

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2022-5-16 15:45:35 | 只看该作者 Only view this author
wolfli369 发表于 2022-5-16 14:32
麻烦问一下,根据钠离子数目怎么计算电势呢

PRL 119, 016801 (2017),程俊老师这篇文章有讲

240

帖子

1

威望

4087

eV
积分
4347

Level 6 (一方通行)

11#
发表于 Post on 2022-5-16 16:34:07 | 只看该作者 Only view this author
fly1992 发表于 2022-5-16 15:45
PRL 119, 016801 (2017),程俊老师这篇文章有讲

收到,谢谢

8

帖子

0

威望

375

eV
积分
383

Level 3 能力者

12#
发表于 Post on 2023-11-7 00:48:49 | 只看该作者 Only view this author
sobereva 发表于 2022-5-10 17:51
1 体系密度、水之间的氢键数,Cu特定距离内O或H数目随时间的变化等,都可以用来判断平衡

2 费米能级和静 ...

想问社长有没有计算周期性体系的金属表面上方H2O中H、O密度的脚本?

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

13#
发表于 Post on 2023-11-7 06:28:29 | 只看该作者 Only view this author
moulihui 发表于 2023-11-7 00:48
想问社长有没有计算周期性体系的金属表面上方H2O中H、O密度的脚本?

可以直接用gromacs的gmx density进行统计
北京科音自然科学研究中心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

8

帖子

0

威望

375

eV
积分
383

Level 3 能力者

14#
发表于 Post on 2023-11-7 09:42:39 | 只看该作者 Only view this author
sobereva 发表于 2023-11-7 06:28
可以直接用gromacs的gmx density进行统计

非常感谢社长指点,我去找找教程

本版积分规则 Credits rule

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

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

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