激发态频率计算出错，谢谢指点

小范范1989

我知道高斯计算激发态频率，不支持二阶解析导数，所以很慢（只是听sob老师经常说起，至于其中的具体原来是不知道的，嘿嘿）。但是我觉着我们的计算机的计算能力还行，所以想算算试试，同时，我记得sob老师给过建议，可以用CIS做激发态的freq。所以我做了一下两个输入文件，同时错误的信息我也给出。希望得到指点，谢谢
1：CIS的方法：
%mem=8000MB
%nprocs=8
#p CIS/3-21g freq nosymm

Title Card Required

0 1
C                  1.95828400    1.74831400   -3.42363200
C                  0.91004100    0.99931100   -2.91657800
C                  1.18654200   -0.13140200   -2.06627900
C                  2.54759000   -0.55414800   -1.88064100

....

错误信息为：
(Enter /home/software/g09/l801.exe)
ExpMin= 1.83D-01 ExpMax= 2.43D+02 ExpMxC= 2.43D+02 IAcc=1 IRadAn=         1 AccDes= 0.00D+00
HarFok:  IExCor=  205 AccDes= 0.00D+00 IRadAn=         1 IDoV=-2
ScaDFX=  1.000000  1.000000  1.000000  1.000000
Largest valence mixing into a core orbital is  3.11D-04
Largest core mixing into a valence orbital is  1.44D-04
Range of M.O.s used for correlation:    52   507
NBasis=   507 NAE=   168 NBE=   168 NFC=    51 NFV=     0
NROrb=    456 NOA=   117 NOB=   117 NVA=   339 NVB=   339
Leave Link  801 at Thu Sep 24 16:30:30 2015, MaxMem= 1048576000 cpu:       0.8
(Enter /home/software/g09/l804.exe)
Closed-shell transformation, MDV=  1048576000 ITran=3 ISComp=1.
Semi-Direct transformation.
ModeAB=           4 MOrb=           117 LenV=    1044118907
LASXX=   2254111704 LTotXX=  2254111704 LenRXX=  4546530261
LTotAB=  2292418557 MaxLAS=  3803464080 LenRXY=           0
NonZer=  6800641965 LenScr= 13654624768 LnRSAI=  3803464080
LnScr1=  7636761088 LExtra=           0 Total=  29641380197
MaxDsk=          -1 SrtSym=           T ITran=            3
DoSDTr:  NPSUse=  8
JobTyp=0 Pass  1:  I=   1 to 117.
Complete sort for first half transformation.
First half transformation complete.
Begin second half transformation for I=  10.
Begin second half transformation for I=  20.
Begin second half transformation for I=  30.
Begin second half transformation for I=  40.
Begin second half transformation for I=  50.
Begin second half transformation for I=  60.
Begin second half transformation for I=  70.
Begin second half transformation for I=  80.
Begin second half transformation for I=  90.
Begin second half transformation for I= 100.
Begin second half transformation for I= 110.
Erroneous write. Write 488685568 instead of 2097152000.
fd = 4
orig len = 8344566744 left = 4150262744
g_write




2：m062x的方法：
%mem=8000MB
%nprocs=8
#p m062x/3-21g freq nosymm td=(singlets,root=1)

Title Card Required

0 1
C                  1.95828400    1.74831400   -3.42363200
C                  0.91004100    0.99931100   -2.91657800
C                  1.18654200   -0.13140200   -2.06627900

...

我记得出错和上面的一样。

谢谢指点。

ps：顺带一个问题，刚刚也是我看计算文件不多，我搜索了一下，没有准确的方式：
我是优化t1  采用的是03的方式，我以前03的方式，都是很快，很方便的。
输入：
%nprocs=8
#p MPWB95/6-31g(d) IOP(3/76=0690003100) opt  scrf(solvent=toluene)

test

0 3
N                 -2.10875100   -2.08639200    1.04900800
C                 -1.86365800   -2.86953800    2.21201400
C                 -3.07740400   -3.02453200    2.91859200
C                 -0.66795200   -3.40879600    2.68460300
C                 -3.09517700   -3.74693100    4.11663500

错误：
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UmPW+HF-B95) =  -2806.27905584     A.U. after  129 cycles
             Convg  =    0.1713D-04             -V/T =  2.0090
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 1.0000 <S**2>= 2.0193 S= 1.0064
<L.S>= 0.000000000000E+00
KE= 2.781298477498D+03 PE=-2.753263034128D+04 EE= 1.143336098716D+04
Annihilation of the first spin contaminant:
S**2 before annihilation     2.0193,   after     2.0002
DiagDN has N=  1128 LTot=  3418 but NE2=     3 cannot use DSYEVD.
Convergence failure -- run terminated.
Error termination via Lnk1e in /home/software/g09/l502.exe at Thu Sep 24 02:27:07 2015.
Job cpu time:  1 days 10 hours 16 minutes 20.2 seconds.
File lengths (MBytes):  RWF=   2095 Int=      0 D2E=      0 Chk=     30 Scr=      1

sobereva

g09中CIS支持解析二阶导数，所以做依赖于二阶导数的振动分析任务快；而TDDFT只支持一阶解析导数，要通过有限差分获得二阶导数，做振动分析耗时是CIS的6N倍的程度（N=原子数）。

1 硬盘不够。确保硬盘剩余空间足够大，或者写上maxdisk来指明还剩多少硬盘空间。或者把CIS改成CIS(direct)再试。

2 视具体输出而定。没有准确信息不好说。

附加问题：老生常谈的不收敛
解决SCF不收敛问题的方法
http://sobereva.com/61

小范范1989

sobereva 发表于 2015-9-24 21:50
g09中CIS支持解析二阶导数，所以做依赖于二阶导数的振动分析任务快；而TDDFT只支持一阶解析导数，要通过有 ...

谢谢sob老师。嘿嘿

石头

你好，请问你是优化激发态时用的 opt td，但是算激发态频率使用CIS吗？这样做可以吗？我用高斯计算激发态频率算很久都没算出来。

sobereva

石头 发表于 2016-7-20 11:50
你好，请问你是优化激发态时用的 opt td，但是算激发态频率使用CIS吗？这样做可以吗？我用高斯计算激发态频 ...

这么做显然不行，频率计算用的级别必须和几何优化严格一致

石头

sobereva 发表于 2016-7-20 13:21
这么做显然不行，频率计算用的级别必须和几何优化严格一致

明白了，谢谢Sob老师，那在用TD计算激发态频率的时候有什么好办法吗？因为分子较大，很长时间都算不出来，很苦恼啊。

sobereva

石头 发表于 2016-7-20 15:21
明白了，谢谢Sob老师，那在用TD计算激发态频率的时候有什么好办法吗？因为分子较大，很长时间都算不出来 ...

没辙，要省时只能降低点基组。
或者等等下个Gaussian大版本，估计不久就会出来，应该会支持TDDFT二阶解析导数。或者买支持TDDFT解析Hessian的Q-Chem，但不是很划算。

石头

sobereva 发表于 2016-7-21 00:19
没辙，要省时只能降低点基组。
或者等等下个Gaussian大版本，估计不久就会出来，应该会支持TDDFT二阶 ...

好的，明白了，非常感谢Sob老师！