计算化学公社

标题: 求助：使用Multiwfn分析NTO时显示unable to locate 'excited state' label [打印本页]

作者
Author: bobo 时间: 2020-5-11 06:15
标题: 求助：使用Multiwfn分析NTO时显示unable to locate 'excited state' label
本帖最后由 bobo 于 2020-5-11 20:05 编辑

老师您好：

问题：基于Gaussian16的输出结果fchk文件，用Multiwfn做NTO分析, 一直显示error： unable to locate 'excited state' label 。
输入文件fchk具体计算信息：gaussian计算的命令行是#p b3lyp/gen TD(singlets, nstates=13) iop(3/33=4) pop=full scf=qc geom=connectivity integral=grid=ultrafine
pseudo=read 计算成功，得到fchk文件。

用fchk文件在Multiwfn中分析一直显示error：unable to locate 'excited state' label，附件是Multiwfn详细的输出文件。

请问老师这种情况该如何解决？

非常感谢

作者
Author: snljty 时间: 2020-5-11 09:33
写iop(3/33=4) pop=full干嘛，一个输出双电子积分一个输出所有轨道布居，和NTO分析没有任何关系，还会让输出文件变得臃肿不堪。
反倒是重要的iop(9/40=4)或者最高的iop(9/40=5)没写，输出的组态系数阈值太大，结果精度会明显不够用。
上传你的Multiwfn的settings.ini，有可能是iloadGaugeom之类的设置的问题。
G16对DFT默认就是Integral=UltrafineGrid，除非你Default.Rou(te)给覆盖了，否则没必要写，而且大多数情况Integral=FineGrid对B3LYP已经足够。

作者
Author: sobereva 时间: 2020-5-11 13:21
按照2L说的，把乱七八糟关键词都去掉重新算，如果还不行，上传压缩后的输出文件

PS：不要轻易用scf=qc，不仅非常慢而且很危险，这里说了
解决SCF不收敛问题的方法
http://sobereva.com/61

作者
Author: bobo 时间: 2020-5-11 15:47

snljty 发表于 2020-5-11 09:33
写iop(3/33=4) pop=full干嘛，一个输出双电子积分一个输出所有轨道布居，和NTO分析没有任何关系，还会让输 ...

好的，非常感谢，我再重新试一下，如果有问题再向您请教。

作者
Author: bobo 时间: 2020-5-11 15:47

sobereva 发表于 2020-5-11 13:21
按照2L说的，把乱七八糟关键词都去掉重新算，如果还不行，上传压缩后的输出文件

PS：不要轻易用scf=qc， ...

嗯嗯，非常感谢，我再重新试一下，如果有问题再向您请教。

作者
Author: bobo 时间: 2020-5-11 15:53
本帖最后由 bobo 于 2020-5-11 15:57 编辑

snljty 发表于 2020-5-11 09:33
写iop(3/33=4) pop=full干嘛，一个输出双电子积分一个输出所有轨道布居，和NTO分析没有任何关系，还会让输 ...

您好，settings.ini 文件已经上传了，请问iloadgaugeom设置有问题的话怎么进行修改？

作者
Author: snljty 时间: 2020-5-11 16:13

bobo 发表于 2020-5-11 15:53
您好，settings.ini 文件已经上传了，请问iloadgaugeom设置有问题的话怎么进行修改？

看来不是这个的问题，请精简一下你的关键词后重试。

作者
Author: bobo 时间: 2020-5-11 18:52
本帖最后由 bobo 于 2020-5-11 19:02 编辑

snljty 发表于 2020-5-11 16:13
看来不是这个的问题，请精简一下你的关键词后重试。

#p b3lyp/gen TD(triplets, nstates=13) iop(9/40=4) pop=full geom=connectivity pseudo=read

您好，改成上述关键词，计算成功后(三线态先算出来了)，用Multiwfn分析还是显示error: unable to locate 'excited state' label

作者
Author: snljty 时间: 2020-5-11 19:07

bobo 发表于 2020-5-11 18:52
#p b3lyp/gen TD(triplets, nstates=13) iop(9/40=4) pop=full geom=connectivity pseudo=read

您好 ...

既然你都会用Multiwfn了，那Pop=Full就没有任何意义了。用Multiwfn作布居分析方便的多。先删掉，还不行的话，建议上传文件请大家帮忙分析。

作者
Author: bobo 时间: 2020-5-11 19:13

snljty 发表于 2020-5-11 19:07
既然你都会用Multiwfn了，那Pop=Full就没有任何意义了。用Multiwfn作布居分析方便的多。先删掉，还不行的 ...

嗯嗯，好的，我再试一下，我正在上传这一次的结果，分子有点大，上传到现在也没上传完成

作者
Author: bobo 时间: 2020-5-11 20:04
本帖最后由 bobo 于 2020-5-11 20:06 编辑

sobereva 发表于 2020-5-11 13:21
按照2L说的，把乱七八糟关键词都去掉重新算，如果还不行，上传压缩后的输出文件

PS：不要轻易用scf=qc， ...

老师，您好，我简化关键词后还是不行，依然是error。我fchk文件压缩之后还有25Mb,网站不让我上传，只能小于9.8Mb才能上传，请问还有其他方法能上传麻烦您帮忙看看吗？

作者
Author: sobereva 时间: 2020-5-11 23:16

bobo 发表于 2020-5-11 20:04
老师，您好，我简化关键词后还是不行，依然是error。我fchk文件压缩之后还有25Mb,网站不让我上传，只能小 ...

传网盘啊

作者
Author: bobo 时间: 2020-5-12 00:22

sobereva 发表于 2020-5-11 23:16
传网盘啊

链接：https://pan.baidu.com/s/1PAJFCh6viRD6X1Quq9qTbA
提取码：6fth
老师，您好，我把fchk和out文件上传到网盘了，麻烦您帮忙看一下。我用Multiwfn分析一直显示error

作者
Author: bobo 时间: 2020-5-12 00:23

snljty 发表于 2020-5-11 19:07
既然你都会用Multiwfn了，那Pop=Full就没有任何意义了。用Multiwfn作布居分析方便的多。先删掉，还不行的 ...

链接：https://pan.baidu.com/s/1PAJFCh6viRD6X1Quq9qTbA
提取码：6fth

您好，老师。我将pop=full删掉了，Multiwfn分析还是显示error.上面是fchk文件，麻烦帮忙看一下问题在哪里。谢谢

作者
Author: sobereva 时间: 2020-5-12 03:37
我用目前最新版本，没有遇到问题，输出信息见附件：
(, 下载次数 Times of downloads: 3)

估计是你的Multiwfn版本太老所致。遇到问题时，总应当先尝试最新版本。

作者
Author: bobo 时间: 2020-5-12 04:57

sobereva 发表于 2020-5-12 03:37
我用目前最新版本，没有遇到问题，输出信息见附件：

感谢老师。看了您的输出信息，我找到原因，应该是out文件也要用的，我每次都是用空白文档让它读取了。非常感谢。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)