高斯优化第一激发态问题

sobereva · 发表于 Post on 2017-1-11 10:05:01

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本身就是默认的，写不写都一样，写了也白写。

ggdh · 发表于 Post on 2017-1-11 10:23:12

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

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

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比较一下激发能和优化后结构的差别就知道了。

ggdh · 发表于 Post on 2017-1-11 10:40:48

本帖最后由 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）。只有那些还在用老版本的人才会遇到这个问题。可以问问楼主用的啥版本.
第二个点待我测下。。测完了发个测试报告上来。

sobereva · 发表于 Post on 2017-1-11 11:35:05

小范范1989 发表于 2017-1-11 08:45
偶，谢谢sob老师。现在已经有卖新的了吗？我刚刚看那个老师回复的网址，我看还是09.谢谢sob老师

没公开卖，但是据代理说，现在订购g09，之后会发货g16

sobereva · 发表于 Post on 2017-1-11 11:37:40

如果是怀疑老版本bug触发"solve for more vectors"，先换成新版本再算一下就完了。
但这个并非只是bug导致的，这种错误本身是有对应的意义的，即新一步得到的态中找不到与上一步被优化的态的对应关系，从而没法对指定的态优化下去。抛开老版本bug导致的情况不谈，增大nstates和减小步长上限都是最正统的解决方法（建议先试减小步长上限的做法，因为当前态数已经不少了）。

小范范1989 · 发表于 Post on 2017-1-11 12:00:57

sobereva 发表于 2017-1-11 11:35
没公开卖，但是据代理说，现在订购g09，之后会发货g16

偶，这样子，谢谢sob老师。

南有乔木 · 发表于 Post on 2017-1-11 21:30:24

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

已经试过了，还是同样的错误，算不出来

南有乔木 · 发表于 Post on 2017-1-11 21:38:36

sobereva 发表于 2017-1-11 01:41
按照liyuanhe的回复处理。SCF一大串纯属多余
SCRF里面写read完全意义不明，应当去掉。没事不要用CPCM，明 ...

由于文献上用的是CPCM模型，不能改，我已经把没用的关键词去掉了，正在算哈，增大了nstates也算不出来

南有乔木 · 发表于 Post on 2017-1-11 21:42:08

sobereva 发表于 2017-1-11 11:37
如果是怀疑老版本bug触发"solve for more vectors"，先换成新版本再算一下就完了。
但这个并非只是bug导致 ...

好的老师，大家提的建议我都会试试的，希望可以算出来

ggdh · 发表于 Post on 2017-1-12 09:26:13

我测试了一下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，基本认为就是数值误差了。
当然一个样本太少，能量一样也许只是巧合。我还得继续测几个不同种类的分子。

南有乔木 · 发表于 Post on 2017-1-13 21:32:26

ggdh 发表于 2017-1-12 09:26
我测试了一下nstates的影响，顺便帮你算完了。我用的g09 E.01不会出错。
关键词分别是这两个。
# td=(sin ...

真的非常感谢，我会用你给的关键词再自己算一下的，我自己试了很多方法，还是出错，希望这次能够顺利算出来啊

南有乔木 · 发表于 Post on 2017-1-13 21:49:57

ggdh 发表于 2017-1-12 09:26
我测试了一下nstates的影响，顺便帮你算完了。我用的g09 E.01不会出错。
关键词分别是这两个。
# td=(sin ...

老师，您发给我的log文件是不是加密了啊，我用UE打开是乱码的

ggdh · 发表于 Post on 2017-1-14 08:10:02

南有乔木发表于 2017-1-13 21:32
真的非常感谢，我会用你给的关键词再自己算一下的，我自己试了很多方法，还是出错，希望这次能够顺利算出 ...

你确定你gaussian版本正确么？
论坛有下载新版本的。
如果实在不想装新版，用td=direct关键词试试。我给你的压缩文件种第二个文件是好的。第一个文件不知道怎么变成乱码了。。。

南有乔木 · 发表于 Post on 2017-1-15 09:25:08

ggdh 发表于 2017-1-14 08:10
你确定你gaussian版本正确么？
论坛有下载新版本的。
如果实在不想装新版，用td=direct关键词试试。我 ...

好的，我自己用你给的关键词计算了，已经算出来了，这个老问题终于解决了，非常感谢。

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

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

评分 Rate

浏览过的版块