计算化学公社

标题: 优化结构总出错 [打印本页]

作者
Author: wx91aaa 时间: 2018-9-30 09:01
标题: 优化结构总出错
本帖最后由 wx91aaa 于 2018-9-30 09:01 编辑

各位老师好！我按照文献里的这个结构（图1）进行计算，方法是  #p opt freq b3lyp/6-31g nosym optcyc=500 scfcyc=500。我的初始结构是这样的（图2）。最终的结果是收敛失败 >>>>>>>>>> Convergence criterion not met. SCF Done:  E(UB3LYP) =  -2461.73071007    A.U. after  501 cycles
         Convg  = 0.3058D-05          -V/T =  2.0023
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 2.5323 S= 1.1680
<L.S>= 0.000000000000E+00
KE= 2.455963348504D+03 PE=-1.187783023688D+04 EE= 3.917609632804D+03
Annihilation of the first spin contaminant:
S**2 before annihilation    2.5323, after    3.1556
Convergence failure -- run terminated.
Error termination via Lnk1e in d:\G09W\l502.exe at Sat Sep 29 19:38:51 2018.
Job cpu time:  0 days  8 hours 22 minutes 36.0 seconds.
File lengths (MBytes):  RWF= 137 Int=    0 D2E=    0 Chk=    15 Scr=    1

我调了好多次结构，都是这样的错误，请问是我最初构建的分子结构就有问题吗？请大神帮忙看看。十分感谢！！！

作者
Author: Adair 时间: 2018-9-30 11:10
建议看老师的博文: 《解决SCF不收敛问题的方法》http://sobereva.com/61 结合自己的体系选择适当的解决方法。另外你的关键词里面 optcyc=500 scfcyc=500对解决问题几乎没用，一般都不能用。

作者
Author: sobereva 时间: 2018-9-30 11:22
你图2里氮原子的位置严重不合理，哪能那样摆，参考吡啶结构。结构严重不合理极容易导致SCF不收敛

另外，就这么点体系还用6-31G完全不能接受，尤其是还牵扯杂原子的情况，仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336

作者
Author: wx91aaa 时间: 2018-9-30 11:24

Adair 发表于 2018-9-30 11:10
建议看老师的博文: 《解决SCF不收敛问题的方法》http://sobereva.com/61 结合自己的体系选择适当的解决方 ...

好的。谢谢！

作者
Author: wx91aaa 时间: 2018-9-30 11:24

sobereva 发表于 2018-9-30 11:22
你图2里氮原子的位置严重不合理，哪能那样摆，参考吡啶结构。结构严重不合理极容易导致SCF不收敛

另外， ...

好的。谢谢sob老师！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)