_{1}

^{*}

We employ a recently amended Born-Oppenheimer (hereafter shortly BO) approximation
[1] to treat inelastic scattering of slow electrons from highly excited Rydberg atoms like e
^{-} + He(1
*s*
*n*
*s*)→He
^{-** }for
*n*
≫ 1. Along these lines we replace the standard BO set of potentials by an evolution operator. In this way we take a momentum-momentum coupling inadvertently disregarded by BO into account. The BO eigenvalue problem is now replaced by an evolution equation. One eigen-evolution has been identified as Wanner channel. That channel describes the diffraction of electron pairs from a potential ridge. That diffraction causes a phase jump of π/2 in the channel evolution. Moreover we present a new conservative attractive force controlling the motion of the electron pair as a whole in the nuclear field whose potential is given by
. The coupling constant
*g* has been calculated. That potential foreign to the standard BO approximation manifests itself by an entirely new series of isolated resonances located slightly below the double ionization threshold. This resonance ensemble compares favorably with experimental data. Further we present an evolution which forces the electron pair to the electrostatically unstable top of the potential ridge. That evolution may be regarded as quantum version of Wannier’s converging trajectory, and manifests itself here as Fresnel distribution.

In 1953 Wannier surprised the atomic community with an entirely unusual threshold ionization law for electron impact ionization of the hydrogen atom [

σ ∝ E 1.127 ⋯ . (1)

In the meantime there exists highly accurate experimental data which confirm that exponent (see [

Wannier’s picture is basically as follows. As long as the scattering electron is still away from the target atom the bound electron performs periodic Bohr orbits. As soon as the incident electron comes closer to the target these Bohr orbits experience a strong distribution due to correlation. At very low energy the electron-electron repulsion forces the three bodies into a collinear configuration (electron-nucleus-electron). Within that collinear configuration the electron pair experiences an unstable equilibrium position when both electrons have equal distances from the nucleus, say r 1 = r 2 . Actually a cut through the potential surface,

V = − Z ( 1 r 1 + 1 r 2 ) + 1 r 12 (2)

shows a repulsive barrier whose top is located at r 1 r 2 = 1 . In (2) r 12 is the electron-electron separation. Wannier has therefore introduced an angular coordinate

α = tan − 1 r 1 r 2 (3)

and studied the Langrage equations in the subspace (R, α) where

R = r 1 2 + r 2 2 (4)

is the so called hyperradius orthogonal to α. That subspace constitutes geometrically a ridge. Unfortunately, the coordinates R, α are neither in the classical Hamilton function nor in its quantum counterpart separable. Nevertheless, Wannier was able to show the existence of two special trajectories during the evolution along R. One trajectory converges to the top of the ridge whereas the other one diverges away from it. It is the aim of the present paper to present a mathematical frame within quantum theory which delivers two kinds of wave propagation which reduce in the classical limit to Wannier’s trajectories.

We do that in two steps. First we introduce an angle different from Wanniers α (see Section 2). Our angle (below denoted by φ) may be obtained by a rotation in the space ℝ 6 such that the better suited angle φ becomes a cyclic coordinate (see for instance [

Wannier [

R = ∑ i = 1 N r i 2 (5)

and complete that generalized coordinate with 3N − 1 hyperpolar angles corresponding to unit vectors on the hypersphere S 3 N − 1 . Each single electron position vector has then the form

r i = R R ^ i ( ω ) (6)

where ω stands for the set of angles, and R ^ is a unit vector depending on these angles. It is then straightforward to see that the Hamiltonian of any atom has in non-relativistic description the structure

H = T R + Λ 2 2 R 2 + C ( ω ) R . (7)

where Λ 2 is the Casimir operator of the rotation group SO_{6} and C ( ω ) is the potential multiplied by R.

In the case of two electrons Wannier pointed out that (7) at zero-energy obeys a similarity principle which describes classical trajectories in which the whole electron cloud performs at a set of constant angles hyperradial breathing motions. This requires of course that the electron cloud has an equilibrium configuration. Actually this is the case for few-electron atoms (see [

The case of two electrons in the field of one nucleus has been worked out in great detail. We follow here Klar et al. [

One of them we chose to describe the relative motion within the collinear configuration (electron-nucleus-electron). Following Sommerfeld [

φ = tan − 1 r 1 2 − r 2 2 2 r 1 r 2 (8)

( 0 ≤ φ ≤ 2 π ) whereas the bending motion may be described by the ratio of moments of inertia,

ψ = 1 2 tan − 1 ( Θ x x Θ y y ) . (9)

Macek [

The equilibrium configuration occurs in terms of our coordinates at ψ = 0 corresponding to the collinear configuration and at φ = π corresponding to equal electron-nucleus separations r 1 = r 2 . Since the Wannier phenomenon occurs only in the subspace ( R , φ ) we fix below the angle ψ = 0 .

Although the standard BO approximation [^{−} delivers indeed good results for singly and doubly excited states. The method fails, however, to describe the Wannier phenomenon. This is not surprising if the hyperradius R is selected as adiabatic coordinate. In the BO approximation the hyperradius R is treated as a constant in the channel expansion.

Actually, the standard BO approximation suffers from a serious mistake. It is usually believed that the validity rests on a small radial momentum in our present problem given by

p = 1 i d d R . (10)

That is not true, neither in the molecular nor in the atomic case. The small quantity for the definition of a collision channel is the radial kinetic energy given by

T R ∝ 1 2 P 2 = − 1 2 d 2 d R 2 (11)

We look now to the Hamiltonian applied to a product wave function, see (6),

H { F ( R ) Φ ( ω ; R ) } = 1 2 ( P 2 F ) Φ + ( P F ) ( P Φ ) + 1 2 F ( P 2 Φ ) + F h a n g Φ = − 1 2 d 2 F d R 2 Φ − d F d R ∂ Φ ∂ R − 1 2 F ∂ 2 Φ ∂ R 2 + F h a n g Φ (12)

where the angular Hamiltonian reads, see [

h a n g = Λ 2 2 R 2 + C ( ω ) R . (13)

BO introduces now mistakenly an adiabatic Hamiltonian disregarding the term PΦ. That corresponds to the assumption of constant values of R. It is obvious that this primitive approximation cannot deliver the Wannier mode since already in Wannier’s classical treatment R is a function of time. We need therefore an improvement which allows for variable values of R.

Recently [

{ ( P F ) P + F h a n g } Φ . (14)

For a finite radial function F ≠ 0 that leads to post-adiabatic eigen-channels defined by

{ − d F d R ∂ ∂ R + h a n g } Φ = W Φ (15)

Before we go more deeply into the new development we review briefly the application of the standard BO approximation to two-electron atoms. After splitting

off a factor of R 5 2 from the wave function the time-independent wave equation shows in hyperspherical coordinates the structure

( − ∂ 2 ∂ R 2 + 15 8 R 2 + Λ 2 2 R 2 + C ( ψ , φ ) R − E ) Ψ ( R , ψ , φ ) = 0 (16)

see [

Ψ ( R , ω ) = F ( R ) Φ ( ω ; R ) (17)

where the channel function Φ solves (15) at constant values of R, i.e.

( 15 8 R 2 + Λ 2 2 R 2 + C ( ψ , φ ) R ) Φ = U Φ . (18)

The separation parameter U constitutes in the case of helium for example an e-He^{+} (nl) potential. The channel function Φ represents standing waves on the hypersurface S 5 . The Wannier mode, however, we expect to be a travelling wave, and therefore is not member of the adiabatic channel family.

The present sketch derivation above suffers from an incorrect separation of variables.

The zero-momentum approach identical to the standard BO approximation delivers for two-electron atoms good results for the lowest channels. But it becomes shaky for increasing excitation, and breaks down near threshold of double escape. Directly at threshold the classical Wannier modes of motion show a finite radial velocity R ˙ ( t ) ≠ 0 , see [^{1}, and expresses the derivative by its logarithmic one, i.e.

F ′ = F ′ F F = d ( log F ) d R F = i P ( log F ) F (19)

The logarithmic derivative, moreover, we use to define collision channels, in our case

P ( log F ) = K ( R ) (20)

K being the hyperradial wavenumber. That trick leads us to the channel equation

{ ( P F ) P + F h } Φ = F { i K ∂ ∂ R + h ( ω ; R ) } Φ = F W ( R ) Φ (21)

which reduces for a finite radial function F ≠ 0 to

{ i K ( R ) ∂ ∂ R + h ( ω ; R ) } Φ = W ( R ) Φ (22)

The radial wavenumber is given by twice the radial kinetic energy, expressed by total energy minus potential energy,

K ( R ) = ± 2 ( E − W ( R ) ) (23)

Equation (22) constitutes therefore a nonlinear eigenvalue equation because the lhs of (22) depends on the eigenvalue W, too. The quantity in curly brackets of (22) resembles a time-dependent wave equation where time has been replaced by radius. Like the time-dependent wave equation or the heat equation or our (22) describes an evolution. Therefore (22) may be rewritten as evolution Equation. To this end we introduce an evolution operator E ( ω | R , 0 ) which describes the evolution from the united atom corresponding to R = 0 to a finite value R ≠ 0 . Let us denote the channel function at R = 0 with Φ ( ω ; 0 ) we put

Φ ( ω ; R ) = E ( ω | R , 0 ) Φ ( ω ; 0 ) (24)

with the initial condition

E ( ω | 0 , 0 ) = 1 (25)

The evolution operator itself satisfies then an equation identical to (22)

( i K d d R + 15 8 R 2 + Λ 2 2 R 2 + C ( ψ , φ ) R ) E ( ω | R , 0 ) = W ( R ) E ( ω | R , 0 ) . (26)

With help of the wavenumber K we are now able to introduce collision channels, for instance by the energy level of the hydrogen-like target, here He^{+}, with energy levels

E n = − Z 2 2 n 2 (27)

In the following, we intent to treat also nuclei with arbitrary charges Z.

The wavenumber K in (23) reads in the present case, see [

K = K n = ± − Z 2 2 n 2 − W ( R ) (28)

The wavenumber in the Wannier spectral range needs a modification. The zero-energy domain is there controlled by a long-range Coulomb force. Note, at threshold E = 0 the Coulomb zone extends to R = ∞ . The Equation (28) needs therefore a modification since we expect W(R) to be short-range potential. We proceed as follows.

Near threshold we expect the electron pair to be located near the ridge top. Therefore we approximate the ridge top by a Taylor expansion,

C ( ψ = 0 , φ ) = − C 0 − C 2 ( φ − π ) 2 + O ( 3 ) . (29)

The coefficients are given by

C 0 = 4 Z − 1 2 (30)

equal to the net charge of the correlated electron pair on the ridge top, and the ridge curvature

C 2 = 12 Z − 1 4 2 (31)

(29) identifies a long-range interaction representing a Coulomb potential in the space ℝ 6 . i.e. the correlated electron pair as a whole with a charge C_{0} is attracted by the nucleus. Thus we have to replace the eigenvalue W in (28) by – C 0 / R . Along these lines and using (28) the Wannier wavenumber becomes

lim n → ∞ K n = K ∞ = 2 C 0 R (32)

Finally, we need the operator Λ 2 . That reads in the frame of our model at fixed ψ = 0 simply

Λ 2 = − 4 ∂ 2 ∂ φ 2 (33)

since the azimuth angle φ is a cyclic coordinate, see [

{ i 2 C 0 R ∂ ∂ R + 15 8 R 2 − 2 R 2 ∂ 2 ∂ φ 2 − C 2 R ( φ − π ) 2 } E ( φ | R , 0 ) = W ( R ) E ( φ | R , 0 ) (34)

The following Section 4 solves that equation (34) and discusses the result.

We solve now the partial differential Equation (34) in the Coulomb zone employing the initial condition E ( φ | 0 , 0 ) = 1 . To this end we disregard the centrifugal term, and try the exponential Ansatz

E = exp { i γ g ( R ) ( φ − π ) 2 } (35)

with g ( 0 ) = 0 and a constant γ ∈ ℂ . Along these lines we conclude immediately g ( R ) = R provided the constant γ satisfies the quadratic equation

8 γ 2 − C 0 2 γ − C 2 = 0 (36)

Equation (36) may be interpreted as description of a potential surface curvature deformation.

The solution (35) delivers for the rhs of (34) 4 i γ R 3 / 2 E . We stress that the Hermiticity of the lhs of (34) excludes an imaginary eigenvalue. We conclude therefore that the diffraction of the correlated two-electron wave from the ridge has caused a phase jump^{2} of π 2 , and forces the rhs of (34) to the interpretation

W exp { i ( π 2 + γ R ( φ − π ) 2 ) } (37)

with the real eigenvalue

W ( R ) = 4 γ R 3 / 2 (38)

provided γ is real. Actually we get the from (35) two real solutions given by

γ 1 , 2 = 2 4 32 ( 4 Z − 1 ± 100 Z − 9 ) (39)

where we have used (30), (31).

The existence of two solutions is in agreement with Wannier’s observation. He got one trajectory converging to the equilibrium point φ = π and one diverging away from it. Our positive γ value describes indeed an evolution converging to the ridge top as may be seen as follows. The evolution (35) constitutes a Fresnel distribution with the property

lim R → ∞ exp { i γ R ( φ − π ) 2 } = i π γ R − 1 4 δ ( φ − π ) E ( φ = π ; R ≫ 1 ) . (40)

Here we have used the definition of a Fresnel Distribution given by

Δ Fresnel = a i π η exp { i a x 2 η } (41)

valid for x real and real a > 0, see [

The diverging solution ( γ < 0 ) leads Wannier’s study to double electron escape. In our quantum description the corresponding solution delivers an attractive potential operating in the ℝ 6 space controlling the motion of the electron pair as a whole. That potential is given by

W ( R ) = − g R 3 / 2 (42)

with the coupling constant

g = 2 4 8 ( 100 Z − 9 − 4 Z − 1 ) (43)

see (35). That result deserves several comments. First, a fractional power of R in (36) is entirely unexpected. At a first glance one would expect the Coulomb law − C 0 / R . But the diffraction pattern of the pair wave from the ridge leads to an ℝ 6 rotational symmetry breakdown and yields the surprising potential (42). The corresponding effective potential

U ( R ) = 15 8 R 2 − g R 3 / 2 (44)

predicts a series of resonances below the double escape threshold. Such levels have been observed a long time ago by Cvejanovic and Read [^{−}. The scattering of slow electrons from highly excited helium e − + He ( 1 s , n s ) → He − ∗ ∗ showed a series of isolated lines below the ionization threshold of helium. A pilot analysis of their observed level spacing shows a decrease slower than in a Rydberg series. We believe that ref [

Shell structure and circular orbits constitute since the early days of quantum mechanics the theoretical frame for multi-electron atoms although these terms stem from the exactly solvable hydrogen atom. Nevertheless models based on that theoretical frame have been surprisingly successful provided correlation is only a small effect. Near thresholds of multiple ionization, however, exist no longer circular orbits. They have been replaced by Wannier modes of motion.

Wannier [

The present paper presents the necessary quantum version, and opens an entirely new field of atomic physics full with surprises. We have shown that the new situation shortly referred to as Wannier phenomenon can be described with help of an amended BO approximation which replaces the set of BO potentials by an evolution operator. That operator delivers in fact an eigen-evolution which we refer to as Wannier channel. That channel delivers the threshold cross section (1), as well as a set of new isolated resonances slightly below threshold. They

manifest themselves as eigenvalues of the potential W ( R ) = − g R 3 / 2 , see (42).

Moreover, we show the Wannier evolution to be identical to a Fresnel distribution converging to a Dirac delta distribution centered at an equilibrium configuration of the electron pair. That corresponds to Wannier’s converging trajectory which describes an electron-electron attraction. That attraction emerges from the diffraction of a dominantly correlated two-electron wave from a potential ridge. In summary we find out that the present work has confirmed our earlier suspiciousness on that topic [

Finally, we remark that this Wannier phenomenon worked out here can be extended to three/four excited electrons [

The author declares no conflicts of interest regarding the publication of this paper.

Klar, H. (2020) Quantum Theory of the Wannier Phenomenon. Journal of Applied Mathematics and Physics, 8, 2416-2426. https://doi.org/10.4236/jamp.2020.811178