计算化学公社

标题: ORCA的错误信息：Expecting C(onstraint) in ScanConstraints [打印本页]

作者
Author: astrozheng 时间: 2016-5-4 17:44
标题: ORCA的错误信息：Expecting C(onstraint) in ScanConstraints
我想要固定两个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的朋友，这个错误信息来自哪？如何避免这个错误信息？

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

作者
Author: astrozheng 时间: 2016-5-5 12:02
非常感谢sob！

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

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

LEAVING ORCA

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

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

作者
Author: sobereva 时间: 2016-5-5 19:50

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

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

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

作者
Author: astrozheng 时间: 2016-5-5 21:51
非常感谢sob。

确实可行。

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

总之，问题得到了解决。

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