24小时热门版块排行榜    

查看: 2761  |  回复: 23
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

fineren

木虫 (正式写手)

[交流] 【求助】优化卟啉的Cu或Fe配合物总是收敛失败 已有10人参与

用混合基组,对金属原子用Lanl2dz,其他原子用3-21G,出错显示如下:
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RB+HF-LYP) =  -10031.4527509     A.U. after  129 cycles
             Convg  =    0.6966D-04             -V/T =  2.0366
             S**2   =   0.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in /local/apl/atx/g03e01/l502.exe

请教计算过这类配合物的大侠,这种情况怎么能够解决?调整初始结构已经试过多次。有没有通过加什么关键词能解决的?谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jjww0001

金虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by fineren at 2011-02-18 08:32:34:


非常感谢,这样多个铜卟啉连在一起的没发现文献,所以就先试了单重态。另外请教一下,我在Gview里画的Cu与卟啉是有配位键的,而且键长在gjf文件里好像都差不多:
  B1             2.04654752
   B2        ...

Gview显示的是否有键对计算没有影响的。recoli   说得很对,如果你还是用s=1的话,如果不是只有1个cu,那就是做成1价的铜了,Gview生成的输入文件最好能亲自检查一下再丢进去算。
11楼2011-02-18 13:33:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答

zhou2009

版主 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
heyo_123(金币+1): 鼓励交流! 2011-02-19 12:21:08
出现Convergence failure -- run terminated,不收敛的原因常仅仅是因为最大循环数不够(默认64),不足以达到预定的收敛指标(某次循环前和循环后的电荷密度差别小于1.00D-06)。

可以通过设置maxcyc=N来增大最大循环步数。

或者采用关键词iop(5/13=1),此时计算到一个默认64次循环、收敛还达不到要求时,计算会在这个优化的基础上自动进行下一轮优化,直至优化达到要求。
3楼2011-02-17 15:47:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jjww0001

金虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
heyo_123(金币+1): 鼓励交流! 2011-02-19 12:21:31
这个地方不收敛一般是初始构型不好,可以先优化配体,金属配合物先用简单的配体优化,再把大配体替换上去。建议基组加上极化基(如3-21g*),你的构型总能量这么大,要有好几十个原子吧
5楼2011-02-17 16:05:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

recoli

金虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zhou2009(金币+1): 2011-02-18 09:38:40
S**2   =   0.0000
自旋多重度设置对了吗?
6楼2011-02-17 16:52:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 293求调剂 +12 zjl的号 2026-03-16 17/850 2026-03-22 16:51 by i_cooler
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-20 3/150 2026-03-22 16:00 by ColorlessPI
[考研] 招08考数学 +3 laoshidan 2026-03-20 7/350 2026-03-22 12:54 by 在风落中
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
信息提示
请填处理意见