In the basic description of electronic systems, each electron is taken to be subject only to the average interaction with the other electrons, i.e. a mean field approximation. The many electron wave function is then constructed from one-electron wavefunctions that are often referred to as ‘orbitals’. A variational optimization of an antisymmetrized product of orbitals gives Hartree-Fock (HF) theory which in some sense appears to be the optimal mean field theory. While molecular structure is quite well predicted by HF and the orbital energies can be used to estimate electronic excitations and photoelectron spectra, the total energy and, thereby, bond energy is not accurate. Kohn-Sham density functional theory (KS-DFT), is much more widely used today. It is also a mean field theory based to some extent on orbitals but gives a better estimate of the total energy, while orbital energies are inaccurate. HF and KS-DFT have several other shortcomings, for example an overemphasis on electron localization by the former and an overemphasis on delocalization by the latter. Is it possible to develop a mean field theory that gives accurate total energy as well as useful orbital energies and a proper balance between localized and delocalized states, thereby providing an optimal mean field theory, more accurate than either KS-DFT and HF? Since our intuition is based strongly on the concept of orbitals and the computational effort of higher level methods increases rapidly with the number of electrons, it seems worthwhile to seek such an optimal mean field theory. A promising approach appears to be an energy functional with explicit orbital density dependence (ODD). This expanded functional form makes it possible to avoid the self-interaction error inherent in KS-DFT while the computational effort still scales as the number of electrons cubed. Functionals formed by applying the Perdew-Zunger self-interaction correction [1] to KS-DFT functionals are examples of ODD functionals. While they do not represent optimal mean field theory, they give encouraging results and could be used to guide the development of such a theory. Various results obtained recently using a fully variational implementation of PZ-SIC applied to gradient dependent KS-DFT functionals and complex valued orbitals will be given in this presentation and some suggestions for further work towards an optimal ODD functional discussed. The calculations presented include atomization energy of molecules [2], localized and delocalized charge states in a cation [3], dipole bound anion [4], Rydberg excited states of molecules [5], electron holes in crystalline oxides [6] and magnetic properties of small transition metal clusters. [1] J.P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). [2] S. Lehtola, E.Ö. Jónsson and H. Jónsson, Journal of Chemical Theory and Computation 12, 4296 (2016). [3] X. Cheng, Y. Zhang, E. Jónsson, H. Jónsson and P.M. Weber, Nature Communications 7, 11013 (2016). [4] Y. Zhang, P. M. Weber and H. Jónsson, J. Phys. Chem. Letters 7, 2068 (2016). [5] Y. Zhang, S. Deb, H. Jónsson and P.M. Weber, J. Phys. Chem. Letters 8, 3740 (2017). [6] H. Gudmundsdóttir, E. Ö. Jónsson and H. Jónsson, New Journal of Physics 17, 083006 (2015).