Useful extremum principle for the variational calculation of matrix elements

Abstract
Variational principles for the estimation of the matrix element Wnn(φn,Wφn) for an arbitrary operator W are of great interest. The variational estimates are constructed from a trial wave function φnt, an approximation to the nth normalized bound-state eigenfunction φn, and of a trial auxiliary function Lt, an approximation to L which satisfies (HEn)L=(WnnW)φnq(φn). Variational-principle applications have been limited by the difficulty of obtaining a reasonable Lt, among other things, one demands that Lt approach L as φnt approaches φn. The equation (HEnt)Lt=q(φnt), where Ent=(φnt,Hφnt), is known not to provide such an Lt. A practical procedure for handling complicated systems given a reasonably accurate Rayleigh-Ritz trial function φnt is called for. This paper provides such a procedure using techniques developed in the establishment of variational bounds on scattering lengths. Given H and φnt, we define Lt by ALt=q(φnt), where A differs from HEn in that the influence of states 1 through n has effectively been "subtracted out"; the operator A is non-negative. A functional M(Ltt) is constructed which is an extremum for Ltt=Lt. Variational parameters contained in

This publication has 16 references indexed in Scilit: