  Multiple energy minima of the LDA+U energy functional are obtained for
$\gamma$-Ce when it is implemented in a full potential, rotationally invariant
scheme including spin-orbit coupling, and different starting local
configurations are chosen. The lowest energy solution leads to a fully spin
polarized 4f state and the lattice constant of $\gamma$-Ce. The higher energy
local minima (additional self-consistent solutions) are shown to be strongly
indicative of crystal electric field and multiplet excitations.
