  A simple procedure to incorporate one-loop quantum electrodynamic (QED)
corrections into the generalized (Gatchina) nonlocal shape-consistent
relativistic pseudopotential model is described. The pseudopotentials for Lu,
Tl, and Ra replacing only inner core shells (with principal quantum numbers
$n\le 3$ for the two former elements and $n\le 4$ for the latter one) are
derived from the solutions of reference atomic SCF problems with the
Dirac-Coulomb-Breit Hamiltonian to which the model Lamb shift operator added.
QED contributions to atomic valence excitation energies evaluated at the SCF
level are demonstrated to exceed the errors introduced by the pseudopotential
approximation itself by an order of magnitude. Pilot applications of the new
model to calculations of excitation energies of two-valence-electron atomic
systems using the intermediate-Hamiltonian relativistic Fock space coupled
cluster method reformulated here for incomplete main model spaces are reported.
Implications for high-accuracy molecular excited state calculations are
discussed.
