## Exchange contribution to the variance of the excitation energy of the electron shell of the daughter atom in doubleβ decay

∗ †

K. S. Tyrin and M. I. Krivoruchenko

National Research Centre 'Kurchatov Institute', 123182 Moscow, Russia

The excitation of electron shell of a daughter atom in a neutrinoless double β -decay causes change in shape of the total energy peak of β electrons at the end of the energy spectrum. The main parameters of the modified distribution are the average energy and variance of the excitation energy of the electron shell. We derive an expression for the variance taking into account exchange effects and make numerical estimates of the average excitation energy and variance based on the non-relativistic Roothaan-Hartree-Fock method and the relativistic Dirac-Hartree-Fock method implemented in the General Relativistic Atomic Structure Package (G rasp 2018). Estimates are made for eleven isotopes, two-neutrino doubleβ decay of which is observed experimentally. The results are determined by the first two negative moments of the electron radii in the parent atom. The values obtained for change in the peak shape can be used to parameterize the energy distribution of β electrons, taking into account the excitation of the electron shell of the daughter atom.

∗ E-mail address: Tyrin KS@nrcki.ru

† E-mail address: mikhail.krivoruchenko@itep.ru

## I. INTRODUCTION

The study of β processes sensitive to mixing, mass, and nature of neutrinos (Dirac/Majorana) is a promising way in the search for extensions of the Standard Model (SM). The neutrinoless doubleβ (0 ν 2 β ) decay is important for testing the conservation of the total lepton number and the nature of neutrinos (for a review see [1]). Observing 0 ν 2 β decay indicates the presence of a Majorana neutrino mass, according to the 'black box theorem' [2, 3]. In the quark sector of the SM, proton decay and neutron-antineutron oscillations play a similar role in testing the total baryon number conservation. The limitation of the period of neutron-antineutron oscillations restricts the neutron Majorana mass to ∆ m &lt; . n 0 8 × 10 -33 m n , where m n = 939 57 MeV/c . 2 is the neutron Dirac mass [4]. The Majorana neutrino mass is limited by a fraction of electronvolts [5].

In 0 ν 2 β decay, the decay energy Q is distributed between two β electrons. Thus, their total energy equals the value of Q . The occurrence of a peak at energy Q against the background of a continuous distribution of electron energy in a two-neutrino doubleβ decay (2 ν 2 β ) is the signature of 0 ν 2 β -decay. The collaborations CUORE [6, 7], EXO [8], KamLAND-Zen [9] and LEGEND [10] are searching for 0 ν 2 β -decay to clarify the Dirac/Majorana nature of neutrinos and estimate their Majorana mass. Multielectron modes can be seen while examining the β electron spectrum of the SuperNEMO collaboration [11, 12].

The limited energy resolution of detector, together with the nuclear recoil and excitation of the parent atom's electron shell, modify the energy spectrum of β electrons. In β processes, a change in the charge of the nucleus acts on the electron shell as a sudden perturbation, forcing the electrons of a daughter ion to move to excited levels (shake-up) or ionization states of a continuous spectrum (shake-off) with a certain probability. The energy probability distribution splits into two terms: the first in which all electrons preserve quantum numbers and the reaction energy is not modified, and the second in which excitation and ionization occur. In the second situation, the effective decay energy decreases, which affects the spectrum of β electrons around the boundary value, sensitive to the Majorana neutrino mass.

Modeling the probability distribution of β electrons in 0 ν 2 β -decay requires estimates of the average excitation energy and variance of the daughter ion's electron shell. The excitation effects are discussed in [13, 14]. The average excitation energy and variance of excitation energy of the electron shell are estimated, taking exchange effects into account. In Section II.B, we provide derivation of the formulas used in [14] to calculate exchange effects. The variance is estimated with the use of the matrix elements of r -1 and r -2 , which account for exchange effects. Section III presents the matrix elements of r -1 and r -2 for eleven isotopes, whose doubleβ decays are observed experimentally. The non-relativistic Roothaan-Hartree-Fock (RHF) method [15] and the software package G rasp 2018 [16-18], based on the relativistic Dirac-Hartree-Fock (DHF) method, are employed for calculations.

## II. EXCITATION ENERGY OF ELECTRON SHELL AND VARIANCE

The paper discusses 0 ν 2 β decay with decay energy Q /greaterorsimilar mc 2 , where m is the electron mass and c is the speed of light. This criterion is fulfilled for all eleven isotopes, whose 2 ν 2 β -decays were observed (see Table I). The condition Q /greaterorsimilar mc 2 indicates that the velocity of β electrons approaches the speed of light and exceeds the typical velocity ∼ αZ 1 / 3 c of electrons in multielectron atoms.

In β -decay problems, an approximation of a sudden perturbation gives adequate precision. For an accurate approximation, the flight duration τ f of β electron through the atomic shell should be smaller than the period τ o of electron's revolution around the nucleus (see, e.g., [30]). The energy of β electrons follows a continuous distribution. When one of the β electrons carries energy alongside the phase space volume boundary, the velocity of the second electron is minimal. In such a situation, the approach of a sudden perturbation of the atomic shell is inequitable. The usefulness of approximation is confined to the middle region of the energy spectrum, where β electron energies are comparable and neutrinos have low energy. In 0 ν 2 β decay, only electrons carry away the decay energy, and the concept of a sudden perturbation has a broader range of validity.

In what follows, the atomic system of units /planckover2pi1 = m = e 2 = 1, c = 137 is employed. The β electron flight duration through the shell, τ f , can be estimated as the ratio of the K-orbit radius a K ∼ 1 /Z to the velocity of the β electron pc /E 2 ∼ c for Q /greaterorsimilar c 2 . The time the electron takes to rotate around the nucleus is τ o ∼ 2 π/Z 2 . The perturbation is sudden if the electron on the shell does not have enough time to perform one circle around the nucleus during the flight of the β electron. The related instantaneous parameter is

$$\xi _ { \text{inst} } = \frac { \tau _ { f } } { \tau _ { o } } \sim \frac { Z } { 2 \pi c }.$$

For a germanium atom of Z = 32, ξ inst ∼ 1 / 25. For outer orbits where the shielded charge of the nucleus is of the

TABLE I. Variance of the excitation energy of the electron shells of the daughter atom for eleven isotopes, the 0 ν 2 β -decays of which are observed experimentally. The second column contains the values of the mass difference, Q , of the neutral atoms involved in the decay. The fourth column presents the average excitation energies calculated by the DHF method using the G rasp 2018 package. The fifth and sixth columns show the exchange contributions to the variance and the total variance obtained in the framework of the non-relativistic RHF method. The last two columns report the exchange contributions and the total variance obtained using the relativistic DHF method implemented into the G rasp 2018.

| Process   | Process   | Process   | Q [keV]      | Ref.   | C DHF [eV]   | ∆ D 1 / 2 RHF [keV]   | D 1 / 2 RHF [keV]   | ∆ D 1 / 2 DHF [keV]   | D 1 / 2 DHF [keV]   |
|-----------|-----------|-----------|--------------|--------|--------------|-----------------------|---------------------|-----------------------|---------------------|
| 48 20 Ca  | →         | 48 22 Ti  | 4267.98(32)  | [19]   | 283          | -0.05                 | 1.61                | -0.05                 | 1.64                |
| 76 32 Ge  | →         | 76 34 Se  | 2039.006(50) | [21]   | 365          | -0.10                 | 2.62                | -0.11                 | 2.77                |
|           |           |           | 2039.061(7)  | [20]   |              |                       |                     |                       |                     |
| 82 34 Se  | →         | 82 36 Kr  | 2997.9(3)    | [23]   | 375          | -0.11                 | 2.79                | -0.12                 | 2.97                |
| 96 40 Zr  | →         | 96 42 Mo  | 3356.097(86) | [24]   | 403          | -0.15                 | 3.29                | -0.15                 | 3.60                |
| 100 42 Mo | →         | 100 44 Ru | 3034.40(17)  | [22]   | 416          | -0.16                 | 3.46                | -0.17                 | 3.83                |
| 116 48 Cd | →         | 116 50 Sn | 2813.50(13)  | [25]   | 448          | -0.20                 | 3.97                | -0.21                 | 4.54                |
| 128 52 Te | →         | 128 54 Xe | 865.87(131)  | [26]   | 468          | -0.21                 | 4.32                | -0.23                 | 5.05                |
| 130 52 Te | →         | 130 54 Xe | 2526.97(23)  | [25]   | 468          | -0.21                 | 4.32                | -0.23                 | 5.05                |
| 136 54 Xe | →         | 136 56 Ba | 2457.83(37)  | [27]   | 475          | -0.22                 | 4.49                | -0.25                 | 5.32                |
| 150 60 Nd | →         | 150 62 Sm | 3371.38(20)  | [28]   | 515          |                       |                     | -0.29                 | 6.20                |
| 238 92 U  | →         | 238 94 Pu | 1437.3       | [29]   | 817          |                       |                     | -0.63                 | 13.95               |

order of one, a substitution of Z → Z eff ∼ 1 should be made, resulting in an additional decline in ξ inst . Thus, for Q /greaterorsimilar c 2 , the concept of suddenness of the perturbation is quite realistic.

## A. Multielectron atoms

In relativistic theory, the Hamiltonian of N electrons in an atom with a nuclear charge Z is written as follows [18]

$$\hat { H } _ { Z, N } & = \sum _ { i = 1 } ^ { N } \left ( \alpha _ { i } \mathbf p _ { i } c + \beta _ { i } c ^ { 2 } - \frac { Z } { | \mathbf x _ { i } | } \right ) + \sum _ { i < j } ^ { N } V _ { \text{CG} / \text{CB} } ( \mathbf x _ { i } - \mathbf x _ { j } ), \\.. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &.. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &.. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &. &.$$

where atomic units are used, α and β are the Dirac gamma matrices, V CG CB / ( x i -x j ) is the gauge-dependent lowest-order relativistic interaction potential of two electrons. Lorenz and Coulomb gauges are usually used to model electromagnetic interaction of electrons in multielectron atoms.

The Coulomb-Gaunt potential refers to the interaction potential of electrons in the Lorenz gauge [31]:

$$V _ { C G } ( x _ { i } - x _ { j } ) = \frac { 1 } { | x _ { i } - x _ { j } | } \left ( 1 - \frac { \alpha _ { i } \cdot \alpha _ { j } } { c ^ { 2 } } \right ). & & ( \text{II.3} ) \\ \text{international as the alaotron volository anorator differential by the second of light } \text{$\text{$the$national$}$}$$

The matrices α i can be interpreted as the electron velocity operators divided by the speed of light c . The potential V CG ( x i -x j ) describes the lowest-order photon-exchange static interaction of two electrons. The second term is responsible for the energy of magnetostatic interaction; it is of the same order ∼ v /c 2 2 as the retardation correction to the one-photon exchange. For the innermost shells v ∼ 1 /r ∼ Z , where r = | x | .

The Coulomb-Breit static potential can be obtained using the Coulomb gauge [32]:

$$V _ { \text{CB} } ( x _ { i } - x _ { j } ) & = \frac { 1 } { | x _ { i } - x _ { j } | } \left ( 1 - \frac { \alpha _ { i } \cdot \alpha _ { j } + \alpha _ { i } \cdot \mathbf n \alpha _ { j } \cdot \mathbf n } { 2 c ^ { 2 } } \right ), \\ \dots \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdot \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdot \quad \cdot \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdclots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \ddots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \dots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \vdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \ \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \ \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots\quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots\ { \mathfrak { m } }$$

where n = ( x i -x j ) / | x i -x j | . V CB ( x i -x j ) is precise to order 1 /c 2 , with the retardation modifications appearing at order 1 /c 4 . The equivalence of potentials (II.3) and (II.4) in the order 1 /c 2 is demonstrated in Ref. [4].

These potentials are used to generate initial approximations for a variety of atomic physics problems. The noncovariant Coulomb gauge is appropriate for directly solving the constraint field equations, allowing a transition from classical electrodynamics to quantum electrodynamics (QED) via canonical quantization. The equivalence of the covariant Lorenz gauge and the non-covariant Coulomb gauge implies relativistic covariance of QED. Quantum electrodynamics of multielectron atoms is a gauge invariant theory in all orders of perturbation theory in relation to atomic energy eigenvalues and electromagnetic form factors [33].

The ground state of N electrons in an atom with a nuclear charge Z is represented by | Z, N 〉 . By the definition

$$\hat { H } _ { Z, N } | Z, N \rangle = E _ { Z, N } | Z, N \rangle, \quad \ \ \ \ \ \ \ ( \text{II.5} ) \\ \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \,.$$

where E Z,N is the total energy of electron shell. The energy spectrum has discrete and continuous parts. The discrete energy eigenstates are normalized by one:

$$\langle Z, N | Z, N \rangle & = 1, & & ( \text{II.6} ) \\. &. &. &.$$

while states with different energy eigenvalues are orthogonal.

## B. Average excitation energy

Doubleβ decay increases the nuclear charge by two units. The Coulomb addition to the Hamiltonian causes a sudden perturbation, transforming the ground state | Z, Z 〉 of the initial Hamiltonian ˆ H Z,Z into a superposition of states with definite energy of the final Hamiltonian:

$$\hat { H } _ { Z + 2, Z } & = \hat { H } _ { Z, Z } - 2 \sum _ { i = 1 } ^ { Z } \frac { 1 } { r _ { i } }, \\ \text{are of $Z$ alone of the parent atom in the Column field of the author name}$$

where r i = | x i | . The average energy of Z electrons of the parent atom in the Coulomb field of the daughter nucleus with a charge of Z +2 is equal to

$$\langle Z, Z | \hat { H } _ { Z + 2, Z } | Z, Z \rangle = E _ { Z, Z } + 2 E _ { Z, Z } ^ { \mathbb { C } } / Z,$$

where E C Z,Z is the average energy of the Coulomb interaction of electrons of the parent atom with the nucleus:

$$E _ { Z, Z } ^ { \text{C} } = - \langle Z, Z | \sum _ { i = 1 } ^ { Z } \frac { Z } { r _ { i } } | Z, Z \rangle. \intertext { t i o n a r n o r n o f t h a l t o f t h a d u b t o r i o n i c a t o r m i n o l h e 1 2 3 1 4 }$$

The average value of the excitation energy of the electron shell of the daughter ion is determined by [13, 14]

$$\mathcal { C } = E _ { Z, Z } + 2 E _ { Z, Z } ^ { \text{C} } / Z - E _ { Z + 2, Z }.$$

Binding energies of neutral atoms are tabulated in Refs. [15, 34-36] and can be calculated with the use of programs designed to model structure of electron shells (see, e.g., [16, 17]). The Coulomb interaction energy E C Z,Z is tabulated in Ref. [36].

## C. Variance of excitation energy taking into account exchange effects

In simple scenarios, the probability distribution functions can be characterized by two parameters, for which it is enough to know only the first two moments. The average of the excitation energy is used as the first, and the variance of excitation energy,

$$\mathcal { D } = \langle Z, Z | \hat { H } _ { Z + 2 } ^ { 2 } | Z, Z \rangle - \langle Z, Z | \hat { H } _ { Z + 2 } | Z, Z \rangle ^ { 2 }, \\. \quad \ \ - \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
.$$

is used as the second. Taking into account Eqs. (II.5) and (II.7) we have

$$\frac { 1 } { 4 } \mathcal { D } & = \langle Z, Z | \sum _ { i = 1 } ^ { z } \sum _ { j = 1 } ^ { z } \frac { 1 } { r _ { i } } \frac { 1 } { r _ { j } } | Z, Z \rangle - \langle Z, Z | \sum _ { i = 1 } ^ { z } \frac { 1 } { r _ { i } } | Z, Z \rangle ^ { 2 }. \\ \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots.$$

Finding the average values of 1 /r i and 1 /r i 1 /r j with the help of the wave function of the parent atom allows to calculate the variance.

It should be noted that Eqs. (II.10) and (II.11) are valid for any form of the two-electron interaction potential V CG CB / ( x i -x j ). When deriving these equations, only the property

$$\frac { \partial } { \partial Z } \hat { H } _ { Z, N } = - \sum _ { i = 1 } ^ { N } \frac { 1 } { | \mathbf x _ { i } | }$$

is used.

Ignoring the exchange contribution, factorization holds for i = j

/negationslash

$$\langle Z, Z | \frac { 1 } { r _ { i } r _ { j } } | Z, Z \rangle \approx \langle Z, Z | \frac { 1 } { r _ { i } } | Z, Z \rangle \langle Z, Z | \frac { 1 } { r _ { j } } | Z, Z \rangle.$$

In this situation, the variance takes the form

$$\frac { 1 } { 4 } \mathcal { D } \approx & \sum _ { i = 1 } ^ { z } \langle Z, Z | \frac { 1 } { r _ { i } ^ { 2 } } | Z, Z \rangle - \sum _ { i = 1 } ^ { z } \langle Z, Z | \frac { 1 } { r _ { i } } | Z, Z \rangle ^ { 2 }. \\ \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \end{pmatrix}$$

To calculate the exchange contribution, it is necessary to know the two-particle density of the electron distribution in the atom. The exchange effects are taken into account in the Hartree-Fock method and its relativistic versions by averaging the operator over the antisymmetrized electron wave function. The Slater determinant,

$$\Psi _ { \alpha _ { 1 } \alpha _ { 2 } \dots \alpha _ { N } } ^ { \sigma _ { 1 } \sigma _ { 2 } \dots \sigma _ { N } } ( x _ { 1 }, x _ { 2 }, \dots, x _ { N } ) = \frac { 1 } { \sqrt { N! } } \varrho ^ { s _ { 1 } s _ { 2 } \dots s _ { N } } \psi _ { \alpha _ { s _ { 1 } } } ^ { \sigma _ { 1 } } ( x _ { 1 } ) \psi _ { \alpha _ { s _ { 2 } } } ^ { \sigma _ { 2 } } ( x _ { 2 } ) \dots \psi _ { \alpha _ { s _ { N } } } ^ { \sigma _ { N } } ( x _ { N } ),$$

represents the shell wave function, where ψ σ α s ( x ) are the wave functions of electrons; index i = 1 , ..., N numbers the spatial coordinates and bispinor indices of the orbitals; and index α represents the quantum numbers of orbitals. In the case under study, α = ( njlm ), where n is the principal quantum number, j is the total angular momentum, l = j ± 1 / 2 ≥ 0 is the orbital angular momentum, and m = -j, . . . , j is the projection of the total angular momentum on the z axis. The Slater determinant is defined as a convolution of the Levi-Civita symbol /epsilon1 s 1 s 2 ...s N = ± 1 with the product of N electron wave functions. The Levi-Civita symbol antisymmetrizes the wave function by summing over N ! permutations. In the one-determinant approximation, a fixed set of quantum numbers ( α α 1 2 . . . α N ) completely describes the state of the shell. The orthonormality condition for the wave functions is expressed as

$$\int d \mathbf x \psi _ { \alpha } ^ { \dagger } ( \mathbf x ) \psi _ { \beta } ( \mathbf x ) & = \delta _ { \alpha \beta }. \\. \quad \dots \dots \dots \dots \dots$$

/negationslash

The calculations use a two-particle density matrix i = k , which is found by integrating the product Ψ σ σ ...σ 1 2 N α α ...α 1 2 N ( x 1 , x 2 , . . . , x N )Ψ ¯ σ σ ...σ 1 ¯ 2 ¯ N ∗ α α ...α 1 2 N (¯ x 1 , ¯ x 2 , . . . , ¯ x N ) over the coordinates for x h = ¯ x h and contracting the bispinor indices σ h and ¯ σ h of orbitals with numbers h other than i and k :

$$\rho _ { \alpha _ { 1 } \alpha _ { 2 } \dots \alpha _ { N } } ^ { \sigma _ { i } r _ { k } \delta _ { i } } ( \mathbf x _ { i }, \mathbf x _ { k } ; \bar { \mathbf x } _ { k }, \bar { \mathbf x } _ { i } ) = \frac { 1 } { N ( N - 1 ) } \left ( \delta ^ { r r ^ { \prime } } \delta ^ { s s ^ { \prime } } - \delta ^ { r s ^ { \prime } } \delta ^ { s r ^ { \prime } } \right ) \psi _ { \alpha _ { r } } ^ { \sigma _ { i } } ( \mathbf x _ { i } ) \psi _ { \alpha _ { s } } ^ { \tau _ { k } } ( \mathbf x _ { k } ) \psi _ { \alpha _ { s ^ { \prime } } } ^ { \tau _ { k } * } ( \bar { \mathbf x } _ { k } ) \psi _ { \alpha _ { s ^ { \prime } } } ^ { \delta _ { i } * } ( \bar { \mathbf x } _ { i } ). \quad \\ \mathfrak { T u } \colon \colon \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \colon = \Delta \end{cases}$$

The indices r, s, r , s ′ ′ = 1 , . . . , N are summed. Assuming y = ¯ y , integrating over the variable y and further contracting the spinor indices τ and ¯, we find a single-particle density matrix: τ

$$\rho _ { \alpha _ { 1 } \alpha _ { 2 } \dots \alpha _ { N } } ^ { \sigma \bar { \sigma } } ( \mathbf x ; \bar { \mathbf x } ) = \frac { 1 } { N } \delta ^ { r r ^ { \prime } } \psi _ { \alpha _ { r } } ^ { \sigma } ( \mathbf x ) \psi _ { \alpha _ { r ^ { \prime } } } ^ { \bar { \sigma } * } ( \bar { \mathbf x } ),$$

which is used to calculate the average values of single-particle operators.

Let's represent the upper and lower components of the electron wave functions as products of the radial and angular parts:

$$\psi _ { \alpha } ( { \mathbf x } ) = \left ( \begin{array} { c } f _ { \varkappa } ( r ) \Omega ^ { l _ { \alpha } } _ { j _ { \alpha } } ( { \mathbf n } ) \\ i g _ { \varkappa } ( r ) \Omega ^ { l _ { \alpha } } _ { j _ { \alpha } } ( { \mathbf n } ) \end{array} \right ).$$

Here f κ ( r ) and g κ ( r ) are real functions, α = ( njlm , ) κ = ( njl ), l ′ = 2 j -l , n = x / | x | , and Ω l jm ( n ) is the spherical spinor, describing a spin-1 / 2 particle. Spherical spinors are normalized by the condition [37]

$$\int d o _ { \mathbf n } \Omega _ { j m } ^ { l } ( \mathbf n ) ^ { \dagger } \Omega _ { j ^ { \prime } m ^ { \prime } } ^ { l ^ { \prime } } ( \mathbf n ) = \delta ^ { l l ^ { \prime } } \delta _ { j j ^ { \prime } } \delta _ { m m ^ { \prime } }. \\ \dots \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

The average value of the sum of the products of single-particle operators O ¯ σσ ( r i ) O ¯ ττ ( r k ) is given by the expression

/negationslash

/negationslash

$$\sum _ { i \neq k } \langle Z, N | \hat { \mathcal { O } } ( r _ { i } ) \hat { \mathcal { O } } ( r _ { k } ) | Z, N \rangle = \sum _ { i \neq k } \int d x _ { i } d x _ { k } \mathcal { O } ^ { \hat { \theta } \sigma } ( r _ { i } ) \mathcal { O } ^ { \hat { \tau } } ( r _ { k } ) \rho ^ { \sigma \tau \hat { \tau } \hat { \sigma } } ( x _ { i }, x _ { k } ; x _ { k } ; x _ { i } ).$$

TABLE II. Matrix elements 〈 κ | r -1 | ¯ κ 〉 and 〈 κ | r -2 | ¯ κ 〉 with κ = ( njl ) and ¯ = (¯ κ njl ) for orbital states of electrons in calcium and germanium atoms. Numerical values are given in atomic units. The computations use non-relativistic Roothaan-Hartree-Fock wave functions [15]. The final two lines depict the diagonal matrix elements of κ = ¯ calculated by the relativistic Dirac-Fock κ method [35]. The higher and lower lines represent the total angular momentum of j = l ± 1 2, respectively. /

| 20 Ca                 | 20 Ca    | 20 Ca       | 20 Ca         | 20 Ca        | 20 Ca   | 20 Ca       | 20 Ca       | 20 Ca       | 20 Ca   | 20 Ca       |
|-----------------------|----------|-------------|---------------|--------------|---------|-------------|-------------|-------------|---------|-------------|
| 〈 κ | r - 1 | ¯ κ 〉 | 1S       | 2S          | 3S            | 4S           |         | 2P          | 3P          | 4P          |         | 3D          |
| 1S                    | 19.54    | 3.458       | 1.158         | 0.275        |         | 3.944       |             |             |         |             |
| 2S                    |          | 4.041       | 0.776         | 0.179        | 2P      |             | -0.634      |             |         |             |
| 3S                    |          |             | 1.189         | 0.150        | 3P      |             | 1.059       |             |         |             |
| 4S                    |          |             |               | 0.300        |         |             |             |             |         |             |
| [35]                  | 19.75    | 4.089       | 1.198         | 0.301        |         | 3.947 3.987 | 1.059 1.066 |             |         |             |
| 〈 κ | r - 2 | ¯ κ 〉 | 1S       | 2S          | 3S            | 4S           |         | 2P          | 3P          | 4P          |         | 3D          |
| 1S                    | 769.4    | 200.4       | 68.52         | 16.35        |         |             |             |             |         |             |
| 2S                    |          | 67.04       | 21.80         | 5.185        | 2P      | 21.49       | -5.373      |             |         |             |
| 3S                    |          |             | 8.395         | 1.921        | 3P      |             | 2.462       |             |         |             |
| 4S                    |          |             |               | 0.535        |         |             |             |             |         |             |
| [35]                  | 794.3    | 70.01       | 8.752         | 0.557        |         | 21.54       | 2.447       |             |         |             |
| 32 Ge                 | 32 Ge    | 32 Ge       | 32 Ge         | 32 Ge        | 32 Ge   | 32 Ge       | 32 Ge       | 32 Ge       | 32 Ge   | 32 Ge       |
| 〈 κ | r - 1 | ¯ κ 〉 | 1S       | 2S          | 3S            | 4S           |         | 2P          | 3P          | 4P          |         | 3D          |
| 1S                    | 31.51    | 5.917 6.950 | -2.203 -1.477 | -0.564       | 2P      | 6.890       | 1.248       | 0.231       |         |             |
| 2S 3S                 |          |             | 2.233         | -0.366 0.310 | 3P      |             | 2.104       | 0.223       | 3D      | 1.793       |
| 4S                    |          |             |               | 0.580        | 4P      |             |             | 0.437       |         |             |
| [35]                  | 32.28    | 7.179       | 2.289         | 0.591        |         | 6.917 7.107 | 2.112 2.156 | 0.441 0.448 |         | 1.783 1.796 |
| 〈 κ | r - 2 | ¯ κ 〉 | 1S 1995. | 2S          | 3S            | 4S           |         | 2P          | 3P          | 4P          |         | 3D          |
| 1S                    |          | 555.1       | -211.8        | -54.40       |         |             |             |             |         |             |
| 2S 3S                 |          | 197.2       | -71.55 30.39  | -18.31 7.508 | 2P 3P   | 64.75       | 18.57 9.613 | 3.527 1.658 | 3D      | 4.278       |
| 4S                    |          |             |               | 2.206        | 4P      |             |             | 0.497       |         |             |
| [35]                  | 2168.    | 221.3       | 34.18         | 2.482        |         | 65.39       | 9.749       | 0.507       |         | 4.233       |

In the case we are interested in, O ¯ σσ ( r ) = δ ¯ σσ O ( r ). Substituting the expression (II.17) into the right side of the equation, followed by integration over the angular variables and taking into account the orthonormality of spherical spinors gives

/negationslash

/negationslash where

$$\sum _ { i \neq k } ( Z, N | \hat { \mathcal { O } } ( r _ { i } ) \hat { \mathcal { O } } ( r _ { k } ) | Z, N ) \ & = \ \sum _ { i \neq k } \int r _ { i } ^ { 2 } d r _ { i } r _ { k } ^ { 2 } d r _ { k } \mathcal { O } ( r _ { i } ) \mathcal { O } ( r _ { k } ) \sum _ { \alpha _ { i } \alpha _ { k } } \frac { 1 } { N ( N - 1 ) } \\ & \times \ \left ( \mathcal { P } _ { \varkappa _ { i } \varkappa _ { i } } ( r _ { i } ) \mathcal { P } _ { \varkappa _ { k } \varkappa _ { k } } ( r _ { k } ) - \delta ^ { l, l _ { k } } \delta _ { j, i j _ { k } } \delta _ { m _ { i } m _ { k } } \mathcal { P } _ { \varkappa _ { i } \varkappa _ { k } } ( r _ { i } ) \mathcal { P } _ { \varkappa _ { k } \varkappa _ { k } } ( r _ { k } ) \right ), \\ \text{where}$$

$$\mathcal { P } _ { \varkappa _ { i } \varkappa _ { k } } ( r ) = f _ { \varkappa _ { i } } ( r ) f _ { \varkappa _ { k } } ( r ) + g _ { \varkappa _ { k } } ( r ) g _ { \varkappa _ { i } } ( r ).$$

The summation by electron states includes the terms κ i = κ k , which correspond to the common shell of electrons i and k . As a consequence of the antisymmetry of the total wave function, the direct and exchange contributions at m i = m k cancel each other. The diagonal terms α i = α k could be excluded from the sum, however, calculations are simplified when the summation over the quantum numbers of two electrons is uncorrelated, i.e., it includes the diagonal.

Let's assume each pair ( jl ) has a maximum of one open shell for some n . When summing over projections m i and m k , one should keep in mind that if m i takes the values from the interval -j, . . . , j , m k can only take some of them. The second line of Eq. (II.19) yields a multiplier of (2 j i +1) N κ k , where N κ ≤ 2 j +1 represents the number of electrons on a completely or partially filled shell κ . For completely filled shells, N κ = 2 j + 1. In the third line, summation by m i removes the Kronecker symbol δ m m i k . Summation by m k yields the multiplier N κ k . The expression below is

TABLE III. Matrix elements 〈 κ | r -1 | ¯ κ 〉 and 〈 κ | r -2 | ¯ κ 〉 for the orbital states of electrons in selenium and zirconium atoms. Notations are identical to those in Table II.

| 34 Se                 | 34 Se   | 34 Se   | 34 Se   | 34 Se   | 34 Se   | 34 Se   | 34 Se   | 34 Se   | 34 Se   | 34 Se   | 34 Se   | 34 Se   |
|-----------------------|---------|---------|---------|---------|---------|---------|---------|---------|---------|---------|---------|---------|
| 〈 κ | r - 1 | ¯ κ 〉 | 1S      | 2S      | 3S      | 4S      | 5S      |         | 2P      | 3P      | 4P      |         | 3D      | 4D      |
| 1S                    | 33.50   | 6.325   | -2.403  | -0.699  |         |         |         |         |         |         |         |         |
| 2S                    |         | 7.434   | -1.610  | -0.452  |         | 2P      | 7.379   | 1.368   | 0.313   |         |         |         |
| 3S                    |         |         | 2.433   | 0.379   |         | 3P      |         | 2.311   | 0.298   | 3D      | 2.040   |         |
| 4S                    |         |         |         | 0.695   |         | 4P      |         |         | 0.563   |         |         |         |
| [35]                  | 34.57   | 7.714   | 2.502   | 0.711   |         |         | 7.412   | 2.322   | 0.561   |         | 2.029   |         |
|                       |         |         |         |         |         |         | 7.645   | 2.375   | 0.573   |         | 2.044   |         |
| 〈 κ | r - 2 | ¯ κ 〉 | 1S      | 2S      | 3S      | 4S      | 5S      |         | 2P      | 3P      | 4P      |         | 3D      | 4D      |
| 1S                    | 2255.   | 631.2   | -245.8  | -71.72  |         |         |         |         |         |         |         |         |
| 2S                    |         | 225.4   | -83.48  | -24.26  |         | 2P      | 74.18   | 21.83   | 5.137   |         |         |         |
| 3S                    |         |         | 36.13   | 10.13   |         | 3P      |         | 11.55   | 2.459   | 3D      | 5.443   |         |
| 4S                    |         |         |         | 3.329   |         | 4P      |         |         | 0.857   |         |         |         |
| [35]                  | 2479.   | 257.1   | 41.26   | 3.809   |         |         | 75.05   | 11.73   | 0.861   |         | 5.395   |         |
|                       |         |         |         |         |         |         | 82.97   | 12.92   | 0.944   |         | 5.483   |         |
| 40 Zr                 | 40 Zr   | 40 Zr   | 40 Zr   | 40 Zr   | 40 Zr   | 40 Zr   | 40 Zr   | 40 Zr   | 40 Zr   | 40 Zr   | 40 Zr   | 40 Zr   |
| 〈 κ | r - 1 | ¯ κ 〉 | 1S      | 2S      | 3S      | 4S      | 5S      |         | 2P      | 3P      | 4P      |         | 3D      | 4D      |
| 1S                    | 39.49   | -7.553  | -3.023  | 1.134   | 0.301   |         |         |         |         |         |         |         |
| 2S                    |         | 8.891   | 2.022   | -0.728  | -0.192  | 2P      | 8.849   | -1.735  | -0.563  |         |         |         |
| 3S                    |         |         | 3.055   | -0.603  | -0.153  | 3P      |         | 2.950   | 0.524   | 3D      | 2.742   | -0.319  |
| 4S                    |         |         |         | 1.067   | 0.140   | 4P      |         |         | 0.951   | 4D      |         | 0.628   |
| 5S                    |         |         |         |         | 0.306   |         |         |         |         |         |         |         |
| [35]                  | 41.26   | 9.368   | 3.177   | 1.098   | 0.317   |         | 8.909   | 2.972   | 0.954   |         | 2.730   | 0.606   |
|                       |         |         |         |         |         |         | 9.307   | 3.064   | 0.978   |         | 2.755   | 0.612   |
| 〈 κ | r - 2 | ¯ κ 〉 | 1S      | 2S      | 3S      | 4S      | 5S      |         | 2P      | 3P      | 4P      |         | 3D      | 4D      |
| 1S                    | 3131.   | -889.0  | -365.2  | 137.7   | 36.48   |         |         |         |         |         |         |         |
| 2S                    |         | 321.8   | 125.7   | -47.15  | -12.49  | 2P      | 106.4   | -33.32  | -11.01  |         |         |         |
| 3S                    |         |         | 57.09   | -20.66  | -5.458  | 3P      |         | 18.61   | 5.617   | 3D      | 9.569   | -1.720  |
| 4S                    |         |         |         | 8.537   | 2.178   | 4P      |         |         | 2.562   | 4D      |         | 0.745   |
| 5S                    |         |         |         |         | 0.657   |         |         |         |         |         |         |         |
| [35]                  | 3577.   | 387.3   | 68.702  | 10.23   | 0.810   |         | 108.2   | 19.06   | 2.616   |         | 9.513   | 0.701   |
|                       |         |         |         |         |         |         | 124.8   | 21.83   | 2.972   |         | 9.714   | 0.719   |

obtained by summing over projections m i and m k , integrating radial variables r i and r k , and summing over i = k :

/negationslash

/negationslash

$$\sum _ { i \neq k } ( Z, N | \hat { \mathcal { O } } ( r _ { i } ) \hat { \mathcal { O } } ( r _ { k } ) | Z, N \rangle = ( \sum _ { \varkappa } \mathcal { N } _ { \varkappa } \langle \varkappa | \mathcal { O } ( r ) | \varkappa \rangle ) ^ { 2 } - \sum _ { \varkappa \bar { \varkappa } } \min ( \mathcal { N } _ { \varkappa }, \mathcal { N } _ { \bar { \varkappa } } ) \langle \varkappa | \mathcal { O } ( r ) | \bar { \varkappa } \rangle ) ^ { 2 },$$

where κ = ( njl ), ¯ = (¯ κ njl ) and

$$\langle \varkappa | \mathcal { O } ( r ) | \bar { \varkappa } \rangle = \int r ^ { 2 } d r \mathcal { O } ( r ) \mathcal { P } _ { \varkappa \bar { \varkappa } } ( r ). \quad & & & & & & & & & & & & & & & & & & &$$

This equation is appropriate for calculating the average values of operators that do not mix the upper and lower components of bispinors. Summing the same terms in Eq. (II.19) cancels the multiplier 1 / N N ( ( -1)). In the second term on the right side of Eq. (II.20), the summation is carried out over the quantum numbers n, n, j ¯ , and l . The averaging is performed over the radial density of electrons in the parent atom N = Z . In non-relativistic theory where g κ ( r ) = 0, the radial part of wave functions f κ ( r ) is exclusively dependent on n and l . For completely filled shells, summation by j = l ± 1 / 2 enables the replacement of | njl 〉 → | nl 〉 and (2 j +1) → 2(2 l +1).

A single-particle density matrix is used to calculate the average value of the sum of single-particle operators ˆ ( O r i ):

$$\sum _ { i } \langle Z, N | \hat { \mathcal { O } } ( r _ { i } ) | Z, N \rangle = \sum _ { \varkappa } \mathcal { N } _ { \varkappa } \langle \varkappa | \mathcal { O } ( r ) | \varkappa \rangle. \quad \quad \quad$$

The variance consists of the sum of the diagonal and off-diagonal components 〈 Z, N |O ˆ ( r i ) ˆ ( O r k ) | Z, N 〉 minus the square of the sum of the components 〈 Z, N |O ˆ ( r i ) | Z, N 〉 . In the case under consideration, O ( r ) = r -1 . The first and

TABLE IV. Matrix elements 〈 κ | r -1 | ¯ κ 〉 and 〈 κ | r -2 | ¯ κ 〉 for the orbital states of electrons in molybdenum and cadmium atoms. Notations are identical to those in Table II.

| 42 Mo                 | 42 Mo   | 42 Mo   | 42 Mo   | 42 Mo       | 42 Mo             | 42 Mo   | 42 Mo       | 42 Mo        | 42 Mo        | 42 Mo   | 42 Mo       | 42 Mo       |
|-----------------------|---------|---------|---------|-------------|-------------------|---------|-------------|--------------|--------------|---------|-------------|-------------|
| 〈 κ | r - 1 | ¯ κ 〉 | 1S      | 2S      | 3S      | 4S          | 5S                |         | 2P          | 3P           | 4P           |         | 3D          | 4D          |
| 1S                    | 41.49   | 7.962   | 3.231   | -1.255      | 0.321 0.204 0.163 | 3P      | 9.339       |              | -0.626       |         |             |             |
| 2S                    |         | 9.378   | 2.160   | -0.803      |                   | 2P      |             | -1.858       |              |         |             |             |
| 3S                    |         |         | 3.264   | -0.665      |                   |         |             | 3.164        | 0.582        | 3D      | 2.970       | -0.361      |
| 4S                    |         |         |         | 1.171       | -0.149            | 4P      |             |              | 1.052        | 4D      |             | 0.714       |
| 5S [35]               |         |         |         |             | 0.327 0.322       |         | 9.412 9.879 | 3.190        | 1.059        |         | 2.958       | 0.695       |
| - 2                   | 43.55   | 9.939   | 3.409   | 1.209       |                   |         |             | 3.300        | 1.089        |         | 2.987       | 0.705       |
| 〈 κ | r | ¯ κ 〉     | 1S      | 2S      | 3S      | 4S          | 5S                |         | 2P          | 3P           | 4P           |         | 3D          | 4D          |
| 1S                    | 3455.   | 984.9   | 410.3   | -160.1      | 40.94             |         |             |              |              |         |             |             |
| 2S                    |         | 357.8   | 141.7   | - 55.02     | 14.06             | 2P      | 118.4       | -37.69       | -13.17       |         |             |             |
| 3S                    |         |         | 65.20   | - 24.42     | 6.223             | 3P      |             | 21.34        | 6.697        | 3D      | 11.17       | -2.120      |
| 4S                    |         |         |         | 10.41       | -2.564            | 4P      |             |              | 3.157        | 4D      |             | 0.965       |
| 5S                    |         |         |         |             | 0.748             |         |             |              |              |         |             |             |
| [35]                  | 4005.   | 439.4   | 80.03   | 12.73       | 0.830             |         | 120.7       | 21.93        | 3.243        |         | 11.11       | 0.930       |
| 48 Cd                 | 48 Cd   | 48 Cd   | 48 Cd   | 48 Cd       | 48 Cd             | 48 Cd   | 48 Cd       | 48 Cd        | 48 Cd        | 48 Cd   | 48 Cd       | 48 Cd       |
| 〈 κ | r - 1 | ¯ κ 〉 | 1S      | 2S      | 3S      | 4S          | 5S                |         | 2P          | 3P           | 4P           |         | 3D          | 4D          |
| 1S                    | 47.48   | 9.194   | -3.856  | -1.617      | 0.3840            |         |             |              |              |         |             |             |
| 2S                    |         | 10.84   | -2.575  | -1.030      | 0.2429            | 2P      | 10.81       | -2.174       | -0.820       |         |             |             |
| 3S                    |         |         | 3.895   | 0.851       | -0.1933           | 3P      |             | 3.774        | 0.756        | 3D      | 3.637       | -0.532      |
| 4S                    |         |         |         | 1.485       | -0.1793           | 4P      |             |              | 1.369        | 4D      |             | 1.063       |
| 5S                    |         |         |         |             | 0.3905            |         |             |              |              |         |             |             |
| [35]                  | 50.64   | 11.72   | 4.127   | 1.554       | 0.4128            |         | 10.93 11.66 | 3.850 4.024  | 1.380 1.432  |         | 3.630 3.676 | 1.049 1.066 |
| 〈 κ | r - 2 | ¯ κ 〉 | 1S      | 2S      | 3S      | 4S          | 5S                |         | 2P          | 3P           | 4P           |         | 3D          | 4D          |
| 1S                    | 4524.   | 1302.   | -561.4  | -236.5      | 56.21             |         |             |              |              |         |             |             |
| 2S                    |         | 477.6   | -195.6  | -81.99      | 19.47             | 2P      | 158.4       | -51.48 29.99 | -20.11 10.44 | 3D      | 16.62       | -3.886      |
| 3S 4S                 |         |         | 92.91   | 37.55 17.16 | -8.893 -3.939     | 3P 4P   |             |              | 5.374        | 4D      |             | 2.098       |
| 5S                    |         |         |         |             | 1.071             |         |             |              |              |         |             |             |
| [35]                  | 5504.   | 672.8   | 121.9   | 22.50       | 1.490             |         | 162.7       | 31.80        | 5.571        |         | 16.55       | 2.058       |
|                       |         |         |         |             |                   |         | 201.4       | 38.96        | 6.775        |         | 17.05       | 2.142       |

third terms are calculated using Eq. (II.22), the second term is given by Eq. (II.20). We remark that the third term is cancelled by the first one on the right side of Eq. (II.20). Finally, the variance of the excitation energy of the electron shell turns out to be

$$\frac { 1 } { 4 } \mathcal { D } = \sum _ { \kappa } \mathcal { N } _ { \kappa } \langle \kappa | r ^ { - 2 } | \kappa \rangle - \sum _ { \kappa \bar { \kappa } } \min ( \mathcal { N } _ { \kappa }, \mathcal { N } _$$

The sum of the diagonal components coincides with the right side of Eq. (II.13). The off-diagonal components in the second term are associated with exchange effects. The exchange effects reduce the variance.

## III. NUMERICAL RESULTS

To calculate the variance with exchange effects, Eq. (II.23) requires knowledge of the matrix elements of r -1 and r -2 . The average values of r -1 in the inner and outer electron orbits are approximately in the ratio of Z : 1. The average values of r -2 correlate roughly as Z 2 : 1. When an atom has multiple partially filled shells, Eq. (II.23) can be used to provide an approximate estimate. Electrons on the inner completely filled shells contribute the most to the variance, hence the accuracy of the estimate is expected to be reasonable.

Wave functions of stationary states are determined up to a phase multiplier and usually selected as real. The offdiagonal matrix elements are then defined up to the sign. The signs of radial wave functions differ in the RHF method and in the program package G rasp 2018. The variance depends on squares of the off-diagonal matrix elements, so the sign arbitrariness has no effect on the results.

TABLE V. Matrix elements 〈 κ | r -1 | ¯ κ 〉 and 〈 κ | r -2 | ¯ κ 〉 for the orbital states of electrons in tellurium and xenon atoms. Notations are identical to those in Table II.

| 52 Te                 | 52 Te   | 52 Te   | 52 Te   | 52 Te        | 52 Te        | 52 Te   | 52 Te       | 52 Te       | 52 Te        | 52 Te        | 52 Te   | 52 Te       | 52 Te       |
|-----------------------|---------|---------|---------|--------------|--------------|---------|-------------|-------------|--------------|--------------|---------|-------------|-------------|
| 〈 κ | r - 1 | ¯ κ 〉 | 1S      | 2S      | 3S      | 4S           | 5S           |         | 2P          | 3P          | 4P           | 5P           |         | 3D          | 4D          |
| 1S                    | 51.47   | -10.02  | -4.273  | 1.889        | 0.639        |         |             |             |              |              |         |             |             |
| 2S                    |         | 11.82   | 2.851   | -1.201       | -0.403       | 2P      | 11.80       | -2.463      | -0.971       | 0.287        |         |             |             |
| 3S                    |         |         | 4.316   | -0.991       | -0.318       | 3P      |             | 4.230       | 0.902        | -0.256       | 3D      | 4.084       | 0.684       |
| 4S                    |         |         |         | 1.722        | 0.290        | 4P      |             |             | 1.616        | -0.244       | 4D      |             | 1.369       |
| 5S                    |         |         |         |              | 0.590        | 5P      |             |             |              | 0.504        |         |             |             |
| [35]                  | 55.57   | 12.96   | 4.627   | 1.815        | 0.602        |         | 11.95 12.91 | 4.294 4.526 | 1.631 1.700  | 0.470 0.493  |         | 4.073 4.133 | 1.354 1.374 |
| 〈 κ | r - 2 | ¯ κ 〉 | 1S      | 2S      | 3S      | 4S           | 5S           |         | 2P          | 3P          | 4P           | 5P           |         | 3D          | 4D          |
| 1S                    | 5316.   | -1539.  | -674.9  | 300.0        | 101.6        |         |             |             |              |              |         |             |             |
| 2S                    |         | 567.1   | 236.4   | -104.5       | -35.35       | 2P      | 188.4       | -63.40      | -26.07       | 7.759        |         |             |             |
| 3S                    |         |         | 114.1   | -48.63       | -16.40       | 3P      |             | 37.75       | 13.91        | -4.107       | 3D      | 20.80       | 5.625       |
| 4S                    |         |         |         | 23.35        | 7.628        | 4P      |             |             | 7.485        | -2.026       | 4D      |             | 3.401       |
| 5S                    |         |         |         |              | 2.841        | 5P      |             |             |              | 0.806        |         |             |             |
| [35]                  | 6713.   | 785.3   | 157.7   | 32.09        | 3.731        |         | 194.6       | 39.51       | 7.809        | 0.719        |         | 20.73       | 3.347       |
| 54 Xe                 | 54 Xe   | 54 Xe   | 54 Xe   | 54 Xe        | 54 Xe        | 54 Xe   | 54 Xe       | 54 Xe       | 54 Xe        | 54 Xe        | 54 Xe   | 54 Xe       | 54 Xe       |
| 〈 κ | r - 1 | ¯ κ 〉 | 1S      | 2S      | 3S      | 4S           | 5S           |         | 2P          | 3P          | 4P           | 5P           |         | 3D          | 4D          |
| 1S                    | 53.47   | 10.43   | -4.482  | 2.027        | 0.712        |         |             |             |              |              |         |             |             |
| 2S                    |         | 12.31   | -2.990  | 1.287        | 0.448        | 2P      | 12.29       | 2.584       | -1.047       | -0.317       |         |             |             |
| 3S                    |         |         | 4.527   | -1.063       | -0.353       | 3P      |             | 4.445       | -0.972       | -0.282       | 3D      | 4.304       | 0.751       |
| 4S                    |         |         |         | 1.843        | 0.321        | 4P      |             |             | 1.741        | 0.267        | 4D      |             | 1.509       |
| 5S                    |         |         |         |              | 0.648        | 5P      |             |             |              | 0.547        |         |             |             |
| [35]                  | 58.11   | 13.61   | 4.884   | 1.950        | 0.681        |         | 12.47 13.55 | 4.519 4.784 | 1.760        | 0.546 0.576  |         | 4.293 4.362 | 1.495 1.517 |
| 〈 κ | r - 2 | ¯ κ 〉 | 1S      | 2S      | 3S      | 4S           | 5S           |         | 2P          | 3P          | 1.841 4P     | 5P           |         | 3D          | 4D          |
| 1S                    | 5735.   | 1665.   | -735.7  | 334.7        | 117.6        |         |             |             |              |              |         |             |             |
| 2S                    |         | 614.8   | -258.2  | 116.8        | 41.01        | 2P      | 204.4       | 69.35       | -29.34       | -8.935       |         |             |             |
| 3S 4S                 |         |         | 125.5   | -54.76 26.85 | -19.16 9.107 | 3P 4P   |             | 41.60       | -15.78 8.687 | -4.761 2.405 | 3D 4D   | 23.06       | 6.526 4.097 |
| 5S                    |         |         |         |              | 3.506        | 5P      |             |             |              | 0.970        |         |             |             |
| [35]                  | 7391.   | 879.8   | 178.5   | 37.89        | 5.020        |         | 211.7 279.6 | 43.73 37.05 | 9.104 11.72  | 0.998 1.305  |         | 22.99 23.89 | 4.048 4.212 |

The RHF method is an excellent tool for solving problems involving the construction of atomic shells. Reference [15] presents a systematic collection of the atomic wave functions through decomposition over a finite set of basis functions. The number of terms in the finite series and the parameters of basis functions are determined independently for each shell. The present paper uses these techniques to evaluate the matrix elements in atoms of the interest. The calculation results are placed in Tables II - V, where calculations of diagonal matrix elements are also compared with the results of the relativistic DHF method [35]. Table I displays the exchange corrections and variance. The exchange effect contributions to D can be estimated at about 10%.

In all examples studied, the square of variance is typically an order of magnitude greater than the average value. This is likely related to the fact that in the doubleβ decay, the predominant contribution to the variance is made by rare excitations of electrons in the inner orbits. Feinberg [38] and Migdal [39] calculate the probability of atom ionization due to shaking caused by knocking an electron from the K shell into the continuum.

The non-relativistic theory has a degeneracy in the total angular momentum j = l ± 1 / 2. For testing purposes, we evaluated the orthonormality of the wave functions of the orbitals against the radial quantum number for each pair ( jl ) and compared them to the known diagonal matrix elements from Ref. [35]. In the relativistic models, the wave functions of orbitals j = l ± 1 / 2 split, resulting in two values for each l ≥ 1. The upper row represents j = l +1 2. / The orthonormality holds with an accuracy of at least 10 -5 . The agreement between the diagonal matrix elements is satisfactory. Deviations increase systematically with increasing Z , owing to relativistic effects. The non-relativistic approach underestimates negative moments of radius because in the relativistic framework the probability density of atomic electrons is shifted to shorter distances. Since the off-diagonal matrix elements are symmetric under the interchange n ↔ ¯, the tables display only the values n n ≤ ¯. n Table I contains the variance data for atoms with atomic numbers Z ≤ 54. The constraint on Z is due to the limitations of the non-relativistic approximation.

TABLE VI. Matrix elements 〈 njl r | -1 | ¯ njl 〉 and 〈 njl r | -2 | ¯ njl 〉 for orbital states of electrons in neodymium and uranium atoms, obtained by the DHF method with the use of the program package G rasp 2018 [17]. Each matrix element except for l = 0 is given two values of the total angular momentum j . The upper and lower lines correspond to j = l ± 1 2. /

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00068v1-figure-c1060097-cf7a-4886-920a-f4ada48cd78d.jpg)

| - 1                                                      | 1S        | 2S              | 3S                | 4S                      | 5S                            | 6S                                  |                   |                | 60 Nd          | 4P                                     |                                                    |                   |                               | 3D                   | 4D                      | 5D          | 6D                   |             | 4F          |                   |
|----------------------------------------------------------|-----------|-----------------|-------------------|-------------------------|-------------------------------|-------------------------------------|-------------------|----------------|----------------|----------------------------------------|----------------------------------------------------|-------------------|-------------------------------|----------------------|-------------------------|-------------|----------------------|-------------|-------------|-------------------|
| 〈 κ | r | ¯ κ 〉 1S 2S 3.524 3S 5.202 4S 5.590 5S 6S 7S | 66.04     | 14.07 15.65     | 6.167 4.074 5.689 | 2.894 1.825 1.439 2.336 | 1.132 0.707 0.533 0.464 0.874 | 0.299 0.186 0.139 0.117 0.112 0.258 | 7S                | 2P 3P 4P 5P 6P | 2P 14.02 15.59 | 3P 3.055 1.497 1.198 1.322 2.105 2.227 | 5P 1.289 0.443 0.393 0.447 0.368 0.405 0.725 0.773 | 6P 0.529          | 3D                            | 4.959 5.055 4D 5D 6D | 0.914 0.943 1.845 1.877 |             | 4F                   | 5F          | 1.291 1.306 | 5F                |
| 92 U                                                     |           |                 |                   |                         |                               |                                     |                   |                |                |                                        |                                                    |                   |                               |                      |                         |             |                      |             |             |                   |
| 〈 κ | r - 1 | ¯ κ 〉                                    | 1S        | 2S              | 3S                | 4S                      | 5S                            | 6S                                  | 7S                |                | 2P             | 3P 4P                                  | 5P                                                 | 6P                |                               | 3D                   | 4D                      | 5D          |                      | 6D          | 4F          | 5F                |
| 1S 2S                                                    | 122.4     | 30.38 30.96     | 14.02 9.205       | 7.179 4.519             | 3.541                         | 1.565 0.973                         | 0.464             | 2P             | 22.78          | 5.335                                  | 2.478 1.150 3.801                                  | 0.445 1.785 0.734 | 3D                            | 8.493 8.924          | 1.778 1.909             | 0.686 0.741 | 0.189 0.204          | 4F          | 3.326 3.384 | 0.419             |
| 3S                                                       |           |                 | 11.28             | 3.348                   | 2.207 1.562                   | 0.682                               | 0.288 0.202       | 3P             | 30.88          | 8.048 9.042 3.130                      | 2.338 1.028 1.401                                  | 0.394             | 4D 0.571                      |                      | 3.663 3.814             | 0.748 0.791 | 0.199 0.210          | 5F          |             | 0.428 0.917 0.929 |
| 4S 5S 6S 7S                                              |           |                 |                   | 4.820                   | 1.273 2.110                   | 0.531 0.459 0.130 0.883             | 0.156 0.125 0.283 | 4P 5P 6P       |                | 11.17                                  | 4.000 0.947 4.709 1.186 1.745 1.999                |                   | 0.346 0.463 0.326 0.407 0.675 | 5D 6D                | 1.519                   | 1.467       | 0.202 0.209 6D 3.519 | 0.435 0.444 |             |                   |
| 〈 κ | r - 2 | ¯ κ 〉 3P 1S 273.7                        | 1S 40946. | 2S 15344. 6502. | 3S 7382. 3072.    | 4S 3819. 1586.          | 5S 1889. 783.9                | 6S 835.5 346.7                      | 7S 247.6 102.7    | 2P             | 2P 712.5 2077. | 4P 915.8                               | 5P 134.7 63.35 462.1 220.4 79.58 37.05             |                   | 0.787 6P 24.61                |                      | 4D                      | 5D          |                      | 4F          |             | 2.803             |
| 177.6 4S 5S                                              |           |                 | 1536.             | 782.5 415.2             | 386.3 202.7                   | 170.8 89.52                         | 50.62 26.53       | 3P             |                | 508.4                                  | 246.3 116.9                                        | 48.24             | 8.119                         |                      | 26.45                   | 9.512       | 2.618                |             |             |                   |
|                                                          |           |                 |                   |                         |                               | 44.82 20.31                         | 13.27 5.951       | 4P             |                |                                        | 48.80 21.09                                        | 25.76             |                               | 5D                   |                         | 4.921       | 1.365                |             |             |                   |
|                                                          |           |                 |                   |                         |                               |                                     |                   | 5P             |                |                                        | 136.2                                              | 4.225             | 6D                            |                      |                         |             |                      |             |             |                   |
| 7S                                                       |           |                 |                   |                         |                               |                                     | 1.829             |                |                |                                        | 11.74                                              |                   |                               |                      |                         |             |                      |             |             |                   |
|                                                          |           |                 |                   |                         |                               |                                     |                   |                |                |                                        | 32.11                                              | 12.80             |                               |                      |                         |             |                      |             |             |                   |
|                                                          |           |                 |                   |                         |                               |                                     |                   | 6P             |                |                                        |                                                    | 1.962             |                               |                      |                         |             |                      |             |             |                   |
|                                                          |           |                 |                   |                         |                               |                                     |                   |                |                |                                        |                                                    | 5.693             |                               |                      |                         |             |                      |             |             |                   |
|                                                          |           |                 |                   |                         |                               |                                     |                   |                |                |                                        |                                                    |                   |                               | 88.84                |                         |             |                      |             |             |                   |
|                                                          |           |                 |                   |                         |                               |                                     |                   |                |                |                                        |                                                    |                   |                               | 3D                   | 31.07                   | 12.61       |                      |             | 13.36       |                   |
| 2S                                                       |           |                 |                   |                         |                               |                                     |                   |                |                |                                        |                                                    |                   | 3D 4D                         | 100.4                | 35.99                   | 14.72       | 4.088                |             | 4F 5F       | 2.685             |
| 3S                                                       |           |                 |                   |                         | 102.5                         |                                     |                   |                |                |                                        |                                                    | 62.64             | 90.99 14.36                   |                      |                         | 8.298       |                      |             | 13.89       |                   |
| 6S                                                       |           |                 |                   |                         |                               |                                     |                   |                |                |                                        |                                                    |                   |                               |                      | 23.45                   |             | 2.290                | 5F          | 1.465       |                   |
|                                                          |           |                 |                   |                         |                               |                                     |                   |                |                |                                        |                                                    |                   |                               |                      |                         | 5.534       | 1.209                |             | 1.517       |                   |
|                                                          |           |                 |                   |                         |                               |                                     |                   |                |                |                                        |                                                    |                   |                               |                      |                         |             |                      | 0.497       |             |                   |
|                                                          |           |                 |                   |                         |                               |                                     |                   |                |                |                                        |                                                    |                   |                               |                      |                         |             |                      | 0.546       |             |                   |

For the relativistic calculation of the quantities entering Eqs. (II.10) and (II.23), the package G rasp 2018 is used. G rasp 2018 is based on the multi-configuration DHF approach to solve the most significant atomic physics problems: identifying energy levels, calculation of radiation decay rates, hyperfine structure, etc. The package allows to compute radial wave functions for any given electron configuration with definite total angular momentum. For each atom from Table I we compute a set of orbital wave functions f α ( r ) and g α ( r ), total energies, and matrix elements 〈 κ | r -1 | ¯ κ 〉 and 〈 κ | r -2 | ¯ κ 〉 needed to calculate the parameters C and D . The values of C are in reasonable agreement with the results of Ref. [14]. The results for ∆ D 1 / 2 DHF and D 1 / 2 DHF are compared with the non-relativistic counterparts ∆ D 1 / 2 RHF and D 1 / 2 RHF . In Table VI we provide the matrix elements for neodymium and uranium atoms obtained with the DHF

method. In this case, the splitting of the levels with j = +1 2 and l / j = l -1 / 2 is significant and relativistic calculation is mandatory.

Knowledge of the parameters C and D allows to model the probability density of the kinetic energy of β electrons. A popular phenomenological distribution in the finite interval [0 , 1] is the beta-distribution [40]. We use the betadistribution to describe the probability density in the process that occurs with the excitation of the shell. The full distribution also takes into account the singular part associated with the probability of the shell remaining in an unexcited state [14]:

$$w ( x ) = K _ { Z } ^ { 2 } \delta ( x ) + ( 1 - K _ { Z } ^ { 2 } ) \frac { \Gamma ( a ) \Gamma ( b ) } { \Gamma ( a + b ) } x ^ { a - 1 } ( 1 - x ) ^ { b - 1 }. & & ( \text{III.1} )$$

The probability is normalized by one, with a &gt; 0 and b &gt; 0. In the given instance, x = /epsilon1/Q ∗ , where /epsilon1 is the excitation energy of the electron shell; Q ∗ = Q -I 2 , where I 2 is the sum of the first two ionization potentials of the daughter atom; K Z is the overlap integral of the electron shells of the parent and daughter atoms. The kinetic energy of β electrons is the difference between the magnitude of Q ∗ and the excitation energy of the electron shell. Since the average excitation energy and variance are known, the parameters a and b can be determined from the equations

$$1$$

$$\int _ { 0 } ^ { 1 } x w ( x ) d x \ & = \ \frac { \mathcal { C } } { Q ^ { * } }, & ( \text{III.2} ) \\ \int _ { 0 } ^ { 1 } x ^ { 2 } w ( x ) d x \ & = \ \frac { \mathcal { D } + \mathcal { C } ^ { 2 } } { Q ^ { * 2 } }. & ( \text{III.3} ) \\ \text{ribution parameters allows for the determination of probability distribution of the}$$

This method for evaluating distribution parameters allows for the determination of probability distribution of the excitation energy. In the decay of germanium, e.g., the excitation energy does not exceed 0 55 keV with 95% proba-. bility.

## IV. CONCLUSION

This work finds the contribution of exchange effects to the variance of the excitation energy of the shell of a daughter atom in 0 ν 2 β -decay. The effect is described by Eq. (II.23), which is derived in Sect. II.B. To compute the variance, one needs to know the matrix elements of the values r -1 and r -2 . These matrix elements were determined for eleven atoms that exhibit 2 ν 2 β -decay. Tables I-VI provide a summary of the results. Table I shows the values of the corrections for exchange effects. The calculations utilized the non-relativistic RHF approach, which is justified for Z /lessorsimilar 54. The results were compared to those obtained with the software package G rasp 2018, which uses the DHF relativistic technique. Both methods are in reasonable agreement at Z /lessorsimilar 54. The relativistic DHF approach is also used to evaluate the average excitation energy of atoms listed in Table I and, in particular, the impact of exchange effects in heavy neodymium and uranium atoms.

The obtained results can be utilized to describe the spectrum of β electrons in 0 ν 2 β -decay, taking into consideration the excitation of the electron shell of the daughter atom.

## ACKNOWLEDGEMENTS

The authors gratefully acknowledge F.F. Karpehsin for discussions.

- [1] S. Bilenky, Introduction to the Physics of Massive and Mixed Neutrinos , 2nd ed., Lecture Notes in Physics, Springer-Verlag, Berlin (2018), v. 947.
- [2] J. Schechter and J. W. F. Valle, Phys. Rev. D 25 , 2951 (1982).
- [3] M. Hirsch, S. Kovalenko, and I. Schmidt, Phys. Lett. B 642 , 106 (2006).
- [4] K. Blaum, S. Eliseev, F. A. Danevich, V. I. Tretyak, Sergey Kovalenko, M. I. Krivoruchenko, Yu. N. Novikov, and J. Suhonen, Rev. Mod. Phys. 92 , 045007 (2020).
- [5] F. ˇ imkovic, Phys.-Uspekhi S 64 , 1238 (2021).
- [6] D.Q. Adams, C. Alduino, K. Alfonso et al. (CUORE Collaboration), Phys. Rev. Lett. 124 , 122501 (2020).
- [7] D.Q. Adams, C. Alduino, K. Alfonso et al. (CUORE Collaboration), Nature 604 , 53 (2022).
- [8] G. Anton, I. Badhrees, P. S. Barbeau et al. (EXO-200 Collaboration), Phys. Rev. Lett. 123 , 161802 (2019).

- [9] S. Abe, S. Asami, M. Eizuka et al. (KamLAND-Zen Collaboration), Phys. Rev. Lett. 130 , 051801 (2023).
- [10] LEGEND collaboration, The Large Enriched Germanium Experiment for Neutrinoless ββ Decay: LEGEND-1000 Preconceptual Design Report , arXiv:2107.11462 (2021).
- [11] R. Arnold, C. Augier, J.D. Baker et al. (NEMO-3 Collaboration), Phys. Rev. D 89, 111101(R) (2014).
- [12] R. Arnold, C. Augier, J.D. Baker et al. (NEMO-3 Collaboration), Phys. Rev. D 92, 072011 (2015).
- [13] M.I. Krivoruchenko, K.S. Tyrin, and F.F. Karpeshin, JETP Lett. 117 , 884 (2023).
- [14] M.I. Krivoruchenko, K.S. Tyrin, and F.F. Karpeshin, JETP Lett. 118 , 470 (2023).
- [15] E. Clementi and C. Roetti, Atomic Data and Nuclear Data Tables 14 , 177-478 (1974).
- [16] K.G. Dyall, I.P. Grant, C.T. Johnson, F.A. Parpia, and E.P. Blummer, Comput. Phys. Commun. 55 , 425-456 (1989).
- [17] C.F. Fischer, G. Gaigalas, P. J¨nsson, J. Biero´, Comp. Phys. Commun. o n 237 , 184 (2019).
- [18] I.P. Grant, Relativistic Quantum Theory of Atoms and Molecules: Theory and Computation , (Springer Science+Business Media, New York 2007).
- [19] A.A. Kwiatkowski, T. Brunner, J.D. Holt, A. Chaudhuri, U. Chowdhury, M. Eibach, J. Engel, A.T. Gallant, A. Grossheim, M. Horoi, A. Lennarz, T.D. Macdonald, M.R. Pearson, B.E. Schultz, M.C. Simon, R.A. Senkov, V.V. Simon, K. Zuber, J. Dilling, Phys. Rev. C 89 , 045502 (2014).
- [20] The GERDA Collaboration, Background-free search for neutrinoless doubleβ decay of 76Ge with GERDA , Nature 544 , 47 - 52 (2017).
- [21] M. Suhonen, I. Bergstr¨m, T. Fritioff, Sz. Nagy, A. Solders and R. Schuch, J. Instrum. o 2 06003, (2007).
- [22] S. Rahaman, V.V. Elomaa, T. Eronen, J. Hakala, A. Jokinen, J. Julin, A. Kankainen, A. Saastamoinen, J. Suhonen, C. Weber, J. Aysto, Phys. Lett. B 662 , 111-116 (2008)
- [23] D.L. Lincoln, J.D. Holt, G. Bollen, M. Brodeur, S. Bustabad, J. Engel, S.J. Novario, M. Redshaw, R. Ringle, and S. Schwarz, Phys. Rev. Lett. 110 , 012501 (2013).
- [24] M. Alanssari, D. Frekers, T. Eronen, L. Canete, J. Dilling, M. Haaranen, J. Hakala, M. Holl, M. Jeˇkovsk´, A. Jokinen, A. s y Kankainen, J. Koponen, A.J. Mayer, I. Moore, D. Nesterenko, I. Pohjalainen, P. Povinec, J. Reinikainen, S. Rinta-Antila, P.C. Srivastava, J. Suhonen, R.I. Thompson, A. Voss, M.E. Wieser, Phys. Rev. Lett. 116 , 072501 (2016).
- [25] S. Rahaman, V.-V. Elomaa, T. Eronen, J. Hakala, A. Jokinen, A. Kankainen, J. Rissanen, J. Suhonen, C. Weber, J. Aysto, Phys. Lett. B 703 , 412-416 (2011).
- [26] N.D. Scielzo, S. Caldwell, G. Savard, J.A. Clark, C.M. Deibel, J. Fallis, S. Gulick, D. Lascar, A.F. Levand, G. Li, J. Mintz, E.B. Norman, K.S. Sharma, M. Sternberg, T. Sun, J. Van Schelt, Phys. Rev. C 80 , 025501 (2009).
- [27] M. Redshaw, E. Wingfield, J. McDaniel, and E.G. Myers, Phys. Rev. Lett. 98 , 053003 (2007).
- [28] V.S. Kolhinen, T. Eronen, D. Gorelov, J. Hakala, A. Jokinen, A. Kankainen, I.D. Moore, J. Rissanen, A. Saastamoinen, J. Suhonen, and J. Ayst¨, Phys. Rev. C ¨ o 82 , 022501 (2010).
- [29] R.B. Firestone, V.S. Shirley, editors; C.M. Baglin, S.Y. Frank Chu, and J. Zipkin, assistant editors. Table of isotopes (Wiley-Interscience, New York 1996)
- [30] L.D. Landau, E.M. Lifschitz, Quantum Mechanics: Non-relativistic Theory. Course of Theoretical Physics. Vol. 3 (3rd ed.), (London: Pergamon1977).
- [31] J.B. Gaunt, The Triplets of Helium , Proc. R. Soc. A 122 , 513 (1929).
- [32] G. Breit, The Effect of Retardation on the Interaction of Two Electrons , Phys. Rev. 34 , 553 (1929).
- [33] M.I. Krivoruchenko, Phys. Atom. Nucl. 86 , 153 (2023).
- [34] C.C. Lu, T.A. Carlson, F.B. Malik, T.C. Tucker, and C.W. Nestor, Jr., At. Data Nucl. Data Tables 3 , 1 (1971).
- [35] J.P. Desclaux, Atomic Data and Nuclear Data Tables, 12 , 31l-406 (1973).
- [36] Ken-Ning Huang, Michio Aoyagi, Mau Hsiung Chen, Bern Grasemann, Hans Mark, Atomic Data and Nuclear Data Tables 18 , 243-29 (1976).
- [37] D.A. Varshalovich, A.N. Moskalev, and V.K. Khersonskii, Quantum Theory of Angular Momentum , (World Scientific Publishing, Singapore 1988).
- [38] E.L. Feinberg, J. Phys. (USSR) 4 , 423 (1941).
- [39] A. Migdal, J. Phys. Acad. Sci. USSR 4 , 449 (1941).
- [40] V.S. Korolyuk, N.I. Portenko, A.V. Skorokhod, and A.F. Turbin, Handbook of Probability Theory and Mathematical Statistics (Nauka, Moscow, 1985) [in Russian].
K. S. Tyrin M. I. Krivoruchenko
Exchange contribution to the variance of the excitation energy of the electron shell of the daughter atom in double-$β$ decay
The excitation of electron shell of a daughter atom in a neutrinoless double
$\beta$-decay causes change in shape of the total energy peak of $\beta$
electrons at the end of the energy spectrum. The main parameters of the
modified distribution are the average energy and variance of the excitation
energy of the electron shell. We derive an expression for the variance taking
into account exchange effects and make numerical estimates of the average
excitation energy and variance based on the non-relativistic
Roothaan-Hartree-Fock method and the relativistic Dirac-Hartree-Fock method
implemented in the General Relativistic Atomic Structure Package (Grasp2018).
Estimates are made for eleven isotopes, two-neutrino double-$\beta$ decay of
which is observed experimentally. The results are determined by the first two
negative moments of the electron radii in the parent atom. The values obtained
for change in the peak shape can be used to parameterize the energy
distribution of $\beta$ electrons, taking into account the excitation of the
electron shell of the daughter atom.