计算化学公社

标题: 在过渡态计算中，出现716错误如何解决 [打印本页]

作者
Author: gesehuapa 时间: 2021-8-30 17:12
标题: 在过渡态计算中，出现716错误如何解决
老师好，我在做过多态计算时候，出现716错误，是有一个夹角大于180度的问题。下面是我的输入文件，和log文件。请教各位老师，我需要做些什么来能解决这个问题呢？万分感谢。

输入文件%nosave
# opt=(calcfc,ts,ef,maxcycle=400) freq b3lyp/genecp

Title Card Required

0 1
Ti
H                1          B1
C                1          B2 2          A1
C                3          B3 1          A2 2          D1 0
H                4          B4 3          A3 1          D2 0
O                1          B5 3          A4 4          D3 0

B1          1.75393709
B2          2.09460710
B3          1.20550440
B4          1.05891832
B5          1.62994081
A1          152.71580457
A2          179.24897336
A3          179.04942098
A4          110.19314191
D1          71.20893838
D2          170.45729945
D3       -176.52049374

H C O       0
aug-cc-pvtz
****
Ti       0
SDD
****

Ti       0
SDD

输出.log文件的末尾
Cartesian Forces:  Max    0.038252470 RMS    0.015964647
Error on Z-matrix line number    4:  angle Alpha is outside the valid range of 0 to 180.
Conversion from Z-matrix to cartesian coordinates failed:
---------------------------------------------------------------------------------------------------
                        Z-MATRIX (ANGSTROMS AND DEGREES)
CD Cent Atom N1    Length/X       N2    Alpha/Y       N3       Beta/Z       J
---------------------------------------------------------------------------------------------------
   1    1  Ti
   2    2  H       1 1.743022(    1)
   3    3  C       1 2.087700(    2)    2  145.029(    6)
   4    4  C       3 1.207573(    3)    1  180.327(    7)    2 57.921( 10)    0
   5    5  H       4 1.059933(    4)    3  179.203(    8)    1  170.951( 11)    0
   6    6  O       1 1.621416(    5)    3  110.900(    9)    4 -183.936( 12)    0
---------------------------------------------------------------------------------------------------
Error termination via Lnk1e in /software/g09.d01/g09/l716.exe at Fri Aug 27 23:03:40 2021.

作者
Author: wzkchem5 时间: 2021-8-30 17:57
本帖最后由 wzkchem5 于 2021-8-30 10:59 编辑

把这个文件用gaussview打开，保存成笛卡尔坐标的输入文件，再试。另外注意：aug-cc-pVTZ基组比SDD大很多，一般金属的基组应该大于等于非金属的基组，否则会出现基组不平衡的问题，而且金属的计算误差一般占主导，一般来说非金属基组加得比金属大等于白加了

作者
Author: gesehuapa 时间: 2021-8-31 09:30
谢谢您的回复。我按照你说的保存了笛卡尔坐标，但是很快就出错了。
输出文件说没有变量。难道过渡态计算只能用内坐标？
输入文件
%nosave
# opt=(calcfc,ts,ef,maxcycle=400) freq b3lyp/genecp

Title Card Required

0 1
Ti                0.63631500 0.38844400 -0.05007800
H                2.12616600 1.01185700 0.63402400
C                -1.42328600 0.01128700 0.00627100
C                -2.61141700 -0.19032800 0.03688100
H                -3.65201800 -0.38470400 0.06287200
O                1.46689200 -1.01233500 0.01823800

H C O       0
aug-cc-pvtz
****
Ti       0
SDD
****

Ti       0
SDD

输出log文件
ERROR IN INITEF. NUMBER OF VARIABLES (  0) **
** INCORRECT (SHOULD BE BETWEEN 1 AND 50) *
Error termination via Lnk1e in /software/g09.d01/g09/l113.exe

作者
Author: zjxitcc 时间: 2021-8-31 09:39
本帖最后由 zjxitcc 于 2021-8-31 09:41 编辑

gesehuapa 发表于 2021-8-31 09:30
谢谢您的回复。我按照你说的保存了笛卡尔坐标，但是很快就出错了。
输出文件说没有变量。难道过渡态计算只 ...

（1）不要写ef；（2）末尾至少空2行，你直接复制粘贴到论坛上，空行数目有时会变，自己确认一下；（3）最好用#p而非#；（4）金属用的SDD，非金属却用aug-cc-pVTZ，内行人看了都会感觉别扭，充其量def2TZVP足矣。如果不懂怎么选基组，看《谈谈量子化学中基组的选择》http://sobereva.com/336

作者
Author: gesehuapa 时间: 2021-9-1 09:45

zjxitcc 发表于 2021-8-31 09:39
（1）不要写ef；（2）末尾至少空2行，你直接复制粘贴到论坛上，空行数目有时会变，自己确认一下；（3）最 ...

谢谢老师的回答，按照您的建议做了修改。这个过渡态问题已经解决。多谢

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