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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】优化卟啉的Cu或Fe配合物总是收敛失败
241/3返回列表123下一页
查看: 2552  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭