Se trata de una cuestión compleja, sobre todo porque la gente suele pensar en términos de una imagen de partícula independiente (es decir, el aufbau llenando orbitales), aunque la función de onda exacta de muchos cuerpos tiene fuertes correlaciones electrón-electrón. Así que permítame reformular su pregunta:
¿Cuál es la relación entre las funciones propias de KS y la función de onda exacta de muchos cuerpos?
Matemáticamente, como dices, las eigenfunciones KS no tienen estrictamente ningún significado físico (que sepamos). Sin embargo, las eigenfunciones KS ofrecen una imagen cualitativa (y a veces cuantitativa) muy útil. La razón es que las funciones propias de KS son una aproximación bastante buena a algo de la teoría de perturbaciones de muchos cuerpos llamado función de onda de cuasipartículas. La función de onda de las cuasipartículas es una propiedad física bien definida de un sistema que, esencialmente, nos dice que si añadimos (o quitamos) un electrón con una cierta cantidad de energía, a dónde irá. Por ejemplo, véase Phys. Rev. B 74, 045102 (2006) .
¿Existen ejemplos de cuando las eigenfunciones KS no dan una buena descripción de las funciones de onda de las cuasipartículas? Bueno, ciertamente hay muchas situaciones en las que las aproximaciones que solemos utilizar en la DFT (como la aproximación de la densidad local) conducen a graves problemas. Sin embargo, no conozco ningún ejemplo en el que alguien haya demostrado que las funciones propias exactas de KS (es decir, las obtenidas con el verdadero funcional de intercambio-correlación) no concuerden, al menos cualitativamente, con las funciones de onda de las cuasipartículas.
Como apunte, todo lo que he dicho arriba se aplica igualmente a las funciones de onda de Hartree-Fock. De hecho, existe una sólida base matemática para interpretar las funciones de onda HF como una aproximación a las funciones de onda de las cuasipartículas. Véase el capítulo 4 de la obra de Fetter Teoría cuántica de los sistemas de muchas partículas .
¿Y los valores propios de KS? En sentido estricto, en general no se corresponden con las energías de ionización (ni con ninguna otra cantidad físicamente útil). La única excepción es el mayor valor propio ocupado, que es exactamente igual a la energía de ionización del sistema. El teorema de Janak nos dice que los demás valores propios están relacionados con la derivada de la energía con respecto a la ocupación de esa función propia:
$$\epsilon_i=\frac{dE}{dn_i}$$
Ver Phys. Rev. B 18, 7165 (1978) y Phys. Rev. B 56, 16021 (1997) . Resulta que, empíricamente, estos valores propios son, no obstante, aproximaciones bastante buenas a los verdaderos niveles de energía del sistema, con algunas salvedades. En particular, las brechas de banda de los sólidos se subestiman sistemáticamente.