**Relativistic Perturbation Theory Formalism to Computing Spectra and Radiation Characteristics: Application to Heavy Elements**

Alexander V. Glushkov, Olga Yu. Khetselius, Andrey A. Svinarenko, Vasily V. Buyadzhi, Valentin B. Ternovsky, Anna A. Kuznetsova and Petr G. Bashkarev

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/intechopen.69102

#### 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 exchangecorrelation 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

© 2017 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## 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 exchangecorrelation 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]).

1. Introduction

132 Recent Studies in Perturbation Theory

atomic, nuclear, and molecular calculations [1–7].

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)

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 exchangecorrelation 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.

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 exchangecorrelation 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 manybody 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 (Tolmachev-Ivanov-Ivanova, 1969–1974). The key algorithm of this approach includes construction of the secular matrix M 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:

$$
\Delta E = \text{Re}\,\Delta E + i \,\text{Im}\,\Delta E \quad \text{Im}\,\Delta E = \,-\,\,\Gamma/2. \tag{1}
$$

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 M. In Refs. [33–37], the ReΔE calculation procedure has been generalized for the case of nearly degenerate states, whose levels form a more or less compact group. Naturally, the matrix M reduces to one term (ΔE) in the case of well-identified and separated energy spectrum. The Gell-Mann and Low formula allow further to obtain the expansion of the M elements into PT series on interelectron interaction and apply the standard Feynman diagrammatic technique. The corresponding PT series is as follows:

$$M = \ \ M^{(0)} + \mathcal{M}^{(1)} + \mathcal{M}^{(2)} + \mathcal{M}^{(3)}.\tag{2}$$

Here, M(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); M(1), M(2),and M(3) are the contributions, which correspond to the one-, two- and three-QP PT diagrams, respectively. The diagonal matrix M(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 M(2) and M(3)contributions. Using the Feynman diagrams technique, the authors [33–38] have in detail analyzed the M(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

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

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 invari-

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 (Tolmachev-Ivanov-Ivanova, 1969–1974). The key algorithm of this approach includes construction of the secular matrix M 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

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 M. In Refs. [33–37], the ReΔE calculation procedure has been generalized for the case of nearly degenerate states, whose levels form a more or less compact group. Naturally, the matrix M reduces to one term (ΔE) in the case of well-identified and separated energy spectrum. The Gell-Mann and Low formula allow further to obtain the expansion of the M elements into PT series on interelectron interaction and apply the standard

<sup>M</sup> <sup>¼</sup> <sup>M</sup>ð0<sup>Þ</sup> <sup>þ</sup> <sup>M</sup>ð1<sup>Þ</sup> <sup>þ</sup> <sup>M</sup>ð2<sup>Þ</sup> <sup>þ</sup> <sup>M</sup>ð3<sup>Þ</sup>

Here, M(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); M(1), M(2),and M(3) are the

Feynman diagrammatic technique. The corresponding PT series is as follows:

ΔΕ ¼ Re ΔE þ i ImΔE ImΔE ¼ � Γ=2: ð1Þ

: ð2Þ

2. Relativistic many-body perturbation theory with optimized

one-quasiparticle zeroth representation

(ion of Hg+).

2.1. General remarks

134 Recent Studies in Perturbation Theory

ance principle fulfilling.

shift can be defined as follows:

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

$$H = \sum\_{i} \{acp\_i - \beta c^2 - Z/r\_i\} + \sum\_{i>j} \exp(i|\omega|r\_{ij})(1 - \alpha\_i \alpha\_j)/r\_{ij\star} \tag{3}$$

where Z is a charge of nucleus, αi, α<sup>j</sup> are the Dirac matrices, ωij is the transition frequency, and c, 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 α<sup>2</sup> (α 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 ρ(r) 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|R) as [14]

$$V\_{mcl}(r|\mathbb{R}) = -\left( (1/r) \int\_0^r dr' r'^2 \rho(r'|\mathbb{R}) + \int\_r^\infty dr' r' \rho(r'|\mathbb{R}).\tag{4}$$

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<sup>0</sup> and perturbation operator can be presented in the standard form as follows [7, 14, 15]:

$$\begin{aligned} H\_0 &= \sum\_i a\_i^+ a\_i E\_i \\ H\_{\text{int}} &= \sum\_{\vec{i}\vec{j}} a\_i^+ a\_{\vec{i}} V\_{\vec{i}\vec{j}} + \frac{1}{2} \sum\_{\vec{i}\vec{j}\vec{k}} V\_{\vec{i}\vec{k}l} a\_i^+ a\_{\vec{j}}^+ a\_k a\_l \\ V\_{\vec{i}\vec{j}} &= \left[ \int \vec{r} \cdot \boldsymbol{\wp}\_i(\vec{r}) [-V\_{\text{MF}}(r)] \cdot \boldsymbol{\wp}(\vec{r}) \right. \\ V\_{\vec{i}\vec{j}\vec{k}} &= \iint d\vec{r}\_1 d\vec{r}\_2 \, \boldsymbol{\wp}(\vec{r}\_1) \, \boldsymbol{\wp}(\vec{r}\_2) \, V(r\_1 r\_2) \, \boldsymbol{\wp}\_k(\vec{r}\_2) \, \boldsymbol{\wp}\_l(\vec{r}\_1), \end{aligned} \tag{5}$$

where ϕðr !Þ are one-electron functions (Dirac bispinors), Ei , one-electron energies, and <sup>V</sup>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 VMF. In our PT version, we use the potential

$$V\_{MF} = V^{DKS}(r) = \left[V\_{Coul}^{D}(r) + V\_X(r) + V\_{\mathcal{C}}(r|a)\right] \tag{6}$$

Further as VX(r) we use the standard Kohn-Sham (KS) exchange potential as follows [8]:

$$V\_X^{KS}(r) = - (\mathbf{1}/\pi) \left[ \mathfrak{J} \pi^2 \rho(r) \right]^{1/3}. \tag{7}$$

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

$$V\_X[\rho(r), r] = \frac{\delta E\_X[\rho(r)]}{\delta \rho(r)},\tag{8}$$

In the relativistic multielectron theory with a Hamiltonian having a transverse vector potential (for describing the photons), one could determine the homogeneous density ρ(r), construct the corresponding exchange energy EX[ρ(r)], and introduce the following exchange potential [16]:

$$V\_X[\rho(r), r] = V\_X^{KS}(r) \cdot \left\{ \frac{3}{2} \ln \frac{[\beta + (\beta^2 + 1)^{1/2}]}{\beta(\beta^2 + 1)^{1/2}} - \frac{1}{2} \right\},\tag{9}$$

where <sup>β</sup> ¼ ½3π<sup>2</sup>ρðrÞ�<sup>1</sup>=<sup>3</sup> =c. The corresponding correlation functional is as follows [16, 17]:

$$V\_{\mathbb{C}}[\rho(r), r] = -0.0333 \cdot b \cdot \ln[1 + 18.3768 \cdot \rho(r)^{1/3}],\tag{10}$$

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 F and G components can be written as [14] follows:

$$\begin{aligned} f' &= - (\chi + |\chi|) f/r - a Z V \mathbf{g} - (a Z E\_{n\chi} + \mathbf{2}/a Z) \mathbf{g}, \\ \mathbf{g'} &= (\chi - |\chi|) \mathbf{g}/r - a Z V f + a Z E\_{n\chi} f. \end{aligned} \tag{11}$$

Here, En<sup>χ</sup> is one-electron energy without the rest energy. The boundary values are defined by the first terms of the Taylor expansion:

$$g = \left( V(0) - E\_{\pi\chi} \right) r\alpha Z / (2\chi + 1); \quad f = 1 \text{ at } \chi < 0,\tag{12a}$$

$$f = \left(V(0) - E\_{n\chi} - 2/\alpha^2 Z^2\right) \alpha Z; \quad g = 1 \text{ at } \chi > 0. \tag{12b}$$

The condition f, g!0 at r!∞ determines the quantified energies of the state Enχ. 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]:

$$\begin{aligned} M\_1^{(2)} &= \{n\_1 l\_1 j\_1 \quad n\_2 l\_2 [j] \, | \, V\_{\text{int}} \, | \, n\_4 l\_4 \dot{j}\_4 \quad n\_3 l\_3 \dot{j}\_3 [l] \} \\ &= P\_1 P\_2 (-1)^{1 + j\_2 + j\_4 + l} [(2j\_1 + 1)(2j\_2 + 1)(2j\_3 + 1)(2j\_4 + 1)]^{1/2} \\ &\quad \times \sum\_{i,k} \sum\_a \left\{ \begin{matrix} j\_i \dot{j}\_k \\ j\_2 j\_1 a \end{matrix} \right\} \left( \delta\_{i,3} \delta\_{k,4} + (-1)^l \delta\_{i,4} \delta\_{k,3} \right) \cdot Q\_{\lambda \nu} \end{aligned} \tag{13}$$

where

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 VMF. In our PT

Further as VX(r) we use the standard Kohn-Sham (KS) exchange potential as follows [8]:

<sup>X</sup> <sup>ð</sup>rÞ ¼ �ð1=πÞ½3π<sup>2</sup>

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

VX½ρðrÞ, r� ¼ <sup>δ</sup>EX½ρðrÞ�

In the relativistic multielectron theory with a Hamiltonian having a transverse vector potential (for describing the photons), one could determine the homogeneous density ρ(r), construct the corresponding exchange energy EX[ρ(r)], and introduce the following exchange potential [16]:

> <sup>X</sup> <sup>ð</sup>rÞ � <sup>3</sup> 2

VC½ρðrÞ, r�¼�0:0333 � b � ln½1 þ 18:3768 � ρðrÞ

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 F and G components can be written as [14] follows:

<sup>0</sup> ¼ �ðχ þ jχjÞf =r � αZVg � ðαZEn<sup>χ</sup> þ 2=αZÞg,

Here, En<sup>χ</sup> is one-electron energy without the rest energy. The boundary values are defined by

Z2 �

The condition f, g!0 at r!∞ determines the quantified energies of the state Enχ. The system of Eq. (11) is numerically solved by the Runge-Kutta method ('Superatom" package is used [7, 13–

<sup>ρ</sup>ðrÞ�<sup>1</sup>=<sup>3</sup>

ln <sup>½</sup><sup>β</sup> þ ðβ<sup>2</sup> <sup>þ</sup> <sup>1</sup><sup>Þ</sup>

<sup>β</sup>ðβ<sup>2</sup> <sup>þ</sup> <sup>1</sup><sup>Þ</sup>

=c. The corresponding correlation functional is as follows [16, 17]:

<sup>g</sup><sup>0</sup> ¼ ð<sup>χ</sup> � jχjÞg=<sup>r</sup> � <sup>α</sup>ZVf <sup>þ</sup> <sup>α</sup>ZEnχ<sup>f</sup> : <sup>ð</sup>11<sup>Þ</sup>

( )

1=2 �

<sup>1</sup>=<sup>2</sup> � <sup>1</sup> 2

1=3

rαZ=ð2χ þ 1Þ; f ¼ 1 at χ < 0, ð12aÞ

αZ; g ¼ 1 at χ > 0: ð12bÞ

CoulðrÞ þ VXðrÞ þ VCðrjaÞ� ð6Þ

: ð7Þ

; ð9Þ

�, ð10Þ

δρðr<sup>Þ</sup> , <sup>ð</sup>8<sup>Þ</sup>

VMF <sup>¼</sup> <sup>V</sup>DKSðrÞ¼½V<sup>D</sup>

VKS

VX½ρðrÞ, r� ¼ <sup>V</sup>KS

f

g ¼ �

f ¼ �

Vð0Þ � En<sup>χ</sup>

�

<sup>V</sup>ð0Þ � En<sup>χ</sup> � <sup>2</sup>=α<sup>2</sup>

the first terms of the Taylor expansion:

version, we use the potential

136 Recent Studies in Perturbation Theory

where <sup>β</sup> ¼ ½3π<sup>2</sup>ρðrÞ�<sup>1</sup>=<sup>3</sup>

23, 34, 36, 47–67]).

$$P\_1 = \begin{cases} 1 & \text{if } \quad n\_1 l\_1 j\_1 \neq n\_2 l\_2 j\_2 \\ \nearrow 2 & \text{if } \quad n\_1 l\_1 j\_1 = n\_2 l\_2 j\_2 \end{cases} \\ P\_2 = \begin{cases} 1 & \text{if } \quad n\_3 l\_3 j\_3 \neq n\_4 l\_4 j\_4 \\ \nmid & \quad n\_3 l\_3 j\_3 = n\_4 l\_4 j\_4 \end{cases} \tag{14}$$

The value of the Q<sup>λ</sup> can be expressed through the radial Slater-like integrals and presented as a sum of the Coulomb and Breit parts: <sup>Q</sup><sup>λ</sup> <sup>¼</sup> <sup>Q</sup>Qul <sup>λ</sup> <sup>þ</sup> <sup>Q</sup>Br <sup>λ</sup> , which corresponds to a partition of the interelectron potential into the Coulomb and Breit ones in the second term of Eq. (1). Let us remind that, for instance, the Coulomb part in Eq. (13) is expressed through the radial integrals and angle coefficients as follows:

$$\begin{split} \mathbf{Q}\_{\boldsymbol{\Lambda}}^{\text{Qui}} &= \frac{1}{Z} \left\{ \mathbf{R}\_{l} (1243) \mathbf{S}\_{\boldsymbol{\Lambda}} (1243) + \mathbf{R}\_{l} (\mathbf{\tilde{1}24\tilde{3}}) \mathbf{S}\_{\boldsymbol{\Lambda}} (\mathbf{\tilde{1}24\tilde{3}}) + \\ &+ \mathbf{R} (1\mathbf{\tilde{2}4}\mathbf{\tilde{3}}) \mathbf{S}\_{\boldsymbol{\Lambda}} (12\tilde{4}\tilde{4}3) + \mathbf{R} (\mathbf{\tilde{1}24\tilde{3}}) \mathbf{S}\_{\boldsymbol{\Lambda}} (\mathbf{\tilde{1}24\tilde{3}}) \right\}. \end{split} \tag{15}$$

In the nonrelativistic limit, there remains only the first term in Eq. (15) depending only on the large component f(r) of the one-electron Dirac functions. For example, its imaginary part is as follows [36]:

$$\begin{aligned} \text{Im } R\_{\lambda}(12; 43) &= \frac{1}{2} (2\lambda + 1) \pi X\_{\lambda}(13) X\_{\lambda}(24) \\ X\_{\lambda}(12) &= \int dr \, r^{3/2} f\_1(r) f\_{\lambda + 1/2}^{(1)} \left( r \alpha Z |\omega| f\_2(r) \right) \end{aligned} \tag{16}$$

The angular coefficient has only a real part:

$$\mathcal{S}\_{\Lambda}(12;43) = \mathcal{S}\_{\Lambda}(13)\mathcal{S}\_{\Lambda}(24) \qquad \qquad \mathcal{S}\_{\Lambda}(13) = \begin{Bmatrix} \lambda l\_1 l\_3 \end{Bmatrix} \begin{pmatrix} j\_1 & j\_3 & \lambda \\ 1 & -\frac{1}{2} & 0 \end{pmatrix} \tag{17}$$

Here, {λl1l3} means that λ, l<sup>1</sup> and l<sup>3</sup> must satisfy the triangle rule and the sum λ + l<sup>1</sup> + l<sup>3</sup> 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 f must be replaced by the small one <sup>g</sup>, and instead of li, <sup>~</sup>li <sup>¼</sup> li � 1 should be taken for ji <sup>&</sup>lt; li and <sup>~</sup>li <sup>¼</sup> li <sup>þ</sup> 1 for ji <sup>&</sup>gt; li. 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]:

$$\begin{aligned} \mathcal{W}\_{\text{pol}}^{\text{f}}(r\_{1}r\_{2}) &= \\ \displaystyle X \left\{ \int \frac{dr' \left(\rho\_{\varepsilon}^{(0)}(r')\right)^{1/3} \Theta(r')}{|r\_{1} - r'| \cdot |r' - r\_{2}|} - \int \frac{dr' \left(\rho\_{\varepsilon}^{(0)}(r')\right)^{1/3} \Theta(r')}{|r\_{1} - r'|} \right\} \frac{dr'' \left(\rho\_{\varepsilon}^{(0)}(r'')\right)^{1/3} \Theta(r'')}{|r'' - r\_{2}|} / \left\langle \left(\rho\_{\varepsilon}^{(0)}\right)^{1/3} \right\rangle \Bigg\} \end{aligned} \tag{18a}$$

$$
\left\langle \left( \rho\_c^{(0)} \right)^{1/3} \right\rangle = \int dr \left( \rho\_c^{(0)}(r) \right)^{1/3} \theta(r) \tag{18b}
$$

$$\theta(r) = \left\{ 1 + [3\pi^2 \cdot \rho\_c^{(0)}(r)]^{2/3} / c^2 \right\}^{1/2} \tag{18c}$$

Figure 1. Some diagrams of the second order, taking into account the exchange and polarization interaction of quasiparticles and electrons of the closed shells core.

Figure 2. Some diagrams of the second order, describing a direct interaction of the two or three external quasiparticles.

where ρ<sup>0</sup> <sup>c</sup> is the core electron density (without account for the quasiparticle), X is the numerical coefficient, and c 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].

In order to take into account the ladder diagrams contributions as well as some of the threequasiparticle 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 <sup>M</sup><sup>~</sup> <sup>ð</sup>1<sup>Þ</sup> and further 2-QP one too; <sup>M</sup><sup>~</sup> <sup>ð</sup>2<sup>Þ</sup> is computed using the PT first-order formulae and the modified radial 1-QP wave functions.

## 2.4. Optimization of the relativistic orbitals basis

small one <sup>g</sup>, and instead of li, <sup>~</sup>li <sup>¼</sup> li � 1 should be taken for ji <sup>&</sup>lt; li and <sup>~</sup>li <sup>¼</sup> li <sup>þ</sup> 1 for ji <sup>&</sup>gt; li. The

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

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

�<sup>1</sup>=<sup>3</sup> � � <sup>8</sup>

ð dr<sup>00</sup> � ρð0<sup>Þ</sup> <sup>c</sup> ðr00Þ

<sup>c</sup> <sup>ð</sup>rÞ�<sup>2</sup>=<sup>3</sup>

<sup>2</sup> n o<sup>1</sup>=<sup>2</sup>

Figure 1. Some diagrams of the second order, taking into account the exchange and polarization interaction of quasipar-

Figure 2. Some diagrams of the second order, describing a direct interaction of the two or three external quasiparticles.

=c

�<sup>1</sup>=<sup>3</sup>

jr<sup>00</sup> � r2j

θðr00Þ

. � ρð0<sup>Þ</sup> c

θðrÞ, ð18bÞ

9 >=

>;

ð18aÞ

ð18cÞ

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

external quasiparticles (or antiscreening in the case, say, of an electron and a vacancy).

jr<sup>1</sup> � r<sup>0</sup> j

> ¼ ð dr � ρð0<sup>Þ</sup> <sup>c</sup> ðrÞ �<sup>1</sup>=<sup>3</sup>

<sup>θ</sup>ðrÞ ¼ <sup>1</sup> þ ½3π<sup>2</sup> � <sup>ρ</sup>ð0<sup>Þ</sup>

operator can be taken in the following form [50]:

j�jr<sup>0</sup> � <sup>r</sup>2<sup>j</sup> �

ð dr<sup>0</sup> � ρð0<sup>Þ</sup> <sup>c</sup> ðr<sup>0</sup> Þ �<sup>1</sup>=<sup>3</sup> θðr<sup>0</sup> Þ

� ρð0<sup>Þ</sup> c �<sup>1</sup>=<sup>3</sup> � �

Vd

X

polðr1r2Þ ¼

138 Recent Studies in Perturbation Theory

ð dr<sup>0</sup> � ρð0<sup>Þ</sup> <sup>c</sup> ðr<sup>0</sup> Þ �<sup>1</sup>=<sup>3</sup> θðr<sup>0</sup> Þ

><

>:

jr<sup>1</sup> � r<sup>0</sup>

ticles and electrons of the closed shells core.

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 gaugedependent multielectron contributions ImΔEninv 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ΔE (the radiation decay width).

In the fourth order of QED PT (the second order of the atomic PT), the diagrams appear, whose contribution to the ImΔEninv accounts for the multielectron exchange-correlation (polarization) effects (diagrams Ad, Aex; 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

Figure 3. B: second other PT diagram contributing the imaginary energy part related to the radiation transitions; Ad and Aex: QED PT fourth (atomic PT second)-order polarization diagrams.

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:

$$\begin{split} \text{Im}\Delta E\_{ninv}(\alpha - s | A\_d) &= -\mathbb{C} \frac{e^2}{4\pi} \int \int \int dr\_1 dr\_2 dr\_3 dr\_4 \sum (\frac{1}{\omega\_{mn} + \omega\_{a\_s}} + 1) \\ \frac{1}{\omega\_{mn} - \omega\_{a\_s}} | \Psi^+\_n(r\_1) \Psi^+\_m(r\_2) \Psi^+\_s(r\_3) \Psi^+\_n(r\_4) (1 - \alpha\_1 a\_2) / r\_{12} . \end{split} \tag{19}$$
 
$$\begin{split} \frac{1}{\omega\_{mn}} [(\alpha\_3 \alpha\_4 - (\alpha\_3 n\_{34})(\alpha\_4 n\_{34})) / r\_{34} \cdot \sin\left[\omega\_{a\_n} (r\_{12} + r\_{34}) + \omega\_{a\_s} \cdot \right. \\ \cos\left[\omega\_{a\_s} (r\_{12} + r\_{34})\right] [(1 + (\alpha\_3 n\_{34})(\alpha\_4 n\_{34}))] | \Psi\_m(r\_3) \Psi\_n(r\_4) \Psi\_n(r\_2) \Psi\_s(r\_1) \end{split} \tag{10}$$

where C is the gauge constant, and f is the boundary of the closed shells.

Realizing a principle of minimization of the functional ImΔEninv, 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:

$$\text{Im}\Delta E = -\frac{e^2}{4\pi} \sum\_{\substack{\alpha > n > f \\ [\alpha < n \le f]}} V\_{\alpha n \alpha n'}^{|\omega\_{\text{un}}|} \tag{20}$$

where X α>n>f is for electron and X α<n ≤ f for vacancy, and Vjωαn<sup>j</sup> <sup>α</sup>nα<sup>n</sup> is determined as follows:

$$\Psi\_{\vec{\eta}\vec{k}}^{|\omega|} = \iint dr\_1 dr\_2 \Psi\_i^\*(r\_1) \Psi\_{\vec{\eta}}^\*(r\_2) \frac{\sin|\omega| r\_{12}}{r\_{12}} (1 - a\_1 a\_2) \Psi\_k^\*(r\_2) \Psi\_l^\*(r\_1) \tag{21}$$

The individual terms of the sum Eq. (21) represent the contributions of different channels and probability, for instance, of the dipole <sup>α</sup>-n transition as <sup>P</sup>α<sup>n</sup><sup>e</sup> <sup>1</sup> <sup>4</sup><sup>π</sup> <sup>V</sup>jωα<sup>n</sup> <sup>j</sup> <sup>α</sup>nα<sup>n</sup>; the probability with accounting for the core polarization correction is <sup>P</sup><sup>α</sup><sup>n</sup> <sup>e</sup> <sup>1</sup> <sup>4</sup><sup>π</sup> � {Vjωα<sup>n</sup> <sup>j</sup> <sup>α</sup>nα<sup>n</sup> þ ðV<sup>d</sup>þex pol Þαnα<sup>n</sup>}. The total probability of a λ-pole transition is usually represented as a sum of the electric PE <sup>λ</sup> and magnetic PM λ

parts. The electric (or magnetic) λ-pole transition γ ! δ connects two states with parities which by λ (or λ + 1) units. In our designations,

$$\begin{array}{ll} P^{E}\_{\lambda}(\gamma \to \delta) = \ 2(2j+1) \mathcal{Q}^{E}\_{\lambda}(\gamma \delta; \gamma \delta) & \mathcal{Q}^{E}\_{\lambda} = \mathcal{Q}^{\mathbb{C}ul}\_{\lambda} + \mathcal{Q}^{Br}\_{\lambda, \ \lambda - 1} + \mathcal{Q}^{Br}\_{\lambda, \ \lambda + 1} \\ P^{M}\_{\lambda}(\gamma \to \delta) = \ 2(2j+1) \mathcal{Q}^{M}\_{\lambda}(\gamma \delta; \gamma \delta) & \mathcal{Q}^{M}\_{\lambda} = \mathcal{Q}^{Br}\_{\lambda, \ \lambda}. \end{array} \tag{22}$$

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: j1j <sup>2</sup>½J� ! j 1j <sup>2</sup>½J�):

$$P(\lambda|j\_1j\_2[\overline{l}] | \, \overline{j}\_1j\_2[\overline{l}]) = (\overline{l}) \begin{Bmatrix} \lambda ... J ... \overline{J} \\ j\_2 ... \overline{j}\_1 ... \overline{j}\_1 \end{Bmatrix} P(\lambda | 1\overline{1}) (\overline{j}\_1) . \tag{23}$$

It is worth noting that all relativistic atomic calculations are usually carried out in the jjcoupling scheme. The transition to the intermediate-coupling scheme is realized by diagonalization of the M matrix, but usually only ReM should be diagonalized. The important simplified moment of the procedure is connected with converting the imaginary part by means of the matrix of eigenvectors {Cmk}, obtained by diagonalization of ReM:

$$\operatorname{Im} M\_{mk} = \sum\_{\vec{\imath}} \mathbb{C}\_{mi}^\* M\_{\vec{\imath}} \operatorname{C}\_{\vec{\imath}k} \tag{24}$$

where Mij are the matrix elements in the jj-coupling scheme, and Mmk in the intermediatecoupling scheme representation. The procedure is correct to terms of the order of Im M/Re 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

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

ðððð

<sup>s</sup> ðr3ÞΨ<sup>þ</sup>

cos ½ωα<sup>n</sup> ðr<sup>12</sup> þ r34Þ�ð1 þ ðα3n34Þðα4n34ÞÞ�}Ψmðr3ÞΨαðr4ÞΨnðr2ÞΨsðr1Þ

Realizing a principle of minimization of the functional ImΔEninv, 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

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

4π

<sup>j</sup> <sup>ð</sup>r2<sup>Þ</sup> sin <sup>j</sup>ωjr<sup>12</sup> r<sup>12</sup>

The individual terms of the sum Eq. (21) represent the contributions of different channels and

X α>n>f ½α<n ≤ f �

for vacancy, and Vjωαn<sup>j</sup> <sup>α</sup>nα<sup>n</sup> is determined as follows:

ð1 � α1α2ÞΨ�

<sup>4</sup><sup>π</sup> � {Vjωα<sup>n</sup> <sup>j</sup> <sup>α</sup>nα<sup>n</sup> þ ðV<sup>d</sup>þex

dr1dr2dr3dr<sup>4</sup>

<sup>n</sup> ðr4Þð1 � α1α2Þ=r12�

<sup>X</sup><sup>ð</sup> <sup>1</sup> ωmn þ ωα<sup>s</sup>

þ

<sup>V</sup>jωαn<sup>j</sup> <sup>α</sup>nα<sup>n</sup>, <sup>ð</sup>20<sup>Þ</sup>

<sup>k</sup> ðr2ÞΨ�

<sup>l</sup> ðr1Þ ð21Þ

<sup>4</sup><sup>π</sup> <sup>V</sup>jωα<sup>n</sup> <sup>j</sup> <sup>α</sup>nα<sup>n</sup>; the probability with

pol Þαnα<sup>n</sup>}. The total prob-

<sup>λ</sup> and magnetic PM

λ

ð19Þ

4π

{½ðα3α<sup>4</sup> � ðα3n34Þðα4n34ÞÞ=r<sup>34</sup> � sin ½ωα<sup>n</sup> ðr<sup>12</sup> þ r34Þ þ ωα<sup>n</sup> �

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

part of electron energy shift, which is defined in the PT lowest order as follows:

ImΔ<sup>E</sup> ¼ � <sup>e</sup><sup>2</sup>

α<n ≤ f

<sup>i</sup> ðr1ÞΨ�

ability of a λ-pole transition is usually represented as a sum of the electric PE

dr1dr2Ψ�

accounting for the core polarization correction is <sup>P</sup><sup>α</sup><sup>n</sup> <sup>e</sup> <sup>1</sup>

probability, for instance, of the dipole <sup>α</sup>-n transition as <sup>P</sup>α<sup>n</sup><sup>e</sup> <sup>1</sup>

mðr2ÞΨ<sup>þ</sup>

where C is the gauge constant, and f is the boundary of the closed shells.

has been calculated, which is as follows:

140 Recent Studies in Perturbation Theory

1 ωmn � ωα<sup>s</sup>

be found in Refs. [7, 14–20, 47–49].

Method of calculation

where X

α>n>f

is for electron and X

V<sup>j</sup>ω<sup>j</sup> ijkl ¼ ðð

ImΔEninvð<sup>α</sup> � <sup>s</sup>jAd޼�<sup>C</sup> <sup>e</sup><sup>2</sup>

<sup>α</sup> ðr1ÞΨ<sup>þ</sup>

ÞΨ<sup>þ</sup>

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<sup>+</sup> states can be treated one- and three-QP states of electrons (6s) and vacancy (5d<sup>1</sup> ) above the core of the closed shells 5d106s2 . 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 (5d107p(P1/2,P3/2)-5d106s(S1/2), 5d107p(P1/2,P3/2)-5d107s(S1/2)), and E2 (5d9 6s<sup>2</sup> (D5/2, D3/2)-5d106s (S1/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).

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


Theoretical data—Hartree-Fock (HF), Dirac-Fock (DF) [31]; relativistic PT with the empirical model potential approximation (RPT-MP) [18]; relativistic PT-RPT-EA (this work); experimental data—Moore (NBS, Washington) [32] (see text).

Table 1. The energies of the 5d<sup>9</sup> 6s<sup>2</sup> (D5/2,D3/2)-5d106s (S1/2), 5d107p(P1/2,P3/2)-5d106s(S1/2), 5d107p(P1/2,P3/2)-5d107s(S1/2), 5d<sup>9</sup> 6s<sup>2</sup> (D5/2,D3/2)-5d106s (S1/2) transitions in Hg+ (Ry).


HF, Hartree-Fock data; DF, Dirac-Fock data; DF (Eexp), DF data using the experimental transitions energies [31]; relativistic perturbation theory with the empirical model potential approximation RPT-MP [18]; relativistic PT-RPT-EA (this work); experimental data—Moore (NBS, Washington) [32] (see text).

Table 2. Probabilities of the transitions 5d107p(P1/2,P3/2)-5d106s(S1/2), 5d107p(P1/2,P3/2)-5d107s(S1/2) in Hg+ (in s<sup>1</sup> ). Relativistic Perturbation Theory Formalism to Computing Spectra and Radiation Characteristics: Application to… http://dx.doi.org/10.5772/intechopen.69102 143


HF, Hartree-Fock data; DF, Dirac-Fock data; DF (Eexp), DF data using the experimental transitions energies [31]; relativistic perturbation theory with the empirical model potential approximation (RPT-MP) [18]; relativistic PT-RPT-EA (this work); experimental data—Moore (NBS, Washington) [32] (see text).

Table 3. The E2 probabilities of the 5d<sup>9</sup> 6s<sup>2</sup> (D5/2,D3/2)-5d106s (S1/2) transition in Hg+ (in s<sup>1</sup> ).

vacancy (5d<sup>1</sup>

142 Recent Studies in Perturbation Theory

Table 1. The energies of the 5d<sup>9</sup>

5d<sup>9</sup> 6s<sup>2</sup> ) above the core of the closed shells 5d106s2

(5d107p(P1/2,P3/2)-5d106s(S1/2), 5d107p(P1/2,P3/2)-5d107s(S1/2)), and E2 (5d9

performed using the modified atomic code "Superatom-ISAN."

(S1/2)) probabilities of the transitions in the spectrum of Hg+

optimized RPT using relativistic energy approach (RPT-EA).

6s<sup>2</sup>

work); experimental data—Moore (NBS, Washington) [32] (see text).

(D5/2,D3/2)-5d106s (S1/2) transitions in Hg+ (Ry).

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

In Tables 1–3, we present the experimental (NIST) [32] and theoretical energies, electric E1

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

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

Method E6s 7P1/2-6S1/2 7P3/2-6S1/2 7P1/2-7S1/2 7P3/2-7S1/2 D3/2-S1/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 Theoretical data—Hartree-Fock (HF), Dirac-Fock (DF) [31]; relativistic PT with the empirical model potential approximation (RPT-MP) [18]; relativistic PT-RPT-EA (this work); experimental data—Moore (NBS, Washington) [32] (see text).

Method 7P3/2-6S1/2 7P1/2-6S1/2 7P3/2-7S1/2 7P1/2-7S1/2 7P3/2-6S1/2 HF 4.75<sup>10</sup><sup>6</sup> 4.75106 3.65<sup>10</sup><sup>7</sup> 3.65<sup>10</sup><sup>7</sup> 3.65<sup>10</sup><sup>7</sup> DF 8.45<sup>10</sup><sup>7</sup> 1.67107 6.89<sup>10</sup><sup>7</sup> 6.89<sup>10</sup><sup>7</sup> 4.71<sup>10</sup><sup>7</sup> DF (Eexp) 1.17108 2.04107 1.10<sup>10</sup><sup>8</sup> 1.10<sup>10</sup><sup>8</sup> 5.52<sup>10</sup><sup>7</sup> RPT-MP 1.49108 2.31107 1.41<sup>10</sup><sup>8</sup> 1.41<sup>10</sup><sup>8</sup> 6.33<sup>10</sup><sup>7</sup> RPT-EA 1.51108 2.33107 1.43<sup>10</sup><sup>8</sup> 1.43<sup>10</sup><sup>8</sup> 6.35<sup>10</sup><sup>7</sup> Exp. 1.53<sup>10</sup><sup>8</sup> 2.35107 1.44<sup>10</sup><sup>8</sup> 1.44<sup>10</sup><sup>8</sup> 6.37<sup>10</sup><sup>7</sup> HF, Hartree-Fock data; DF, Dirac-Fock data; DF (Eexp), DF data using the experimental transitions energies [31]; relativistic perturbation theory with the empirical model potential approximation RPT-MP [18]; relativistic PT-RPT-EA (this

Table 2. Probabilities of the transitions 5d107p(P1/2,P3/2)-5d106s(S1/2), 5d107p(P1/2,P3/2)-5d107s(S1/2) in Hg+ (in s<sup>1</sup>

(D5/2,D3/2)-5d106s (S1/2), 5d107p(P1/2,P3/2)-5d106s(S1/2), 5d107p(P1/2,P3/2)-5d107s(S1/2),

. The interaction "quasiparticle core"

6s<sup>2</sup> (D5/2, D3/2)-5d106s

).

. The theoretical results are

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].

## Author details

Alexander V. Glushkov\*, Olga Yu. Khetselius, Andrey A. Svinarenko, Vasily V. Buyadzhi, Valentin B. Ternovsky, Anna A. Kuznetsova and Petr G. Bashkarev

\*Address all correspondence to: glushkovav@gmail.com

Odessa State Environmental University, Odessa, Ukraine

## References


[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

Author details

144 Recent Studies in Perturbation Theory

References

2003. p. 680

University Press; 2007. p. 780

Physics. 2007;79:291–328

Alexander V. Glushkov\*, Olga Yu. Khetselius, Andrey A. Svinarenko, Vasily V. Buyadzhi,

[1] Wilson S. Handbook on Molecular Physics and Quantum Chemistry. Chichester: Wiley;

[2] Dyall KG, Faegri R, Jr. Introduction to Relativistic Quantum Theory. Oxford: Oxford

[3] Bartlett RJ, Musiał MN. Coupled-cluster theory in quantum chemistry. Review of Modern

[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

[6] Parpia FA, Froese-Fischer C, Grant IP. Generalized relativistic atomic structure package:

[7] Glushkov AV. Relativistic Quantum Theory. Relativistic Quantum Mechanics of Atomic

[8] Kohn W, Sham S. Quantum density oscillations in an inhomogeneous electron gas.

[9] Johnson WR, Sapistein J, Blundell SA. Atomic structure calculations associated with PNC

[10] Desclaux JP. A multiconfiguration relativistic DIRAC-FOCK program. Computer Physics

[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

Strong External Fields. Vol. 1. New York, NY: Kluwer; 1997. pp. 101–108

GRASP. Computer Physics Communications. 1996;94(2):249–270

experiments in atomic Cs. Physica Scripta T. 1993;46:184–192

quantum electrodynamics. Physica Scripta. 1989;39:682–688

Systems. Odessa: Astroprint; 2008. p. 700

Physical Review A. 1965;137:1697–1710

Communications. 1975;9(1):31–45

Valentin B. Ternovsky, Anna A. Kuznetsova and Petr G. Bashkarev

\*Address all correspondence to: glushkovav@gmail.com

Odessa State Environmental University, Odessa, Ukraine


[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

[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

Berlin: Springer; 2015. pp. 54–76

146 Recent Studies in Perturbation Theory

Review A. 2009;79:022512

ions. Physical Review A. 2007;75:062501

Optics and Spectroscopy. 1989;67:16–22

2005;71:052506

NIST; 2004

[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.

[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 Theo-

[24] Safronova UI, Johnson W, Derevianko A. Relativistic many-body calculations of energy levels, hyperfine constants, electric-dipole matrix elements, static polarizabilities for

[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

[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.

[28] Yerokhin V, Artemyev A, Shabaev VM. QED treatment of electron correlation in Li-like

[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-

[31] Ostrovsky VN, Sheynerman SA. Radiation transitions in external shells of the ion Hg+.

[32] Martin W. NIST Spectra Database, Version 2.0 (http://physics.nist.gov.asd). Washington:

[33] Ivanov LN, Ivanova EP. Extrapolation of atomic ion energies by model potential method:

nonconservation amplitude in cesium. Physical Reviews A. 2001;63:044103

Na-like spectra. Atomic Data and Nuclear Data Tables. 1979;24(2):95–121

retical Chemistry and Physics. Vol. 29. Berlin: Springer; 2015. pp. 197–218

alkali-metal atoms. Physical Review A. 1999;60:4476–4486


[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

[49] Glushkov AV. Negative ions of inert gases. JETP Letters. 1992;55(2):97–100

theory. Journal of Structural Chemistry. 1990;31(4):529–533

for molecules. Soviet Physics Journal. 1991;34(10):871–876

USSR Academy of Sciences (ISAN); 1992. AS-N2:1–10

Journal of Quantum Chemistry. 2004;99(5):936–939

tum optics. European Physics Journal. 2008;T160:195–204

copy of the USSR Academy of Sciences (ISAN); 1992. AS-N3:1–12

of Sciences (ISAN); 1992. AS-N1:1–16

Journal. 1990;33(1):1–8

148 Recent Studies in Perturbation Theory

AS-N4:1–16

L379-L396

[50] Glushkov AV. Relativistic polarization potential of a many-electron atom. Soviet Physics

[51] Glushkov AV. Correction for exchange and correlation effects in multielectron system

[52] Glushkov AV. Relativistic multiconfiguration time-dependent self-consistent-field theory

[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.

[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

[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

[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 Spectros-

[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):

[58] Glushkov AV, Ambrosov SV, Ignatenko AV, Korchevsky DA. DC strong field Stark effect for non-hydrogenic atoms: New consistent quantum mechanical approach. International

[59] Glushkov AV, Malinovskaya SV, Chernyakova YuG, Svinarenko AA. Cooperative laserelectron-nuclear processes: QED calculation of electron satellites spectra for multi-charged

[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 quan-

[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

ion in laser field. International Journal of Quantum Chemistry. 2004;99(5):889–893

