24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2799  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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的回帖

levyleo

金虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
Gaussian计算过渡金属收敛有时候会比较困难,ADF算过渡金属比较好
我太渺小了。
19楼2011-03-03 11:35:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

linzhongaiguo

铜虫 (正式写手)

试一试 scf=xqc


heyo_123(金币+1): 鼓励交流! 2011-02-19 12:21:15
scf=xqc
DFT
2楼2011-02-17 15:00:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

fineren

木虫 (正式写手)

非常感谢二位
4楼2011-02-17 15:51:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fineren

木虫 (正式写手)

引用回帖:
Originally posted by jjww0001 at 2011-02-17 16:05:21:
这个地方不收敛一般是初始构型不好,可以先优化配体,金属配合物先用简单的配体优化,再把大配体替换上去。建议基组加上极化基(如3-21g*),你的构型总能量这么大,要有好几十个原子吧

初始构型是用以前算的Zn卟啉,换掉Zn原子得来的。有四个卟啉,体系很大
7楼2011-02-17 19:49:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fineren

木虫 (正式写手)

引用回帖:
Originally posted by recoli at 2011-02-17 16:52:22:
S**2   =   0.0000
自旋多重度设置对了吗?

不知道自旋多重度应该是多少,先用1计算的。4个卟啉的2价Cu配体,连在一起成一个封闭的环,您觉得多重度应该是多少呢?
8楼2011-02-17 19:50:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

recoli

金虫 (正式写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zhou2009(金币+3): 2011-02-18 09:37:51
引用回帖:
Originally posted by fineren at 2011-02-17 19:50:46:


不知道自旋多重度应该是多少,先用1计算的。4个卟啉的2价Cu配体,连在一起成一个封闭的环,您觉得多重度应该是多少呢?

二价锌的d电子是3d10,闭壳层,没有问题,但是旁边的铜不一样,二价铜是3d9所以是开壳层,自旋多重度应该是2。这是一个铜离子的情况。如果是多个铜离子的情况我也说不上来了,建议不要上来就算,先查查文献找点依据再计算。
9楼2011-02-17 22:15:59
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 proton00 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 283求调剂 086004考英二数二 +17 那个噜子 2026-04-10 18/900 2026-04-11 16:27 by 明月此时有
[考研] 266求调剂,一志愿哈工程电子信息,本科获多项国奖和省奖 +7 lumine1 2026-04-06 7/350 2026-04-11 15:31 by 学zh
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 化学工程调剂289 +43 yang婷 2026-04-07 49/2450 2026-04-11 10:37 by 紫曦紫棋
[考研] 电气工程专硕320求调剂 +4 小麻子111 2026-04-10 4/200 2026-04-11 10:19 by Delta2012
[考研] 346,工科求调剂 +3 moser233 2026-04-09 3/150 2026-04-11 10:04 by zhq0425
[考研] 调剂 +5 青灯不负 2026-04-09 5/250 2026-04-11 09:59 by zhq0425
[考研] 中药学调剂 初试324 +4 洋甘菊、 2026-04-10 6/300 2026-04-11 09:41 by gong120082
[考研] 314求调剂 +23 wakeluofu 2026-04-09 24/1200 2026-04-10 15:31 by MOF_Catal
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
[考研] 070300化学 求调剂 +13 73372112 2026-04-08 13/650 2026-04-09 20:22 by maddjdld
[考研] 生物学学硕,初试351分,求调剂 +4 …~、王…~ 2026-04-08 5/250 2026-04-08 21:49 by limeifeng
[考研] 一志愿南京航空航天大学 材料与化工329分求调剂 +11 Mr. Z 2026-04-05 12/600 2026-04-08 16:15 by luoyongfeng
[考研] 264求调剂 +11 麦小叮当 2026-04-07 11/550 2026-04-08 16:05 by 一只好果子?
[考研] 336求调剂,一志愿中科大 +9 墨彧 yuyu 2026-04-06 9/450 2026-04-08 11:24 by 想读书的菌菌
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 304求调剂 +4 luoye0105 2026-04-05 4/200 2026-04-06 21:05 by 木子君1218
[考研] 一志愿北交大材料工程总分358求调剂 +10 cs0106 2026-04-05 12/600 2026-04-06 19:41 by 无际的草原
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
信息提示
请填处理意见