24小时热门版块排行榜    

查看: 2759  |  回复: 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的回帖

fineren

木虫 (正式写手)

引用回帖:
Originally posted by recoli at 2011-02-17 22:15:59:


二价锌的d电子是3d10,闭壳层,没有问题,但是旁边的铜不一样,二价铜是3d9所以是开壳层,自旋多重度应该是2。这是一个铜离子的情况。如果是多个铜离子的情况我也说不上来了,建议不要上来就算,先查查文献找 ...

非常感谢,这样多个铜卟啉连在一起的没发现文献,所以就先试了单重态。另外请教一下,我在Gview里画的Cu与卟啉是有配位键的,而且键长在gjf文件里好像都差不多:
  B1             2.04654752
   B2             2.04654752
   B3             2.04654751
   B4             2.04654752
但是我保存成gjf文件后再用Gview打开,Cu原子都不成键,这样对计算有没有影响呢?
10楼2011-02-18 08:32:34
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学调剂 +4 Surekei 2026-03-21 4/200 2026-03-22 07:18 by ilovexiaobin
[考研] 求调剂院校信息 +4 CX 330 2026-03-21 4/200 2026-03-21 23:48 by ms629
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 317求调剂 +9 申子申申 2026-03-19 15/750 2026-03-21 17:31 by 学员8dgXkO
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 材料工程(专)一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 初始318分求调剂(有工作经验) +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 307求调剂 +10 冷笙123 2026-03-17 10/500 2026-03-21 01:54 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 材料专硕英一数二306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 281求调剂(0805) +14 烟汐忆海 2026-03-16 25/1250 2026-03-20 15:47 by yuncha
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
信息提示
请填处理意见