计算化学公社

标题: Gaussian16计算振动分辨的电子光谱出错?ERROR: Low progression after class 2. [打印本页]

作者
Author:
lindlar    时间: 2021-5-30 22:57
标题: Gaussian16计算振动分辨的电子光谱出错?ERROR: Low progression after class 2.
本帖最后由 lindlar 于 2021-5-30 23:17 编辑

试了好几次,都是出现下面的错误,不知道什么原因?
==================================================
              Calculations of Band Intensities
     ==================================================

-- To: vibronic fundamental state --
    Spectrum progression:    0.56%

-- To: single overtones --
    Spectrum progression:    3.48%

-- To: combinations of  2 simultaneously excited modes --
    Spectrum progression:   11.03%

ERROR: Low progression after class 2. Total convergence = 11.0%.
        The vibronic spectrum will likely be unreliable. Stopping.

计算方法参考:http://sobereva.com/223
第一步:# opt td=(root=1)  freq=savenm cam-b3lyp/6-311g(d,p) scrf=(solvent=n,n-dimethylacetamide)
第二步:# opt freq=(fc,emission) cam-b3lyp/6-311g(d,p) scrf=(solvent=n,n-dimethylacetamide)  
              。。。。
              _exc.chk(第一步中freq=saveNM任务的chk文件)


作者
Author:
sobereva    时间: 2021-5-31 05:10
论坛首页google框搜你的报错提示,诸如Low progression after class ,看以前的讨论。有问题先搜再问
作者
Author:
lindlar    时间: 2021-6-1 20:58
本帖最后由 lindlar 于 2021-6-1 21:00 编辑

不好意思,就是没有找到明确的解决方法,所以才问的。
作者
Author:
wzkchem5    时间: 2021-6-1 23:50
lindlar 发表于 2021-6-1 13:58
不好意思,就是没有找到明确的解决方法,所以才问的。

你把你搜到的帖子一一列举一下,并且告诉我们每个帖子分别缺哪方面的信息,这样我们才好针对性回答。
作者
Author:
kay    时间: 2021-12-30 10:36
楼主这个问题解决了吗?我也遇到此类问题,搜了帖子确实没有找到解决问题的办法,我的分子结构差异不算大。不知有没有好的解决方法?
作者
Author:
wzkchem5    时间: 2021-12-30 17:16
kay 发表于 2021-12-30 03:36
楼主这个问题解决了吗?我也遇到此类问题,搜了帖子确实没有找到解决问题的办法,我的分子结构差异不算大。 ...

还是得学会问问题。显然不会完全搜不到相关的帖子,只是搜到的帖子相关性不够而已,那你当然要告诉我们你已经搜到的帖子讲了什么、还缺哪些方面没有说清楚,这样我们才知道哪些方面需要给你澄清。
作者
Author:
Huangff    时间: 2024-3-26 09:02
lindlar 发表于 2021-6-1 20:58
不好意思,就是没有找到明确的解决方法,所以才问的。

您解决了吗

作者
Author:
lindlar    时间: 2024-3-26 20:03
没有解决
作者
Author:
wzkchem5    时间: 2024-3-26 22:50
本帖最后由 wzkchem5 于 2025-5-22 10:25 编辑

最近我这边有合作者也遇到了类似的问题。我个人虽然不太了解高斯计算振动分辨光谱的方法,但是大体来说,高斯默认是一个一个地把vibronic跃迁算出来的(https://doi.org/10.1063/1.2437197)。虽然如果提前已知一个峰的强度很弱,高斯可以不算这个峰,但是即便如此,需要算的峰的数目应该也是随体系大小指数增加的。所以对于足够大的体系,高斯会发现即使把按理说贡献最大的几类振动能级(没有任何振动模被激发的跃迁,只有一个振动模被激发的跃迁,有两个振动模被激发的跃迁)都算出来,也只算出了总振动分辨光谱的一小部分,也就是spectrum progression远小于100%(至于怎么知道算了总光谱的多大部分,可能是因为理论上可以求得整个振动分辨光谱的总吸收强度,然后把所有已经算的振动分辨峰的吸收强度加起来,和理论上的总吸收强度比较?)。此时高斯认为继续计算下去要么需要极长的时间才能算完,要么只能做极其粗糙的近似,导致结果不可靠,因此就放弃计算了。以上是我的个人理解,不一定对。
解决方法是改用不需要一个一个把vibronic跃迁算出来的方法,例如ORCA的ESD模块使用的路径积分方法,计算量应该是与体系大小的三次方成正比,而非指数增加。也可用MOMAP计算,需要用高斯、Q-Chem、BDF等程序提供基态和激发态的Hessian,原理类似。这样的缺点在于无法直接指认每个振动分辨峰是哪些振动导致的,也就是只能画出振动分辨光谱的包络线而没有底下表示vibronic跃迁的竖线;但有可能可以通过设置多个tcutfreq值、分别跑ESD的方式迂回解决(http://bbs.keinsci.com/forum.php ... 50634&fromuid=21811)。不过这种做法似乎没有文献先例,所以建议谨慎使用。
我们这边最近成功用ORCA把一个高斯算出来total convergence连1%都不到、因而报错的算例给算出来了,结果和实验符合不错,近期会投稿。这篇文章估计也会用上用tcutfreq指认振动分辨光谱的方法,如果审稿人对这种做法没有意见,这篇文章将来就可以作为文献先例,供其他想用ORCA算振动光谱、又想指认振动分辨峰的人来引用

EDIT(2025.5.22):这篇文章去年发出来了(https://doi.org/10.3390/molecules29102377),我们算的振动光谱在SI里,没有用到用tcutfreq指认振动分辨光谱的方法

作者
Author:
Huangff    时间: 2024-5-20 16:13
您解决了吗?

作者
Author:
wzkchem5    时间: 2024-7-13 02:07
Huangff 发表于 2024-5-20 09:13
您解决了吗?

你楼上我已经说了解决方法了




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