gaussian计算Zn原子单点能出现无具体报错信息的L914错误

Stepheee

使用gaussian 16 C.01计算单重态Zn原子的单点能，出现L914错误，但是不同于论坛上的报错总结，在输出文件中没有看到具体的报错信息。

输入文件如下：

1. %nprocshared=32
   
2. %mem=40GB
   
3. %chk=ZnH_2_Zn_1_m06l.chk
   
4. #p m06l/6-311++g(3df,3pd) stable scf=fermi
   
5. 
   
6. Title Card Required
   
7. 
   
8. 0 1
   
9. Zn 0.0 0.0 0.0
   
10. 
    
11. 
    

复制代码

输出文件末尾段如下：

1. Iteration     5 Dimension   101 NMult    98 NNew      3
   
2. CISAX will form     3 AO SS matrices at one time.
   
3. NMat=     3 NSing=     3 JSym2X=-1.
   
4. Vector  26 SDiff=3.58D-01 TDiff=1.00D+00 assumed to be singlet.
   
5. Vector  27 SDiff=8.53D-01 TDiff=3.13D-02 assumed to be triplet.
   
6. Vector  28 SDiff=5.13D-01 TDiff=3.82D-01 assumed to be triplet.
   
7. Vector  29 SDiff=8.92D-01 TDiff=2.05D-01 assumed to be triplet.
   
8. Vector  30 SDiff=9.30D-01 TDiff=5.30D-04 assumed to be triplet.
   
9. Vector  31 SDiff=4.25D-04 TDiff=1.13D+00 assumed to be singlet.
   
10. Vector  32 SDiff=1.09D-04 TDiff=1.08D+00 assumed to be singlet.
    
11. Vector  35 SDiff=1.86D-04 TDiff=1.14D+00 assumed to be singlet.
    
12. Orthogonalize current: ORTVEC returned         100 roots
    
13. Error termination via Lnk1e in /es01/paratera/sce3678/software-sce3678/g16/l914.exe at Fri Sep 16 13:30:07 2022.
    
14. Job cpu time:       0 days  0 hours  7 minutes 26.4 seconds.
    
15. Elapsed time:       0 days  0 hours  0 minutes 14.7 seconds.
    
16. File lengths (MBytes):  RWF=     89 Int=      0 D2E=      0 Chk=      2 Scr=      1

复制代码

完整的输出文件在附件中。

zjxitcc

七尺贱 发表于 2022-9-20 13:18
借楼问一下，电子亲和能不加弥散算不准确，那么def2-QZVPP去计算电子亲合能还需要加弥散吗？

这只能算了对比才知道。

七尺贱

zjxitcc 发表于 2022-9-19 13:08
Zn原子是闭壳层，SCF很好收敛，没有必要加scf=fermi。另外这个基组有点浪费计算资源了，用def2-TZVP更好， ...

借楼问一下，电子亲和能不加弥散算不准确，那么def2-QZVPP去计算电子亲合能还需要加弥散吗？

Stepheee

LittlePupil 发表于 2022-9-20 00:13
一看C.03吓我一跳结果是C.01...
看样子是程序bug.

不好意思，写错了

LittlePupil

一看C.03吓我一跳结果是C.01...
看样子是程序bug.

Stepheee

hebrewsnabla 发表于 2022-9-19 14:18
我用G16B算了一遍，没有报错。

老师好，那请问为什么会出现这种情况呢，和我用的超算集群有关系吗？

hebrewsnabla

Stepheee 发表于 2022-9-19 14:04
老师您好，按照您说的修改输入文件以后依旧会如上报错，输入文件的关键字是：

我用G16B算了一遍，没有报错。

zjxitcc

Stepheee 发表于 2022-9-19 14:04
老师您好，按照您说的修改输入文件以后依旧会如上报错，输入文件的关键字是：

波函数稳定性没必要在QZ级别下做，过于耗时，合理的不同基组下波函数稳定性情况是一样的。在def2-TZVP下算单点、检验波函数稳定性，然后由def2-QZVP读入即可，例如

1. %chk=Zn_BLYP.chk
   
2. %mem=16GB
   
3. %nprocshared=16
   
4. #p BLYP/def2TZVP nosymm stable
   
5. 
   
6. title
   
7. 
   
8. 0 1
   
9. Zn   0.0   0.0   0.0
   
10. 
    
11. --Link1--
    
12. %chk=Zn_BLYP.chk
    
13. %mem=16GB
    
14. %nprocshared=16
    
15. #p BLYP/def2QZVP nosymm guess=read geom=allcheck
    

复制代码

Stepheee

zjxitcc 发表于 2022-9-19 13:08
Zn原子是闭壳层，SCF很好收敛，没有必要加scf=fermi。另外这个基组有点浪费计算资源了，用def2-TZVP更好， ...

老师您好，按照您说的修改输入文件以后依旧会如上报错，输入文件的关键字是：

1. #p blyp/def2qzvp stable

复制代码

Stepheee

zjxitcc 发表于 2022-9-19 13:08
Zn原子是闭壳层，SCF很好收敛，没有必要加scf=fermi。另外这个基组有点浪费计算资源了，用def2-TZVP更好， ...

老师好，基组是为了和文献一致，目前在复现他的数据

zjxitcc

Zn原子是闭壳层，SCF很好收敛，没有必要加scf=fermi。另外这个基组有点浪费计算资源了，用def2-TZVP更好，要是觉得不够还可以用def2-QZVP，比什么辣眼睛的6-311++g(3df,3pd)好多了。有时候不一定要解决报错，让输入文件更干净合理一些，计算就正常了。