计算化学公社

标题: CCSD(T)计算出现L502错误 [打印本页]

作者
Author: 氢原子 时间: 2019-5-4 14:22
标题: CCSD(T)计算出现L502错误
大家好，我计算一个含过渡金属的配合物时，先用M06L优化好了结构，再用CCSD(T)计算其三重态的单点能，一直出现l502错误。还提示自旋污染，请问怎么解决呢？

作者
Author: zjxitcc 时间: 2019-5-4 14:32
（1）如果不懂自旋污染啥意思，就不要管它；
（2）改正“Lxxx（数字）错误”这种不良习惯，而是直接说有效报错信息，更为直接、便于他人分析；
（3）这些问题老一套，早就有答案合辑了，提问前应该先搜索，去里面搜错误《Gaussian FAQ；新手求助报错时的注意事项》http://bbs.keinsci.com/thread-4829-1-1.html

作者
Author: 我本是个娃娃 时间: 2019-5-4 15:56
read it carefully

http://sobereva.com/61

作者
Author: 量化小菜鸡 时间: 2019-5-4 16:21
自旋污染你用UHF或多或少一定有，假如我理解的对的话。能量不一定被影响太多，但是波函数的形状会被影响。关键看你要干嘛。

作者
Author: sobereva 时间: 2019-5-4 16:40
SCF不收敛和自旋污染有没有以及多大毫无直接联系

作者
Author: 氢原子 时间: 2019-5-4 17:06

sobereva 发表于 2019-5-4 16:40
SCF不收敛和自旋污染有没有以及多大毫无直接联系

我已经尝试了减小基组、conver=6的方法，还是收敛失败。现在试scf=qc, 但是担心这种体系算CCSD(T)加了qc算不动啊。而且我只是算个能量而已，基于优化好的结构也这么难收敛吗？

作者
Author: sobereva 时间: 2019-5-4 17:14

氢原子发表于 2019-5-4 17:06
我已经尝试了减小基组、conver=6的方法，还是收敛失败。现在试scf=qc, 但是担心这种体系算CCSD(T)加了qc ...

当前给配体加弥散函数莫名其妙，本来CCSD(T)算这样的体系就费劲，还加这种对当前体系结果改进基本没意义的徒增耗时的东西，SCF难收敛也等于是自找的

没事别乱加弥散函数，这点下文写得非常清楚
谈谈弥散函数和“月份”基组
http://sobereva.com/119

Pople基组根本就不适合给CCSD(T)用

conver=6明显不能给CCSD(T)用

配体改成def-TZVP，若CCSD(T)还算不动就改用双杂化，或者ORCA跑DLPNO-CCSD(T)
甭优先想着scf=qc，又慢又容易收敛到不稳定波函数，此文明确说了
解决SCF不收敛问题的方法
http://sobereva.com/61

作者
Author: 氢原子 时间: 2019-5-4 18:12

sobereva 发表于 2019-5-4 17:14
当前给配体加弥散函数莫名其妙，本来CCSD(T)算这样的体系就费劲，还加这种对当前体系结果改进基本没意义 ...

老师，我配体换了def-TZVP还是收敛失败，是否可以出尝试SCF=vshift=300 这个体系都含有过渡金属

作者
Author: sobereva 时间: 2019-5-5 06:42

氢原子发表于 2019-5-4 18:12
老师，我配体换了def-TZVP还是收敛失败，是否可以出尝试SCF=vshift=300 这个体系都含有过渡金属

可。博文提到的都可以尝试

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