Please, correct me if i am wrong. I have some high symmetric molecule with transition metal. I am thinking about CASSCF for the molecule. I did DFT calculations and found all orbitals which were proposed by group-theoretical analysis and used them as starting orbitals for following calculations. The active space is too large: 18e19o. Since i am going to use cistep=guga for CASSCF i decided to run cityp=guga for my active sapce (iexcit=2, nstate=16, wstate=1,-0) for my large active space to estimate the population of the orbitals. I am interesting mainly in ground state, but of course low laying excited states also will be considered. After cityp=guga run the population of the natural orbital in atomic basis was analyzed. I found that population of two the lowest orbital and two the highest orbitals within my active significantly differ from the other orbitals. Does it mean i can remove two valence and two virtual orbitals from the active space and thus contract my active space 18e19o->14e15o for the casscf run? Does this methodology make sense?
[ This message was edited on Sat May 18 '13 at 4:14am by the author ]