For my localization scheme, I need to get MO coefficients for the full molecule, and for atomic fragments. The problem is that for the atomic fragments like carbon, oxygen, …, the occupied p orbitals are not spherically averaged. In the code DIRAC, one can specify fractional electron number in those orbitals (like 2 electrons in 6 spin-orbitals for carbon) to get this spherical averaging. It is called “Open shell SCF calculation using Average-of-Configuration”. Basically, instead of optimizing the energy for a single open-shell state, it optimizes the energy for a limited set of open-shell states.
Is there anything like this in psi4 ?
I am reopening this question, i.e. is there any way to get fractional occupation to have spherical symmetry of atomic densities (using Restricted-HF). As @hokru said, this is possible with UHF but the function is disabled in RHF.
I guess it should be too much of a problem to add it in the restricted case as well ?
Problem is that when I do calculations on atoms, I get MO coefficients (well, should be called orthonormal AO coefficients I guess as it is for a single atom) that are not unique and I can get really different results depending on which architecture my code is running.