计算化学公社

标题: AuCl2和BBr3之间的弱相互作用能计算错误 [打印本页]

作者
Author: nzh19981019 时间: 2022-1-19 20:36
标题: AuCl2和BBr3之间的弱相互作用能计算错误
各位老师大家晚上好，我在计算AuCl2和BBr3之间的弱相互作用时，得到的作用能为正值，但是从他们的结合能为负值，结合距离也小于范德华半径之和，大家可以帮我看看是什么原因造成的吗，是我的自旋多重度不合适吗？我的输入文件和输出文件在下面了，谢谢。
(, 下载次数 Times of downloads: 20) (, 下载次数 Times of downloads: 10)

作者
Author: liyuanhe211 时间: 2022-1-20 03:42
这啥意思：“得到的作用能为正值，但是从他们的结合能为负值”，分别什么能

作者
Author: nzh19981019 时间: 2022-1-20 11:49

liyuanhe211 发表于 2022-1-20 03:42
这啥意思：“得到的作用能为正值，但是从他们的结合能为负值”，分别什么能

您好，作用能指的是二聚体的能量与二聚体中两个单体对应结构能量加和的差值，结合能是二聚体与两个没有形成二聚体前单体结构对应的能量加和的差值。因为两个单体形成二聚体时会变形来降低能量，作用能考虑了变形后降低的能量，结合能没有考虑单体的变形，所以作用能一般会比结合能更负，但是我的这个结果作用能确是正值，结合能为负值。作用能为正值意味着它们间没有弱相互作用，但是我其他的数据都显示存在弱作用，所以我感觉作用能的值不对，想问问您为什么会这样，是不是自旋多重度选择的不合理导致的。

作者
Author: 冰释之川 时间: 2022-1-20 11:50

liyuanhe211 发表于 2022-1-20 03:42
这啥意思：“得到的作用能为正值，但是从他们的结合能为负值”，分别什么能

我猜：作用能是指用电子能量算的；结合能是指优化两个片段后的能量？

作者
Author: nzh19981019 时间: 2022-1-20 13:29

冰释之川发表于 2022-1-20 11:50
我猜：作用能是指用电子能量算的；结合能是指优化两个片段后的能量？

作用能是通过高斯直接算出来的，结合能是二聚体的能量与两个片段没形成二聚体前单独结构优化能量加和的差值。不知道是否和您表述的一样。

作者
Author: sobereva 时间: 2022-1-20 23:10
别用那么诡异的自定义基组方法
仔细看
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60
算这个直接都用def2-QZVPP就完了，这么大基组下也不用考虑BSSE问题

这俩之间的作用已经是化学键范畴了，用counterpoise有害无益，看
计算化学键键能时考虑BSSE不仅是多余的甚至是有害的
http://sobereva.com/381（http://bbs.keinsci.com/thread-6047-1-1.html）

另外强烈不建议用MP2算这种体系，本来MP2算闭壳层体系都很烂，对于开壳层体系更是巨垃圾，仔细看
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

作者
Author: nzh19981019 时间: 2022-1-22 12:15

sobereva 发表于 2022-1-20 23:10
别用那么诡异的自定义基组方法
仔细看
详解Gaussian中混合基组、自定义基组和赝势基组的输入

好的，谢谢您。

作者
Author: 火星哥 时间: 2022-6-18 14:51
您好，想请教您，我再计算烷基自由基和Br阴离子的相互作用能，您的计算方法适合我吗？

作者
Author: nzh19981019 时间: 2022-11-10 09:58

火星哥发表于 2022-6-18 14:51
您好，想请教您，我再计算烷基自由基和Br阴离子的相互作用能，您的计算方法适合我吗？

如果自旋污染不大的话是可以的

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