A coupled two frequency Hill\'s equation is solved. Analytically approximate solution correct up-to first order is derived using modified Lindstedt-Poincare perturbation method. For a wide range of controlling parameters, we compare the numerical and analytical solutions. The solution is the first step towards developing a comprehensive understanding of the electrodynamics of charged particles in a combinational ion trap utilizing both electrostatic DC and RF fields along with a constant static magnetic field with prospects of confining antimatter such as anti-hydrogen for a reasonably long durations of time.
Introduction
Conclusion
For particle trajectory in x-y plane, the analytical approximate solution correct up-to first order, is derived for the coupled two frequency Hill\'s equation using modified Lindstedt-Poincare method. The analytical solution matches well with the numerical solution obtained by numerical simulating the system of coupled differential equations given in Eq. (3). The analytical solution has a limited number of harmonic terms, i.e., (?±?_1 ) and (?±??_1 ) terms, whereas the numerical solution encompasses the effect of all the harmonic terms which make up the complete solution. Therefore, the matching is observed for some range of controlling parameters only.
If the order of the analytical solution is increased, the range of operating parameters for which the two solutions coincide will widen. However, the derivation of such higher order terms will be mathematically challenging. It is important to see that the solution described by Eq. 16 and Eq. 17 will blow up when ??1.
To obtain single frequency solutions one can simply substitute q_r=0 and keep ? away from unity. In most of the practical settings [7,15], the value of ? is substantially higher than unity, a regime wherein the analytical solutions are a good match to the numerical solutions. Experience guides us that analytical solution correct up-to first and second order are usually sufficient to provide deeper insights to both individual particle as well as collective dynamics inside the trap [17-19].
The relevance of an analytical solution cannot be understated when one must study the collective dynamics inside such combinational traps. Since the fields are spatially linear in this set up, one must see if a distribution function can be constructed for the particles by the method of inversion [17]. It is well known that RF heating on account of applied RF fields will increase the temperature of the charged particles.
The analytical tracking of temperature variation for each species inside such a trap is therefore important [18-19]. Temperature can be evaluated as the second order moment of the distribution function. To the best of my knowledge, such analytical work on collective dynamics for combinational traps has not been undertaken. Going ahead in this direction will require us to choose some operating parameters for stable configuration. The analytical expressions for particle dynamics derived in this work assumes importance as a vital starting point.
Imperfections in electrode geometry of the trap introduce deviations from the quadrupole potential. It would be interesting to see if analytical solutions can be derived for particle dynamics in such a scenario. Study of nonlinear resonances, deviation in the values of secular frequencies, changes in the stability regimes of the dynamics are all very interesting problems that could be taken up as future work.
References
[1] W. Magnus, S. Winkler, Hill\'s Equations, Dover, New York (1979), p. 4
[2] G. W. Hill, On the part of the motion of lunar perigee which is a function of the mean motions of the sun and moon, Acta Math., 8, (1886), p. 1
[3] Edmond de Hoffmann, Vincent Stroobant, Mass spectroscopy:Principles and Applications, 2nd Edition, John Wiley & Sons Ltd, p. 65
[4] Karl W. Boer, Udo W. Pohl, Qantum Mechanics of electrons in crystal, Semiconductor physics, Springer, (2018)
[5] Christopher Gerry, Peter Knight, Introduction to Quantum Optics, Cambridge university press, (2004)
[6] W. Paul, Electromagnetic traps for charged and neutral particles, Rev. Mod. Phys., 69, (1990), p. 531
[7] N. Leefer, K. Krimmel, W. Bertsche, et al., Investigation of two-frequency Paul traps for anti-hydrogen production, Hyperfine Interact, 238, (2017), p. 12
[8] D. Offenberg et. al., Translation cooling and storage of protonated proteins in an ion trap at sub-kelvin temperatures, Phys. Rev. A, 78, (2008), p. 0061401
[9] G. B. Anderson et.al., Trapped anti hydrogen, Nature, 468, (2010), p. 673
[10] G. B. Anderson et. al., Confinement of antihydrogen for 1,000 seconds, Nature Phys., 7, (2011), p. 263401
[11] C. H. Story et.al., First laser controlled antihydrogen production, Phys. Rev. Lett., 93, (2004), p. 263401
[12] G. Gabrielse, Trapped antihydrogen in its ground state, Phys. Rev. Lett., 108, (2012), p. 113002
[13] G. Werth, V. N. Gheorghe, F. G. Major, Charged Particle Traps II, Springer Series on Atomic, Optical, and Plasma Physics, Vol. 54, Springer Berlin Heidelberg, Berlin, Heidelberg (2009), p. 354.
[14] J. Walz, S. Ross, C. Zimmermann, L. Ricci, M. Prevedelli, T. W. Hansch, Combined Trap with the Potential for anti-hydrogen production, Physical Review Letters 75, (1995), p. 3257
[15] C. J. Foot, D. Trypogeorgos, E. Bentinea, A. Gardnerc, M.Kellerc, Twofrequency operation of a Paul trap to optimize confinement of two species of ions, International Journal of Mass Spectrometry, 430, (2018), p. 117
[16] P. Amore and A. Aranda, Presenting a new method for the solution of nonlinear problems, Phys. Lett. A, 316, (2003), p. 218
[17] D. J. Berkeland, J. D. Miller, J. C. Bergquist, W. M. Itano, D. J. Wineland, Minimization of ion micro-motion in a Paul trap, Journal of Applied Physics, 83(10), (1998), p. 5025
[18] V. Saxena and K. Shah, Plasma dynamics in dual frequency Paul trap using Tsallis distribution, IEEE Transactions on Plasma Science, 46(3), (2018), p. 474
[19] V. Saxena and K. Shah, Time evolution of Tsallis distribution in Paul trap, IEEE Transactions on Plasma Science, 45(6), (2017), p. 918