计算化学公社

标题: 过渡金属配合物SCF不收敛问题求助 [打印本页]

作者
Author: sun666 时间: 2021-12-3 14:03
标题: 过渡金属配合物SCF不收敛问题求助
本帖最后由 sun666 于 2021-12-3 14:10 编辑

在对过渡金属配合物优化计算单点能时看到文献对配体使用的6-311g(d) 对过渡金属用的Stuttgart小核赝势。因此采用了同样的方法基组，但是优化过渡金属配合物出现了SCF不收敛问题。输入输出文件在附件
有几个问题请教
1.我对Co原子直接用的默认SDD是小核赝势吗，跟文文献中的“Stuttgart小核赝势”是否一样（如何区别用的为大核还是小核）
2.自选多重都我设置为4不知是否合理

请各位老师指点如何解决不收敛的问题

作者
Author: wzkchem5 时间: 2021-12-3 15:43
按照http://sobereva.com/61的方法解决，如果每个方法都尝试多次，没有遗漏哪怕一种方法，还是不收敛，再来问，问的时候详细描述这篇文章里的哪些方法让收敛变好了，哪些让收敛变差了，这样我们好知道应该加什么关键词

作者
Author: sun666 时间: 2021-12-4 13:35

wzkchem5 发表于 2021-12-3 15:43
按照http://sobereva.com/61的方法解决，如果每个方法都尝试多次，没有遗漏哪怕一种方法，还是不收敛，再来 ...

谢谢老师。
尝试了以下办法：
（1）SCF=vshift=400
（2）SCF=conver=6、
(3) 对轻原子的基组换为6-31g(d),Co依然用SDD
这三种依然没有任何效果，SCF仍不收敛
（4）gencp换成pbe1pbe/def2svp方法基组：结构优化15步后依然SCF不收敛，但是最后一步结构优化Maximum Force 和RMS Force已经为YES

作者
Author: wzkchem5 时间: 2021-12-4 16:53

sun666 发表于 2021-12-4 06:35
谢谢老师。
尝试了以下办法：
（1）SCF=vshift=400

继续试，那个帖子里还有很多种方法。
另外，我觉得你要算的应该是溶液相反应吧，如果是的话加上溶剂模型再试

作者
Author: sun666 时间: 2021-12-10 21:52
本帖最后由 sun666 于 2021-12-10 21:56 编辑

作者
Author: mfdsrax2 时间: 2021-12-10 23:00

sun666 发表于 2021-12-10 21:52
**** 本内容被作者隐藏 ****

500分，想看一下帖子都不行吗

作者
Author: wzkchem5 时间: 2021-12-10 23:24

sun666 发表于 2021-12-10 14:52
**** 本内容被作者隐藏 ****

水分子的方向摆放不合理。建议查一下过渡金属阳离子和水分子的配合物里，水分子的取向一般是什么样的，有一个大致概念以后再做计算

作者
Author: sun666 时间: 2021-12-11 00:38

wzkchem5 发表于 2021-12-10 23:24
水分子的方向摆放不合理。建议查一下过渡金属阳离子和水分子的配合物里，水分子的取向一般是什么样的，有 ...

这个结构是真空中优化的结果，做了震动分析也没有虚频。您的意思是这个结构不适合做为庚烷溶剂环境下的初始结构吗。水分子是否会与P=O形成氢键呢

作者
Author: 喵星大佬 时间: 2021-12-11 00:44
你到底是SCF不收敛还是结构优化不收敛。。。。还有这玩意为啥自旋多重度是4，这不是闭壳层的嘛

作者
Author: wzkchem5 时间: 2021-12-11 03:16

sun666 发表于 2021-12-10 17:38
这个结构是真空中优化的结果，做了震动分析也没有虚频。您的意思是这个结构不适合做为庚烷溶剂环境下的初 ...

我刚没意识到是非极性溶剂里，那还有一定可能。不过还是按你楼下说的，先确定电荷、自旋多重度没错再看

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