Yuanyaya 发表于 2021-1-13 17:40 6-31+g(d,p)算(超)极化率烂得一塌糊涂。用aug-cc-pVDZ还勉勉强强 多少态能收敛和体系关系很大,算到几百个态才收敛都有可能。 |
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 |
Graph1.jpg (110.76 KB, 下载次数 Times of downloads: 32)
sobereva 发表于 2019-4-17 13:10 了解了,感谢老师的解答! |
cornerj 发表于 2019-4-17 07:43 显然用你做电子激发计算任务产生的fch文件 |
sobereva 发表于 2019-4-17 00:45 明白这个操作了,感谢老师深夜解答! 有个莫名其妙的问题,优化和SOS一低一高俩计算级别 都产生了fchk文件, 想看SOS中激发态对应的轨道跃迁,用哪一个fchk?这俩fchk包含的轨道信息有没有区别? |
cornerj 发表于 2019-4-16 12:52 使用Multiwfn观看分子轨道 http://sobereva.com/269(http://bbs.keinsci.com/thread-462-1-1.html) 这里说了,主功能0里面在文本框里输入负号就代表选择beta轨道,正值就代表alpha |
sobereva 发表于 2019-4-16 00:08 尝试了下,输出特别迅速!谢谢老师! 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 请问这两种序号是存在对应关系的嘛?谢谢! |
cornerj 发表于 2019-4-15 21:52 轨道跃迁信息肯定有,不然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 |
sobereva 发表于 2019-4-15 09:43 哦哦那我老老实实算个200试试看。谢谢老师! 还有个小问题,从图上看出来第3,4个激发态有很明显的正、负贡献, 怎么用Multiwfn导出这几个激发态对应的轨道跃迁组成?我在TD的输出文件里找不到轨道跃迁信息、 |
| 只能再增加,没办法。这是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.