计算化学公社

 找回密码 Forget password
 注册 Register

态求和计算150states还是没有呈现收敛趋势,需要继续增加态数嘛?

查看数: 7759 | 评论数: 11 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2019-4-14 19:52

正文摘要:

请问态求和算β0值,计算150states,用Multiwfn导出来来的β0与态数的关系图,还是没有明显的收敛趋势,还有必要继续增加态数吗? 使用计算水平CAM-B3LYP/aug-cc-pVDZ,体系原子数60多个,再增加态数计算耗时 ...

回复 Reply

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

6-31+g(d,p)算(超)极化率烂得一塌糊涂。用aug-cc-pVDZ还勉勉强强
多少态能收敛和体系关系很大,算到几百个态才收敛都有可能。
biogon 发表于 Post on 2021-1-13 19:39:51
Yuanyaya 发表于 2021-1-13 17:40
sob老师您好,我看了您使用Multiwfn基于完全态求和(SOS)方法计算极化率和超极化率(http://sobereva.com/23 ...

不收敛就加多态数目继续算
Yuanyaya 发表于 Post on 2021-1-13 17:40:39
本帖最后由 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

Graph1.jpg (110.76 KB, 下载次数 Times of downloads: 32)

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

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

显然用你做电子激发计算任务产生的fch文件
cornerj 发表于 Post on 2019-4-17 07:43:27
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包含的轨道信息有没有区别?
sobereva 发表于 Post on 2019-4-17 00:45:47
cornerj 发表于 2019-4-16 12:52
尝试了下,输出特别迅速!谢谢老师!

Some MO transitions sorted by absolute contributions:

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

这里说了,主功能0里面在文本框里输入负号就代表选择beta轨道,正值就代表alpha
cornerj 发表于 Post on 2019-4-16 12:52:20
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
   请问这两种序号是存在对应关系的嘛?谢谢!


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

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

在Multiwfn里也可以直接输出某个电子激发中主要的轨道跃迁,见
电子激发任务中轨道跃迁贡献的计算
http://sobereva.com/230
cornerj 发表于 Post on 2019-4-15 21:52:15
sobereva 发表于 2019-4-15 09:43
只能再增加,没办法。这是SOS方法的缺点。估计到了两百多个态总应该能收敛了

哦哦那我老老实实算个200试试看。谢谢老师!
还有个小问题,从图上看出来第3,4个激发态有很明显的正、负贡献,
怎么用Multiwfn导出这几个激发态对应的轨道跃迁组成?我在TD的输出文件里找不到轨道跃迁信息、
sobereva 发表于 Post on 2019-4-15 09:43:02
只能再增加,没办法。这是SOS方法的缺点。估计到了两百多个态总应该能收敛了

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-24 08:30 , Processed in 0.197080 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list