计算化学公社

标题: 对配合物进行优化时出现问题 [打印本页]

作者
Author: 2238663 时间: 2020-7-3 09:46
标题: 对配合物进行优化时出现问题
本帖最后由 2238663 于 2020-11-4 11:18 编辑

在对配合物进行构型优化时，采用CASSCF方法和ANO基组。出现了如下情况，不知如何解决？
(, 下载次数 Times of downloads: 45)

作者
Author: zjxitcc 时间: 2020-7-3 10:19
本帖最后由 zjxitcc 于 2020-7-3 10:29 编辑

截图所示没有报错信息，是正常的输出内容。

PS：CASSCF计算别是简单直接写了关键词/方法，不能像DFT那样直接算，如果不熟悉的话先看看ORCA的CASSCF tutorial学习一下（可在ORCA官方论坛下载）

作者
Author: 2238663 时间: 2020-7-3 12:36

zjxitcc 发表于 2020-7-3 10:19
截图所示没有报错信息，是正常的输出内容。

PS：CASSCF计算别是简单直接写了关键词/方法，不能像DFT那样 ...

虽说正常但是我是用win系统跑到这里报了2070的错 link died

作者
Author: snljty 时间: 2020-7-3 13:01

2238663 发表于 2020-7-3 12:36
虽说正常但是我是用win系统跑到这里报了2070的错 link died

不要说“2070”，因为这个报错信息没有任何价值。Windows下所有Gaussian报错，全都显示2070。提问前可以先看看《Gaussian 常见报错及解决方法；新手求助报错时的注意事项》，http://bbs.keinsci.com/thread-4829-1-1.html

作者
Author: zjxitcc 时间: 2020-7-3 13:24
600+基函数的CASSCF，不知道你活性空间取多大。感觉这是个很艰难的任务，强烈建议你在足够了解CASSCF计算技巧的基础上再做计算。

假设你是32位的高斯，那内存有限制，根本没戏。就算是64位高斯，这在win下算效率太低，不知猴年马月，而且你的内存也未必跟得上。

更理想的情况是换到linux g16，那还得看你CASSCF初始轨道好不好，否则也是白谈。

作者
Author: 2238663 时间: 2020-7-3 14:28

snljty 发表于 2020-7-3 13:01
不要说“2070”，因为这个报错信息没有任何价值。Windows下所有Gaussian报错，全都显示2070。提问前可以 ...

收到！谢谢你

作者
Author: 2238663 时间: 2020-7-3 14:37

zjxitcc 发表于 2020-7-3 13:24
600+基函数的CASSCF，不知道你活性空间取多大。感觉这是个很艰难的任务，强烈建议你在足够了解CASSCF计算技 ...

是计算一个单核的稀土离子配合物活性空间仅位于稀土金属离子周围看文献上对稀土化合物计算采用CASSCF的的基组进行计算所以我也就用这个方法了在出现这个错误以后我试过用更为低级的方法比如DFT的b3lyp和PBEPBE 都出现了相同的现象
所以我也估计就是得在linux系统上才可以新手一枚非常感谢你的帮助

作者
Author: zjxitcc 时间: 2020-7-3 15:11

2238663 发表于 2020-7-3 14:37
是计算一个单核的稀土离子配合物活性空间仅位于稀土金属离子周围看文献上对稀土化合物计算采用CASSCF的 ...

Win10的linux子系统也不错，安装方便。不用装双系统啥的。

作者
Author: 2238663 时间: 2020-7-3 15:22

zjxitcc 发表于 2020-7-3 15:11
Win10的linux子系统也不错，安装方便。不用装双系统啥的。

组里配的电脑没有到到了可以搞搞我把方法降到b3lyp 对稀土离子用SDD 对其他轻原子用了6-31g* 粗算一个结果结果它跑起来了请问这样可行吗

作者
Author: 喵星大佬 时间: 2020-7-3 15:31
你这个体系做CASSCF也太大了吧。。。。641个基函数
不知道活性空间选了多少个轨道
感觉选一个合适的泛函去跑还比较现实一点

作者
Author: zjxitcc 时间: 2020-7-3 15:46

2238663 发表于 2020-7-3 15:22
组里配的电脑没有到到了可以搞搞我把方法降到b3lyp 对稀土离子用SDD 对其他轻原子用了6-31g* 粗算一个 ...

从描述上看是恰当的。但是没见到输入文件，没法判断你是否写对了输入文件

写错了跑起来也是白搭

作者
Author: 2238663 时间: 2020-7-3 15:50

喵星大佬发表于 2020-7-3 15:31
你这个体系做CASSCF也太大了吧。。。。641个基函数
不知道活性空间选了多少个轨道
感觉选一个合适的泛函 ...

嗯你说得对但是我不知道怎么选择合适的方法去优化和频率分析因为文献里基本上都是只介绍了直接计算性质用CASSCF……

作者
Author: 2238663 时间: 2020-7-3 15:50
本帖最后由 2238663 于 2020-11-4 11:18 编辑

zjxitcc 发表于 2020-7-3 15:46
从描述上看是恰当的。但是没见到输入文件，没法判断你是否写对了输入文件

写错了跑起来也是白搭

%chk=G:\five.chk
# opt=calcfc freq b3lyp/genecp

作者
Author: zjxitcc 时间: 2020-7-3 16:05

2238663 发表于 2020-7-3 15:50
%chk=G:\five.chk
# opt=calcfc freq b3lyp/genecp

看起来是OK的。

作者
Author: 喵星大佬 时间: 2020-7-3 16:11
本帖最后由喵星大佬于 2020-7-3 16:24 编辑

2238663 发表于 2020-7-3 15:50
%chk=G:\five.chk
# opt=calcfc freq b3lyp/genecp

我竟然愣了一下Dy是什么元素。。。。
看里面有没有金属-金属键，有的话用TPSS可能会好点，多参考性质强的体系SCAN、M06L、MN15L比较好，泛函选择可参考http://sobereva.com/272
杂化泛函算得金属-金属键长会偏差比较大。不过如果用纯泛函的话，ORCA倒是可以算得飞起，结合Gaussian的优化
具体见sob老师的帖子http://bbs.keinsci.com/thread-10141-1-1.html

作者
Author: 2238663 时间: 2020-7-3 16:16

喵星大佬发表于 2020-7-3 16:11
我竟然愣了一下Dy是什么元素。。。。
看里面有没有金属-金属键，有的话用TPSS可能会好点
杂化泛函算得 ...

好滴！！非常感谢！！

作者
Author: 2238663 时间: 2020-7-3 16:16

zjxitcc 发表于 2020-7-3 16:05
看起来是OK的。

非常感谢！！

作者
Author: 2238663 时间: 2020-7-3 20:35

zjxitcc 发表于 2020-7-3 16:05
看起来是OK的。

No special actions if energy rises.
Defaulting to unpruned grid for atomic number  66.
Defaulting to unpruned grid for atomic number  66.
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did    9 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did    5 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did    7 forward-backward iterations
EnCoef did 26 forward-backward iterations
EnCoef did 25 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did 50 forward-backward iterations
EnCoef did 25 forward-backward iterations
EnCoef did 38 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
Rare condition: small coef for last iteration:  0.858D-15
Restarting incremental Fock formation.
EnCoef did 52 forward-backward iterations
Rare condition: small coef for last iteration:  0.000D+00
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did 19 forward-backward iterations
Rare condition: small coef for last iteration:  0.333D-15
EnCoef did 98 forward-backward iterations
Rare condition: small coef for last iteration:  0.000D+00
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did 12 forward-backward iterations
Rare condition: small coef for last iteration:  0.111D-15
EnCoef did 96 forward-backward iterations
Rare condition: small coef for last iteration:  0.000D+00
EnCoef did 100 forward-backward iterations
EnCoef did 32 forward-backward iterations
Rare condition: small coef for last iteration: -0.888D-15
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did 24 forward-backward iterations
Rare condition: small coef for last iteration: -0.111D-14
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
Restarting incremental Fock formation.
EnCoef did 24 forward-backward iterations
Rare condition: small coef for last iteration: -0.222D-15
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did 24 forward-backward iterations
Rare condition: small coef for last iteration:  0.133D-14
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did 24 forward-backward iterations
Rare condition: small coef for last iteration:  0.333D-15
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did 24 forward-backward iterations
Rare condition: small coef for last iteration:  0.999D-15
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did 24 forward-backward iterations
Rare condition: small coef for last iteration: -0.222D-15
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
Restarting incremental Fock formation.
EnCoef did 24 forward-backward iterations
Rare condition: small coef for last iteration:  0.666D-15
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did 24 forward-backward iterations
Rare condition: small coef for last iteration:  0.222D-15
EnCoef did 100 forward-backward iterations
EnCoef did 100 forward-backward iterations
EnCoef did 24 forward-backward iterations
Rare condition: small coef for last iteration:  0.133D-14

出现这种情况知道不是错误但是是不是意味着无法收敛啊

作者
Author: sobereva 时间: 2020-7-4 04:22
你这个计算槽点太多。那文章作者9成可能是molcas的用户，而且还是CASSCF方面的资深研究者，你拿Gaussian生搬硬套显然不行
ANO-RCC-VDZP这种完全广义收缩基组在Gaussian里慢得要死
Gaussian的CASSCF做得不怎么样，收敛性差、速度慢，你拿它优化本来就不好搞的稀土配合物，几乎没什么希望
当前问题直接用合适的泛函做opt freq就完事了，你用CASSCF完全是瞎折腾

18L的信息明显是出现了SCF难收敛的情况。记得加上#P监控收敛趋势，并且仔细看
解决SCF不收敛问题的方法
http://sobereva.com/61

含稀土元素的体系是出了名的SCF难收敛，得非常耐心折腾。并且可以考虑用SCF更容易收敛的SDD大核版本的赝势+赝势基组。

作者
Author: 2238663 时间: 2020-7-4 09:33

sobereva 发表于 2020-7-4 04:22
你这个计算槽点太多。那文章作者9成可能是molcas的用户，而且还是CASSCF方面的资深研究者，你拿Gaussian生 ...

啊啊啊啊啊谢谢Sob老师

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