计算化学公社

 找回密码 Forget password
 注册 Register

Cu配合物优化不收敛问题求助

查看数: 8865 | 评论数: 6 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2018-6-3 13:43

正文摘要:

前辈老师们好,我想询问一个Cu配合物优化时选用基组的困扰。 背景:我按照chemcomm上一篇2010年的实验文章的化合物进行优化,结构如图1所示。文献DOI: 10.1039/c0cc01864c。该文献也有一点计算,但是对XRD得到的结 ...

回复 Reply

luktian 发表于 Post on 2018-6-4 09:41:38
liyuanhe211 发表于 2018-6-4 08:59
应当分清是优化不收敛还是SCF不收敛。你的L1.log是SCF不收敛,不是优化不收敛。
加上弥散函数之后SCF容易 ...

感谢前辈您的指点,我之前写帖子的时候给混淆了SCF不收敛和优化不收敛;
有关弥散函数的部分我去尝试一下,感谢您的指路;
关于电荷,前面一位前辈也帮我指出了,现在改正后正在计算,感谢您的建议
liyuanhe211 发表于 Post on 2018-6-4 08:59:26
应当分清是优化不收敛还是SCF不收敛。你的L1.log是SCF不收敛,不是优化不收敛。
加上弥散函数之后SCF容易不收敛是正常合理的,看sobereva.com/61(尤其是其中与弥散有关的选项)解决。
另外你做的是Cu(0),文献中是Cu(I),确认是否是自己想要的。
  1. Coeff:      0.103D+00-0.418D-01 0.103D-01 0.251D-02-0.126D+00 0.297D+00
  2. Coeff:      0.896D+00-0.780D+00-0.390D+01 0.369D+01 0.444D+00-0.157D+00
  3. Coeff:     -0.633D+00 0.168D+00 0.499D+01-0.543D+01 0.958D+00-0.293D+01
  4. Coeff:      0.177D+01 0.167D+01
  5. Gap=     0.045 Goal=   None    Shift=    0.000
  6. Gap=     0.176 Goal=   None    Shift=    0.000
  7. RMSDP=1.00D-03 MaxDP=1.46D-01 DE= 7.72D-04 OVMax= 8.49D-02

  8. >>>>>>>>>> Convergence criterion not met.
  9. SCF Done:  E(UwB97XD) =  -2640.91781552     A.U. after  129 cycles
  10.             NFock=128  Conv=0.10D-02     -V/T= 2.2075
  11. <Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7659 S= 0.5079
  12. <L.S>= 0.000000000000E+00
  13. KE= 2.187074654595D+03 PE=-1.900821066278D+04 EE= 7.627967488886D+03
  14. Annihilation of the first spin contaminant:
  15. S**2 before annihilation     0.7659,   after     0.7502
  16. Convergence failure -- run terminated.
  17. Error termination via Lnk1e in /home/software/g16/l502.exe at Tue May 29 02:31:32 2018.
  18. Job cpu time:      54 days  4 hours 16 minutes 40.5 seconds.
  19. Elapsed time:       1 days 16 hours 54 minutes 18.9 seconds.
  20. File lengths (MBytes):  RWF=   1766 Int=      0 D2E=      0 Chk=    146 Scr=      1
复制代码


luktian 发表于 Post on 2018-6-3 22:07:33
h840473807 发表于 2018-6-3 21:24
你这个片段整体应该是带一个正电的,所以应该是偶数个电子,电荷和自旋多重度设为1 1
之后再进行优化,如 ...

前辈您好,感谢您的指导,我之前没注意设置正电荷与多重度的设置,现在我设为了1 1,用混合基组重新尝试。
h840473807 发表于 Post on 2018-6-3 21:24:11
你这个片段整体应该是带一个正电的,所以应该是偶数个电子,电荷和自旋多重度设为1 1
之后再进行优化,如果出现几何优化不收敛,再用http://sobereva.com/164 量子化学计算中帮助几何优化收敛的常用方法 这篇博文里的方法解决
不同计算等级下有的等级容易收敛,有的等级不容易收敛是很正常的事情。但你这个分子全用赝势基组肯定是不合适的(C,H,O,N应当用全电子基组),混合基组挺合适,如果嫌麻烦,可以用def2-SVP优化。
关于基组选取的看以下博文http://sobereva.com/336谈谈量子化学中基组的选择
http://sobereva.com/373谈谈赝势基组的选用

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2 我很赞同

查看全部评分 View all ratings

luktian 发表于 Post on 2018-6-3 14:40:00
让你变成回忆 发表于 2018-6-3 14:19
你把你的所有文件打包嘛。这么多~

您好前辈,我把6个附件打包上传在最后了
让你变成回忆 发表于 Post on 2018-6-3 14:19:23
你把你的所有文件打包嘛。这么多~

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-23 16:17 , Processed in 0.223411 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list