On Invariant Vectors in the Presence of Electric and Magnetic Fields

1. Introduction

The subject of invariants in physics, even restricted to the narrower field of atomic physics, is far too broad for any complete survey. The aim of the present article is accordingly rather more limited: attempting to review a selection of works, starting from the Kepler problem and then considering the cases of electric and magnetic fields.

The German mathematician Noether showed, in 1918, that any invariance in a continuous group of transformations, of spatio-temporal nature or not, with n parameters, implies the existence of n conservation laws [1,2]. As Curie’s principle (the symmetries of the causes are to be found in the effects) and Wigner’s theorem (specifying how physical symmetries such as rotations, translations, as well as charge, parity and time reversal symmetries manifest themselves on the Hilbert space of states), Noether’s theorem allows one to exploit the knowledge of the continuous symmetries of a system. Noether showed that, conversely, a conservation law is never accidental but derives from the existence of a continuous symmetry. It is worth mentioning that the importance of Noether’s contribution to science was only recognized around 1960.

In the case of mechanics problems, “dynamical” invariants are revealed, which cannot be associated with any continuous space-time symmetry and are therefore linked to the particular analytical properties of the interactions responsible for the motion. The first example of dynamical symmetry was discovered in 1710 by Hermann, who demonstrated [3,4] with the help of Bernoulli [5], that the motion of a mass in a Keplerian potential is characterized by two invariants: the angular momentum, which determines the plane of motion, and a second vector - called today the Laplace-Runge-Lenz vector - which fixes the orientation of the major axis of the elliptical orbit in its plane, and whose modulus is proportional to the eccentricity of the orbit. Laplace exploited the existence of this vector [6], which was rediscovered by Hamilton [7], described by Runge in a didactic work [8], and finally used by Lenz [9] and Pauli [10] to find the quantum properties of the hydrogen atom [11], characterized by a specific degeneracy of the energy levels called Coulomb degeneracy. In 1873, Bertrand showed that the orbit of a particle placed in central potential does not close on itself except if the potential is Keplerian or harmonic [12]: in the first case, there is a unique privileged direction in the plane of the orbit, which is precisely that of the Laplace vector; in the second case, there are two privileged directions, which are perpendicular and constitute the proper directions of the invariant tensor. According to Noether’s theorem, a dynamical symmetry is associated with the Laplace vector: it is that of the $SO\left(4\right)$ group (with six parameters) of rotations in four-dimensional space, obtained by Fock [13] and Bargmann [14]. The symmetry of the isotropic harmonic oscillator in three dimensions was discovered in 1940 by Jauch and Hill [15] and is that of the $SU\left(3\right)$ group (with eight parameters) of unitary matrices $3\times 3$ of determinant + 1. More than twenty years later, Redmond found an invariant in the presence of an electric field [16].

There are two main versions of the so-called two-center problem. The first one is the motion of an electron in the field of two stationary Coulomb centers of charges Z and $Z^{\prime }$ separated by a distance R, which is one of the most fundamental problems in quantum mechanics. The second one is the motion of a planet in the gravitational field of two stationary stars of generally different masses, which is one of the most important problems in celestial mechanics. Kryukov and Oks [17,18] derived a supergeneralized Runge-Lenz vector whose projection on the internuclear or interstellar axis does conserve.

A dyon is a hypothetical particle with both electric and magnetic charges. When its electric charge is zero, the particle is usually referred to as a magnetic monopole. Many grand unified theories predict the existence of both magnetic monopoles and dyons. Dyons were first proposed as a phenomenological alternative to quarks [19]. Schwinger extended the Dirac quantization condition [20] to the dyon and predicted the existence of a particle with the properties of the $J/\psi$ meson prior to its discovery in 1974. There is much interest in the problem of bound states of a fermion in the field of a fixed Dirac monopole.

In Section 2, we introduce the Runge-Lenz vector in classical mechanics and discuss the Bertrand theorem. In Section 3, in the framework of the $SO\left(4\right)$ rotation group, the quantum-mechanical analog to the Lenz vector is obtained and applied to the determination of the energy levels of the hydrogen atom and their degeneracy. In Section 4, the analogy between the Kepler and harmonic potentials is discussed and the Laplace invariant tensor is introduced. The Redmond invariant in the presence of an electric field is introduced in Section 5 and the Kryukov-Oks supergeneralized Runge-Lenz vector for the two-center problem is presented in Section 6. The Landau-Avron-Sivardière approach of the case of a magnetic field is discussed in Section 7, as well as the treatment of the hydrogen atom in crossed electric and magnetic fields. Finally, the dyon-dyon system is briefly evoked in Section 8 in the framework of the MIC (McIntosh-Cisneros)-Kepler problem. The linear Stark effect in the MIC-Kepler problem describing the interaction of charged particle with Dirac’s dyon is considered.

2. The Kepler problem: the Runge-Lenz vector and Bertrand’s theorem

Let us consider a mass m subject to the (attractive) Kepler potential of a mass M. The reduced mass being $\mu =mM/(m+M)$, the potential reads

$$U\left(r\right)=-\frac{\alpha }{r}=-G\frac{Mm}{r}=-G\frac{\mu (m+M)}{r}=-\kappa \frac{\mu }{r}.$$

(1)

with $\alpha =\kappa \mu >0$ and G the Cavendish gravitational constant. One thus has $\kappa =G(m+M)$. The potential has spherical symmetry and is independent of time. Therefore, kinetic energy and angular momentum are invariant with time. We also know that, whatever the initial conditions, the orbit is an ellipse for the bound states (negative energy), a parabola for the zero-energy states and a hyperbola for the positive-energy states.

2.1. The Runge-Lenz, eccentricity and Hamilton vectors

As mentioned in the introduction, the so-called Runge-Lenz vector was discovered by Hermann, transmitted by Bernoulli, used intensively by Laplace, exhumed by Runge at the beginning of the 20^th century and used by Lenz and Pauli in quantum mechanics in order to study the hydrogen atom. Although the vector should therefore be called Hermann-Bernoulli-Laplace-Runge-Lenz vector [4,5,21], we will in the following, for simplicity and because it is the most common way it is referred to in the literature, retain the denomination “Runge-Lenz” vector. The Runge-Lenz vector reads [22,23]:

$$\mathbf{A}=\mathbf{v}\wedge \mathbf{L}-\alpha \frac{\mathbf{r}}{r}=\mathbf{v}\wedge \mathbf{L}-\alpha {\mathbf{u}}_r,$$

(2)

where $\mathbf{v}$ is the velocity of our system (mass m), $\mathbf{L}=\mu \mathbf{r}\wedge \mathbf{v}$ the angular momentum, and ${\mathbf{u}}_r=\mathbf{r}/r$ the radial unit vector. Throughout the paper, the symbol ∧ denotes the vector product:

$$\left(\begin{array}{c}{a}_1\\ a_2\\ a_3\end{array}\right)\wedge \left(\begin{array}{c}{b}_1\\ b_2\\ b_3\end{array}\right)=\left(\begin{array}{c}{a}_2{b}_3-a_3{b}_2\\ a_3{b}_1-a_1{b}_3\\ a_1{b}_2-a_2{b}_1\end{array}\right).$$

(3)

The equation of motion is

$$\frac{d\mathbf{v}}{dt}=-\frac{\alpha }{\mu r^2}\frac{\mathbf{r}}{r}=-\frac{\alpha }{\mu r^2}{\mathbf{u}}_r,$$

(4)

and one has

$$\frac{d}{dt}\frac{\mathbf{r}}{r}=\frac{1}{L}\mathbf{L}\wedge {\mathbf{u}}_r\frac{d\theta }{dt}=\frac{\mathbf{L}\wedge {\mathbf{u}}_r}{\mu r^2}.$$

(5)

It is also possible to obtain it more directly. Using the invariance of angular momentum

$$L=\mu r^2\frac{d\theta }{dt}$$

(6)

which gives

$$\frac{d\mathbf{v}}{dt}+\frac{\alpha }{L}\frac{d\theta }{dt}{\mathbf{u}}_r=0$$

(7)

and taking the vector product of the latter equality by $\mathbf{L}$, one has

$$\mathbf{L}\wedge \frac{\mathbf{v}}{dt}+\alpha \frac{d{\mathbf{u}}_r}{dt}=0,$$

(8)

giving directly the invariance of the Runge-Lenz vector. The projection of the Runge-Lenz vector on the big axis of the ellipse gives the equation of the conic. Calculating the scalar product $\mathbf{A}.\mathbf{r}$, one gets

$$Arcos\theta =\mu \left[\mathbf{v}\wedge (\mathbf{r}\wedge \mathbf{v})\right]-\alpha r=\mu \left[r^2{v}^2-{(\mathbf{r}.\mathbf{v})}^2\right]-\alpha r=\mu C^2-\alpha r,$$

(9)

where we have used

$$\mathbf{v}\wedge \left(\mathbf{r}\wedge \mathbf{v}\right)=\left(\mathbf{v}·\mathbf{v}\right)\mathbf{r}-\left(\mathbf{r}·\mathbf{v}\right)\mathbf{v}.$$

(10)

The polar equation of the orbit is then [24,25]:

$$r=\frac{p}{1+ecos\theta },$$

(11)

with $p=C^2/k=L^2/\left(\kappa {\mu }^2\right)$, $C=L/\mu$ being the areal constant (twice the areal velocity, Kepler’s second law) and $e=A/\left(k\mu \right)=A/\alpha$ the eccentricity of the ellipse. The energy of the motion can be obtained by calculating the norm of the Runge-Lenz vector

$$E=\frac{\mu }{2}\frac{A^2-{\alpha }^2}{L^2}.$$

(12)

Using the definitions of p and e, as well as the relation $p=a(1-e^2)$, we deduce

$$E=-\frac{\kappa \mu }{2a}=-\frac{\alpha }{2a}$$

(13)

without using the Binet formulas [26] (see Section 2.2). There are other invariant vectors, which are related to the Runge-Lenz vector.

In celestial mechanics, the eccentricity vector of a Kepler orbit [27,28] is the dimensionless vector with direction pointing from apoapsis to periapsis (respectively farthest and nearest points in the orbit of a planetary body about its primary body to the center of force) and with magnitude equal to the orbit’s scalar eccentricity. For Kepler orbits, the eccentricity vector is a constant of motion. It is of great interest in the analysis of almost circular orbits, as perturbing (non-Keplerian) forces on an actual orbit will cause the osculating eccentricity vector to change continuously as opposed to the eccentricity and argument of periapsis parameters for which eccentricity zero (circular orbit) corresponds to a singularity. The eccentricity vector is

$$\mathbf{e}=\frac{\mathbf{v}\wedge \mathbf{L}}{\alpha }-\frac{\mathbf{r}}{\left|\mathbf{r}\right|},$$

(14)

which follows immediately from the vector identity (10). One has

$$\mathbf{e}=\frac{\mathbf{A}}{\mu \kappa }=\frac{1}{\kappa }(\mathbf{p}\wedge \mathbf{L})-{\mathbf{u}}_r,$$

(15)

with $\mathbf{p}=\mu \mathbf{v}$. The Runge-Lenz vector $\mathbf{A}$ is directed from the centre of force towards the periapsis, making it easy to visualize (Lenz and Pauli called it the “Achsenvektor”) [29]. However, its physical dimension (SI unit kg m${ }^3$ s${ }^{-2}$) does not endow it with any familiar meaning; the alternative invariant $\mathbf{A}/m$ is hardly better. If $\mathbf{e}$ is used instead of $\mathbf{A}$ in the derivation of the orbit equation, a simple equation appears, with an immediate physical meaning. This is why it is often considered that e should be preferred over A for didactic applications. In the literature, the only “eccentricity-like” vector referred to seems to be the one defined by Hamilton, $\mathbf{L}\wedge \mathbf{A}/L^2$: it is co-linear to $\mathbf{e}$, but has the dimension of a velocity or of a momentum according to the alternative definition of [30]. It does not coincide with a particular value of v or $\mu v$, and thus has no special significance. More importantly, while equations of motion have been obtained with all vectors, each one has its own mathematical simplicity (complexity) and physical meaning. Finally, the eccentricity vector naturally arises in solving the dynamical equation at an elementary level, without resorting to Lagrangian formalism or group theory. The eccentricity vector also only arises in the case of an inverse-square law of force.

The vector product of $\mathbf{L}$ and $\mathbf{A}$ is an invariant. It is thus possible to deduce the invariant Hamilton $\mathbf{S}$ vector

$$\mathbf{S}=\frac{1}{L^2}\mathbf{L}\wedge \mathbf{A}=\mathbf{v}-\frac{\alpha }{L^2}\mathbf{L}\wedge {\mathbf{u}}_r.$$

(16)

The Laplace vector indicates the direction of the periapsis. In the same way, the Hamilton vector is parallel to the velocity of the system at the periapsis [31,32,33,34].

2.2. The Bertrand theorem

The Bertrand theorem [21,35,36,37,38,39,40,41,42,43,44,45] states that the orbit is closed only if the potential is in $1/r$ or harmonic $r^2$. In the Keplerian case, there are two privileged directions: the one relating the attractive center to the periapsis of the orbit, and the perpendicular direction. The two directions play different roles. This justifies the existence of the vector $\mathbf{A}$. As mentioned above, the vector $\mathbf{A}$ is parallel to the big axis of the ellipse and indicating the direction of the periapsis. In the case of an harmonic ($r^2$) potential, there are also two privileged directions. The Binet formulas are, with $u=1/r$:

• First Binet formula:

  $$v^2=C^2\left[u^2+{\left(\frac{du}{d\theta }\right)}^2\right],$$

  (17)
• Second Binet formula:

  $$a_r=-C^2{u}^2\left[u+\frac{d^{2}u}{d{\theta }^2}\right],$$

  (18)

where $a_r$ is the acceleration. Let us consider a force proportional to $r^n$ [46]:

$$a_r=-\kappa r^n$$

(19)

and a small perturbation $ϵ$ to the circular trajectory of radius $r_0$ (corresponding to $d^{2}u/d{\theta }^2=0$):

$$C^2/r_0=\kappa r_0^{n+2}.$$

(20)

One has

$$\frac{1}{r}=\frac{1}{r_0}(1+ϵ),$$

(21)

where $ϵ\ll 1$. Using the second Binet formula, one obtains

$$\frac{d^2ϵ}{d{\theta }^2}=-(n+3)(ϵ+\beta {ϵ}^2+\gamma {ϵ}^3)$$

(22)

where we have set $2\beta =-(n+2)$ and $6\gamma =(n+2)(n+1)$. We have to solve an equation describing a non-linear oscillation of variable $ϵ$. For that purpose, we can resort to a perturbative expansion of the solution and evaluate the apsidal angle of the orbits in the neighborhood of the circular one. At first order of approximation, in the vicinity of the circular trajectory, one has

$$\frac{d^2{ϵ}_1}{d{\theta }^2}=-(n+3){ϵ}_1$$

(23)

and thus

$${ϵ}_1=Mcos\left(p\theta \right),$$

(24)

with $p^2=n+3$. At the next (second) order, we get

$$\frac{d^2{ϵ}_2}{d{\theta }^2}=-(n+3)\left\{{ϵ}_2-\beta \frac{M^2}{2}\left[1+cos\left(2p\theta \right)\right]\right\}$$

(25)

yielding

$${ϵ}_2=Mcos\left(p\theta \right)+M^2(G+Acos\left(2p\theta \right)),$$

(26)

with $G=(n+2)/4$ and $A=-(n+2)/12$. Finally, at the third order, we look for solutions of the type

$${ϵ}_3=Mcos\left(q\theta \right)+M^2\left[G+Acos\left(2p\theta \right)\right]+M^{3}Bcos\left(3p\theta \right)$$

(27)

satisfying the equation if and only if

$$q^2=(-n+3)\left[1-\frac{M^2}{12}(n-1)(n+2)\right]$$

(28)

where we have set

$$G=(n+2)/4$$

(29)

as well as

$$B=(n+2)(n+3)/96.$$

(30)

The solution at third order, around the circular trajectory, is of the form $ϵ=Mcos\left(q\theta \right)+O\left(M^2\right)$, where

$$q=\sqrt{n+3}\left[1-\frac{M^2}{24}(n+2)(n-1)\right]$$

(31)

and the apsidal angle, equal to $\pi /q$ can not be commensurable with $\pi$$\forall M\ll 1$, except if $n=1$ or $n=-2$. Therefore, the orbits of a point subject to a central force proportional to the inverse of the power n of the distance between the active center and the point are closed if and only if $n=1$ or $n=-2$.

3. The Runge-Lenz vector in quantum mechanics

3.1. Construction of the $SO\left(4\right)$ rotation group - Quantum mechanical analog to the Runge-Lenz vector

The resemblance between the Kepler problem of planetary motion and the charged particle in Coulombic potential motivates us to understand symmetries possessed by such systems. The Hamiltonian can be written formally as

$$H=\frac{{\mathbf{p}}^2}{2\mu }-\frac{\alpha }{r},$$

(32)

where $\mu$ is the reduced mass $Mm/(M+m)$, $\alpha =GMm=G\mu (M+m)=\kappa \mu$ (for Newton potential, see above) and $\alpha =Z{e}^2$ (for Coulomb potential with $4\pi {ϵ}_0$ set to unity). With $\mathbf{r}$, $\mathbf{p}$, and $\mathbf{L}$ as already established quantum mechanical operators, the classical Lenz vector $\mathbf{A}$ can be translated directly into quantum mechanics (vector $\mathbf{M}$). There is only a bit of subtlety involved in its construction. We must take care in defining the cross product (see Equation (3). We notice that $\mathbf{p}\wedge \mathbf{L}\ne -\mathbf{L}\wedge \mathbf{p}$. The classical Runge-Lenz operator is not hermitian. Following the derivation of Ref. [47], we redefine $\mathbf{M}$ as a symmetric average. The quantum-mechanical operator $\mathbf{M}$, extension of the classical Runge-Lenz vector, defined as follows

$$\mathbf{M}=\frac{1}{2\mu }\left(\mathbf{p}\wedge \mathbf{L}-\mathbf{L}\wedge \mathbf{p}\right)-\alpha {\mathbf{u}}_r.$$

(33)

As expected, this vector commutes with the Hamiltonian, $\left[\mathbf{M},H\right]=0$ and defines a conserved quantity in the quantum mechanical analog of the classical Kepler problem. We find that

$$\mathbf{L}.\mathbf{M}=\mathbf{M}.\mathbf{L}=0,$$

(34)

and

$${\mathbf{M}}^2=\frac{2H}{\mu }({\mathbf{L}}^2+{\hslash }^2)+{\alpha }^2=0.$$

(35)

From here, we look to understand how this new commuting operator relates to the physical system and all the existing conserved quantities. We examine the system in light of well-established geometrical symmetries.

3.2. The algebra of the $\mathbf{M}$, $\mathbf{L}$ generators

The three components of $\mathbf{M}$, like the three components of $\mathbf{L}$, can be treated as generators of some infinitesimal transformations. With the goal of exploring the algebraic structure of these new generators, we work out their commutation relations (there are thirty-six in all). Some are known; indeed, $\left[L_i,L_j\right]$ are given by

$$\left[L_i,L_j\right]=i\hslash {ϵ}_{ijk}{L}_k$$

(36)

for $i,j,k=1,2,3$, where ${ϵ}_{ijk}$ is the Levi-Civita symbol

$${ϵ}_{ijk}=\left|\begin{array}{ccc}{\delta }_{i1}& {\delta }_{i2}& {\delta }_{i3}\\ {\delta }_{j1}& {\delta }_{j2}& {\delta }_{j3}\\ {\delta }_{k1}& {\delta }_{k2}& {\delta }_{k3}\end{array}\right|$$

(37)

and ${\delta }_{ij}$ the Kronecker symbol. The others are

$$\begin{array}{ccc}\hfill \left[M_i,L_i\right]& =& 0\hfill \\ \hfill \left[M_i,L_j\right]& =& i\hslash {ϵ}_{ijk}{M}_k,\hfill \end{array}$$

(38)

and

$$\left[M_i,M_j\right]=-\frac{2i\hslash }{\mu }{ϵ}_{ijk}H{L}_k.$$

(39)

The first set of these commutators (38) shows that $\mathbf{M}$ is a vector. $M_i$ commutes with $L_i$ since $L_i$ generates rotation about the i axis, which has no effect on a vector pointed along that axis. On the other hand, if j is perpendicular to i, $L_j$ induces a change in the direction of $M_i$. This is characterized by the second commutation relation in Equation (38), which implies that the direction of this change in $M_j$ should be orthogonal to both $M_j$ and $L_i$, in this case, along $-M_k$. The relation in Equation (39) breaks the closed algebra of the L and M operators together, since H is involved. We now make a useful substitution. Considering a degenerate subspace (manifold) of H with energy eigenvalue E, we can replace the Hamiltonian with the resulting energy eigenvalue, and define ${\mathbf{M}}^{\prime }$ such that

$${\mathbf{M}}^{\prime }={\left(-\frac{\mu }{2E}\right)}^{1/2}\mathbf{M}.$$

(40)

The commutation relations remain unchanged with M replaced everywhere by $M^{\prime }$, with the exception of Equation (39), which now becomes

$$\left[M_i^{\prime },M_j^{\prime }\right]=i{ϵ}_{ijk}{L}_k,$$

(41)

signature of a closed algebraic system.

3.3. Identifying the closed algebra of $\mathbf{L}$ and ${\mathbf{M}}^{\prime }$ with the $SO\left(4\right)$ group

We can relabel $L_x$, $L_y$, and $L_z$ into $L_1$, $L_2$, and $L_3$. Substituting 1, 2 and 3 for x, y and z allows for a natural extension beyond three dimensions. The fact that the commutators $\left[r_i,p_j\right]$ are zero only when $i\ne j$ implies that the ${\mathbf{M}}^{\prime }$ matrices are merely generalized angular momentum operators in the fourth dimension. In four dimensions, we have $\left(\begin{array}{c}4\\ 2\end{array}\right)=6$ orthogonal angular momentum generators, where $\left(\begin{array}{c}n\\ p\end{array}\right)=n!/p\S (n-p)!$ is the usual binomial coefficient. Let us make the following associations:

$$\begin{array}{ccc}\hfill M_x^{\prime }& =& L_{14},M_y^{\prime }=L_{24},M_z^{\prime }=L_{34},\hfill \\ \hfill L_x& =& L_{23},L_y=L_{31},L_z=L_{12},\hfill \end{array}$$

(42)

which can be summarized as

$$M_i^{\prime }=L_{i4},and{L}_i={ϵ}_{ijk}{L}_{jk}$$

(43)

We add a fictitious coordinate $\omega$, so that the rotation in four-dimensional space is given by

$$\left(\begin{array}{c}{x}^{\prime }\\ y^{\prime }\\ z^{\prime }\\ {\omega }^{\prime }\end{array}\right)=exp\left[\sum _{i<j}{\theta }_{ij}{L}_{ij}\right]\left(\begin{array}{c}x\\ y\\ z\\ \omega \end{array}\right).$$

In order to obtain the energy, we can exploit the $SO\left(4\right)$ symmetry. It is easy to verify that the commutation relations in Section 4.3 still hold. In four dimensions, a degree of decoupling emerges: the maximal commuting set of operators includes two angular momentum operators at the same time: $L_{12}$ with $L_{34}$, etc. Such a decoupling stems from the property that the two operators generate rotations in orthogonal, non-intersecting planes, defined by completely disjoint sets of vectors. This fact alludes to the structure of the group generated by the operators: the rotational group in four dimensions, $SO\left(4\right)$, of rank 2, which contains, as a subgroup, the rank 1 $SO\left(3\right)$ rotational group. We have therefore generalized the closed algebra of $\mathbf{L}$ and $\mathbf{M}$ operators. The remarkable insight here is that a dynamical symmetry (as physically evidenced by eccentricities in closed orbits) is a mere artifact of rotational symmetry in a higher dimension.

3.4. Energy levels of hydrogen atom

The partial decoupling of angular momentum operators in four dimensions hints at a decomposition of the $SO\left(4\right)$ group generated by the operators L and M[48,49]. Rank 2 $SO\left(4\right)$ groups by rule can be broken down into two completely decoupled rank 1 $SU\left(2\right)$ groups:

$$SO\left(4\right)\to SU\left(2\right)\times SU\left(2\right),$$

(44)

with three non-commuting generators each, yielding $3+3=6$ generators. We make then the following choice of basis

$$\mathbf{I}=\frac{1}{2}(\mathbf{L}+{\mathbf{M}}^{\prime })$$

(45)

and

$$\mathbf{K}=\frac{1}{2}(\mathbf{L}-{\mathbf{M}}^{\prime })$$

(46)

with

$$\begin{array}{ccc}\hfill \left[I_i,I_j\right]& =& i\hslash {ϵ}_{ijk}{I}_k,\hfill \end{array}$$

$$\begin{array}{ccc}\hfill \left[K_i,K_j\right]& =& i\hslash {ϵ}_{ijk}{K}_k,\hfill \\ \hfill \left[\mathbf{I},\mathbf{K}\right]& =& 0,\hfill \end{array}$$

(47)

$$\begin{array}{ccc}\hfill \left[\mathbf{I},H\right]& =& \left[\mathbf{K},H\right]=0.\hfill \end{array}$$

(48)

Since both $\mathbf{I}$ and $\mathbf{K}$ constitute an $SU\left(2\right)$ known algebra, it is isomorphic to the group generated by angular momentum in three dimensions. We use our findings in Section 4.2 to get their possible eigenvalues. We can write states which are simultaneous eigenstates of H, $I_z$ and $K_z$:

$$|E;i,i_z,k,k_z\rangle .$$

(49)

and one has

$$I_z|E;i,i_z,k,k_z\rangle =i_z\hslash |E;i,i_z,k,k_z\rangle$$

(50)

for $i,k=0,1/2,1,\cdots$. The corresponding Casimir operators, ${\mathbf{I}}^2$ and ${\mathbf{K}}^2$, are respectively

$${\mathbf{I}}^2=\frac{1}{4}{\left(\mathbf{L}+{\mathbf{M}}^{\prime }\right)}^2$$

(51)

and

$${\mathbf{K}}^2=\frac{1}{4}{\left(\mathbf{L}-{\mathbf{M}}^{\prime }\right)}^2.$$

(52)

It is important to remember that Casimir operators are distinct from generators. They are actually functions of generators and although they commute with all the generators of a group, are not considered to be in the set of commuting generators itself. A linear combination of these operators is diagonal in the basis of their common eigenvectors. The Casimir operators ${\mathbf{I}}^2$ and ${\mathbf{K}}^2$ commute with all the generators $(I_x,I_y,I_z,H,K_x,K_y,K_z)$. One has

$${\mathbf{I}}^2|E;i,i_z,k,k_z\rangle =i(i+1){\hslash }^2|E;i,i_z,k,k_z\rangle$$

(53)

and

$${\mathbf{K}}^2|E;i,i_z,k,k_z\rangle =k(k+1){\hslash }^2|E;i,i_z,k,k_z\rangle .$$

(54)

From $\mathbf{I}$ and $\mathbf{K}$ we can construct operators

$$C_1={\mathbf{I}}^2-{\mathbf{K}}^2$$

(55)

and

$$C_2={\mathbf{I}}^2+{\mathbf{K}}^2,$$

(56)

which commute with all $so\left(4\right)$ generators and H. $C_1$ yields a constraint on the eigenvalues through

$$\begin{array}{ccc}\hfill C_1& =& {\mathbf{I}}^2-{\mathbf{K}}^2\hfill \\ & =& \frac{1}{4}({\mathbf{L}}^2+{{\mathbf{M}}^{\prime }}^2+\mathbf{L}.{\mathbf{M}}^{\prime }+{\mathbf{M}}^{\prime }.\mathbf{L})\hfill \\ & =& \mathbf{L}.{\mathbf{M}}^{\prime }=0.\hfill \end{array}$$

(57)

The last line used the relation from Equations (34) and (35). Actually, $C_1=0$ means that i and k are the same. This brings additional information (not given by $SU\left(2\right)$). We find the eigenvalues for operator $C_2$ with ease using Equations (53) and (54), yielding

$$C_2=2{\mathbf{K}}^2$$

(58)

and

$$C_2|E;i,i_z,k,k_z\rangle =2k(k+1){\hslash }^2|E;i,i_z,k,k_z\rangle .$$

(59)

Combining Equations (34), (35), (40), (51) and (52), we derive a constant expression for $C_2$:

$$\begin{array}{ccc}\hfill C_2& =& {\mathbf{I}}^2+{\mathbf{K}}^2=\frac{1}{2}\left({\mathbf{L}}^2+{{\mathbf{M}}^{\prime }}^2\right)=\frac{1}{2}\left({\mathbf{L}}^2-\frac{\mu }{2E}{{\mathbf{M}}^{\prime }}^2\right)\hfill \\ & =& \frac{1}{2}\left[{\mathbf{L}}^2-\frac{\mu }{2E}\left(\frac{2E}{\mu }({\mathbf{L}}^2+{\hslash }^2)+{\alpha }^2\right)\right]\hfill \\ & =& -\frac{{\hslash }^2}{2}-\frac{\mu {\alpha }^2}{4E}.\hfill \end{array}$$

(60)

Equating Equations (58) and (59) we obtain our permissible energy levels for a hydrogenic atom:

$$E_k=-\mu \frac{{\alpha }^2}{2{\hslash }^2{(2k+1)}^2}$$

(61)

for $k=0,1/2,1,\cdots$. Recalling that $\alpha =Z{e}^2$ and setting $2k+1=n$, for $n=1,2,3,\cdots$ we recover the well-known expression

$$E_n=-\mu \frac{{\alpha }^2}{2{\hslash }^2{n}^2}$$

(62)

for $n=1,2,3,\cdots$. The elegance of this method stems from the fact that, in addition to the hydrogen energy spectrum, the relevant constraints on various quantum numbers have emerged naturally as well. For instance, consider $\mathbf{L}=\mathbf{I}+\mathbf{K}$ from the initial definition of $\mathbf{I}$ and $\mathbf{K}$. Since $\mathbf{I}=\mathbf{K}$ , we see that by the triangular rule of vector addition, L can be at most $k+k=n-1$, and its values extend to $|k-k|=0$ in integer steps, as expected. In addition, following the algebra of $SU\left(2\right)$ groups, $I_z$ and $K_z$, analogous to $L_z$ would each have $2k+1=n$ independent eigenvalues and eigenstates. The total degeneracy of an energy level, $n\times n=n^2$, is therefore also recovered [50].

4. Analogy with the harmonic potential: invariant tensor

Let us summarize the analogies between the Kepler and the harmonic oscillator problems [51]:

• The orbit is closed (Bertrand’s theorem).
• The motion is periodic (degeneracy).
• There is a Laplace invariant.
• The motion has a dynamical symmetry higher than the rotation geometrical symmetry.
• For each value of the energy, one can find an infinity of classical orbits of different shapes.
• The separation of variables applies in several systems of coordinates.
• The quantum states are degenerate.

4.1. The Laplace tensor

We have seen that the Runge-Lenz vector is an invariant of the Kepler motion. According the the reciprocal of the Noether theorem [1], a continuous symmetry can be associated to it. Such a symmetry was discovered by Fock and Bargmann in 1935: the Coulomb potential has an internal symmetry (or dynamical symmetry), isomorphic (for the bound states) to the one of the continuous group $SO\left(5\right)$ of the rotations in the four-dimensional space. Such a group has six parameters and there are indeed six constants of motion: the three components of the angular momentum and the three components of the Runge-Lenz vector. Actually, only five of these constants are independent, since the vectors $\mathbf{L}$ and $\mathbf{A}$ are orthogonal. The energy is not a constant of motion independent of the previous ones, since it can be calculated as a function of the invariants of $A^2$ and $L^2$, which are Casimir invariants of the $SO\left(4\right)$ group.

The harmonic potential also has a dynamical symmetry, isomorphic to the one of the continuous group $SU\left(3\right)$ of $3\times 3$ matrices with determinant equal to one. This group has eight parameters. The eight corresponding constants of motion are the three components of the angular momentum and the five components of the symmetric Laplace tensor. Five of them are independent since $\mathbf{L}$ is an eigenvector of the Laplace tensor. Let us consider a mass $\mu$ in an isotropic harmonic potential

$$\mathcal{U}\left(r\right)=\frac{1}{2}k{r}^2.$$

(63)

The equation of motion is

$$\mu \frac{d\mathbf{v}}{dt}=-k\mathbf{r}.$$

(64)

It is possible to define an invariant of the motion [52] as

$$\overline{\overline{\mathrm{A}}}=\frac{k}{2}\mathbf{r}.\mathbf{r}+\frac{\mu }{2}\mathbf{v}.\mathbf{v}.$$

(65)

Let us set ${\omega }^2=k/\mu$. The Laplace tensor reads [53,54,55]:

$$\overline{\overline{\mathrm{A}}}=\left(\begin{array}{cc}{\omega }^2{x}^2+v_x^2& {\omega }^{2}xy+v_x{v}_y\\ {\omega }^{2}xy+v_x{v}_y& {\omega }^2{y}^2+v_y^2\end{array}\right),$$

(66)

where x and y are the coordinates of position operator $\mathbf{r}$ and $v_x$ and $v_y$ the coordinates of velocity $\mathbf{v}$. One has

$$\mathrm{Tr}\left(\overline{\overline{\mathrm{A}}}\right)={\lambda }_1+{\lambda }_2=E$$

(67)

as well as

$$\mathrm{Det}\left(\overline{\overline{\mathrm{A}}}\right)={\lambda }_1{\lambda }_2=\frac{1}{4}{\omega }^2{L}^2$$

(68)

yielding

$${\lambda }_1=\frac{1}{2}\left(E+\sqrt{E^2-{\omega }^2{L}^2}\right)$$

(69)

and

$${\lambda }_2=\frac{1}{2}\left(E-\sqrt{E^2-{\omega }^2{L}^2}\right)$$

(70)

and finally

$$E^2={\left|\overline{\overline{\mathrm{A}}}\right|}^2+\frac{k}{2m}{L}^2.$$

(71)

We have

$$\begin{array}{ccc}\hfill \mathbf{r}\overline{\overline{\mathrm{A}}}\mathbf{r}& =& \left(\begin{array}{cc}rcos\theta & rsin\theta \end{array}\right)\left(\begin{array}{cc}{\lambda }_1& 0\\ 0& {\lambda }_2\end{array}\right)\left(\begin{array}{c}rcos\theta \\ rsin\theta \end{array}\right)\hfill \\ & =& {\lambda }_1{r}^2{cos}^2\theta +{\lambda }_2{r}^2{sin}^2\theta \hfill \end{array}$$

(72)

and

$$\begin{array}{ccc}\hfill \mathbf{r}\overline{\overline{\mathrm{A}}}\mathbf{r}& =& \frac{k}{2}{r}^4+\frac{\mu }{2}{\left(\mathbf{r}.\mathbf{v}\right)}^2\hfill \\ & =& \frac{k}{2}{r}^4+\frac{\mu }{2}\left(r^2{v}^2-C^2\right).\hfill \end{array}$$

(73)

Eliminating v using the energy conservation gives

$$\frac{1}{r^2}=\frac{2{\lambda }_2}{\mu C^2}{cos}^2\theta +\frac{2{\lambda }_1}{\mu C^2}{sin}^2\theta .$$

(74)

If both eigenvalues are positive, the orbit is an ellipse with the semi-minor and semi-major axis given by

$$a^2=\frac{\mu C^2}{2{\lambda }_2},b^2=\frac{\mu C^2}{2{\lambda }_1}.$$

(75)

It is worth mentioning that the transformations of the group $SO\left(4\right)$ or $SU\left(3\right)$ enable one to change from an orbit of energy E to degenerated orbits of the same energy [56,57].

Analogies between the Kepler and harmonic motions were already known by Newton [58]. A mathematical transform from the one to the other has been introduced in 1965 by Kustaanheimo and Stiefel [59].

4.2. The Kustaanheimo-Stiefel transformation

The Kustaanheimo-Stiefel transformation is a particular case of the surjective application

$${\mathbb{R}}^4\to {\mathbb{R}}^3:(u_0,u_1,u_2,u_3)⟼(x_0,x_2,x_3),$$

(76)

defined by

$$\begin{array}{c}{x}_0=u_0^2-p{u}_1^2+t{u}_2^2-pt{u}_3^2,\hfill \\ x_2=2(u_0{u}_2+p{u}_1{u}_3),\hfill \\ x_3=2(u_0{u}_3+u_1{u}_2),\hfill \end{array}$$

(77)

subject to the constraint $\omega =0$ with

$$\omega =2(u_{1}d{u}_0-u_{0}d{u}_1+t{u}_{3}d{u}_2-t{u}_{2}d{u}_3),$$

(78)

where the parameters p and t can take the values $\pm 1$. The transformation introduced by Kustaanheimo and Stiefel [60] corresponds to $p=t=-1$, while the cases $p=-t=-1$ and $p=-t=1$ (or $p=t=1$) correspond to two other (inequivalent) transformations.

The applications of the transformations (77) and (78) range from number theory to physics (classical and quantum mechanics, gauge theories, etc.). The Hurwitz matrix corresponding to (77) and (78) is connected to the problem of factoring the sum of four squared numbers. In classical mechanics, the Kustaanheimo–Stiefel transformation is used for the regularization of the Kepler problem [60]. In quantum mechanics, the latter transformation enables one to transform the Schrödinger equation for the three-dimensional hydrogen atom (in an electromagnetic field) into a Schrödinger equation for a four-dimensional isotropic harmonic oscillator (with quartic and sextic anharmonic terms) subject to a constraint [61,62,63,64,65]. Conversely, the mappings (77) and (78) may be used in some dimensional reduction process for converting a dynamical system in ${\mathbb{R}}^4$ or ${\mathbb{R}}^2\times {\mathbb{R}}^2$ into a dynamical system in ${\mathbb{R}}^3$.

It is worth mentioning that Yoshida proposed a derivation of the Kustaanheimo–Stiefel in parabolic coordinates. Through the latter derivation, it becomes clearer where and how the additional dimension is introduced and what the bilinear relation means [66].

4.3. Integrability and similarity

A system with N degrees of freedom has at maximum $2N-1$ independent constants of motion. If there are N constants, the system is integrable. If there are $N+n$ ones, it is superintegrable and n times degenerated. If $n=N-1$, the system is totally degenerated, it is the case of the Keplerian and harmonic motions [67].

Prince and Eliezer have hown that the Laplace invariants are linked to the fact that the Kepler and harmonic motions are invariant by similarity, which explains respectively the third Kepler law and the isochronism of harmonic oscillations [68].

5. Invariant in the presence of an electric field: the Redmond invariant

Let us introduce, in the presence of an electric field $\mathbf{F}$, the so-called Redmond vector [16]:

$$\mathcal{R}=\mathbf{A}+\frac{1}{2}\left(\mathbf{r}\wedge \mathbf{F}\right)\wedge \mathbf{r}$$

(79)

and its projection on the z axis

$${\mathcal{R}}_z=A_z-\frac{1}{2}\left(x^2+y^2\right)F.$$

(80)

The Redmond vector is not an invariant [69]. However, the quantity

$$\mathcal{R}.\mathbf{F}=\mathbf{A}.\mathbf{F}+\frac{1}{2}\left[(\mathbf{r}\wedge \mathbf{F})\wedge \mathbf{r}\right].\mathbf{F}$$

(81)

which, using the relation

$${\displaystyle \mathbf{a}·(\mathbf{b}\wedge \mathbf{c})=\mathbf{b}·(\mathbf{c}\wedge \mathbf{a})-\mathbf{c}·(\mathbf{a}\wedge \mathbf{b})}$$

(82)

is equal to

$$\mathcal{R}.\mathbf{F}=\mathbf{A}.\mathbf{F}-\frac{1}{2}{\left(\mathbf{r}\wedge \mathbf{F}\right)}^2,$$

(83)

is an invariant. Indeed, for an electron, we have (since in atomic units the force $\mathbf{f}$ is related to the field $\mathbf{F}$ by $\mathbf{f}=-\mathbf{F}$):

$$\frac{d\mathbf{L}}{dt}=\mathbf{r}\wedge \mathbf{f}=-\mathbf{r}\wedge \mathbf{F}$$

(84)

and thus

$$\frac{d\mathbf{A}}{dt}=-\left[\mathbf{F}\wedge (\mathbf{r}\wedge \mathbf{v})+\mathbf{v}\wedge (\mathbf{r}\wedge \mathbf{F})\right]$$

(85)

i.e.

$$\frac{d\mathbf{A}}{dt}.\mathbf{F}=-\left[(\mathbf{F}.\mathbf{r})(\mathbf{F}.\mathbf{v})-\mathbf{r}.\mathbf{v}{F}^2\right].$$

(86)

We have also

$$\frac{1}{2}\frac{d}{dt}{(\mathbf{r}\wedge \mathbf{F})}^2=(\mathbf{r}\wedge \mathbf{F}).(\mathbf{v}\wedge \mathbf{F})$$

(87)

and using formula

$$(\mathbf{a}\wedge (\mathbf{b})\wedge (\mathbf{c}\wedge \mathbf{d})=(\mathbf{a}.\mathbf{c}).(\mathbf{b}.\mathbf{d})-(\mathbf{a}.\mathbf{d}).(\mathbf{b}\mathbf{c}),$$

(88)

we obtain zero [16]. The surface of the orbit precesses around the direction of the field, its angle with the field oscillates around a given average value.

This invariant is interesting or the study of the classical Stark effect, i.e. the motion of an electron in an electric field [70].

6. Oks’ supergeneralized Runge-Lenz vector for the two-center problem

Let us consider the problem of two Coulomb centers of charges Z and $Z^{\prime }$ separated by a distance R. This is an example of the so-called “two-center problem”. In the limit of large R, the problem of two Coulomb centers reduces to the problem of a hydrogenic ion in the nuclear charge Z in the uniform electric field $F=Z^{\prime }/R^2$, which is another fundamental problem of quantum mechanics (its analog in celestial physics is the one-center Kepler problem in a uniform gravitational field).

This simpler quantum has also a dynamical symmetry (which is linked to the fact that this problem allows the separation of variables in the parabolic coordinates). The corresponding integral of the motion in that case is the projection of a generalized Runge-Lenz vector on the internuclear axis. This generalized Runge-Lenz vector is the Redmond vector [16] (see Equation (79). After Redmond introduced it for this asymptotic case of the two-center problem, the challenge was to find a supergeneralization of the Runge-Lenz vector in the general (i.e. non necessarily asymptotic) case.

Throughout this section, we use atomic units $m=\hslash =e=1$ to conform with the convention of the authors. Gurarie, assuming a charge $Z_1$ at $z=a$ and a charge $Z_2$ at $z=-a$ [71], considered the following Hamiltonian:

$$H=\frac{p^2}{2}-\frac{Z_1}{r_1}-\frac{Z_2}{r_2},$$

(89)

and found the invariant vector

$${\mathbf{A}}^{\left(G\right)}=\mathbf{p}\wedge \mathbf{L}-Z_1\frac{\mathbf{r}-a{\mathbf{u}}_{\alpha }}{|\mathbf{r}-a{\mathbf{u}}_{\alpha }|}-Z_2\frac{\mathbf{r}+a{\mathbf{u}}_{\alpha }}{|\mathbf{r}+a{\mathbf{u}}_{\alpha }|},$$

(90)

where ${\mathbf{u}}_{\alpha }$ is the unit vector from $Z_1$ to $Z_2$. It can be checked that, at large R, the Redmond vector is recovered.

Krivchenko and Liberman considered the Hamiltonian [72]:

$$H=\frac{p^2}{2}-\frac{Z_1}{r_1}-\frac{Z_2}{r_2}+\frac{Z_2}{R}$$

(91)

with $Z_1$ at the origin and $Z_2$ at $z=R$. They obtained (using parabolic coordinates) the following vector

$${\mathbf{A}}^{\left(KL\right)}=\mathbf{p}\wedge \mathbf{L}-Z_1\frac{{\mathbf{r}}_1}{|{\mathbf{r}}_1|}+Z_2\frac{{\mathbf{r}}_2}{|{\mathbf{r}}_2|}+Z_2{\mathbf{u}}_{\alpha }$$

(92)

but its projection on the axis is not preserved. This led Kryukov and Oks to introduce the supergeneralized Runge-Lenz vector for the two-center Coulomb problem [17]:

$${\mathbf{A}}^{\left(KO\right)}=\mathbf{p}\wedge \mathbf{L}-\frac{L^2}{R}{\mathbf{u}}_{\alpha }-Z\frac{\mathbf{r}}{\left|\mathbf{r}\right|}-Z^{\prime }\frac{\mathbf{R}-\mathbf{r}}{|\mathbf{R}-\mathbf{r}|}+Z^{\prime }{\mathbf{u}}_{\alpha },$$

(93)

the z axis going from Z to $Z^{\prime }$.

As an application, Sanders and Oks used the robust perturbation theory [18,73] for calculating quadrupole corrections to the wave-functions. This is important for Stark effect, and the asymmetry of hydrogenic spectral lines in plasmas. If for the perturbed quantum system there is an operator $\mathbf{A}$ that commutes with the Hamiltonian H and the parts of these operators ${\mathbf{A}}_{\mathbf{0}}$ and $H_0$, characterizing the unperturbed quantum system, also commute, then the perturbation theory can be constructed in terms of the perturbation ($\mathbf{A}-{\mathbf{A}}_{\mathbf{0}}$) to the operator ${\mathbf{A}}_{\mathbf{0}}$, rather than in terms of the perturbation (H - $H_0$) to the operator $H_0$. For calculating corrections to the wave functions (which are common for both ${\mathbf{A}}_{\mathbf{0}}$ and $H_0$), the advantage is that the eigenvalues of the operator ${\mathbf{A}}_{\mathbf{0}}$ are typically non-degenerate (which is not the case for the eigenvalues of the operator $H_0$). Therefore, for calculating the first order corrections to the wave functions, it is sufficient to use the first order of the non-degenerate perturbation theory with respect to the perturbation ($\mathbf{A}-{\mathbf{A}}_{\mathbf{0}}$) and it would not involve infinite summations. In distinction, for calculating the same corrections in terms of the perturbation ($H-H_0$), one would have to proceed to the second order of the degenerate perturbation theory, involving infinite summations. Sanders and Oks chose the projection $A_z$ of the super-generalized Runge-Lenz vector, derived by Kryukov and Oks [17], on the axis connecting the nucleus of the hydrogenic atom/ion with the perturbing ion. The operator of the unperturbed projection $A_z^{\left(0\right)}$ has the well-known eigenvalues $q/n$. According to Equation (12) from [18], the first non-vanishing term of the expansion of the operator ($A_z-A_z^{\left(0\right)}$) in terms of the small parameter $n^2/R$ (here and below we use atomic units) is $-L^2/R$. Then, the corrections to the wave functions are given by

$$-\frac{1}{R}\frac{{\left(L^2\right)}_{n{q}_{p}m}^{n{q}_p^{\prime }m}}{A_{z,\alpha }^{\left(0\right)}-A_{z,{\alpha }^{\prime }}^{\left(0\right)}}=-\frac{n}{R}\frac{{\left(L^2\right)}_{n{q}_{p}m}^{n{q}_p^{\prime }m}}{(q_p-q_p^{\prime })}$$

(94)

where $q_p-q_p^{\prime }=\pm 2$ with $q_p=n_1-n_2$ and $q_p^{\prime }=n_1^{\prime }-n_2^{\prime }$ (not to be confused with the charge in Section 7.1). Sholin et al. have shown that [74]:

$$\langle n_1+1,n_2-1,m|L^2|n_1{n}_{2}m\rangle =-\hslash {\left[n_2(n-n_2)(n_1+1)(n-n_1-1)\right]}^{1/2}$$

(95)

and

$$\langle n_1-1,n_2+1,m|L^2|n_1{n}_{2}m\rangle =-\hslash {\left[n_1(n-n_1)(n_2+1)(n-n_2-1)\right]}^{1/2}.$$

(96)

The relation between parabolic and spherical wave-functions reads

$$|n_1{n}_2{m\rangle =(-1)}^{\frac{2n_2+\left|m\right|-m}{2}}|j_1{j}_2{m}_1{m}_2\rangle$$

(97)

where $j_1$ and $j_2$ are the angular momenta associated to operators

$${\mathbf{J}}_{\mathbf{1}}=\frac{1}{2}\left(\mathbf{L}+\mathbf{M}\right)$$

(98)

and

$${\mathbf{J}}_{\mathbf{2}}=\frac{1}{2}\left(\mathbf{L}-\mathbf{M}\right)$$

(99)

and one has

$$j_1=j_2=j=\frac{n-1}{2},m_1=\frac{m+n_2-n_1}{2}and{m}_2=\frac{m+n_1-n_2}{2}.$$

(100)

One defines then

$$J_{1\pm }=J_{1x}\pm i{J}_{1y},and,J_{2\pm }=J_{2x}\pm i{J}_{2y},$$

(101)

as well as

$$M_{+}=2J_{1+}-L_{+}=L_{+}-2J_{2+}$$

(102)

with $L_{\pm }=L_x\pm i{L}_y$ and

$$M_{-}=2J_{1-}-L_{-}=L_{-}-2J_{2-},$$

(103)

yielding

$$\mathbf{M}.\mathbf{M}=\frac{1}{2}\left(M^{+}{M}^{-}+M^{-}{M}^{+}\right)+M_z{M}_z.$$

(104)

Sholin et al. obtained

$$\langle n_1-1,n_2,m+1|M_{+}|n_1{n}_{2}m\rangle =\hslash \sqrt{n_1(n-n_1)},$$

(105)

$$\langle n_1,n_2-1,m+1|M_{+}|n_1{n}_{2}m\rangle =\hslash \sqrt{n_2(n-n_2)},$$

(106)

$$\langle n_1+1,n_2,m-1|M_{-}|n_1{n}_{2}m\rangle =\hslash \sqrt{(n_1+1)(n-n_1-1)},$$

(107)

and

$$\langle n_1,n_2+1,m-1|M_{-}|n_1{n}_{2}m\rangle =\hslash \sqrt{(n_2+1)(n-n_2-1)}.$$

(108)

We have also

$$\langle n_1-1,n_2,m+1|L_{+}|n_1{n}_{2}m\rangle =\hslash \sqrt{n_1(n-n_1)},$$

(109)

$$\langle n_1,n_2-1,m+1|L_{+}|n_1{n}_{2}m\rangle =-\hslash \sqrt{n_2(n-n_2)},$$

(110)

$$\langle n_1+1,n_2,m-1|L_{-}|n_1{n}_{2}m\rangle =\hslash \sqrt{(n_1+1)(n-n_1-1)},$$

(111)

$$\langle n_1,n_2+1,m-1|L_{-}|n_1{n}_{2}m\rangle =-\hslash \sqrt{(n_2+1)(n-n_2-1)},$$

(112)

$$\langle n_1,n_2,m|\mathbf{L}.\mathbf{L}|n_1{n}_{2}m\rangle =\frac{{\hslash }^2}{2}\left[n^2-{(n_1-n_2)}^2+{\left|m\right|}^2-1\right],$$

(113)

$$\langle n_1+1,n_2-1,m|\mathbf{L}.\mathbf{L}|n_1{n}_{2}m\rangle =-{\hslash }^2\sqrt{n_2(n-n_2)(n_1+1)(n-n_1-1)},$$

(114)

$$\langle n_1-1,n_2+1,m|\mathbf{L}.\mathbf{L}|n_1{n}_{2}m\rangle =-{\hslash }^2\sqrt{n_1(n-n_1)(n_2+1)(n-n_2-1)},$$

(115)

and

$$\langle n_1,n_2,m|\mathbf{L}.\mathbf{L}+\mathbf{M}.\mathbf{M}|n_1{n}_{2}m\rangle ={\hslash }^2\left(n^2-1\right).$$

(116)

As mentioned by Sanders and Oks, the non-diagonal matrix elements of the operator $L_x$ in parabolic coordinates were reproduced in Gavrilenko’s article [75]. In addition, the non-diagonal matrix elements of the operators $L_{\pm }$ can also be deduced from their proportionality (with the n-manifold) to the non-diagonal matrix elements of the operators $(x+iy)$. As a consequence, the non-diagonal matrix elements of the operator $L^2$ in parabolic coordinates can be obtained using their similar proportionality to the non-diagonal matrix elements of the operator $(x^2+y^2)$. The latter matrix elements were calculated, for instance, by Clark [76]. It is important to mention that Sanders and Oks [73] provided corrections to the Sholin tables from Ref. [77], which were used by many authors to calculate the asymmetry of hydrogenic spectral lines in plasmas.

7. Invariants in the presence of a magnetic field

7.1. The Landau-Avron-Sivardière approach

In the case of an external magnetic field, much less is known. Liberman does not provide any invariant in that case, but obtains series expansions, revealing some underlying symmetry [78]. The equation of motion of a charge q in the uniform field $\mathbf{B}$ is

$$\mu \frac{d\mathbf{v}}{dt}=q\mathbf{v}\wedge \mathbf{B}.$$

(117)

Integrating with respect to time, we get

$$\mu \mathbf{v}=q\mathbf{r}\wedge \mathbf{B}+\mathbf{k},$$

(118)

where $\mathbf{k}=\mu {\mathbf{v}}_0-\mathbf{B}\wedge {\mathbf{r}}_0$ is the time-invariant Landau vector (${\mathbf{r}}_0=\mathbf{r}(t=0)$ and ${\mathbf{v}}_0=\mathbf{v}(t=0)$). Combining Equations (117) and (118), one gets

$$\frac{d^2\mathbf{r}}{d{t}^2}+\frac{q^2{B}^2}{{\mu }^2}\mathbf{r}=\frac{q}{{\mu }^2}\mathbf{k}\wedge \mathbf{B},$$

(119)

which means that the motion of q is harmonic: the center C of the elliptic orbit is given by $\mathbf{OC}=(\mathbf{k}\wedge \mathbf{B})/\left(q{B}^2\right)$ and the orbit is in fact a circle since, according to Equation (117), the tangential acceleration is zero and the normal acceleration is constant. However [79], it is not necessary to know how to solve Equation (119). Multiplying Equation (118) by $\mathbf{B}$, we get where $\mathbf{v}$ is known from Equation (117). The nature of the motion is then found easily. Let us choose a new origin C with $\mathbf{OC}=\mathbf{a}$ such that the new Landau vector $\mathbf{K}$ is zero. We have, setting $\mathbf{r}=\mathbf{a}+\mathbf{R}$ and $\mathbf{v}=\mathbf{V}=d\mathbf{R}/dt$:

$$\mu \mathbf{V}=q\mathbf{R}\wedge \mathbf{V}.$$

(120)

and combining Equations (118) and (120):

$$\mathbf{a}=\frac{1}{q{B}^2}(\mathbf{k}\wedge \mathbf{B})={\mathbf{r}}_0+\frac{\mu }{q{B}^2}{\mathbf{v}}_0\wedge \mathbf{B}$$

(121)

and

$$\mathbf{a}\wedge \mathbf{k}=\frac{\mu v^2}{q{B}^2}\mathbf{B}+\mathbf{L},$$

(122)

where $\mathbf{L}=\mu \mathbf{r}\wedge \mathbf{v}+q{r}^2\mathbf{B}/2$ is the angular momentum $\mathbf{r}\wedge \mathbf{p}$, $\mathbf{p}=\mu \mathbf{v}+\mathbf{A}$ and $\mathbf{A}=(\mathbf{B}\wedge \mathbf{r})/2$ is the vector potential in the symmetric gauge. The invariance of $\mathbf{L}$ is a consequence of the rotational symmetry around $\mathbf{B}$ [80,81] and it is easy to show that

$$L=\frac{qB}{2}(a^2-R^2).$$

(123)

In the case of the motion of a charge in the field of a magnetic dipole, the motion may occur in the plane of symmetry of the dipole. If the motion is bound, the orbit fills a circular crown which axis carry the dipole: no privileged direction exists in the plane and no dynamical invariant can be found.

In the case of the motion of a charge in a magnetic monopole, there is a privileged direction in space different from the one of the cone axis: the direction of the periapsis. Such a direction is the one of a dynamical invariant which could be determined by Kerner [82,83,84]. The corresponding dynamical symmetry is the one of the $SO(4,2)$ group of matrices acting in the four-dimensional space and conserving the pseudo-norm $x^2+y^2-z^2-t^2$.

7.2. The Zeeman effect

The symmetry of the quadratic Zeeman effect was investigated by different authors [76,85,86,87]. In Ref. [88] on the motion of a charge in a magnetic field, Sivardière considers the Zeeman effect for the harmonic oscillator, and shows that the motion is epicyclic.

In a very interesting work, Yuzbashyan et al., motivated by the renewed interest in non-trivial degeneracies of a simple spin Hamiltonian, proposed a general approach to extract all non-trivial symmetries from the energy spectrum [89].

7.3. The spectrum of the hydrogen atom in electric and magnetic fields

The effect of a magnetic field on spectral line shapes, with or without an additional electric field, is an important topic of atomic physics (see for example the non-exhaustive list of references: [90,91,92,93,94,95,96,97,98,99]), for astrophysical or magnetic-fusion applications, both from the theoretical side and for the interpretation of Z-pinch experiments for instance.

A consistent treatment of the splitting of hydrogen atom in crossed weak electric and magnetic fields was given by Demkov et al. [100]. In that work, only the first-order perturbation theory was considered. The resulting formula for the energy corrections valid for arbitrary mutual orientation of the fields $\mathbf{F}$ and $\mathbf{B}$ was [101]:

$$E^{\left(1\right)}=n^{\prime }{e}_1+n^{\prime \prime }{e}_2,$$

(124)

where $n^{\prime }$ and $n^{\prime \prime }$ belong to the set $\left\{-j,-j+1,\cdots ,j\right\}$, with $j=(n-1)/2$ and $e_1$ and $e_2$ are absolute values (in atomic units) of the operator:

$${\mathbf{e}}_{\mathbf{1}}=\frac{\mathbf{B}}{2c}-\frac{3}{2}n\mathbf{F}$$

(125)

where c is the speed of light in vacuum, and

$${\mathbf{e}}_{\mathbf{2}}=\frac{\mathbf{B}}{2c}+\frac{3}{2}n\mathbf{F}.$$

(126)

The zero-order wave-functions are eigenvectors of $I_{1\alpha }$ and $I_{2\alpha }$:

$$I_{1\alpha }{\psi }_{n{n}^{\prime }{n}^{\prime \prime }}=n^{\prime }{\psi }_{n{n}^{\prime }{n}^{\prime \prime }}$$

(127)

and

$$I_{2\alpha }{\psi }_{n{n}^{\prime }{n}^{\prime \prime }}=n^{\prime \prime }{\psi }_{n{n}^{\prime }{n}^{\prime \prime }}$$

(128)

where

$$I_{i\alpha }=\frac{{\mathbf{I}}_{\mathbf{i}}.{\mathbf{e}}_{\mathbf{i}}}{e_i},$$

(129)

where the operators ${\mathbf{I}}_{\mathbf{1}},{\mathbf{I}}_{\mathbf{2}}$ are connected to $\mathbf{L}$ and $\mathbf{M}$ via

$${\mathbf{I}}_{\mathbf{1}}=\frac{\mathbf{L}+\mathbf{M}}{2}$$

(130)

and

$${\mathbf{I}}_{\mathbf{1}}=\frac{\mathbf{L}-\mathbf{M}}{2},$$

(131)

exactly as Equations (45) and (46). In general, first-order formula (124) completely removes the degeneracy. Second-order corrections lead only to small energy shifts. The latter was calculated by Solov’ev [102] using the dynamical symmetry group $O(4,2)$ of atomic hydrogen. The problem gets much more involved if the two fields are perpendicular ($\mathbf{F}\perp \mathbf{B}$, when $e_1=e_2=e$. The first-order correction $E^{\left(1\right)}$ depends then only on the set of quantum numbers $n^{\prime }+n=q$, implying a residual $n-$fold $(n-q)$ degeneracy. The second-order perturbation theory removes the degeneracy:

$$\begin{array}{ccc}\hfill E& =& -\frac{1}{2n^2}+E^{\left(1\right)}+\frac{F^2{n}^4}{16}\left[3q^2-17n^2-19-\frac{6}{1+{\gamma }^2}(n^2-3q^2-1)\right]\hfill \\ & & +\frac{H^2{n}^2}{16c^2}(3n^2+1-q^2+ϵ),\hfill \end{array}$$

(132)

where $\gamma =3ncF/H$ and $ϵ$ is the eigenvalue of the operator

$$b{\left(I_{1\alpha }-I_{2\alpha }\right)}^2-16I_{1\beta }{I}_{2\beta }.$$

(133)

In the above formula, $I_{1\beta }$ is the projection of operator ${\mathbf{I}}_{\mathbf{i}}$ on a direction belonging to the plane (${\mathbf{e}}_1$, ${\mathbf{e}}_2$) and orthogonal to ${\mathbf{e}}_i$. The parameter b is given by the formula

$$b={\gamma }^2-1-\frac{2}{1+{\gamma }^2}.$$

(134)

In fact, b grows monotonously when the ration $F/H$ is increased. The extreme values $b=-3$ and $b=\infty$ correspond to purely magnetic and electric perturbations (quadratic Zeeman and Stark effects). In these limits, the quantum number q turns into the magnetic quantum number m (in the case where $F=0$) or into $n_1-n_2$, where $n_1$ and $n_2$ are parabolic quantum numbers (in the case where $H=0$).

7.4. Magnetic monopole

The existence of isolated magnetic monopoles was first suggested by Dirac in 1931 [20] and could explain the quantization of the charge. Since there is still no known experimental or observational evidence that magnetic monopoles exist, they remain an active field of research. The magnetic field of the particle of charge q in $\mathbf{r}=r{\mathbf{u}}_{\mathbf{r}}$ is spherical and reads [84]

$$\mathbf{B}=g\frac{\mathbf{r}}{r^3}$$

(135)

and the equation of motion is

$$\mu \frac{d\mathbf{v}}{dt}=\delta \frac{\mathbf{L}}{\mu r^3},$$

(136)

where $\delta =-qg$, $\mathbf{L}=\mu \mathbf{r}\wedge \mathbf{v}$ is not an invariant, but the Poincaré vector $\mathbf{J}=\mathbf{L}+\alpha {\mathbf{u}}_r$ is one [103].

8. The charge-dyon system

A dyon is a hypothetical particle in 4-dimensional theories with both electric and magnetic charges. A dyon with a zero electric charge is usually referred to as a magnetic monopole. Many grand unified theories predict the existence of both magnetic monopoles and dyons. Dyons were introduced by Schwinger in 1969 as a phenomenological alternative to quarks [19]. He extended the Dirac quantization condition [20] to the dyon and used the model to predict the existence of a particle with the properties of the $J/\psi$ meson prior to its discovery in 1974. There is much interest in the problem of bound states of a fermion in the field of a fixed Dirac monopole or in a ’t Hooft-Polyakov monopole [104,105]. In theoretical physics, the ’t Hooft–Polyakov monopole is a topological soliton similar to the Dirac monopole but without the Dirac string. The allowed charges of dyons are restricted by the Dirac quantization condition. This states in particular that their magnetic charge must be integral, and that their electric charges must all be equal modulo one. It was shown that constant homogeneous electric field completely removes the degeneracy of the energy levels on orbital quantum number [106]. An important property is that the charge-dyon system possesses hidden symmetry, the corresponding group being $SO\left(4\right)$ for the discrete spectrum, and $SO(3,1)$ for the continuous spectrum. It was shown that the problem of the charge-dyon system due to hidden symmetry can be factorized not only in spherical but also in parabolic coordinates [106].

8.1. Classical approach of the charged monopole problem

In the case of the so-called “augmented charged monopole problem”, the monopole creates the potential

$$\frac{{\delta }^2}{2\mu r^2}.$$

(137)

Such an additive modification of the Coulomb potential by a repulsive inverse-square potential was introduced a long time ago [107,108] in order to enhance the dynamical symmetry. Thus, in addition to the Lorentz force, the monopole creates a central force $\mathbf{f}$ deriving from the potential $U\left(r\right)$, and the equation of motion reads

$$\mathbf{f}=\mu \frac{d\mathbf{v}}{dt}=\delta \frac{\mathbf{L}}{\mu r^3}-\frac{dU}{dr}{\mathbf{u}}_r.$$

(138)

Taking the vector product of the latter equation with $\mathbf{J}=\mathbf{L}+\delta {\mathbf{u}}_r$ and using

$$\frac{d{\mathbf{u}}_r}{dt}=\frac{1}{\mu r^2}\mathbf{L}\wedge {\mathbf{u}}_r,$$

(139)

we obtain

$$\frac{d\mathbf{v}}{dt}\wedge \mathbf{J}-\left(r^2\frac{dU}{dr}+\frac{{\delta }^2}{mr}\right)\frac{d{\mathbf{u}}_r}{dt}=\mathbf{0}.$$

(140)

If the coefficient of $d{\mathbf{u}}_r/dt$ is equal to a constant $\gamma$, that is if

$$U\left(r\right)=-\frac{\gamma }{r}+\frac{{\delta }^2}{2\mu r^2},$$

(141)

then a direct integration of Equation (140) provides the following invariant Laplace-like vector

$$\mathbf{A}=\mathbf{v}\wedge \mathbf{J}-\gamma {\mathbf{u}}_r=\mathbf{v}\wedge \mathbf{L}-\delta \frac{\mathbf{L}}{\mu r}-\gamma {\mathbf{u}}_r.$$

(142)

8.2. Quantum-mechanical MIC-Kepler problem

The integrable system MIC-Kepler was constructed by Zwanziger [107] and rediscovered by McIntosh and Cisneros [108]. This system is described by the Hamiltonian

$$H_0=\frac{{\hslash }^2}{2\mu }{\left(i\nabla +s\mathbf{A}\right)}^2+\frac{{\hslash }^2{s}^2}{2\mu r^2}-\frac{\gamma }{r},$$

(143)

where

$$\nabla \wedge \mathbf{A}=\frac{\mathbf{r}}{r^3}.$$

(144)

Its distinctive peculiarity is the hidden symmetry given by the following constants of motion

$$I=\frac{\hslash }{2\mu }\left[\left(i\nabla +s\mathbf{A}\right)\wedge \mathbf{J}-\mathbf{J}\wedge \left(i\nabla +s\mathbf{A}\right)\right]+\gamma \frac{\mathbf{r}}{r},$$

(145)

with

$$\mathbf{J}=-\hslash \left(i\nabla +s\mathbf{A}\right)\wedge \mathbf{r}+\frac{\hslash s\mathbf{r}}{r}.$$

(146)

These constants of motion, together with the Hamiltonian, form the quadratic symmetry algebra of the Coulomb problem. The operator $\mathbf{J}$ defines the angular momentum of the system, while the operator $\mathbf{I}$ is the analog of the Runge-Lenz vector. For the negative values of energy the constants of motion form the $so\left(4\right)$ algebra, whereas for positive values of energy, they build the $so(3,1)$ one. Due to the hidden symmetry the MIC-Kepler problem could be factorized in few coordinate systems, e.g., in the spherical and parabolic ones. Hence, the MIC-Kepler system is a generalization of the Coulomb problem in the presence of Dirac monopole. It can be considered as a one-parameter deformation family of the standard Kepler problem with the property that it retains its dynamical symmetries. Physically, the deformation parameter $\delta$ can be interpreted as the magnetic charge of the particle at rest and measures the pitch of the cone on which the trajectory lies. Actually, its genuine mathematical interpretation is as a cohomology class [109] of the symplectic structure [110]. The global symmetry group of the problem is either $SO\left(4\right)$, $E\left(3\right)$ or $SO(3,1)$, depending on whether the energy is negative, zero or positive. The motion obeys the three Kepler laws. The monopole number s satisfies the Dirac’s charge quantization rule, $s=0,\pm 1/2,\pm 1,\cdots$. The MIC-Kepler system could be constructed by the reduction of the four-dimensional isotropic oscillator by the use of the so-called Kustaanheimo-Stiefel transformation (see Section 4.2) both on classical and quantum mechanical levels [111,112,113,114]. In a similar way, reducing the two- and eight- dimensional isotropic oscillator, one can obtain the two- [115,116] and five-dimensional [117,118,119] analogs of MIC-Kepler system. An infinitely thin solenoid providing the system by the spin 1/2, plays the role of monopole in two-dimensional case, whereas in the five-dimensional case this role is performed by the $SU\left(2\right)$ Yang monopole [120], endowing the system by the isospin. All the above-mentioned systems have Coulomb symmetries and can be solved in spherical and parabolic coordinates both in discrete and continuous parts of energy spectra [121,122,123]. Finally, it is worth mentioning that there are generalizations of MIC-Kepler systems on three-dimensional sphere [124].

9. Conclusion

In this review, we discussed the importance of invariant vectors in atomic physics, such as the Runge-Lenz vector, the Redmond vector in the presence of an electric field, and the Landau-Avron-Sivardière vector when the system is subject to a magnetic field. The supergeneralized vector for the two-center problem is probably the most important recent result on the subject. The existence of constants of motion in the charge-dyon system is also briefly introduced. It is worth mentioning that Yoshida [125] discussed the main known ways of generalizing the Runge-Lenz vector [126,127]. De Castro Moreira generalized the Landau vector in some physical situations where it is possible to define a “potential” linear momentum [128]. From a more general viewpoint the author discussed the conditions for the introduction of the concept of a “generalized” linear (angular) momentum in analogy with the concept of mechanical energy. Oks reviewed classical studies of the oscillatory-precessional motion of an electron in the field of an electric dipole (the latter representing the polar molecule) with or without external magnetic or electric fields [129]. In a future work, we plan to discuss the case of the anisotropic harmonic oscillator, and especially its connections with the effect of an external electric field on the energy levels of an hydrogen atom [130,131].