计算化学公社

标题: 关于金属离子和卟啉结合能计算结果异常的求助 [打印本页]

作者
Author: imasen 时间: 2024-2-29 00:03
标题: 关于金属离子和卟啉结合能计算结果异常的求助
本帖最后由 imasen 于 2024-2-29 00:28 编辑

我在计算TCPP卟啉和AuI的复合物时，采用ORCA5.0.3在TPSSh-D3/def2-SVP级别进行了几何优化。然后在wB97M-V/def2-TZVP级别下分别进行了TCPP(采用单独优化后的结构)、AuI-TCPP和AuI的单点计算，得到的结合能有-500多kcal/mol。考虑到没有加上弥散函数，我又改用Gaussian16在M06-2X-D3/ma-TZVP级别下做了结合能计算，并用counterpoise方法做了BSSE校正，得到的结果仍然有-600多kcal/mol，结果如下图，Gaussian计算结果添加于附件。请问我在计算中忽略什么部分或者进行了什么不当操作吗？为什么会得到这么大的结合能？
(, 下载次数 Times of downloads: 6)
(, 下载次数 Times of downloads: 4)
(, 下载次数 Times of downloads: 11)

(, 下载次数 Times of downloads: 6)

作者
Author: sobereva 时间: 2024-2-29 00:16
“wB97M-V/def2-TZVP级别下分别进行了单独优化” 意义不明，怎么又优化一遍
把Gaussian的counterpoise任务的输入输出文件上传
量子化学领域表示结合能要用kJ/mol或kcal/mol，不要用非主流的eV

当前算的是化学键作用能，带显著负电荷的区域又不在相互作用位点上，弥散跟当前问题没关系。

作者
Author: imasen 时间: 2024-2-29 00:31

sobereva 发表于 2024-2-29 00:16
“wB97M-V/def2-TZVP级别下分别进行了单独优化” 意义不明，怎么又优化一遍
把Gaussian的counterpoise任务 ...

sob老师您好，我已经修改问题描述，并添加附件。我指的是wB97M-V/def2-TZVP下计算单点的卟啉结构是单独优化后的而不是金属卟啉中的。当前结构下不做BSSE校正也仍有-452.02kcal/mol，这个结果合理么？

作者
Author: sobereva 时间: 2024-2-29 00:40
你的描述里说AuI，但是结构里根本没有I

这个数量级并不离谱，一个普通单键100 kcal/mol数量级的程度，当前牵扯4个键的相互作用，而且当前还是正负离子的作用

注意根据实际研究条件，确认卟啉衍生物部分用的净电荷对不对。另外考虑羧基的质子化状态选择是否得当。

作者
Author: imasen 时间: 2024-2-29 00:43

sobereva 发表于 2024-2-29 00:40
你的描述里说AuI，但是结构里根本没有I

这个数量级并不离谱，一个普通单键100 kcal/mol数量级的程度，当 ...

感谢老师的解答，应该是Au(I)，是我书写不当

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