计算化学公社

 找回密码 Forget password
 注册 Register

高斯TD运行出错

查看数: 16429 | 评论数: 11 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2015-9-22 19:09

正文摘要:

小弟,用cam-b3lyp/3-21g 来对激发态进行优化,遇到这样的运行错误,不知道是哪个地方出错了,还望大神指点迷津,谢谢~~ SCF Done:  E(RCAM-B3LYP) =  -2359.09050787      ...

回复 Reply

sobereva 发表于 Post on 2015-9-28 21:06:04
xamiracle407 发表于 2015-9-28 09:00
sob老师,激发态优化不如基态优化容易,又出了问题了,(图片,输入输出文件附上) 关键词:#p opt td=fu ...


这都是数值巧合
还是按博文里的方法尝试加上各种帮助SCF收敛的关键词,比如scf(noincfock,novaracc) int=ultrafine或者scf=xqc或者vshift之类的再试,或者尝试其它泛函或基组。
xamiracle407 发表于 Post on 2015-9-26 21:42:30
sobereva 发表于 2015-9-26 21:00
对你的体系nosymm没意义。scf=diis本来就是默认的没必要写。
做TDDFT激发态振动分析非常耗时,确定是 ...

好的,感谢sob老师的宝贵建议,我再做一下
sobereva 发表于 Post on 2015-9-26 21:00:59
xamiracle407 发表于 2015-9-26 18:40
笔误,是关键词#p opt freq td=full b3lyp/6-31g(d) scrf=(solvent=water) nosymm scf=diis ,这样子可行 ...


对你的体系nosymm没意义。scf=diis本来就是默认的没必要写。
做TDDFT激发态振动分析非常耗时,确定是否真的有必要做。建议把opt和freq分成两步做。
其它的没问题。
xamiracle407 发表于 Post on 2015-9-26 18:40:57
sobereva 发表于 2015-9-26 18:27
老生常谈的SCF不收敛
解决SCF不收敛问题的方法
http://sobereva.com/61

笔误,是关键词#p opt freq td=full b3lyp/6-31g(d) scrf=(solvent=water) nosymm scf=diis ,这样子可行吗,老师?
xamiracle407 发表于 Post on 2015-9-26 18:37:10
sobereva 发表于 2015-9-26 18:27
老生常谈的SCF不收敛
解决SCF不收敛问题的方法
http://sobereva.com/61

SOB老师,输入文件关键词改成如下:#p opt freq td=full b3lyp/6-11g(d) scrf=(solvent=water) nosymm scf=diis ,这样子可行吗
sobereva 发表于 Post on 2015-9-26 18:27:31
老生常谈的SCF不收敛
解决SCF不收敛问题的方法
http://sobereva.com/61

建议先尝试换泛函,CAM-B3LYP对这样的体系并没什么好处。
xamiracle407 发表于 Post on 2015-9-26 08:49:50
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

RLH_21_26.gjf~

2.89 KB, 下载次数 Times of downloads: 18

RLH_21_26.LOG

5.25 MB, 下载次数 Times of downloads: 4

xamiracle407 发表于 Post on 2015-9-22 19:28:56
sobereva 发表于 2015-9-22 19:24
这说明内核轨道和价层轨道之间有明显混合,冻核近似此时可能不合理。最简单的做法就是不用冻核近似,即写成 ...

谢谢老师,我再试一下
xamiracle407 发表于 Post on 2015-9-22 19:28:29
谢谢老师,我在试一下
sobereva 发表于 Post on 2015-9-22 19:24:53
这说明内核轨道和价层轨道之间有明显混合,冻核近似此时可能不合理。最简单的做法就是不用冻核近似,即写成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.

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