Calculation of the relativistic Bethe logarithm in the two-center problem

We present a variational approach to evaluate relativistic corrections of order α² to the Bethe logarithm for the ground electronic state of the Coulomb two-center problem. That allows us to estimate the radiative contribution at mα⁷ order in molecular-like three-body systems such as hydrogen molecular ions H₂⁺ and HD⁺ or antiprotonic helium atoms. While we get ten significant digits for the nonrelativistic Bethe logarithm, calculation of the relativistic corrections is much more involved, especially for small values of bond length R. We were able to achieve a level of three to four significant digits starting from R=0.2 bohr, which will allow us to reach 10^-10 relative uncertainty on transition frequencies.