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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】优化卟啉的Cu或Fe配合物总是收敛失败
121/2返回列表12下一页
查看: 2624  |  回复: 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的回帖

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

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

fineren

木虫 (正式写手)
引用回帖:
Originally posted by jjww0001 at 2011-02-18 13:33:47:

Gview显示的是否有键对计算没有影响的。recoli   说得很对，如果你还是用s=1的话，如果不是只有1个cu，那就是做成1价的铜了，Gview生成的输入文件最好能亲自检查一下再丢进去算。

谢谢，输入文件我也检查不出什么问题。请问，依你看，4个铜卟啉连在一起，自旋多重度应该选多少为好？按4个单电子算吗？

[ Last edited by fineren on 2011-2-18 at 19:11 ]
赞一下 回复此楼
12楼2011-02-18 17:00:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fineren

木虫 (正式写手)
引用回帖:
Originally posted by jjww0001 at 2011-02-18 11:18:47:
一个铜应该有一个未成对电子，你有4个铜，如果相互间离得比较远，显然不会有键合，应该有4个未成对电子，当然如果铜离得较近的话，也有可能未成对电子数变少。

Cu之间是没有成键的

[ Last edited by fineren on 2011-2-18 at 20:21 ]
赞一下 回复此楼
14楼2011-02-18 19:36:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fineren

木虫 (正式写手)
引用回帖:
Originally posted by caipingl at 2011-02-18 22:02:23:
2＋Cu，可能的自旋多重度都要考虑一下；此外，对于第一过渡系列也可以用全电子基组优化（不一定是过渡金属就要用赝势基组的）；当然泛函也可能会有影响的。

2价的Cu最外层是3d9，也是只有一个单电子吧？请问可能是自旋多重度都有哪些呢？全电子基组优化耗时太大，因为是4个铜卟啉。
赞一下 回复此楼
16楼2011-02-19 10:03:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fineren

木虫 (正式写手)
引用回帖:
Originally posted by caipingl at 2011-02-19 17:25:48:
高自旋：五重态，四个Cu单电子同为alpha或者beta，S=5；低自旋：单重态，四个Cu单电子两个为alpha两个为beta，S=1；还有可能是三重态，即三个为alpha一个为beta（或反之），S=3。

非常感谢
回复此楼
18楼2011-02-20 10:11:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fineren

木虫 (正式写手)
引用回帖:
Originally posted by levyleo at 2011-03-03 11:35:02:
Gaussian计算过渡金属收敛有时候会比较困难，ADF算过渡金属比较好

换软件？没买啊
回复此楼
21楼2011-03-05 12:10:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 proton00 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭