ORCA的错误信息：Expecting C(onstraint) in ScanConstraints

astrozheng

我想要固定两个dihedral angles， 然后计算其能量。

输入的inp文件如下：
# avogadro generated ORCA input file
# Basic Mode
#
! HF 6-31G*
! Opt TightSCF
%geom
      Constraints
        {D 6 8 9  10 -180.0  C}  ## phi angle
        {D 8 9 10 13 -180.0  C}  ## phi angle
      end
end

% pal nprocs 8
end

* xyz 0 1
   C       -2.98900        1.32900        0.54500
   H       -3.78200        1.68100       -0.09600
   H       -2.47200        2.17500        0.98400
   H       -3.41400        0.73400        1.34600
   C       -2.03400        0.48600       -0.27500
   O       -2.22200        0.20100       -1.44200
   N       -0.92500        0.07200        0.41100
   H       -0.86700        0.24300        1.39300
   C        0.01800       -0.88300       -0.16400
   H       -0.25300       -0.96600       -1.20800
   C       -0.07700       -2.25700        0.50800
   H       -1.06400       -2.67500        0.35600
   H        0.11000       -2.17200        1.57200

......（一部分xyz数据省略）
*

然后遇到了下面的错误信息：
Your calculation utilizes the basis: 6-31G
Cite in your paper:
H - He: W.J. Hehre, R. Ditchfield and J.A. Pople, J. Chem. Phys. 56,
Li - Ne: 2257 (1972).  Note: Li and B come from J.D. Dill and J.A.
         Pople, J. Chem. Phys. 62, 2921 (1975).
Na - Ar: M.M. Francl, W.J. Pietro, W.J. Hehre, J.S. Binkley, M.S. Gordon,
         D.J. DeFrees and J.A. Pople, J. Chem. Phys. 77, 3654 (1982)
K  - Zn: V. Rassolov, J.A. Pople, M. Ratner and T.L. Windus, J. Chem. Phys.
         (accepted, 1998)
Note: He and Ne are unpublished basis sets taken from the Gaussian program


Error - Expecting C(onstraint) in ScanConstraints


Error in [GEOM] block - Scan_Constraints - Line 10 (180.0)


LEAVING ORCA



请问各位了解ORCA的朋友，这个错误信息来自哪？如何避免这个错误信息？

astrozheng

非常感谢sob。

确实可行。

一方面，dihedral的角度要不为零，但可以大于180度。
另外，一句后面的##的内容需要删掉。 但是，有时候在有的语句后面的comment却不会报错，这就有点奇怪了。

总之，问题得到了解决。

sobereva

astrozheng 发表于 2016-5-5 12:02
非常感谢sob！

我试了下，用20度的二面角来算，还是有错误信息。

把## phi 去掉，从提示看读到了phi字符。

astrozheng

非常感谢sob！

我试了下，用20度的二面角来算，还是有错误信息。

Error in [GEOM] block - Scan_Constraints - Line 10 (PHI)

LEAVING ORCA


但是考虑到如果直接生成20度的二面角的三维坐标，会让分子结构存在一些问题，而且也需要其他软件来做，就会有些麻烦。

在ORCA的网站找到的manual也没有讲太多constraint的细节，好忧伤。

sobereva

二面角不能直接给负值，否则会报这种错
你可以先把坐标摆成你要限制的值，然后这么设定约束的时候不写初值，默认就会用坐标里当前的值。