Gaussian wB97XD算nics出现l301错误是什么原因？

aaq2800

这个nics计算用wB97XD就不行，换m062x就没问题，只是不好收敛。但结构搜索和结构优化，其他一些性质的计算都是用wB97XD。gaussian09w D01版本
输入文件
%chk=SrD18d1.chk
%mem=512MB
%nproc=1
#wB97XD/gen NMR

SrD18d1

0 1   
C       -3.461081   -1.220535   -0.002191
C        0.679707   -3.531581    0.001958
C        3.471088   -1.199078   -0.001387
C       -2.957265   -2.387954    0.000465
C        2.973715   -2.369255    0.003452
C       -1.812423   -3.037662   -0.002038
C       -0.656881   -3.537368    0.000968
C       -2.971257    2.368798    0.002549
C       -3.920979   -0.012406    0.001441
C        1.831940   -3.024414    0.001046
C       -3.468297    1.198616    0.000029
C       -1.831255    3.026993    0.002043
C       -0.679109    3.534391    0.001333
C        2.956188    2.384896    0.000559
C        3.462670    1.218938   -0.001750
C        3.924898    0.011447    0.002197
C        0.657522    3.537665    0.000125
C        1.812101    3.035811   -0.002380
Sr       -0.006399    1.376525   -0.001582
Bq       0.000282    0.000067    1.000211
Bq       0.000282    0.000067    2.000211

C Sr 0
def2TZVP
****
Bq 0
LanL2DZ
****

Bq 0
LanL2DZ


输出文件

Cannot cope with ghost atoms or with translation vectors.
Stoichiometry    C18Sr
Framework group  C1[X(C18Sr)]
NUMDOF-- NAT=    19 NAtoms=    21
Deg. of freedom     0
Full point group                 C1      NOp   1
Rotational constants (GHZ):      0.3367100      0.3263531      0.1657254
General basis read from cards:  (5D, 7F)
Warning:  center  20 has no basis functions!
Warning:  center  21 has no basis functions!
   591 basis functions,   976 primitive gaussians,   684 cartesian basis functions
    73 alpha electrons       73 beta electrons
       nuclear repulsion energy      1417.3874886977 Hartrees.
NAtoms=   21 NActive=   21 NUniq=   21 SFac= 1.00D+00 NAtFMM=   60 NAOKFM=F Big=F
Integral buffers will be    262144 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned off.
R6DC6: No C6 coefficient available for IA=  0
Error termination via Lnk1e in C:\G09W\l301.exe at Wed Apr 26 08:36:18 2023.
Job cpu time:       0 days  0 hours  0 minutes  0.0 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      1

paramecium86

这是因为色散矫正不支持加到bq原子上。 加上iop(3/124=4) 就行了。 社长也是有帖子的。 http://bbs.keinsci.com/thread-15540-1-1.html

aaq2800

paramecium86 发表于 2023-4-26 08:53
这是因为色散矫正不支持加到bq原子上。 加上iop(3/124=4) 就行了。 社长也是有帖子的。 http://bbs.keinsci ...

好的，多谢多谢

aaq2800

paramecium86 发表于 2023-4-26 08:53
这是因为色散矫正不支持加到bq原子上。 加上iop(3/124=4) 就行了。 社长也是有帖子的。 http://bbs.keinsci ...

但是又出现这种情况
Deg. of freedom     0
Full point group                 C1      NOp   1
Rotational constants (GHZ):      0.8335839      0.6434643      0.3631596
Standard basis: def2TZVP (5D, 7F)
Atomic number out of range in Df2TZV.
Error termination via Lnk1e in C:\G09W\l301.exe at Wed Apr 26 08:59:15 2023.
Job cpu time:       0 days  0 hours  0 minutes  0.0 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      4 Scr=      1

sobereva

算NICS的目的，鬼原子根本不需要定义基组

aaq2800

sobereva 发表于 2023-4-26 09:31
算NICS的目的，鬼原子根本不需要定义基组

谢谢社长，我再试试

mfdsrax2

Sr已经是第五周期了，def2-TZVP是赝势基组，需要读取赝势