Two-loop vacuum polarization in a Coulomb field

The leading-order two-loop vacuum-polarization potential, linear in the Coulomb field of a nucleus, was first derived in the seminal 1955 work of Källén and Sabry [ , 1 (1955)]. The higher-order two-loop vacuum-polarization corrections, however, have remained unknown until now. In this work, we compute Coulomb corrections to the Källén-Sabry potential, specifically those involving three, five, and seven Coulomb interactions inside the vacuum-polarization loop. The potentials are evaluated in momentum space and subsequently used to calculate one-electron energy shifts. Our results reduce the theoretical uncertainty of the two-loop vacuum-polarization contribution to transition energies, which is required for next-generation tests of bound-state QED in heavy one- and few-electron ions as well as for the determination of nuclear charge radii.