计算化学公社

标题: 利用molpro计算H2S阳离子2B1平衡位置的频率出错，请大家帮忙，谢谢 [打印本页]

作者
Author: yzsong 时间: 2016-4-12 11:01
标题: 利用molpro计算H2S阳离子2B1平衡位置的频率出错，请大家帮忙，谢谢
本帖最后由 yzsong 于 2016-4-12 11:02 编辑

大家好，我在利用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

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

作者
Author: yzsong 时间: 2016-4-14 08:13

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

谢谢，我试试。

作者
Author: xujian 时间: 2016-8-5 15:15
请问前辈，您用molpro做过几何优化吗？

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