计算化学公社

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

[波函数分析求助] 做基于力场的能量分解时发现总能量差于高斯计算得到的结果不一致?请问结果能用吗?

[复制链接 Copy URL]

97

帖子

0

威望

531

eV
积分
628

Level 4 (黑子)

本帖最后由 wangxiaofeng 于 2022-3-29 18:11 编辑

图中单位都为kcal/mol,通过高斯得到的电子能差为正值,仅为1.12,而基于力场的能量分解分析得到总能量差为负值,而且差值还极大,这种现象正常吗?

202203291801371800..png (17.53 KB, 下载次数 Times of downloads: 53)

202203291801371800..png

97

帖子

0

威望

531

eV
积分
628

Level 4 (黑子)

20#
 楼主 Author| 发表于 Post on 2022-4-2 10:45:06 | 只看该作者 Only view this author
sobereva 发表于 2022-4-1 23:58
那就显然属于隔三个键以上的情况

好的,明白。谢谢老师!

6万

帖子

99

威望

6万

eV
积分
125147

管理员

公社社长

19#
发表于 Post on 2022-4-1 23:58:41 | 只看该作者 Only view this author
wangxiaofeng 发表于 2022-4-1 15:50
但是对于两部分没有化学键相连相连的体系,怎么判断它们之间有没有间隔三个化学键呢?

那就显然属于隔三个键以上的情况
北京科音自然科学研究中心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

97

帖子

0

威望

531

eV
积分
628

Level 4 (黑子)

18#
 楼主 Author| 发表于 Post on 2022-4-1 15:50:51 | 只看该作者 Only view this author
sobereva 发表于 2022-4-1 12:00
你的图中黄色区域和另一个部分有一个化学键相连,显然不属于EDA-FF适用范畴

EDA-FF是靠分子力场的非键 ...

但是对于两部分没有化学键相连相连的体系,怎么判断它们之间有没有间隔三个化学键呢?

6万

帖子

99

威望

6万

eV
积分
125147

管理员

公社社长

17#
发表于 Post on 2022-4-1 12:00:04 | 只看该作者 Only view this author
wangxiaofeng 发表于 2022-4-1 11:29
好的,谢谢。您的意思我理解了。但是我对“片段间至少要隔三个化学键(相对于键连关系而言,而不是原子间 ...

你的图中黄色区域和另一个部分有一个化学键相连,显然不属于EDA-FF适用范畴

EDA-FF是靠分子力场的非键相互作用项描述片段间相互作用,如果隔<=3个键,那两个片段间的作用就不是光靠非键相互作用描述了。
北京科音自然科学研究中心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

97

帖子

0

威望

531

eV
积分
628

Level 4 (黑子)

16#
 楼主 Author| 发表于 Post on 2022-4-1 11:29:27 | 只看该作者 Only view this author
sobereva 发表于 2022-4-1 00:53
这种情况明显不适合EDA-FF
使用Multiwfn做基于分子力场的能量分解分析
http://sobereva.com/442(http: ...

好的,谢谢。您的意思我理解了。但是我对“片段间至少要隔三个化学键(相对于键连关系而言,而不是原子间距而言)”这句话不是很理解,能麻烦您具体解释一下吗?

6万

帖子

99

威望

6万

eV
积分
125147

管理员

公社社长

15#
发表于 Post on 2022-4-1 00:53:42 | 只看该作者 Only view this author
wangxiaofeng 发表于 2022-3-31 15:41
老师,您好!体系分为黄色选中部分和另外未标注一的部分,麻烦您了

这种情况明显不适合EDA-FF
使用Multiwfn做基于分子力场的能量分解分析
http://sobereva.com/442http://bbs.keinsci.com/thread-10907-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

97

帖子

0

威望

531

eV
积分
628

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2022-3-31 15:41:02 | 只看该作者 Only view this author
本帖最后由 wangxiaofeng 于 2022-3-31 15:42 编辑
sobereva 发表于 2022-3-31 14:58
直接把你算的体系的截图发出来,在上面标注片段是怎么划分的,一看就明白了
光靠猜没法给你准确的回复

老师,您好!体系分为黄色选中部分和另外未标注一的部分,麻烦您了

202203311541254058..png (91.83 KB, 下载次数 Times of downloads: 49)

202203311541254058..png

6万

帖子

99

威望

6万

eV
积分
125147

管理员

公社社长

13#
发表于 Post on 2022-3-31 14:58:00 | 只看该作者 Only view this author
直接把你算的体系的截图发出来,在上面标注片段是怎么划分的,一看就明白了
光靠猜没法给你准确的回复
北京科音自然科学研究中心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

97

帖子

0

威望

531

eV
积分
628

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2022-3-31 11:52:17 | 只看该作者 Only view this author
本帖最后由 wangxiaofeng 于 2022-3-31 12:01 编辑
agent99 发表于 2022-3-31 02:51
有过渡金属,还是过渡态,普通的力场怎么可能可靠?除非你的力场是专门针对这类体系开发的

请问什么力场适合过渡金属的有机体系呢?目前Multiwfn只有这两种力场可以选择 1.UFF van der Waals model 2: AMBER99 & GAFF van der Waals model

269

帖子

1

威望

3030

eV
积分
3319

Level 5 (御坂)

11#
发表于 Post on 2022-3-31 02:51:00 | 只看该作者 Only view this author
有过渡金属,还是过渡态,普通的力场怎么可能可靠?除非你的力场是专门针对这类体系开发的

108

帖子

0

威望

3946

eV
积分
4054

Level 6 (一方通行)

10#
发表于 Post on 2022-3-30 21:57:01 | 只看该作者 Only view this author
wangxiaofeng 发表于 2022-3-30 21:36
请问是需要提供分子结构吗?还有一个问题EDA-FF分析得到的数据其中total这项数据说的是两个片段间的相互 ...

想问的是你想算什么问题?

97

帖子

0

威望

531

eV
积分
628

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2022-3-30 21:36:12 | 只看该作者 Only view this author
Kangtor 发表于 2022-3-30 18:49
实际上,你到底算了什么内容是描述的不够的,没法给你具体的建议。卢老师也提出了这一点。

我也是猜测 ...

请问是需要提供分子结构吗?还有一个问题EDA-FF分析得到的数据其中total这项数据说的是两个片段间的相互作用的能量而不是整个结构式的能量吧?

202203302136267837..png (45.75 KB, 下载次数 Times of downloads: 54)

202203302136267837..png

108

帖子

0

威望

3946

eV
积分
4054

Level 6 (一方通行)

8#
发表于 Post on 2022-3-30 18:49:15 | 只看该作者 Only view this author
wangxiaofeng 发表于 2022-3-30 17:32
请问键接项能量从哪里得到?

实际上,你到底算了什么内容是描述的不够的,没法给你具体的建议。卢老师也提出了这一点。

我也是猜测一下

97

帖子

0

威望

531

eV
积分
628

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2022-3-30 17:32:30 | 只看该作者 Only view this author
Kangtor 发表于 2022-3-30 16:40
你加上键接项的能量再来比较

力场开发是在一定特定的化学环境的,和量化计算结果有一定差异很正常

请问键接项能量从哪里得到?

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

GMT+8, 2026-2-22 14:34 , Processed in 0.179918 second(s), 24 queries , Gzip On.

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