[molpro-user] Geometry specification with Cs symmetry

[Jayashree]

Hi all,
I am trying to obtain two states of the methoxy radical - one belonging to
A' and the other belonging to A" symmetry. The optimized geometry
corresponding to A' symmetry is a true minimum and that of A" is a saddle
point. I have the same starting geometry specification for the two states,
but the symmetry of the open shell orbital and thereby the final
wavefunction specification is different in the two cases. The optimisation
is successful for A' state, but in the case of A" state, I get the following
error:

 *Attempt to include non-existent symmetry X
 Because of noorient option, program will
 not align to your requested symmetry elements
 Try again with different symmetry request
 This error can also happen in geometry optimizations that specify
  symmetry-breaking internal coordinates in terms of active variables
 User-specified symmetry elements: X
 Symmetry elements:
 ? Error
 ? Symmetry error
 ? The problem occurs in zmatrix.f:zmat_evaluate*

I believe that in a zmatrix input, MOLPRO automatically reorients the
molecule for the requested symmetry. If I have not specified the noorient
option in my input(below), why does MOLPRO assume it and thereby not
reorient?

*Input for A" state-*
*geomtype={zmat}
 geometry={Angstrom;
           q1;
           o,q1,q1o;
           h1,q1,q1h1,o,h1q1o;
           h2,q1,q1h2,o,h2q1o,h1,dih1;
           h3,q1,q1h3,o,h3q1o,h1,dih2;
           c,q1,q1c,o,180.0,h1,0.0}
 q1o=   0.7481460595672777      ,
 q1c=   0.7277119994327224      ,
 q1h1=    1.479915970144683      ,
 q1h2=    1.479915970144683      ,
 q1h3=    1.479915970144683      ,
 h1q1o=    147.9783986150496      ,
 h2q1o=    130.4352205872918      ,
 h3q1o=    130.4352205872918      ,
 dih1=    120.0000000000000      ,
 dih2=    240.0000000000000
 dummy,q1
{uhf;occ,7,2;open,2.2;wf,17,2,1}
{casscf;occ,7,2;closed,5,1;wf,17,2,1;maxiter,40}
{optg;root,2;active,q1h1,q1h2,q1h3}
*
Any suggestion would be helpful.
Thanks,
-- 
Jayashree
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20090318/b67e61ad/attachment.html>