计算化学公社

标题: Gaussian wB97XD算nics出现l301错误是什么原因？ [打印本页]

作者
Author: aaq2800 时间: 2023-4-26 08:42
标题: Gaussian wB97XD算nics出现l301错误是什么原因？
本帖最后由 aaq2800 于 2023-4-26 09:53 编辑

这个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

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

作者
Author: aaq2800 时间: 2023-4-26 08:54

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

好的，多谢多谢

作者
Author: aaq2800 时间: 2023-4-26 09:00

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

作者
Author: sobereva 时间: 2023-4-26 09:31
算NICS的目的，鬼原子根本不需要定义基组

作者
Author: aaq2800 时间: 2023-4-26 09:43

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

谢谢社长，我再试试

作者
Author: mfdsrax2 时间: 2023-4-27 09:53
Sr已经是第五周期了，def2-TZVP是赝势基组，需要读取赝势

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