Correct small-truncated excited state wave functions obtained via minimization principle for excited states compared/ opposed to Hylleraas-Undheim and McDonald higher ``roots''

We demonstrate that, if a truncated expansion of a wave function is large, then the standard excited states computational method, of optimizing one ``root'' of a secular equation, according to the theorem of Hylleraas, Undheim and McDonald (HUM), tends to the correct excited wave function, comparable to that obtained via our proposed minimization principle for excited states [J. Comput. Meth. Sci. Eng. 8, 277(2008)] (independent of orthogonality to lower lying approximants). However, if a truncated expansion of a wave function is small - that would be desirable for large systems - then the HUM-based methods may lead to an incorrect wave function - despite the correct energy (: according to the HUM theorem) whereas our method leads to correct, reliable, albeit small truncated wave functions. The demonstration is done in He excited states, using truncated series ``small'' expansions both in Hylleraas coordinates, and via standard configuration-interaction truncated ``small'' expansions, in comparison with corresponding ``large'' expansions. Beyond that, we give some examples of linear combinations of Hamiltonian eigenfunctions that have the energy of the 1^st excited state, albeit they are orthogonal to it, demonstrating that the correct energy is not a criterion of correctness of the wave function.