高斯优化第一激发态问题

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

sobereva · 发表于 Post on 2018-9-21 20:43:48

zhangwei2713 发表于 2018-9-21 10:30
前辈你好！我看你的小*虫帖子上面写的是加上关键词td=direct后优化成功的，你也是试了在G09A版本上面去掉 ...

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

zhangwei2713 · 发表于 Post on 2018-9-21 10:30:08

ggdh 发表于 2017-1-11 10:23
去掉root=1确实可以，因为当时g09刚出的时候我就遇到了这个问题。首先想到的就是根据提示增加激发态数量 ...

前辈你好！我看你的小*虫帖子上面写的是加上关键词td=direct后优化成功的，你也是试了在G09A版本上面去掉root=1得到的了优化的激发态吗？

sobereva · 发表于 Post on 2017-1-17 11:12:11

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

新版本手册还没完整给出，现在还没说。

caofan_success · 发表于 Post on 2017-1-16 17:04:28

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

sob老师啥时候发文讲解一下新版高斯的各种变化

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

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

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

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

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

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

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

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

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

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

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

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

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-11 21:42:08

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

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

南有乔木 · 发表于 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:30:24

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

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

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

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

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

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

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

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

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

评分 Rate

浏览过的版块