计算化学公社

标题: 关于使用sobEDA方法对CO和N2进行能量分解的几点疑问 [打印本页]

作者
Author: duke_yzu 时间: 2024-10-30 22:16
标题: 关于使用sobEDA方法对CO和N2进行能量分解的几点疑问
关于使用sobEDA方法对 CO 和 N2 做能量分解（刚接触EDA，对这方面不是太懂）：

得到的结果虽然定性上是正确的：CO键更强。
但各项贡献与文献中有较大出入：文献中认为CO键更强原因在于较弱的泡利排斥，然而我计算的结果完全相反。根据我的计算结果来看反而是CO轨道间的相互作用更占优。

很可能是我的计算有问题，但不清楚问题出在哪里，希望有做过这方面工作的老师能帮忙看一下问题出在哪里。
（也考虑过片段的自旋多重度问题，但结果类似）
附件的压缩包提供了此能量分解相关的所有输入输出文件以及脚本文件，希望各位老师帮忙看一下，非常感谢！！

作者
Author: JCenter 时间: 2024-10-30 22:40
sobEDA是做片段间的能量分解解析的，你是想讨论N2与N2间的弱相互作用能、CO与CO间和弱相互作用能，还是N2与CO之间的呢，描述不清晰。此外，这也与你的团簇构型有关，不知道你的团簇结构从何而来，理应使用卢老师开发的molclus程序，寻找能量最低的团簇。《使用molclus程序做团簇构型搜索和分子构象搜索》（http://bbs.keinsci.com/thread-577-1-1.html）

作者
Author: duke_yzu 时间: 2024-10-31 08:49
你好，我讨论的是化学键，即CO中的C-O键，和N2中N-N键。不是讨论这两个分子间的弱相互作用。这两个分子已在B3LYP-D3(BJ)/6-31G*级别下进行了优化（如附图中所示）。

作者
Author: sobereva 时间: 2024-10-31 16:54
分解结果直接取决于片段状态的选取，看sobEDA原文里CH2片段用三重态和单重态时候区别
试图对照文献时片段状态的选取也必须严格一致，计算级别也应当尽可能统一

作者
Author: duke_yzu 时间: 2024-11-6 15:36
感谢各位老师提供的建议。参照文献中的片段自旋多重度选取，以及文献中提供的CO键长。重新用sobEDA方法做了能量分解，得到了与文献中大体一致的结果。 (, 下载次数 Times of downloads: 3)

作者
Author: duke_yzu 时间: 2024-11-6 15:46
采用完全相同的方法，只是片段自旋多重度不同，能量分解的结果却差别极大，对此不能很好理解。希望各位老师能指点一下 (, 下载次数 Times of downloads: 3)

作者
Author: zjxitcc 时间: 2024-11-6 16:13
本帖最后由 zjxitcc 于 2024-11-6 16:15 编辑

duke_yzu 发表于 2024-11-6 15:46
采用完全相同的方法，只是片段自旋多重度不同，能量分解的结果却差别极大，对此不能很好理解。希望各位老师 ...

C原子/O原子基态是三重态，N原子基态是四重态，划分片段时要注意尊重化学事实。

一种EDA方法对片段电荷和自旋采用不同的划分方式时，结果不同，这是全世界EDA方法面临的共同问题。EDA不是获得一个能量分解“标准答案”，事实上，能量分解没有标准答案，是不可观测量。EDA是当用户有明确的划分方式和相应的能量分解要求时，给用户提供一个便捷的工具。如果用户自己没想清楚、没定义清楚问题，那不应该开展EDA计算，而是应该先自己定义清楚问题（我这不是骂人，而是建议用户应该首先想明白他/她想要算什么）。

作者
Author: duke_yzu 时间: 2024-11-6 16:41
感谢老师指点！

作者
Author: gaojiahao 时间: 前天 15:56
各位老师好。sobEDA可以分析每个Π键的占比吗？类似与所给附件的分析。

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