Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential

We prove that a neutral atom in mean-field approximation has O(4) symmetry and this fact explains the empirical [n+l,n]-rule or Madelung rule which describes effectively periods, structure and other properties of the Mendeleev table of chemical elements.

Gespeichert in:
Bibliographische Detailangaben
Datum:2017
1. Verfasser: Belokolos, E.D.
Format: Artikel
Sprache:English
Veröffentlicht: Інститут математики НАН України 2017
Schriftenreihe:Symmetry, Integrability and Geometry: Methods and Applications
Online Zugang:http://dspace.nbuv.gov.ua/handle/123456789/148640
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential / E.D. Belokolos // Symmetry, Integrability and Geometry: Methods and Applications. — 2017. — Т. 13. — Бібліогр.: 34 назв. — англ.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-148640
record_format dspace
spelling irk-123456789-1486402019-02-19T01:31:48Z Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential Belokolos, E.D. We prove that a neutral atom in mean-field approximation has O(4) symmetry and this fact explains the empirical [n+l,n]-rule or Madelung rule which describes effectively periods, structure and other properties of the Mendeleev table of chemical elements. 2017 Article Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential / E.D. Belokolos // Symmetry, Integrability and Geometry: Methods and Applications. — 2017. — Т. 13. — Бібліогр.: 34 назв. — англ. 1815-0659 2010 Mathematics Subject Classification: 81Q05; 81V45 DOI:10.3842/SIGMA.2017.038 http://dspace.nbuv.gov.ua/handle/123456789/148640 en Symmetry, Integrability and Geometry: Methods and Applications Інститут математики НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
description We prove that a neutral atom in mean-field approximation has O(4) symmetry and this fact explains the empirical [n+l,n]-rule or Madelung rule which describes effectively periods, structure and other properties of the Mendeleev table of chemical elements.
format Article
author Belokolos, E.D.
spellingShingle Belokolos, E.D.
Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential
Symmetry, Integrability and Geometry: Methods and Applications
author_facet Belokolos, E.D.
author_sort Belokolos, E.D.
title Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential
title_short Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential
title_full Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential
title_fullStr Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential
title_full_unstemmed Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential
title_sort mendeleev table: a proof of madelung rule and atomic tietz potential
publisher Інститут математики НАН України
publishDate 2017
url http://dspace.nbuv.gov.ua/handle/123456789/148640
citation_txt Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential / E.D. Belokolos // Symmetry, Integrability and Geometry: Methods and Applications. — 2017. — Т. 13. — Бібліогр.: 34 назв. — англ.
series Symmetry, Integrability and Geometry: Methods and Applications
work_keys_str_mv AT belokolosed mendeleevtableaproofofmadelungruleandatomictietzpotential
first_indexed 2025-07-12T19:51:17Z
last_indexed 2025-07-12T19:51:17Z
_version_ 1837472034938945536
fulltext Symmetry, Integrability and Geometry: Methods and Applications SIGMA 13 (2017), 038, 15 pages Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential Eugene D. BELOKOLOS Department of Theoretical Physics, Institute of Magnetism, National Academy of Sciences of Ukraine, 36-b Vernadsky Blvd., Kyiv, 252142, Ukraine E-mail: bel@imag.kiev.ua Received February 27, 2017, in final form May 22, 2017; Published online June 07, 2017 https://doi.org/10.3842/SIGMA.2017.038 Abstract. We prove that a neutral atom in mean-field approximation has O(4) symmetry and this fact explains the empirical [n+l, n]-rule or Madelung rule which describes effectively periods, structure and other properties of the Mendeleev table of chemical elements. Key words: Madelung rule; Mendeleev periodic system of elements; Tietz potential 2010 Mathematics Subject Classification: 81Q05; 81V45 1 Introduction In 1869 D.I. Mendeleev discovered a periodic dependence of chemical element properties on Z with periods 2, 8, 8, 18, 18, 32, . . . , (1.1) where Z is a number of electrons in the chemical element atom1. With creation of quantum mechanics physicists tried to explain the Mendeleev empirical law in terms of the one-particle quantum numbers n, l, m, σ, where n = nr+l+1 is a principal quantum number, nr is a radial quantum number, l is an orbital quantum number, m is a magnetic quantum number, −l ≤ m ≤ l, σ = ±1/2 is a projection of electron spin. For example, energy levels of the Hydrogen atom, which has O(4) symmetry, depend on the principal quantum number n only and are degenerate in other quantum numbers: l, 0 ≤ l ≤ n − 1 (so called “accidental” degeneracy); m, −l ≤ m ≤ l; σ, σ = ±1/2. It is a so called [n, l]-rule. According to it the energy spectrum of atom consists of electron shells, enumerated by the principal quantum number n and having the degeneracy Nn = n∑ l=0 2(2l + 1) = 2n2. This formula gives the following set of periods 2, 8, 18, 32, . . . . (1.2) Comparing the period sequences (1.1) and (1.2), we can explain only the first two periods of the Mendeleev table. Although we see a remarkable fact that (1) the period lengths have cardi- nalities that correspond to the Hydrogen degeneracy dimensions, and (2) the same cardinalities always occur in pairs, except for the very first one. In this paper we show that these similarities are not accidental. 1See https://en.wikipedia.org/wiki/Periodic_table. mailto:bel@imag.kiev.ua https://doi.org/10.3842/SIGMA.2017.038 https://en.wikipedia.org/wiki/Periodic_table 2 E.D. Belokolos 2 The Madelung [n + l, n]-rule We get an exact expression for the Mendeleev periods and other properties of the Mendeleev table with the Madelung [n+ l, n]-rule [22]). The [n + l, n]-rule asserts: with growth of atomic charge Z the electrons fill up in atom consecutively the one-particle states with the least possible value of the quantum number n+ l; and, for a given value n + l, the electrons fill up states with the least possible value of the quantum number n. Here it is reasonable to introduce the quantum number M = n+ l = nr + 2l + 1. The [n+l, n]-rule is in fact an algorithm for consecutive building-up of atoms2. It describes the states for each of about 5000 electrons of the atoms to corresponding elements of the periodic system (indeed, now we have 118 chemical elements, and for approximately 100 elements we know every state of electron configuration, that is 100∑ 1 Z = 5050 states) and predicts correctly the real electron configurations of all elements of periodic system with small number of exclusions. There exist only 19 elements (Cr, Cu, Nb, Mo, Ru, Rh, Pd, Ag, La, Ce, Gd, Pt, Au, Ac, Th, Pa, U, Np, Cm), whose electron configurations differ from the configurations predicted with the [n+ l, n]-rule [24]. The [n+ l, n]-rule allows to obtain relations between the order number Z of chemical element and quantum numbers M , n, l of the appropriate electron configuration. Let us designate the order number Z of chemical element in which electrons from the (n+ l)- subgroups, (n, l)-subgroups, n-subgroups, l-subgroups or nr-subgroups appear for the first time by symbols Zn+l, Zn,l, Zn, Zl, Znr . Let us designate an order number Z of a chemical element in which the (n+ l)-subgroups, (n, l)- subgroups, n-subgroups, l-subgroups or nr-subgroups are filled up completely for the first time by symbols n+lZ, n,lZ, nZ, lZ, nrZ. Then the [n+ l, n]-rule with help of four actions of arithmetic and well known formulas k=n∑ k=1 k = (1/2)n(n+ 1), k=n∑ k=1 k2 = (1/6)n(n+ 1)(2n+ 1), leads to the Klechkovski–Hakala formulas [11, 15] Zn+l = K(n+ l) + 1, n+lZ = K(n+ l + 1), Zn,l = K(n+ l + 1)− 2(l + 1)2 + 1, n,lZ = K(n+ l + 1)− 2l2, Zn = K(n+ 1)− 1, nZ = K(2n)− 2(n− 1)2 = (1/6) [ (2n− 1)3 + 11(2n− 1) ] , Zl = K(2l + 1) + 1 = (1/6) [ (2l + 1)3 + (5− 2l) ] , Znr = K(nr + 2)− 1. (2.1) 2See https://en.wikipedia.org/wiki/Aufbau_principle. https://en.wikipedia.org/wiki/Aufbau_principle Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential 3 Here K(x) is the following function K(x) = (1/6)x [ x2 + 2− 3µ(x) ] , x ∈ N, µ(x) = x mod (2) = { 1, x is odd, 0, x is even. Relations for Zn+l and n+lZ are exact. Other relations are exact up to exclusions pointed out above. For example, for l = 0, 1, . . ., we have Zl = 1, 5, 21, 57, . . ., that is, according to the Madelung rule, the onset of the 4f block starts with La (Z = 57). But, according to the IUPAC data3 on electron configurations in atoms, the onset of the 4f block starts with Ce (Z = 58). The [n+ l, l]-rule defines essential characteristics of the Mendeleev periodic system. In the Mendeleev periodic system, the table rows are enumerated by the Mendeleev number M which is a linear function of the Madelung number M [15]: M = M − 1 + δl,0 = n+ l − 1 + δl,0, where δj,k is the Kronecker delta. It is easy to show that the number of elements in the Mendeleev M-th period of the periodic system according to the [n+ l, n]-rule is equal to LM = K(M + 2)−K(M + 1) = 2 ([ M 2 ] + 1 )2 , where [x] is an integer part of the real number x. Numbers LM, M = 1, 2, . . . form the sequence 2, 8, 8, 18, 18, 32, 32, . . . , which coincides with empirical lengths (1.1) of the periods of the system of elements. According to the [n+ l, n]-rule, the number ZM for the initial element of the Mendeleev M-th period is equal to ZM = K(M + 1)− 1, and the number MZ for the final element of the Mendeleev M-th period is equal to MZ = K(M + 2)− 2. The sequence ZM = 1, 3, 11, 19, 37, 55, 87, . . . corresponds to alkaline metals, and the sequence MZ = 2, 10, 18, 36, 54, 86, . . . corresponds to noble gases. Thus, the empirical [n + l, n]-rule is a very efficient method to explain periodic table and atomic properties. Research to justify it continues up to now (see, e.g., [2, 13, 14, 26, 27]. But after 80 years of studies we have not yet good understanding for it. Standard textbooks on quantum mechanics even do not mention this rule. In the present paper we give a theoretical basis for the [n+ l, n]-rule. A preliminary version of this paper was published in 2002 [6]. Further we use everywhere the atomic units ~ = e = m = 1, where ~ is the reduced Plank constant, −e is the electron charge, and m is the electron mass. 3See https://iupac.org/what-we-do/periodic-table-of-elements/. https://iupac.org/what-we-do/periodic-table-of-elements/ 4 E.D. Belokolos 3 Atomic potential in the mean-field and semi-classical approximations The Hamiltonian of a free neutral atom is HA = Z∑ k=1 ( −1 2 ∆k − Z ~rk ) + 1 2 Z∑ i,k=1 1 ~rik = Z∑ k=1 ( −1 2 ∆k + v(~rk) ) , v(~rk) = −Z ~rk + 1 2 Z∑ i=1 1 ~rik , where v(~r) describes the electron interactions with the atomic nucleus and other electrons. Since it is good approximation to describe the electrons in an atom in terms of the electron configuration, i.e., the electron distribution on one-particle states, n, l, m, σ, therefore we may consider the electron interaction v(~r) in the mean-field approximation with the Hamiltonian H = Z∑ k=1 ( −1 2 ∆k + V (~rk) ) , where V (~r) is the mean-field atomic potential. Since we consider a one-particle angular momen- tum, quantum number l as a good quantum number, this potential has to be central, V (~r) = V (r). For the Schrödinger equation HΨ = EΨ we look for the ground state solution in the Slater determinant form Ψ = det ||ψj(rk)||, where one-particle wave function ψk(r) satisfies the one-particle Schrödinger equation describing an electron in a 3-dimensional central potential V (r),( −1 2 ∆ + V (r) ) ψj(r) = Ejψj(r). In the semiclassical approximation, solutions of this equation have the form ψ(r) = exp(iS(r)), where S(r) is the classical action. 3.1 The Thomas–Fermi atomic potential Since in the case under consideration we have three integrals of motion (energy E, angular momentum L, and its projection Lz) our problem is integrable by quadratures due to the Arnold–Liouville theorem, and the radial action looks as follows S(r) = 1 π ∫ r [ 2(E − V (r))− I2θ r2 ]1/2 dr, Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential 5 where the radial action variable is Ir = 1 π ∫ rmax rmin [ 2(E − V (r))− I2θ r2 ]1/2 dr. In the semiclassical approximation we should change action variables by quantum numbers, Ir = nr(E, l), Iθ = l. In this way we come to the Bohr–Sommerfeld quantization rule nr(E, l) = 1 π ∫ r+(E,l) r−(E,l) [ 2E − 2V (r)− l2 r2 ]1/2 dr. Let us set in the latter expression E = 0, then nr(0, l) = √ 2 π ∫ r+(0,l) r−(0,l) [ −V (r)− l2 2r2 ]1/2 dr. Thus the number of bound states N in the atom is N = lmax∑ l=0 nr(0, l)2(2l + 1). Approximately a value N looks as follows N ' ∫ lmax l=0 nr(0, l)2(2l + 1)dl = √ 2 π ∫ r+ r− ∫ lmax l=0 [ −V (r)− l2 2r2 ]1/2 2(2l + 1)dldr = √ 2 π 2 3 ∫ r+ r− [−V (r)]3/2 4r2dr = 2 √ 2 3π2 ∫ r+ r− [−V (r)]3/2 4πr2dr = ∫ r+ r− ρ(r)4πr2dr. This is the well-known asymptotic formula [29] for the number of bound states in a central potential with the density of bound states ρ(r) equal to ρ(r) = 2 √ 2 3π2 [−V (r)]3/2. The electrostatic potential of atomic electrons φ(r) = −V (r) satisfies the Poisson equation ∆φ(r) = −4πρ, and therefore ∆φ(r) = 8 √ 2 3π φ3/2(r), φ(r)r = Z, r → 0, φ(r) = 0, r →∞. We can present the Thomas–Fermi potential φ(r) in such a way φ(r) = Z r χ(x), x = r R0 , R0 = bZ−1/3, b = 1 2 ( 3π 4 )2/3 ' 0.885. For large Z the ground state energy ETF(Z) of the Thomas–Fermi atom is the asymptotics for a ground state energy EHF(Z) of the Hartree–Fock atom [20, 21]: lim Z→∞ EHF(Z)/ETF(Z) = 1. 6 E.D. Belokolos The Thomas–Fermi potential does not depend on quantum numbers although naturally it should do. Nevertheless, let us calculate the Zl for the Thomas–Fermi potential [8]. The effective potential ul(r) is ul(r) = −φ(r) + (l + 1/2)2 2r2 = −(l + 1/2)2 2r2 [ζlxχ(x)− 1], where ζl = 2ZR0 (l + 1/2)2 = 2b Z2/3 (l + 1/2)2 . Conditions for the first appearance of the energy level with certain l are ul(r) = 0, u′l(r) = 0 or, which is the same, ul(x) = 0, u′l(x) = 0. It is equivalent to equations ζlxχ(x)− 1 = 0, d dx (xχ(x)) = 0, which have a solution x0 ' 2.104, χ(x0) ' 0.231, x0χ(x0) ' 0.486, χ′(x0) ' −0.110, ζl = [x0χ(x0)] −1 ' 2.056. This means that ZTF l (2l + 1)3 = ( ζl 8b )3/2 = 0.155, and hence ZTF l ' 0.155(2l + 1)3. In the general case ZTFl is not integer and we have to write ZTF l = ∣∣0.155(2l + 1)3 ∣∣, where |x| means the integer which is the closest to the real x. If we change in this formula the coefficient 0.155 to 0.169 ' 1/6 we get a better agreement of the values ZTFl with that in the Mendeleev table (see [12] and [17, Section 73]). In this case the r.h.s. of the formula will coincide with the first summand of the Klechkovski–Hakala expression for Zl (2.1) obtained with the [n+ l]-rule. This example shows that the Thomas–Fermi theory yields only some approximation to the [n+ l, n]-rule. According to the above discussion the electrons in the atom in the mean-field approximation interact by means of central atomic potential, having geometrical O(3) symmetry. The Thomas– Fermi theory takes into account the O(3) symmetry of the atomic Hamiltonian with central potential. We do the next step. The Hamiltonian with any central potential has dynamical O(4) symmetry [4, 10, 25], similar to the Hydrogen atom [5, 9, 28]. Therefore, the atomic Hamiltonian must have an additional integral of motion in involution with respect to integrals of motion corresponding to O(3) symmetry. Further we shall show that such additional integral of motion in involution does exist for the atomic Hamiltonian and leads to [n+ l, n]-rule. Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential 7 3.2 The Tietz atomic potential We begin with certain basic facts on symmetry and integrability in the classical Hamiltonian systems (see [3, Sections 49–51], [16, Section 52] and [33]). Let us assume that on a symplectic 2n-dimensional manifold there are n functions in involu- tion F1, . . . , Fn, {Fi, Fj} = 0, i, j = 1, . . . , n. Consider a level set of functions Fi Mf = {x : Fi(x) = fi, i = 1, . . . , n}. Suppose that n functions Fi are independent on Mf (i.e., the n 1-forms dFi are linearly inde- pendent in every point of Mf ). Then 1. Mf is a smooth manifold, invariant with respect to the phase flow with the Hamiltonian H = F1. 2. If the manifold Mf is compact and connected, then it is diffeomorphic to n-dimensional torus Tn = {(φ1, . . . , φn) (mod 2π)}. 3. Phase flow with the Hamiltonian H defines on Mf a quasi-periodic movement with angular variables (φ1, . . . , φn), dφi dt = ωi(f), φi(t) = φi(0) + ωit, i = 1, . . . , n. Instead of functions F = (F1, . . . , Fn) it is possible to define new functions I = (I1(F ), . . . , In(F )) which are called action variables and which together with angle variables form in the neighborhood of manifold Mf the canonical system of action-angle coordinates. 4. Canonical equations with the Hamiltonian function are integrable in quadratures. 5. If frequencies ω = (ω1, . . . , ωn) are degenerated, i.e., if there exists such an integer-valued vector k = (k1, . . . , kn) ∈ Zn that (k, ω) = 0, then there appears one more single-valued function Fn+1 which is in involution with func- tions F1, . . . , Fn. Now let us go back to our problem: an electron in the central atomic potential. In this case the 3-dimensional electron movement is reduced to 2-dimensional one in the plane perpendicular to the orbital momentum vector. We describe this movement in the semiclassical approximation by the action S = ∫ r √ 2[E − V (r)]− I2θ r2 dr with radial (r, Ir) and orbital (θ, Iθ) action-angle variables, where Ir = 1 π ∫ rmax rmin √ 2[E − V (r)]− I2θ r2 dr. 8 E.D. Belokolos And so in our problem we have a quasi-periodic movement on 2-dimensional torus which is described by the Fourier series G(t) = ∑ l1∈Z ∑ l2∈Z Gl1,l2 exp(l1φr + l2φθ) = ∑ l1∈Z ∑ l2∈Z Gl1,l2 exp[(l1ωr + l2ωθ)t] with 2 basic frequencies ωr = ∂E ∂Ir , ωθ = ∂E ∂Iθ . In the general case these frequencies are independent. But under certain circumstances, as we have, when new additional integral appears, they may become degenerate (or commensurate, or resonance), qωr = pωθ, q, p ∈ Z, g.c.d.(p, q) = 1, where g.c.d.(p, q) means the greatest common divisor of the integers p and q. This condition leads to important consequences. 1. Since the latter relation is equivalent to q ∂E ∂Ir = p ∂E ∂Iθ it means that energy E depends on pIr + qIθ: E = E(pIr + qIθ). And due to the Bohr–Sommerfeld semiclassical quantization rule Ir → nr, Iθ → l, we have E = E(pnr + ql). 2. We have the additional independent integral of motion. 3. The canonical variables of the problem are separated in several systems of coordinates. For example, for the Kepler–Coulomb potential 1/r we have p = q = 1, energy E depends on principal quantum number n = nr + l+ 1, the canonical variables of the problem are separated in polar and parabolic coordinates. Let us study the similar situation for the atomic potential. Consider a system of equations Ir(E) = 1 π ∫ r+(E) r−(E) [ 2E − 2V (r)− I2θ r2 ]1/2 dr, Mp,q = pIr + qIθ, q, p ∈ Z, g.c.d.(p, q) = 1, the first of which is the standard relation between action variables Ir, Iθ, and energy E in central potential and the second one is relation between Ir, Iθ arising as result of frequencies degeneracy. Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential 9 We shall study equations at the energy E = 0 at which new bound states appear from continuous spectrum, when Z grows. Then differentiating equations with respect to Iθ we come to the equation πα = ∫ r+ r− [ − 2V (r)r2 − I2θ ]−1/2 Iθ dr r , where α = −∂Ir ∂Iθ = q p ∈ Q, q, p ∈ Z, g.c.d.(p, q) = 1. From this equation we deduce an expression for atomic potential V (r). A similar integral equation arising in problem for generalized tautochrone (isochrone) curve was solved by Abel [1]. We shall use his method in the form presented in [16, Section 12]. Theorem 1. Equation πα = ∫ r+ r− [ −2V (r)r2 − L2 ]−1/2 L dr r (3.1) has a solution Vα,β,R(r) = − β r2[(r/R)1/α + (R/r)1/α]2 , where β and R are certain constants, and its deformations. If the potential Vα,β,R(r) coincides at small r with the Kepler–Coulomb potential Vα,β,R = −Z/r, r → 0, then we have α = q/p = 2, β = ZR, and the potential takes the following form V (r) ≡ V2,ZR,R(r) = − Z r(1 + (r/R))2 . Remark 1. According to the theorem in the neighborhood E ' 0 we have E = E(n+ l) and this fact proves the first part of the [n + l, n]-rule. A second part of the [n + l, n]-rule is a consequence of the oscillation theorem. Proof. Let us rewrite equation (3.1) in the form πα = ∫ r+ r− [ w(x)− L2 ]−1/2 Ldx, where x = ln(r/R), r = R exp(x), w(x) = −2V (r)r2 at r = R exp(x), w0 = maxw(x), 10 E.D. Belokolos are new variables and R is a parameter. We assume that w(x) is the one-well potential and therefore the inverse function is two-valued, i.e., the values w are reached in two points x−(w) and x+(w). We shall assume also that x−(w) ≤ x+(w) and at w0 we have x−(w0) = x+(w0). As a result we obtain πα = ∫ w0 L2 ( w − L2 )−1/2(dx+ dw − dx− dw ) Ldw. Multiplying this equality by ( L2 − w1 )−1/2 2dL and integrating from (w1) 1/2 to (w0) 1/2 we get 2πα ∫ (w0)1/2 (w1)1/2 ( L2 − w1 )−1/2 dL = ∫ (w0)1/2 (w1)1/2 2LdL ∫ w0 L2 ( dx+ dw − dx− dw )[( w − L2 )( L2 − w1 )]−1/2 dw = ∫ w0 w1 dw ( dx+ dw − dx− dw )∫ (w)1/2 (w1)1/2 [( w2 − L2 )( L2 − w1 )]−1/2 2LdL. Since ∫ (w0)1/2 (w1)1/2 ( L2 − w1 )−1/2 dL = Arcosh(w0/w1) 1/2,∫ (w)1/2 (w1)1/2 [( w2 − L2 )( L2 − w1 )]−1/2 2LdL = π, this equality acquires the following form 2αArcosh(w0/w1) 1/2 = (x+(w)− x−(w)) ∣∣∣w0 w1 . Taking into account that x−(w0) = x+(w0) and setting w1 = w we come to the equality x+(w)− x−(w) = 2αArcosh(w0/w)1/2. This equality defines only the difference x−(w)− x+(w) of two functions x−(w) and x+(w), any of which remains actually undefined. It means that there exists the infinite set of potentials which satisfy the latter equation and differ by deformations which do not change the difference of two values of x corresponding to one value of w. Among these potentials there is a symmetric potential with the property x+(w) = −x−(w) ≡ x(w). For this potential in this case w(x) = w0/ cosh2(x/α), or, in the previous notations, Vα,β,R(r) = − β r2[(r/R)1/α + (R/r)1/α]2 , where β = 2w0. These potentials are degenerate for arbitrary value of the parameter α. In order for the degenerate potential to coincide with Coulomb potential at small r Vα,β,R = −Z/r, r → 0, we must have α = q/p = 2, β = ZR. Then this potential takes the following form V (r) ≡ V2,ZR,R(r) = − Z r(1 + (r/R))2 . � Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential 11 In this potential the total number of bound states N is equal to the total number of elect- rons Z if N = 8 √ 2 3π ∫ [−V (r)]3/2dr = ( 9/2R3 ) = Z, i.e., R = (9/2Z)1/3 ' 1.651Z−1/3, and we have finally V (r|Z) = − Z r(1 + (r/R(Z)))2 . For the given electron configuration we must use the Klechkovski–Hakala formulas for Z. In this case the potential will depend on quantum numbers. We can present this potential as the sum V (x) = − (Z/R) x(1 + x)2 = −Z R [ 1 x − 1 (x+ 1) − 1 (x+ 1)2 ] , x = r R , where the first summand describes the Coulomb attraction of the atomic nucleus for the single electron and two other summands describe a nucleus screening by the other electrons. T. Tietz [30] proposed the potential V (r) = − Z r[1 + (r/R)]2 as a good rational approximation to the Thomas–Fermi potential and used it for calculation of various atomic properties and explanation of the periodic system elements [31] (see also [34]). Further we shall call this potential as the Tietz atomic potential. Due to the proximity of the Tietz and Thomas–Fermi potentials it is very likely that at large Z a ground state energy of the Tietz atom, ET(Z), is an asymptotics for the ground state energy of the Hartree–Fock atom, EHF(Z): lim Z→∞ EHF(Z)/ET(Z) = 1. Yu.N. Demkov and V.N. Ostrovski [7] pointed out that this potential is a particular case of a so called focussing (in other words “degenerate”) potentials studied in connection to certain problems of optics by J.C. Maxwell [23] and V. Lenz [19]. The Tietz potential is a rational function, and thus we can easily do various calculations with it. For example, we can calculate the atomic spectrum. 4 Semiclassical atomic spectrum for the Tietz atomic potential The Bohr–Sommerfeld semiclassical condition of quantization for the Tietz atomic potential is nr = 1 π ∫ r+ r− [ 2E + 2Z r(1 + (r/R))2 − (l + 1/2)2 r2 ]1/2 dr. In the scaled quantum numbers x = r R , ε = 2ER2 (l + 1/2)2 , νr = nr l + 1/2 , ηl = 2ZR (l + 1/2)2 12 E.D. Belokolos the Bohr–Sommerfeld equation takes the form νr = 1 π ∫ x+ x− [ ε+ ηl x(1 + x)2 − 1 x2 ]1/2 dx = 1 π ∫ x+ x− √ P4(x) x(1 + x) dx, where P4(x) is a 4-th degree polynomial P4(x) = εx2(1 + x)2 + ηlx− (1 + x)2, and boundaries of integration are real non-negative zeros of this polynomial. Therefore scaled radial quantum number νr is a period of the elliptic integral with scaled energy ε and scaled charge ηl as the parameters. We can obtain an atomic spectrum in the semiclassical approxi- mation ε = f(νr, ηl), by means of inversion of the elliptic integral. Let us study this problem in the framework of perturbation theory in the energy ε. We present the Bohr–Sommerfeld semiclassical equation in the form νr = 1 π ∫ x+ x− [ ε+ ηlx− (1 + x)2 x2(1 + x)2 ]1/2 dx = 1 π ∫ x+ x− [ ε+ (x+ − x)(x− x−) x2(1 + x)2 ]1/2 dx, where ηlx− (1 + x)2 = (x+ − x)(x− x−), x2± − (ηl − 2)x± + 1 = 0, x± = (ηl − 2)± √ (ηl − 2)2 − 4 2 = (ηl − 2)± √ ηl(ηl − 4) 2 , −(ηl − 4) > ε > −1 4 ηl(ηl − 4). At the first approximation we have νr ' J0 + εJ1 = 1 π ∫ x+ x− [(x+ − x)(x− x−)]1/2 dx x(1 + x) + ε 1 2π ∫ x+ x− x(1 + x) [(x+ − x)(x− x−)]1/2 dx, where J0 = 1 π ∫ x+ x− [(x+ − x)(x− x−)]1/2 dx x(1 + x) = √ ηl − 2, J1 = 1 2π ∫ x+ x− x(1 + x) [(x+ − x)(x− x−)]1/2 dx = 1 16 ( 3η2l − 8ηl ) . The energy spectrum in the semiclassical theory is εn,l = −J0 − νr J1 = −16 √ ηl − 2− νr( 3η2l − 8ηl ) = −16 (l + (1/2)) √ ηl − 2l − 1− nr (l + (1/2)) ( 3η2l − 8ηl ) = −16 (l + (1/2)) √ ηl − (n+ l) (l + (1/2)) ( 3η2l − 8ηl ) = −16 √ 2ZR−M (l + (1/2)) ( 3η2l − 8ηl ) , where we have used the Madelung number M = n+ l. Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential 13 Since εM,l = 2EM,lR 2 (l + (1/2))2 , we have the following expression for the atomic spectrum in the atomic units EM,l = −8 ( √ 2ZR−M)(l + (1/2))( 3η2l − 8ηl ) R2 , ηl = 2ZR (l + (1/2))2 . If √ 2ZR = M then EM,l = 0, and all states (M, l) are degenerate with respect to l. Thus at energy E = 0 there appears the full set of energy levels (n, l) with M = n + l. Since R = (9/2Z)1/3 the equality M = √ 2ZR is equivalent to Z = (1/6)M3. We can write down a complete perturbation series νr = ∞∑ k=0 εkJk. (4.1) Since all integrals Jk = ck π ∫ x+ x− { x2(1 + x)2 [(x+ − x)(x− x−)] }k−(1/2) dx, ck = Γ(3/2) Γ(k + 1)Γ((3/2)− k) , k ≥ 2, are divergent we should take their regularized values. i.e., valeur principale. Inverting series (4.1) by means of the Bürmann–Lagrange theorem we can obtain an expression for atomic spectrum in the semiclassical approximation ε = f(νr, ηl). We plans to compare this atomic spectrum (both eigenvalues and eigenfunctions) with those presented in articles [18, 32]. 5 Conclusion For the Mendeleev periodic system of elements we have proved the empirical (n+l, n)-rule which explains very efficiently the structure and properties of chemical elements. In order to prove this rule we are forced to build for the Hamiltonian describing an electron in central atomic potential one more integral of motion in involution in addition to energy and orbital integral. In this case the atomic potential appears to be the Tietz potential. For the Tietz potential we have calculated the atomic energy spectrum in the semiclassical approximation. In addition, we are going to study spectrum of the Schrödinger operator with the Tietz potential, when the radial part of wave function satisfies the confluent Heun equation and compare energy levels, obtained in such a way, with the NIST Atomic Spectra Database. 14 E.D. Belokolos For the atom in the paper we have studied the non-relativistic atomic Hamiltonian. We plan to consider also the relativistic Dirac Hamiltonian. It is interesting also to study isospectral many-well deformation of the Tietz potential. Questions studied in the paper are related only to the ground state of atoms. But according to experimental data the weakly excited atomic states also follow to [n+ l, n]-rule [15]. It may be of interest for physical chemistry, especially for understanding chemical reactions rules. We hope to do that in the near future. References [1] Abel N.H., Solution de quelques problémes à l’aide d’intégrales définies, in Œuvres completes de Niels Henrik Abel, Editors L. Sylow, S. Lie, Christiania, Norway, 1881, 11–27. [2] Allen L.C., Knight E.T., The Löwdin challenge: origin of the n + l, n (Madelung) rule for filling the orbital configurations of the periodic table, Int. J. Quantum Chem. 90 (2002), 80–88. [3] Arnol’d V.I., Mathematical methods of classical mechanics, 3rd ed., Nauka, Moscow, 1989. [4] Bacry H., Ruegg H., Souriau J.M., Dynamical groups and spherical potentials in classical mechanics, Comm. Math. Phys. 3 (1966), 323–333. [5] Bargmann V., Zur Theorie des Wasserstoffatoms. Bemerkungen zur gleichnamigen Arbeit von V. Fock, Z. Phys. 99 (1936), 576–582. [6] Belokolos E.D., The integrability and the structure of atom, J. Math. Phys. Anal. Geometry 9 (2002), 339–351. [7] Demkov Y.N., Ostrovskii V.N., Internal symmetry of the Maxwell “fish eye” problem and the Fock group for the hydrogen atom, Sov. Phys. JETP 33 (1971), 1083–1087. [8] Fermi E., Über die Anwendung der statistischen Methode auf die Probleme des Atombaues, in Quantenthe- orie und Chemie, Editor H. Falkenhagen, Leipziger Vorträge, S. Hirzel-Verlag, Leipzig, 1928, 95–111. [9] Fock V., Zur Theorie des Wasserstoffatoms, Z. Phys. 98 (1935), 145–154. [10] Fradkin D.M., Existence of the dynamic symmeries O(4) and SU3 for all classical central potential problems, Progr. Theoret. Phys. 37 (1967), 798–812. [11] Hakala R., The periodic law in mathematical form, J. Phys. Chem. 56 (1952), 178–181. [12] Ivanenko D., Larin S., Theory of periodic system of elements, Dokl. Akad. Nauk USSR 88 (1953), 45–48. [13] Kibler M.R., On a group-theoretical approach to the periodic table of chemical elements, quant-ph/0408104. [14] Kibler M.R., From the Mendeleev periodic table to particle physics and back to the periodic table, quant-ph/0611287. [15] Klechkovski V.M., The distribution of atomic electrons and the filling rule for (n + l)-groups, Atomizdat, Moscow, 1968. [16] Landau L.D., Lifshits E.M., Theoretical physics, Vol. 1, Mechanics, Nauka, Moscow, 1973. [17] Landau L.D., Lifshits E.M., Theoretical physics, Vol. 3, Quantum mechanics: nonrelativistic theory, Nauka, Moscow, 1974. [18] Latter R., Atomic energy levels for the Thomas–Fermi and Thomas–Fermi–Dirac potentials, Phys. Rev. 99 (1955), 510–519. [19] Lenz W., Zur Theorie der optischen Abbildung, in Probleme der mathematischen Physik, Leipzig, 1928, 198–207. [20] Lieb E.H., Thomas–Fermi and related theories of atoms and molecules, Rev. Modern Phys. 53 (1981), 603–641. [21] Lieb E.H., Simon B., The Thomas–Fermi theory of atoms, molecules and solids, Adv. Math. 23 (1977), 22–116. [22] Madelung E., Die mathematischen Hilfsmittel des Physikers, Die Grundlehren der Mathematischen Wis- senschaften, Vol. 4, Springer-Verlag, Berlin, 1936. [23] Maxwell J.C., Solutions of problems, in The Scientific Papers of James Clerk Maxwell, Editor W.D. Niven, Dover Publications, New York, 1952, 76–79. https://doi.org/10.1007/BF01645086 https://doi.org/10.1007/BF01645086 https://doi.org/10.1007/BF01645086 https://doi.org/10.1007/BF01338811 https://doi.org/10.1007/BF01336904 https://doi.org/10.1143/PTP.37.798 https://doi.org/10.1021/j150494a006 http://arxiv.org/abs/quant-ph/0408104 http://arxiv.org/abs/quant-ph/0611287 https://doi.org/10.1103/PhysRev.99.510 https://doi.org/10.1103/RevModPhys.53.603 https://doi.org/10.1016/0001-8708(77)90108-6 https://doi.org/10.1007/978-3-662-02177-4 https://doi.org/10.1007/978-3-662-02177-4 Mendeleev Table: a Proof of Madelung Rule and Atomic Tietz Potential 15 [24] Meek T.I., Allen L.C., Configuration irregularities: deviations from the Madelung rule and inversion of orbital energy levels, Chem. Phys. Lett. 362 (2002), 362–364. [25] Mukunda N., Dynamical symmeries and classical mechanics, Phys. Rev. 155 (1967), 1383–1386. [26] Ostrovsky V.N., Dynamic symmetry of atomic potential, J. Phys. B: At. Mol. Phys. 14 (1981), 4425–4439. [27] Ostrovsky V.N., What and how physics contributes to understanding the periodic law?, Found. Chem. 3 (2001), 145–182. [28] Pauli W., Über das Wasserstoffspektrum vom Standpunkt der neuen Quantenmechanik, Z. Phys. 36 (1926), 336–363. [29] Reed M., Simon B., Methods of modern mathematical physics. IV. Analysis of operators, Academic Press, New York – London, 1978. [30] Tietz T., Über eine Approximation der Fermischen Verteilungsfunktion, Ann. Physik 450 (1955), 186–188. [31] Tietz T., Elektronengruppen im periodischen System der Elemente in der statistischen Theorie des Atoms, Ann. Physik 460 (1960), 237–240. [32] Tietz T., Über Eigenwerte und Eigenfunktionen der Schrödinger-Gleichung für das Thomas–Fermische Po- tential, Acta Phys. Acad. Sci. Hungar. 11 (1960), 391–400. [33] Torrielli A., Classical integrability, J. Phys. A: Math. Theor. 49 (2016), 323001, 31 pages, arXiv:1606.02946. [34] Wong D.P., Theoretical justification of Madelung’s rule, J. Chem. Educ. 56 (1979), 714–717. https://doi.org/10.1016/S0009-2614(02)00919-3 https://doi.org/10.1103/PhysRev.155.1383 https://doi.org/10.1088/0022-3700/14/23/008 https://doi.org/10.1023/A:1011476405933 https://doi.org/10.1007/BF01450175 https://doi.org/10.1002/andp.19554500309 https://doi.org/10.1002/andp.19604600313 https://doi.org/10.1007/BF03160433 https://doi.org/10.1088/1751-8113/49/32/323001 http://arxiv.org/abs/1606.02946 https://doi.org/10.1021/ed056p714 1 Introduction 2 The Madelung [n+l,n]-rule 3 Atomic potential in the mean-field and semi-classical approximations 3.1 The Thomas–Fermi atomic potential 3.2 The Tietz atomic potential 4 Semiclassical atomic spectrum for the Tietz atomic potential 5 Conclusion References