计算化学公社

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

[综合交流] 利用径向分布函数(RDF)可以计算小分子在固体表面吸附的自由能吗?

[复制链接 Copy URL]

47

帖子

0

威望

284

eV
积分
331

Level 3 能力者

通过分子模拟计算自由能的方法有很多种,但是其中操作最方便的好像是利用径向分布函数(RDF)计算?


使用RDF计算小分子在固体表面吸附过程的自由能时可行的吗?因为注意到有的帖子里说计算化学反应的自由能时尽量使用增强采样的方法,但是计算吸附过程的自由能的话,是不是直接用RDF计算就好?与伞状采样等增强采样计算自由能的结果存在误差吗?

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

2#
发表于 Post on 2023-7-7 10:49:34 | 只看该作者 Only view this author
不要泛化rdf的用处
rdf能讨论的问题极其有限,然而却经常被胡用乱用
先想清楚表面吸附对应的路径是什么,比如考察的坐标是垂直于表面方向的,显然这和rdf所表现的粒子间径向坐标根本不对应,应当用垂直于表面方向的密度分布来折算PMF
北京科音自然科学研究中心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

47

帖子

0

威望

284

eV
积分
331

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-7-7 14:42:28 | 只看该作者 Only view this author
sobereva 发表于 2023-7-7 10:49
不要泛化rdf的用处
rdf能讨论的问题极其有限,然而却经常被胡用乱用
先想清楚表面吸附对应的路径是什么, ...

谢谢站长提醒!利用密度分布来折算PMF的公式或者相关文章,您能给个链接吗?之前一直没有看到这方面的文章

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

4#
发表于 Post on 2023-7-8 09:21:16 | 只看该作者 Only view this author
sigusigu 发表于 2023-7-7 14:42
谢谢站长提醒!利用密度分布来折算PMF的公式或者相关文章,您能给个链接吗?之前一直没有看到这方面的文 ...

参考
浅谈PCA与g_covar+g_anaeig+ddtdp+sigmaplot做自由能面图的方法
http://sobereva.com/73

另参考北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/workshop/KFP_content.html)这页ppt,其中g(r)换成垂直于表面的密度分布
北京科音自然科学研究中心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

47

帖子

0

威望

284

eV
积分
331

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-7-11 14:04:26 | 只看该作者 Only view this author
sobereva 发表于 2023-7-8 09:21
参考
浅谈PCA与g_covar+g_anaeig+ddtdp+sigmaplot做自由能面图的方法
http://sobereva.com/73

站长 感谢您的回复。 我是输出了水分子Z轴的数密度和质量密度,这个数值可以直接代入替换g(r)吗?总感觉这个数值不算概率。或者说将密度数值归一化处理一下?还是可以写成ρi/ρ0(前者为该点密度,后者为平均密度)?

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

6#
发表于 Post on 2023-7-11 14:30:52 | 只看该作者 Only view this author
sigusigu 发表于 2023-7-11 14:04
站长 感谢您的回复。 我是输出了水分子Z轴的数密度和质量密度,这个数值可以直接代入替换g(r)吗?总感 ...

可以
你感兴趣的是不同位置的相对自由能,即便你额外给密度加个分母项,在求差时也抵消了
北京科音自然科学研究中心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

47

帖子

0

威望

284

eV
积分
331

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2023-7-11 15:38:10 | 只看该作者 Only view this author
sobereva 发表于 2023-7-11 14:30
可以
你感兴趣的是不同位置的相对自由能,即便你额外给密度加个分母项,在求差时也抵消了

感谢站长 我去处理一下试试!

47

帖子

0

威望

284

eV
积分
331

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2023-7-28 16:02:08 | 只看该作者 Only view this author
sobereva 发表于 2023-7-11 14:30
可以
你感兴趣的是不同位置的相对自由能,即便你额外给密度加个分母项,在求差时也抵消了

站长,我最近跑了几个case,并且计算得到了PMF曲线。想请教一个问题:PMF既然是一个自由能随自变量变化的曲线(我这里就是Z轴坐标),那么我想计算的水分子从A点到B点的自由能变化,就是这两个点的PMF值作差吗?

6万

帖子

99

威望

5万

eV
积分
124694

管理员

公社社长

9#
发表于 Post on 2023-7-28 22:32:37 | 只看该作者 Only view this author
sigusigu 发表于 2023-7-28 16:02
站长,我最近跑了几个case,并且计算得到了PMF曲线。想请教一个问题:PMF既然是一个自由能随自变量变化的 ...

北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-1-25 07:32 , Processed in 0.165426 second(s), 23 queries , Gzip On.

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