计算化学公社

标题: 能量分解疑问 [打印本页]

作者
Author: dy19930508 时间: 2016-8-14 21:50
标题: 能量分解疑问
各位老师好：
   我计算一个煤分子（约70个原子）吸附一个水分子时，它们之间以氢键作用吸附在一起。利用GAMESS做能量分解，先尝试分别利用BLYP/6-31+G(d,p)和BLYP-D3/6-31+G(d,p)进行计算，得到如下两张图中的结果。左图为只使用BLYP时得到的能量，右图为使用BLYP-D3得到的能量。
   为什么加入D3修正后electrostatic能量会变成正值并且这么大，polarization能量绝对值也会增大，而色散能几乎没有变化呢？
   我现在想到的是会不会是因为煤分子是柔性结构，其片段之间存在范德华作用，使得加入D3修正后引起能量的较大变化，这个变化值大于煤分子与水分子之间的相互作用，从而产生如右图所示的结果。
   还是由于方法和基组选择正确呢？
   谢谢各位老师！

作者
Author: sobereva 时间: 2016-8-15 07:00
DFT-D3校正不应直接用于LMOEDA，否则物理意义会错乱
应该在不用D3校正下做LMOEDA，然后自己算出DFT-D3对相互作用能的校正，作为色散作用能（也是GKS-EDA的划分方式）
而对于BLYP这样完全描述不了色散作用的泛函，当前输出的“色散能”绝对不是物理上的色散能，而应当作为相关能一项处理。

作者
Author: dy19930508 时间: 2016-8-15 08:57

sobereva 发表于 2016-8-15 07:00
DFT-D3校正不应直接用于LMOEDA，否则物理意义会错乱
应该在不用D3校正下做LMOEDA，然后自己算出DFT-D3对相 ...

原来是这样啊，谢谢sob老师！

作者
Author: dy19930508 时间: 2016-8-15 09:35

sobereva 发表于 2016-8-15 07:00
DFT-D3校正不应直接用于LMOEDA，否则物理意义会错乱
应该在不用D3校正下做LMOEDA，然后自己算出DFT-D3对相 ...

sob老师您好，如果我使用M06-2X/6-31+G(d,p)方法和基组得出的色散能，对于我这个70个原子的体系是不是近似准确的呢？

作者
Author: sobereva 时间: 2016-8-15 09:49

dy19930508 发表于 2016-8-15 09:35
sob老师您好，如果我使用M06-2X/6-31+G(d,p)方法和基组得出的色散能，对于我这个70个原子的体系是不是近 ...

用这个也得不到“色散能”，因为此时LMOEDA给出的色散项其实里面掺进去了很多并不是用来表现色散作用的电子相关作用能。

作者
Author: dy19930508 时间: 2016-8-15 10:45

sobereva 发表于 2016-8-15 09:49
用这个也得不到“色散能”，因为此时LMOEDA给出的色散项其实里面掺进去了很多并不是用来表现色散作用的 ...

sob老师，如果我用ORCA，BLYP-D3计算出色散校正能
这部分能量就是色散能，还是要加上利用GAMESS(BLYP方法)计算出的“色散能”才是真正的色散能呢？
这部分能量加上GAMESS(BLYP方法)计算出的相互作用能是不是更为准确真实的相互作用能？
十分感谢sob老师！

作者
Author: sobereva 时间: 2016-8-15 12:35

dy19930508 发表于 2016-8-15 10:45
sob老师，如果我用ORCA，BLYP-D3计算出色散校正能
这部分能量就是色散能，还是要加上利用GAMESS(BLYP方 ...

色散部分相互作用能你用grimme的DFT-D3程序来算就完了，这里写了
DFT-D色散校正的使用
http://sobereva.com/210

加到BLYP算的总相互作用能上就是BLYP-D3级别的相互作用能，比较真实。

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