Review of Classical Analytical Results for the Motion of a Rydberg Electron around a Polar Molecule under Magnetic or Electric Fields of Arbitrary Strengths in Axially Symmetric Configurations

We review 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. The focus is on the most recent studies. In one study (at zero external field), it was shown that, generally, the oscillations being in the meridional direction and the precession being along parallels of latitude can take place on the same time scale—contrary to the statement from the previous literature. In another study, it was shown that a magnetic field enables new ranges of the bound oscillatory-precessional motion of the Rydberg electron and that in one of the new ranges, the period of the θ-oscillations has the non-monotonic dependence on primary parameter of the system. This is a counterintuitive result. In yet another study, it was shown that under the electric field there are two equilibrium circular states of a positive energy and one equilibrium state of a negative energy. The existence of the equilibrium state of the negative energy is a counterintuitive result since at the absence of the field, the bound state was possible only for the zero energy. Thus, it is a counterintuitive result that in this case the electric field can play the role of a stabilizing factor.


Introduction
The motion of an electron in the field of a polar molecule is one of the most fundamental problems in atomic and molecular physics-on par with the problem of a hydrogen atom. If the dimension of the molecule is much smaller than its separation from the electron, the molecule can be considered as the point-like dipole. Obviously, such system possesses the axial symmetry-just from the geometrical considerations. Fermi and Teller [1] were the first to study a mathematically an equivalent problem in 1947: the motion of a low-energy muon in the vicinity of a hydrogen atom. In the subsequent years, there were lots of other relevant studies-see, e.g., paper [2] and references therein.
When the electron is only weakly bound to the dipole, i.e., the electron is in the Rydberg state, the classical treatment is appropriate. For the case of the finite electric dipole, the following two works have been published. In paper [3], the authors provided an analytical treatment of this problem by using the elliptical coordinates (ξ, η, φ). They demonstrated that the bound motion is possible in some finite region in terms of the coordinates ξ and η for any finite dipole. The region has the shape of a torus of the revolution with respect to the axis of the dipole.
In paper [4], the authors analyzed the motion of a negative charge around a finite dipole by using the cylindrical coordinates. First, they investigated circular orbits where the negative charge revolves in a plane orthogonal to the dipole. (A particular application described in detail in paper [4] was for the motion of a muon in the field of a finite dipole consisting of a stationary proton and electron, but their mathematical results are also valid for the motion of an electron around any finite dipole.) The authors of paper [4] demonstrated the existence of stable circular orbits under the condition that the dipole size exceeds some threshold value, which coincided with the corresponding threshold value obtained in quantum-mechanics-based investigations of this system-see, e.g., paper [2] and references therein. Then, they removed the assumption of circular orbits and showed the existence of stable conic-helical orbits of the negative charge in such systems. The orbit is located on the surface of a cone (specifically, on a right frustum of the cone), the cone axis coinciding with the dipole axis.
It is important to emphasize the following. An electron around a finite dipole is a particular case of the general problem: a charge moving around two stationary Coulomb centers. This general problem possesses a higher than geometrical symmetry-the algebraic symmetry. The manifestation of this algebraic symmetry is the existence of an additional integral of the motion: the projection of the super-generalized Runge-Lenz vector on the geometrical axis of symmetry, the super-generalized Runge-Lenz vector being derived in paper [5].
For the situation where the dipole can be considered point-like (i.e., where the separation of the electron from the dipole is much greater than the dipole size, the corresponding classical motion was first studied in paper [6] in 1968. The author of paper [6] demonstrated that only for the zero total energy of the electron the bound motion can happen, and it occurs on the surface of a sphere. Specific analytical results were developed in paper [6] to the form of integrals, but no integration has been performed and no qualitative description of the classical motion was given (see some details below).
Later, in paper [7], the author considered on semicircular meridional orbits on a sphere. He assumed the light particle to be charged positively and noted that in this limiting case, the situation is mathematically identical to a pendulum.
Then, in paper [8], the author concentrated on the following two types of orbits. One of them was the same as considered in paper [7]. Another one was a circular orbit coinciding with some specific parallel of latitude. The author of paper [8] noted that in the general case, the motion is a combination of relatively large oscillations in terms of the polar angle θ and a relatively slow precession along the azimuthal angle ϕ.
The most detailed classical analytical solution for the motion of an electron in the field of a point-like dipole was presented in paper [9]. The author of paper [9] showed that, in the general case, the oscillations in the θ-direction and the precession in the ϕ-precessions can actually happen on the same time scale-consequently, the opposite statement from paper [8] was erroneous.
The author of paper [9] obtained numerous specific analytical results for a variety of particular cases. He also deduced a condition, under which the electron motion (being in general conditionally periodic) can turn out to be periodic, manifested by a closed trajectory of the electron. Details are presented in Section 2.
In Section 3, we present a more complicated problem where the system is subjected to a magnetic field B parallel or antiparallel to the dipole. This configuration retains the axial symmetry. Analytical results are obtained for an arbitrary strength of the magnetic field. The magnetic field enables additional ranges of the bound motion of the electron. In one of the additional ranges of the oscillatory-precessional motion, it turns out that the dependence of the period of the meridional oscillations on the primary parameter of the system is non-monotonic, which is a counterintuitive result.
In Section 4, we present the study of the situation where instead of the magnetic field, there is an electric field F parallel or antiparallel to the electric dipole. This configuration also retains the axial symmetry. Again, analytical results are obtained for an arbitrary strength of the electric field. The electric field makes possible one equilibrium state of a negative energy and two equilibrium states of a positive energy. While at the absence of the electric field the bound state existed only for the zero energy, the possibility of the equilibrium state of the negative energy constitutes a counterintuitive result. Indeed, in this case, the electric field is a stabilizing factor, while, usually, the electric field promotes the ionization of atomic systems. We provide analytical expressions for the equilibrium radius of the electron orbit, the polar angle of the orbit, and the electron revolution frequency.
Since the frequency of spectral lines emitted by the Rydberg electron coincides with the revolution frequency, then these analytical results make the experimental determination of the energy of the equilibrium states possible.

Classical Non-Circular Orbits at the Absence of External Fields
Following paper [9], the energy of the electron in the field of the point-like electric dipole (of the dipole moment D) located at the origin of the spherical polar coordinate reference frame (the z-axis being along the dipole) can be expressed as: where e and m are the absolute value of the electron charge and the electron mass, respectively. The energy does not depend on ϕ, so that the corresponding generalized momentum p ϕ (which is the z-projection of the angular momentum M z of the electron) is conserved: Of course, this is the consequence of the axial symmetry of the configuration under the consideration, and it enables the separation of the θand φ-motions.
In paper [6], it was demonstrated that the bound motion can occur only if E = 0 and r = const. In other words, the electron is confined to a sphere: only θ and φ remain the dynamical variables.
The differential equation for the θ-motion derived in paper [6] is as follows: Here: x = cosθ, From Equation (2) follows the differential equation for the φ-motion: So, the sequence of tasks can be to find x(t) from Equation (3) and then obtain φ (t) from Equation (5).
For the bound motion to be possible, the following polynomial in Equation (3) should be positive in some segment of the variable x within the interval from −1 to 1, as noted in paper [6]. The polynomial has a maximum equal to −K + 2/3 3/2 at x = 1/3 1/2 and a negative minimum equal to −K − 2/3 3/2 at x = −1/3 1/2 . The necessary condition for the existence of the bound motion is for the maximum to be positive. This condition translates into the following requirement [6]: K < 2/3 3/2 = K max (7) or equivalently (according to Equation (4)): Under condition (7), the x-motion (i.e., the θ-motion) could be confined between two positive turning points. In other words, the bound motion can happen only in the upper hemisphere [6].
The results below have been obtained in paper [9]. The positive turning points x 2 and x 3 (x 3 > x 2 ) of the motion in the x-coordinate are the real roots of the cubic equation: The explicit expressions for these roots are the following: In Figure 1, there are plots of both x 2 (K) and x 3 (K). Specifically, x 2 (K) is the lower part of the double-valued curve, and x 3 (K) is the upper part of the double-valued curve. These two parts intersect at K = K max = 2/3 3/2 , where x 2 (2/3 3/2 ) = x 3 (2/3 3/2 ) = 1/3 1/2 . After introducing a scaled dimensionless time τ as follows: Equation (3) can be rewritten as: Both the plus and minus signs should be considered in Equation (13)-contrary to the choice of only the minus sign in Equation (12) of paper [6]. Now we present important particular cases. For K = K max = 2/3 3/2 , the electron executes a circular orbit representing the parallel of latitude that corresponds to cosθ = 1/3 1/2 , so that θ = 0.9553 and rad = 54.74 degrees. The angular velocity of the electron is constant (see Equation (5)): so that the period is: According to Equation (4), p φ = (2KmeD) 1/2 , meaning that for K = K max = 2/3 3/2 , one has p φ = 2(meD) 1 Next, in the case of K << 1, Equations (9) and (10) can be rewritten as follows: In the particular case of K = 0 (that is, p φ = 0), the electron oscillates along a semicircle. The semicircle is positioned in a meridional plane in the upper hemisphere. In this particular case, Equation (13) can be integrated analytically. The result is the explicit dependence of the scaled time τ on x (i.e., on cosθ): F(α, q) being the incomplete elliptic integral of the first kind. The motion occurs in the range of x from 0 to 1. In this range, the right side of Equation (18) is real-despite the presence of the imaginary unit i in Equation (18). For x << 1, Equation (18) yields: (We note that there was a misprint in the corresponding Equation (19) from paper [9].) Figure 2 shows the time evolution of the electron as it moves from a point on the equator (x = 0 at τ = 0) to the north pole (x = 1 at τ = 2.622), then changes its geographical longitude of in the upper hemisphere by 180 degrees and moves to the equator (x = 0 at τ = 5.244), then back to the north pole (x = 1 at τ = 7.866), then changes again its geographical longitude of in the upper hemisphere by 180 degrees and moves to the equator (x = 0 at τ = 10.488). (We note that in paper [9], the corresponding time evolution was not described accurately.) in the usual units. For the general case of an arbitrary K from the segment [0, K max = 2/3 3/2 ], the scaled time depends on x as follows: the lower limit of the integration being the smaller of the two turning points. In Figure 3, the time evolution versus x during one period of the θ-motion, as x changes from x 2 (K) to x 3 (K) and back to x 2 (K) is presented for K = 0.3 (dotted line), K = 0.2 (dashed line), and K = 0.1 (solid line). The integral in Equation (22) can be calculated analytically for K << 1 by expanding it in Taylor series up to the terms~K: where: In Equation (25), 2 F 1 (a, b, c, z) is the hypergeometric function. The accuracy of this approximate analytical result for τ(x, K) is demonstrated in Figure 4   The scaled period T θ of the θ-motion can be represented as follows (in units of mr 4 /(2eD)): The scaled period T θ versus K is displayed in Figure 5. An explicit analytical result for the scaled period of the θ-motion can be obtained for K << 1, as follows: where functions f and g are defined by Equations (24) and (25), respectively.
Next is the analysis of the φ-motion. Equation (5) can be represented as follows: Then, the dependence of the angular variable ϕ on the angular variable x = cosθ can be derived by using the relation between dt and dx from Equation (3): (29) Figure 6 shows φ versus x during one half-period of the θ-motion for K = 0.1 (solid line), K = 0.2 (dashed line), and K = 0.3 (dotted line). One can see that the curve φ(x) becomes steeper and the change of φ over one half-period of the θ-motion slightly increases with the growth of the parameter K. The following explicit analytical result for φ(K, x) can be derived for K << 1: where: In Equation (31), F(α, q) is the incomplete elliptic integral of the first kind. In Equation (30), the fact that x 2 (K) ≈ K for K << 1 has been used.
The accuracy of the approximate analytical result for φ(K, x) from Equation (30) is demonstrated in Figure 7 for K = 0.02 (solid line), where it is compared to the corresponding exact result from Equation (29) (dashed line). The accuracy of the analytical result is very good, as shown. During one half-period of the θ-motion, the angular variable φ changes by the following amount: In the case where K is close to zero, one has ∆φ ≈ π/2, while for the case where K is close to K max = 2/3 3/2 , one has ∆ϕ ≈ π/2 1/2 , so that: The quantity ∆φ monotonically increases between these two limits as K increases.
In the general case, the precession of the electron along parallels of latitude is combined with the oscillations in the upper hemisphere in the meridional direction. In general, the amount ∆ϕ(K), by which the angle φ advances during one period of the θ-oscillation, is not equal to nπ/m, where n and m are relatively small integers. Therefore, in the general case, the electron trajectory is not a closed curve: The motion is conditionally periodic. However, in some special cases, where  Figure 8 shows the intersections of the curve ∆φ(K)/π with the horizontal lines, corresponding to n/m, for the above three pairs. Below are the details for all three cases. For n = 2, m = 3, the angle φ completes 2 full circles during 3 periods of the θ-oscillations. This occurs for K = 0.2103. The corresponding frequency of the θ-oscillations is 1.2703 in units of 2eD/(mr 4 ). In this case, the intensity of the spectral line emitted at the frequency of 2.5406 in units of 2eD/(mr 4 ), which is the double-frequency of the θ-oscillations, would be enhanced (compared to the intensity of the spectral line at the primary frequency of the θ-oscillations) due to the resonance with the triple-frequency of the ϕ-precession. For n = 3, m = 5, the angle φ completes 3 full circles during 5 periods of the θ-oscillations. This occurs for K = 0.0632. The corresponding frequency of the θ-oscillations is 1.2223 in units of 2eD/(mr 4 ). In this case, the intensity of the spectral line emitted at the frequency 3.6669 in units of 2eD/(mr 4 ), which is the triple-frequency of the θ-oscillations, would be enhanced (compared to the intensity of the second harmonic of the frequency of the θ-oscillations) due to the resonance with the fifth harmonic of the frequency of the ϕ-precession.
For n = 4, m = 7, the angle φ completes 4 full circles during 7 periods of the θ-oscillations. This occurs for K = 0.0310. The corresponding frequency of the θ-oscillations is 1.2105 in units of 2eD/(mr 4 ). In this case, the intensity of the spectral line emitted at the frequency 4.8420 in units of 2eD/(mr 4 ), which is the quadruple-frequency of the θ-oscillations, would be enhanced (compared to the intensity of the third harmonic of the frequency of the θ-oscillations) due to the resonance with the seventh harmonic of the frequency of the ϕ-precession.
Observing the enhanced intensity of the spectral lines at the frequencies indicated above would be the manifestation of the states of the truly periodic motion of the electron around the polar molecule.
Thus, it is seen that the φ-precessions can actually occur on the same time scale as the θ-oscillations. It is true also in the general case. Therefore, the statement from [8] that "in general the motion should consist of large oscillations with respect to the polar angle θ combined with a slow precession about the z-axis" does not seem to be correct.

Effects of a Magnetic Field of an Arbitrary Strength
In this section, we consider the same system as in Section 2 but under a uniform magnetic field B along the z-axis. Thus, the axial symmetry is preserved. The analytical results presented below were obtained for the arbitrary strength of the magnetic field.
In Equation (35): is the scaled magnetic field: Ω > 0 or Ω < 0 if B is parallel or antiparallel to D, respectively.
So, the two angular variables become separated.
Obviously, the left side of the first formula in Equation (48) is non-negatively defined, so that the polynomial y(x, K) should be positive or zero. On denoting the following: Equation (49) becomes: From Equation (52), it is easy to find out that the parameter K can be negative for the values of u from the interval (−2, −2/3). (The negative values of u correspond to the antiparallel configuration of vectors B and D.) This is a new range of motion (with respect to the coordinate θ) that was impossible at the zero field.
There are three cases for the x-motion (we remind the x = cosθ): A.
In the case of K < −2/3 3/2 , the x-motion is in the segment [−1, 1]. Obviously, now there are additional ranges of the x-motion compared to the no-field situation where the x-motion was allowed only in the segment [x 2 (K), x 3 (K)].
Let us now obtain results for some particular cases.
On introducing the following notations:  (49), it follows that K = 0 for the following two values of the scaled (to the frequency dimension) magnetic field:

K = 0. From Equation
, Here, the most interesting situation is the one where Ω = Ω 1 . The electron moves in a semicircular path in a meridional plane in the northern hemisphere for any p φ , while in the no-field case, this trajectory was valid only for the zero value of p φ .
For K = 0, in Section 2, the explicit dependence of the scaled time τ on x = cosθ was obtained, given by Equation (18). For the range 0 < x < 1 where the motion occurs, the right side of Equation (18) yields real values. We remind the reader that for the subcase of x << 1, Equation (18) is simplified to Equation (19). Now, we provide some further details (in addition to those given by Equations (22), (26) and (29)) for the general situation where K is arbitrary. Figure 9 displays the angular variable φ versus x (during the half-period of the x-motion) for K = 0.15 (solid line) and K = 0.3 (dashed line).  Table 1 presents the scaled period T θ (in units of mr 4 /(2eD)) of the x-motion versus the parameter K in the negative range of K. It is a counterintuitive result that the period of the x-motion depends on K non-monotonically.   Further details on the effects of the magnetic field on this system can be found in paper [10].

Effects of an Electric Field of an Arbitrary Strength
In this Section we consider the same system as in Section 2, but under a uniform magnetic field F parallel or antiparallel to the dipole. Thus, again the axial symmetry is preserved. The analytical results presented below are obtained for the arbitrary strength of the electric field.
The Lagrangian has the form L = m[(dr/dt) 2 + r 2 (dθ/dt) 2 + r 2 sin 2 θ (dφ/dt) 2 ]/2 + eDcosθ/r 2 − eFr cosθ, where m and e are the mass and the absolute value of the electron charge, respectively. The quantity F can be positive or negative depending on whether the electric field is parallel or antiparallel to the dipole. From one of the Lagrange equations of the following motion: where p φ is the projection of the angular momentum on the dipole From Equation (63), we find the following: dϕ/dt = p φ /(mr 2 sin 2 θ).
so that:  (72), the latter equation would depend only on r and θ. Together with Equation (70), it would constitute the system of two equations for the two unknowns (r and θ) in the general case.
For obtaining explicit analytical results, from now on we consider the situation where the electron revolves around the dipole in a circular orbit (concentric with the dipole) perpendicular to the dipole axis, so that r = const and θ = const. Then, from Equation (70) follows the interrelation between these two constant values: From Equation (64), it follows that: where the constant angular velocity is denoted by ω. From Equation (67) one obtains: mrω 2 sin 2 θ − (2eD/r 3 + eFr) cosθ = 0.
On substituting cosθ from Equation (73) in Equation (77), we find: On substituting cosθ from Equation (73) in Equation (76), we obtain: After equating to each other the right sides of Equations (78) and (79) and doing simple transformations, we obtain the following equation containing only one unknown quantity-the quantity r r 5 − br 2 + 3ab = 0, where we introduced the following notations: It is important to note that while the quantity a from Equation (81) is always nonnegative, the quantity b can be positive (if F > 0) or negative (if F < 0). Equation (80) can be also represented as follows: where we denoted the following: Physically, E s is the scaled dimensionless absolute value of the energy. While in Equation (80), the left side has a slightly simpler form than in Equation (82), the latter equation explicitly contains the scaled dimensionless absolute value of the energy E s , which turns out to be the physical parameter controlling (together with the parameter b) various solutions of Equation (82), as will be shown below.
We seek the solution of Equation (82) in the following form: On substituting Equation (84) in Equation (82), we obtain: Let us first consider the situation where b > 0, i.e., F > 0, meaning that the electric field points in the same direction as the dipole. Then, Equation (85) takes the following form: According to one of the theorems of algebra-namely, Descartes's rule of signs-for b > 0, the polynomial P(c) from Equation (86) has either two or zero positive roots (we are not interested in negative roots-they are unphysical since it should be r > 0). By equating the derivative of P(c) to zero, it is easy to find after a simple analysis that at c = 0 the polynomial has a local maximum P max = (4/3) E s 2 and that it has the absolute minimum at c = (2/5) 1/3 : P min = (2/5) 5/3 − (2/5) 2/3 + 4E s 2 /3 = 4E s 2 /3 − 0.326.
Obviously, for the existence of two different positive roots, there should be P min < 0, leading to the following requirement: At E s < E s,max , the two positive roots merge together at c = (2/5) 1/3 = 0.737-see Figure 11. It is well known that there are no exact analytical solutions for the roots of polynomials of the 5th (or higher) degree. However, we can find approximate analytical solutions for the case where E s << 1, as follows.
One of the positive roots corresponds to the situation where in Equation (86) the first term is much smaller than the other two terms. The truncated equation has the following form: −c 2 + 4E s 2 /3 ≈ 0, yielding the root: so that: r 1 ≈ 2E s |b| 1/3 /3 1/2 = 2|E|/(3 1/2 e|F|. It is easy to verify that at r = r 1 , the first term in Equation (86) is indeed much smaller than the two other terms, given that E s << 1.
On substituting r 1 from Equation (91) in Equation (78), we find the following expression for the square of the electron revolution frequency (we remind that we still consider the situation where F > 0): Obviously, ω 1 2 should be positive. Consequently, this solution has the physical meaning only for E > 0. Then, by substituting r 1 from Equation (91) in Equation (73), we obtain the corresponding value of cosθ: The second positive root corresponds to the situation where, in Equation (86), the third term is much smaller than the two other terms. The truncated equation has the form yielding c 2 ≈ 1. It is seen that the third term in Equation (86) is indeed much smaller than the two other terms since E s << 1. For obtaining a non-vanishing result for the square of the electron revolution frequency, one has to go to the next approximation. By substituting c = (1 + α), where |α| << 1, in Equation (86), we find the following: so that: On substituting r 2 from Equation (96) in Equation (78), we obtain: Since ω 2 2 should be positive, this second solution also has the physical meaning only for E > 0. Then, by substituting r 2 from Equation (96) in Equation (73), we obtain the corresponding value of cosθ: Now, we proceed to the situation where b < 0, i.e., F < 0. Then, Equation (85) takes the following form: According to Descartes's rule of signs, this equation has one and only one positive root for any E s . This is illustrated in Figure 12. For E s >> 1 or E s << 1, we can find an approximate analytical result for the root of Equation (99), as follows. We start from the case of E s >> 1. In this case, the second term in Equation (99) is much smaller than the two other terms (see the proof below). The truncated equation has the following form: yielding: For this value of c 3 , the ratio of the second term in Equation (99) to, e.g., the first term is 1/c 3 3~1 /E s 6/5 << 1, so that the neglect of the second term is justified. With the c 3 from Equation (101), we obtain: Obviously, in Equation (103) |cosθ 3 | >> 1 (since, in this this case, E s >> 1), so that this formal mathematical solution does not have a physical meaning.
(The values of c 4 and r 4 coincide with the values of c 1 and r 1 , respectively). It is again easy to verify that at r = r 4 , the first term in Equation (99) is indeed much smaller than the two other terms, given that E s << 1.
The expression for the square of the electron revolution frequency ω 4 2 in this case differs by the sign from ω 1 2 in Equation (92) because now F < 0: Since ω 4 2 should be positive, this solution has the physical meaning only for E < 0. Then, by substituting r 4 from Equation (106) in Equation (73), we obtain the corresponding value of cosθ: It should be emphasized that the results given by Equations (106)-(108) for F < 0 are quite remarkable. They correspond to the bound states of negative energies (E < 0). We remind readers that in the absence of the electric field (the situation studied in papers [6][7][8][9]), as well as in the presence of the magnetic field along the electric dipole (the situation presented in Section 3 of the current review), the only possible bound states were of zero energy. Thus, the electric field can stabilize the system into bound states of negative energies. This is a counterintuitive result. Indeed, while a magnetic field typically inhibits the ionization of atomic systems, the electric field promotes the ionization and is thus considered as a destabilizing factor. However, the solution presented by Equations (106)-(108) for F < 0 (i.e., for the electric field antiparallel to the dipole) demonstrates that the electric field can work as the stabilizing factor for the Rydberg electron in the field of a point-like electric dipole.
We note that the frequency of the electron revolution found above would manifest as the frequency of spectral lines emitted by the system. By measuring the frequency of the spectral lines and knowing the electric dipole moment D and the field strength F, one can deduce the energy E-see, e.g., Equation (107).

Conclusions
In Section 2 of this review, we followed paper [9] to present the analysis of the classical bound motion of a Rydberg electron around a polar molecule. It was demonstrated that, generally, both the θ-oscillations and the φ-precessions can actually occur on the same time scale, so that the statement to the contrary from work [8] seems to be incorrect.
Furthermore, the relation between the dynamical variable φ and the dynamical variable θ was derived in the form of a one-fold integral in the general situation. In addition, an explicit analytical result for φ(θ) was obtained for the situation where the parameter K, which is the scaled square of the projection p φ of the angular momentum on the axis of the electric dipole, is relatively small.
For the case of K = 0, corresponding to the motion of the electron along a meridional semicircle crossing the North pole, an explicit analytical result for the dependence of time t on θ was obtained. For the case of K = K max = 2/3 3/2 , where the electron goes in a circle, the circle being the parallel of latitude corresponding to θ = 0.9553 rad = 54.74 degrees, an explicit analytical result was derived for the period of the revolution.
Furthermore, in the general case, the time evolution of the dynamical variable θ and the period T θ of the θ-oscillations was derived in the form of a one-fold integral. This was complemented by the corresponding explicit analytical expressions for the case of K << 1.
In addition, in the general case, the change of the angular variable ϕ during one half-period of the θ-motion was derived in the form of a one-fold integral.
Finally, it was shown in detail at what values of the parameter K, the conditionally periodic motion of the electron would become truly periodic. A general condition for the truly periodic motion of the electron was obtained and three examples were provided of the corresponding values of the parameter K.
From the fundamental point of view, deriving all of the above analytical results was possible due to the axial symmetry of the system.
In Section 3, the same system from Section 2 was analyzed under a magnetic field parallel or antiparallel to the electric dipole, so that the axial symmetry was preserved. The analytical results were derived in paper [10] for the arbitrary strength of the magnetic field.
It was demonstrated that under the magnetic field, new ranges of the bound motion of the Rydberg electron open up. Especially important is that, as the magnetic field opened up several new ranges of the parameter K, in the range 0 < K < 2/3 3/2 , the period of the θ-oscillations has the non-monotonic dependence on K. This is a counterintuitive result.
Next, while at B = 0, the oscillations of the electron along a semicircle located in a meridional plane in the Northern hemisphere were possible only for the zero projection of the angular momentum p φ on the electric dipole; the magnetic field makes this motion possible for any value of p φ -provided that Ω = −2p ϕ /(mr 2 ), where Ω = eB/(mc) is the scaled magnetic field.
In addition, while at B = 0, the circular trajectory of the electron along the parallel of latitude at cosθ = 1/3 1/2 was possible only for two values of the projection of the angular momentum corresponding to p φ /(meD) 1/2 = ±2/3 3/4 = ±0.8774, a proper choice of the magnetic field strength enables this kind of motion for any value of p φ -see Equation (59).
In Section 4, the same system from Section 2 was analyzed under an electric field F parallel or antiparallel to the electric dipole, so that the axial symmetry was preserved. The analytical results were derived for the arbitrary strength of the electric field.
First, some results were obtained for the general case. Then, the focus was on the circular orbits of the electron, the orbits being concentric with the dipole axis and the orbital plane being perpendicular to the dipole axis. Analytical results were obtained for the revolution frequency, as well as for the polar angle and the radius of the orbit-for three different equilibrium states of the system: two states of a positive energy (i.e., when the electric field and the dipole are parallel) and one state of the negative energy (i.e., when the electric field and the dipole are antiparallel).
The existence of the equilibrium state of the negative energy is a counterintuitive result. Indeed, at F = 0, the equilibrium state was possible only for the zero energy. This result is counterintuitive because it defies the usual expectation is that an electric field promotes the ionization: In the above situation, the electric field was found to be the stabilizing factor.
By observing the spectral lines emitted by the system and measuring their frequency, it is possible to deduce the energy by using our analytical results for the electron revolution frequency.
All of the above classical analytical results, obtained for the arbitrary strength of the magnetic and electric fields, provide a physical insight in the complicated dynamics of a Rydberg electron around a polar molecule. These analytical results, being valid for the arbitrarily large strength of the magnetic and electric fields, cannot be obtained analytically by the perturbation theory. The statement that the analytical results are obtained for the arbitrary strength of the magnetic and electric fields refers to the system "an electron plus an electric dipole". When the dipole is represented by a polar molecule, the following clarification is in order. The external electric field induces an additional dipole moment in the molecule. In this case, in the corresponding formulas from Section 4, the dipole moment D should include the contributions from both the built-in dipole moment and the induced one.