计算化学公社

标题: 新手求助L502计算出错 [打印本页]

作者
Author: lxl123456 时间: 2021-9-23 21:57
标题: 新手求助L502计算出错
%mem=1GB
# opt td=(nstates=5) cam-b3lyp/6-31g(d) IOp(9/40=4) scrf=(solvent=chloroform)

H

0 1
C                -2.50221000 1.63121600 -0.80871900
C                -1.68182600 0.65161400 -0.30024700
C                -0.57694200 1.21919900 0.45464800
C                -0.62875400 2.58935100 0.49752000
S                -1.98238800 3.22298100 -0.37874100
C                2.59030000 -0.37694600 1.70231400
C                1.89473800 0.47848300 0.87852000
C                0.46280600 0.41346600 1.13883400
C                0.15319400 -0.50385100 2.11005200
S                1.55997300 -1.27880700 2.75651900
C                2.48615700 1.39003600 -0.11077600
C                -1.80417100 -0.79026300 -0.56233600
C                3.65473100 1.32921600 -0.78503500
C                -2.89696100 -1.56386900 -0.73482500
C                -4.33411500 -1.26728200 -0.58015500
C                4.68206100 0.27139500 -0.86619500
C                -4.84089700 -0.34870700 0.35872600
C                -6.21305500 -0.14848100 0.49395600
C                -7.11524900 -0.86004100 -0.30295900
C                -6.63038100 -1.78660100 -1.22801400
C                -5.25704500 -1.99374700 -1.35564200
C                6.03246300 0.64670700 -0.98967600
C                7.03942000 -0.31046100 -1.11389400
C                6.71512700 -1.66829900 -1.14082200
C                5.37531800 -2.05696400 -1.04721400
C                4.37089500 -1.10028100 -0.91254100
H                -3.37696600 1.49999100 -1.42981100
H                0.03881700 3.26031000 1.02119500
H                3.65930200 -0.52773000 1.74180100
H                -0.82825900 -0.75176000 2.49055400
H                1.88300600 2.27267000 -0.30814500
H                -0.84532700 -1.29671800 -0.66105700
H                3.89091800 2.20972000 -1.38290000
H                -2.69648200 -2.59238300 -1.03494700
H                -4.15375900 0.20151600 0.99258400
H                -6.58091800 0.56123400 1.23046400
H                -8.18468400 -0.70044000 -0.19564000
H                -7.32098500 -2.35448200 -1.84590000
H                -4.88709700 -2.72565700 -2.07010100
H                6.29112500 1.70308800 -0.98175700
H                8.07599600 0.00541100 -1.19723200
H                7.49622100 -2.41649600 -1.24450300
H                5.11052400 -3.11038200 -1.08741700
H                3.33360100 -1.41401300 -0.85477300
我的输入文件
Error on total polarization charges =  0.01209
SCF Done:  E(RCAM-B3LYP) =  -1721.01565900    A.U. after  129 cycles
         NFock=128  Conv=0.16D-04    -V/T= 2.0062
Convergence failure -- run terminated.
Error termination via Lnk1e in D:\app\GAUSSIAN\G09W\l502.exe at Wed Sep 22 19:23:57 2021.
Job cpu time:    2 days 21 hours  8 minutes 39.0 seconds.
File lengths (MBytes):  RWF= 749 Int=    0 D2E=    0 Chk=    52 Scr=    1
输出文件的最后
其他两个基团都能计算完，就这个基团

作者
Author: wzkchem5 时间: 2021-9-23 22:09
结构优化的最后一帧的激发能是不是很接近0？如果是的话，说明结构优化到了S0/S1锥形交叉点，也就是S1平衡结构不存在

作者
Author: lxl123456 时间: 2021-9-23 22:59

wzkchem5 发表于 2021-9-23 22:09
结构优化的最后一帧的激发能是不是很接近0？如果是的话，说明结构优化到了S0/S1锥形交叉点，也就是S1平衡结 ...

Rare condition: small coef for last iteration: 0.000D+00
Rare condition: small coef for last iteration: 0.000D+00
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
是这个吗？纯新手

作者
Author: wzkchem5 时间: 2021-9-23 23:51

lxl123456 发表于 2021-9-23 15:59
Rare condition: small coef for last iteration: 0.000D+00
Rare condition: small coef for last it ...

不是，我问的是TDDFT激发能
如果不会看TDDFT激发能，先用这个论坛上面的搜索框，搜怎么看

作者
Author: sobereva 时间: 2021-9-24 08:30
老生常谈的问题，有问题先搜公社首页google框，找不到答案再发帖问
解决SCF不收敛问题的方法
http://sobereva.com/61

作者
Author: lxl123456 时间: 2021-9-25 07:55

sobereva 发表于 2021-9-24 08:30
老生常谈的问题，有问题先搜公社首页google框，找不到答案再发帖问
解决SCF不收敛问题的方法
http://sobe ...

老师，是用第六条的方法解决吗？

作者
Author: sobereva 时间: 2021-9-25 08:20

lxl123456 发表于 2021-9-25 07:55
老师，是用第六条的方法解决吗？

把文章开头的字仔细看了

作者
Author: lxl123456 时间: 2021-9-25 10:28

wzkchem5 发表于 2021-9-23 23:51
不是，我问的是TDDFT激发能
如果不会看TDDFT激发能，先用这个论坛上面的搜索框，搜怎么看

Excited State 5:    Singlet-A    3.3590 eV  369.11 nm  f=0.0156  <S**2>=0.000
   39 ->134       -0.00014
   39 ->135       -0.00015
最后一个不是0

作者
Author: wzkchem5 时间: 2021-9-25 16:21

lxl123456 发表于 2021-9-25 03:28
Excited State 5: Singlet-A 3.3590 eV 369.11 nm f=0.0156 =0.000
39 ->134 ...

应该看第一个激发能。看激发能的目的是为了看是否到了S0/S1锥形交叉点，当且仅当S1的激发能接近0时，可以认为接近S0/S1锥形交叉点。和S5激发能没有任何关系

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