计算化学公社

标题: 用Gaussian的counterpoise=2关键字计算两分子相互作用能，需要进行零点振动能校正吗？ [打印本页]

作者
Author: LSmall 时间: 2018-9-7 11:22
标题: 用Gaussian的counterpoise=2关键字计算两分子相互作用能，需要进行零点振动能校正吗？
各位前辈好！用Gaussian的counterpoise=2关键字计算A、B两分子的相互作用能，直接读取complexation energy = -6.56 kcal/mole (corrected)就是A、B分子在特定构型下的相互作用能吗？是否需要进行零点振动能校正吗？怎么计算这个零点振动能呢？目前使用的优化、频率验证方法是 m06-2x/6-311g(d,p)empiricaldispersion=gd3 opt freq ；能量计算方法是m06-2x/6-311+g(d,p) EmpiricalDispersion=GD3 counterpoise=2。感谢解答！

作者
Author: sobereva 时间: 2018-9-8 03:42
是
理论上研究结合能不用考虑热力学校正量（包括ZPE）

作者
Author: LSmall 时间: 2018-9-8 08:42
嗯嗯，明白了，谢谢Sob老师解惑！

作者
Author: granvia 时间: 2018-9-8 12:16

sobereva 发表于 2018-9-8 03:42
是
理论上研究结合能不用考虑热力学校正量（包括ZPE）

具体也看目的吧。比如氢键结合的二聚体，好像带零点能校正和不带零点校正的能量，实验上都可以测的

作者
Author: sobereva 时间: 2018-9-9 01:53

granvia 发表于 2018-9-8 12:16
具体也看目的吧。比如氢键结合的二聚体，好像带零点能校正和不带零点校正的能量，实验上都可以测的

当前没说和实验对比，跟实验对比显然得考虑热力学校正量
纯理论研究结合能不需要考虑ZPE或其它热力学校正量

作者
Author: TonyAugustine 时间: 2021-1-14 05:17

sobereva 发表于 2018-9-9 01:53
当前没说和实验对比，跟实验对比显然得考虑热力学校正量
纯理论研究结合能不需要考虑ZPE或其它热力学校 ...

可以请教一下sob老师，在这种情况下，如果想算零点能怎么算呢？
因为我在计算的时候也用了counterpoise，然而freq计算在使用counterpoise和不使用counterpoise是不一样的。而counterpoise的计算只给出了最终构型的零点能校正，没有monomer的零点能矫正。

作者
Author: wzkchem5 时间: 2021-1-14 10:37

TonyAugustine 发表于 2021-1-14 05:17
可以请教一下sob老师，在这种情况下，如果想算零点能怎么算呢？
因为我在计算的时候也用了counterpoise ...

建议opt、freq时不要做counterpoise，实在怕BSSE可以把基组加大一点，或者用orca的gCP方法校正。单点能计算用counterpoise就行了。

作者
Author: TonyAugustine 时间: 2021-1-14 23:39

wzkchem5 发表于 2021-1-14 10:37
建议opt、freq时不要做counterpoise，实在怕BSSE可以把基组加大一点，或者用orca的gCP方法校正。单点能计 ...

我可以理解为用counterpoise算的单点能，再去加上用freq去算的ZPC和thermal correction，去得到最后的ZPE和thermal energy吗

作者
Author: 喵星大佬 时间: 2021-1-15 01:04
你这两种结合能定义都不一样，参考态也不一样。考虑countpoise的情况，显然用的是结合之后片段的几何结构，这时候单算片段的话根本就不在极小点，频率和热力学校正量也就显得莫名其妙。类似于ALMO-EDA的那个文献里提的“绝热结合能”和“垂直结合能”的关系。如果你要算考虑热力学校正量的结合能，那片段的参考态结构也应该是各自的极小点结构，这时候也压根不存在counterpoise的问题。

作者
Author: sobereva 时间: 2021-1-15 07:54

TonyAugustine 发表于 2021-1-14 05:17
可以请教一下sob老师，在这种情况下，如果想算零点能怎么算呢？
因为我在计算的时候也用了counterpoise ...

类似的问题别在两个帖子里问
http://bbs.keinsci.com/thread-21294-1-1.html

作者
Author: zfy524 时间: 2022-8-7 13:25

喵星大佬发表于 2021-1-15 01:04
你这两种结合能定义都不一样，参考态也不一样。考虑countpoise的情况，显然用的是结合之后片段的几何结构， ...

老师你说的ALMO-EDA是哪篇文献啊

作者
Author: 喵星大佬 时间: 2022-8-7 15:37

zfy524 发表于 2022-8-7 13:25
老师你说的ALMO-EDA是哪篇文献啊

这篇

作者
Author: zfy524 时间: 2022-8-7 16:42

喵星大佬发表于 2022-8-7 15:37
这篇

谢谢老师

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