Next Article in Journal
Electron-Impact Excitation of the λ190.8 nm and λ179.9 nm Intercombination Lines in the Tl+ Ion
Next Article in Special Issue
Electron Correlations in Sequential Two-Photon Double Ionization of an Ar Atom
Previous Article in Journal
Domainal Investigation of a Quartz-Fluorite Composite Using Spectroscopic Techniques
Previous Article in Special Issue
Spectroscopic Peculiarities at Ionization of Excited 2p5(2PJf)3s[K]0,1,2 States of Ne: Cooper Minima and Autoionizing Resonances
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fock Expansion for Two-Electron Atoms: High-Order Angular Coefficients

by
Evgeny Z. Liverts
1,* and
Rajmund Krivec
2
1
Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel
2
Department of Theoretical Physics, J. Stefan Institute, 1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Atoms 2022, 10(4), 135; https://doi.org/10.3390/atoms10040135
Submission received: 23 September 2022 / Revised: 1 November 2022 / Accepted: 2 November 2022 / Published: 7 November 2022

Abstract

:
The Fock expansion, which describes the properties of two-electron atoms near the nucleus, is studied. The angular Fock coefficients ψ k , p ( α , θ ) with the maximum possible value of subscript p are calculated on examples of the coefficients with 5 k 10 . The presented technique makes it possible to calculate such angular coefficients for any arbitrarily large k. The mentioned coefficients being leading in the logarithmic power series representing the Fock expansion, they may be indispensable for the development of simple methods for calculating the helium-like electronic structure. The theoretical results obtained are verified by other suitable methods. The Wolfram Mathematica is used extensively.

1. Introduction

The properties of a two-electron atomic (helium-like) system with an infinitely massive nucleus of charge Z and nonrelativistic energy E are defined by the wave function (WF) Ψ ( r 1 , r 2 , r 12 ) , where r 1 and r 2 are the electron–nucleus distances, and r 12 is the distance between the electrons. The behavior of the ground state WF in the vicinity of the nucleus located at the origin is determined by the Fock expansion [1]
Ψ ¯ ( r 1 , r 2 , r 12 ) Ψ ( R , α , θ ) = k = 0 R k p = 0 [ k / 2 ] ψ k , p ( α , θ ) ln p R ,
where the hyperspherical coordinates R , α and θ are defined by the relations:
R = r 1 2 + r 2 2 , α = 2 arctan r 2 r 1 , θ = arccos r 1 2 + r 2 2 r 12 2 2 r 1 r 2 .
The convergence of expansion (1) was proven in Ref. [2]. The angular Fock coefficients (AFC) ψ k , p satisfy the Fock recurrence relation (FRR)
Λ 2 k ( k + 4 ) ψ k , p ( α , θ ) = h k , p ( α , θ )
with the RHS of the form [3,4]:
h k , p = 2 ( k + 2 ) ( p + 1 ) ψ k , p + 1 + ( p + 1 ) ( p + 2 ) ψ k , p + 2 2 V ψ k 1 , p + 2 E ψ k 2 , p .
The dimensionless Coulomb potential representing the electron–electron and electron–nucleus interactions is
V R r 12 Z R r 1 + R r 2 = 1 ξ 2 Z η sin α ,
where we have introduced the important (in what follows) angular quantities:
ξ = 1 sin α cos θ , η = 1 + sin α .
The hyperspherical angular momentum operator, projected on S states, is defined as
Λ 2 = 4 2 α 2 + 2 cot α α + 1 sin 2 α 2 θ 2 + cot θ θ .
It is clear that all circumnuclear features of the two-electron atoms (ions) are defined by the Fock expansion (1). There are a large number of methods for calculating the electronic structure of the two-electron atomic systems. An excellent review on this topic can be found in Refs. [3,5,6,7,8]. However, we know only one technique that correctly represents the WF Ψ ( r 1 , r 2 , r 12 ) near the nucleus. It is the so-called correlation function hyperspherical harmonic method (CFHHM) [9,10,11]. The expansion in hyperspherical harmonics (HHs) provides the correct representation of the AFCs. However, the HH expansion is known to converge very slowly. Although this method makes it possible to increase the convergence of the HH expansion, a sufficiently good accuracy requires a large HHs’ basis size, which, in turn, creates great computational difficulties.
Thus, it would be extremely useful to develop a much simpler method for calculating the WF with correct behavior near the nucleus. In this regard, we would like to emphasize the following important peculiarities of the Fock expansion (FE). It follows from definition (1) that the FE can be split into individual power series (lines) associated with definite power of ln R . In other words, the FE can be represented in the form:
Ψ = ( ln R ) 0 ψ 0 , 0 + R ψ 1 , 0 + R 2 ψ 2 , 0 + . . . + ( ln R ) 1 R 2 ψ 2 , 1 + R ψ 3 , 1 + R 2 ψ 4 , 1 + . . . + ( ln R ) 2 R 4 ψ 4 , 2 + R ψ 5 , 2 + R 2 ψ 6 , 2 + . . . + ( ln R ) 3 R 6 ψ 6 , 3 + R ψ 7 , 3 + R 2 ψ 8 , 3 + . . . + ( ln R ) 4 R 8 ψ 8 , 4 + R ψ 9 , 4 + R 2 ψ 10 , 4 + . . . + . . .
It is seen that the leading term of each line represents the product ( ln R ) k / 2 R k ψ k , k / 2 ( α , θ ) with even k. The first AFCs ( ψ 0 , 0 = 1 ) corresponding to k = 0 , 2 , 4 are well-known (see, e.g., [3,4]):
ψ 1 , 0 = 1 2 ξ Z η ,
ψ 2 , 1 = Z ( π 2 ) 3 π ( 1 ξ 2 ) ,
ψ 3 , 1 = Z ( π 2 ) 36 π 6 Z η ( 1 ξ 2 ) + ξ ( 5 ξ 2 6 ) ,
ψ 4 , 2 = Z 2 ( π 2 ) ( 5 π 14 ) 540 π Y 40 ( α , θ ) + 2 Y 42 ( α , θ ) .
The normalized HHs are
Y 40 ( α , θ ) = π 3 / 2 ( 4 cos 2 α 1 ) , Y 42 ( α , θ ) = 2 2 π 3 / 2 sin 2 α P 2 ( cos θ ) ,
where the P n ( x ) denote the Legendre polynomials.
In this paper, we present the theoretical calculations of the AFCs ψ 5 , 2 ( α , θ ) , ψ 6 , 3 ( α , θ ) , ψ 7 , 3 ( α , θ ) , ψ 8 , 4 ( α , θ ) and ψ 9 , 4 ( α , θ ) included into the k = 4 , k = 6 and k = 8 “lines” of the expansion (8), and also the AFC ψ 10 , 5 ( α , θ ) representing the leading term of the k = 10 “line”. It is important to note that all mentioned angular coefficients represent the AFCs ψ k , p with the maximum possible p for a given k.

2. Derivation of the Angular Fock Coefficient ψ 5 , 2 ( α , θ )

The FRR (3) and (4) for k = 5 and p = 2 reduce to the form
Λ 2 45 ψ 5 , 2 ( α , θ ) = h 5 , 2 ( α , θ ) ,
where
h 5 , 2 ( α , θ ) = 2 V ψ 4 , 2 ( α , θ ) .
Using Equations (5), (6), (12) and (13), it is convenient to represent the RHS of Equation (14) in the form
h 5 , 2 ( α , θ ) = Z 2 ( π 2 ) ( 5 π 14 ) 270 π 3 π 3 / 2 ( 2 h 1 + h 2 ) 2 Z ( h 3 + 2 h 4 ) ,
where
h 1 = ( 1 ξ 2 ) 2 ξ , h 2 = cos ( 2 α ) ξ , h 3 = η Y 40 ( α , θ ) sin α , h 4 = η Y 42 ( α , θ ) sin α .
Accordingly, we obtain the solution of Equation (14) in the identical form
ψ 5 , 2 ( α , θ ) = Z 2 ( π 2 ) ( 5 π 14 ) 270 π 3 π 3 / 2 ( 2 f 1 + f 2 ) 2 Z ( f 3 + 2 f 4 ) ,
where the AFC-components f i satisfy the individual Fock recurrence relations (IFRRs)
Λ 2 45 f i = h i . ( i = 1 , 2 , 3 , 4 )
We sequentially find solutions to each of the IFRRs (19) using various methods presented in Ref. [4].

2.1. Solution of the IFRR ( Λ 2 45 ) f 3 = η Y 40 / sin α

Moving from simpler to more complex solutions, let us start with IFRR
Λ 2 45 f 3 = h 3 .
The RHS h 3 h 3 ( α ) represents the function of only one angle variable α . It was shown [4] that the solution of the corresponding IFRR (20) reduces to the solution g ( ρ ) = f 3 ( α ) of the inhomogeneous differential equation
( ρ 2 + 1 ) 2 g ( ρ ) + 2 ρ 1 ( ρ 2 + 1 ) g ( ρ ) + 45 g ( ρ ) = h ( ρ ) ,
where ρ = tan ( α / 2 ) , and
h ( ρ ) h 3 ( α ) = ( 4 cos 2 α 1 ) 1 + sin α π 3 / 2 sin α = ( ρ + 1 ) ( 3 ρ 4 10 ρ 2 + 3 ) 2 π 3 / 2 ρ ( ρ 2 + 1 ) 3 / 2 .
For convenience, we solve the Equation (21) with the RHS h ( ρ ) not containing the multiplier π 3 / 2 . The final solution f 3 will be multiplied by this factor.
Using the method of variation of parameters, one obtains [4] the particular solution of Equation (21) in the form
g ( p ) ( ρ ) = v 50 ( ρ ) u 50 ( ρ ) h ( ρ ) d ρ ( ρ 2 + 1 ) 2 W 0 ( ρ ) u 50 ( ρ ) v 50 ( ρ ) h ( ρ ) d ρ ( ρ 2 + 1 ) 2 W 0 ( ρ ) ,
where
W 0 ( ρ ) = ( ρ 2 + 1 ) / ρ 2 .
The independent solutions of the homogeneous equation associated with Equation (21) are [4]:
u 50 ( ρ ) = ( ρ 2 + 1 ) 9 / 2 ρ 2 F 1 4 , 7 2 ; 1 2 ; ρ 2 = 1 7 ρ 2 ( 3 5 ρ 2 + ρ 4 ) ρ ( ρ 2 + 1 ) 5 / 2 ,
v 50 ( ρ ) = ( ρ 2 + 1 ) 9 / 2 2 F 1 4 , 9 2 ; 3 2 ; ρ 2 = 1 35 ρ 2 + 21 ρ 4 ρ 6 7 ( ρ 2 + 1 ) 5 / 2 ,
where 2 F 1 ( . . . ) is the Gaussian hypergeometric function. The substitution of Equations (24)–(26) into the general representation (23) yields
g ( p ) ( ρ ) = 7 ρ ρ 5 ρ ( 3 ρ 4 ) ( 3 ρ + 5 ) 24 + 23 23 420 ρ ( ρ 2 + 1 ) 5 / 2 .
The general solution of the inhomogeneous equation can be expressed as the sum of the general solution of the associated homogeneous (complementary) equation and the particular solution of the inhomogeneous equation, whence
g ( ρ ) = g ( p ) ( ρ ) + c u u 50 ( ρ ) + c v v 50 ( ρ ) ,
where the coefficients c u and c v are currently undetermined. To choose these coefficients, it is necessary to determine the behavior of all independent solutions on the boundaries of the domain [ 0 , ] . We easily obtain:
g ( p ) ( ρ ) = ρ 0 23 420 ρ 23 60 + 451 ρ 840 + O ( ρ 2 ) ,
g ( p ) ( ρ ) = ρ 3 4 ρ 1 4 ρ 2 + 85 24 ρ 3 + O ( ρ 4 ) ,
u 50 ( ρ ) = ρ 0 1 ρ 47 ρ 2 + O ( ρ 3 ) ,
u 50 ( ρ ) = ρ 7 + 105 2 ρ 2 + O ( ρ 4 ) ,
v 50 ( ρ ) = ρ 0 1 15 ρ 2 2 + O ( ρ 4 ) ,
v 50 ( ρ ) = ρ ρ 7 + 47 14 ρ + O ( ρ 3 ) .
It is seen that the particular solution g ( p ) ( ρ ) is divergent at ρ = 0 , whereas the solutions of the homogeneous equation associated with Equation (21) are divergent, at ρ = 0 and ρ = for u 50 ( ρ ) and v 50 ( ρ ) , respectively. Thus, to avoid the divergence on the whole range of definition, one should set c u = 23 / 420 and c v = 0 in the general solution (28). Then, the final physical solution becomes
f 3 ( α ) = ( ρ + 1 ) ( 23 ρ 4 + 22 ρ 3 122 ρ 2 + 22 ρ + 23 ) 60 π 3 / 2 ( ρ 2 + 1 ) 5 / 2 = = 1 60 π 3 / 2 11 sin α + 21 cos ( 2 α ) + 2 1 + sin α .

2.2. Solution of the IFRR ( Λ 2 45 ) f 4 = η Y 42 / sin α

It was shown in Ref. [4] that the solution f 4 f 4 ( α , θ ) of the IFRR
Λ 2 45 f 4 = 2 π 3 / 2 η 2 sin α P 2 ( cos θ )
can be found in the form
f 4 = 2 2 π 3 / 2 sin 2 α P 2 ( cos θ ) g 4 ( ρ ) ,
where the function g 4 ( ρ ) satisfies the equation
( ρ 2 + 1 ) 2 g 4 ( ρ ) + 2 ρ 1 [ 1 + ρ 2 + 2 ( 1 ρ 4 ) ] g 4 ( ρ ) + 13 g 4 ( ρ ) = h 4 ( ρ ) ,
with
h 4 ( ρ ) = η sin α = ( ρ + 1 ) ρ 2 + 1 2 ρ .
Using the method of variation of parameters, one obtains [4] the particular solution of Equation (38) in the form
g 4 ( p ) ( ρ ) = v 52 ( ρ ) u 52 ( ρ ) h 4 ( ρ ) d ρ ( ρ 2 + 1 ) 2 W 2 ( ρ ) u 52 ( ρ ) v 52 ( ρ ) h 4 ( ρ ) d ρ ( ρ 2 + 1 ) 2 W 2 ( ρ ) ,
where
W 2 ( ρ ) = 5 ρ ρ 2 + 1 ρ 5 .
The independent solutions of the homogeneous equation associated with Equation (38) are:
u 52 ( ρ ) = ( ρ 2 + 1 ) 13 / 2 ρ 5 2 F 1 4 , 3 2 ; 3 2 ; ρ 2 = 1 + 11 ρ 2 + 99 ρ 4 231 ρ 6 ρ 5 ρ 2 + 1 ,
v 52 ( ρ ) = ( ρ 2 + 1 ) 13 / 2 2 F 1 4 , 13 2 ; 7 2 ; ρ 2 = 231 99 ρ 2 11 ρ 4 ρ 6 231 ρ 2 + 1 .
Thus, the particular solution (40) reduces to the form:
g 4 ( p ) ( ρ ) = 11 ( 21 ρ 5 + 9 ρ 4 + ρ 2 ) + 1 2772 ρ 5 ρ 2 + 1 .
Considering the series expansions for g 4 ( p ) ( ρ ) , u 52 ( ρ ) and v 52 ( ρ ) on the boundaries of the range of definition ( ρ [ 0 , ] ), it can be shown that function
g 4 ( p ) ( ρ ) + 1 2772 u 52 ( ρ ) = ρ + 1 12 ρ 2 + 1 = 1 12 1 + sin α
represents the physical (finite) solution of Equation (38), and hence we finally obtain:
f 4 ( α , θ ) = 2 6 π 3 / 2 ( sin α ) 2 1 + sin α P 2 ( cos θ ) .

2.3. Solution of the IFRR ( Λ 2 45 ) f 1 = ( 1 ξ 2 ) 2 / ξ

It is important to note that the RHS h 1 (see Equation (17)) of the IFRR
Λ 2 45 f 1 = h 1 ( ξ )
is a function of ξ (only) defined by Equation (6). For this case [4,12], the solution of Equation (47) coincides with the solution of the inhomogeneous differential equation
( ξ 2 2 ) f 1 ( ξ ) + ξ 1 ( 5 ξ 2 4 ) f 1 ( ξ ) 45 f 1 ( ξ ) = h 1 ( ξ ) .
A particular solution of Equation (48) can be found by the method of variation of parameters in the form [12]
f 1 ( p ) ( ξ ) = 1 7 2 u 5 ( ξ ) v 5 ( ξ ) w ( ξ ) d ξ v 5 ( ξ ) u 5 ( ξ ) w ( ξ ) d ξ ,
where
w ( ξ ) = h 1 ( ξ ) ξ 2 2 ξ 2 .
The linearly independent solutions of the homogeneous equation associated with Equation (48) are defined by the relations
u 5 ( ξ ) = P 13 / 2 1 / 2 ( ξ / 2 ) ξ 2 ξ 2 4 = 2 1 / 4 ( 8 ξ 6 28 ξ 4 + 28 ξ 2 7 ) π ( 2 ξ 2 ) ,
v 5 ( ξ ) = Q 13 / 2 1 / 2 ( ξ / 2 ) ξ 2 ξ 2 4 = π ( 8 ξ 6 + 20 ξ 4 12 ξ 2 + 1 ) 2 3 / 4 ξ .
where P ν μ ( x ) and Q ν μ ( x ) are the associated Legendre functions of the first and second kind, respectively. A substitution of the representations (50)–(52) into (49) yields the particular solution:
f 1 ( p ) ( ξ ) = 1 60 ξ ( 13 ξ 4 30 ξ 2 + 15 ) .
It can be verified that the particular solution f 1 ( p ) ( ξ ) is finite on the whole range of definition ( ξ [ 0 , 2 ] ), whereas the solutions of the homogeneous equation associated with Equation (48) are divergent at ξ = 2 ( α = π / 2 , θ = π ) and ξ = 0 ( α = π / 2 , θ = 0 ) for u 5 ( ξ ) and v 5 ( ξ ) , respectively. Thus, we can conclude that the final physical solution of the IFFR (47) coincides with the particular solution (53), whence
f 1 = 1 60 ξ ( 13 ξ 4 30 ξ 2 + 15 ) = 1 60 1 sin α cos θ sin α cos θ ( 4 + 13 sin α cos θ ) 2 .

2.4. Solution of the IFRR ( Λ 2 45 ) f 2 = cos ( 2 α ) / ξ

To solve the IFRR
Λ 2 45 f 2 = h 2 ,
with the RHS h 2 defined by Equation (17), first of all, it is necessary to recall Sack’s representation [13] (see also [3,4]) for ξ ν with ν = 1 :
ξ 1 = l = 0 P l ( cos θ ) sin α 2 l F l ( ρ ) ,
where
F l ( ρ ) = 2 F 1 l 2 + 1 4 , l 2 + 3 4 ; l + 3 2 ; 4 ρ 2 ( ρ 2 + 1 ) 2 = F l ( ρ ) 0 ρ 1 F l ( 1 / ρ ) ρ 1
with
F l ( ρ ) = ( ρ 2 + 1 ) l + 1 2 .
This enables us to present the RHS of Equation (55) in the form
h 2 cos ( 2 α ) ξ = l = 0 P l ( cos θ ) ( sin α ) l h l ( ρ ) ,
where
h l ( ρ ) = 2 l F l ( ρ ) cos ( 2 α ) = 2 l F l ( ρ ) 1 8 ρ 2 ( ρ 2 + 1 ) 2 .
In turn, it was shown in Ref. [4] that in the case where the RHS is determined by Equation (59), the solution of the corresponding IFRR (55) can be found in the form
f 2 ( α , θ ) = l = 0 P l ( cos θ ) ( sin α ) l σ l ( ρ ) ,
where the function σ l ( ρ ) satisfies the inhomogeneous differential equation
( ρ 2 + 1 ) 2 σ l ( ρ ) + 2 ρ 1 [ 1 + ρ 2 + l ( 1 ρ 4 ) ] σ l ( ρ ) + ( 5 2 l ) ( 9 + 2 l ) σ l ( ρ ) = h l ( ρ ) .
The linearly independent solutions of the homogeneous equation associated with Equation (62) are:
u 5 l ( ρ ) = ( ρ 2 + 1 ) l + 9 / 2 ρ 2 l + 1 2 F 1 4 , 7 2 l ; 1 2 l ; ρ 2 = ( ρ 2 + 1 ) l 5 / 2 ρ 2 l + 1 × ρ 6 1 + 120 2 l 5 120 2 l 3 + 24 2 l 1 + 3 ρ 4 1 + 40 2 l 3 24 2 l 1 + 3 ρ 2 1 + 8 2 l 1 + 1 ,
v 5 l ( ρ ) = ( ρ 2 + 1 ) l + 9 / 2 2 F 1 4 , 9 2 + l ; 3 2 + l ; ρ 2 = ( ρ 2 + 1 ) l 5 / 2 × ρ 6 1 24 2 l + 3 + 120 2 l + 5 120 2 l + 7 + 3 ρ 4 1 + 24 2 l + 3 40 2 l + 5 + 3 ρ 2 1 8 2 l + 3 + 1 .
The method of variation of parameters enables us to obtain the particular solution of the inhomogeneous differential Equation (62) in the form
σ l ( p ) ( ρ ) = v 5 l ( ρ ) u 5 l ( ρ ) h l ( ρ ) d ρ ( ρ 2 + 1 ) 2 W l ( ρ ) u 5 l ( ρ ) v 5 l ( ρ ) h l ( ρ ) d ρ ( ρ 2 + 1 ) 2 W l ( ρ ) ,
where
W l ( ρ ) = 2 l + 1 ρ ρ 2 + 1 ρ 2 l + 1 .
Note that due to different representations for the function F l ( ρ ) (see Equation (57)) at values of ρ less and greater than 1, we obtain special representations for a particular solution in these two regions:
σ l ( 0 ) ( ρ ) = ( ρ 2 + 1 ) l 5 / 2 2 l + 1 ( 2 l 1 ) ( 2 l 3 ) ( 2 l 3 ) ρ 4 4 ( l 2 ) ρ 2 4 l 2 + 4 l 27 3 ( 2 l 5 ) , 0 ρ 1
σ l ( 1 ) ( ρ ) = ρ 2 l 1 ( ρ 2 + 1 ) l 5 / 2 2 l + 1 ( 2 l + 3 ) ( 2 l + 5 ) 4 l 2 + 4 l 27 3 ( 2 l + 7 ) + 4 ( l + 3 ) ρ 2 ( 2 l + 5 ) ρ 4 . ρ 1
It can be verified that both functions (67) and (68) have no singularities on their domains of definition. On the other hand, function u 5 l ( ρ ) is singular at ρ = 0 , whereas v 5 l ( ρ ) is singular at ρ = . This means that one should search the general solution of Equation (62) in the form:
σ l ( ρ ) = σ l ( 0 ) ( ρ ) + c 5 l ( v ) v 5 l ( ρ ) , 0 ρ 1
σ l ( ρ ) = σ l ( 1 ) ( ρ ) + c 5 l ( u ) u 5 l ( ρ ) . ρ 1
Note that two coefficients c 5 l ( v ) and c 5 l ( u ) are presently undetermined. To calculate them, we need to find two equations relating these coefficients. The first equation is quite obvious. It follows from the condition that the representations (69) and (70) are coincident at the common point ρ = 1 , that is
σ l ( 0 ) ( 1 ) + c 5 l ( v ) v 5 l ( 1 ) = σ l ( 1 ) ( 1 ) + c 5 l ( u ) u 5 l ( 1 ) .
This relationship reduces to the first desired equation:
c 5 l ( u ) ( 2 l + 3 ) ( 2 l + 5 ) ( 2 l + 7 ) = c 5 l ( v ) ( 2 l 5 ) ( 2 l 3 ) ( 2 l 1 ) + 27 4 l 4 l 2 3 × 2 l + 1 .
It can be verified that h l ( 0 ) = h l ( ) = 2 l . It follows from these relations that σ l ( 0 ) ( 0 ) + c 5 l ( v ) v 5 l ( 0 ) = σ l ( 1 ) ( ) + c 5 l ( u ) u 5 l ( ) . It can be assumed that the last equation represents the second desired equation. However, this assumption turns out to be false, because it again leads to Equation (72).
We propose the following method to find the second desired equation. Recall that any suitable function of the angles α and θ may be expanded into HHs since they form a complete set:
f ( α , θ ) = n = 0 ( 2 ) l = 0 n / 2 f n , l Y n , l ( α , θ ) ,
where (see, e.g., [3])
f n , l = f ( α , θ ) Y n , l ( α , θ ) d Ω
with
d Ω = π 2 sin 2 α sin θ d α d θ . α [ 0 , π ] , θ [ 0 , π ]
For the function f ( α , θ ) = f 2 ( α , θ ) represented by Equation (61), the expansion coefficient with n = 2 l becomes
f 2 l , l = π 2 0 π 0 π f 2 ( α , θ ) Y 2 l , l ( α , θ ) sin 2 α sin θ d α d θ = 2 π 2 N 2 l , l 2 l + 1 K 0 ( l ) + c l ( v ) K v ( l ) + c l ( u ) K u ( l ) ,
where
K 0 ( l ) = 0 π / 2 ( sin α ) 2 l + 2 σ l ( 0 ) ( ρ ) d α + π / 2 π ( sin α ) 2 l + 2 σ l ( ) ( ρ ) d α = = 2 24 l 3 100 l 2 + 198 l + 249 ( 2 l 5 ) ( 2 l 3 ) ( 2 l + 3 ) ( 2 l + 5 ) ( 2 l + 7 ) ( 2 l + 9 ) ,
K v ( l ) = 0 π / 2 ( sin α ) 2 l + 2 v 5 l ( ρ ) d α = 2 l + 3 / 2 ( 2 l 1 ) ( 2 l + 5 ) ( 2 l + 9 ) ,
K u ( l ) = π / 2 π ( sin α ) 2 l + 2 u 5 l ( ρ ) d α = 2 l + 3 / 2 ( 2 l + 3 ) ( 2 l + 7 ) ( 2 l 5 ) ( 2 l 3 ) ( 2 l + 9 ) .
To derive the results (76)–(79), we use the representation Y 2 l , l ( α , θ ) = N 2 l , l sin l α P l ( cos θ ) for the particular case of the HHs, and the orthogonality property for the Legendre polynomials. It should be noted that the explicit form of the normalization constant N 2 l , l is not required.
On the other hand, expanding f 2 ( α , θ ) in HHs, and inserting this expansion into the LHS of the IFRR (55), we obtain
Λ 2 45 f 2 ( α , θ ) = n = 0 ( 2 ) l = 0 n / 2 f n , l n ( n + 4 ) 45 Y n , l ( α , θ ) .
To derive the last equation, we use the fact that Y n , l ( α , θ ) is an eigenfunction of the operator Λ 2 with an eigenvalue equal to n ( n + 4 ) , that is
Λ 2 Y n , l ( α , θ ) = n ( n + 4 ) Y n , l ( α , θ ) .
The HH expansion of the RHS of Equation (55) is
h 2 ( α , θ ) = n = 0 ( 2 ) l = 0 n / 2 h n , l Y n , l ( α , θ ) .
Hence,
f 2 l , l = h 2 l , l 4 l ( l + 2 ) 45 .
Using again Sack’s representation (56) and (57) and Equation (74), we obtain the expansion coefficient h 2 l , l in explicit form:
h 2 l , l = π 2 0 π 0 π h 2 ( α , θ ) Y 2 l , l ( α , θ ) sin 2 α sin θ d α d θ = 2 7 / 2 π 2 N 2 l , l ( 4 l 2 + 24 l + 19 ) ( 2 l + 1 ) ( 2 l + 3 ) ( 2 l + 5 ) ( 2 l + 7 ) .
Thus, inserting (84) into the RHS of Equation (83) and equating the result to the RHS of Equation (76) we obtain the desired second equation in the form:
c 5 l ( u ) ( 2 l + 3 ) ( 2 l + 5 ) ( 2 l + 7 ) = c 5 l ( v ) ( 2 l 5 ) ( 2 l 3 ) ( 2 l 1 ) 2 l 1 ( 2 l + 1 ) .
Solving the system of two linear Equations (72) and (85) gives the desired coefficients:
c 5 l ( u ) = 2 l 1 ( l + 4 ) ( 2 l 3 ) 3 ( 2 l + 3 ) ( 2 l + 5 ) ( 2 l + 7 ) , c 5 l ( v ) = 2 l 1 ( l 3 ) ( 2 l + 5 ) 3 ( 2 l 5 ) ( 2 l 3 ) ( 2 l 1 ) .
It should be noted that the method described above for calculating the coefficients c 5 l ( u ) and c 5 l ( v ) is very reliable, but quite complex. A much simpler method is based on the statement that the point ρ = 1 represents the match point for the functions defined by Equations (69) and (70). This means that not only these functions, but also their first (at least) derivatives must coincide at this point. Thus, the second required equation relating the coefficients c 5 l ( u ) and c 5 l ( v ) is:
d σ l ( 0 ) ( ρ ) d ρ | ρ = 1 + c 5 l ( v ) d v 5 l ( ρ ) d ρ | ρ = 1 = d σ l ( 1 ) ( ρ ) d ρ | ρ = 1 + c 5 l ( u ) d u 5 l ( ρ ) d ρ | ρ = 1 .
The solution of the system of two Equations (72) and (87) again gives the coefficients defined by Equation (86). Substituting these coefficients into the representations (69) and (70), we finally obtain:
f 2 ( α , θ ) = 1 6 l = 0 ζ l ( ρ ) P l ( cos θ ) ( 2 l 1 ) ( 2 l + 3 ) ,
where
ζ l ( ρ ) = χ l ( ρ ) , 0 ρ 1 χ l ( 1 / ρ ) , ρ 1
with
χ l ( ρ ) = ρ l ( ρ 2 + 1 ) 5 / 2 ( l 3 ) ( 2 l 1 ) ρ 6 2 l + 7 + 9 l ρ 4 9 ( l + 1 ) ρ 2 ( l + 4 ) ( 2 l + 3 ) 2 l 5 .
It is clear that only the function (69) is required for calculating the function χ l ( ρ ) . Thus, in fact, we need to calculate only one coefficient c 5 l ( v ) to define this function. In this regard, it is important to emphasize that the representation (89) reflects the fact that the WF of a two-electron atomic system must preserve its parity when interchanging the electrons. For the singlet S-states (which include the ground state) this means that the AFC and/or its component preserves its form (including the sign) under the transformation α π α . For the AFC-component f 2 ( α , θ ) , represented by the series expansion (61), this property corresponds (in terms of variable ρ ) to the relationship:
σ l ( 0 ) ( ρ 1 ) + c 5 l ( v ) v 5 l ( ρ 1 ) = σ l ( 1 ) ( ρ ) + c 5 l ( u ) u 5 l ( ρ ) .
The elimination of the RHSs between Equations (71) and (91) for ρ = 1 yields the identity, whereas the use of Equation (87) instead of Equation (71) yields the required equation:
d σ l ( 0 ) ( ρ ) d ρ | ρ = 1 + c 5 l ( v ) d v 5 l ( ρ ) d ρ | ρ = 1 = d σ l ( 0 ) ( ρ 1 ) d ρ | ρ = 1 + c 5 l ( v ) d v 5 l ( ρ 1 ) d ρ | ρ = 1 .
The solution of the last equation gives the coefficients c 5 l ( v ) presented by Equation (86). Note that the coefficient c 5 l ( u ) can then be calculated by the use of Equation (71) if needed.
In the general case, we cannot sum the infinite series (88) to obtain the function f 2 ( α , θ ) in an explicit closed form. However, this can be done for some special angles α and/or θ . For example, it is worth noting that the angles θ = 0 , π correspond to the collinear configuration [14] of the two-electron atomic system in question. For these cases we obtain
f 2 ( α , 0 ) = ± ( ρ 1 ) ( 12 ρ 4 13 ρ 3 88 ρ 2 13 ρ + 12 ) 90 ( ρ 2 + 1 ) 5 / 2 ,
f 2 ( α , π ) = ( ρ + 1 ) ( 12 ρ 4 + 13 ρ 3 88 ρ 2 + 13 ρ + 12 ) 90 ( ρ 2 + 1 ) 5 / 2 .
Sign “+” in Equation (93) corresponds to 0 α π / 2 ( 0 ρ 1 ), whereas “−” to π / 2 α π ( ρ 1 ). The list of special θ -angles can be supplemented with an intermediate angle θ = π / 2 :
f 2 α , π 2 = 2 ( ρ 4 3 ρ 2 + 1 ) 15 ( ρ 2 + 1 ) 2 .
It is worth noting that for the important cases of the nucleus–electron and electron–electron coalescence, Equation (88), respectively, reduces to:
f 2 0 , θ = 2 15 , f 2 π 2 , 0 = 0 .
To derive the results (93)–(96) we used the relationships:
P n ( 0 ) = π Γ 1 1 n 2 Γ 1 n 2 + 1 , P n ( 1 ) = 1 , P n ( 1 ) = ( 1 ) n ,
where Γ ( x ) is the gamma function.

3. Derivation of the Angular Fock Coefficient ψ 6 , 3 ( α , θ )

We start this section by considering the FRR (3) and (4) for k = 6 and p = 2 :
Λ 2 60 ψ 6 , 2 = 48 ψ 6 , 3 2 V ψ 5 , 2 + 2 E ψ 4 , 2 .
Next, let us expand each function in Equation (98) into HHs, using Equation (73). This gives
ψ k , p = n = 0 ( 2 ) l = 0 n / 2 c n l ( k p ) Y n , l ( α , θ ) ,
with { k , p } = { 6 , 3 } , { 6 , 2 } , { 4 , 2 } and
V ψ 5 , 2 = n = 0 ( 2 ) l = 0 n / 2 f n l Y n , l ( α , θ ) ,
where the dimensionless potential V is defined by Equation (5), whereas the expansion coefficient f n l can be calculated by the formula
f n l = V ψ 5 , 2 Y n , l ( α , θ ) d Ω ,
according to Equation (74).
It follows from Equation (4) that h k , k / 2 = 0 for even k. Using additionally Equation (81), we can conclude that the AFC ψ k , k / 2 (with even k) represents the linear combination of the HHs, Y k , l ( α , θ ) . Hence,
c n l ( 63 ) = 0 f o r n 6 ,
c n l ( 42 ) = 0 f o r n 4 .
Equating the coefficients for the HHs, Y 6 , l ( α , θ ) in both sides of Equation (98), we obtain:
0 = 48 c 6 l ( 63 ) 2 f 6 l .
Hence, (using additionally Equation (101)),
c 6 l ( 63 ) = 1 24 V ψ 5 , 2 Y 6 , l ( α , θ ) d Ω .
Note that the LHS of Equation (104) equals zero, because ( Δ 2 60 ) Y 6 , l = 0 as follows from Equation (81).
Thus, according to Equations (99) and (102), the AFC ψ 6 , 3 ( α , θ ) represents a linear combination of four HHs, Y 6 , l ( α , θ ) with l = 0 , 1 , 2 , 3 . The contribution of each HH is determined by the coefficient c 6 l ( 63 ) given by Equation (105). However, it is easy to prove that only the coefficients with odd values of l are nonzero for ψ 6 , 3 ( α , θ ) . Indeed, this has already been mentioned in Section 2.4 that the WF of a two-electron atom/ion must preserve its parity when interchanging the electrons. For the singlet S-states this means that only the HHs, which preserve the sign under transformation α π α , differ from zero in the expansion of the WF, and hence in the expansion of any AFC. In turn, it is easy to show that only Y n , l ( α , θ ) with even values of ( n / 2 l ) satisfy the above property. Hence, the AFC in question, becomes
ψ 6 , 3 ( α , θ ) = a 61 Y 6 , 1 ( α , θ ) + a 63 Y 6 , 3 ( α , θ ) ,
where we denoted a 6 l c 6 l ( 63 ) ( l = 1 , 3 ) for convenience and simplicity, and where the normalized HHs are
Y 6 , 1 ( α , θ ) = 2 sin α + 3 sin ( 3 α ) cos θ π 3 / 2 5 , Y 6 , 3 ( α , θ ) = 8 sin 3 α P 3 ( cos θ ) π 3 / 2 5 .
Using Formula (105) and taking into account the representations (5) and (18) for the dimensionless potential V and the AFC ψ 5 , 2 , respectively, we can represent the desired coefficients in the form
a 6 l = ( π 2 ) ( 5 π 14 ) 6480 I l , 4 Z 4 + I l , 3 Z 3 + I l , 2 Z 2 , ( l = 1 , 3 )
where
I l , 4 = 4 π 3 / 2 0 π 0 π f 3 ( α , θ ) + 2 f 4 ( α , θ ) η Y 6 , l ( α , θ ) sin α sin θ d α d θ ,
I l , 3 = 2 0 π 0 π 3 η 2 f 1 ( α , θ ) + f 2 ( α , θ ) sin α + π 3 / 2 f 3 ( α , θ ) + 2 f 4 ( α , θ ) ξ × × Y 6 , l ( α , θ ) sin 2 α sin θ d α d θ ,
I l , 2 = 3 0 π 0 π 2 f 1 ( α , θ ) + f 2 ( α , θ ) ξ 1 Y 6 , l ( α , θ ) sin 2 α sin θ d α d θ .
To calculate the integral (110), it is useful to separate the contributions which include the functions f 1 , f 3 and f 4 , represented by the explicit closed expressions, and the function f 2 , represented by the infinite series (88). We obtain:
I l , 3 = I l , 3 ( 134 ) 6 I l , 3 ( 2 ) ,
where
I l , 3 ( 134 ) = 2 0 π 0 π 6 η f 1 ( α , θ ) sin α + π 3 / 2 f 3 ( α , θ ) + 2 f 4 ( α , θ ) ξ Y 6 , l ( α , θ ) sin 2 α sin θ d α d θ ,
I l , 3 ( 2 ) = 0 π 0 π f 2 ( α , θ ) η Y 6 , l ( α , θ ) sin α sin θ d α d θ .
The integrals (113) can be taken in an explicit (closed) form that gives:
I 1 , 3 ( 134 ) = 3 ( 45 π 122 ) 35 π 3 / 2 5 , I 3 , 3 ( 134 ) = 245 π 816 70 π 3 / 2 5 .
The problem of calculating the integrals (114) is that the corresponding integrands contain the function f 2 ( α , θ ) represented by the infinite series (88). Fortunately, using the orthogonality relationship for the Legendre polynomials, we can get these integrals also in explicit form. Changing the order of summation and integration, we easily obtain:
I 1 , 3 ( 2 ) = π 3 / 2 3 5 l = 0 0 π [ sin α + 3 sin ( 3 α ) ] η ζ l ( ρ ) ( 2 l 1 ) ( 2 l + 3 ) sin α d α 0 π P l ( cos θ ) cos θ sin θ d θ = 2 π 3 / 2 45 5 0 π [ sin α + 3 sin ( 3 α ) ] η ζ 1 ( ρ ) sin α d α = 7 π + 22 210 π 3 / 2 5 ,
I 3 , 3 ( 2 ) = 4 π 3 / 2 3 5 l = 0 0 π η ζ l ( ρ ) ( 2 l 1 ) ( 2 l + 3 ) sin 4 α d α 0 π P l ( cos θ ) P 3 ( cos θ ) sin θ d θ = 8 π 3 / 2 945 5 0 π η ζ 3 ( ρ ) sin 4 α d α = 3 π 32 180 π 3 / 2 5 .
Recall that η η ( α ) is defined by Equation (6) and ρ = tan ( α / 2 ) .
It can be shown (using fairly long nontrivial derivations) that the integrals I l , 2 and I l , 4 vanish both for l = 1 and l = 3 . This means that (according to the representations (106) and (108)) the AFC, ψ 6 , 3 ( α , θ ) is proportional to the third power of the nucleus charge Z (only), which is in full agreement with Formula (13) from Ref. [4].
Thus, combining the results of this section, we obtain the nonzero coefficients a 6 , l in the simple final form:
a 61 = ( π 2 ) ( 5 π 14 ) ( 32 π 97 ) 56700 π 3 / 2 5 Z 3 ,
a 63 = ( π 2 ) ( 5 π 14 ) ( 357 π 1112 ) 680400 π 3 / 2 5 Z 3 .

4. Derivation of the Angular Fock Coefficients ψ 7 , 3 ( α , θ ) and ψ 8 , 4 ( α , θ )

In Section 2 and Section 3 we detailed the derivation of the AFCs ψ 5 , 2 ( α , θ ) and ψ 6 , 3 ( α , θ ) , respectively. Therefore, for the AFCs ψ 7 , 3 ( α , θ ) and ψ 8 , 4 ( α , θ ) , we give only abbreviated derivations, and include extended explanations only in case of significant differences.

4.1. The AFC ψ 7 , 3 ( α , θ )

The FRR (3) and (4) for k = 7 and p = 3 reduces to the form
Λ 2 77 ψ 7 , 3 ( α , θ ) = h 7 , 3 ( α , θ ) ,
where
h 7 , 3 ( α , θ ) = 2 V ψ 6 , 3 ( α , θ ) .
Using Equations (106), (107) and (5) the RHS of Equation (120) can be represented in the form:
h 7 , 3 ( α , θ ) = ( π 2 ) ( 5 π 14 ) Z 3 340200 5 π 3 / 2 h ¯ 1 + h ¯ 2 5 π 3 / 2 2 Z 12 ( 32 π 97 ) h ¯ 3 + ( 357 π 1112 ) h ¯ 4 ,
where
h ¯ 1 = 20 ξ 1 12 ( 32 π 97 ) ( 1 ξ 2 ) + ( 357 π 1112 ) ( 1 ξ 2 ) 3 ,
h ¯ 2 = 60 ( 688 255 π ) ξ 1 sin 3 α cos θ ,
h ¯ 3 = η ( sin α ) 1 Y 6 , 1 ( α , θ ) , h ¯ 4 = η ( sin α ) 1 Y 6 , 3 ( α , θ ) .
Accordingly, the solution of the FRR (120) can be found in the form:
ψ 7 , 3 ( α , θ ) = ( π 2 ) ( 5 π 14 ) Z 3 340200 5 π 3 / 2 f ¯ 1 + f ¯ 2 5 π 3 / 2 2 Z 12 ( 32 π 97 ) f ¯ 3 + ( 357 π 1112 ) f ¯ 4 ,
where the AFC components f ¯ i satisfy the IFRRs
Λ 2 77 f ¯ i = h ¯ i . ( i = 1 , 2 , 3 , 4 )
Note that the components h ¯ i of the RHS h 7 , 3 of the FRR (120) for the AFC ψ 7 , 3 are reasonably close to the components h i of the RHS h 5 , 2 of the FRR (14). Therefore, we only briefly dwell on the conclusions of the corresponding results, as we noted earlier.
It is seen from Equation (123) that the RHS h ¯ 1 is a function of a single variable ξ defined by Equation (6). The solution of the corresponding IFRR was described in Section IV of Ref. [4] (see also Section II of Ref. [12]) and illustrated (among others) in Section 2.3 of the current article. Thus, following the technique mentioned above, we obtain:
f ¯ 1 = 41437 π 12 74342 7 ξ 7 + 36476 35588 π 3 ξ 5 + + 5 2 ( 4931 π 15156 ) ξ 3 + 5 ( 2276 741 π ) ξ .
It can be verified that the RHSs h ¯ 3 and h ¯ 4 represent functions of the form f ( α ) P l ( cos θ ) with l equals 1 and 3, respectively. The solution of the corresponding IFRR was described in Section V of Ref. [4] and illustrated in Section 2.1 and Section 2.2 of the current article. This enables us to obtain:
f ¯ 3 = ρ ( 1 + ρ ) 29 + ρ 16 + ρ [ ρ ( 16 + 29 ρ ) 114 ] cos θ 9 5 π 3 / 2 ( ρ 2 + 1 ) 7 / 2 ,
f ¯ 4 = sin 3 α 1 + sin α 2 5 π 3 / 2 P 3 ( cos θ ) .
Recall that the ρ variable was defined previously in Section 2.1.
The RHS h ¯ 2 represented by Equation (124) is slightly more complicated than h 2 discussed in Section 2.4. In this regard, it would be useful to clarify two points.
First, using representation (56) for ξ 1 , we can rewrite Equation (124) in the form:
h ¯ 2 = 60 ( 688 225 π ) h ˜ 2
where
h ˜ 2 = l = 0 2 l ( sin α ) l + 3 F l ( ρ ) cos θ P l ( cos θ ) ,
and where F l ( ρ ) is defined by Equations (57) and (58). In order to apply the solution of the corresponding IFRR by the method described in Section 2.4 (see also [4]), the θ -dependent l-component in the series expansion of h ˜ 2 must be pure P l ( cos θ ) . To solve the problem, one could use the general formula representing the Clebsch–Gordan series for a product of two spherical harmonics. However, in our simple case, it is easier to use the recurrence relation for the Legendre polynomials
( l + 1 ) P l + 1 ( x ) ( 2 l + 1 ) x P l ( x ) + l P l 1 ( x ) = 0 ,
which enables us to represent h ˜ 2 in the desired form:
h ˜ 2 = l = 0 h ¯ l ( ρ ) ( sin α ) l P l ( cos θ ) ,
where
h ¯ l ( ρ ) = l 2 l 1 ( 2 l 1 ) sin 2 α F l 1 ( ρ ) + l + 1 2 l + 1 ( 2 l + 3 ) sin 4 α F l + 1 ( ρ ) .
The second point is related to the calculation of the coefficient
h ˜ 2 l , l = π 2 0 π 0 π h ˜ 2 ( α , θ ) Y 2 l , l ( α , θ ) sin 2 α sin θ d α d θ
in the HH expansion of h ˜ 2 (see the corresponding Equation (84) for calculation of ψ 5 , 2 ( α , θ ) ). Of course, we can use representation (134) and (135) and then apply the orthogonality condition for the Legendre polynomials. However, the simpler way is to use the original representation (132) taking into account that cos θ P 1 ( cos θ ) . In this case, we can apply the well-known formula for the integral of three Legendre polynomials
1 1 P l ( x ) P L ( x ) P l ( x ) d x = 2 l L l 0 0 0 2 ,
where the RHS represents twice the square of the Wigner 3-j symbol.
Thus, applying the methodologies outlined in Section 2.4, and given the above features, one obtains
f ¯ 2 = 1 48 l = 0 ζ ¯ l ( ρ ) P l ( cos θ ) ( 2 l 1 ) ( 2 l + 3 ) ,
where
ζ ¯ l ( ρ ) = χ ¯ l ( ρ ) , 0 ρ 1 χ ¯ l ( 1 / ρ ) , ρ 1
with
χ ¯ l ( ρ ) = ρ l ( ρ 2 + 1 ) 7 / 2 { ( 32 l 2 + 26 l 25 ) ρ 6 2 l + 5 ( 2 l 1 ) ρ 2 2 l + 9 + 4 + + 1 2 l 3 6 ( 84 l 2 + 84 l 95 ) ρ 4 2 l + 5 ( 32 l 2 + 38 l 19 ) 2 l + 3 2 l 7 + 4 ρ 2 } .
Recall that the component f ˘ 2 in the RHS of Equation (126) is equal to 60 ( 688 225 π ) f ¯ 2 according to representation (131).
As in the case of the AFC ψ 5 , 2 ( α , θ ) , there are combinations of special hyperspherical angles α and θ for which the component f ¯ 2 f ¯ 2 ( α , θ ) of the AFC ψ 7 , 3 ( α , θ ) can be obtained in closed form. In particular, one obtains:
f ¯ 2 ( α , 0 ) = ( ρ 1 ) ( 95 ρ 6 + 1166 ρ 5 1879 ρ 4 8844 ρ 3 1879 ρ 2 + 1166 ρ + 95 ) 5040 ( ρ 2 + 1 ) 7 / 2 ,
f ¯ 2 ( α , π ) = ( ρ + 1 ) ( 95 ρ 6 1166 ρ 5 1879 ρ 4 + 8844 ρ 3 1879 ρ 2 1166 ρ + 95 ) 5040 ( ρ 2 + 1 ) 7 / 2 ,
f ¯ 2 ( α , π 2 ) = 19 ρ 4 + 10 ρ 2 + 19 1008 ( ρ 2 + 1 ) 2 .
Sign “−” in Equation (141) corresponds to 0 α π / 2 ( 0 ρ 1 ), whereas “+” corresponds to π / 2 α π ( ρ 1 ).
For the important cases of the nucleus–electron and electron–electron coalescence, representation (138)–(140) is simplified to:
f ¯ 2 0 , θ = 19 1008 , f ¯ 2 π 2 , 0 = 0 .

4.2. The AFC ψ 8 , 4 ( α , θ )

Having at our disposal the AFC ψ 7 , 3 ψ 7 , 3 ( α , θ ) , we can calculate the AFC ψ 8 , 4 ψ 8 , 4 ( α , θ ) using the FRR (3) and (4) for k = 8 and p = 3 :
Λ 2 96 ψ 8 , 3 = 80 ψ 8 , 4 2 V ψ 7 , 3 + 2 E ψ 6 , 3 .
It follows from Equation (81) and the FRR (3) and (4) for k = 8 and p = 4 that the AFC ψ 8 , 4 is a linear combination of the HHs, Y 8 , l Y 8 , l ( α , θ ) . Moreover, given that only Y n , l ( α , θ ) with even values of n / 2 l are suitable for singlet S-states, we obtain:
ψ 8 , 4 = a 80 Y 8 , 0 + a 82 Y 8 , 2 + a 84 Y 8 , 4 .
For further derivations, it is advisable to represent the HHs in the form
Y 8 , l ( α , θ ) = y 8 l ( α ) P l ( cos θ ) ,
where
y 80 ( α ) = π 3 / 2 2 cos ( 4 α ) + 2 cos ( 2 α ) + 1 ,
y 82 ( α ) = 2 π 3 / 2 10 7 sin 2 α 4 cos ( 2 α ) + 3 ,
y 84 ( α ) = 8 π 3 / 2 2 7 sin 4 α .
It was found in Section 3 that ψ 6 , 3 ψ 6 , 3 ( α , θ ) is the linear combination of the HHs Y 6 , l ( α , θ ) . Thus, expanding each function of Equation (145) in HHs, and equating the coefficients for Y 8 , l , we obtain (see the corresponding result (105) for a 6 l )
a 8 l = 1 40 V ψ 7 , 3 Y 8 , l d Ω ,
where the potential V is defined by Equation (5). When deriving the last equation, it was taken into account that ( Δ 2 96 ) Y 8 , l = 0 , as follows from Equation (81).
A direct substitution of the representations (5), (126) and (147) into the RHS of Equation (151) yields:
a 8 l = ( π 2 ) ( 5 π 14 ) Z 3 13608000 5 π 3 / 2 × × 1 ξ 2 Z η sin α f ˘ 1 + f ˘ 2 5 π 3 / 2 2 Z 12 ( 32 π 97 ) f ˘ 3 + ( 357 π 1112 ) f ˘ 4 Y 8 , l ( α , θ ) d Ω .
It follows from Equation (13) of Ref. [4] that only the coefficients at Z 4 are nonzero on the RHS of the last equation. Hence, Equation (152) reduces to the form:
a 8 l = ( π 2 ) ( 5 π 14 ) Z 4 6804000 5 π 3 / 2 S 1 l + S 2 l 5 π 3 / 2 + 12 ( 32 π 97 ) S 3 l + ( 357 π 1112 ) S 4 l ,
where
S 1 l = π 2 0 π 0 π f ˘ 1 ( ξ ) Y 8 , l ( α , θ ) η sin α sin θ d α d θ ,
S 2 l = π 2 0 π 0 π f ˘ 2 ( α , θ ) Y 8 , l ( α , θ ) η sin α sin θ d α d θ = = 5 π 2 ( 688 225 π ) ( 2 l 1 ) ( 2 l + 1 ) ( 2 l + 3 ) 0 π / 2 χ ¯ l ( ρ ) y 8 l ( α ) η sin α d α ,
S n l = π 2 0 π 0 π f ˘ n ( α , θ ) Y 8 , l ( α , θ ) ξ 1 sin 2 α sin θ d α d θ . ( n = 3 , 4 )
The identifiers ξ and η are defined by Equation (6), whereas functions χ ¯ l ( ρ ) can be calculated by Formula (140). When deriving Equation (155), we applied the orthogonality condition for the Legendre polynomials. Fortunately, all integrals (154)–(156) can be taken in closed form. Thus, by collecting these results and substituting them into the RHS of Equation (153), we finally obtain the desired coefficients in the form:
a 8 l = Z 4 ( π 2 ) ( 5 π 14 ) π 5 / 2 b 8 l ,
with
b 80 = π ( 150339 π 927292 ) + 1430792 19289340000 , b 82 = π ( 751965 π 4654046 ) + 7200976 1928934000 70 , b 84 = π ( 3190317 π 19828996 ) + 30802176 25719120000 14 .

5. Results and Discussions

The angular Fock coefficients ψ k , p ψ k , p ( α , θ ) with the maximum possible value of subscript p were calculated on examples of the coefficients with 5 k 10 . The results obtained in Section 2, Section 3 and Section 4 are summarized in Appendix A and Appendix B. The AFCs ψ 9 , 4 and ψ 10 , 5 are presented in Appendix C without derivations. To find the latter AFCs, the methods described in the main sections were used. The presented technique makes it possible to calculate such AFCs for any arbitrarily large k. These coefficients are leading in the logarithmic power series representing the Fock expansion (see Equation (8)). As such, they may be indispensable for the development of simple methods for calculating the helium-like electronic structure.
The proposed technique, as well as the final results, are quite complex. Therefore, both require verification. We are aware of two ways for the above-mentioned verification. The first one is to use the Green’s function (GF) approach (see Ref. [1] and also Ref. [15], Section 4) which enables us to calculate (at least, numerically) the AFC (or its component) by the following integral representation:
ψ k , p ( α , θ ) = 1 8 π 0 π d α sin 2 α 0 π d θ sin θ h k , p ( α , θ ) 0 π cos k 2 + 1 ω sin ω ( 1 λ ) d φ ,
where ω is an angle defined by the relation
cos ω = cos α cos α + sin α sin α cos θ cos θ + sin θ sin θ cos φ ,
whereas
λ = 0 k o d d ω / π k e v e n .
For even k and maximum value of p = k / 2 , the RHS h k , k / 2 of the FRR (3) equals zero. This implies that the GF formula (159) cannot be applied in this case. Hence, only the AFCs ψ k , p with odd values of k (and maximum p) can be verified with the GF method. Thus, numerically calculating (for various combinations of angles α and θ ) the triple integrals (159) representing the AFCs ψ 5 , 2 ( α , θ ) , ψ 7 , 3 ( α , θ ) and ψ 9 , 4 ( α , θ ) , we verified that the representations obtained for them in Section 2 and Section 4.1 and in Appendix C were correct.
The second verification method considered, covering all possible combinations of angles, being quite complex, is the only method known to us that correctly displays the WF near the nucleus. This is the CFHH method mentioned in the Introduction. It is based on decomposing the full WF into a form
Ψ CFHH ( r 1 , r 2 , r 12 ) = exp f ( r 1 , r 2 , r 12 ) Φ CFHH ( R , α , θ ) ,
where the so-called correlation function f can be taken in a simple linear form
f ( r 1 , r 2 , r 12 ) = c 1 r 1 + c 2 r 2 + c 12 r 12 .
The so-called “cusp parametrization”
c 1 = c 2 = Z , c 12 = 1 / 2
is used as a rule. For a small enough hyperspherical radius R, the function Φ is represented as
Φ CFHH ( R , α , θ ) = 1 d 0 , 0 ( α , θ ) k = 0 K ( 2 κ R ) k p = 0 [ k / 2 ] d k , p ( α , θ ) ln p ( 2 κ R ) ,
where κ = 2 E , and functions d k , p ( α , θ ) are expanded in N (basis size) HHs. It follows from representation (165) that the AFCs ψ k , p ( α , θ ) can be expressed in terms of the functions d k , p ( α , θ ) calculated by the CFHHM. For example, for the AFCs in question, one obtains:
ψ k , k / 2 CFHH ( α , θ ) = ( 2 κ ) k d k , k / 2 ( α , θ ) d 0 , 0 ( α , θ ) .
We calculated all AFCs discussed in this article using CFHHM with K = 18 and N = 1600 . The angles 0 α π and 0 θ π with step π / 6 were considered. The relative difference | 1 ψ k , p ( α , θ ) / ψ k , p CFHH ( α , θ ) | was less than 10 7 for all examined cases, including 1 Z 5 . This indicates that all our theoretical calculations were correct.

Author Contributions

Conceptualization, E.Z.L.; Methodology, E.Z.L.; Software, E.Z.L. and R.K.; Validation, R.K.; Writing—original draft, E.Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Summarizing the results of Section 2, one obtains.
ψ 5 , 2 ( α , θ ) = Z 2 ( π 2 ) ( 5 π 14 ) 270 π 3 π 3 / 2 ( 2 f 1 + f 2 ) 2 Z ( f 3 + 2 f 4 ) ,
where
f 1 = 1 60 1 sin α cos θ sin α cos θ ( 4 + 13 sin α cos θ ) 2 ,
f 3 ( α ) = 1 60 π 3 / 2 11 sin α + 21 cos ( 2 α ) + 2 1 + sin α ,
f 4 ( α , θ ) = 2 6 π 3 / 2 ( sin α ) 2 1 + sin α P 2 ( cos θ ) ,
f 2 ( α , θ ) = 1 6 l = 0 ζ l ( ρ ) P l ( cos θ ) ( 2 l 1 ) ( 2 l + 3 ) .
The ζ function is defined as follows:
ζ l ( ρ ) = χ l ( ρ ) , 0 ρ 1 χ l ( 1 / ρ ) , ρ 1
where
χ l ( ρ ) = ρ l ( ρ 2 + 1 ) 5 / 2 ( l 3 ) ( 2 l 1 ) ρ 6 2 l + 7 + 9 l ρ 4 9 ( l + 1 ) ρ 2 ( l + 4 ) ( 2 l + 3 ) 2 l 5 .
Recall that ρ = tan ( α / 2 ) , and special cases of the function f 2 f 2 ( α , θ ) , when they can be obtained in closed form, are represented by Equations (93)–(96).
Summarizing the results of Section 3, one obtains:
ψ 6 , 3 ( α , θ ) = Z 3 ( π 2 ) ( 5 π 14 ) 56700 π 3 / 2 5 ( 97 32 π ) Y 6 , 1 ( α , θ ) + ( 1112 357 π ) 12 Y 6 , 3 ( α , θ ) ,
where Y n , l ( α , θ ) are the normalized hyperspherical harmonics.

Appendix B

Summarizing the results of Section 4.1, one obtains:
ψ 7 , 3 ( α , θ ) = ( π 2 ) ( 5 π 14 ) Z 3 340200 5 π 3 / 2 f ¯ 1 + f ¯ 2 5 π 3 / 2 2 Z 12 ( 32 π 97 ) f ¯ 3 + ( 357 π 1112 ) f ¯ 4 ,
where
f ¯ 1 = 41437 π 12 74342 7 ξ 7 + 36476 35588 π 3 ξ 5 + + 5 2 ( 4931 π 15156 ) ξ 3 + 5 ( 2276 741 π ) ξ .
f ¯ 3 = ρ ( 1 + ρ ) 29 + ρ 16 + ρ [ ρ ( 16 + 29 ρ ) 114 ] cos θ 9 5 π 3 / 2 ( ρ 2 + 1 ) 7 / 2 ,
f ¯ 4 = sin 3 α 1 + sin α 2 5 π 3 / 2 P 3 ( cos θ ) ,
f ¯ 2 = 1 48 l = 0 ζ ¯ l ( ρ ) P l ( cos θ ) ( 2 l 1 ) ( 2 l + 3 ) .
The ζ ¯ function is defined as follows:
ζ ¯ l ( ρ ) = χ ¯ l ( ρ ) , 0 ρ 1 χ ¯ l ( 1 / ρ ) , ρ 1 ,
where
χ ¯ l ( ρ ) = ρ l ( ρ 2 + 1 ) 7 / 2 { ( 32 l 2 + 26 l 25 ) ρ 6 2 l + 5 ( 2 l 1 ) ρ 2 2 l + 9 + 4 + + 1 2 l 3 6 ( 84 l 2 + 84 l 95 ) ρ 4 2 l + 5 ( 32 l 2 + 38 l 19 ) 2 l + 3 2 l 7 + 4 ρ 2 } .
Recall that variable ξ is defined by Equation (6), and special cases of the function f ¯ 2 f ¯ 2 ( α , θ ) , when they can be obtained in closed form, are represented by Equations (141)–(144).
Summarizing the results of Section 4.2, one obtains:
ψ 8 , 4 = Z 4 ( π 2 ) ( 5 π 14 ) π 5 / 2 b 80 Y 8 , 0 ( α , θ ) + b 82 Y 8 , 2 ( α , θ ) + b 84 Y 8 , 4 ( α , θ ) ,
where
b 80 = π ( 150339 π 927292 ) + 1430792 19289340000 , b 82 = π ( 751965 π 4654046 ) + 7200976 1928934000 70 , b 84 = π ( 3190317 π 19828996 ) + 30802176 25719120000 14 .

Appendix C

In Section 2 and Section 3, the AFCs ψ 5 , 2 ( α , θ ) and ψ 6 , 3 ( α , θ ) were calculated with detailed derivations. In Section 4, the AFCs ψ 7 , 3 ( α , θ ) and ψ 8 , 4 ( α , θ ) were presented with a very brief derivations. The corresponding results were summarized in Appendix A and Appendix B. The current Appendix presents the AFCs ψ 9 , 4 ψ 9 , 4 ( α , θ ) and ψ 10 , 5 ψ 10 , 5 ( α , θ ) without derivations. To find the latter AFCs, the methods described in the main sections were used.
So, the first AFC under consideration can be represented as:
ψ 9 , 4 = 2 Z 4 2 Z X 1 ( α , θ ) X 2 ( α , θ ) ,
where
X 1 ( α , θ ) = a ˇ 80 f ˇ 1 + a ˇ 82 f ˇ 2 + a ˇ 84 f ˇ 3 ,
X 2 ( α , θ ) = 35 π 3 / 2 2 7 a ˇ 84 f ˇ 4 + ( π 2 ) ( 5 π 14 ) 123451776000 π 4 c 5 f ˇ 5 + c 6 f ˇ 6 + c 7 f ˇ 7 + 16 ( c 8 f ˇ 8 + c 9 f ˇ 9 ) .
Here, a ˇ 8 l = Z 4 a 8 l , where the coefficients a 8 l are defined by Equations (157) and (158), whereas the other coefficients are:
c 5 = π ( 29757524 4780401 π ) 46286848 , c 6 = π ( 9581100 π 59458928 ) + 92239360 , c 7 = π ( 28060 + 10149 π ) 167168 , c 8 = 9 π ( 134543 π 828732 ) + 11488128 , c 9 = π ( 4804833 π 29773780 ) + 46119680 .
The functions f ˇ i f ˇ i ( α , θ ) are:
f ˇ 1 = ( ρ + 1 ) ( 563 ρ 8 + 1012 ρ 7 8932 ρ 6 3668 ρ 5 + 23954 ρ 4 3668 ρ 3 8932 ρ 2 + 1012 ρ + 563 ) 1260 π 3 / 2 ( ρ 2 + 1 ) 9 / 2 ,
f ˇ 2 = 8 π 3 / 2 10 7 ρ 2 ( 1 + ρ ) ( 126 + 49 ρ 424 ρ 2 + 49 ρ 3 + 126 ρ 4 ) 300 ( ρ 2 + 1 ) 9 / 2 P 2 ( cos θ ) ,
f ˇ 3 = 2 5 π 3 / 2 2 ( 1 + sin α ) 7 sin 4 α P 2 ( cos θ ) ,
f ˇ 4 = ξ ( 315 1680 ξ 2 + 2814 ξ 4 1854 ξ 6 + 419 ξ 8 ) 1260 ,
f ˇ 5 = ξ 60 ( 2 ξ 2 3 ) ( 2 ξ 2 1 ) ( 4 ξ 4 10 ξ 2 + 5 ) .
The remaining f ˇ functions are represented by series:
f ˇ j = ( ρ 2 + 1 ) 9 / 2 k j l = 0 ρ l ζ ˇ j l ( ρ ) ( 2 l 1 ) ( 2 l + 3 ) P l ( cos θ ) , ( j = 6 , 7 , 8 , 9 )
where
k 6 = 6 , k 7 = 60 , k 8 = 24 , k 9 = 40 ,
and the corresponding ζ ˇ functions are:
ζ ˇ 6 l ( ρ ) = ( 2 l 15 ) ( 2 l 1 ) ( l + 1 ) ρ 10 ( 2 l + 7 ) ( 2 l + 11 ) + ( 22 l 2 5 l 12 ) ρ 8 ( 2 l + 7 ) + 10 ( 2 l 2 + 11 l + 3 ) ρ 6 2 l + 7 10 ( 2 l 2 7 l 6 ) ρ 4 2 l 5 ( 22 l 2 + 49 l + 15 ) ρ 2 2 l 5 l ( 2 l + 3 ) ( 2 l + 17 ) ( 2 l 9 ) ( 2 l 5 ) ,
ζ ˇ 7 l ( ρ ) = ( 2 l 1 ) ( 4 l 2 + 160 l 189 ) ρ 10 ( 2 l + 7 ) ( 2 l + 11 ) + 35 ( 4 l 2 + 40 l 9 ) ρ 8 ( 2 l + 7 ) 350 ( 4 l 2 + 16 l + 3 ) ρ 6 2 l + 7 + + 350 ( 4 l 2 8 l 9 ) ρ 4 2 l 5 35 ( 4 l 2 32 l 45 ) ρ 2 2 l 5 ( 2 l + 3 ) ( 4 l 2 152 l 345 ) ( 2 l 9 ) ( 2 l 5 ) ,
ζ ˇ 8 l ( ρ ) = ( 2 l 1 ) ( 56 l 3 + 250 l 2 + 338 l + 171 ) ρ 10 ( 2 l + 5 ) ( 2 l + 7 ) ( 2 l + 11 ) ( 136 l 3 + 314 l 2 110 l 153 ) ρ 8 ( 2 l + 5 ) ( 2 l + 7 ) 2 ( 80 l 4 + 652 l 3 + 566 l 2 1824 l 873 ) ρ 6 ( 2 l 3 ) ( 2 l + 5 ) ( 2 l + 7 ) + 2 ( 80 l 4 332 l 3 910 l 2 + 1320 l + 945 ) ρ 4 ( 2 l 5 ) ( 2 l 3 ) ( 2 l + 5 ) + + ( 136 l 3 + 94 l 2 330 l 135 ) ρ 2 ( 2 l 5 ) ( 2 l 3 ) + ( 2 l + 3 ) ( 56 l 3 82 l 2 + 6 l 27 ) ( 2 l 9 ) ( 2 l 5 ) ( 2 l 3 ) ,
ζ ˇ 9 l ( ρ ) = ( 2 l 1 ) ( 24 l 3 150 l 2 670 l 439 ) ρ 10 ( 2 l + 5 ) ( 2 l + 7 ) ( 2 l + 11 ) + 5 ( 72 l 3 + 162 l 2 70 l 103 ) ρ 8 ( 2 l + 5 ) ( 2 l + 7 ) + + 10 ( 16 l 4 + 220 l 3 + 222 l 2 804 l 423 ) ρ 6 ( 2 l 3 ) ( 2 l + 5 ) ( 2 l + 7 ) 10 ( 16 l 4 156 l 3 342 l 2 + 652 l + 399 ) ρ 4 ( 2 l 5 ) ( 2 l 3 ) ( 2 l + 5 ) 5 ( 72 l 3 + 54 l 2 178 l 57 ) ρ 2 ( 2 l 5 ) ( 2 l 3 ) ( 2 l + 3 ) ( 24 l 3 + 222 l 2 298 l 57 ) ( 2 l 9 ) ( 2 l 5 ) ( 2 l 3 ) ,
It is important to emphasize that the representations (A27)–(A32) are valid only for 0 ρ 1 . For values ρ > 1 , one should replace ρ with 1 / ρ , which is equivalent to simply redefining ρ as cot ( α / 2 ) .
The second AFC under consideration is of the form:
ψ 10 , 5 = Z 5 ( π 2 ) ( 5 π 14 ) π 7 / 2 b 10 , 1 Y 10 , 1 ( α , θ ) + b 10 , 3 Y 10 , 3 ( α , θ ) + b 10 , 5 Y 10 , 5 ( α , θ ) ,
where
b 10 , 1 = π [ 3 π ( 6840010557 π 63828704998 ) + 595609133656 ] 617517605744 401025378600000 105 ,
b 10 , 3 = π [ π ( 9194460432 π 85833963053 ) + 267084629592 ] 277009842768 100256344650000 30 ,
b 10 , 5 = π [ π ( 622341848670 π 5812646794643 ) + 18095537797140 ] 18776793358080 10025634465000000 42 ,
and Y 10 , l ( α , θ ) with l = 1 , 3 , 5 are the normalized HHs.

References

  1. Fock, V.A. On the Schrödinger Equation of the Helium Atom. Izv. Akad. Nauk SSSR Ser. Fiz. 1954, 18, 161–174. [Google Scholar]
  2. Morgan, J.D., III. Convergence properties of Fock’s expansion for S-state eigenfunctions of the helium atom. Theor. Chim. Acta 1986, 69, 181–223. [Google Scholar] [CrossRef]
  3. Abbott, P.C.; Maslen, E.N. Coordinate systems and analytic expansions for three-body atomic wavefunctions: I. Partial summation for the Fock expansion in hyperspherical coordinates. J. Phys. A Math. Gen. 1987, 20, 2043–2075. [Google Scholar] [CrossRef]
  4. Liverts, E.Z.; Barnea, N. Angular Fock coefficients. Refinement and further development. Phys. Rev. A 2015, 92, 042512. [Google Scholar] [CrossRef]
  5. Nakashima, H.; Nakatsuji, H. Solving the Schrödinger equation for helium atom and its isoelectronic ions with the free iterative complement interaction (ICI) method. J. Chem. Phys. 2007, 127, 224104. [Google Scholar] [CrossRef] [Green Version]
  6. Rodriguez, K.V.; Gasaneo, G.; Mitnik, D.M. Accurate and simple wavefunctions for the helium isoelectronic sequence with correct cusp conditions. J. Phys. B 2007, 40, 3923–3939. [Google Scholar] [CrossRef]
  7. Forrey, R.C. Compact representation of helium wave functions in perimetric and hyperspherical coordinates. Phys. Rev. A 2004, 69, 022504. [Google Scholar] [CrossRef] [Green Version]
  8. Drake, G.W.F. (Ed.) High Precision Calculations for Helium. In Atomic, Molecular, and Optical Physics Handbook; AIP Press: New York, NY, USA, 1996; Section 11. [Google Scholar]
  9. Haftel, M.I.; Mandelzweig, V.B. Exact Solution of Coupled Equations and the Hyperspherical Formalism: Calculation of Expectation Values and Wavefunctions of Three Coulomb-Bound Particles. Ann. Phys. 1983, 150, 48–91. [Google Scholar] [CrossRef]
  10. Haftel, M.I.; Mandelzweig, V.B. Fast Convergent Hyperspherical Harmonic Expansion for Three-Body Systems. Ann. Phys. 1989, 189, 29–52. [Google Scholar] [CrossRef]
  11. Haftel, M.I.; Krivec, R.; Mandelzweig, V.B. Power Series Solution of Coupled Differential Equations in One Variable. J. Comp. Phys. 1996, 123, 149–161. [Google Scholar] [CrossRef] [Green Version]
  12. Liverts, E.Z. Analytic calculation of the edge components of the angular Fock coefficients. Phys. Rev. A 2016, 94, 022504. [Google Scholar] [CrossRef]
  13. Sack, R.A. Generalization of Laplace’s expansion to arbitrary powers and functions of the distance between two points. J. Math. Phys. 1964, 5, 245–251. [Google Scholar] [CrossRef]
  14. Liverts, E.Z.; Krivec, R.; Barnea, N. Collinear configuration of the helium atom two-electron ions. Ann. Phys. 2020, 422, 168306. [Google Scholar] [CrossRef]
  15. Liverts, E.Z. Co-spherical electronic configuration of the helium-like atomic systems. Ann. Phys. 2022, 436, 168669. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liverts, E.Z.; Krivec, R. Fock Expansion for Two-Electron Atoms: High-Order Angular Coefficients. Atoms 2022, 10, 135. https://doi.org/10.3390/atoms10040135

AMA Style

Liverts EZ, Krivec R. Fock Expansion for Two-Electron Atoms: High-Order Angular Coefficients. Atoms. 2022; 10(4):135. https://doi.org/10.3390/atoms10040135

Chicago/Turabian Style

Liverts, Evgeny Z., and Rajmund Krivec. 2022. "Fock Expansion for Two-Electron Atoms: High-Order Angular Coefficients" Atoms 10, no. 4: 135. https://doi.org/10.3390/atoms10040135

APA Style

Liverts, E. Z., & Krivec, R. (2022). Fock Expansion for Two-Electron Atoms: High-Order Angular Coefficients. Atoms, 10(4), 135. https://doi.org/10.3390/atoms10040135

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop