关于金属离子和卟啉结合能计算结果异常的求助

imasen · 发表于 Post on 2024-2-29 00:03:52

本帖最后由 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计算结果添加于附件。请问我在计算中忽略什么部分或者进行了什么不当操作吗？为什么会得到这么大的结合能？

AuI_TCPP.gjf (5.25 KB, 下载次数 Times of downloads: 11)

AuI_TCPP.log (416.32 KB, 下载次数 Times of downloads: 6)

sobereva · 发表于 Post on 2024-2-29 00:16:16

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

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

imasen · 发表于 Post on 2024-2-29 00:31:55

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

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

sobereva · 发表于 Post on 2024-2-29 00:40:16

你的描述里说AuI，但是结构里根本没有I

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

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

imasen · 发表于 Post on 2024-2-29 00:43:45

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

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 关于金属离子和卟啉结合能计算结果异常的求助