We propose expanding the space of variations in traditional variational calculations for the energy by considering the wave function ψ to be a functional of a set of functions χ : ψ = ψ [χ], rather than a function. A constrained search in a subspace over all functions χ such that the functional ψ[χ] satisfies a sum rule or leads to a physical observable is then performed. An upper bound to the energy is subsequently obtained by variational minimization. The rigorous construction of such a constrained-search-variational wave function functional is demonstrated.One of the mostly extensively employed and accurate approximation methods in quantum mechanics is the variational principle for the energy. Consider a quantum mechanical system with Hamiltonian operatorĤ. The ground state eigenenergies E and eigenfunctions Ψ for this system satisfy the Schrödinger equationĤΨ = EΨ. . In application of the variational principle, however, the space of variations is limited by the choice of form of the function chosen for the approximate wave function. For example, if Gaussian or Slater-type orbitals are employed, then the variational space is limited to such functions. In this paper we propose the idea of overcoming this limitation by expanding the space over which the variations are performed. This then allows for a greater flexibility for the structure of the approximate wave function. We demonstrate the idea of expansion of the variational space by example.We expand the space of variations by considering the approximate wave function to be a functional of the set of functions χ: ψ = ψ[χ], rather than a function. The space of variations is expanded because the functional ψ[χ] can be adjusted through the function χ to reproduce any well behaved function. However, this space of variations is still too large for practical purposes, and so we consider a subset of this space. In addition to the function ψ being of a particular analytical form and dependent on the variational parameters c i , the functions χ are chosen such that the functional ψ[χ] satisfies a constraint. Examples of such constraints on the wave function functional ψ[χ] are those of normalization or the satisfaction of the Fermi-Coulomb hole charge sum rule, or the requirement that it lead to observables such as the electron density, the diamagnetic susceptibility, nuclear magnetic constant or any other physical property of interest. A constrained-search over all functions χ such that ψ[χ] satisfies a particular condition is then performed. With the functional ψ[χ] thus determined, the functional I[ψ [χ]] is then minimized with respect to the parameters c i . In this manner both a particular system property of interest as well as the energy are obtained accurately, the latter being a consequence of the variational principle. We refer to this way of determining an approximate wave function as the constrained-search-variational method.As an example of the method we consider its application to the ground state of the Helium atom. In atomic units e =h = m = 1, the non...