The energies of the 5d^{9}6s^{2}(D_{5/2},D_{3/2})-5d^{10}6s (S_{1/2}), 5d^{10}7p(P_{1/2},P_{3/2})-5d^{10}6s(S_{1/2}), 5d^{10}7p(P_{1/2},P_{3/2})-5d^{10}7s(S_{1/2}), 5d^{9}6s^{2}(D_{5/2},D_{3/2})-5d^{10}6s (S_{1/2}) transitions in Hg^{+} (Ry).

## Abstract

Fundamentals of gauge-invariant relativistic many-body perturbation theory (PT) with optimized ab initio zeroth approximation in theory of relativistic multi-electron systems are presented. The problem of construction of optimal one-electron representation is directly linked with a problem of the correct accounting for multielectron exchange-correlation effects and gauge-invariance principle fulfilling in atomic calculations. New approach to construction of optimal PT zeroth approximation is based on accurate treating the lowest order multielectron effects, in particular, the gauge-dependent radiative contribution for the certain class of photon propagator (for instance, the Coulomb, Feynman, Babushkin ones) gauge. This value is considered to be a typical representative of important multielectron exchange-correlation effects, whose minimization is a reasonable criteria in the searching for optimal PT one-electron orbital basis. This procedure derives an undoubted profit in the routine many-body calculations as it provides the way of refinement of the atomic characteristics calculations, based on the “first principles”. The relativistic density-functional approximation is taken as the zeroth one. There have taken into account all exchange-correlation corrections of the second order and dominated classes of the higher orders diagrams (polarization interaction, quasiparticles screening, etc.). New form of multi-electron polarization functional is used. As illustration, the results of computing energies, transition probabilities for some heavy ions are presented.

### Keywords

- relativistic many-body perturbation theory
- density-functional approximation
- exchange-correlation effects
- radiative transitions
- oscillator strengths
- heavy atoms

## 1. Introduction

Perturbation theory (PT) formalism has a long history in studying different multielectron (more generally, multifermion) systems, including different atomic, molecular, and nuclear properties. Really, one should say about formalism of the many-body PT as, a rule, usually it applies to studying different properties of the multiparticle systems, for instance, ionization and excitation energies, spectra, electron exchange-correlation energies, hyperfine structure, radiative and autoionization decay rates (transition probabilities, oscillator, and lines strengths), as well as the influence of an external electromagnetic fields. In the last few decades, the PT methods have been refined with a sophisticated and comprehensive approach of more correct treatment of the exchange-correlation effects, electron-nuclear dynamics, and so on [1–44]. Rephrasing the known interesting quote by Bartlett and Musiał [3, 4] and earlier by Wilson, one could say that * the PT methods are an emerging computational area that is sixty years ahead of lattice gauge theory… and a rich source of new ideas and new approaches to the computation of many fermion systems*. The old multibody quantum theoretical approaches often take place, which have been primarily developed in a theory of a superfluity and/or a superconductivity, and generally speaking in a theory of solids, became the powerful tools for developing new conceptions in many-body (multielectron) atomic, nuclear, and molecular calculations [1–7].

A number of the PT versions include a synthesis of cluster expansions, Brueckner’s summation of ladder diagrams, the summation of ring diagrams Gell-Mann, and an infinite-order generalization of manybody PT (Kelly, 1969; Ivanov-Tolmachev, 1969–1974, Bartlett and Silver, 1974–1976, etc.; see review in Ref. [7]). Using quantum-field methods in atomic and molecular theory allowed obtaining a very powerful approach for the correct treatment of the exchange-correlation effects in many-electron systems. In this context, it is useful to remind about such sophisticated methods as a coupled-cluster theory, the Green-functions method, configuration interaction methods, and so on. Only with this property are applications to solids or the electron gas possible, and, even for small atoms and molecules, its effects are numerically quite essential. When relativistic effects became essential in the studied multielectron (fermion) system, naturally it is necessary to formulate a formalism of the relativistic many-body PT. In the first attempts, an account for the relativistic effects had been reduced to treating the Darwin, mass-velocity, and spin-orbit effects, which have to be added to the nonrelativistic solution and provide different approximations lying between the Schrödinger equation and the four-component Dirac equation [2, 6, 7]. Among recent developments in this field, special attention should be given to two very general and important computer systems for relativistic and QED calculations of atomic and molecular properties developed in the Oxford, Troitsk, and other groups (known as ”GRASP,” ”Dirac,” ”BERTHA,” ”QED,” “Superatom,” etc.; Ref. [1–13] and references therein). For example, a new relativistic molecular structure theory within the QED framework with accounting of the electron correlation and higher-order QED effects has been formulated and further realized as the BERTHA program. The master system of equations includes the so-called Dirac-Hartree-Fock-Breit self-consistent field equations. The useful overview of the relativistic electronic structure theory is presented in Refs. [2, 7] from the QED point of view. The next important step is an adequate taking into account the QED corrections. This topic has been a subject of intensive theoretical and experimental interest.

Hitherto, most many-body PT studies concerned atoms with a simple electron-shell structure, namely atoms of the inert gases and atoms and ions with a single electron (or hole) above (or inside) the closed shells core. The fundamental limitation to extend the many-body procedure beyond such simple atomic systems arises from the complexity of any perturbation expansion if more than just one or two * effective* particles appear in the derivation (see detailed analysis in Ref. [5]). In order to overcome this difficulty, a number of different efficient techniques of complex expansions were developed. As a rule, the atomic PT expansions are expressed in terms of the Feynman-Goldstone diagrams in Rayleigh-Schrödinger PT formulation. Above the most popular and known versions of the PT formalism, one should mention formally exact relativistic many-body PT with the model zeroth approximation by Ivanova-Ivanov et al., relativistic PT with the Hartree-Fock (HF) or Dirac-Fock (DF) zeroth approximations by Johnson et al., Flambaum-Dzuba et al., Safronova and Safronova et al., Khetselius et al., and so on [9–38]).

The searching for the optimal one-electron zeroth representation is one of the oldest in the theory of multielectron atoms and, respectively, in the formulation of the effective PT formalism. Two decades ago, Davidson had pointed the principal disadvantages of the traditional representation based on the self-consistent field approach and suggested the optimal “natural orbitals” representation [11]. Nevertheless, there remain insurmountable computational difficulties in the realization of the Davidson program (see, e.g., Refs. [11, 12]). One of the simplified recipes represents, for example, a density functional theory (DFT) formalism [8]. Unfortunately, this approach does not provide a regular refinement procedure in the case of the complicated atom with few quasiparticles (QPs) (electrons or vacancies above a core of the closed electronic shells). The problem of construction of the optimal one-electron representation is tightly linked with the problem of the correct accounting for the multielectron exchange-correlation effects. In Refs. [47, 48], the PT lowest-order multielectron effects, in particular, the gauge-dependent radiative contribution (gauge-noninvariant) for the certain class of the photon propagator gauge is treated. This value is considered to be the typical representative of the multielectron exchange-correlation effects contribution. New fundamental idea has been proposed in Refs. [47, 48] in order to construct the optimal PT one-electron basis and is in minimization of the gauge-noninvariant contribution into a radiation width of atomic level. Such an approach allows to determine an effectiveness of accounting of the multielectron exchange-correlation effects and provides the practical way of the refinement of the atomic characteristics calculations, based on the “first principles.” Really, the known standard criterion of the multielectron computing quality in atomic spectroscopy is linked with a closeness of the atomic level radiation width values, calculated using two alternative forms of the transition operator (the “length” and the “velocity” forms). It is of special interest to verify the compatibility of the new optimization principle with the other requirements conditioning a “good” one-electron representation. We suppose that this point should be obligatory in formulation of the effective, optimal PT formalism.

In this chapter, we present the theoretical fundamentals of the gauge-invariant relativistic many-body PT with using the optimized one-QP representation in the theory of relativistic multielectron systems [21–23, 47, 48]. All exchange-correlation corrections of the second-order and dominated classes of the higher-orders diagrams (polarization interaction, QPs screening, etc.) [47–67] have been taken into account. As illustration of application of the presented PT formalism, we list the results of computing energies, transition probabilities (oscillator strengths) in some heavy atoms (ion of Hg+).

## 2. Relativistic many-body perturbation theory with optimized one-quasiparticle zeroth representation

### 2.1. General remarks

Our relativistic PT version is constructed on the same principles as the known formally exact PT with model zeroth approximation by Ivanova-Ivanov et al. [33–47]; however, there a few principal points, where our formalism differs from this known theory. At first, this is another definition of the zeroth approximation, namely within the relativistic DFT one [14–17, 19–22]. Second, this is an implementation of the principally new approach to construction of the optimized one-QP representation, which allows correctly to take into account a gauge invariance principle fulfilling.

In nonrelativistic theory of multielectron atoms, a powerful field approach for computing the electron energy shift Δ* E* of the degenerate states is known, which are usually present in the dense spectra of the complex relativistic atomic multielectron systems

**. The key algorithm of this approach includes construction of the secular matrix**(Tolmachev-Ivanov-Ivanova, 1969–1974)

*using the known Gell-Mann and Low adiabatic formula and its further diagonalization. The analogous approach using the Gell-Mann and Low formula with an electrodynamic scattering matrix has been developed in a theory of the relativistic atom [33–36]; however, the*M

*matrix elements in the relativistic representation are complex; the corresponding imaginary parts determine the values of radiation widths. According to Ref. [34], the total electron energy shift can be defined as follows:*M

Here, Γ is a radiation width of the atomic level (or a possibility ** P** of the radiation decay or transition: P = Γ. Within the general framework, the corresponding energies of a nondegenerated excited states and their radiation decay amplitudes can be determined by means of the computing and further diagonalization of the matrix

*. In Refs. [33–37], the ReΔ*M

*calculation procedure has been generalized for the case of nearly degenerate states, whose levels form a more or less compact group. Naturally, the matrix*E

*reduces to one term (Δ*M

*) in the case of well-identified and separated energy spectrum. The Gell-Mann and Low formula allow further to obtain the expansion of the*E

*elements into PT series on interelectron interaction and apply the standard Feynman diagrammatic technique. The corresponding PT series is as follows:*M

Here, ^{(0)} is the contribution of the PT all-orders vacuum diagrams (in fact, this is a real matrix, which determines only the general atomic levels shift); ^{(1)}, ^{(2)},and ^{(3)} are the contributions, which correspond to the one-, two- and three-QP PT diagrams, respectively. The diagonal matrix ^{(1)} can be easily calculated as it represents a sum of the one-QP contributions. Generally speaking, computing all the one-QP diagrams contributions within the PT formalism is the most simple procedure. The more complicated problem is computing the ^{(2)} and ^{(3)}contributions. Using the Feynman diagrams technique, the authors [33–38] have in detail analyzed the ^{(2)} contributions. Naturally, the fundamental point of the whole consideration is the definition of the PT zeroth approximation.

### 2.2. The perturbation theory zeroth approximation

We will describe an atomic multielectron system by the relativistic Dirac Hamiltonian (the atomic units are used) as follows [14, 15]:

where * Z* is a charge of nucleus,

α

_{i},

α

_{j}are the Dirac matrices,

ω

_{ij}is the transition frequency, and

*, a light velocity. The interelectron interaction potential (second term in Eq. (3)) takes into account the retarding effect and magnetic interaction in the lowest order on parameter*c

α

^{2}(

*is the fine structure constant). Let us note that in order to account for the nuclear finite size effect (in the zeroth approximation), one could describe a charge distribution in the atomic nucleus*α

*(*ρ

*) by the Gaussian or Fermi (another variant is relativistic mean-field theory of a nucleus) functions and write the Coulomb potential for the spherically symmetric nuclear density*r

*(*ρ

*) as [14]*r|R

Here, * R* is a nuclear radius. According to the known Ivanova-Ivanov et al. method of differential equations [33–36], computing the potential (20) can be reduced to solving the system of the differential equations. By the way, this method is used by us in further under computing the PT first- and second-order corrections. The zeroth-order Hamiltonian

H

_{0}and perturbation operator can be presented in the standard form as follows [7, 14, 15]:

where _{i}, one-electron energies, and _{MF} is the central field self-consistent potential of the Coulomb type. The latter can be taken in the form of the usual Dirac-Fock potential or even any appropriate model potential, which imitates an effect of the electron subsystem. Let us remind that in the relativistic PT by Ivanova-Ivanov et al., the consistent model (as a rule, empirical) potential was taken as _{MF.} In our PT version, we use the potential

Further as _{X}(* r*) we use the standard Kohn-Sham (KS) exchange potential as follows [8]:

The standard definition of the exchange potential in the density-functional theory is as follows:

In the relativistic multielectron theory with a Hamiltonian having a transverse vector potential (for describing the photons), one could determine the homogeneous density * ρ*(

*), construct the corresponding exchange energy*r

E

_{X}[

*(*ρ

*)], and introduce the following exchange potential [16]:*r

where

where * b* is the optimization parameter (for details, see below and Refs. [16–19, 47–49] too). Naturally, potential (6) is subtracted from the interelectron potential in Eq. (3) in the perturbation operator. The Dirac equations for

*and*F

*components can be written as [14] follows:*G

Here, _{nχ} is one-electron energy without the rest energy. The boundary values are defined by the first terms of the Taylor expansion:

The condition * f, g*→0 at

*→∞ determines the quantified energies of the state*r

E

_{nχ}. The system of Eq. (11) is numerically solved by the Runge-Kutta method (‘Superatom” package is used [7, 13–23, 34, 36, 47–67]).

### 2.3. The perturbation theory first- and second-orders corrections: correlation effects

In the PT first order, one should determine the matrix elements of the PT operator with the relativistic Coulomb-Breit potential, which are the contributions of the following type [36]:

where

The value of the _{λ} can be expressed through the radial Slater-like integrals and presented as a sum of the Coulomb and Breit parts:

In the nonrelativistic limit, there remains only the first term in Eq. (15) depending only on the large component * f*(

*) of the one-electron Dirac functions. For example, its imaginary part is as follows [36]:*r

The angular coefficient has only a real part:

Here, {_{1}_{3}} means that _{1} and _{3} must satisfy the triangle rule and the sum * λ* +

l

_{1}+

l

_{3}must be an even number. The rest terms in Eq. (16) include the small components of the Dirac functions. The tilde in Eq. (13) designates that the large radial component

*must be replaced by the small one*f

*, and instead of*g

j

_{i}<

l

_{i}and

j

_{i}>

l

_{i}. The Breit (magnetic) part can be expressed by the similar way (see details in Refs. [13–16]).

Then, exchange-correlation effects can be treated within the PT formalism as effects of the second and higher PT orders. In the second order, one should especially note the polarization and ladder diagrams. In Figures 1 and 2, we list some important diagrams of the second order describing the effects of the polarization interaction of quasiparticles and screening of the external quasiparticles (or antiscreening in the case, say, of an electron and a vacancy).

The polarization diagrams take into account the quasiparticle interaction through the polarizable core, and the ladder diagrams account for the immediate quasiparticle interaction. An effective approach to accounting the polarization contributions is in adding the effective two-QP polarizable operator into the first-order matrix elements. The corresponding polarization operator can be taken in the following form [50]:

where * X* is the numerical coefficient, and

*is the light velocity. The similar approximate potential representation has been received for the exchange polarization interaction of quasiparticles (see details in Refs. [7, 14–19]). The polarization potential Eqs. (18a)–(18c) generalizes the corresponding nonrelativistic operator, which has been derived in Ref. [36].*c

In order to take into account the ladder diagrams contributions as well as some of the three-quasiparticle diagram contributions in all PT orders, we use the special procedure, which includes a modification of the mean-filled potential, which describes the effects of screening (antiscreening) of the core potential of each QP by the others (see details in Refs. [7, 14–19, 33–38]). Introduction of the additional screening potential into the Dirac equations for the large and small components changes the 1-QP energies and orbitals. It results in the corresponding modification of the diagonal 1-QP matrix

### 2.4. Optimization of the relativistic orbitals basis

In order to obtain a precise description of the spectral characteristics of multielectron atomic systems, within the PT framework one should generate the optimized relativistic orbitals basis (see “Introduction” section) [1–7, 9–15]. The powerful ab initio approach to construction of the optimized PT basis has been developed in Ref. [48] and reduced to consistent treating gauge-dependent multielectron contributions ImΔ_{ninv} of the lowest relativistic PT corrections to the atomic level radiation width and their further functional minimization.

For simplicity, let us consider now the one-quasiparticle atomic system (i.e., atomic system with one electron or vacancy above a core of the closed electronic shells). The multiquasiparticle case does not contain principally new moments. In the PT lowest, second order for the Δ* E*, there is only one-quasiparticle Feynman diagram B (see Figure 3), contributing the ImΔ

*(the radiation decay width).*E

In the fourth order of QED PT (the second order of the atomic PT), the diagrams appear, whose contribution to the ImΔ_{ninv} accounts for the multielectron exchange-correlation (polarization) effects (diagrams _{d}, _{ex}; Figure 3). This multielectron contribution is dependent on the photon propagator gauge (the gauge-noninvariant contribution). Let us remind about the known criterion of the correctness of the atomic-computing radiation transition probabilities using the alternative forms for the transition operator (“length” and “velocity” transition operator forms). Their closeness of the “length” and “length” transition probabilities values confirms the correctness of the relativistic orbitals basis construction. Correspondingly, their noncoincidence is provided by multielectron by their nature and gauge-noninvariant terms.

In Ref. [48], the gauge-noninvariant contribution to an imaginary part of the electron energy has been calculated, which is as follows:

where * C* is the gauge constant, and

*is the boundary of the closed shells.*f

Realizing a principle of minimization of the functional ImΔ_{ninv,} one could obtain the Dirac-Kohn-Sham (DKS)-like equations for an electron density. Their numerical solution allows to obtain the optimized basis of the one-QP relativistic orbitals. The corresponding procedure is described in detail, for example, in Refs. [18–23]. All details of the presented PT formalism can be found in Refs. [7, 14–20, 47–49].

### 2.5. Radiation decay probability as an imaginary part of the electron energy shift. Method of calculation

The method of computing the radiation decay (transition probabilities, oscillator strengths) probabilities within the relativistic energy approach is presented in, for instance, Refs. [16–19, 33–35, 47, 48]. Here, we only note that a probability is directly linked with the imaginary part of electron energy shift, which is defined in the PT lowest order as follows:

where

The individual terms of the sum Eq. (21) represent the contributions of different channels and probability, for instance, of the dipole α-n transition as * λ*-pole transition is usually represented as a sum of the electric

*-pole transition*λ

*→*γ

*connects two states with parities which by*δ

*(or*λ

*+ 1) units. In our designations,*λ

In a case of the two-quasiparticle states (for instance, the excited atomic state is treated as a state with the two QP: electron and vacancy above the closed shells core), the corresponding probability has the following form (say, transition:

It is worth noting that all relativistic atomic calculations are usually carried out in the * jj*-coupling scheme. The transition to the intermediate-coupling scheme is realized by diagonalization of the

*matrix, but usually only Re*M

*should be diagonalized. The important simplified moment of the procedure is connected with converting the imaginary part by means of the matrix of eigenvectors {*M

C

_{mk}}, obtained by diagonalization of Re

*:*M

where _{ij} are the matrix elements in the * jj*-coupling scheme, and

M

_{mk}in the intermediate-coupling scheme representation. The procedure is correct to terms of the order of Im

*/Re*M

*.*M

In conclusion, let us also underline that the tedious procedure of phase convention in calculating the matrix elements of different operators is avoided in the energy approach, although the final formulae, certainly, must coincide with the formulae obtained using the traditional amplitude quantum-mechanical method. All other details can be found in Refs. [7, 16–19, 33–36, 47–50].

## 3. Some results and conclusions

As illustration of the application of the above presented formalism, we present the results of computing energies, transition probabilities (oscillator strengths) in the heavy multielectron ion of Hg^{+}. A great interest to studying similar systems (Hg) is explained by the importance of the corresponding data, for instance, for laser effect studying. The collision of atoms of the Mendeleev table second raw with ions of helium (other inert gases) leads to creating ions in the excited states which is important for creating the inverse populations and laser effect. The available literature data on radiative characteristics are definitely insufficient. An account of the relativistic and correlation effects has a critical role in the cited systems as the studied transitions occur in the external shells in a strong field of atom with large Z. Within the relativistic PT, the Hg^{+} states can be treated one- and three-QP states of electrons (6s) and vacancy (5d^{−1}) above the core of the closed shells 5d^{10}6s^{2}. The interaction “quasiparticle core” is described by the potential (6). The polarization interaction of the quasiparticles through the core is described by the two-particle effective potential Eqs. (18a)–(18c). All calculations are performed using the modified atomic code “Superatom-ISAN.”

In Tables 1–3, we present the experimental (NIST) [32] and theoretical energies, electric E1 (5d^{10}7p(P_{1/2},P_{3/2})-5d^{10}6s(S_{1/2}), 5d^{10}7p(P_{1/2},P_{3/2})-5d^{10}7s(S_{1/2})), and E2 (5d^{9}6s^{2} (D_{5/2}, D_{3/2})-5d^{10}6s (S_{1/2})) probabilities of the transitions in the spectrum of Hg^{+}. The theoretical results are obtained within the Hartree-Fock, Dirac-Fock methods by Ostrovsky-Sheynerman, relativistic PT theory with the empirical model potential zeroth approximation (RPT-MP) [18, 31], and our optimized RPT using relativistic energy approach (RPT-EA).

Method | _{6s} | 7P_{1/2}-6S_{1/2} | 7P_{3/2}-6S_{1/2} | 7P_{1/2}-7S_{1/2} | 7P_{3/2}-7S_{1/2} | D_{3/2}-S_{1/2} |
---|---|---|---|---|---|---|

HF | −1.07 | 0.721 | 0.721 | 0.095 | 0.095 | 0.863 |

DF | −1.277 | 0.904 | 0.922 | 0.109 | 0.127 | 0.608 |

RPT-MP | −1.377 | 0.986 | 1.019 | 0.114 | 0.147 | 0.462 |

RPT-EA | −1.378 | 0.987 | 1.020 | 0.115 | 0.148 | 0.462 |

Exp. | −1.378 | 0.987 | 1.020 | 0.115 | 0.148 | 0.461 |

Method | 7P_{3/2}-6S_{1/2} | 7P_{1/2}-6S_{1/2} | 7P_{3/2}-7S_{1/2} | 7P_{1/2}-7S_{1/2} | 7P_{3/2}-6S_{1/2} |
---|---|---|---|---|---|

HF | 4.75×10^{6} | 4.75×10^{6} | 3.65×10^{7} | 3.65×10^{7} | 3.65×10^{7} |

DF | 8.45×10^{7} | 1.67×10^{7} | 6.89×10^{7} | 6.89×10^{7} | 4.71×10^{7} |

DF (_{exp}) | 1.17×10^{8} | 2.04×10^{7} | 1.10×10^{8} | 1.10×10^{8} | 5.52×10^{7} |

RPT-MP | 1.49×10^{8} | 2.31×10^{7} | 1.41×10^{8} | 1.41×10^{8} | 6.33×10^{7} |

RPT-EA | 1.51×10^{8} | 2.33×10^{7} | 1.43×10^{8} | 1.43×10^{8} | 6.35×10^{7} |

Exp. | 1.53×10^{8} | 2.35×10^{7} | 1.44×10^{8} | 1.44×10^{8} | 6.37×10^{7} |

Method | D_{3/2}-S_{1/2} | D_{5/2}-S_{1/2} |
---|---|---|

HF | 1360 | 1360 |

DF | 257.0 | 77.4 |

DF (_{exp}) | 63.9 | 13.3 |

RPT-MP | 54.54 | 11.8 |

RPT-EA | 54.52 (0.2%) | 11.7 (0.2%) |

Exp. | 53.5 ± 2.0 | 11.6 ± 0.4 |

The standard HF and DF approaches in the single-configuration approximations do not allow to obtain very accurate results. Using the empirical transition energies significantly improve the theoretical results as in fact it means an account for very important interparticle correlations effects. In our approach, the corresponding exchange-correlation effects (the polarization interaction of the QPs, mutual screening and anti-screening corrections, etc.) are taken into account more accurately. The core polarization correction to the transition probability is of great importance as it changes significantly the probability value (~15–40%). It should be also noted that the gauge-noninvariant contribution to radiation width is very small (0.2%; see Table 2 in the line “EA”) that means equivalence of the calculation results in the standard amplitude approach with using the length and velocity forms for transition operator. From the other side, this is an evidence of the successful choice of the PT zeroth approximation and accurate account of the multi-particle correlation effects.

We have presented the fundamentals of the new relativistic many-body PT formalism with construction of the optimized one-QP representation in the theory of relativistic multielectron systems. The relativistic density-functional approximation with the Kohn-Sham potential is taken as the zeroth one and all exchange-correlation corrections of the second-order and dominated classes of the higher-orders diagrams (polarization interaction, QPs screening, etc.) have been taken into account. In order to reach the corresponding optimization, we have used a procedure of the accurate treating of the PT lowest-order multielectron effects, in particular, the gauge-dependent radiative contribution for the certain class of the photon propagator gauge. The corresponding contribution is considered to be the typical representative of the important multielectron exchange-correlation effects, whose minimization is reasonable criteria in the searching for the optimal PT one-electron basis. This procedure derives an undoubted profit in the routine many-body calculations as it provides the way of the refinement of the atomic (molecular) characteristics calculations, based on the “first principles.” The presented relativistic PT formalism can be further generalized, in particular, by the way of accounting for the radiation, QED (the Lamb shift self-energy and vacuum polarization corrections, for instance in the effective Uhling-Serber approximation with account for the Källen-Sabry and Wichmann-Kroll corrections), and nuclear (the Bohr-Weisskopf and Breit-Rosenthal-Crawford-Schawlow effects, nuclear finite size correction, magnetic moment distribution, etc.) effects [13–23].

## References

- 1.
Wilson S. Handbook on Molecular Physics and Quantum Chemistry. Chichester: Wiley; 2003. p. 680 - 2.
Dyall KG, Faegri R, Jr. Introduction to Relativistic Quantum Theory. Oxford: Oxford University Press; 2007. p. 780 - 3.
Bartlett RJ, Musiał MN. Coupled-cluster theory in quantum chemistry. Review of Modern Physics. 2007; 79 :291–328 - 4.
Shavitt I, Bartlett RJ. Many-Body Methods in Chemistry and Physics: MBPT and Coupled-Cluster Theory. New York, NY: Cambridge University Press; 2009. p. 390 - 5.
Fritzsche S, Fricke B, Sepp W-D. Computer-algebraic derivation of atomic Feynman–Goldstone expansions. In: Schweizer W, Schmelcher P, editors. Atoms and Molecules in Strong External Fields. Vol. 1. New York, NY: Kluwer; 1997. pp. 101–108 - 6.
Parpia FA, Froese-Fischer C, Grant IP. Generalized relativistic atomic structure package: GRASP. Computer Physics Communications. 1996; 94 (2):249–270 - 7.
Glushkov AV. Relativistic Quantum Theory. Relativistic Quantum Mechanics of Atomic Systems. Odessa: Astroprint; 2008. p. 700 - 8.
Kohn W, Sham S. Quantum density oscillations in an inhomogeneous electron gas. Physical Review A. 1965; 137 :1697–1710 - 9.
Johnson WR, Sapistein J, Blundell SA. Atomic structure calculations associated with PNC experiments in atomic Cs. Physica Scripta T. 1993; 46 :184–192 - 10.
Desclaux JP. A multiconfiguration relativistic DIRAC-FOCK program. Computer Physics Communications. 1975; 9 (1):31–45 - 11.
Feller D, Davidson ER. An approximation to frozen natural orbitals through use of Hartree-Fock exchange potential. Journal of Chemical Physics. 1981; 74 (7):3977–3989 - 12.
Dietz K, Hess BA. Single particle orbitals for configuration interaction derived from quantum electrodynamics. Physica Scripta. 1989; 39 :682–688 - 13.
Glushkov AV, Khetselius OYu, Malinovskaya SV. New laser-electron nuclear effects in the nuclear γ transition spectra in atomic and molecular systems. In: Wilson S, Grout PJ, Maruani J, Delgado-Barrio G, Piecuch P, editors. Frontiers in Quantum Systems in Chemistry and Physics. Series: Progress in Theoretical Chemistry and Physics. Vol. 18. Berlin: Springer; 2008. pp. 523–542 - 14.
Glushkov AV, Ambrosov SV, Loboda AV, Gurnitskaya EP, Khetselius OYu. QED calculation of heavy multicharged ions with account for the correlation, radiative and nuclear effects. In: Julien J-P, Maruani J, Mayou D, Wilson S, Delgado-Barrio G, editors. Recent Advances in Theoretical Physics and Chemistry Systems. Series: Progress in Theoretical Chemistry and Physics. Vol. 15. Berlin: Springer; 2006. pp. 285–300 - 15.
Glushkov AV, Khetselius OYu, Gurnitskaya EP, Loboda AV, Florko TA, Sukharev DE, Lovett L. Gauge-invariant QED perturbation theory approach to calculating nuclear electric quadrupole moments, hyperfine structure constants for heavy atoms and ions. In: Wilson S, Grout PJ, Maruani J, Delgado-Barrio G, Piecuch P, editors. Frontiers in Quantum Systems in Chemistry and Physics. Series: Progress in Theoretical Chemistry and Physics. Vol. 18. Berlin: Springer; 2008. pp. 501–522 - 16.
Glushkov AV, Khetselius OYu, Lovett L. Electron-β-nuclear spectroscopy of atoms and molecules and chemical environment effect on the β-decay parameters. In: Piecuch P, Maruani J, Delgado-Barrio G, Wilson S, editors. In: Advances in the Theory of Atomic and Molecular Systems Dynamics, Spectroscopy, Clusters, and Nanostructures. Series: Progress in Theoretical Chemistry and Physics. Vol. 20. Berlin: Springer; 2010. pp. 125–152 - 17.
Glushkov AV, Khetselius OYu, Svinarenko AA. Relativistic theory of cooperative muon-gamma-nuclear processes: Negative muon capture and metastable nucleus discharge. In: Hoggan P, Brandas E, Maruani J, Delgado-Barrio G, Piecuch P, editors. Advances in the Theory of Quantum Systems in Chemistry and Physics. Series: Progress in Theoretical Chemistry and Physics. Vol. 22. Berlin: Springer; 2011. pp. 51–70 - 18.
Glushkov AV. Advanced relativistic energy approach to radiative decay processes in multielectron atoms and multicharged ions. In: Nishikawa K, Maruani J, Brandas E, Delgado-Barrio G, Piecuch P, editors. Quantum Systems in Chemistry and Physics: Progress in Methods and Applications. Series: Progress in Theoretical Chemistry and Physics. Vol. 26. Berlin: Springer; 2012. pp. 231–254 - 19.
Khetselius OYu. Relativistic energy approach to cooperative electron-gamma-nuclear processes: NEET Effect. In: Nishikawa K, Maruani J, Brandas E, Delgado-Barrio G, Piecuch P, editors. Quantum Systems in Chemistry and Physics: Progress in Methods and Applications. Series: Progress in Theoretical Chemistry and Physics. Vol. 26. Berlin: Springer; 2012. pp. 217–230 - 20.
Glushkov AV, Khetselius OYu, Svinarenko AA, Prepelitsa GP. Energy approach to atoms in a laser field and quantum dynamics with laser pulses of different shape. In: Duarte FJ, editor. Coherence and Ultrashort Pulsed Emission Rijeka: InTech; 2010. pp. 159–186 - 21.
Glushkov AV. Operator perturbation theory for atomic systems in a strong DC electric field. In: Hotokka M, Maruani J, Brandas E, Delgado-Barrio G, editors. Advances in Quantum Methods and Applications in Chemistry, Physics, and Biology. Series: Progress in Theoretical Chemistry and Physics. Vol. 27. Berlin: Springer; 2013. pp. 161–178 - 22.
Khetselius OYu. Optimized perturbation theory to calculating the hyperfine line shift and broadening for heavy atoms in the buffer gas. In: Nascimento M, Maruani J, Brandas E, Delgado-Barrio G, Piecuch P, editors. Frontiers in Quantum Methods and Applications in Chemistry and Physics. Series: Progress in Theoretical Chemistry and Physics. Vol. 29. Berlin: Springer; 2015. pp. 54–76 - 23.
Glushkov AV, Svinarenko AA, Khetselius OYu, Buyadzhi VV, Florko TA, Shakhman AN. Relativistic quantum chemistry: Advanced approach to construction of the Green's function of the Dirac equation with complex energy and mean-field nuclear potential. In: Nascimento M, Maruani J, Brandas E, Delgado-Barrio G, Piecuch P, editors. Frontiers in Quantum Methods and Applications in Chemistry and Physics. Series: Progress in Theoretical Chemistry and Physics. Vol. 29. Berlin: Springer; 2015. pp. 197–218 - 24.
Safronova UI, Johnson W, Derevianko A. Relativistic many-body calculations of energy levels, hyperfine constants, electric-dipole matrix elements, static polarizabilities for alkali-metal atoms. Physical Review A. 1999; 60 :4476–4486 - 25.
Safronova MS, Safronova UI. Critically evaluated theoretical energies, lifetimes, hyperfine constants, and multipole polarizabilities in 87Rb. Physical Review A. 2011; 83 :052508 - 26.
Safronova UI, Safronova MS. Third-order relativistic many-body calculations of energies, transition rates, hyperfine constants, and blackbody radiation shift in 171Yb+. Physical Review A. 2009; 79 :022512 - 27.
Safronova MS, Johnson WR, Safronova UI. Excitation energies, hyperfine constants, E1,E2,M1 transition rates and lifetimes of 6s2nl states in TlI and PbII . Physical Review A. 2005; 71 :052506 - 28.
Yerokhin V, Artemyev A, Shabaev VM. QED treatment of electron correlation in Li-like ions. Physical Review A. 2007; 75 :062501 - 29.
Dzuba VA, Flambaum VV, Silvestrov PG, Sushkov DE. Many-body perturbation theory calculations in atoms with open shells. Physical Reviews A. 1991; 44 :2828–2831 - 30.
Dzuba VA, Harabati C, Johnson WR, Safronova MS. Breit correction to the parity-nonconservation amplitude in cesium. Physical Reviews A. 2001; 63 :044103 - 31.
Ostrovsky VN, Sheynerman SA. Radiation transitions in external shells of the ion Hg+. Optics and Spectroscopy. 1989; 67 :16–22 - 32.
Martin W. NIST Spectra Database, Version 2.0 (http://physics.nist.gov.asd). Washington: NIST; 2004 - 33.
Ivanov LN, Ivanova EP. Extrapolation of atomic ion energies by model potential method: Na-like spectra. Atomic Data and Nuclear Data Tables. 1979; 24 (2):95–121 - 34.
Driker MN, Ivanova EP, Ivanov LN, Shestakov AF. Relativistic calculation of spectra of 2-2 transitions in O-and F-like atomic ions. Journal of Quantitative Spectroscopy and Radiative Transfer. 1982; 28 (6):531–535 - 35.
Ivanov LN, Letokhov VS. Spectroscopy of autoionization resonances in heavy elements atoms. Com.Modern Phys. D.: Atomic and Molecular Physics. 1985; 4 :169–184 - 36.
Ivanova EP, Ivanov LN, Glushkov AV, Kramida AE. High-order corrections in relativistic perturbation theory with model zero approximation. Physica Scripta. 1985; 32 (4):512–524 - 37.
Ivanova EP, Glushkov AV. Theoretical study of multicharged ions spectra of F-, Ne-isoelectronic sequences. Journal of Quantitative Spectroscopy & Radiative Transfer. 1986; 36 (N2):127–145 - 38.
Bekov GI, Vidolova-Angelova EP, Ivanov LN, Letokhov VS, Mishin VI. Laser spectroscopy of narrow doubly excited autoionizing states of ytterbium atoms. Soviet Physics JETP. 1981; 53 (3):441–447 - 39.
Vidolova-Angelova E, Ivanov LN, Ivanova EP, Angelov DA. Relativistic perturbation method for studying radiation decay of highly-excited many-electron atoms: Tm. Journal of Physics B: Atomic, Molecular and Optical Physics. 1986; 19 :2053–2069 - 40.
Vidolova-Angelova E, Ivanov LN, Angelov DA. Autoionization decay of excited Rydberg Tm states. Journal of Physics B: Atomic, Molecular and Optical Physics. 1988; 21 :3877–3890 - 41.
Ivanova EP, Ivanov LN, Aglitsky EV. Modern trends in spectroscopy of multi-charged ions. Physics Reports. 1988; 166 (N6):315–390 - 42.
Vidolova-Angelova E, Ivanov LN. Autoionizing Rydberg states of thulium. Re-orientation decay due to monopole interaction. Journal of Physics B: Atomic, Molecular and Optical Physics. 1991; 24 :4147–4158 - 43.
Ivanova EP, Gulov AV. Theoretical investigation of the neon isoelectronic sequence. Atomic Data And Nuclear Data Tables. 1991; 49 :1–64 - 44.
Ivanov LN, Ivanova EP, Knight L. Energy approach to consistent QED theory for calculation of electron-collision strengths. Physics Review A. 1993; 48 (9):4365–4374 - 45.
Ivanov LN, Ivanova EP. Sturm orbitals method in calculation of physical characteristics of radiation of the atoms. JETP. 1996; 110 (6):483–498 - 46.
Ivanova EP, Zinoviev NA. The possibility of X-ray lasers based on inner-shell transitions of Ne-like ions. Physics Letters A. 2001; 274 :239–246 - 47.
Glushkov AV, Ivanov LN, Ivanova, EP. Relativistic decay of excited atomic states. Generalized energy approach. In: Autoionization Phenomena in Atoms. Moscow: Moscow University Press; 1986. pp. 152–164 - 48.
Glushkov AV, Ivanov LN. Radiation decay of atomic states: Atomic residue and gauge non-invariant contributions. Physics Letters A. 1992; 170 (1):33–37 - 49.
Glushkov AV. Negative ions of inert gases. JETP Letters. 1992; 55 (2):97–100 - 50.
Glushkov AV. Relativistic polarization potential of a many-electron atom. Soviet Physics Journal. 1990; 33 (1):1–8 - 51.
Glushkov AV. Correction for exchange and correlation effects in multielectron system theory. Journal of Structural Chemistry. 1990; 31 (4):529–533 - 52.
Glushkov AV. Relativistic multiconfiguration time-dependent self-consistent-field theory for molecules. Soviet Physics Journal. 1991; 34 (10):871–876 - 53.
Glushkov AV, Ivanov LN, Letokhov VS. Nuclear Quantum Optics. Moscow-Troitsk: Preprint of Institute for Spectroscopy of the USSR Academy of Sciences (ISAN); 1991. AS-N4 :1–16 - 54.
Glushkov AV, Ivanov LN. DC strong field Stark effect: Consistent quantum mechanical approach. Moscow-Troitsk: Preprint of Institute for Spectroscopy of the USSR Academy of Sciences (ISAN); 1992. AS-N1 :1–16 - 55.
Glushkov AV, Ivanov LN. A broadening of the thulium atom autoionization resonances in a weak electric field. Moscow-Troitsk: Preprint of Institute for Spectroscopy of the USSR Academy of Sciences (ISAN); 1992. AS-N2 :1–10 - 56.
Glushkov AV, Ivanov LN. Shift and deformation of radiation atomic lines in the laser emission field. Multiphoton processes. Moscow-Troitsk: Preprint of Institute for Spectroscopy of the USSR Academy of Sciences (ISAN); 1992. AS-N3 :1–12 - 57.
Glushkov AV, Ivanov LN. DC strong-field Stark-effect: Consistent quantum-mechanical approach. Journal of Physics B: Atomic, Molecular and Optical Physics. 1993; 26 (16):L379-L396 - 58.
Glushkov AV, Ambrosov SV, Ignatenko AV, Korchevsky DA. DC strong field Stark effect for non-hydrogenic atoms: New consistent quantum mechanical approach. International Journal of Quantum Chemistry. 2004; 99 (5):936–939 - 59.
Glushkov AV, Malinovskaya SV, Chernyakova YuG, Svinarenko AA. Cooperative laser-electron-nuclear processes: QED calculation of electron satellites spectra for multi-charged ion in laser field. International Journal of Quantum Chemistry. 2004; 99 (5):889–893 - 60.
Glushkov AV, Ambrosov SV, Loboda AV, Chernyakova YuG, Svinarenko AV, Khetselius OYu. QED calculation of the superheavy elements ions: Energy levels, radiative corrections and hyperfine structure for different nuclear models. Nuclear Physics A. 2004; 734S :21–25 - 61.
Glushkov AV, Khetselius OYu, Malinovskaya SV. Optics and spectroscopy of cooperative laser-electron nuclear processes in atomic and molecular systems—New trend in quantum optics. European Physics Journal. 2008; T160 :195–204 - 62.
Glushkov AV, Ambrosov SV, Loboda AV, Gurnitskaya EP, Prepelitsa GP. Consistent QED approach to calculation of electron-collision excitation cross-sections and strengths: Ne-like ions. International Journal of Quantum Chemistry. 2005; 104 (1):562–569 - 63.
Glushkov AV, Malinovskaya SV, Loboda AV, Gurnitskaya EP, Korchevsky DA. Diagnostics of the collisionally pumped plasma and search of the optimal plasma parameters of x-ray lasing: Calculation of electron-collision strengths and rate coefficients for Ne-like plasma. Journal of Physics: Conference Series. 2005; 178 :188–198 - 64.
Glushkov AV, Loboda AV, Gurnitskaya EP, Svinarenko AA. QED theory of radiation emission and absorption lines for atoms in a strong laser field. Physica Scripta. 2009; T.135 :014022 - 65.
Malinovskaya SV, Glushkov AV, Khetselius OY, Svinarenko AA, Mischenko EV, Florko TA. Optimized perturbation theory scheme for calculating the interatomic potentials and hyperfine lines shift for heavy atoms in the buffer inert gas. International Journal of Quantum Chemistry. 2009; 109 (14):3325–3329 - 66.
Malinovskaya SV, Glushkov AV, Khetselius OY, Lopatkin YM, Loboda AV, Nikola LV, Svinarenko AA, Perelygina TB. Generalized energy approach for calculating electron collision cross-sections for multicharged ions in a plasma: Debye shielding model. International Journal of Quantum Chemistry. 2011; 111 (2):288–296 - 67.
Glushkov AV, Malinovskaya SV, Khetselius OYu, Loboda AV, Sukharev DE, Lovett L. Green's function method in quantum chemistry: New numerical algorithm for the Dirac equation with complex energy and Fermi-model nuclear potential . International Journal of Quantum Chemistry. 2009; 109 (8):1717–1727