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

计算化学公社

 找回密码
 现在注册!
查看: 1001|回复: 31

[Gaussian/gview] 高斯优化第一激发态问题

[复制链接]

23

帖子

0

威望

31

eV
积分
54

Level 2 能力者

发表于 2017-1-10 15:58:55 | 显示全部楼层 |阅读模式
用高斯优化第一激发态时一直出错,关键是其他分子都能正常计算,就这一个分子,死活算不出来,大家帮我看一下是出了什么问题,麻烦大家了 X59-fasheopt.com (6.09 KB, 下载次数: 23)

863

帖子

0

威望

2252

eV
积分
3115

Level 5 (御坂)

发表于 2017-1-10 16:12:37 | 显示全部楼层
我的感觉哈:
#p td=(singlets,nstates=10,root=1) opt B3P86/6-31g* scf=(maxcycle=300,conver=7,xqc) scrf=(CPCM,read,solvent=dichloromethane)
1:你的分子不小,优化激发态,表头越简单越好:opt td scrf=(CPCM,read,solvent=dichloromethane)
2:可以尝试先采用小机组,比如说3-21g等,先找到差不多结构,然后再在此结构上换为你要的机组。
3:我平日也计算激发态,特别是频率,那就更慢了。
PS:不知道新的高斯支不支持解析求解激发态频率。如果能的话,我们想买一个linux版本的。(今年预算没做,但是需要的话,可以买)大约多少钱呢?谢谢。
有梦想,才有实现。奋斗吧,骄傲的少年。

23

帖子

0

威望

31

eV
积分
54

Level 2 能力者

 楼主| 发表于 2017-1-10 16:21:37 | 显示全部楼层
小范范1989 发表于 2017-1-10 16:12
我的感觉哈:
#p td=(singlets,nstates=10,root=1) opt B3P86/6-31g* scf=(maxcycle=300,conver=7,xqc) sc ...

非常感谢您的回答,主要是我后边还算了其他比这个分子大的激发态,都是没有任何问题的,就这一个不行,我也试了很多方法,最后结果都是这样的输出文件,实在计算不出来
新的高斯是可以求解的,至于价钱,因为都是我们老板添置的,我也不是很清楚诶。

253

帖子

3

威望

1103

eV
积分
1416

科音成员

发表于 2017-1-10 21:30:57 | 显示全部楼层
小范范1989 发表于 2017-1-10 16:12
我的感觉哈:
#p td=(singlets,nstates=10,root=1) opt B3P86/6-31g* scf=(maxcycle=300,conver=7,xqc) sc ...

高斯价格官网上都有、
一花一世界,一叶一追寻。一曲一场叹,一生为一人

253

帖子

3

威望

1103

eV
积分
1416

科音成员

发表于 2017-1-10 21:35:38 | 显示全部楼层
scf=(maxcycle=300,
conver=7,xqc)去掉试试。
一花一世界,一叶一追寻。一曲一场叹,一生为一人

23

帖子

0

威望

31

eV
积分
54

Level 2 能力者

 楼主| 发表于 2017-1-10 22:00:46 | 显示全部楼层
chrischen1128 发表于 2017-1-10 21:35
scf=(maxcycle=300,
conver=7,xqc)去掉试试。

好的,我试试哈,谢谢哈

23

帖子

0

威望

31

eV
积分
54

Level 2 能力者

 楼主| 发表于 2017-1-10 22:24:23 | 显示全部楼层
小范范1989 发表于 2017-1-10 16:12
我的感觉哈:
#p td=(singlets,nstates=10,root=1) opt B3P86/6-31g* scf=(maxcycle=300,conver=7,xqc) sc ...

这里有各个版本的价格http://www.gaussian.com/g_prod/prix.htm

299

帖子

6

威望

902

eV
积分
1321

Level 4 (黑子)

发表于 2017-1-10 22:54:20 | 显示全部楼层
优化激发态 首先把nstates=10去掉。算这么多态,却只优化第一态。坑爹呢你这是。

1816

帖子

19

威望

5874

eV
积分
8070

Level 6 (一方通行)

发表于 2017-1-10 23:10:05 | 显示全部楼层
本帖最后由 liyuanhe211 于 2017-1-10 23:12 编辑
ggdh 发表于 2017-1-10 22:54
优化激发态 首先把nstates=10去掉。算这么多态,却只优化第一态。坑爹呢你这是。

你这建议才坑爹。。。
  1. Root      2 not converged, maximum delta is    0.535926330535614
  2. No map to state      3
  3. No map to state      4
  4. You need to solve for more vectors in order to follow this state.
  5. Error termination via Lnk1e in /public/software/g09/l914.exe at Sun Jan  8 21:04:36 2017.
  6. Job cpu time:  5 days  0 hours  6 minutes 41.1 seconds.
  7. File lengths (MBytes):  RWF=   3460 Int=      0 D2E=      0 Chk=    413 Scr=      1
复制代码
他输出末尾如图,显然是在优化过程中发生了激发态顺序的翻转之类,应该“solve for more vectors”,即增加nstates的数量,你要求将nstates=10去掉相当于使用了nstates=3,与问题(至少是与其提示)相反。

评分

参与人数 3eV +8 收起 理由
哈哈宝 + 2 我很赞同
冰释之川 + 5 我很赞同
sobereva + 1 我很赞同

查看全部评分

1816

帖子

19

威望

5874

eV
积分
8070

Level 6 (一方通行)

发表于 2017-1-10 23:11:39 | 显示全部楼层
以当前优化的最后一个结构为初猜结构(如果合理),适当增加nstates的数量,然后对照下文去掉你输入文件中用的一堆垃圾关键词:
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

评分

参与人数 1eV +1 收起 理由
sobereva + 1

查看全部评分

9938

帖子

16

威望

1万

eV
积分
22815

管理员

公社社长

发表于 2017-1-11 01:41:36 | 显示全部楼层
按照liyuanhe的回复处理。SCF一大串纯属多余
SCRF里面写read完全意义不明,应当去掉。没事不要用CPCM,明明默认的IEFPCM比它更好

如果增大了nstates还不行,把步长上限减小后再试。nstates太小或步长太大都有可能造成此类报错。
思想家公社的门口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!

9938

帖子

16

威望

1万

eV
积分
22815

管理员

公社社长

发表于 2017-1-11 01:42:36 | 显示全部楼层
小范范1989 发表于 2017-1-10 16:12
我的感觉哈:
#p td=(singlets,nstates=10,root=1) opt B3P86/6-31g* scf=(maxcycle=300,conver=7,xqc) sc ...

肯定支持,这是新版本的关键性重点,其它相对来说倒都是次要
思想家公社的门口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!

863

帖子

0

威望

2252

eV
积分
3115

Level 5 (御坂)

发表于 2017-1-11 08:43:11 | 显示全部楼层

恩,谢谢你。我看还是09.我们有09.我们想买一个最新的,看看能不能支持解析求解激发态频率计算
有梦想,才有实现。奋斗吧,骄傲的少年。

863

帖子

0

威望

2252

eV
积分
3115

Level 5 (御坂)

发表于 2017-1-11 08:45:06 | 显示全部楼层
sobereva 发表于 2017-1-11 01:42
肯定支持,这是新版本的关键性重点,其它相对来说倒都是次要

偶,谢谢sob老师。现在已经有卖新的了吗?我刚刚看那个老师回复的网址,我看还是09.谢谢sob老师
有梦想,才有实现。奋斗吧,骄傲的少年。

299

帖子

6

威望

902

eV
积分
1321

Level 4 (黑子)

发表于 2017-1-11 09:35:28 | 显示全部楼层
本帖最后由 ggdh 于 2017-1-11 09:43 编辑
liyuanhe211 发表于 2017-1-10 23:10
你这建议才坑爹。。。
他输出末尾如图,显然是在优化过程中发生了激发态顺序的翻转之类,应该“solve fo ...

这个“solve for more vectors”是老版本g09的bug,应该是把root=1去掉就能解决的。。难道nstates=10优化的第一激发态和nstates=3优化的第一激发态不一样么?
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2017-8-17 03:23 , Processed in 0.107192 second(s), 28 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

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