应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】优化卟啉的Cu或Fe配合物总是收敛失败
241/3返回列表123下一页
查看: 2802  |  回复: 23
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

请教计算过这类配合物的大侠，这种情况怎么能够解决？调整初始结构已经试过多次。有没有通过加什么关键词能解决的？谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-02-17 14:44:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

jjww0001

金虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
heyo_123(金币+1): 鼓励交流！ 2011-02-19 12:21:31
这个地方不收敛一般是初始构型不好，可以先优化配体，金属配合物先用简单的配体优化，再把大配体替换上去。建议基组加上极化基(如3－21g*)，你的构型总能量这么大，要有好几十个原子吧
赞一下(3人) 回复此楼
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
自旋多重度设置对了吗？
赞一下(3人) 回复此楼
6楼2011-02-17 16:52:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

levyleo

金虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
Gaussian计算过渡金属收敛有时候会比较困难，ADF算过渡金属比较好
赞一下(2人) 回复此楼
我太渺小了。
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次循环、收敛还达不到要求时，计算会在这个优化的基础上自动进行下一轮优化，直至优化达到要求。
赞一下(2人) 回复此楼
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。这是一个铜离子的情况。如果是多个铜离子的情况我也说不上来了，建议不要上来就算，先查查文献找点依据再计算。
赞一下(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 的主题更新
241/3返回列表123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 267求调剂 +5 一个叶子zsw 2026-04-05 6/300 2026-04-11 19:46 by laoshidan
[考研] 268求调剂 +3 你好tg 2026-04-05 5/250 2026-04-11 16:45 by Delta2012
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 一志愿211生物学280分 求调剂 +7 李rien 2026-04-05 7/350 2026-04-11 11:16 by 逆水乘风
[考研] 283求调剂，工科！ +12 苏打水7777 2026-04-08 12/600 2026-04-11 10:28 by 逆水乘风
[考研] 080100力学316求调剂 +8 L_Hairui 2026-04-07 8/400 2026-04-11 10:00 by zhq0425
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 材料调剂 +5 hzhahg 2026-04-06 5/250 2026-04-10 10:10 by may_新宇
[考研] 一志愿双非085400电子信息344 求调剂，对材料和化学方向也感兴趣 +8 无情的小羊 2026-04-09 9/450 2026-04-10 09:30 by 松花缸1201
[考研] 材料299专硕求调剂 +10 +21 2026-04-09 10/500 2026-04-09 17:34 by 1753564080
[考研] 求调剂材料科学与工程一志愿985初试365分 +5 材化李可 2026-04-08 5/250 2026-04-09 17:00 by Lilly_Li
[考研] 085400电子信息类（川大控制工程）求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 349学科化学045106求调剂，化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 生物学308分求调剂（一志愿华东师大） +13 相信必会光芒万� 2026-04-06 16/800 2026-04-09 13:54 by 徐良白眉大侠
[考研] 085801 总分275 本科新能源 求调剂 +8 bradoner 2026-04-08 9/450 2026-04-09 13:43 by only周
[考研] 一志愿华南师范大学0702物理学305调剂 +4 念常安 2026-04-07 6/300 2026-04-08 22:53 by bljnqdcc
[考研] 机械工程313分找工科调剂 +3 双一流本科机械 2026-04-08 3/150 2026-04-08 20:41 by 土木硕士招生
[考研] 287求调剂 +6 Fnhc 2026-04-07 6/300 2026-04-08 10:05 by xingguangj
[考研] 302分求调剂 一志愿安徽大学085601 +12 zyx上岸！ 2026-04-04 12/600 2026-04-07 02:09 by BruceLiu320
[考研] 324求调剂 +3 k可乐 2026-04-05 4/200 2026-04-06 09:54 by 蓝云思雨
信息提示
关闭
请填处理意见
关闭