计算化学公社

标题: 高斯计算单点能L801 L804报错问题 [打印本页]

作者
Author: 桃子量子 时间: 2023-5-31 11:56
标题: 高斯计算单点能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

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

作者
Author: zjxitcc 时间: 2023-5-31 12:06
g_write是硬盘被写满了，要清理硬盘；不是说写个关键词再算一次企图希望它不写满。当然，不排除你硬盘不够大，即使全清理了拿来算CCSD(T)也还不够用。

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

作者
Author: sobereva 时间: 2023-6-1 05:01
(, 下载次数 Times of downloads: 12)

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

作者
Author: 桃子量子 时间: 2023-6-2 09:18

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

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

作者
Author: 桃子量子 时间: 2023-6-2 09:19

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

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

作者
Author: zjxitcc 时间: 2023-6-2 09:59

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

一定

作者
Author: 桃子量子 时间: 2023-6-2 10:25

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

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

作者
Author: 桃子量子 时间: 2023-6-2 10:28

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

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

作者
Author: zjxitcc 时间: 2023-6-2 10:42
本帖最后由 zjxitcc 于 2023-6-2 11:00 编辑

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

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

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

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

作者
Author: 桃子量子 时间: 2023-6-2 10:54

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

可以用aug ...

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

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