计算化学公社

标题: 优化阳离子结构一直L502报错 [打印本页]

作者
Author: 盐城人790 时间: 2017-5-17 14:45
标题: 优化阳离子结构一直L502报错
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出错信息

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

作者
Author: 盐城人790 时间: 2017-5-17 14:48

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

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

作者
Author: xujian 时间: 2017-5-17 14:56
试试scf=xqc

作者
Author: 盐城人790 时间: 2017-5-17 15:00

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

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

作者
Author: 我本是个娃娃 时间: 2017-5-17 15:05
如果是做有机分子的几何结构，可以尝试PBE0

作者
Author: xujian 时间: 2017-5-17 15:26

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

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

作者
Author: 盐城人790 时间: 2017-5-17 15:33

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

好的谢谢你

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

作者
Author: 盐城人790 时间: 2017-5-17 19:41

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

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

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

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

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

作者
Author: 盐城人790 时间: 2017-5-18 08:34

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

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

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

作者
Author: 盐城人790 时间: 2017-5-19 08:30

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

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

作者
Author: Flower 时间: 2017-5-19 08:33
不是的哦，是把中性分子的chk复制到阳离子要优化的chk,在阳离子输入文件里面加guess=check.然后提交作业就行了。

作者
Author: 盐城人790 时间: 2017-5-19 08:47

Flower 发表于 2017-5-19 08:33
不是的哦，是把中性分子的chk复制到阳离子要优化的chk,在阳离子输入文件里面加guess=check.然后提交作业就 ...

阳离子的chk 文件不是提交任务时才会出现的哈？请问是如何将中性的chk，复制到阳离子要优化的chkl里面的哈？谢谢

作者
Author: 盐城人790 时间: 2017-5-19 09:13

Flower 发表于 2017-5-19 08:33
不是的哦，是把中性分子的chk复制到阳离子要优化的chk,在阳离子输入文件里面加guess=check.然后提交作业就 ...

你好具体的方法我明白了非常感谢你

作者
Author: Flower 时间: 2017-5-19 12:18
不客气。会了就好。

作者
Author: 风玲儿 时间: 2017-9-20 09:18

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

老师，对于含有过渡金属的体系，用的B3LYP泛函，尝试了老师您说的几种方法，依然出现l502错误

作者
Author: sobereva 时间: 2017-9-20 11:31

风玲儿发表于 2017-9-20 09:18
老师，对于含有过渡金属的体系，用的B3LYP泛函，尝试了老师您说的几种方法，依然出现l502错误

所有能解决SCF不收敛的办法没有没被我那篇博文囊括、提到的，只能说还没恰当使用它们

作者
Author: new-hand 时间: 2022-3-31 20:29

Flower 发表于 2017-5-19 08:33
不是的哦，是把中性分子的chk复制到阳离子要优化的chk,在阳离子输入文件里面加guess=check.然后提交作业就 ...

不好意思，不是很明白你说的chk是要复制到哪里，可以截图解释一下吗？谢谢！

作者
Author: hhx521 时间: 2023-11-9 08:28

new-hand 发表于 2022-3-31 20:29
不好意思，不是很明白你说的chk是要复制到哪里，可以截图解释一下吗？谢谢！

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