高斯计算单点能L801 L804报错问题

桃子量子

使用CCSD(T)/cc-pvtz计算单点能报错，
如下所示，
PCMQM: There is more than one cavity, setting IComp=0.
Error on total polarization charges =  0.06888
SCF Done:  E(RHF) =  -1432.92974520     A.U. after   12 cycles
            NFock= 12  Conv=0.61D-08     -V/T= 2.0008
KE= 1.431754411466D+03 PE=-6.118278208422D+03 EE= 1.916389221833D+03
SMD-CDS (non-electrostatic) energy       (kcal/mol) =       6.22
(included in total energy above)
Leave Link  502 at Wed May 31 10:53:07 2023, MaxMem=   536870912 cpu:      1244.9
(Enter /data/home/liuzzgroup/g09/g09/l801.exe)
Windowed orbitals will be sorted by symmetry type.
ExpMin= 1.03D-01 ExpMax= 3.74D+05 ExpMxC= 1.18D+03 IAcc=3 IRadAn=         5 AccDes= 0.00D+00
HarFok:  IExCor=  205 AccDes= 0.00D+00 IRadAn=         5 IDoV=-2 UseB2=F ITyADJ=14
ICtDFT= 12500011 ScaDFX=  1.000000  1.000000  1.000000  1.000000
Largest valence mixing into a core orbital is  5.36D-04
Largest core mixing into a valence orbital is  8.79D-05
Range of M.O.s used for correlation:    19   438
NBasis=   438 NAE=    62 NBE=    62 NFC=    18 NFV=     0
NROrb=    420 NOA=    44 NOB=    44 NVA=   376 NVB=   376
Leave Link  801 at Wed May 31 10:53:08 2023, MaxMem=   536870912 cpu:        22.4
(Enter /data/home/liuzzgroup/g09/g09/l804.exe)
Closed-shell transformation, MDV=   536870912 ITran=4 ISComp=1.
Semi-Direct transformation.
ModeAB=           4 MOrb=            44 LenV=     532500122
LASXX=   1522847491 LTotXX=  1522847491 LenRXX=  3056649871
LTotAB=  1533802380 MaxLAS=  1980039600 LenRXY=           0
NonZer=  4579497362 LenScr=  9194915328 LnRSAI=  1980039600
LnScr1=  2981709824 LExtra=           0 Total=  17213314623
MaxDsk=          -1 SrtSym=           T ITran=            4
DoSDTr:  NPSUse= 28
JobTyp=0 Pass  1:  I=   1 to  44.
Complete sort for first half transformation.
First half transformation complete.
Begin second half transformation for I=  10.
Erroneous write. Write -1 instead of 1355776.
fd = 4
orig len = 8384512 left = 1355776
g_write

输入文件如下
%nprocshared=28
%mem=4GB
%chk=C3F7SO3-2-C2-7-IRC-1-E.chk
#p ccsd(t)/cc-pvtz scrf=(smd,solvent=water) nosymm


C3F7SO3-2-C2-7-IRC-1-E


-2 1
C                 -2.05855100    0.31435000   -0.05854500
C                 -0.71337000   -0.38912500   -0.34797500
C                  0.50902900    0.51294500   -0.07535300
S                  2.09876700   -0.40643800    0.14705000
O                  3.08796800    0.65478800    0.05083600
O                  1.97729800   -1.00664500    1.46417600
O                  2.09499800   -1.34630800   -0.96347400
F                 -0.71324300   -0.74856800   -1.63776800
F                 -2.05509400    1.53651000   -0.58229500
F                 -0.65864600   -1.49070000    0.41005500
F                  0.64207400    1.35551400   -1.11212700
F                 -2.26740600    0.40562000    1.24578500
F                  0.27225900    1.24346200    1.02992000
F                 -3.04907800   -0.37610200   -0.60326300
H                 -3.33500200   -2.73923400    1.40927900
看了报错的解析帖子，加入了full，仍然报错

报错如下，
PCMQM: There is more than one cavity, setting IComp=0.
Error on total polarization charges =  0.06888
SCF Done:  E(RHF) =  -1432.92974520     A.U. after   12 cycles
            NFock= 12  Conv=0.61D-08     -V/T= 2.0008
KE= 1.431754411466D+03 PE=-6.118278208422D+03 EE= 1.916389221833D+03
SMD-CDS (non-electrostatic) energy       (kcal/mol) =       6.22
(included in total energy above)
Leave Link  502 at Wed May 31 11:43:20 2023, MaxMem=   536870912 cpu:      1274.5
(Enter /data/home/liuzzgroup/g09/g09/l801.exe)
Windowed orbitals will be sorted by symmetry type.
Range of M.O.s used for correlation:     1   438
NBasis=   438 NAE=    62 NBE=    62 NFC=     0 NFV=     0
NROrb=    438 NOA=    62 NOB=    62 NVA=   376 NVB=   376
Leave Link  801 at Wed May 31 11:43:21 2023, MaxMem=   536870912 cpu:         0.4
(Enter /data/home/liuzzgroup/g09/g09/l804.exe)
Closed-shell transformation, MDV=   536870912 ITran=4 ISComp=1.
Semi-Direct transformation.
ModeAB=           4 MOrb=            62 LenV=     532271288
LASXX=   2366132722 LTotXX=  2366132722 LenRXX=  4749606751
LTotAB=  2383474029 MaxLAS=  2909629620 LenRXY=           0
NonZer=  7115739473 LenScr= 14287292928 LnRSAI=  2909629620
LnScr1=  5842082816 LExtra=           0 Total=  27788612115
MaxDsk=          -1 SrtSym=           T ITran=            4
DoSDTr:  NPSUse= 28
JobTyp=0 Pass  1:  I=   1 to  62.
Complete sort for first half transformation.
Erroneous write. Write -1 instead of 1782749344.
fd = 4
orig len = 4248942120 left = 1782749344
g_write
输入文件如下
%nprocshared=28
%mem=4GB
%chk=C3F7SO3-2-C2-7-IRC-1-E.chk
#p ccsd(t,full)/cc-pvtz scrf=(smd,solvent=water) nosymm


C3F7SO3-2-C2-7-IRC-1-E


-2 1
C                 -2.05855100    0.31435000   -0.05854500
C                 -0.71337000   -0.38912500   -0.34797500
C                  0.50902900    0.51294500   -0.07535300
S                  2.09876700   -0.40643800    0.14705000
O                  3.08796800    0.65478800    0.05083600
O                  1.97729800   -1.00664500    1.46417600
O                  2.09499800   -1.34630800   -0.96347400
F                 -0.71324300   -0.74856800   -1.63776800
F                 -2.05509400    1.53651000   -0.58229500
F                 -0.65864600   -1.49070000    0.41005500
F                  0.64207400    1.35551400   -1.11212700
F                 -2.26740600    0.40562000    1.24578500
F                  0.27225900    1.24346200    1.02992000
F                 -3.04907800   -0.37610200   -0.60326300
H                 -3.33500200   -2.73923400    1.40927900

请教各位老师应该怎么修改才能正确计算呢

zjxitcc

g_write是硬盘被写满了，要清理硬盘；不是说写个关键词再算一次 企图希望它不写满。当然，不排除你硬盘不够大，即使全清理了拿来算CCSD(T)也还不够用。

另外提醒一下，当前体系带-2负电荷，用cc-pVTZ计算没啥意义，基组一定要加弥散函数，可以给所有原子都用上弥散函数，也可以只给阴离子加上弥散函数。

sobereva

此外，不要做出
%nprocshared=28
%mem=4GB
这么荒唐的设定，尤其还是高级别后HF计算，好好看下文
Gaussian的安装方法及运行时的相关问题
http://sobereva.com/439（http://bbs.keinsci.com/thread-10814-1-1.html）

桃子量子

sobereva 发表于 2023-6-1 05:01
此外，不要做出
%nprocshared=28
%mem=4GB

好的老师，我的内存设的太小了对叭，核数大的话，内存应该相应增大，
老师我这么算是可以是不

桃子量子

zjxitcc 发表于 2023-5-31 12:06
g_write是硬盘被写满了，要清理硬盘；不是说写个关键词再算一次 企图希望它不写满。当然，不排除你硬盘不够 ...

好的，谢谢老师，我再下去好好看看，对于阴离子，一定要加弥散函数吗

zjxitcc

桃子量子 发表于 2023-6-2 09:19
好的，谢谢老师，我再下去好好看看，对于阴离子，一定要加弥散函数吗

一定

桃子量子

zjxitcc 发表于 2023-6-2 09:59
一定

计算的时候加上弥散函数，是不是就不可以使用CCSD(T)/cc-pvtz计算单点能了呢

桃子量子

zjxitcc 发表于 2023-6-2 09:59
一定

加上aug-，是不是就可以了呢，使用ccsd算单点能，
这样的话计算时间会大大增加，如果不使用CCSD去计算，请教是否还有别的基组可以代替呢

zjxitcc

桃子量子 发表于 2023-6-2 10:28
加上aug-，是不是就可以了呢，使用ccsd算单点能，
这样的话计算时间会大大增加，如果不使用CCSD去计算， ...

我看你在一开始提问时讲的是CCSD(T)，现在咋又说CCSD，方法名称要严谨，一个字符都不能漏。

可以用aug-cc-pVTZ，还可以降低弥散函数的使用量，见http://sobereva.com/119。另外，在2L已经说过了，可以只给涉及阴离子的原子加上弥散函数。

如果嫌慢，可以换用ORCA做DLPNO-CCSD(T)计算，可以使用与高斯中一样的基组。顺带一提，不要回复我说目前只能用高斯（因为这句话不成立，ORCA是免费程序，谁都可以下载使用）

桃子量子

zjxitcc 发表于 2023-6-2 10:42
我看你在一开始提问时讲的是CCSD(T)，现在咋又说CCSD，方法名称要严谨，一个字符都不能漏。

可以用aug ...

抱歉老师，不严谨了，好的谢谢老师