请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 157|回复: 17

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

[复制链接]

20

帖子

0

威望

175

eV
积分
195

Level 3 能力者

发表于 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出错信息



这是用int=ultrafine 计算出错的信息

这是用int=ultrafine 计算出错的信息

2021

帖子

8

威望

2343

eV
积分
4524

Level 6 (一方通行)

首席卖萌官

发表于 2017-5-17 14:47:23 | 显示全部楼层
有过渡金属原子?那就换M06或者B3LYP泛函

20

帖子

0

威望

175

eV
积分
195

Level 3 能力者

 楼主| 发表于 2017-5-17 14:48:32 | 显示全部楼层
我本是个娃娃 发表于 2017-5-17 14:47
有过渡金属原子?那就换M06或者B3LYP泛函

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

345

帖子

0

威望

744

eV
积分
1089

Level 4 (黑子)

发表于 2017-5-17 14:56:45 | 显示全部楼层
试试scf=xqc

20

帖子

0

威望

175

eV
积分
195

Level 3 能力者

 楼主| 发表于 2017-5-17 15:00:48 | 显示全部楼层

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

2021

帖子

8

威望

2343

eV
积分
4524

Level 6 (一方通行)

首席卖萌官

发表于 2017-5-17 15:05:16 | 显示全部楼层
如果是做有机分子的几何结构,可以尝试PBE0

345

帖子

0

威望

744

eV
积分
1089

Level 4 (黑子)

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

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

20

帖子

0

威望

175

eV
积分
195

Level 3 能力者

 楼主| 发表于 2017-5-17 15:33:01 | 显示全部楼层
我本是个娃娃 发表于 2017-5-17 15:05
如果是做有机分子的几何结构,可以尝试PBE0

好的  谢谢你

38

帖子

0

威望

422

eV
积分
460

Level 3 能力者

发表于 2017-5-17 18:50:05 | 显示全部楼层
会不会是阳离子平衡结构与中性分子相差过大导致的?
尝试一下把scf和opt的收敛标准都loose下来,先初步优化下?

20

帖子

0

威望

175

eV
积分
195

Level 3 能力者

 楼主| 发表于 2017-5-17 19:41:24 | 显示全部楼层
Jasminer 发表于 2017-5-17 18:50
会不会是阳离子平衡结构与中性分子相差过大导致的?
尝试一下把scf和opt的收敛标准都loose下来,先初步优 ...

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

9119

帖子

16

威望

1万

eV
积分
20756

管理员

公社社长

发表于 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优化
思想家公社的门口Blog:http://sobereva.com
北京科音自然科学研究中心:http://www.keinsci.com
Multiwfn主页:http://sobereva.com/multiwfn
计算化学公社:http://bbs.keinsci.com
思想家公社QQ群1号:18616395(2000人已满),思想家公社QQ群2号:466017436(讨论计算化学为主,加入时必须注明研究方向)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

20

帖子

0

威望

175

eV
积分
195

Level 3 能力者

 楼主| 发表于 7 天前 | 显示全部楼层
sobereva 发表于 2017-5-17 21:54
你没用弥散函数,此时int=acc2e=12 SCF=novaracc SCF=noincfock基本不会解决问题
stable关键词和解决SCF不 ...

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

28

帖子

0

威望

179

eV
积分
207

Level 3 能力者

发表于 6 天前 | 显示全部楼层
我有时候会直接读取中性分子的波函数,把它作为阳离子的初猜。记得必须是同一个分子。不知道对你有没有用。

20

帖子

0

威望

175

eV
积分
195

Level 3 能力者

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

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

28

帖子

0

威望

179

eV
积分
207

Level 3 能力者

发表于 6 天前 | 显示全部楼层
不是的哦,是把中性分子的chk复制到阳离子要优化的chk,在阳离子输入文件里面加guess=check.然后提交作业就行了。
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

Archiver|手机版|小黑屋|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2017-5-25 22:29 , Processed in 0.105954 second(s), 24 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表