计算化学公社

标题: 激发态频率计算出错，谢谢指点 [打印本页]

作者
Author: 小范范1989 时间: 2015-9-24 20:54
标题: 激发态频率计算出错，谢谢指点
本帖最后由小范范1989 于 2015-9-24 21:19 编辑

我知道高斯计算激发态频率，不支持二阶解析导数，所以很慢（只是听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

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

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

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

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

作者
Author: 小范范1989 时间: 2015-9-25 07:39

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

谢谢sob老师。嘿嘿

作者
Author: 石头 时间: 2016-7-20 11:50
你好，请问你是优化激发态时用的 opt td，但是算激发态频率使用CIS吗？这样做可以吗？我用高斯计算激发态频率算很久都没算出来。

作者
Author: sobereva 时间: 2016-7-20 13:21

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

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

作者
Author: 石头 时间: 2016-7-20 15:21

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

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

作者
Author: sobereva 时间: 2016-7-21 00:19

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

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

作者
Author: 石头 时间: 2016-7-21 08:35

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

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

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