<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>