Let's say we want to run a CASSCF calculation of a radical, for example the cation of a neutral closed-shell molecule (therefore, we calculate an open-shell radical cation). MOLCAS, as well as probably every other program package out there, needs starting orbitals when performing a CASSCF calculation. For a neutral species, these are normally assumed to be the regular old Hartree-Fock orbitals. Usually, when dealing with a cation, it seems that the Hartree-Fock-SCF calculation is performed for the neutral species or the di-cationic species to yield the necessary HF starting-orbitals. They are in these cases the result of closed-shell calculations. MOLCAS also permits the use of a UHF calculation as the basis for a subsequent CASSCF calculation. Then, the natural orbitals of the unrestricted calculation are used.
Is the usage of the natural orbitals of an unrestricted HF calculation reasonable as a basis for a subsequent CASSCF calculation? If not, why? When is it a good idea / when is it a bad idea? What are the differences in using natural UHF orbitals in contrast to the SCF orbitals of the neutral or di-cationic species? What are the pitfalls?
Edit: I am aware (or at least it is my understanding) that conceptually it is no problem to perform the CASSCF using the natural UHF orbitals, as it could for example also be performed using only the alpha orbitals of the UHF calculation. My question is therefore aimed at finding out in which cases it is generally a good or bad idea to use the natural UHF orbitals.