Hi all,<br>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:<br>
<br> <i>Attempt to include non-existent symmetry X<br> Because of noorient option, program will<br> not align to your requested symmetry elements<br> Try again with different symmetry request<br> This error can also happen in geometry optimizations that specify<br>
symmetry-breaking internal coordinates in terms of active variables<br> User-specified symmetry elements: X<br> Symmetry elements:<br> ? Error<br> ? Symmetry error<br> ? The problem occurs in zmatrix.f:zmat_evaluate</i><br>
<br>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?<br>
<br><u>Input for A" state-</u><br><i>geomtype={zmat}<br> geometry={Angstrom;<br> q1;<br> o,q1,q1o;<br> h1,q1,q1h1,o,h1q1o;<br> h2,q1,q1h2,o,h2q1o,h1,dih1;<br> h3,q1,q1h3,o,h3q1o,h1,dih2;<br>
c,q1,q1c,o,180.0,h1,0.0}<br> q1o= 0.7481460595672777 ,<br> q1c= 0.7277119994327224 ,<br> q1h1= 1.479915970144683 ,<br> q1h2= 1.479915970144683 ,<br> q1h3= 1.479915970144683 ,<br>
h1q1o= 147.9783986150496 ,<br> h2q1o= 130.4352205872918 ,<br> h3q1o= 130.4352205872918 ,<br> dih1= 120.0000000000000 ,<br> dih2= 240.0000000000000<br> dummy,q1<br>{uhf;occ,7,2;open,2.2;wf,17,2,1}<br>
{casscf;occ,7,2;closed,5,1;wf,17,2,1;maxiter,40}<br>{optg;root,2;active,q1h1,q1h2,q1h3}<br></i><br clear="all">Any suggestion would be helpful.<br>Thanks,<br>-- <br>Jayashree<br>