计算化学公社

标题: 求助稀土(Tb)作为中心及有机物的配合物SCF不收敛 [打印本页]

作者
Author: zhangjinkang 时间: 2019-10-8 10:45
标题: 求助稀土(Tb)作为中心及有机物的配合物SCF不收敛
新手求助各位老师：我自己按照http://sobereva.com/61提供的方法尝试收敛，这些方法都不能收敛，我计算的是稀土(Tb)作为中心及有机物的配合物。计算的结果都是迭代128次后直接显示不能收敛，用QC方法显示L508错误。麻烦各位老师给看看。以下附件是我提供的输入输出文件。 (, 下载次数 Times of downloads: 33) (, 下载次数 Times of downloads: 8)

作者
Author: Accelerator 时间: 2019-10-8 11:13
L508是no lower point found?
如果是这个的话我也想知道怎么解决遇到这种情况都只好修改初始构型/计算水平了

作者
Author: zhangjinkang 时间: 2019-10-8 11:47

Accelerator 发表于 2019-10-8 11:13
L508是no lower point found?
如果是这个的话我也想知道怎么解决遇到这种情况都只好修改初始构型/计算水 ...

L508是因为我使用了QC方法，不用这个方法的话就都是显示L501的错误，不管用哪种方法都不能收敛。您好，那对于初始构型您有什么好的建议吗？谢谢您

作者
Author: sobereva 时间: 2019-10-8 16:32
你都具体尝试过什么？诸如文中说的非常、极度、超级重要的尝试其它计算级别收敛的波函数当初猜波函数的做法你都尝试了哪些级别？
不可能显示L501错误，报错也是L502
我博文里明确强调没事别考虑用scf=qc
你这体系酰胺部分结构一看就扭曲得厉害，博文里明确说了结构扭曲容易导致SCF难收敛

作者
Author: zhangjinkang 时间: 2019-10-8 18:00

sobereva 发表于 2019-10-8 16:32
你都具体尝试过什么？诸如文中说的非常、极度、超级重要的尝试其它计算级别收敛的波函数当初猜波函数的做法 ...

老师，不好意思我打错了是显示L502错误。老师您说的酰胺部分的扭曲可以理解成就是画图的时候尽量大部分的键在一个平面上吗？我刚接触这个计算差不多半个月，您发的那些东西我也在学习中。以后还请您多多指教。谢谢老师的帮助。

作者
Author: zjxitcc 时间: 2019-10-8 19:28

zhangjinkang 发表于 2019-10-8 18:00
老师，不好意思我打错了是显示L502错误。老师您说的酰胺部分的扭曲可以理解成就是画图的时候尽量大部分的 ...

多看几遍博文，会受益良多，总有解决办法。另外，应抛弃“Lxxx错误”这种无用低效说法，直接说有效报错信息。

作者
Author: sobereva 时间: 2019-10-9 08:38

zhangjinkang 发表于 2019-10-8 18:00
老师，不好意思我打错了是显示L502错误。老师您说的酰胺部分的扭曲可以理解成就是画图的时候尽量大部分的 ...

你先单独优化一下配体，以消除严重不合理的某些局部严重变形，再拿它当配体搭建配合物结构优化整个配合物
另外，你这个体系有很长的烃链，用B3LYP还不加D3是搞笑的，仔细看
谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413（http://bbs.keinsci.com/thread-9772-1-1.html）
DFT-D色散校正的使用
http://sobereva.com/210

另外建议先把烃链用甲基替换，这样试解决SCF不收敛的做法的时候能节约不少时间

作者
Author: zhangjinkang 时间: 2019-10-9 08:53

sobereva 发表于 2019-10-9 08:38
你先单独优化一下配体，以消除严重不合理的某些局部严重变形，再拿它当配体搭建配合物结构优化整个配合物 ...

好的，谢谢老师，我先按您提供的办法解决，真是太感谢您了

作者
Author: Harrypotter 时间: 2019-10-9 09:25
你这是正三价的Tb吧，正三价的Tb外面是有6个f单电子的，而你用的SDD基组默认是小核基组和赝势，并没有将f电子考虑进内核中，所以你的自旋多重度应该设置为7而不是1。Tb和Dy配合物就是比其他稀土配合物难优化，以我个人经验来说，在确定配体结构正确的前提下，先用大核基组优化，然后以得到的结构为初猜再用小核基组优化，期间还是可能遇到不收敛的问题，可以先降低收敛限，优化出一步再提升一步地慢慢“磨”，直到最后。
PS1：这篇文章不可不读（http://sobereva.com/373）
PS2：配体上的烷基链可以先不加，等拿到主体的稳定结构以后再加上去，这样算的快。

作者
Author: zhangjinkang 时间: 2019-10-9 18:35

Harrypotter 发表于 2019-10-9 09:25
你这是正三价的Tb吧，正三价的Tb外面是有6个f单电子的，而你用的SDD基组默认是小核基组和赝势，并没有将f电 ...

谢谢大神，您真是说到点上了，谢谢谢谢

作者
Author: zhangjinkang 时间: 2019-10-10 10:54

Harrypotter 发表于 2019-10-9 09:25
你这是正三价的Tb吧，正三价的Tb外面是有6个f单电子的，而你用的SDD基组默认是小核基组和赝势，并没有将f电 ...

老师，您说用大核基组赝势，我查资料对于我这个体系来说是不是得用SDD的大核基组赝势，那这个基组应该怎么输入了，有点懵了。您能给详细说说吗？谢谢谢谢您。

作者
Author: Harrypotter 时间: 2019-10-10 15:13
赝势基组的使用以及输入在sob老师这篇文章里说的很详细http://sobereva.com/60，搭配PS1中的文章（http://sobereva.com/373）将会无往不利。大核就用SDD的就行，从哪里找博文里也有说。

作者
Author: wang7344412 时间: 2019-10-17 15:29

Harrypotter 发表于 2019-10-9 09:25
你这是正三价的Tb吧，正三价的Tb外面是有6个f单电子的，而你用的SDD基组默认是小核基组和赝势，并没有将f电 ...

你好，我有个问题想请教，直接用大核基组可以么？为什么还要用小核继续优化，什么情况下需要继续用小核心继续优化啊？

作者
Author: Harrypotter 时间: 2019-10-17 21:41

wang7344412 发表于 2019-10-17 15:29
你好，我有个问题想请教，直接用大核基组可以么？为什么还要用小核继续优化，什么情况下需要继续用小核心 ...

这得看你的计算需求以及是否具有“强迫症”。大核优化的结构也是可以的，但是如果接下来的计算的性质涉及到f电子，那接下来必须把基组换成小核。此文http://sobereva.com/156的第2、3部分不得不看。

作者
Author: wang7344412 时间: 2019-10-18 09:07

Harrypotter 发表于 2019-10-17 21:41
这得看你的计算需求以及是否具有“强迫症”。大核优化的结构也是可以的，但是如果接下来的计算的性质涉及 ...

恩，明白了。谢谢。

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