计算化学公社

标题: 高斯优化第一激发态问题 [打印本页]
作者
Author: 南有乔木    时间: 2017-1-10 15:58
标题: 高斯优化第一激发态问题
用高斯优化第一激发态时一直出错，关键是其他分子都能正常计算，就这一个分子，死活算不出来，大家帮我看一下是出了什么问题，麻烦大家了 (, 下载次数 Times of downloads: 56) (, 下载次数 Times of downloads: 36)


作者
Author: 小范范1989    时间: 2017-1-10 16:12
我的感觉哈：
#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版本的。（今年预算没做，但是需要的话，可以买）大约多少钱呢？谢谢。
作者
Author: 南有乔木    时间: 2017-1-10 16:21
小范范1989 发表于 2017-1-10 16:12
我的感觉哈：
#p td=(singlets,nstates=10,root=1) opt B3P86/6-31g* scf=(maxcycle=300,conver=7,xqc) sc ...

非常感谢您的回答，主要是我后边还算了其他比这个分子大的激发态，都是没有任何问题的，就这一个不行，我也试了很多方法，最后结果都是这样的输出文件，实在计算不出来。
新的高斯是可以求解的，至于价钱，因为都是我们老板添置的，我也不是很清楚诶。
作者
Author: chrischen1128    时间: 2017-1-10 21:30
小范范1989 发表于 2017-1-10 16:12
我的感觉哈：
#p td=(singlets,nstates=10,root=1) opt B3P86/6-31g* scf=(maxcycle=300,conver=7,xqc) sc ...

高斯价格官网上都有、
作者
Author: chrischen1128    时间: 2017-1-10 21:35
scf=(maxcycle=300,
conver=7,xqc)去掉试试。
作者
Author: 南有乔木    时间: 2017-1-10 22:00
chrischen1128 发表于 2017-1-10 21:35
scf=(maxcycle=300,
conver=7,xqc)去掉试试。

好的，我试试哈，谢谢哈
作者
Author: 南有乔木    时间: 2017-1-10 22:24
小范范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
作者
Author: ggdh    时间: 2017-1-10 22:54
优化激发态 首先把nstates=10去掉。算这么多态，却只优化第一态。坑爹呢你这是。
作者
Author: liyuanhe211    时间: 2017-1-10 23:10
本帖最后由 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，与问题（至少是与其提示）相反。
作者
Author: liyuanhe211    时间: 2017-1-10 23:11
以当前优化的最后一个结构为初猜结构（如果合理），适当增加nstates的数量，然后对照下文去掉你输入文件中用的一堆垃圾关键词：
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331
作者
Author: sobereva    时间: 2017-1-11 01:41
按照liyuanhe的回复处理。SCF一大串纯属多余
SCRF里面写read完全意义不明，应当去掉。没事不要用CPCM，明明默认的IEFPCM比它更好

如果增大了nstates还不行，把步长上限减小后再试。nstates太小或步长太大都有可能造成此类报错。
作者
Author: sobereva    时间: 2017-1-11 01:42
小范范1989 发表于 2017-1-10 16:12
我的感觉哈：
#p td=(singlets,nstates=10,root=1) opt B3P86/6-31g* scf=(maxcycle=300,conver=7,xqc) sc ...

肯定支持，这是新版本的关键性重点，其它相对来说倒都是次要
作者
Author: 小范范1989    时间: 2017-1-11 08:43
南有乔木 发表于 2017-1-10 22:24
这里有各个版本的价格http://www.gaussian.com/g_prod/prix.htm

恩，谢谢你。我看还是09.我们有09.我们想买一个最新的，看看能不能支持解析求解激发态频率计算
作者
Author: 小范范1989    时间: 2017-1-11 08:45
sobereva 发表于 2017-1-11 01:42
肯定支持，这是新版本的关键性重点，其它相对来说倒都是次要

偶，谢谢sob老师。现在已经有卖新的了吗？我刚刚看那个老师回复的网址，我看还是09.谢谢sob老师
作者
Author: ggdh    时间: 2017-1-11 09:35
本帖最后由 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优化的第一激发态不一样么？
作者
Author: sobereva    时间: 2017-1-11 10:05
ggdh 发表于 2017-1-11 09:35
这个“solve for more vectors”是老版本g09的bug，应该是把root=1去掉就能解决的。。难道nstates=10优化 ...


http://bbs.keinsci.com/thread-299-1-1.html
此贴8L有高斯官方回复

我未听说过有root=1去掉能解决的情况。root=1本身就是默认的，写不写都一样，写了也白写。
作者
Author: ggdh    时间: 2017-1-11 10:23
sobereva 发表于 2017-1-11 10:05
http://bbs.keinsci.com/thread-299-1-1.html
此贴8L有高斯官方回复

去掉root=1确实可以，因为当时g09刚出的时候我就遇到了这个问题。首先想到的就是根据提示增加激发态数量。但是发现没有用。后来去掉了root=1居然就行了。所以印象深刻。
网上一搜，居然还发现了我当年发的贴
http://muchong.com/html/200910/1633849.html
作者
Author: liyuanhe211    时间: 2017-1-11 10:23
ggdh 发表于 2017-1-11 09:35
这个“solve for more vectors”是老版本g09的bug，应该是把root=1去掉就能解决的。。难道nstates=10优化 ...

root=1去掉之后还是root=1，本来就是默认的，自己看看去掉、不去掉的Gaussian IOp路径有没有区别就知道了。再说我要优化第二激发态怎么办？也把root=2去掉？！

nstates不同当然不一样。激发态迭代有数值问题，关心第x个态，当然要计算比它多出一定数量的态，如x+5，才能保证x是准确的。经验上如果激发态不太“密集”确实一般x+5就够了，但如果分子比较复杂，就需要多很多态才能算准。

再不信你自己去找个有大共轭的片段、中间用个双键连上，用nstates=1和nstates=10比较一下激发能和优化后结构的差别就知道了。
作者
Author: ggdh    时间: 2017-1-11 10:40
本帖最后由 ggdh 于 2017-1-11 10:50 编辑
liyuanhe211 发表于 2017-1-11 10:23
root=1去掉之后还是root=1，本来就是默认的，自己看看去掉、不去掉的Gaussian IOp路径有没有区别就知道了 ...

你这都是理论上的，理论上我也知道root=1 加不加都一样。当年我是有血泪经验的。加root=1就出错，不加就能算，所以我才印象深刻。而且我也用这个方法解救过不少人。
不过现在争论这个意义不大，因为g09 c01版本以后这个bug去掉了（据说是一个内存分配的bug）。只有那些还在用老版本的人才会遇到这个问题。可以问问楼主用的啥版本.
第二个点待我测下。。测完了发个测试报告上来。
作者
Author: sobereva    时间: 2017-1-11 11:35
小范范1989 发表于 2017-1-11 08:45
偶，谢谢sob老师。现在已经有卖新的了吗？我刚刚看那个老师回复的网址，我看还是09.谢谢sob老师

没公开卖，但是据代理说，现在订购g09，之后会发货g16
作者
Author: sobereva    时间: 2017-1-11 11:37
如果是怀疑老版本bug触发"solve for more vectors"，先换成新版本再算一下就完了。
但这个并非只是bug导致的，这种错误本身是有对应的意义的，即新一步得到的态中找不到与上一步被优化的态的对应关系，从而没法对指定的态优化下去。抛开老版本bug导致的情况不谈，增大nstates和减小步长上限都是最正统的解决方法（建议先试减小步长上限的做法，因为当前态数已经不少了）。
作者
Author: 小范范1989    时间: 2017-1-11 12:00
sobereva 发表于 2017-1-11 11:35
没公开卖，但是据代理说，现在订购g09，之后会发货g16

偶，这样子，谢谢sob老师。
作者
Author: 南有乔木    时间: 2017-1-11 21:30
chrischen1128 发表于 2017-1-10 21:35
scf=(maxcycle=300,
conver=7,xqc)去掉试试。

已经试过了，还是同样的错误，算不出来
作者
Author: 南有乔木    时间: 2017-1-11 21:38
sobereva 发表于 2017-1-11 01:41
按照liyuanhe的回复处理。SCF一大串纯属多余
SCRF里面写read完全意义不明，应当去掉。没事不要用CPCM，明 ...

由于文献上用的是CPCM模型，不能改，我已经把没用的关键词去掉了，正在算哈，增大了nstates也算不出来
作者
Author: 南有乔木    时间: 2017-1-11 21:42
sobereva 发表于 2017-1-11 11:37
如果是怀疑老版本bug触发"solve for more vectors"，先换成新版本再算一下就完了。
但这个并非只是bug导致 ...

好的老师，大家提的建议我都会试试的，希望可以算出来
作者
Author: ggdh    时间: 2017-1-12 09:26
我测试了一下nstates的影响，顺便帮你算完了。我用的g09 E.01不会出错。
关键词分别是这两个。
# td=(singlets,nstates=10,root=1) opt B3P86/6-31g* scrf=(CPCM,read,solvent=dichloromethane)
# td=(singlets,nstates=1,root=1) opt B3P86/6-31g* scrf=(CPCM,read,solvent=dichloromethane)
最后的能量是：
X59-fasheopt_n10.log: Total Energy, E(TD-HF/TD-KS) =  -2537.90857326
X59-fasheopt_n1.log: Total Energy, E(TD-HF/TD-KS) =  -2537.90857333
可以看到只算一个态和算10个态的结果相差只有0.00000007Ha，基本认为就是数值误差了。
当然一个样本太少，能量一样也许只是巧合。我还得继续测几个不同种类的分子。
作者
Author: 南有乔木    时间: 2017-1-13 21:32
ggdh 发表于 2017-1-12 09:26
我测试了一下nstates的影响，顺便帮你算完了。我用的g09 E.01不会出错。
关键词分别是这两个。
# td=(sin ...

真的非常感谢，我会用你给的关键词再自己算一下的，我自己试了很多方法，还是出错，希望这次能够顺利算出来啊
作者
Author: 南有乔木    时间: 2017-1-13 21:49
ggdh 发表于 2017-1-12 09:26
我测试了一下nstates的影响，顺便帮你算完了。我用的g09 E.01不会出错。
关键词分别是这两个。
# td=(sin ...

老师，您发给我的log文件是不是加密了啊，我用UE打开是乱码的
作者
Author: ggdh    时间: 2017-1-14 08:10
南有乔木 发表于 2017-1-13 21:32
真的非常感谢，我会用你给的关键词再自己算一下的，我自己试了很多方法，还是出错，希望这次能够顺利算出 ...

你确定你gaussian版本正确么？
论坛有下载新版本的。
如果实在不想装新版，用td=direct关键词试试。我给你的压缩文件种第二个文件是好的。第一个文件不知道怎么变成乱码了。。。
作者
Author: 南有乔木    时间: 2017-1-15 09:25
ggdh 发表于 2017-1-14 08:10
你确定你gaussian版本正确么？
论坛有下载新版本的。
如果实在不想装新版，用td=direct关键词试试。我 ...

好的，我自己用你给的关键词计算了，已经算出来了，这个老问题终于解决了，非常感谢。
作者
Author: caofan_success    时间: 2017-1-16 17:04
sobereva 发表于 2017-1-11 01:42
肯定支持，这是新版本的关键性重点，其它相对来说倒都是次要

sob老师啥时候发文讲解一下新版高斯的各种变化
作者
Author: sobereva    时间: 2017-1-17 11:12
caofan_success 发表于 2017-1-16 17:04
sob老师啥时候发文讲解一下新版高斯的各种变化

新版本手册还没完整给出，现在还没说。
作者
Author: zhangwei2713    时间: 2018-9-21 10:30
ggdh 发表于 2017-1-11 10:23
去掉root=1确实可以，因为当时g09刚出的时候我就遇到了这个问题。首先想到的就是根据提示增加激发态数量 ...

前辈你好！我看你的小*虫帖子上面写的是加上关键词td=direct后优化成功的，你也是试了在G09A版本上面去掉root=1得到的了优化的激发态吗？
作者
Author: sobereva    时间: 2018-9-21 20:43
zhangwei2713 发表于 2018-9-21 10:30
前辈你好！我看你的小*虫帖子上面写的是加上关键词td=direct后优化成功的，你也是试了在G09A版本上面去掉 ...

TD本身默认就用direct，加这个没意义



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3