高斯优化第一激发态问题

南有乔木 · 发表于 Post on 2017-1-10 15:58:55

用高斯优化第一激发态时一直出错，关键是其他分子都能正常计算，就这一个分子，死活算不出来，大家帮我看一下是出了什么问题，麻烦大家了

X59-fasheopt.com (6.09 KB, 下载次数 Times of downloads: 56)

X59-fasheopt.log (546.27 KB, 下载次数 Times of downloads: 36)

小范范1989 · 发表于 Post on 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版本的。（今年预算没做，但是需要的话，可以买）大约多少钱呢？谢谢。

南有乔木 · 发表于 Post on 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 ...

非常感谢您的回答，主要是我后边还算了其他比这个分子大的激发态，都是没有任何问题的，就这一个不行，我也试了很多方法，最后结果都是这样的输出文件，实在计算不出来

。
新的高斯是可以求解的，至于价钱，因为都是我们老板添置的，我也不是很清楚诶。

chrischen1128 · 发表于 Post on 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 ...

高斯价格官网上都有、

chrischen1128 · 发表于 Post on 2017-1-10 21:35:38

scf=(maxcycle=300,
conver=7,xqc)去掉试试。

南有乔木 · 发表于 Post on 2017-1-10 22:00:46

chrischen1128 发表于 2017-1-10 21:35
scf=(maxcycle=300,
conver=7,xqc)去掉试试。

好的，我试试哈，谢谢哈

南有乔木 · 发表于 Post on 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

ggdh · 发表于 Post on 2017-1-10 22:54:20

优化激发态首先把nstates=10去掉。算这么多态，却只优化第一态。坑爹呢你这是。

liyuanhe211 · 发表于 Post on 2017-1-10 23:10:05

本帖最后由 liyuanhe211 于 2017-1-10 23:12 编辑

ggdh 发表于 2017-1-10 22:54
优化激发态首先把nstates=10去掉。算这么多态，却只优化第一态。坑爹呢你这是。

你这建议才坑爹。。。

Root 2 not converged, maximum delta is 0.535926330535614
No map to state 3
No map to state 4
You need to solve for more vectors in order to follow this state.
Error termination via Lnk1e in /public/software/g09/l914.exe at Sun Jan 8 21:04:36 2017.
Job cpu time: 5 days 0 hours 6 minutes 41.1 seconds.
File lengths (MBytes): RWF= 3460 Int= 0 D2E= 0 Chk= 413 Scr= 1

复制代码

他输出末尾如图，显然是在优化过程中发生了激发态顺序的翻转之类，应该“solve for more vectors”，即增加nstates的数量，你要求将nstates=10去掉相当于使用了nstates=3，与问题（至少是与其提示）相反。

liyuanhe211 · 发表于 Post on 2017-1-10 23:11:39

以当前优化的最后一个结构为初猜结构（如果合理），适当增加nstates的数量，然后对照下文去掉你输入文件中用的一堆垃圾关键词：

常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

sobereva · 发表于 Post on 2017-1-11 01:41:36

按照liyuanhe的回复处理。SCF一大串纯属多余
SCRF里面写read完全意义不明，应当去掉。没事不要用CPCM，明明默认的IEFPCM比它更好

如果增大了nstates还不行，把步长上限减小后再试。nstates太小或步长太大都有可能造成此类报错。

sobereva · 发表于 Post on 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 ...

肯定支持，这是新版本的关键性重点，其它相对来说倒都是次要

小范范1989 · 发表于 Post on 2017-1-11 08:43:11

南有乔木发表于 2017-1-10 22:24
这里有各个版本的价格http://www.gaussian.com/g_prod/prix.htm

恩，谢谢你。我看还是09.我们有09.我们想买一个最新的，看看能不能支持解析求解激发态频率计算

小范范1989 · 发表于 Post on 2017-1-11 08:45:06

sobereva 发表于 2017-1-11 01:42
肯定支持，这是新版本的关键性重点，其它相对来说倒都是次要

偶，谢谢sob老师。现在已经有卖新的了吗？我刚刚看那个老师回复的网址，我看还是09.谢谢sob老师

ggdh · 发表于 Post on 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优化的第一激发态不一样么？

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

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

评分 Rate

评分 Rate

浏览过的版块