24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2801  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿厦大0856,306求调剂 +11 Bblinging 2026-04-11 11/550 2026-04-11 19:48 by lqspecial
[考研] 284求调剂 +10 archer.. 2026-04-09 11/550 2026-04-11 19:45 by laoshidan
[考研] 305求调剂 +6 玛卡巴卡boom 2026-04-11 6/300 2026-04-11 19:00 by wutongshun
[考研] 327求调剂 +4 拾光任染 2026-04-06 4/200 2026-04-11 12:01 by Delta2012
[考研] 求调剂 +6 archer.. 2026-04-09 8/400 2026-04-11 10:55 by zhq0425
[考研] 085402通信工程调剂,有4项学科竞赛国奖(电赛国二),硕士研究生调剂自荐信。 +5 m永o不v言o弃m 2026-04-09 5/250 2026-04-11 09:33 by zhq0425
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 284求调剂 +19 梵@@ 2026-04-06 21/1050 2026-04-10 21:12 by zhouxiaoyu
[考研] 22408 366分,本科211,一志愿西工大 +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 282,电气工程专业,求调剂,不挑专业 +9 jggshjkkm 2026-04-10 9/450 2026-04-10 14:55 by 逆水乘风
[考研] 一志愿华东师范生物学326分,求调剂 +8 刘墨墨 2026-04-09 8/400 2026-04-10 12:00 by pengliang8036
[考研] 调剂 +24 不逢春 2026-04-07 25/1250 2026-04-10 09:39 by 690616278
[考研] 085501机械英二77总分294求调剂,接受跨专业学习 +6 守法公民亓纪 2026-04-08 6/300 2026-04-09 15:55 by wp06
[考研] 1U盾记得记得就 +9 sanjin020722 2026-04-08 10/500 2026-04-09 14:11 by 诗与自由
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 265求调剂 +19 小木虫085600 2026-04-06 21/1050 2026-04-08 10:38 by 逆水乘风
[考研] 085602调剂 初试总分335 +3 19123253302 2026-04-06 3/150 2026-04-07 18:00 by jp9609
[考研] 324求调剂 +3 k可乐 2026-04-05 4/200 2026-04-06 09:54 by 蓝云思雨
信息提示
请填处理意见