计算化学公社

标题: 高斯TD运行出错 [打印本页]

作者
Author: xamiracle407 时间: 2015-9-22 19:09
标题: 高斯TD运行出错
小弟，用cam-b3lyp/3-21g 来对激发态进行优化，遇到这样的运行错误，不知道是哪个地方出错了，还望大神指点迷津，谢谢~~
SCF Done:  E(RCAM-B3LYP) =  -2359.09050787    A.U. after 19 cycles
         Convg  = 0.3457D-08          -V/T =  2.0062
KE= 2.344485029694D+03 PE=-1.144423071194D+04 EE= 3.801274037537D+03
Leave Link  502 at Tue Sep 22 16:02:51 2015, MaxMem=16106127360 cpu: 1472.8
(Enter /home/zy/g09/l801.exe)
ExpMin= 2.44D-02 ExpMax= 1.21D+03 ExpMxC= 1.21D+03 IAcc=3 IRadAn=       5 AccDes= 0.00D+00
HarFok:  IExCor=  205 AccDes= 0.00D+00 IRadAn=       5 IDoV=-2
ScaDFX=  1.000000  1.000000  1.000000  1.000000
Largest valence mixing into a core orbital is  4.07D-01
Largest core mixing into a valence orbital is  3.93D-01
Excessive mixing of frozen core and valence orbitals.
Error termination via Lnk1e in /home/zy/g09/l801.exe at Tue Sep 22 16:02:52 2015.
Job cpu time:  1 days 10 hours  5 minutes  3.8 seconds.
File lengths (MBytes):  RWF= 682 Int=    0 D2E=    0 Chk=    31 Scr=

作者
Author: sobereva 时间: 2015-9-22 19:24
这说明内核轨道和价层轨道之间有明显混合，冻核近似此时可能不合理。最简单的做法就是不用冻核近似，即写成TD(full)，计算量会相应地增加。

作者
Author: xamiracle407 时间: 2015-9-22 19:28
谢谢老师，我在试一下

作者
Author: xamiracle407 时间: 2015-9-22 19:28

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

谢谢老师，我再试一下

作者
Author: xamiracle407 时间: 2015-9-26 08:49

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

作者
Author: sobereva 时间: 2015-9-26 18:27
老生常谈的SCF不收敛
解决SCF不收敛问题的方法
http://sobereva.com/61

建议先尝试换泛函，CAM-B3LYP对这样的体系并没什么好处。

作者
Author: xamiracle407 时间: 2015-9-26 18:37

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 ,这样子可行吗

作者
Author: xamiracle407 时间: 2015-9-26 18:40

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 ,这样子可行吗，老师？

作者
Author: sobereva 时间: 2015-9-26 21:00

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分成两步做。
其它的没问题。

作者
Author: xamiracle407 时间: 2015-9-26 21:42

sobereva 发表于 2015-9-26 21:00
对你的体系nosymm没意义。scf=diis本来就是默认的没必要写。
做TDDFT激发态振动分析非常耗时，确定是 ...

好的，感谢sob老师的宝贵建议，我再做一下

作者
Author: xamiracle407 时间: 2015-9-28 09:00

sobereva 发表于 2015-9-26 21:00
对你的体系nosymm没意义。scf=diis本来就是默认的没必要写。
做TDDFT激发态振动分析非常耗时，确定是 ...

sob老师，激发态优化不如基态优化容易，又出了问题了，（图片，输入输出文件附上）关键词：#p opt td=full b3lyp/6-31g(d) scrf=(solvent=water) nosymm geom=connectivity maxdisk=200GB
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RB3LYP) =  -2372.80347919    A.U. after  129 cycles
         Convg  = 0.4710D-03          -V/T =  2.0069
KE= 2.356564680829D+03 PE=-1.161040295842D+04 EE= 3.874222937388D+03
Convergence failure -- run terminated.
Error termination via Lnk1e in /home/zy/g09/l502.exe at Sun Sep 27 20:08:46 2015.
Job cpu time:  9 days 21 hours 46 minutes 37.8 seconds.
File lengths (MBytes):  RWF= 1443 Int=    0 D2E=    0 Chk=    50 Scr=    1

作者
Author: sobereva 时间: 2015-9-28 21:06

xamiracle407 发表于 2015-9-28 09:00
sob老师，激发态优化不如基态优化容易，又出了问题了，（图片，输入输出文件附上）关键词：#p opt td=fu ...

这都是数值巧合
还是按博文里的方法尝试加上各种帮助SCF收敛的关键词，比如scf(noincfock,novaracc) int=ultrafine或者scf=xqc或者vshift之类的再试，或者尝试其它泛函或基组。

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