<div dir="ltr"><div><div><div><div><div>Dear Molpro Experts,<br><br></div>I have used following input to get CASPT3 energies and population analysis:<br><br>***,PSB3<br>memory, 450, M<br>SET, CHARGE=1<br>gexpec, dm<br>gdirect<br>geomtyp=xyz<br>symmetry,nosym<br>geometry={<br>14<br>PSB3<br>C         -0.0002566337        0.0012806335       -0.0045289738<br>C         -0.0030717618       -0.0004852300        1.3489960614<br>C          1.2442813204       -0.0012856474        2.0756133424<br>C          1.4475005421       -0.0039130467        3.4356099462<br>C          0.3842489326       -0.0058656017        4.3728298332<br>N          0.5682361422       -0.0095509528        5.6709166926<br>H          2.4654050705       -0.0048405031        3.8177998498<br>H          2.1426826325       -0.0002301120        1.4591434928<br>H         -0.6530291517       -0.0042819065        4.0474832878<br>H         -0.9514175607       -0.0016008345        1.8796761931<br>H         -0.2127208483       -0.0110877687        6.3195543342<br>H          1.5004228539       -0.0110165037        6.0776308342<br>H         -0.9240107645        0.0015068704       -0.5729784019<br>H          0.9255382266        0.0022226034       -0.5733504920<br>}<br>basis=6-31G*<br>rhf<br>put, molden, PSB3.scf.molden;<br>{multi<br>  rotate, 19.1, 20.1, 0<br>  rotate, 25.1, 28.1, 0<br>  occ, 25, 0<br>  closed, 19, 0<br>  wf, 44, 1, 0<br>  state, 3<br>}<br>put, molden, PSB3.rasscf.molden;<br>{rs3, shift=0.1, ipea=0.25<br>  state, 3<br>}<br>pop<br>!<br><br></div>however I got population analysis for one state only. Is it possible to get it<br></div>for ground and two first excited states in one run?<br><br></div>Best regards,<br><br></div>Dawid Grabarek<br><br></div>