计算化学公社

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

[综合交流] 静电能占总相互作用能的比例应该怎么算?

[复制链接 Copy URL]

103

帖子

0

威望

629

eV
积分
732

Level 4 (黑子)

各位老师好,我想计算不同离子液体中静电能占总能量的比例,以此来探究静电能对其性质的影响。计算的时候用静电能除以总能量来算发现根本不对,我看到一个表格,不知道它里面的比例是怎么算出来的,请问老师,这个是有什么不一样公式或算法吗?我暂时在网上还没有查到。。。。谢谢老师!

VF$GTO0LK29PY5B8M5I)__3.png (136.95 KB, 下载次数 Times of downloads: 22)

VF$GTO0LK29PY5B8M5I)__3.png

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

2#
发表于 Post on 2016-7-13 19:53:54 | 只看该作者 Only view this author
做SAPT能量分解,此文第10节提了
Multiwfn支持的弱相互作用的分析方法概览
http://sobereva.com/252

不是占总能量,而是占总相互作用能。已经把你的标题改了。
北京科音自然科学研究中心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

103

帖子

0

威望

629

eV
积分
732

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2016-7-14 10:32:19 | 只看该作者 Only view this author
sobereva 发表于 2016-7-13 19:53
做SAPT能量分解,此文第10节提了
Multiwfn支持的弱相互作用的分析方法概览
http://sobereva.com/252

谢谢老师,我现在是可以计算出各项能量了,但是不会算静电能占总相互作用能的比例。比如那个表格中第一个化合物的静电能是-1.2,总能量是-1.4,那么EElst(%)怎么会是13%呢?我查资料没有找到,发邮件问作者一直没回,自己试了一些的算法都不对。。。。所以想问问老师这是不是有什么特别的算法?谢谢老师!

239

帖子

0

威望

2340

eV
积分
2579

Level 5 (御坂)

4#
发表于 Post on 2016-7-14 12:55:59 | 只看该作者 Only view this author
浅原问问 发表于 2016-7-14 10:32
谢谢老师,我现在是可以计算出各项能量了,但是不会算静电能占总相互作用能的比例。比如那个表格中第一个 ...
原文Elst(%)是怎么定义的?
十室之邑,必有忠信如丘焉,不如丘之好学也。

103

帖子

0

威望

629

eV
积分
732

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2016-7-14 15:47:50 | 只看该作者 Only view this author
brothers 发表于 2016-7-14 12:55
原文Elst(%)是怎么定义的?

应该就是静电能占总能量的比例。。。

3L3]FFNP~80PQLPK830TO9G.png (45.79 KB, 下载次数 Times of downloads: 32)

3L3]FFNP~80PQLPK830TO9G.png

239

帖子

0

威望

2340

eV
积分
2579

Level 5 (御坂)

6#
发表于 Post on 2016-7-14 18:42:46 | 只看该作者 Only view this author
本帖最后由 brothers 于 2016-7-14 20:02 编辑
浅原问问 发表于 2016-7-14 15:47
应该就是静电能占总能量的比例。。。

搞不清楚,一个简单的比例按说很好算。
没用过sapt,但是感觉直接拿静电能除以总能量好像没什么意义啊,比较靠谱的是静电能在所有吸引作用中的比例?
表格里的IP1的Eelst%不是75%吗...正文里怎么又成了77%了,作者好粗心啊

太不靠谱了,第一行应该是把数据搞错了。其他各行都是 静电能/所有吸引作用%
十室之邑,必有忠信如丘焉,不如丘之好学也。

44

帖子

0

威望

2624

eV
积分
2668

Level 5 (御坂)

7#
发表于 Post on 2016-7-14 18:53:36 | 只看该作者 Only view this author
静电能在所有吸引作用中的比例,就是这么算的。
http://pymol.chenzhaoqiang.com/index.html

103

帖子

0

威望

629

eV
积分
732

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2016-7-14 20:53:18 | 只看该作者 Only view this author
brothers 发表于 2016-7-14 18:42
搞不清楚,一个简单的比例按说很好算。
没用过sapt,但是感觉直接拿静电能除以总能量好像没什么意义啊, ...

您好,我不是太明白您说的所有吸引作用是指什么?如果是总相互作用能的话,sapt中,把总相互作用能分成了四类,就是表格中的E Elst ,E Ind,E Disp,E Exch,这四类能量。这四类能量再加上误差能量δE int HF ,就是E SAPT了,这样理解的话,还是得不到表格中的比例啊。。。

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

9#
发表于 Post on 2016-7-14 21:44:42 | 只看该作者 Only view this author
浅原问问 发表于 2016-7-14 20:53
您好,我不是太明白您说的所有吸引作用是指什么?如果是总相互作用能的话,sapt中,把总相互作用能分成了 ...


你用相互作用能的某个分量值除以总相互作用能就完了
北京科音自然科学研究中心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

103

帖子

0

威望

629

eV
积分
732

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2016-7-15 09:11:07 | 只看该作者 Only view this author
sobereva 发表于 2016-7-14 21:44
你用相互作用能的某个分量值除以总相互作用能就完了

老师我本来也想这么算的,但是您看表格中第一行数据,总相互作用E SAPT是 -1.4,分量静电能E Elst是 -1.2,相除根本得不到13%啊,而且第二行里头,E SAPT才 -26.3,但是E Elst 却是 -40.4,这样除的话,分量值占总能量的150%还多了。。。。我也试了按绝对值算,可还是对不上。。。。

239

帖子

0

威望

2340

eV
积分
2579

Level 5 (御坂)

11#
发表于 Post on 2016-7-15 09:41:51 | 只看该作者 Only view this author
本帖最后由 brothers 于 2016-7-15 09:43 编辑
浅原问问 发表于 2016-7-15 09:11
老师我本来也想这么算的,但是您看表格中第一行数据,总相互作用E SAPT是 -1.4,分量静电能E Elst是 -1.2 ...

负值的都是吸引作用。第二行数据,40.4/(40.4+35.0+16.2+10.4)=40%。其余类此
十室之邑,必有忠信如丘焉,不如丘之好学也。

103

帖子

0

威望

629

eV
积分
732

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2016-7-15 11:01:59 | 只看该作者 Only view this author
brothers 发表于 2016-7-15 09:41
负值的都是吸引作用。第二行数据,40.4/(40.4+35.0+16.2+10.4)=40%。其余类此

哦。。。。谢谢您!除了第一行的数据其他的都能对上了,难道是数据小的时候有别的算法?或者是作者算错了?我再去找一些文献看看,太谢谢您啦!

239

帖子

0

威望

2340

eV
积分
2579

Level 5 (御坂)

13#
发表于 Post on 2016-7-15 13:42:30 | 只看该作者 Only view this author
浅原问问 发表于 2016-7-15 11:01
哦。。。。谢谢您!除了第一行的数据其他的都能对上了,难道是数据小的时候有别的算法?或者是作者算错了 ...
应该是他搞错了,正文IP1的比例那不也错了~
十室之邑,必有忠信如丘焉,不如丘之好学也。

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 13:14 , Processed in 0.160365 second(s), 23 queries , Gzip On.

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