频率计算时SCF不收敛

yukio07943 · 发表于 Post on 2020-12-8 19:13:03

题主对一个结构优化完成后，想要计算抛除结构中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中的方法，均不奏效，以下是计算能量的输入输出文件，请各位老师多多指点

北大-陶豫 · 发表于 Post on 2020-12-8 19:27:30

（1）你想算的话，直接删了F离子，然后再单独算F离子，然后加起来，应该会容易收敛得多。
（2）这篇文章应该会很有用：http://sobereva.com/82

yukio07943 · 发表于 Post on 2020-12-8 19:33:35

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

多谢指教，我明白了些

zjxitcc · 发表于 Post on 2020-12-8 19:38:36

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

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

yukio07943 · 发表于 Post on 2020-12-8 19:57:57

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

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

zjxitcc · 发表于 Post on 2020-12-8 20:00:15

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

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

北大-陶豫 · 发表于 Post on 2020-12-8 20:29:20

我已经给你算好了。感觉主要是靠 SCF=XQC 抑制了振荡……你之后的任务可以读我给你算好的波函数。

北大-陶豫 · 发表于 Post on 2020-12-8 20:38:41

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

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

北大-陶豫 · 发表于 Post on 2020-12-8 21:02:18

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

biogon · 发表于 Post on 2020-12-8 21:23:20

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

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

sobereva · 发表于 Post on 2020-12-8 23:20:45

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

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

搭建的结构异想天开

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 频率计算时SCF不收敛

浏览过的版块