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

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

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

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

浅原问问 发表于 2016-7-15 09:11
老师我本来也想这么算的，但是您看表格中第一行数据，总相互作用E SAPT是 -1.4，分量静电能E Elst是 -1.2 ...

负值的都是吸引作用。第二行数据，40.4/(40.4+35.0+16.2+10.4)=40%。其余类此

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

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

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

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

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

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

静电能在所有吸引作用中的比例，就是这么算的。

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

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

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

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

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

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

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

做SAPT能量分解，此文第10节提了
Multiwfn支持的弱相互作用的分析方法概览
http://sobereva.com/252

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