计算化学公社

标题: 询问一个关于过渡金属盐的计算方法和基组选择问题 [打印本页]

作者
Author: soulggl 时间: 2020-6-26 17:19
标题: 询问一个关于过渡金属盐的计算方法和基组选择问题
sobe老师，
我最近计算一些MnF2+H2O的体系，我用了cam-b3lyp（也做了d3矫正，但是几乎没有效果），基组都采用的是6-311++g(d,p) ，还算接近实验值，大概是理论值2.46 eV，实验值2.29 eV这样，相差0.1-0.2 eV，为了更好的计算结果，我尝试了b3lyp/LC-wpbe/PBE0//6-311++g(d,p) ，这几个结果差不多都偏离的更远，不如cam-b3lyp，另外CCSD(T)/aug-cc-pVTZ的结果也不好，阴离子结构T1值0.048，这个我不太会处理，就又尝试了CCSD(T,T1diag,full)，结果也不好。
您看这样的体系有没有什么更合适的计算方法推荐，谢谢您。

作者
Author: pwzhou 时间: 2020-6-26 19:10
你这个2.46eV具体是什么都没说清楚

作者
Author: sobereva 时间: 2020-6-27 00:19
完全看不明白你到底要算什么，测试那些方法目的是什么
如果是计算相互作用能，能量都是用kJ/mol或kcal/mol表示的，量化里没有用eV表示的
CAM-B3LYP算这种过渡金属配合物体系完全是莫名其妙，别瞎用。该用什么仔细看
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

b3lyp/LC-wpbe/PBE0//6-311++g(d,p) 这种表示更是严重错误

也不知道怎么冒出来一个“阴离子结构”。

不要瞎算。

作者
Author: soulggl 时间: 2020-6-28 17:52
sobe老师，我算的是VDE值，所以优化了阴离子结构，然后用阴离子结构的能量和中性结构能量做差减，用这个理论VDE值的结果和实验值做对比，eV的单位是我自己转换过来的。
我优化结构的时候分别算了这几种方法b3lyp/LC-wpbe/PBE0，基组都采用了6-311++g(d,p)
文献中我找到的算过渡金属盐类的有王来生老师算的CuBr2和CuCl2，用的是b3lyp，我试了一下结果并不好，所以又算了其他几种方法，但是感觉都不理想，所以想问您一下，对于MnF2+H2O这种体系有没有什么建议。

作者
Author: soulggl 时间: 2020-6-28 17:53

pwzhou 发表于 2020-6-26 19:10
你这个2.46eV具体是什么都没说清楚

垂直脱附能VDE值

作者
Author: sobereva 时间: 2020-6-29 01:23

soulggl 发表于 2020-6-28 17:52
sobe老师，我算的是VDE值，所以优化了阴离子结构，然后用阴离子结构的能量和中性结构能量做差减，用这个理 ...

6-311++G**太垃圾。Pople系列基组的弥散函数质量很低
尝试aug-cc-pVTZ。LC-wPBE这种泛函没考虑的价值
算得动就用CCSD(T)或DLPNO-CCSD(T)，或者双杂化泛函，诸如B2PLYP-D3(BJ)

作者
Author: soulggl 时间: 2020-6-29 18:57

sobereva 发表于 2020-6-29 01:23
6-311++G**太垃圾。Pople系列基组的弥散函数质量很低
尝试aug-cc-pVTZ。LC-wPBE这种泛函没考虑的价值
...

谢谢sob老师，我们用的是Gaussian，采用CCSD(T)有点太高了，我们可能算不动，我试试B2PLYP-D3(BJ)/aug-cc-pVTZ，看看结果怎样再跟您说。

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