在过渡态计算中，出现716错误如何解决

gesehuapa

老师好，我在做过多态计算时候，出现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.

gesehuapa

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

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

zjxitcc

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

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

gesehuapa

谢谢您的回复。我按照你说的保存了笛卡尔坐标，但是很快就出错了。
输出文件说没有变量。难道过渡态计算只能用内坐标？
输入文件
%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

wzkchem5

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