计算化学公社

标题: 求助！Unrecognized atomic symbol Title ;L301错误。确定波函数稳定性时 [打印本页]

作者
Author: Scx 时间: 2018-1-6 16:17
标题: 求助！Unrecognized atomic symbol Title ;L301错误。确定波函数稳定性时
求键能过程中，确定波函数稳定性时，总是link301错误，
输出：Center    Atomic    Atomic          Coordinates (Angstroms)
Number    Number    Type          X          Y          Z
---------------------------------------------------------------------
   1       47          0       0.000000 0.000000 0.000000
---------------------------------------------------------------------
General basis read from cards:  (5D, 7F)
Unrecognized atomic symbol Title
Error termination via Lnk1e in /home/sunchuanxin/g09/l301.exe at Sun Jan  7 01:05:29 2018.
Job cpu time:    0 days  0 hours  0 minutes  4.0 seconds.
File lengths (MBytes):  RWF=    27 Int=    0 D2E=    0 Chk=    2 Scr=    2
输入：
%chk=Ag-M2-I0-B2PLYP-stable.chk
%mem=8GB
%nproc=16
# B2PLYP/genecp stable=opt guess=read geom=allcheck IOp(3/125=0360003600,3/78=0640006400,3/76=0350006500,3/77=1000010000,5/33=1) empiricaldispersion=gd3bj int(grid=ultrafine)

Title Card Required

0 2
Ag             -0.56306310 0.87837837 0.00000000

Ag    0
S 2 1.00
   19.000000000          -0.16600104141
   15.428199933          0.35665095918
S 1 1.00
   6.0553507268          1.0000000
S 1 1.00
   1.4162368935          1.0000000
S 1 1.00
   0.61758635858       1.0000000
S 1 1.00
   0.10474197431       1.0000000
S 1 1.00
   0.37685106264E-01          1.0000000
P 4 1.00
   13.188180180          0.66928737147E-01
   7.7952789138       -0.24735235409
   2.0351571912          0.49154280216
   0.98093914842       0.49741609006
P 1 1.00
   0.44451179958       1.0000000
P 1 1.00
   0.13000000000       1.0000000
P 1 1.00
   0.41200000000E-01          1.0000000
D 4 1.00
   25.784397351          0.35645063082E-02
   11.396636755          -0.12984262784E-01
   2.7345581361          0.24108826548
   1.1873583605          0.42412330744
D 1 1.00
   0.47316910566       1.0000000
D 1 1.00
   0.16746017986       1.0000000
F 1 1.00
   1.3971100             1.0000000
****

Ag    0
Ag-ECP    3    28
f-ul potential
  2
2    14.2200000          -33.68992012
2    7.1100000          -5.53112021
s-ul potential
  4
2    13.1300000          255.13936452
2    6.5100000          36.86612154
2    14.2200000          33.68992012
2    7.1100000             5.53112021
p-ul potential
  4
2    11.7400000          182.18186871
2    6.2000000          30.35775148
2    14.2200000          33.68992012
2    7.1100000             5.53112021
d-ul potential
  4
2    10.2100000          73.71926087
2    4.3800000          12.50211712
2    14.2200000          33.68992012
2    7.1100000             5.53112021

谢谢！

作者
Author: sobereva 时间: 2018-1-6 16:21
去掉geom=allcheck

作者
Author: jiangning198511 时间: 2018-1-6 16:27
去掉 guess=read geom=allcheck 试试

作者
Author: Scx 时间: 2018-1-6 17:12

sobereva 发表于 2018-1-6 16:21
去掉geom=allcheck

谢谢老师！问题解决了。

作者
Author: Scx 时间: 2018-1-6 17:12

jiangning198511 发表于 2018-1-6 16:27
去掉 guess=read geom=allcheck 试试

谢谢前辈！问题解决了

作者
Author: Scx 时间: 2018-1-6 17:33

sobereva 发表于 2018-1-6 16:21
去掉geom=allcheck

请教老师，在算完波函数稳定性后，继续下一步优化时，又出现了这样的错误是怎么回事儿呢？
输入文件
%chk=Ag-M2-I0-B2PLYP-stable-opt.chk
%mem=8GB
%nproc=16
# B2PLYP/genecp freq(noraman) opt IOp(3/125=0360003600,3/78=0640006400,3/76=0350006500,3/77=1000010000,5/33=1) empiricaldispersion=gd3bj int(grid=ultrafine)
本来加着guess=read geom=allcheck 出错后，依次去掉试了试，仍然是这样。
查看的输出文件
tail Ag-M2-I0-B2PLYP-stable-opt.log
         Y1          0.00000 0.00000
         Z1          0.00000 0.00000 0.00000
ITU=  0  0
   Eigenvalues ---
En-DIIS/RFO-DIIS IScMMF=       0 using points:    2 1
RFO step:  Lambda=-5.96046448D-09.
RFO step:  ModMin= 1 Lambda=-5.96046448D-09 EMin= 0.00000000D+00
Using NR instead of RFO step for point 1.
RFO step:  ModMin= 1 Lambda=-5.96046448D-09 EMin= 0.00000000D+00
Using NR instead of RFO step for point 2.

作者
Author: sobereva 时间: 2018-1-6 18:34

Scx 发表于 2018-1-6 17:33
请教老师，在算完波函数稳定性后，继续下一步优化时，又出现了这样的错误是怎么回事儿呢？
输入文件
%c ...

没看到报错提示
noraman完全多余

作者
Author: Scx 时间: 2018-1-16 10:06

sobereva 发表于 2018-1-6 18:34
没看到报错提示
noraman完全多余

谢谢老师！noraman的确是多余，我去掉之后，问题仍没有解决，查看输出文件，发现这是计算到中途自己停止了。我试着把核数和内存降低一些还是不行。恳请老师给指导一下！谢谢老师！
输出文件：
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Using GEDIIS/GDIIS optimizer.
Search for a local minimum.
Step number 2 out of a maximum of 20
All quantities printed in internal units (Hartrees-Bohrs-Radians)
Mixed Optimization -- En-DIIS/RFO-DIIS
Swapping is turned off.
Update second derivatives using D2CorN and points 2
The second derivative matrix:
                        X1       Y1       Z1
         X1          0.00000
         Y1          0.00000 0.00000
         Z1          0.00000 0.00000 0.00000
ITU=  0  0
   Eigenvalues ---
En-DIIS/RFO-DIIS IScMMF=       0 using points:    2 1
RFO step:  Lambda=-5.96046448D-09.
RFO step:  ModMin= 1 Lambda=-5.96046448D-09 EMin= 0.00000000D+00
Using NR instead of RFO step for point 1.
RFO step:  ModMin= 1 Lambda=-5.96046448D-09 EMin= 0.00000000D+00
Using NR instead of RFO step for point 2.

输入文件：
%chk=Ag-M2-I0-B2PLYP-stable-opt.chk
%mem=2GB
%nproc=4
# B2PLYP/genecp empiricaldispersion=gd3bj freq guess=read opt IOp(3/125=0360003600,3/78=0640006400,3/76=0350006500,3/77=1000010000,5/33=1)

Title Card Required

0 2
Ag             -0.56306310 0.87837837 0.00000000

Ag    0
S 2 1.00
   19.000000000          -0.16600104141
   15.428199933          0.35665095918
S 1 1.00
   6.0553507268          1.0000000
S 1 1.00
   1.4162368935          1.0000000
S 1 1.00
   0.61758635858       1.0000000
S 1 1.00
   0.10474197431       1.0000000
S 1 1.00
   0.37685106264E-01          1.0000000
P 4 1.00
   13.188180180          0.66928737147E-01
   7.7952789138       -0.24735235409
   2.0351571912          0.49154280216
   0.98093914842       0.49741609006
P 1 1.00
   0.44451179958       1.0000000
P 1 1.00
   0.13000000000       1.0000000
P 1 1.00
   0.41200000000E-01          1.0000000
D 4 1.00
   25.784397351          0.35645063082E-02
   11.396636755          -0.12984262784E-01
   2.7345581361          0.24108826548
   1.1873583605          0.42412330744
D 1 1.00
   0.47316910566       1.0000000
D 1 1.00
   0.16746017986       1.0000000
F 1 1.00
   1.3971100             1.0000000
****

Ag    0
Ag-ECP    3    28
f-ul potential
  2
2    14.2200000          -33.68992012
2    7.1100000          -5.53112021
s-ul potential
  4
2    13.1300000          255.13936452
2    6.5100000          36.86612154
2    14.2200000          33.68992012
2    7.1100000             5.53112021
p-ul potential
  4
2    11.7400000          182.18186871
2    6.2000000          30.35775148
2    14.2200000          33.68992012
2    7.1100000             5.53112021
d-ul potential
  4
2    10.2100000          73.71926087
2    4.3800000          12.50211712
2    14.2200000          33.68992012
2    7.1100000             5.53112021

作者
Author: sobereva 时间: 2018-1-16 10:12

Scx 发表于 2018-1-16 10:06
谢谢老师！noraman的确是多余，我去掉之后，问题仍没有解决，查看输出文件，发现这是计算到中途自己停止 ...

这是基本常识问题
一个原子怎么优化？优化什么？

作者
Author: Scx 时间: 2018-1-16 10:37

sobereva 发表于 2018-1-16 10:12
这是基本常识问题
一个原子怎么优化？优化什么？

优化频率，没有虚频。在进行键能计算的过程中，对每一个分子，原子的计算都必须确保经历了这三个步骤（第一步确定基态构型，第二步确定分子波函数稳定性，第三步再次进行opt的优化，使其在矫正了波函数之后进一步做能量的优化计算），并最终没有虚频，对称性正确才可以接着计算。这个是师兄师姐留下的资料，这样算不对吗？老师

作者
Author: sobereva 时间: 2018-1-16 15:15

Scx 发表于 2018-1-16 10:37
优化频率，没有虚频。在进行键能计算的过程中，对每一个分子，原子的计算都必须确保经历了这三个步骤（第 ...

“师兄师姐留下的资料”这类东西，总是错误百出。很多人严重的知识误区、对概念和程序的误解、关键词的滥用等等，全都是网上乱七八糟地方道听途说，或者师兄师姐传下来的。

一个原子怎么可能优化？根本都没有要被优化的变量，程序不报错才怪。
做什么事都绝对不要盲目效仿，一定要搞清楚原理，明白了为什么，再去计算，懂了原理自然能保证计算不犯错、出了错能马上明白是因为什么。

作者
Author: Scx 时间: 2018-1-16 16:15

sobereva 发表于 2018-1-16 15:15
“师兄师姐留下的资料”这类东西，总是错误百出。很多人严重的知识误区、对概念和程序的误解、关键词的 ...

嗯嗯，谢谢您的指导！

作者
Author: lao7 时间: 2019-10-30 17:06

Scx 发表于 2018-1-16 16:15
嗯嗯，谢谢您的指导！

单个原子无法优化。优化指的是原子与原子之间形成的结构。没有结构，无从优化。即使做了，也不收敛，报错！

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