计算化学公社

标题: ETS-NOCV和单点能的结果不匹配 [打印本页]

作者
Author: dfsdfv 时间: yesterday 23:42
标题: ETS-NOCV和单点能的结果不匹配
老师们好，我最近尝试用Orca计算了一些Pd(II)配合物(净电荷1，自旋多重度1)的单点能，并进行了ETS-NOCV分析，遇到一些问题比较费解，希望请教各位老师
1. 首先尝试了用下面的关键词计算了两个结构相近配合物以及结构片段的单点能：
! wB97M-V def2-TZVP SP SMD(Acetonitrile) RIJCOSX def2/J TightSCF

电荷和自旋多重度：
配体：0 1
Pd 1 1 (连有一个氢)

2. 计算完成后，将算得到的gbw文件转换为了Molden文件，计算出的结果如下：

	电子能量	总相互作用能
配合物A	-1592.680325485 hartree	-228.87 kcal/mol
配合物B	-1592.675424644 hartree	-230.07 kcal/mol
差值 (A-B)	-3.07 kcal/mol	1.2 kcal/mol

从单点能的结果来看，配合物A是比B要更稳定的，这也和实验结果吻合，但当尝试用ETS-NOCV计算出总的轨道作用反而是B更稳定，想请教各位老师如何理解这个结果，非常感谢！

作者
Author: sobereva 时间: 10 hour ago
ETS-NOCV只能用于真空下，不要带着溶剂模型产生用于Multiwfn做ETS-NOCV的波函数文件。先解决了这个硬伤再说
注意看下文弄清楚分析方法的关系
谈谈分子间结合能的构成以及分解分析思想

http://sobereva.com/733（http://bbs.keinsci.com/thread-50890-1-1.html）

如果你只需要得到溶剂下的相互作用能，而不需要做更进一步的分析，直接用溶剂模型下算的整体能量减片段能量就完了

作者
Author: dfsdfv 时间: halfhour ago

sobereva 发表于 2026-8-17 00:50
ETS-NOCV只能用于真空下，不要带着溶剂模型产生用于Multiwfn做ETS-NOCV的波函数文件。先解决了这个硬伤再说 ...

感谢老师回复，我尝试用在真空情况下做了计算，但是计算出来的总相互作用能依然是配合物B更高。从NOCV轨道来看，配体对金属中心的donation确实是A比B强1.8 kcal/mol，但是金属d轨道对配体的反馈则情况相反，B大约比A强了接近4 kcal/mol，请问是否是计算方法的原因，非常感谢！

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