计算化学公社

标题: MP2/6-311G(d,p)结构优化失败案例 [打印本页]

作者
Author: 尚艳磊 时间: 2017-12-5 14:33
标题: MP2/6-311G(d,p)结构优化失败案例
我找到了过渡态，并且通过IRC得到了PC和RC。PC的结构优化可以完成，但是RC的结构优化一直就是失败。用的都是MP2/6-311G（d,p）基组。
我贴一下失败的描述哈。
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UHF) =  -358.186283030    A.U. after  129 cycles
         Convg  = 0.4801D-05          -V/T =  2.0009
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7951 S= 0.5223
<L.S>= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation    0.7951, after    0.7513
Convergence failure -- run terminated.
Error termination via Lnk1e in E:\Shang\Gaussian\G09W\l502.exe at Fri Nov 24 22:01:57 2017.
Job cpu time:  0 days  0 hours 55 minutes  2.0 seconds.
File lengths (MBytes):  RWF= 5396 Int=    0 D2E=    0 Chk=    4 Scr=    1

是的，收敛达不到。我试过了scf=qc，降低收敛精度。都是无法正常结束的。用其他的密度泛函，B3LYP/6-311G(d,p)可以优化。
我看到过一些文献，说有的情况就是优化不了。
想请教一下大家，这种情况到底是什么导致的？
另外6-311和6-311G有什么区别吗
谢谢大家！

作者
Author: 我本是个娃娃 时间: 2017-12-5 15:24
M06-2X对积分格点很敏感，君可以添加int=ultrafine尝试一下

作者
Author: sobereva 时间: 2017-12-5 17:03
用MP2是极为不划算的，特别是开壳层体系，常规MP2表现渣，完全费力不讨好
如果死活非要用MP2，解决方法看
解决SCF不收敛问题的方法
http://sobereva.com/61

“说有的情况就是优化不了”此说法明显错误，“优化不了”完全意义不明。只要在当前计算级别，势能面存在对应的极小点，就不可能“优化不了”，碰到数值问题尝试相应解决就完了。

6-311是6-311G的错误写法

作者
Author: 尚艳磊 时间: 2017-12-6 09:04

我本是个娃娃发表于 2017-12-5 15:24
M06-2X对积分格点很敏感，君可以添加int=ultrafine尝试一下

谢谢大家~

作者
Author: 尚艳磊 时间: 2017-12-6 09:05

sobereva 发表于 2017-12-5 17:03
用MP2是极为不划算的，特别是开壳层体系，常规MP2表现渣，完全费力不讨好
如果死活非要用MP2，解决方法看
...

谢谢sob老师，我试一下。因为这半个月尝试了很多方法都不行。难过至极，主要是MP2算了很多了，就那么一两个算不出来

作者
Author: gfcjyb 时间: 2017-12-7 12:33
6-311G里的G是指高斯型函数。。。6-311G表示内层用一个STO-6G，价层用STO-3G，STO-1G和STO-1G三个基来描述。

作者
Author: sobereva 时间: 2017-12-7 17:55

gfcjyb 发表于 2017-12-7 12:33
6-311G里的G是指高斯型函数。。。6-311G表示内层用一个STO-6G，价层用STO-3G，STO-1G和STO-1G三个基来描述 ...

“6-311G表示内层用一个STO-6G，价层用STO-3G，STO-1G和STO-1G三个基来描述”这种说法是错误的，记得某本低水平的中文书里也这么写。你看过6-311G的原文就会知道，构建思想和STO-nG完全不同。
正确描述是：每个内层原子轨道靠一个收缩度为6的基函数描述，每个价层原子轨道靠一个收缩度为3的基函数和两个未收缩的基函数描述。

作者
Author: 尚艳磊 时间: 2017-12-7 18:26
我昨天脑子一抽微调了一下RC的结构，其实就是微调了一个键长，只做opt竟然收敛了。
但是计算频率的时候有虚频-3.12，说明没有到最优结构嘛。
我在之前优化的基础上继续优化结构又不能收敛了。。。

痛并快乐。。。

作者
Author: 尚艳磊 时间: 2017-12-8 20:58

sobereva 发表于 2017-12-7 17:55
“6-311G表示内层用一个STO-6G，价层用STO-3G，STO-1G和STO-1G三个基来描述”这种说法是错误的，记得某 ...

sob 老师您好！
我试过了您介绍的关于SCF不收敛的问题解决方法。还是不不能够收敛，我观察了能量也是逐渐平缓的（scf=qc），最后都超过步数了还是没结束。

我有一个疑问，我觉得就像您说的势能面存在对应的极小点，应该不可能优化不了的。密度泛函的方法确实成功了，但是MP2的方法却不行。也不是一个很大的体系，13个原子，双分子反应的RC。这到底会是什么原因呢？

谢谢sob老师之前的答疑，虽然没有解决，但也了解到了促进SCF收敛的一些方法及其原因。感觉自己掌握的知识还不够多，希望前辈们能够多多指点迷津。

作者
Author: sobereva 时间: 2017-12-9 01:25

尚艳磊发表于 2017-12-8 20:58
sob 老师您好！
我试过了您介绍的关于SCF不收敛的问题解决方法。还是不不能够收敛，我观察了能量也 ...

不要问“优化不了”，要交代优化具体出现了什么问题，什么导致的“优化不了”这个现象

那个帖子已经把所有可能解决SCF不收敛的方法全都归纳了，没解决说明还没有结合实际情况恰当使用帖子里的方法。

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