<div dir="ltr"><div><div><div>Dear All,<br><br></div>I have performed geometry optimization with XMS-CASPT2 method using following input:<br><span style="font-family:monospace,monospace"><br>***,PSB3<br>memory, 600, M<br>SET, CHARGE=0<br>gexpec, dm<br>gdirect<br>geomtyp=xyz<br>symmetry,nosym<br>geometry={<br>6<br>eten<br>C     0.000000     0.000000    -0.000893<br>C     0.000000     0.000000     1.335893<br>H     0.923249     0.000000     1.905412<br>H    -0.923249     0.000000     1.905412<br>H    -0.923249     0.000000    -0.570412<br>H     0.923249     0.000000    -0.570412<br>}<br>basis=6-31G*<br>rhf<br>put, molden, PSB3.scf.molden;<br>{multi<br>  occ, 9, 0<br>  closed, 7, 0<br>  wf, 16, 1, 0<br>  state, 2<br>  weight, 1, 1<br>}<br>put, molden, PSB3.rasscf.molden;<br>{rs2, xms=2, mix=2, root=1;<br>}<br>optg<br>!<br><br></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">However when I change root=1 to root=2 I still get optimization of the ground state. Why is <br></font></span></div><div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">that?<br><br></font></span></div><div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">Best regards,<br></font></span></div><div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">Dawid Grabarek<br></font></span></div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif"></font></span><div><div><br></div></div></div>