优化阳离子结构一直L502报错

盐城人790 · 发表于 Post on 2017-5-17 14:45:51

SOB 老师好：我用 # opt freq m062x/6-31g(d,p) 的方法优化一个中性分子使其带一个正电荷多重度为2 试过SOB 老师博文里的 int=ultrafine int=acc2e=12 SCF=novaracc SCF=noincfock 以及改变初猜， SCF=vshift=400 解决不收敛的方法都一会就L502报错，同一个中性分子不带电荷优化可以正常结束。请问老师这种情况下如何才能继续优化。（分子里具有较多数量的氧）
下图是用int=ultrafine 计算时的L502出错信息

我本是个娃娃 · 发表于 Post on 2017-5-17 14:47:23

有过渡金属原子？那就换M06或者B3LYP泛函

盐城人790 · 发表于 Post on 2017-5-17 14:48:32

我本是个娃娃发表于 2017-5-17 14:47
有过渡金属原子？那就换M06或者B3LYP泛函

你好分子里面只有C H O 没有过渡金属呢不知为何一直在报错谢谢你的解答

xujian · 发表于 Post on 2017-5-17 14:56:45

试试scf=xqc

盐城人790 · 发表于 Post on 2017-5-17 15:00:48

xujian 发表于 2017-5-17 14:56
试试scf=xqc

这种方法我也在试还挺慢的请问SOB 老师博文里说的加stable关键词看波函数的稳定性具体看哪个能量值哈？

我本是个娃娃 · 发表于 Post on 2017-5-17 15:05:16

如果是做有机分子的几何结构，可以尝试PBE0

xujian · 发表于 Post on 2017-5-17 15:26:50

盐城人790 发表于 2017-5-17 15:00
这种方法我也在试还挺慢的请问SOB 老师博文里说的加stable关键词看波函数的稳定性具体看哪个能量值哈？

我猜想这种方法优化好之后，可以保存优化好的结构，去掉关键词再优化一次，看是否可行，仅为个人臆测

盐城人790 · 发表于 Post on 2017-5-17 15:33:01

我本是个娃娃发表于 2017-5-17 15:05
如果是做有机分子的几何结构，可以尝试PBE0

好的谢谢你

Jasminer · 发表于 Post on 2017-5-17 18:50:05

会不会是阳离子平衡结构与中性分子相差过大导致的？
尝试一下把scf和opt的收敛标准都loose下来，先初步优化下？

盐城人790 · 发表于 Post on 2017-5-17 19:41:24

Jasminer 发表于 2017-5-17 18:50
会不会是阳离子平衡结构与中性分子相差过大导致的？
尝试一下把scf和opt的收敛标准都loose下来，先初步优 ...

采用 opt=loose 和SCF=loose 的指令吗？好的我试试看谢谢你的解答

sobereva · 发表于 Post on 2017-5-17 21:54:32

你没用弥散函数，此时int=acc2e=12 SCF=novaracc SCF=noincfock基本不会解决问题
stable关键词和解决SCF不收敛完全无关，此文根本没说要加stable
解决SCF不收敛问题的方法
http://sobereva.com/61

非过渡金属体系用vshift也没什么用

先用其它泛函，比如B3LYP优化，之后用最后的结构，读取其收敛的波函数作为初猜再继续用M06-2X优化

盐城人790 · 发表于 Post on 2017-5-18 08:34:15

sobereva 发表于 2017-5-17 21:54
你没用弥散函数，此时int=acc2e=12 SCF=novaracc SCF=noincfock基本不会解决问题
stable关键词和解决SCF不 ...

好的非常感谢SOB 老师的解答我试试看

Flower · 发表于 Post on 2017-5-19 08:28:05

我有时候会直接读取中性分子的波函数，把它作为阳离子的初猜。记得必须是同一个分子。不知道对你有没有用。

盐城人790 · 发表于 Post on 2017-5-19 08:30:37

Flower 发表于 2017-5-19 08:28
我有时候会直接读取中性分子的波函数，把它作为阳离子的初猜。记得必须是同一个分子。不知道对你有没有用。

好的就是将中性分子的out 文件保存成gjf 改变电荷多重度然后加上 guess=read 吗？谢谢你

Flower · 发表于 Post on 2017-5-19 08:33:14

不是的哦，是把中性分子的chk复制到阳离子要优化的chk,在阳离子输入文件里面加guess=check.然后提交作业就行了。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 优化阳离子结构一直L502报错