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

计算化学公社

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

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

[复制链接]

21

帖子

0

威望

211

eV
积分
232

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 计算出错的信息

2104

帖子

8

威望

2538

eV
积分
4802

Level 6 (一方通行)

首席卖萌官

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

21

帖子

0

威望

211

eV
积分
232

Level 3 能力者

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

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

411

帖子

0

威望

890

eV
积分
1301

Level 4 (黑子)

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

21

帖子

0

威望

211

eV
积分
232

Level 3 能力者

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

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

2104

帖子

8

威望

2538

eV
积分
4802

Level 6 (一方通行)

首席卖萌官

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

411

帖子

0

威望

890

eV
积分
1301

Level 4 (黑子)

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

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

21

帖子

0

威望

211

eV
积分
232

Level 3 能力者

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

好的  谢谢你

82

帖子

0

威望

547

eV
积分
629

Level 4 (黑子)

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

21

帖子

0

威望

211

eV
积分
232

Level 3 能力者

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

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

9760

帖子

16

威望

1万

eV
积分
22360

管理员

公社社长

发表于 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,2号:466017436。用于讨论理论、计算化学,两个群讨论内容相同,可加入任意其一但不可都加入,申请信息必须注明研究方向,否则一概不批。

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

21

帖子

0

威望

211

eV
积分
232

Level 3 能力者

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

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

30

帖子

0

威望

217

eV
积分
247

Level 3 能力者

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

21

帖子

0

威望

211

eV
积分
232

Level 3 能力者

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

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

30

帖子

0

威望

217

eV
积分
247

Level 3 能力者

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

本版积分规则

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

GMT+8, 2017-7-24 18:39 , Processed in 0.109336 second(s), 24 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

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