This work further develops the calculation of QED effects in a finite Gaussian basis. We focus on the nonlinear $\alpha {\left(Z\alpha \right)}^{n\ge 3}$ contribution to the vacuum polarization density, computing the energy shift of $1s_{1/2}$ states in hydrogenlike ions. Our goal is to improve the numerical computations to achieve a precision comparable to that of Green's function methods reported in the literature. To do so, an analytic expression for the linear contribution to the vacuum polarization density is derived using Riesz projectors. Alternative formulations of the vacuum polarization density and their relation are discussed. The convergence of the finite Gaussian basis scheme is investigated, and the numerical difficulties that arise are characterized. In particular, an error analysis is performed to assess the method's robustness to numerical noise. We then report a strategy for computing the energy shift with sufficient precision to enable a sensible extrapolation of the partial-wave expansion. A key feature of the procedure is the use of even-tempered basis sets, allowing extrapolation towards the complete basis-set limit.

Close

PrevNext

View figure in article

• Relativistic & quantum electrodynamic effects in atoms, molecules,& ions