请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 168|回复: 11

[新手求助] 求助稀土(Tb)作为中心及有机物的配合物SCF不收敛

[复制链接]

9

帖子

0

威望

23

eV
积分
32

Level 2 能力者

发表于 6 天前 | 显示全部楼层 |阅读模式
新手求助各位老师:我自己按照http://sobereva.com/61提供的方法尝试收敛,这些方法都不能收敛,我计算的是稀土(Tb)作为中心及有机物的配合物。计算的结果都是迭代128次后直接显示不能收敛,用QC方法显示L508错误。麻烦各位老师给看看。以下附件是我提供的输入输出文件。 A1.gjf (6.39 KB, 下载次数: 9)

257

帖子

2

威望

2085

eV
积分
2382

Level 5 (御坂)

BSJ Institute

发表于 6 天前 | 显示全部楼层
L508是no lower point found?
如果是这个的话我也想知道怎么解决 遇到这种情况都只好修改初始构型/计算水平了

9

帖子

0

威望

23

eV
积分
32

Level 2 能力者

 楼主| 发表于 6 天前 | 显示全部楼层
Accelerator 发表于 2019-10-8 11:13
L508是no lower point found?
如果是这个的话我也想知道怎么解决 遇到这种情况都只好修改初始构型/计算水 ...

L508是因为我使用了QC方法,不用这个方法的话就都是显示L501的错误,不管用哪种方法都不能收敛。您好,那对于初始构型您有什么好的建议吗?谢谢您

2万

帖子

25

威望

2万

eV
积分
46654

管理员

公社社长

发表于 6 天前 | 显示全部楼层
你都具体尝试过什么?诸如文中说的非常、极度、超级重要的尝试其它计算级别收敛的波函数当初猜波函数的做法你都尝试了哪些级别?
不可能显示L501错误,报错也是L502
我博文里明确强调没事别考虑用scf=qc
你这体系酰胺部分结构一看就扭曲得厉害,博文里明确说了结构扭曲容易导致SCF难收敛
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

9

帖子

0

威望

23

eV
积分
32

Level 2 能力者

 楼主| 发表于 6 天前 | 显示全部楼层
sobereva 发表于 2019-10-8 16:32
你都具体尝试过什么?诸如文中说的非常、极度、超级重要的尝试其它计算级别收敛的波函数当初猜波函数的做法 ...

老师,不好意思我打错了是显示L502错误。老师您说的酰胺部分的扭曲可以理解成就是画图的时候尽量大部分的键在一个平面上吗?我刚接触这个计算差不多半个月,您发的那些东西我也在学习中。以后还请您多多指教。谢谢老师的帮助。

473

帖子

0

威望

2434

eV
积分
2907

Level 5 (御坂)

发表于 6 天前 | 显示全部楼层
zhangjinkang 发表于 2019-10-8 18:00
老师,不好意思我打错了是显示L502错误。老师您说的酰胺部分的扭曲可以理解成就是画图的时候尽量大部分的 ...

多看几遍博文,会受益良多,总有解决办法。另外,应抛弃“Lxxx错误”这种无用低效说法,直接说有效报错信息。

2万

帖子

25

威望

2万

eV
积分
46654

管理员

公社社长

发表于 5 天前 | 显示全部楼层
zhangjinkang 发表于 2019-10-8 18:00
老师,不好意思我打错了是显示L502错误。老师您说的酰胺部分的扭曲可以理解成就是画图的时候尽量大部分的 ...

你先单独优化一下配体,以消除严重不合理的某些局部严重变形,再拿它当配体搭建配合物结构优化整个配合物
另外,你这个体系有很长的烃链,用B3LYP还不加D3是搞笑的,仔细看
谈谈“计算时是否需要加DFT-D3色散校正?”
http://sobereva.com/413http://bbs.keinsci.com/thread-9772-1-1.html
DFT-D色散校正的使用
http://sobereva.com/210

另外建议先把烃链用甲基替换,这样试解决SCF不收敛的做法的时候能节约不少时间
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

9

帖子

0

威望

23

eV
积分
32

Level 2 能力者

 楼主| 发表于 5 天前 | 显示全部楼层
sobereva 发表于 2019-10-9 08:38
你先单独优化一下配体,以消除严重不合理的某些局部严重变形,再拿它当配体搭建配合物结构优化整个配合物 ...

好的,谢谢老师,我先按您提供的办法解决,真是太感谢您了

8

帖子

1

威望

729

eV
积分
757

Level 4 (黑子)

发表于 5 天前 | 显示全部楼层
你这是正三价的Tb吧,正三价的Tb外面是有6个f单电子的,而你用的SDD基组默认是小核基组和赝势,并没有将f电子考虑进内核中,所以你的自旋多重度应该设置为7而不是1。Tb和Dy配合物就是比其他稀土配合物难优化,以我个人经验来说,在确定配体结构正确的前提下,先用大核基组优化,然后以得到的结构为初猜再用小核基组优化,期间还是可能遇到不收敛的问题,可以先降低收敛限,优化出一步再提升一步地慢慢“磨”,直到最后。
PS1:这篇文章不可不读(http://sobereva.com/373
PS2:配体上的烷基链可以先不加,等拿到主体的稳定结构以后再加上去,这样算的快。

9

帖子

0

威望

23

eV
积分
32

Level 2 能力者

 楼主| 发表于 5 天前 | 显示全部楼层
Harrypotter 发表于 2019-10-9 09:25
你这是正三价的Tb吧,正三价的Tb外面是有6个f单电子的,而你用的SDD基组默认是小核基组和赝势,并没有将f电 ...

谢谢大神,您真是说到点上了,谢谢谢谢

9

帖子

0

威望

23

eV
积分
32

Level 2 能力者

 楼主| 发表于 4 天前 | 显示全部楼层
Harrypotter 发表于 2019-10-9 09:25
你这是正三价的Tb吧,正三价的Tb外面是有6个f单电子的,而你用的SDD基组默认是小核基组和赝势,并没有将f电 ...

老师,您说用大核基组赝势,我查资料对于我这个体系来说是不是得用SDD的大核基组赝势,那这个基组应该怎么输入了,有点懵了。您能给详细说说吗?谢谢谢谢您。

8

帖子

1

威望

729

eV
积分
757

Level 4 (黑子)

发表于 4 天前 | 显示全部楼层
赝势基组的使用以及输入在sob老师这篇文章里说的很详细http://sobereva.com/60,搭配PS1中的文章(http://sobereva.com/373)将会无往不利。大核就用SDD的就行,从哪里找博文里也有说。
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-10-14 22:19 , Processed in 0.164156 second(s), 26 queries .

快速回复 返回顶部 返回列表