高斯优化激发态出错

小范范1989

高斯优化s1态出错，我也查过，说加上direct。但是我也加了，还是出错，不知道怎么改？谢谢

错误为

小范范1989

我优化t1也是出错，
输入文件：
%chk=test.chk
%mem=2000MB
%nprocs=2
#p b3lyp/6-31* opt td(direct,triplets,nstates=6,root=1)

test

0 1
N                  1.55801000    0.19783500   -0.79778700
C                  1.09085500   -0.95453600   -1.48926800
C                  2.14146600   -1.89758300   -1.54718200
C                 -0.13547400   -1.20956500   -2.10922100


错误
NMat=     8 NSing=     0.
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  1 NPass=   1 NMax=     8.
CISAX will form     8 AO SS matrices at one time.
NMat=     8 NSing=     0.
Root      1 has converged.
Root      2 has converged.
Root      3 has converged.
Root      4 has converged.
Root      5 not converged, maximum delta is    0.000011281594252
Root      6 not converged, maximum delta is    0.000006103042082
Excitation Energies [eV] at current iteration:
Root      1 :     1.916757579771638   Change is   -0.000000000000918
Root      2 :     2.060965374014521   Change is   -0.000000000000801
Root      3 :     2.576906613973796   Change is   -0.000000000107453
Root      4 :     2.603692086423475   Change is   -0.000000000049297
Root      5 :     2.632340173082165   Change is   -0.000000002306130
Root      6 :     2.637826619814619   Change is   -0.000000000652924
Iteration    13 Dimension   136 NMult   136
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  1 NPass=   1 NMax=     4.
CISAX will form     4 AO SS matrices at one time.
NMat=     4 NSing=     0.
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  1 NPass=   1 NMax=     4.
CISAX will form     4 AO SS matrices at one time.
NMat=     4 NSing=     0.

小范范1989

输入表头
%mem=2000MB
%nprocs=2
#p b3lyp/6-31* opt td(direct,triplets,nstates=6,root=1)

test

0 1
N                  1.50647400    1.80483400    0.04491200
C                  1.22289100    1.69315200    1.41808500

输出
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  1 NPass=   1 NMax=    12.
CISAX will form    12 AO SS matrices at one time.
NMat=    12 NSing=     0.
Root      1 not converged, maximum delta is    0.298200414274617
No map to state      2
No map to state      3
No map to state      4
New state      5 was old state      2
Root      5 not converged, maximum delta is    0.672121386957955
No map to state      6
You need to solve for more vectors in order to follow this state.
Error termination via Lnk1e in /home/software/g09/l914.exe at Sat Nov 29 21:23:51 2014.
Job cpu time:  3 days  1 hours  9 minutes 25.1 seconds.
File lengths (MBytes):  RWF=   2932 Int=      0 D2E=      0 Chk=    321 Scr=      1

小范范1989

#p b3lyp/6-31* opt

test

0 3
N                  1.90437600   -0.13917600   -0.93457000
C                  1.38502800   -1.42638600   -0.65539200

错误：
NTrRot=    -1 NTRed=   652 NAtoms=   108 NSkip=   334 IsLin=F
Error in internal coordinate system.
Error termination via Lnk1e in /home/software/g09/l103.exe at Fri Nov 28 21:28:37 2014.
Job cpu time:  0 days  5 hours  2 minutes 45.9 seconds.
File lengths (MBytes):  RWF=   1195 Int=      0 D2E=      0 Chk=     30 Scr=      1


对于错误，我已经看过帖子，说加上OPT=Cartesian，可以解决，但是，我还是出错。

sobereva

都是同一类问题，请发在同一个帖子里，不要分成四个单独发，否则看帖和回答都比较麻烦。我已经进行合并了

sobereva

1、3L：看此贴http://bbs.keinsci.com/forum.php?mod=viewthread&tid=299
2L：输出中并没有出错信息。如果还在算就继续等着，如果已经停了就重算，还不行就换版本或者改方法，或者直接给高斯中文客服发邮件
4L：改用opt=cartesian

小范范1989

sobereva 发表于 2014-11-30 11:52
都是同一类问题，请发在同一个帖子里，不要分成四个单独发，否则看帖和回答都比较麻烦。我已经进行合并了

嘿嘿，麻烦sob了。我以后会注意的。:lol

小范范1989

sobereva 发表于 2014-11-30 12:09
1、3L：看此贴http://bbs.keinsci.com/forum.php?mod=viewthread&tid=299
2L：输出中并没有出错信息。如果 ...

sob你好，我出现这样的错误，我按照你说的，设置nstates改大，但是还是出错，不管是单重态还是三重态。谢谢sob
输入%chk=test.chk
%mem=2GB
%nprocs=2
#p b3lyp/6-31* opt td(direct,triplets,nstates=20,root=1)

test

0 1
N                  1.50647400    1.80483400    0.04491200
C                  1.22289100    1.69315200    1.41808500
C                  2.22262000    0.88887400    2.02102800

输出
New state      1 was old state      3
Root      1 not converged, maximum delta is    0.571065717764176
New state      2 was old state     10
Root      2 not converged, maximum delta is    0.624581561156494
New state      3 was old state      1
Root      3 not converged, maximum delta is    0.020512414872001
New state      4 was old state      2
Root      4 not converged, maximum delta is    0.289690504679012
New state      5 was old state      4
Root      5 not converged, maximum delta is    0.436164850671245
New state      6 was old state      5
Root      6 not converged, maximum delta is    0.513899117437082
No map to state      7
No map to state      8
No map to state      9
No map to state     10
New state     11 was old state      6
Root     11 not converged, maximum delta is    0.455626045840689
No map to state     12
New state     13 was old state      8
Root     13 not converged, maximum delta is    0.309748276085905
No map to state     14
You need to solve for more vectors in order to follow this state.
Error termination via Lnk1e in /home/software/g09/l914.exe at Sun Dec  7 05:49:12 2014.
Job cpu time:  9 days 11 hours 37 minutes 46.1 seconds.
File lengths (MBytes):  RWF=   6110 Int=      0 D2E=      0 Chk=    844 Scr=      1

sobereva

又是三重态又是root=1，直接用0 3计算就行了，不需要TDDFT，又快又不会出现这些问题

小范范1989

sobereva 发表于 2014-12-7 15:43
又是三重态又是root=1，直接用0 3计算就行了，不需要TDDFT，又快又不会出现这些问题

嗯，好的，谢谢sob。我用两种方式都做了，确实0  3 出来的快多了，也没有错误。
但是，我优化s1的时候还是这样的错误。
我也是按照，nstates扩大，都到20了。
%chk=test.chk
%mem=2GB
%nprocs=2
#p b3lyp/6-31* opt=(MaxStep=10) td(nstates=20,singlets,root=1,direct)

test

0 1
N                  1.55801000    0.19783500   -0.79778700
C                  1.09085500   -0.95453600   -1.48926800

还是说you need tosolve.....

sobereva

nstates再设大一些。maxstep不仅设小，还应该加上notrust以免自动在优化中更新置信半径，否则步长仍可能超过maxstep的值。

小范范1989

sobereva 发表于 2014-12-7 16:35
nstates再设大一些。maxstep不仅设小，还应该加上notrust以免自动在优化中更新置信半径，否则步长仍可能超 ...

谢谢sob。thanks:lol

小范范1989

sobereva 发表于 2014-12-7 16:35
nstates再设大一些。maxstep不仅设小，还应该加上notrust以免自动在优化中更新置信半径，否则步长仍可能超 ...

sob，再问您个问题。还是激发态的问题。谢谢你
优化第二激发态后计算的结果发现第二激发态的能量比第一激发态能量还低，正常吗，下面分别是优化第一和第二激发态的结果。
                                
优化激发态（root=2）：
Excited State   1:      Singlet-A      2.4155 eV  513.29 nm  f=0.3072  <S**2>=0.000
     152 ->156         0.10882
     153 ->156         0.10052
     154 ->156         0.29460
     155 ->156        -0.62377

Excited State   2:      Singlet-A      1.6508 eV  751.06 nm  f=0.0014  <S**2>=0.000
     153 ->156         0.21356
     154 ->156         0.58364
     155 ->156         0.31673

Excited State   3:      Singlet-A      2.6231 eV  472.67 nm  f=0.1146  <S**2>=0.000
     148 ->156        -0.13364
     153 ->156         0.65309
     154 ->156        -0.22779
，这个出来的第二个激发能比第一个还小。这是怎么理解？优化的是第二激发态。
下面是我的理解
在基于s2的结构做td，出来的第一个td，也就是2.4155ev，指的是，从基态势能面，s2的结构处，直接td到s2势能面的最低处 的激发能。
而这个1.6508ev，指的是，基态势能面，s2的结构处，直接td，到s1势能面上的激发能（注意，不是s1时能面的的最低处，而是，垂直激发）。

sobereva

这是因为你选定对root=2进行优化，后来势能面发生交叉，能量高低顺序发生改变，目前的结构下曾经的root=2已经不是激发能第二高的态了，但是优化一直是对这个态来优化的。

flyskyzff

小范范1989 发表于 2014-12-7 15:56
嗯，好的，谢谢sob。我用两种方式都做了，确实0  3 出来的快多了，也没有错误。
但是，我优化s1的时候还 ...

请问优化t1,你这里说的“直接0 3计算就行了，不需要TDDFT”的意思是不是route部分只要# b3lyp/6-31* opt，然后电荷为0 ，自旋为3就可以了？