Determination of wave-function functionals: The constrained-search variational method

In a recent paper [Phys. Rev. Lett. 93, 130401 (2004)], we proposed the idea of expanding the space of variations in variational calculations of the energy by considering the approximate wave function {psi} to be a functional of functions {chi}, {psi}={psi}[{chi}], rather than a function. A constrained search is first performed over all functions {chi} such that the wave-function functional {psi}[{chi}] satisfies a physical constraint or leads to the known value of an observable. A rigorous upper bound to the energy is then obtained via the variational principle. In this paper we generalize the constrained-search variational method, applicable to both ground and excited states, to the determination of arbitrary Hermitian single-particle operators as applied to two-electron atomic and ionic systems. We construct analytical three-parameter ground-state functionals for the H{sup -} ion and the He atom through the constraint of normalization. We present the results for the total energy E, the expectations of the single-particle operators W={sigma}{sub i}r{sub i}{sup n}, n=-2,-1,1,2, W={sigma}{sub i}{delta}(r{sub i}), and W={sigma}{sub i}{delta}(r{sub i}-r), the structure of the nonlocal Coulomb hole charge {rho}{sub c}(rr{sup '}), and the expectations of the two particle operators u{sup 2},u,1/u,1/u{sup 2}, where u= vertical barr{sub i}-r{sub j} vertical bar. The results for all the expectation values are remarkably accurate when compared with the 1078-parameter wave function of Pekeris, and other wave functions that are not functionals. We conclude by describing our current work on how the constrained-search variational method in conjunction with quantal density-functional theory is being applied to the many-electron case.