We present a new perturbation theory for quantum mechanical energy eigenstates when the potential equals the sum of two localized, but not necessarily weak potentials V1(r) and V2(r), with the distance L between the respective centers of the two taken to be quite large. It is assumed that complete eigenfunctions of the local Hamiltonians (i.e., in the presence of V1(r) or V2(r) only) are available as inputs to our perturbation theory. If the two local Hamiltonians have degenerate bound-state energy levels, a systematic extension of the molecular orbital theory (or the tight-binding approximation) follows from our formalism. Our approach can be viewed as a systematic adaptation of the multiple scattering theory to the problem of bound states.
The partial-wave Schrödinger equation with real boundary conditions is recast as an equation for the probability density. When a small additional potential is included, the changes in the bound-state energy eigenvalues are obtained, up to third order in the perturbation, purely in terms of the perturbing potential and the unperturbed probability density. Although the approach is different, our results are equivalent to those derived by Bender (Bender C M 1978 Advanced Mathematical Methods for Scientists and Engineers (New York: McGraw-Hill) p 330). Knowledge of neither the unperturbed energy spectrum nor the wavefunctions of excited states is required. Evaluations of the second-order energy shift are given for some soluble S -wave problems.