计算化学公社

标题: 这个体系dmol3收敛了，为何gaussian还没收敛呢？ [打印本页]

作者
Author: zyj19831206 时间: 2018-9-13 17:35
标题: 这个体系dmol3收敛了，为何gaussian还没收敛呢？
一个弱相互作用的体系，用dmol3计算收敛了，但是用gaussian无法收敛，希望大家帮我看下哪里有问题。。

作者
Author: sobereva 时间: 2018-9-13 18:38
dmol3能收敛说明不了什么，收敛限，用的基组、泛函等都不一样
你这体系也不算弱相互作用体系，里面的分子的氮都和Fe有配位作用
另外，有的铁，比如Fe34，配位方式看起来不对。结构不合理的话，电子结构会很诡异，导致SCF也难以收敛

作者
Author: zyj19831206 时间: 2018-9-13 18:43

sobereva 发表于 2018-9-13 18:38
dmol3能收敛说明不了什么，收敛限，用的基组、泛函等都不一样
你这体系也不算弱相互作用体系，里面的分子 ...

那社长，这个结构是否用gaussian就无法收敛？

作者
Author: liyuanhe211 时间: 2018-9-13 21:54
本帖最后由 liyuanhe211 于 2018-9-13 21:55 编辑

zyj19831206 发表于 2018-9-13 18:43
那社长，这个结构是否用gaussian就无法收敛？

否。有的是解决“收敛”的关键词，不过你是什么不收敛都没搞清。。。另外你确定体系的基态是doublet吗？

作者
Author: zyj19831206 时间: 2018-9-13 22:23

liyuanhe211 发表于 2018-9-13 21:54
否。有的是解决“收敛”的关键词，不过你是什么不收敛都没搞清。。。另外你确定体系的基态是doublet吗？
...

感觉这个gjf再次打开的时候，Fe-O键都断开了，需要自己手动调节atom bond parameters才行，是不是就是因为这个造成了很多悬挂键才会导致不收敛？

作者
Author: liyuanhe211 时间: 2018-9-14 00:13
本帖最后由 liyuanhe211 于 2018-9-14 01:13 编辑

zyj19831206 发表于 2018-9-13 22:23
感觉这个gjf再次打开的时候，Fe-O键都断开了，需要自己手动调节atom bond parameters才行，是不是就是因 ...

量子化学的计算结果完全不受键级影响，问这个问题说明对量化的基本原理都有很大欠缺。铁配合物是相对不容易算准的物种，不像普通有机物一样从别人那儿抄个方法来，随便搞搞，计算机吐出来个数就是凑合能用的。不明白咋回事儿可能收敛完了都是个没有物理意义的数。而你这还是个大的、不“紧”的团簇，应该十分小心。

作者
Author: zyj19831206 时间: 2018-9-14 08:22

liyuanhe211 发表于 2018-9-14 00:13
量子化学的计算结果完全不受键级影响，问这个问题说明对量化的基本原理都有很大欠缺。铁配合物是相对不容 ...

抱歉，的确不是量化出身，所以基础一般，Fe原子是有点多。

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