Classical Dynamics of Muonic-Electronic Helium and Heliumlike Ions: The Allowance for the Eccentricity of the Muon and Nucleus Orbits

In the previous paper by one of us (hereafter paper I), the author considered Rydberg states of the muonic-electronic helium atom or helium-like ion and used the fact that the muon motion occurs much more rapidly than the electron motion. Assuming that the muon and nucleus orbits are circular, he ap-plied the analytical method based on separating rapid and slow subsystems. He showed that the electron moves in an effective potential that is mathematically equivalent to the potential of a satellite orbiting an oblate planet like the Earth. He also showed that the “unperturbed” elliptical orbit of the electron engages in two precessions simultaneously: the precession of the electron orbit in the plane of the orbit and the precession of the orbital plane of the electron around the axis perpendicular to the plane of the muon and nuclear orbits. The problem remained whether or not the allowance for the ellipticity of the orbit could significantly change the results. In the present paper, we address this problem: we study how the allowance for a relatively low eccentricity ε of the muon and nucleus orbits affects the motion of the electron. We derive an additional, ε-dependent term in the effective potential for the motion of the electron. We show analytically that in the particular case of the planar geometry (where the electron orbit is in the plane of the muon and nucleus orbits), it leads to an additional contribution to the frequency of the precession of the electron orbit. We demonstrate that this additional, ε-depen-dent contribution to the precession frequency of the electron orbit can reach the same order of magnitude as the primary, ε-independent contribution to the precession frequency. Therefore, the results of our paper seem to be important not only qualitatively, but also quantitatively.

In particular, in paper [6] the author considered Rydberg states of the muonic-electronic helium atom or helium-like ion and used the fact that the muon motion occurs much more rapidly than the electron motion. Therefore, he applied the analytical method centered on separating rapid and slow subsystems. He showed that the electron moves in an effective potential that is mathematically equivalent to the potential of a satellite orbiting an oblate planet (the Earth satellite being an example).
Further, in paper [6] it was shown that the "unperturbed" elliptical orbit of the electron engages in the following two precessions simultaneously: the precession of the electron orbit in the plane of the orbit and the precession of the orbital plane of the electron around the axis of symmetry of the muonic orbit. Despite these two precessions, the elliptical orbit of the Rydberg electron does not change its shape. The fact that the area of the elliptical orbit is conserved manifests the conservation of the square of the electron angular momentum. Thus, the system has higher than geometrical symmetry: indeed, from the geometrical symmetry (which is axial) followed only the conservation of the projection of the angular momentum on the axis of symmetry. This was a counterintuitive result of general physical interest.
In paper [6] the muon and nucleus orbits were considered circular. The problem remained whether or not the allowance for the ellipticity of the orbit could significantly change the results. In the present paper, we address this problem: we consider a more general situation where the muon and nucleus orbits are elliptical. We study how the allowance for a relatively low eccentricity ε of the muon and nucleus orbits affects the motion of the electron. We derive an additional, ε-dependent term in the effective potential for the motion of the electron.
We show analytically that in the particular case of the planar geometry (where the electron orbit is in the plane of the muon and nucleus orbits), it leads to an additional contribution to the frequency of the precession of the electron orbit.
We demonstrate that this additional, ε-dependent contribution to the precession frequency of the electron orbit can reach the same order of magnitude as the primary, ε-independent contribution to the precession frequency.

New Results
As in paper [6], in this study we analyze the system consisting of a muon, an their principal quantum numbers are 1 n µ  (for the muon) and 1 e n  (for the electron). The electron is much further way from the nucleus than the muondue to the large difference in the masses of the two leptons.
In this study we consider the case of low-eccentricity orbits of the muon and the nucleus. For a Coulomb potential: (where α = Ze 2 , e being the electron charge), the equation of the motion in the orbital plane is: Here: where ε is the eccentricity, (r, φ) are the polar coordinates, E is the energy, L in the angular momentum, and m r is the reduced mass of the subsystem "particle-Coulomb center". For low-eccentricity orbits ( 1 ε  ), we have: where r 0 is the radius of the circular orbit for ε = 0. As the case ε = 0 was analyzed in paper [6], the Rydberg electron perceives the rapid subsystem (the nucleus and the muon) as two uniformly charged rings of radii R nucl0 and R μ0 , where  (m nucl is the mass of the nucleus). The effective potential for the Rydberg electron in that case was: In our case of 1 ε  , the quantities R μ and R nucl are the following functions of time: is the frequency of revolution of the muon and of the nucleus about their center of mass, and is the reduced mass of the pair "nucleus-muon". Substituting Equation (6) where W is defined in Equation (9) Therefore, the complete effective potential in this case is: Next, we consider the orbits of the electron in the plane of the orbits of the muon and the nucleus, i.e., θ = π/2. It is easy to check, by looking at the first and second derivatives of the effective potential from Equation (14) with respect to θ, that this position corresponds to the stable equilibrium. In this case, the effective potential will be a Coulomb potential with two terms of the 1/r n -perturbation (from now on we will use the atomic units The calculation of the 1/r n -perturbation for the Kepler problem can be found, e.g., in [14] (the treatment for the cases n = 2 and n = 3 can be found also in the textbook [15]). For the Coulomb potential -α/r perturbed by the potential -β/r n+1 , the orbit undergoes a precession with the perihelion advance: ( ) where m is the reduced mass of the pair "particle-coulomb center", L is the angular momentum of the particle, ε e is the eccentricity of its orbit, and p = L 2 /(mα), in our case the particle being the electron. The ratio of the precession frequency to the Kepler frequency is the perihelion advance scaled by 2π, so for the second term in Equation (15) we obtain the following ratio of the precession frequency in the plane of the orbit to the Kepler frequency: The same result could be obtained by using Equation (1.7.10) from book [16], where the potential corresponds to the gravitational potential of the oblate Earth and is mathematically equivalent to Equation (13)  (the same treatment was used in paper [6]). Substituting the corresponding quantities into Equation (1.7.10) from book [16] and considering the case of θ = π/2, we obtain the same result as in Equation (17) (with e = 1 in the atomic units).
As an example, Figure 1 shows the plot of the ratio ( ) 1 pip K ω ω versus R μ0 for L = 3 and Z = 2. It is seen that in these ranges of R μ0 , the precession frequency Applying Equation (16) for the second perturbing term (the last term in Equation (15)), which corresponds to the effective potential due to the low eccentricity of the muon-nucleus orbits, we obtain the following additional contribution to the precession frequency (scaled by the Kepler frequency ω K of the electron):  is the absolute value of the scaled dimensionless energy of the electron. For the bounded motion of the electron, the eccentricity of its orbit is ε e = (1 − E s ) 1/2 and 0 < E s ≤ 1 (E s is the squared ratio of the semi-minor axis to the semi-major axis of the unperturbed elliptical orbit). Figure 2 shows the plot of the function f(E s ) from Equation (20). As E s increases from 0 to 1, f(E s ) monotonically decreases from 429 to 64.
As an example, Figure 3 shows the plot of the ratio  Figure 4 shows the plot of K 21 versus S = R μ0 /L 2 and E s (defined in Equation (21)) for Z = 2 and ε = 0.02. Figure 5 presents the plot of the ratio K 21 versus S and Z for ε = 0.02 and       ω precession frequency, caused by a relatively small eccentricity of the muon orbit (and of the nucleus orbit), can reach the same order of magnitude as the primary contribution ( ) 1 pip ω . This was the sole purpose of the above illustrations.

Conclusion
We considered a situation where the muon and nucleus orbits in the "nucleusmuon-electron" system are elliptical-the situation is more general compared to paper [6], where the muon and nucleus orbits were set to be circular. For the case where the eccentricity ε of the muon and nucleus orbits is relatively small, we obtained an additional, ε-dependent term in the effective potential for the motion of the electron. By analytical calculations, we demonstrated that in the particular case of the planar geometry (where the electron orbit is in the plane of the muon and nucleus orbits), it leads to an additional contribution to the fre-quency of the precession of the electron orbit. We showed that this additional, ε-dependent contribution to the precession frequency of the electron orbit can reach the same order of magnitude as the primary, ε-independent contribution to the precession frequency. Thus, the results of our paper seem to be important not only qualitatively, but also quantitatively.