计算化学公社

标题: 求助计算CASSCF收敛不了的问题 [打印本页]

作者
Author: Abandon-fmt 时间: 2019-6-18 08:51
标题: 求助计算CASSCF收敛不了的问题
按照高斯手册用CASSCF计算Conical Intersections，在计算激发态时，遇到收敛失败的报错，应该怎么看log的输出文件的信息来修改输入文件，可以使收敛成功？
输入：
%mem=20GB
%nprocshared=32
%chk=1.chk
%nosave
#p CAS(10,10,NRoot=2) 6-31G(d) Guess(Read) Pop(NaturalOrbital) Test Geom=Check

3

0 1

报错：

WARNING! : large rotation  I  J  =                      70
                     69
Step scaled by    0.2578534549246134
ITN= 64 MaxIt= 64 E= -998.4507551722 DE= 1.38D-02 Acc= 1.00D-08 Lan= 47
Convergence failure -- run terminated.

作者
Author: sobereva 时间: 2019-6-18 18:23
CASSCF的计算远没有你想的那么简单，SCF收敛难度也明显高于普通HF计算
仔细检查活性空间选取。同时算多个态一般都得用态平均。若迭代过程能量有收敛迹象增大迭代次数上限

作者
Author: 函数与激情 时间: 2019-6-18 20:05
高斯的CASSCF实在很不给力，大多分子默认的64圈数SCF收敛都吃力，你这距离收敛限差的还很远，可以减小活化空间同时增大圈数试试，G03的SCF收敛限比较松一点，也可以试试。

作者
Author: Abandon-fmt 时间: 2019-6-19 08:14

sobereva 发表于 2019-6-18 18:23
CASSCF的计算远没有你想的那么简单，SCF收敛难度也明显高于普通HF计算
仔细检查活性空间选取。同时算多个 ...

好的，谢谢sob老师！

作者
Author: Abandon-fmt 时间: 2019-6-19 08:14

函数与激情发表于 2019-6-18 20:05
高斯的CASSCF实在很不给力，大多分子默认的64圈数SCF收敛都吃力，你这距离收敛限差的还很远，可以减小活化 ...

好的好的，谢谢

作者
Author: Abandon-fmt 时间: 2019-7-4 14:55
本帖最后由 Abandon-fmt 于 2019-7-4 15:30 编辑

函数与激情发表于 2019-6-18 20:05
高斯的CASSCF实在很不给力，大多分子默认的64圈数SCF收敛都吃力，你这距离收敛限差的还很远，可以减小活化 ...

请问增大圈数的关键词这两个IOp(5/7=x)，scf(maxcycle=500)任选其一么？

作者
Author: Abandon-fmt 时间: 2019-7-4 15:40

sobereva 发表于 2019-6-18 18:23
CASSCF的计算远没有你想的那么简单，SCF收敛难度也明显高于普通HF计算
仔细检查活性空间选取。同时算多个 ...

请问sob老师，一般默认收敛的值是多少？可以把收敛值改大么？可以的话，关键词应该怎么写呢？

作者
Author: sobereva 时间: 2019-7-4 18:19

Abandon-fmt 发表于 2019-7-4 14:55
请问增大圈数的关键词这两个IOp(5/7=x)，scf(maxcycle=500)任选其 ...

用后者。有关键词就不要用IOp
设500没意义，超过一百多还不收敛再增大也白搭，继续算只会白浪费时间

作者
Author: sobereva 时间: 2019-7-4 18:19

Abandon-fmt 发表于 2019-7-4 15:40
请问sob老师，一般默认收敛的值是多少？可以把收敛值改大么？可以的话，关键词应该怎么写呢？

scf=conver=x

作者
Author: Abandon-fmt 时间: 2019-7-5 08:37

sobereva 发表于 2019-7-4 18:19
scf=conver=x

谢谢sob老师，我看到之前的帖子，算CASSCF的时候,G16中设成6是不是没有意义？后续精度需要8？

作者
Author: sobereva 时间: 2019-7-5 13:10

Abandon-fmt 发表于 2019-7-5 08:37
谢谢sob老师，我看到之前的帖子，算CASSCF的时候,G16中设成6是不是没有意义？后续精度需要8？

那是对单行列式SCF的情况而言的

作者
Author: Abandon-fmt 时间: 2019-7-5 13:58

sobereva 发表于 2019-7-5 13:10
那是对单行列式SCF的情况而言的

好的，谢谢老师！

作者
Author: Abandon-fmt 时间: 2019-7-5 14:01

sobereva 发表于 2019-7-5 13:10
那是对单行列式SCF的情况而言的

再请问sob老师，如果想用小基组收敛的波函数做为中/大基组的初猜，关键词应该加什么呢？

作者
Author: sobereva 时间: 2019-7-5 17:04

Abandon-fmt 发表于 2019-7-5 14:01
再请问sob老师，如果想用小基组收敛的波函数做为中/大基组的初猜，关键词应该加什么呢？

小基组下CASSCF算完了，chk里就有此时优化出来的轨道了。改成大基组，guess=read即可

作者
Author: Abandon-fmt 时间: 2019-7-5 19:06

sobereva 发表于 2019-7-5 17:04
小基组下CASSCF算完了，chk里就有此时优化出来的轨道了。改成大基组，guess=read即可

好的呢，谢谢老师，我去试试

作者
Author: Abandon-fmt 时间: 2019-7-8 13:16

sobereva 发表于 2019-7-5 17:04
小基组下CASSCF算完了，chk里就有此时优化出来的轨道了。改成大基组，guess=read即可

谢谢老师，用小基组的结果做初猜能算出来！

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