捐赠本论坛 (Donate this site)
开启辅助访问 切换到宽版

计算化学公社

 找回密码 Forget password
 注册 Register
计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 对配合物进行优化时出现问题
12下一页 Next page
返回列表 Return 发新帖
Views: 17654|回复 Reply: 19
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 对配合物进行优化时出现问题

[复制链接 Copy URL]
2238663

33

帖子

0

威望

384

eV
积分
417

Level 3 能力者
跳转到指定楼层 Go to specific reply
楼主
发表于 Post on 2020-7-3 09:46:35 | 只看该作者 Only view this author |只看大图 View big picture 回帖奖励 |正序浏览 |阅读模式 Reading model
本帖最后由 2238663 于 2020-11-4 11:18 编辑

在对配合物进行构型优化时，采用CASSCF方法和ANO基组。出现了如下情况，不知如何解决？


回复 Reply

举报 Report

2238663

33

帖子

0

威望

384

eV
积分
417

Level 3 能力者
20#
 楼主 Author| 发表于 Post on 2020-7-4 09:33:05 | 只看该作者 Only view this author
sobereva 发表于 2020-7-4 04:22
你这个计算槽点太多。那文章作者9成可能是molcas的用户，而且还是CASSCF方面的资深研究者，你拿Gaussian生 ...

啊啊啊啊啊谢谢Sob老师
回复 Reply

举报 Report

sobereva

6万

帖子

99

威望

6万

eV
积分
125136

管理员

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

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

含稀土元素的体系是出了名的SCF难收敛，得非常耐心折腾。并且可以考虑用SCF更容易收敛的SDD大核版本的赝势+赝势基组。
北京科音自然科学研究中心（http://www.keinsci.com）致力于计算化学的发展和传播，长期开办极高质量的各种计算化学类培训：初级量子化学培训班、中级量子化学培训班、高级量子化学培训班、量子化学波函数分析与Multiwfn程序培训班、分子动力学与GROMACS培训班、CP2K第一性原理计算培训班，内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路！培训各种常见问题见《北京科音办的培训班FAQ》。
欢迎加入北京科音微信公众号获取北京科音培训的最新消息，并避免错过网上有价值的计算化学文章！
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群（群号见此链接），合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔，提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog：http://sobereva.com（发布大量原创计算化学相关博文）
Multiwfn主页：http://sobereva.com/multiwfn（十分强大、极为流行的量子化学波函数分析程序）
Google Scholar：https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate：https://www.researchgate.net/profile/Tian_Lu
回复 Reply

举报 Report

2238663

33

帖子

0

威望

384

eV
积分
417

Level 3 能力者
18#
 楼主 Author| 发表于 Post on 2020-7-3 20:35:52 | 只看该作者 Only view this author
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

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

举报 Report

2238663

33

帖子

0

威望

384

eV
积分
417

Level 3 能力者
17#
 楼主 Author| 发表于 Post on 2020-7-3 16:16:51 | 只看该作者 Only view this author
zjxitcc 发表于 2020-7-3 16:05
看起来是OK的。

非常感谢！！
回复 Reply

举报 Report

2238663

33

帖子

0

威望

384

eV
积分
417

Level 3 能力者
16#
 楼主 Author| 发表于 Post on 2020-7-3 16:16:25 | 只看该作者 Only view this author
喵星大佬 发表于 2020-7-3 16:11
我竟然愣了一下Dy是什么元素。。。。
看里面有没有金属-金属键，有的话用TPSS可能会好点
杂化泛函算得 ...

好滴！！非常感谢！！
回复 Reply

举报 Report

喵星大佬

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人
15#
发表于 Post on 2020-7-3 16:11:27 | 只看该作者 Only view this author
本帖最后由 喵星大佬 于 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
回复 Reply

举报 Report

zjxitcc

4289

帖子

4

威望

9540

eV
积分
13909

Level 6 (一方通行)

MOKIT开发者
14#
发表于 Post on 2020-7-3 16:05:43 | 只看该作者 Only view this author
2238663 发表于 2020-7-3 15:50
%chk=G:\five.chk
# opt=calcfc freq b3lyp/genecp

看起来是OK的。
自动做多参考态计算的程序MOKIT
回复 Reply

举报 Report

2238663

33

帖子

0

威望

384

eV
积分
417

Level 3 能力者
13#
 楼主 Author| 发表于 Post on 2020-7-3 15:50:51 | 只看该作者 Only view this author
本帖最后由 2238663 于 2020-11-4 11:18 编辑
zjxitcc 发表于 2020-7-3 15:46
从描述上看是恰当的。但是没见到输入文件，没法判断你是否写对了输入文件

写错了跑起来也是白搭

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






回复 Reply

举报 Report

2238663

33

帖子

0

威望

384

eV
积分
417

Level 3 能力者
12#
 楼主 Author| 发表于 Post on 2020-7-3 15:50:26 | 只看该作者 Only view this author
喵星大佬 发表于 2020-7-3 15:31
你这个体系做CASSCF也太大了吧。。。。641个基函数
不知道活性空间选了多少个轨道
感觉选一个合适的泛函 ...

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

举报 Report

zjxitcc

4289

帖子

4

威望

9540

eV
积分
13909

Level 6 (一方通行)

MOKIT开发者
11#
发表于 Post on 2020-7-3 15:46:36 | 只看该作者 Only view this author
2238663 发表于 2020-7-3 15:22
组里配的电脑没有到 到了可以搞搞 我把方法降到b3lyp 对稀土离子用SDD 对其他轻原子用了6-31g* 粗算一个 ...

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

写错了跑起来也是白搭
自动做多参考态计算的程序MOKIT
回复 Reply

举报 Report

喵星大佬

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人
10#
发表于 Post on 2020-7-3 15:31:32 | 只看该作者 Only view this author
你这个体系做CASSCF也太大了吧。。。。641个基函数
不知道活性空间选了多少个轨道
感觉选一个合适的泛函去跑还比较现实一点
回复 Reply

举报 Report

2238663

33

帖子

0

威望

384

eV
积分
417

Level 3 能力者
9#
 楼主 Author| 发表于 Post on 2020-7-3 15:22:30 | 只看该作者 Only view this author
zjxitcc 发表于 2020-7-3 15:11
Win10的linux子系统也不错，安装方便。不用装双系统啥的。

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

举报 Report

zjxitcc

4289

帖子

4

威望

9540

eV
积分
13909

Level 6 (一方通行)

MOKIT开发者
8#
发表于 Post on 2020-7-3 15:11:33 | 只看该作者 Only view this author
2238663 发表于 2020-7-3 14:37
是计算一个单核的稀土离子配合物 活性空间仅位于稀土金属离子周围 看文献上对稀土化合物计算采用CASSCF的 ...

Win10的linux子系统也不错，安装方便。不用装双系统啥的。
自动做多参考态计算的程序MOKIT
回复 Reply

举报 Report

2238663

33

帖子

0

威望

384

eV
积分
417

Level 3 能力者
7#
 楼主 Author| 发表于 Post on 2020-7-3 14:37:01 | 只看该作者 Only view this author
zjxitcc 发表于 2020-7-3 13:24
600+基函数的CASSCF，不知道你活性空间取多大。感觉这是个很艰难的任务，强烈建议你在足够了解CASSCF计算技 ...

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

举报 Report

12下一页 Next page
返回列表 Return 发新帖

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-20 12:18 , Processed in 0.199340 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list