  The possibility to use functionals of occupation numbers and natural orbitals
for interacting fermions is discussed as an alternative to multi-reference
energy density functional method. An illustration based on the two-level Lipkin
model is discussed.
