计算化学公社

标题: 频率计算时SCF不收敛 [打印本页]

作者
Author: yukio07943 时间: 2020-12-8 19:13
标题: 频率计算时SCF不收敛
题主对一个结构优化完成后，想要计算抛除结构中F离子后的体系能量，在把F拉到一定距离后用freq关键词通入计算，结果却显示SCF不收敛输出文件中报错语句为
Convergence failure -- run terminated.
Error termination via Lnk1e in /opt/soft/gauss/g16c01_avx2/g16/l502.exe at Tue Dec 8 12:39:50 2020.
之后尝试了http://sobereva.com/61中的方法，均不奏效，以下是计算能量的输入输出文件，请各位老师多多指点

作者
Author: 北大-陶豫 时间: 2020-12-8 19:27
（1）你想算的话，直接删了F离子，然后再单独算F离子，然后加起来，应该会容易收敛得多。
（2）这篇文章应该会很有用：http://sobereva.com/82

作者
Author: yukio07943 时间: 2020-12-8 19:33

北大-陶豫发表于 2020-12-8 19:27
（1）你想算的话，直接删了F离子，然后再单独算F离子，然后加起来，应该会容易收敛得多。
（2）这篇文章应 ...

多谢指教，我明白了些

作者
Author: zjxitcc 时间: 2020-12-8 19:38

yukio07943 发表于 2020-12-8 19:33
多谢指教，我明白了些

单独算负离子的话，由于是F-负离子，要加弥散函数，不能仅用6-31G(d)。为各物质统一起见，可以给所有的F加上弥散函数。

作者
Author: yukio07943 时间: 2020-12-8 19:57

zjxitcc 发表于 2020-12-8 19:38
单独算负离子的话，由于是F-负离子，要加弥散函数，不能仅用6-31G(d)。为各物质统一起见，可以给所有的F ...

谢谢老师指点，还想问如果只单独求F离子能量时添加弥散，此前没有用弥散，对结果影响大吗

作者
Author: zjxitcc 时间: 2020-12-8 20:00

yukio07943 发表于 2020-12-8 19:57
谢谢老师指点，还想问如果只单独求F离子能量时添加弥散，此前没有用弥散，对结果影响大吗

不懂“此前没有用弥散”是什么意思。你没法统一啊，复合物没有弥散、两个单体其中一个没有弥散，另一个单体（即F-）有弥散，这样复合物-（单体1+单体2）这个式子都不成立了，不能做减法，因为基组不一样。

作者
Author: 北大-陶豫 时间: 2020-12-8 20:29
我已经给你算好了。感觉主要是靠 SCF=XQC 抑制了振荡……你之后的任务可以读我给你算好的波函数。

作者
Author: 北大-陶豫 时间: 2020-12-8 20:38

北大-陶豫发表于 2020-12-8 20:29
我已经给你算好了。感觉主要是靠 SCF=XQC 抑制了振荡……你之后的任务可以读我给你算好的波函数。

诶我用 GView 怎么看到 F 和 Na 电荷数都是 0 ……

作者
Author: 北大-陶豫 时间: 2020-12-8 21:02
我又试了一下：
（1）将钠离子和分子整个作为一个 Fragment，然后设置总电荷为+1，好像没什么用，该不收敛还是不收敛。
（2）将基组减小为 STO-3G，片段收敛的波函数作初猜仍然振荡。
我认为以后处理这种体系比较好的方法是：低基组下 SCF=XQC（用低基组降低计算量），而后逐步加大基组。

作者
Author: biogon 时间: 2020-12-8 21:23

北大-陶豫发表于 2020-12-8 21:02
我又试了一下：
（1）将钠离子和分子整个作为一个 Fragment，然后设置总电荷为+1，好像没什么用，该不收敛 ...

这个初始结构那么摆会导致scf非常难收敛，不知道这么摆的意图是什么，显然是不够合理的

作者
Author: sobereva 时间: 2020-12-8 23:20
初始结构摆得毫无化学意义，F-摆得那么远完全不知道想做什么，就算SCF能收敛，优化到平衡结构也需要非常多步数才能让F-逐渐靠过来

越离谱的结构SCF往往越不容易收敛，这就是我博文里明确说过的情况

搭建的结构异想天开

当前问题完全没有超越博文的范畴

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