计算化学公社

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

[波函数分析求助] EDA-FF结果数据怎么才能和VMD原子着色图对应起来

[复制链接 Copy URL]

39

帖子

0

威望

701

eV
积分
740

Level 4 (黑子)

本帖最后由 zhutuziqin 于 2019-12-14 18:22 编辑

各位老师大家好 本人在用Sobereva老师帖子 http://sobereva.com/442  对五个研究对象的七分子团簇结构进行弱相互作用分析
然后在分析的时候 遇到如下问题:

如表格一为各研究对象(用12345表示) 7个片段总的静电互斥色散作用能(单位KJ/mol)等   
  


情况1:其中2和3研究对象的 各片段间总的静电作用比色散作用大
然后氢键(键长都设为3.5埃 键角都设为30度)的成键的片段能够和片段间静电作用由大到小排序能够对应起来
比如下图的研究对象C5H4N4O 有8条氢键(标红的为氢键成键的片段) 然后能依次对应于6个静电作用最高的片段(包括两个片段间有2条氢键的情况)





情况2:如表格一中 1,4和5研究对象总的色散作用大比静电作用大 然后氢键成键的片段 并不能和静电作用由大到小的片段一一对应起来
比如说 下图的研究对象C5H4N4   氢键有5条(标红的为氢键成键的片段)但是它对应的成键的片段的静电相互作用的数据有些并不是最大的  甚至静电作用大的并没有成键




针对以上情况 有2个问题:第1:色散作用比静电作用大的体系 氢键成键的片段不能和静电相互作用的大小一一对应起来吗?  如果确实不能对应的话 该用哪一方面的知识解释这种情况呢?
第2:我检查了所有分析的步骤 原子电荷类型和电荷等 确保没有问题 在分析以上五个研究对象的时候也确保其他外在条件一致  得出以上表一的结果是正确的吗?
  希望老师们提出宝贵意见 衷心感谢



6万

帖子

99

威望

5万

eV
积分
120205

管理员

公社社长

2#
发表于 Post on 2019-12-14 11:01:24 | 只看该作者 Only view this author
注意恰当用逗号和句号。

1 没有你的体系具体结构说不清楚,也不清楚情况1和情况2到底对应什么情况。两个情况都给个复合物结构图,一看就清楚了
直接从结构上就能判断色散作用是否能起到重要作用

2 我不清楚你的表1是怎么算的
北京科音自然科学研究中心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

39

帖子

0

威望

701

eV
积分
740

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2019-12-14 22:58:19 | 只看该作者 Only view this author
sobereva 发表于 2019-12-14 11:01
注意恰当用逗号和句号。

1 没有你的体系具体结构说不清楚,也不清楚情况1和情况2到底对应什么情况。两个 ...

1.sobereva老师我更新了帖子插入了具体情况1,2的图片,麻烦您再看看具体的情况。

2.表1的数据是这样计算得来的:比如说表1的研究对象1的静电相互作用(-273.69KJ/mol)是将各片段间的静电相互作用加起来的(f1-f2+f1-f3+...+f6-f7)。

6万

帖子

99

威望

5万

eV
积分
120205

管理员

公社社长

4#
发表于 Post on 2019-12-15 12:30:46 | 只看该作者 Only view this author
VMD里用30度阈值太严了。VMD的键角定义和一般氢键键角定义不同,平时说的30度在VMD里应当对应43度

“研究对象”指什么?7聚体的不同构型?那两个图是其中两种构型?
北京科音自然科学研究中心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

39

帖子

0

威望

701

eV
积分
740

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2019-12-15 13:00:13 | 只看该作者 Only view this author
sobereva 发表于 2019-12-15 12:30
VMD里用30度阈值太严了。VMD的键角定义和一般氢键键角定义不同,平时说的30度在VMD里应当对应43度

“研 ...

“研究对象”就是五钟不同的物质,情况1具体的研究对象是次黄嘌呤,情况2的研究对象是指嘌呤。每种研究对象只有一个唯一的构型(是从晶体结构里面抠取的,来自于您的帖子:http://bbs.keinsci.com/thread-11390-1-1.html)。

6万

帖子

99

威望

5万

eV
积分
120205

管理员

公社社长

6#
发表于 Post on 2019-12-16 12:22:28 | 只看该作者 Only view this author
不能抠出来团簇后就直接用,因为晶体是测不准氢的位置的,至少也得用半经验方法(能用DFT当然更好)把氢的位置重新优化后再去做这个分析。

情况2 当中Frag 1 -- Frag 7静电作用才那么小,这是最可疑的一项,应当仔细看这两个结构、考察两个片段间原子之间的静电相互作用,弄清楚到底是怎么回事。这个解决了,其它的应该也就自然解决了。
北京科音自然科学研究中心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

39

帖子

0

威望

701

eV
积分
740

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2019-12-17 10:37:20 | 只看该作者 Only view this author
sobereva 发表于 2019-12-16 12:22
不能抠出来团簇后就直接用,因为晶体是测不准氢的位置的,至少也得用半经验方法(能用DFT当然更好)把氢的 ...

我扣出来的团簇结构有做全几何优化,拿优化后的结构直接做的EDA-FF分析。Sobereva老师您的意思是把抠出来的团簇结构,只对氢原子做限制性优化(优化方法我看了您的帖子:http://bbs.keinsci.com/thread-9022-1-1.html),然后再做EDA-FF分析吗?

您说的单独再考察两个片段间原子之间的静电相互作用,我下来之后再试一试,就是这方面理论知识太欠缺,深入一点的东西根本不知道该从哪儿下手该怎么解释,自己找了很多资料也没有什么可用有效的,老师您有推荐的学习资料吗??

6万

帖子

99

威望

5万

eV
积分
120205

管理员

公社社长

8#
发表于 Post on 2019-12-19 18:17:26 | 只看该作者 Only view this author
zhutuziqin 发表于 2019-12-17 10:37
我扣出来的团簇结构有做全几何优化,拿优化后的结构直接做的EDA-FF分析。Sobereva老师您的意思是把抠出来 ...

抠出来重新对所有原子优化,结构相对于晶体状态可能会发生很大变化。如果你的目的是分析晶体状态中的单体间相互作用,就不能对所有原子重新优化。

EDA-FF里静电作用计算非常简单,就是中学学过的点电荷间作用的库仑公式。片段间总静电作用就是两个片段间每一对原子静电作用的加和,因此通过考察关键原子间的相互作用就能帮助你了解片段间相互作用的细节,这没有任何难理解的。
北京科音自然科学研究中心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

39

帖子

0

威望

701

eV
积分
740

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2019-12-21 12:27:30 | 只看该作者 Only view this author
sobereva 发表于 2019-12-19 18:17
抠出来重新对所有原子优化,结构相对于晶体状态可能会发生很大变化。如果你的目的是分析晶体状态中的单体 ...

谢谢老师的指导,我明白您的意思了,我下来先思考一下自己EDA-FF分析的目的到底是什么,然后再作调整。

本版积分规则 Credits rule

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

GMT+8, 2025-8-18 07:08 , Processed in 0.162254 second(s), 23 queries , Gzip On.

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