利用molpro计算H2S阳离子2B1平衡位置的频率出错，请大家帮忙，谢谢

yzsong

大家好，我在利用molpro计算H2S阳离子2B1（自旋多重度是2）平衡位置的频率出错，请大家帮忙解答。
我的计算步骤是：
1. 在C2v对称性下优化H2S阳离子2B1，得到平衡位置的结构；
2. 利用优化得到的平衡位置的结构来计算频率；
但是总是得到下面的错误，这个错误好像是提示计算频率时不能用对称性，但是如果不用对称性的话，怎么确定该分子的2B1态呢？我的输入文件已上传，请有该方面经验的帮忙解决这个问题，非常感谢。


PROGRAM * HESSIAN

*** Long output written to logfile /tmp/10607.node00/HSH-2B1-freq-AVQdZ.log ***


Computing numerical hessian using default procedure for command MRCI

Symmetry elements for unique displacements: X  Y

Using no symmetry in wavefunction calculations

Running default procedure: CASSCF000  MRCI

method=MRCI  nstsym= 1  istsym(1)=2  isyref=2

Numerical hessians for multireference methods can only be computed without symmetry.
Use NOSYM!
? Error
? Illegal use of symmetry
? The problem occurs in numhess_mppx

sobereva

反正2B1是基态，结构稍微位移破坏了对称性时，正常情况也会收敛到这个态上。也可以比较一下计算数值Hessian在每一次有限差分过程中得到的能量值，因为位移很小，应该和平衡结构下2B1态能量相差较小，这也说明是在同一个势能面上。

yzsong

sobereva 发表于 2016-4-14 04:14
反正2B1是基态，结构稍微位移破坏了对称性时，正常情况也会收敛到这个态上。也可以比较一下计算数值Hessian ...

谢谢，我试试。

xujian

请问前辈，您用molpro做过几何优化吗？