Abstract
A general computational scheme for the (non-relativistic) Bethe logarithm is developed opening the route to `routine' evaluation of the leading-order quantum electrodynamics correction (QED) relevant for spectroscopic applications for small polyatomic and polyelectronic molecular systems. The implementation relies on Schwartz' method and minimization of a Hylleraas functional. In relation with electronically excited states, a projection technique is considered, which ensures positive definiteness of the functional over the entire parameter (photon momentum) range. Using this implementation, the Bethe logarithm is converged to a relative precision better than 1:10³ for selected electronic states of the two-electron H₂ and H₃⁺, and the three-electron He₂⁺ and H+H₂ molecular systems. The present work focuses at nuclear configurations near the local minimum of the potential energy surface, but the computations can be repeated also for other structures.