计算化学公社

标题: gaussian计算Zn原子单点能出现无具体报错信息的L914错误 [打印本页]
作者
Author: Stepheee    时间: 2022-9-19 10:43
标题: gaussian计算Zn原子单点能出现无具体报错信息的L914错误
本帖最后由 Stepheee 于 2022-9-20 12:44 编辑

使用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

  6. Title Card Required

  8. 0 1
  9. Zn 0.0 0.0 0.0


复制代码
输出文件末尾段如下：
  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
复制代码
完整的输出文件在附件中。



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

老师好，基组是为了和文献一致，目前在复现他的数据
作者
Author: Stepheee    时间: 2022-9-19 14:04
zjxitcc 发表于 2022-9-19 13:08
Zn原子是闭壳层，SCF很好收敛，没有必要加scf=fermi。另外这个基组有点浪费计算资源了，用def2-TZVP更好， ...

老师您好，按照您说的修改输入文件以后依旧会如上报错，输入文件的关键字是：
  1. #p blyp/def2qzvp stable
复制代码
作者
Author: zjxitcc    时间: 2022-9-19 14:14
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

  6. title

  8. 0 1
  9. Zn   0.0   0.0   0.0

  11. --Link1--
  12. %chk=Zn_BLYP.chk
  13. %mem=16GB
  14. %nprocshared=16
  15. #p BLYP/def2QZVP nosymm guess=read geom=allcheck
复制代码

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

我用G16B算了一遍，没有报错。
作者
Author: Stepheee    时间: 2022-9-19 15:15
hebrewsnabla 发表于 2022-9-19 14:18
我用G16B算了一遍，没有报错。

老师好，那请问为什么会出现这种情况呢，和我用的超算集群有关系吗？
作者
Author: LittlePupil    时间: 2022-9-20 00:13
一看C.03吓我一跳结果是C.01...
看样子是程序bug.
作者
Author: Stepheee    时间: 2022-9-20 12:44
LittlePupil 发表于 2022-9-20 00:13
一看C.03吓我一跳结果是C.01...
看样子是程序bug.

不好意思，写错了
作者
Author: 七尺贱    时间: 2022-9-20 13:18
zjxitcc 发表于 2022-9-19 13:08
Zn原子是闭壳层，SCF很好收敛，没有必要加scf=fermi。另外这个基组有点浪费计算资源了，用def2-TZVP更好， ...

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

这只能算了对比才知道。



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3