l401报错Inaccurate quadrature in CalDSu

hujingwen

在计算时遇到了L401（U2OF7-1.log)和L502(U2OF10-a.log)的报错，报错原因显示的都是Inaccurate quadrature in CalDSu.输入文件U2OF7-1如下：
%chk=U2OF7-1.chk
%nprocshared=8
%mem=4GB
# opt freq=raman b3lyp/genecp

Title Card Required

-1 3
U                  2.07117850    0.34422403    0.00000000
F                  0.09118185    0.34738265    0.00181089
F                  2.73414377    1.27550414    1.61666150
F                  2.72819866   -1.52358785   -0.00180887
O                  2.75785706    1.31530840   -1.68198327
U                  3.45786917    2.30524877   -3.39662641
F                  5.32462348    1.97190103   -2.82697070
F                  2.52764569    3.95636581   -2.82313116
F                  2.52133835    0.98747948   -4.53977737
F                  4.11788060    3.23862112   -5.01328994

U 0
SDD
****
F O 0
6-311+g*
****

U 0
SDD

输出文件U2OF7-1报错如下：
Spurious integrated density or basis function:
NE=  135 NElCor=    0 El error=2.16D-04 rel=8.46D-07 Tolerance=1.00D-03
Shell    74     absolute error=3.39D-02              Tolerance=1.20D-02
Shell    74       signed error=3.39D-02              Tolerance=1.00D-01
Inaccurate quadrature in CalDSu.
Error termination via Lnk1e in /home/hujingwen/Program/g09/l401.exe at Fri Mar 19 14:30:17 2021.
Job cpu time:  0 days  0 hours  0 minutes 55.2 seconds.
File lengths (MBytes):  RWF=     62 Int=      0 D2E=      0 Chk=      2 Scr=      2

g09 job started: Fri Mar 19 14:30:00 CST 2021
g09 job finished: Fri Mar 19 14:30:20 CST 2021
Elapsed time: 20


输入文件U2OF10-a如下：
%chk=U2OF10-a.chk
%nprocshared=8
%mem=4GB
# opt freq=raman b3lyp/genecp int=ultrafine SCF=novaracc

Title Card Required

-4 3
U                  0.00000000    0.00000000    2.24095700
F                 -0.12250000    2.30038700    2.50125800
F                  2.30038700    0.12250000    2.50125800
F                  0.00000000    0.00000000    4.62385600
F                  0.12250000   -2.30038700    2.50125800
F                 -2.30038700   -0.12250000    2.50125800
U                  0.00000000    0.00000000   -2.24095700
F                  2.30038700   -0.12250000   -2.50125800
F                  0.12250000    2.30038700   -2.50125800
F                  0.00000000    0.00000000   -4.62385600
F                 -2.30038700    0.12250000   -2.50125800
F                 -0.12250000   -2.30038700   -2.50125800
O                  0.00000000    0.00000000    0.00000000

U 0
SDD
****
F O 0
6-311+g*
****

U 0
SDD

输出文件U2OF10-a报错如下：

Density matrix breaks symmetry, PCut= 1.00D-04
Density has only Abelian symmetry.
Density matrix breaks symmetry, PCut= 1.00D-04
Density has only Abelian symmetry.
Spurious integrated density or basis function:
NE=  166 NElCor=    0 El error=6.37D-01 rel=2.26D-03 Tolerance=1.00D-03
Shell    80     absolute error=1.02D-02              Tolerance=1.20D-02
Shell   115       signed error=1.08D-05              Tolerance=1.00D-01
Inaccurate quadrature in CalDSu.
Error termination via Lnk1e in /home/hujingwen/Program/g09/l502.exe at Thu Mar 18 13:20:32 2021.
Job cpu time:  0 days  0 hours  3 minutes 56.6 seconds.
File lengths (MBytes):  RWF=    226 Int=      0 D2E=      0 Chk=      5 Scr=      2

g09 job started: Thu Mar 18 13:19:58 CST 2021
g09 job finished: Thu Mar 18 13:20:34 CST 2021
Elapsed time: 36


在采用了L401报错的那个ppt的解决方法以及结构调整后仍然无法解决问题
求助大家如何解决

sobereva

初始结构很诡异。U-F、U-O都不应当会那么长
另外也尝试int=superfine

描述时不要用什么“L401报错的那个ppt”这种其他人都看不懂的说法，明确说清楚到底都尝试过什么

paramecium86

我尝试用你贴上来的原始的输入文件在我这边试了试（改到40核心 20GB内存）。我试了 G16 A03 ，C01， G09 D01 都没遇到你这个错误。（计算了20分钟左右最终都是 “常见”的SCF不收敛错误）感觉会不会不是关键词的问题。 记得以前论坛里有大佬发过一个帖子提到这错误可能跟chk文件出了问题有关系

hujingwen

sobereva 发表于 2021-3-19 18:32
初始结构很诡异。U-F、U-O都不应当会那么长
另外也尝试int=superfine

谢谢sob老师，但是改了结构和关键词依旧会报错。我还尝试了先用小基组算，但是也还是不行。请问关于铀的配合物是否需要一些特别的手段呢？

hujingwen

paramecium86 发表于 2021-3-19 18:40
我尝试用你贴上来的原始的输入文件在我这边试了试（改到40核心 20GB内存）。我试了 G16 A03 ，C01， G09 D0 ...

我试了好多次了，最近还有一个别的铀的配合物的结构，依旧报相似的错误/(ㄒoㄒ)/~~

zjxitcc

hujingwen 发表于 2021-3-27 13:02
我试了好多次了，最近还有一个别的铀的配合物的结构，依旧报相似的错误/(ㄒoㄒ)/~~

输入文件里用#p而非#。然后还有什么错误信息再贴上来，大家分析一下就知道了。“相似的错误”有点宽泛了。

paramecium86

hujingwen 发表于 2021-3-27 13:02
我试了好多次了，最近还有一个别的铀的配合物的结构，依旧报相似的错误/(ㄒoㄒ)/~~

你有没有试过换一个计算的平台？至少我用你贴上来的输入文件 换了几台电脑 几个不同的版本，都还没遇见跟你一样的报错。

hujingwen

zjxitcc 发表于 2021-3-27 15:09
输入文件里用#p而非#。然后还有什么错误信息再贴上来，大家分析一下就知道了。“相似的错误”有点宽泛了 ...

Defaulting to unpruned grid for atomic number  92.
NRdTot=   11209 NPtTot=     3072740 NUsed=     3152841 NTot=     3152873
NSgBfM=   618   619   619   619   619 NAtAll=   125   125.
Leave Link  302 at Wed Mar 31 22:04:48 2021, MaxMem= 1073741824 cpu:      98.0
(Enter /home/hujingwen/Program/g09/l303.exe)
DipDrv:  MaxL=1.
Leave Link  303 at Wed Mar 31 22:04:50 2021, MaxMem= 1073741824 cpu:       2.2
(Enter /home/hujingwen/Program/g09/l401.exe)
Harris functional with IExCor=  402 diagonalized for initial guess.
ExpMin= 5.00D-03 ExpMax= 1.53D+03 ExpMxC= 3.22D+02 IAcc=3 IRadAn=         5 AccDes= 0.00D+00
HarFok:  IExCor=  402 AccDes= 0.00D+00 IRadAn=         5 IDoV= 1
ScaDFX=  1.000000  1.000000  1.000000  1.000000
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Defaulting to unpruned grid for atomic number  92.
Spurious integrated density or basis function:
NE=  440 NElCor=    0 El error=2.99D-05 rel=5.98D-08 Tolerance=1.00D-03
Shell   325     absolute error=1.70D-02              Tolerance=1.20D-02
Shell   325       signed error=1.70D-02              Tolerance=1.00D-01
Inaccurate quadrature in CalDSu.
Error termination via Lnk1e in /home/hujingwen/Program/g09/l401.exe at Wed Mar 31 22:05:22 2021.
Job cpu time:  0 days  0 hours  3 minutes 17.8 seconds.
File lengths (MBytes):  RWF=    115 Int=      0 D2E=      0 Chk=      2 Scr=      2

g09 job started: Wed Mar 31 22:04:02 CST 2021
g09 job finished: Wed Mar 31 22:05:25 CST 2021
Elapsed time: 83

zjxitcc

hujingwen 发表于 2021-4-1 10:28
Defaulting to unpruned grid for atomic number  92.
NRdTot=   11209 NPtTot=     3072740 NUsed=   ...

这次的输入文件呢？输入输出要结合着看