<div dir="ltr">I am trying to locate the conical intersection for the methoxy molecule at the C3v geometry. I have used nosym option in geometry input, as it gives me a lower energy for the molecule at the end of a geometry optimisation than without. With the optimised UHF function I performed a CASSCF energy calculation. Now I would like to do a conical intersection optimisation using this wavefunction and geometry. The trouble is, since I have used the nosym option, I am unable to specify the symmetry of the two states that I need to calculate, to do a state-averaged calculation. As far as the manual goes, the only way that I can give CPMCSCF command is with specification of the symmetry corresponding to each electronic state. Is there a way around this problem?<br>
Following is my input - <br><br>***,methoxy calculation.<br>memory,128,m<br>file,2,<a href="http://met_cas_aug.wf">met_cas_aug.wf</a>,new;<br>punch,met_cas_aug.pun<br><br>gprint,basis,orbitals,civector<br><br>basis=aug-cc-pVDZ<br>
geomtype={zmat}<br>geometry={nosym,Angstrom;<br>          c;<br>          o,c,oc;<br>          h1,c,h1c,o,h1co;<br>          h2,c,h1c,o,h1co,h1,120.0;<br>          h3,c,h1c,o,h1co,h1,240.0}<br>oc=1.410000, h1c=1.100000, h1co=109.100000<br>
<br>{uhf} !uhf input w/o symmetry<br>put,molden,met_uhf_aug.molden<br><br> {optg,space=zmat}<br> put,molden,met_opt_aug.molden<br><br> {casscf; occ,9;closed,2;<br>  start,orbital=2200.2;wf,17,1,1<br>  CPMCSCF,NACM,1.1,2.1,accu=1.0d-7,record=5100.1 !If I give in 1.1,2.1, molpro terminates<br>
  CPMCSCF,GRAD,1.1,spin=0.5,accu=1.0d-7,record=5101.1<br>  CPMCSCF,GRAD,2.1,spin=0.5,accu=1.0d-7,record=5102.1}<br> put,molden,met_cas_aug.molden<br><br>  {Force<br>  SAMC,5100.1<br>  CONICAL,6100.1}<br><br>  {Force<br>  SAMC,5101.1<br>
  CONICAL,6100.1}<br><br>  {Force<br>  SAMC,5102.1<br>  CONICAL,6100.1}<br><br> optg,startcmd=multi<br> put,molden,met_cas_aug_con.molden<br><br>---<br clear="all">Thanks,<br>-- <br>Jayashree<br>
</div>