计算化学公社

标题: 态求和计算150states还是没有呈现收敛趋势，需要继续增加态数嘛？ [打印本页]

作者
Author: cornerj 时间: 2019-4-14 19:52
标题: 态求和计算150states还是没有呈现收敛趋势，需要继续增加态数嘛？
(, 下载次数 Times of downloads: 54)
请问态求和算β0值，计算150states，用Multiwfn导出来来的β0与态数的关系图，还是没有明显的收敛趋势，还有必要继续增加态数吗？

使用计算水平CAM-B3LYP/aug-cc-pVDZ，体系原子数60多个，再增加态数计算耗时太大估计算不动了。请问该如何处理这种情况？

作者
Author: sobereva 时间: 2019-4-15 09:43
只能再增加，没办法。这是SOS方法的缺点。估计到了两百多个态总应该能收敛了

作者
Author: cornerj 时间: 2019-4-15 21:52

sobereva 发表于 2019-4-15 09:43
只能再增加，没办法。这是SOS方法的缺点。估计到了两百多个态总应该能收敛了

哦哦那我老老实实算个200试试看。谢谢老师！
还有个小问题，从图上看出来第3,4个激发态有很明显的正、负贡献，
怎么用Multiwfn导出这几个激发态对应的轨道跃迁组成？我在TD的输出文件里找不到轨道跃迁信息、

作者
Author: sobereva 时间: 2019-4-16 00:08

cornerj 发表于 2019-4-15 21:52
哦哦那我老老实实算个200试试看。谢谢老师！
还有个小问题，从图上看出来第3,4个激发态有很明显的正、负 ...

轨道跃迁信息肯定有，不然Multiwfn都没法产生用于做SOS的文件
如果你嫌输出量太大不好找主要的，可以用这个程序处理
简化用了IOp(9/40=x)的Gaussian的CIS/TDDFT任务的输出文件的程序：simpIOp940
http://sobereva.com/467（http://bbs.keinsci.com/thread-12188-1-1.html）

在Multiwfn里也可以直接输出某个电子激发中主要的轨道跃迁，见
电子激发任务中轨道跃迁贡献的计算
http://sobereva.com/230

作者
Author: cornerj 时间: 2019-4-16 12:52

sobereva 发表于 2019-4-16 00:08
轨道跃迁信息肯定有，不然Multiwfn都没法产生用于做SOS的文件
如果你嫌输出量太大不好找主要的，可以用 ...

尝试了下，输出特别迅速！谢谢老师！

Some MO transitions sorted by absolute contributions:
#  43327 197B ->  208B Coeff.: 0.53203 Contri.: 28.3056%
#  23226 210A ->  213A Coeff.:  -0.37279 Contri.: 13.8972%

比如开壳层体系，对应的轨道跃迁197B ->  208B  210A ->  213A
我想知道这两个跃迁涉及到的分子轨道 197B 208B 210A 213A 各自的分子轨道图是哪几个？
因为用Multiwfn导出来的分子轨道的序号：206 is HOMO of alpha spin, orbital  1045 is HOMO of beta spin
请问这两种序号是存在对应关系的嘛？谢谢！

作者
Author: sobereva 时间: 2019-4-17 00:45

cornerj 发表于 2019-4-16 12:52
尝试了下，输出特别迅速！谢谢老师！

Some MO transitions sorted by absolute contributions:

使用Multiwfn观看分子轨道
http://sobereva.com/269（http://bbs.keinsci.com/thread-462-1-1.html）

这里说了，主功能0里面在文本框里输入负号就代表选择beta轨道，正值就代表alpha

作者
Author: cornerj 时间: 2019-4-17 07:43

sobereva 发表于 2019-4-17 00:45
使用Multiwfn观看分子轨道
http://sobereva.com/269（http://bbs.keinsci.com/thread-462-1-1.html）

...

明白这个操作了，感谢老师深夜解答！
有个莫名其妙的问题，优化和SOS一低一高俩计算级别都产生了fchk文件，
想看SOS中激发态对应的轨道跃迁，用哪一个fchk？这俩fchk包含的轨道信息有没有区别？

作者
Author: sobereva 时间: 2019-4-17 13:10

cornerj 发表于 2019-4-17 07:43
明白这个操作了，感谢老师深夜解答！
有个莫名其妙的问题，优化和SOS一低一高俩计算级别都产生了fchk文 ...

显然用你做电子激发计算任务产生的fch文件

作者
Author: cornerj 时间: 2019-4-17 13:17

sobereva 发表于 2019-4-17 13:10
显然用你做电子激发计算任务产生的fch文件

了解了，感谢老师的解答！

作者
Author: Yuanyaya 时间: 2021-1-13 17:40
本帖最后由 Yuanyaya 于 2021-1-13 17:47 编辑

sob老师您好，我看了您使用Multiwfn基于完全态求和(SOS)方法计算极化率和超极化率（http://sobereva.com/232）的帖子，用到自己体系分析γ的时候，有时在100个态左右就收敛的很好，有时会出现算200个态仍不收敛的情况（如图），请问对于不收敛的可以怎么修改呢？
计算td用的cam-b3lyp/genecp em=gd3 td(nstates=200) IOp(9/40=5) nosymm，C H 用6-31+g(d,p)，Sn 用 Lanl2DZ，大概90个原子。C:\Users\Yuanyaya\Desktop

作者
Author: biogon 时间: 2021-1-13 19:39

Yuanyaya 发表于 2021-1-13 17:40
sob老师您好，我看了您使用Multiwfn基于完全态求和(SOS)方法计算极化率和超极化率（http://sobereva.com/23 ...

不收敛就加多态数目继续算

作者
Author: sobereva 时间: 2021-1-16 08:02

Yuanyaya 发表于 2021-1-13 17:40
sob老师您好，我看了您使用Multiwfn基于完全态求和(SOS)方法计算极化率和超极化率（http://sobereva.com/23 ...

6-31+g(d,p)算（超）极化率烂得一塌糊涂。用aug-cc-pVDZ还勉勉强强
多少态能收敛和体系关系很大，算到几百个态才收敛都有可能。

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