求助过渡金属与有机小分子配位的binding enthalpies的数值范围

Jenny_Zhao

请问各位老师，我在计算比较不同金属离子对同一有机小分子配位的Binding Enthalpies时，用TPSSh-D3(BJ)/6-31G**/SMD(CH2Cl2)进行opt和freq，revDSD-PBEP86-D3(BJ)/def2-TZVP/SMD(CH2Cl2)计算单点。Fe事先通过计算确定了是五重态，Cu、Zn根据常识分别设为二重态和单重态，也都检验了波函数的稳定性。但是我原本想用来比较的revDSD-PBEP86-D3(BJ)/def2-TZVP级别下的ΔH数值(绿色)Cu、Zn似乎有些离谱（几百和几十千卡的差别），而TPSSh-D3(BJ)/6-31G**级别下的ΔE、ΔH(黄色)都还在一个水平。请问这些绿色部分是可以接受的binding enthalpies数值范围，还是我选的方法基组不合适，又或是我忽略了哪些细节？请各位老师指点

hebrewsnabla

你说的检验了波函数稳定性，是指每个配合物单点的波函数稳定性都检验了么

wzkchem5

检查是不是理论级别不统一。差这么多的话，我感觉最大可能是因为写错关键词、读错文件等，导致理论级别不统一

Jenny_Zhao

hebrewsnabla 发表于 2024-1-30 20:24
你说的检验了波函数稳定性，是指每个配合物单点的波函数稳定性都检验了么

之前只做了结构优化级别下的stable=opt，昨天看到您的回复又检验了单点的波函数稳定性，其中Fe，Zn都是稳定的，但Cu的计算报错了，情况和http://bbs.keinsci.com/thread-32832-1-1.html这篇帖子一样，但我没看懂二楼老师的指导，暂时还没解决问题

hebrewsnabla

Jenny_Zhao 发表于 2024-1-31 11:41
之前只做了结构优化级别下的stable=opt，昨天看到您的回复又检验了单点的波函数稳定性，其中Fe，Zn都是稳 ...

高斯不支持直接做双杂化泛函的稳定性分析。间接的方法见 http://bbs.keinsci.com/thread-29232-1-1.html

Jenny_Zhao

wzkchem5 发表于 2024-1-30 21:55
检查是不是理论级别不统一。差这么多的话，我感觉最大可能是因为写错关键词、读错文件等，导致理论级别不统 ...

谢谢老师的提醒，我把输入文件都检查了，没有问题。最后检查发现是读取Zn，Cu数据的时候误读了双杂化泛函的"HF="后面的数值，实在是低级错误，抱歉耽误了老师们的时间。
这也是Sob老师在http://sobereva.com/488这篇博文3.4节强调了的问题，引以为戒。

Jenny_Zhao

hebrewsnabla 发表于 2024-1-31 11:59
高斯不支持直接做双杂化泛函的稳定性分析。间接的方法见 http://bbs.keinsci.com/thread-29232-1-1.html

谢谢老师。不过有个小疑问，为什么Fe、Zn体系的单点做stable=opt的时候没有因为高斯不支持双杂化泛函的稳定性分析而报错

hebrewsnabla

Jenny_Zhao 发表于 2024-1-31 12:52
谢谢老师。不过有个小疑问，为什么Fe、Zn体系的单点做stable=opt的时候没有因为高斯不支持双杂化泛函的稳 ...

要看你关键词怎么写了。