Abstract:In this work, we propose an efficient and accurate computational method to evaluate the many-potential $\alpha\left(Z\alpha\right)^{n\ge3}$ vacuum polarization density of hydrogen-like atoms within the finite-basis approximation of the Dirac equation. To prove the performance of our computational method, we choose to work with the one-electron $_{\,\,\,92}^{238}\text{U}$ atom. In summary, we find that compliance with charge conjugation symmetry is a priori required to obtain physical results that are in line with our knowledge of the analytical problem. We also note that the final numerical results are found to be in excellent agreement with previous formal analytical (and numerical) evaluations that are limited to a few simple nuclear distribution models. Our technique can be efficiently implemented and evaluated in codes that solve the radial Dirac equation in the finite basis set framework and allows the use of arbitrary (radial) nuclear charge distribution. The obtained numerical results of the non-perturbative vacuum polarization density automatically account for the extended nuclear size effect. This method is hence of special importance for atomic Dirac problems whose analytical Green's functions expressions are not at hand or have relatively complicated analytical forms. Furthermore, we propose a vacuum polarization density formula that forces compliance with charge conjugation symmetry and can be used in cases where the relativistic basis violates this symmetry, as is the case in most relativistic basis set programs. In addition, we have shown that vector components of the vacuum polarization four-current vanish in the case where the Dirac Hamiltonian is symmetric under time-reversal symmetry.

Finite-Basis-Set Approach to the Two-Center Heteronuclear \\ Dirac Problem

Ab initio QED calculations in diatomic quasimolecules

Calculating the many-potential vacuum polarization density of the Dirac equation in the finite-basis approximation

Solving One-Electron Systems in a Novel Gaussian-Sinc Mixed Basis Set

Analytical Structure Matching and Very Precise Approach to the Coulombic Quantum Three-Body Problem

Theoretical Investigation of the Spectroscopic Constants for the Ground-State Diatomic Species Cu2, Ag2, and Au2

Finite Space Complete Basis Method: Precision Computation of High-Resolution Spectrum Near Ionization Threshold

Three-Center Tight-Binding together with Multipolar Auxiliary Functions

Two-center harmonic oscillator basis for Skyrme-DFT calculations (I): formalism and Proof of Principle

Light one-electron quasi-molecular ions within the finite-basis-set method for the two-center Dirac equation

Computational method for general multicenter electronic structure calculations

Correlation Consistent Gaussian Basis Sets for H, B–Ne with Dirac–Fock AREP Pseudopotentials: Applications in Quantum Monte Carlo Calculations

Precision of Finite-Difference Representation in 3D Coordinate-Space Hartree-Fock-Bogoliubov Calculations

The critical notes to the solutions of the "puzzle" of hyperfine structure in 82+ and 80+ $^{209}$Bi ions

Investigating the basis set convergence of diagrammatically decomposed coupled-cluster correlation energy contributions for the uniform electron gas

A complete basis set model chemistry. VII. Use of the minimum population localization method

Odd nuclei and quasiparticle excitations within the Barcelona Catania Paris Madrid energy density functional

The Nuclear Magnetic Moment of 208bi and Its Relevance for a Test of Bound-State Strong-Field QED

High-precision solution of the Dirac Equation for the hydrogen molecular ion using a basis-set expansion