xamiracle407 发表于 2015-9-28 09:00 这都是数值巧合 还是按博文里的方法尝试加上各种帮助SCF收敛的关键词,比如scf(noincfock,novaracc) int=ultrafine或者scf=xqc或者vshift之类的再试,或者尝试其它泛函或基组。 |
sobereva 发表于 2015-9-26 21:00 好的,感谢sob老师的宝贵建议,我再做一下 ![]() |
xamiracle407 发表于 2015-9-26 18:40 对你的体系nosymm没意义。scf=diis本来就是默认的没必要写。 做TDDFT激发态振动分析非常耗时,确定是否真的有必要做。建议把opt和freq分成两步做。 其它的没问题。 |
sobereva 发表于 2015-9-26 18:27 笔误,是关键词#p opt freq td=full b3lyp/6-31g(d) scrf=(solvent=water) nosymm scf=diis ,这样子可行吗,老师? |
sobereva 发表于 2015-9-26 18:27 SOB老师,输入文件关键词改成如下:#p opt freq td=full b3lyp/6-11g(d) scrf=(solvent=water) nosymm scf=diis ,这样子可行吗 |
老生常谈的SCF不收敛 解决SCF不收敛问题的方法 http://sobereva.com/61 建议先尝试换泛函,CAM-B3LYP对这样的体系并没什么好处。 |
sobereva 发表于 2015-9-22 19:24 Sob老师,按照您说的,最后还是出问题了,输入输出文件已上传,麻烦老师,再指点迷津~~Error on total polarization charges = 0.01271 E= -2372.19856922817 Delta-E= -0.000206195058 Rises=F Damp=F DIIS: error= 6.59D-04 at cycle 128 NSaved= 20. NSaved=20 IEnMin= 8 EnMin= -2372.20440193226 IErMin= 4 ErrMin= 2.41D-04 ErrMax= 6.59D-04 EMaxC= 1.00D-01 BMatC= 5.09D-04 BMatP= 2.36D-05 IDIUse=3 WtCom= 9.98D-01 WtEn= 2.41D-03 Matrix for removal 19 Erem= -2372.19836303312 Crem= 0.000D+00 Coeff-Com: 0.886D-02-0.379D-01-0.181D+00-0.385D+00 0.767D+00 0.799D+01 Coeff-Com: -0.434D+01-0.513D+01 0.283D+01-0.346D+01 0.114D+01 0.166D+01 Coeff-Com: -0.473D+01 0.434D+01 0.285D-01 0.135D+01-0.642D+00-0.146D+01 Coeff-Com: 0.794D+00 0.443D+00 Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+00 Coeff-En: 0.000D+00 0.100D+01 Coeff: 0.884D-02-0.378D-01-0.181D+00-0.384D+00 0.765D+00 0.797D+01 Coeff: -0.433D+01-0.512D+01 0.283D+01-0.345D+01 0.114D+01 0.165D+01 Coeff: -0.472D+01 0.433D+01 0.284D-01 0.135D+01-0.640D+00-0.145D+01 Coeff: 0.793D+00 0.444D+00 Gap= 0.157 Goal= None Shift= 0.000 RMSDP=8.44D-06 MaxDP=4.42D-04 DE=-2.06D-04 OVMax= 2.71D-03 >>>>>>>>>> Convergence criterion not met. SCF Done: E(RCAM-B3LYP) = -2372.19856923 A.U. after 129 cycles Convg = 0.8438D-05 -V/T = 2.0067 KE= 2.356481944463D+03 PE=-1.167371682379D+04 EE= 3.906361056201D+03 Convergence failure -- run terminated. Error termination via Lnk1e in /home/zy/g09/l502.exe at Fri Sep 25 21:06:20 2015. Job cpu time: 70 days 15 hours 28 minutes 34.0 seconds. File lengths (MBytes): RWF= 1482 Int= 0 D2E= 0 Chk= 86 Scr= 1 |
2.89 KB, 下载次数 Times of downloads: 18
5.25 MB, 下载次数 Times of downloads: 4
sobereva 发表于 2015-9-22 19:24 谢谢老师,我再试一下 ![]() |
谢谢老师,我在试一下![]() |
这说明内核轨道和价层轨道之间有明显混合,冻核近似此时可能不合理。最简单的做法就是不用冻核近似,即写成TD(full),计算量会相应地增加。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-15 15:37 , Processed in 0.254893 second(s), 25 queries , Gzip On.