计算化学公社

标题: 静电能占总相互作用能的比例应该怎么算？ [打印本页]

作者
Author: 浅原问问 时间: 2016-7-13 19:06
标题: 静电能占总相互作用能的比例应该怎么算？
各位老师好，我想计算不同离子液体中静电能占总能量的比例，以此来探究静电能对其性质的影响。计算的时候用静电能除以总能量来算发现根本不对，我看到一个表格，不知道它里面的比例是怎么算出来的，请问老师，这个是有什么不一样公式或算法吗？我暂时在网上还没有查到。。。。谢谢老师！

作者
Author: sobereva 时间: 2016-7-13 19:53
做SAPT能量分解，此文第10节提了
Multiwfn支持的弱相互作用的分析方法概览
http://sobereva.com/252

不是占总能量，而是占总相互作用能。已经把你的标题改了。

作者
Author: 浅原问问 时间: 2016-7-14 10:32

sobereva 发表于 2016-7-13 19:53
做SAPT能量分解，此文第10节提了
Multiwfn支持的弱相互作用的分析方法概览
http://sobereva.com/252

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

作者
Author: brothers 时间: 2016-7-14 12:55

浅原问问发表于 2016-7-14 10:32
谢谢老师，我现在是可以计算出各项能量了，但是不会算静电能占总相互作用能的比例。比如那个表格中第一个 ...

原文Elst(%)是怎么定义的？

作者
Author: 浅原问问 时间: 2016-7-14 15:47

brothers 发表于 2016-7-14 12:55
原文Elst(%)是怎么定义的？

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

作者
Author: brothers 时间: 2016-7-14 18:42
本帖最后由 brothers 于 2016-7-14 20:02 编辑

浅原问问发表于 2016-7-14 15:47
应该就是静电能占总能量的比例。。。

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

太不靠谱了，第一行应该是把数据搞错了。其他各行都是静电能/所有吸引作用%

作者
Author: 数据挖掘 时间: 2016-7-14 18:53
静电能在所有吸引作用中的比例，就是这么算的。

作者
Author: 浅原问问 时间: 2016-7-14 20:53

brothers 发表于 2016-7-14 18:42
搞不清楚，一个简单的比例按说很好算。
没用过sapt，但是感觉直接拿静电能除以总能量好像没什么意义啊， ...

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

作者
Author: sobereva 时间: 2016-7-14 21:44

浅原问问发表于 2016-7-14 20:53
您好，我不是太明白您说的所有吸引作用是指什么？如果是总相互作用能的话，sapt中，把总相互作用能分成了 ...

你用相互作用能的某个分量值除以总相互作用能就完了

作者
Author: 浅原问问 时间: 2016-7-15 09:11

sobereva 发表于 2016-7-14 21:44
你用相互作用能的某个分量值除以总相互作用能就完了

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

作者
Author: brothers 时间: 2016-7-15 09:41
本帖最后由 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%。其余类此

作者
Author: 浅原问问 时间: 2016-7-15 11:01

brothers 发表于 2016-7-15 09:41
负值的都是吸引作用。第二行数据，40.4/(40.4+35.0+16.2+10.4)=40%。其余类此

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

作者
Author: brothers 时间: 2016-7-15 13:42

浅原问问发表于 2016-7-15 11:01
哦。。。。谢谢您！除了第一行的数据其他的都能对上了，难道是数据小的时候有别的算法？或者是作者算错了 ...

应该是他搞错了，正文IP1的比例那不也错了~

欢迎光临计算化学公社 (http://bbs.keinsci.com/)