Analysis of the fractal intrinsic quality in the ionization of Rydberg helium and lithium atoms
Zhang Yanhui1, †, Xu Xiulan1, Kang Lisha1, Cai Xiangji2, Tang Xu1
College of Physics and Electronics, Shandong Normal University, Jinan 250014, China
School of Physics, Shandong University, Jinan 250100, China

 

† Corresponding author. E-mail: yhzhang@sdnu.edu.cn

Abstract

We study the fractal rhythm in the ionization of Rydberg helium and lithium atoms in an electric field by using the semiclassical method. The fractal structures present a nested relationship layer by layer in the initial launch angles of the ionized electrons versus the escape time, which is defined as the rhythm attractor, and exhibit similar rhythm endings. The gradually enhanced chaotic regions of the escape time plots tend to broaden as the scaled energy increases. In addition, the fractal rhythm changes synchronously with the oscillations of the kinetic energy spectrum. We note that the intrinsic quality of the fractal rhythm is closely related to the kinetic energy distribution, that is, the inherent dynamic properties of the Hamiltonian equations have an impact on the fractal regularities. In addition, different ionizing closed trajectories exhibit iterate properties and the inherent beauty of symmetry. Our results and analysis can not only reveal new laws in the ionization of Rydberg atoms, but also promote the establishment of the dynamic mechanism of fractals.

1. Introduction

In recent years, nonlinear science, especially chaos theory, has revealed the universality of the fractal structure in physics. Chaotic phenomena exist in many research fields, such as the escape dynamics of particles in mesoscopic devices, the ionization of high excited atoms, and so on.[17] The interaction between lasers and matter is an important means to explore new laws of microstructures.[813] The emergence and development of atomic physics in extreme strong field conditions can not only reveal new physical phenomena and effects, promote the progress of the ultra-fast laser, but also enable people to truly understand molecules, atoms, and their internal properties.[1418]

With the advancement of laser technology and semiclassical theory, considerable attention has been paid to photoionization microscopy and the ionization dynamics of atoms in external fields.[2,6,7,1925] The chaotic ionization of hydrogen in electric and magnetic fields is described by the Poincaré map.[2629] Recently, with the development of chaos and fractal theory in physics, consideration has been given to fractal dynamics in many fields, for example, the fractal dynamics in quantum and semiclassical transport.[12,3035] The escape dynamics of trapped ultracold atoms and Bose–Einstein condensates always exhibit fractal features.[3638] What is more, diverse attractors exist in nonlinear systems, such as the limit cycle attractor, torus attractor, Lorenz attractor, strange attractor, and so on, which all reflect the diversity of the nature and the flexuosity in the process of evolution. The ionization of Rydberg atoms is one of the most interesting fields in chaos and fractal physics. It is notable that, for the ionization of multi-electron Rydberg atoms in external fields, the intrinsic quality of fractal remains unrealized. Hence, it is a challenge to find the correspondence between the fractal mechanism and the chaotic escape laws of electrons in ionization systems.

In the present paper, we study the intrinsic quality of the fractal in the ionization of Rydberg helium and lithium atoms in an electric field and define the rhythm attractor to analyze the fractal structures within the framework of the semiclassical closed orbit theory.[8,9] Moreover, we study the influence of the potential energy distribution on the chaotic extent of the ionized electrons and calculate the kinetic energy spectra, which show excellent synchronism to the fractal structures. The symmetrical and iterative properties of the typical closed trajectories are discussed to reveal the dynamic mechanism of fractals within the ionization process of Rydberg atoms.

The paper is organized as follows. In Section 2, we briefly review the ionization process and the potential energy distribution. In Section 3, we discuss the relationship between the fractal rhythm of the escape time plots and the kinetic energy spectra, and the iterative and symmetrical features of the classical closed trajectories for both helium and lithium atoms. Section 4 presents the conclusions of the present paper.

2. Semiclassical closed-orbit theory
2.1. Ionization mechanism and the potential

Due to the fact that Rydberg atoms are very sensitive to external fields, we consider helium Rydberg atoms (with only one excited electron) and lithium atoms placed in an electric field directed along the z-axis, respectively. Within the framework of the closed-orbit theory, the outgoing wave packet which propagates away from the atomic core can be described as the evolution of ensemble electron trajectories.[8,9,13] In cylindrical coordinates (ρ,z), the ionization schematic view of the Rydberg system is depicted in Fig. 1. When a laser pulse strikes the atom, the excited electron moves outward from the nucleus in all directions, following the classical trajectories and being turned around by the external field, which are shown by the black solid lines in Fig. 1. Eventually the escaping electrons form a concentric wave packet interference pattern on the detector plane.

Fig. 1. Schematic of some classical trajectories in the process of the ionization of a Rydberg atom which is placed in an electric field F (or magnetic field B) that is oriented along the z axis and perpendicular to the detector plane. If the ionized electron absorbs a photon, the electron moves away from the nucleus with its initial outgoing angle, and the waves propagate outward in all directions. Some trajectories head directly downhill, other trajectories head uphill initially and then return to the nucleus. Thus, the outgoing wave can be described as an ensemble of classical trajectories.

The Hamiltonian of highly excited electrons in helium atoms with an electric field oriented along the z-axis can be written as where F is the electric field strength, and the z component of the angular momentum is taken to be zero. is the Hüpper model potential written as[4,19] where Z = 2 is the nuclear charge, and a = 1/3 is the parameter of the model potential determined by the range of the ion core for the helium Rydberg atom.[19]

By the scaling law , , , and , the transformed Hamiltonian can be written as The Coulomb singularity of the above Hamiltonian can be eliminated through adopting semiparabolic coordinates (u,v) and conjugate momenta which are respectively given by The new scaled time variable is , and the regularized Hamiltonian can be given by introducing an effective form Atomic units are used here and in the following. Clearly, the model potential is controlled by the scaled energy ε. By using analogous scaling laws[2729], for the lithium Rydberg atom in electric fields, the regularized Hamiltonian of the electron has a similar form.[19]

In order to show how the electrons escape in the ionization process, we depict the potential energy distribution of the Rydberg helium and lithium atoms at the same scaled energy in Figs. 2(a) and 2(b), respectively. The atoms are located at the center of the potential u = v = 0, the potential fields at the same coordinate range exhibit different structures. There are two obvious saddle points on the left and right sides of the counter plot of helium in Fig. 2(a), while the saddle points of lithium are not so obvious, which indicates that the potential well of lithium is relatively deeper and the core scattering effect is stronger than that of helium. Therefore, the ionized electrons escaping from the potentials will show different dynamical and chaotic features. The form of the Hamiltonian given in Eq. (5) is symmetric with respect to the coordinates u and v, which denotes that the electronic dynamics for the Rydberg atoms in the external fields will present unique properties, as discussed in the following.

Fig. 2. Contour plots of the potential energy distribution of Rydberg (a) helium and (b) lithium atoms. The ionized electrons begin at the nucleus located at u = v = 0, and eventually pass over the saddle region, and escape from the well.
2.2. Fractal and chaotic attractors

Three successive magnifications of the escape time plots (icicle map) for the helium atom are shown in Figs. 3(a)3(c). Figure 3(b) is the magnification of the angle interval between 1.4865 and 1.51737 of Fig. 3(a). Figure 3(c) is the magnification between 1.4946 and 1.50175 of Fig. 3(b). By analyzing the three different colored lines, the nested relationship layer by layer from Fig. 3(a) to Fig. 3(c) depicts the rhythm properties in the escape time plots. The fractal dimension of the self-similarity structure can reflect the chaotic properties of the ionized electrons, which is calculated within the framework of the box-counting method.[20] We define the phenomena in the ionization systems as the rhythm attractor with the fractal dimension which expresses the irregularity of the rhythm structures. These curves exhibit the same rhythm at the endings of the plots, which are called “rhythm endings” denoting the universality in the ionization process.

Fig. 3. Fractal self-similarity structures in the process of the ionization of a Rydberg helium atom placed in an electric field F = 19 V/cm. Three successive magnifications of the escape time plots: panel (b) is the magnification of the assigned angle interval of the fractal structures in panel (a), and panel (c) is the magnification of the assigned angle interval of the fractal structures in panel (b).
3. Results and discussion
3.1. Fractal and chaotic structures of the escape time plots

In order to identify the rhythm properties in the ionization process, we compute the Hamiltonian motion of the Rydberg helium and lithium atoms based on the canonical equations, which are written as With the initial conditions , , , and of the lithium atom, we perform iterative calculations, obtain and search the ionization trajectories in any outgoing angles. The iterative calculation has been widely used in researching the transport and ionization problems in the semiclassical field.

In Fig. 4, we compare the escape times and the initial launch angle ranges at different scaled energies for the lithium atom.

Fig. 4. Escape time plots for different energies ( ) of lithium atoms in an electric field F = 19 V/cm. The ionization process shows chaotic characteristics and the dotted rectangles exhibit the same rhythm endings. The chaotic complexity regions become narrower with the decrease of the scaled energy.

We use the escape time curve, which is called the icicle map,[2729] to provide more detailed insight into the connection between the fractal structure and the escape processes of the ionized electrons. As shown in Figs. 4(a)4(e), with the decrease of the scaled energy, the fractal icicle distribution shifts toward larger initial launch angles with the maximum escape time 50, and the chaotic complexity regions become narrower at the same time. This implies that the ionized probability decreases when the scaled energy decreases, and the electrons of larger ejection angles will be easier to reach the detector plane exhibiting more chaotic escape properties. These escape time plots contain the rhythm attractors and the dotted rectangles exhibit the same rhythm endings. Moreover, there are fractal self-similarity structures in the various scale ranges.

Considering the core scattering effect on the chaotic escape of electrons for the lithium atom, we further study the ionization of the helium Rydberg atom in a pure electric field at the scaled energy ε = −1.3. As shown in Figs. 5(a) and 5(b), the escape time plots display the same fractal self-similarity characteristics. Some electron trajectories head directly downhill to the detector plane within a short time, while other trajectories are turned around by the influence of the electric field. However, the chaotic regions in Fig. 5(b) show more sophisticated structures than those in Fig. 5(a), which suggests that the core scattering effect causes the enhancement of the ionization probability for the lithium atom under the same external conditions. The position of the fractal self-similarity region A in Fig. 5(a) is located in a finite interval and region B in as shown by the dashed line regions, respectively.

Fig. 5. For Rydberg atoms in an electric field F = 19 V/cm, the escape time that electrons take to strike the detector is plotted as a function of the initial launch angle at the same scaled energy ε = −1.3. (a) The escape time plot of the helium atom, chaotic region A is the typical fractal self-similarity region. (b) The escape time plot of the lithium atom, chaotic region B is the typical fractal self-similarity region. The dotted rectangles exhibit the same rhythm endings.

We next calculate the kinetic energy spectrum of families of ionizing trajectories. Figures 6(a) and 6(b) correspond to region A of Fig. 5(a). Figures 6(c) and 6(d) correspond to region B of Fig. 5(b). Clearly, the fractal rhythm in Figs. 6(a) and 6(c) and the kinetic energy oscillations in Figs. 6(b) and 6(d) exhibit excellent synchronous variations, respectively. The regions marked by the red lines from a1 to e1 in Fig. 6(a) correspond exactly to the chaotic positions of the escape time plot, the non-chaotic region between a1 and b1 corresponds to the non-chaotic kinetic energy oscillation. The regions marked by the blue lines from a2 to g2 in Fig. 6(c) show that the oscillations of the kinetic energy spectrum of Fig. 6(d) are changing synchronously with the fractal rhythm of different ionized electrons of Fig. 6(d). The intrinsic quality of the fractal rhythm is closely related to the energy distribution, that is, the inherent dynamical properties of the Hamiltonian affect the fractal regularities.

Fig. 6. (a) The fractal structures and (b) the kinetic energy spectra of the helium atoms; the kinetic energy oscillations exhibit excellent synchronous variations with the fractal rhythm as shown in the regions marked by red lines from a1 to e1. (c) The fractal structures and (d) the kinetic energy spectra for the lithium atoms; the kinetic energy oscillations exhibit excellent synchronous variations with the fractal rhythm as shown in the regions marked by red lines from a2 to g2.

Moreover, figures 6(c) and 6(d) of the lithium atom illustrate enhanced chaotic properties compared to those of the helium atom shown in Figs. 6(a) and 6(b), implying that the electronic dynamics is more complex, which is caused by the deeper ionization potential well than that of the helium atom. Since one point presents one classical trajectory in the escape time plots, and the plots show fractal regularity structures related to the kinetic energy spectrum, we shall see that classical trajectories in this system must contain certain properties, which will be discussed next.

3.2. Symmetry analysis of classical ionizing electron trajectories

To get a deeper insight into the fractal and the classical trajectories, we calculate the closed trajectories of the ionized electrons of the helium and lithium atoms.

In principle, all the initial outgoing angles should be examined, but we find that the Hamiltonian has a symmetry about z and the ionized electrons escape from the atom in all directions. Therefore, the trajectories can be obtained to research the fractal rhythm in the ionization process. We calculate the fractal dimensions of the escape time plots at different scaled energies for both the helium and the lithium atoms to study the escape behavior of the ionized electrons. For the helium atom, the fractal dimensions of the escape time plots are at different scaled energies ( ), respectively. The results indicate that the fractal dimensions decrease with the scaled energy, at the same time, the ionization of the helium atom shows enhanced chaotic characteristics. A set of ionizing closed trajectories of the helium and lithium atoms at the scaled energy ε = −1.3 are shown in Fig. 7. For the helium atom, from Fig. 7(a) to Fig. 7(c) we find that the three trajectories have similar iterative rules and present interesting rhythm properties. Moreover, the evolutional trends in Figs. 7(d)7(f) tend to complexity step by step with the increasing number of crossings of the z-axis. The fractal dimensions of the escape time plots for the lithium atom in Figs. 4(a)4(e) are , respectively. The fractal dimensions decrease with the decreasing scaled energy, moreover, the ionization of the lithium atom exhibits enhanced chaotic properties compared with that of the helium atom. In Fig. 7, the ionizing closed trajectories show perfect symmetry and iterative characteristics for the lithium atom, just like those of the helium atom. With the fractal dimensions calculated above, at the scaled energy ε = −1.3, Df of the lithium atom is larger than that of the helium atom, which demonstrates that the trajectories of the lithium atom are more complex than those of helium. Because of the stronger core effects of the lithium atom, the dynamical motion of the ionized electrons is iterated numerically in the regularized coordinates, and iterate more and more times from Fig. 7(g) to Fig. 7(l). With a small variation of the initial launch angles, the trajectories show great disparities,[26,27] implying that the ionization of the Rydberg atom satisfies the basic rule of the chaotic system: sensitivity to the initial conditions.

Fig. 7. Important typical types of closed trajectories of (a)–(f) helium and (g)–(l) lithium atoms in an electric field. Different trajectories are calculated for the two atoms, the panels (a)–(l) show the similar iterative rules and interesting rhythm properties.

Such rules of closed trajectories and fractal structures occur in many angles, which indicates the scale-free features in the ionization process. The symmetrical shapes of the closed trajectories and the rhythm reflect that the Hamiltonian equations influence the intrinsic dynamics of the ionization systems.

4. Conclusion

In summary, we have used the semiclassical closed-orbit theory to calculate the ionization fractal rhythm of Rydberg helium and lithium atoms in an electric field. The fractal structures exhibit the nested relationship layer by layer which is defined as the rhythm attractor and present similar rhythm endings in the initial launch angles versus the escape time. Due to the energy distributions of the potentials being distinct, the scattering effects of the helium and lithium atoms show different dynamical impacts on the electrons. In addition, the ionization chaotic regions of the escape time plots tend to broaden as the scaled energy increases for the Rydberg atoms. Oscillations in the kinetic energy spectra of the helium and lithium atoms show excellent correspondence with the fractal structures, which illustrates that the intrinsic dynamic properties of the Hamiltonian equations are closely related to the fractal rhythm. The typical ionizing closed trajectories have been calculated and show iterate properties and perfect symmetry. We believe this work will provide more research methods for revealing new laws in the ionization of Rydberg atoms.

Reference
[1] Yang G C Mao J M Du M L 1999 Phys. Rev. 59 2053
[2] Zhao L B Delos J B 2010 Phys. Rev. 81 053417
[3] Zhao L B Delos J B 2010 Phys. Rev. 81 053418
[4] Hüpper B Main J Wunner G 1996 Phys. Rev. 53 744
[5] Lankhuijzen G M Noordam L D 1996 Phys. Rev. Lett. 76 1784
[6] Zhou H Li H Y Gao S Zhang Y H Jia Z M Lin S L 2008 Chin. Phys. 17 4428
[7] Wang D H 2011 Chin. Phys. 20 013403
[8] Du M L 1989 Phys. Rev. 40 4983
[9] Du M L Delos J B 1988 Phys. Rev. 38 1896
[10] Du M L Delos J B 1987 Phys. Rev. Lett. 58 1731
[11] Fabian L Iva Březinová Joachim B Florian L 2013 Phys. Rev. 88 022916
[12] Zhao G P Zhang Y H Cai X J Xu X L Kang L S 2016 Physica 84 10
[13] Yang B C Delos J B Du M L 2014 Phys. Rev. 89 013417
[14] Blondel C Delsart C Dulieu F 1996 Phys. Rev. Lett. 77 3755
[15] Nicole C Offerhaus H L Vrakking M J J Lépine F Bordas C 2002 Phys. Rev. Lett. 88 133001
[16] Lépine F Bordas C Nicole C Vrakking M J J 2004 Phys. Rev. 70 033417
[17] Nicole C Sluimer I Rosca-Pruna F Warntjes M Vrakking M Bordas C Texier F Robicheaux F 2000 Phys. Rev. Lett. 85 4024
[18] Ivanov I A Nam C H Kim K T 2016 Phys. Rev. 93 043404
[19] Hüpper B Main Jörg Wunner G 1995 Phys. Rev. Lett. 74 2650
[20] Xu X L Zhang Y H Cai X J Zhao G P Kang L S 2016 Chin. Phys. 25 110301
[21] Liu F L Zhao L B 2017 Phys. Rev. 95 043428
[22] Deng M Gao W Lu R Delos J B You L Liu H P 2016 Phys. Rev. 93 063411
[23] Stodolna A S Lépine F Bergeman T Robicheaux F Gijsbertsen A Jungmann J H Bordas C Vrakking M J J 2014 Phys. Rev. Lett. 113 103002
[24] Chen X Wu Y Zhang J T 2017 Phys. Rev. 95 013402
[25] Wang L Yang H F Liu X J Liu H P Zhan M S Delos J B 2010 Phys. Rev. 82 022514
[26] Mitchell K A Handley J P Tighe B Flower A Delos J B 2004 Phys. Rev. 70 043407
[27] Mitchell K A Handley J P Tighe B Flower A Delos J B 2004 Phys. Rev. Lett. 92 073001
[28] Mitchell K A Handley J P Tighe B Delos J B Knudson S K 2003 Chaos 13 880
[29] Mitchell K A Handley J P Delos J B Knudson S K 2003 Chaos 13 892
[30] Kotimäki V Räsänen E Hennig H Heller E J 2013 Phys. Rev. 88 022913
[31] Yi X H Bi J J Liu R Wang C K Jiao Y Li Z L 2016 Chin. Phys. 25 128503
[32] Byrd T A Delos J B 2014 Phys. Rev. 89 022907
[33] Cai X J Zhang Y H Li Z L Jiang G H Yang Q N Xu X Y 2013 Chin. Phys. 22 020501
[34] Mackay R S Meiss J D Percival I C 1984 Physica 13 55
[35] Zhang Y H Cai X J Li Z L Jiang G H Yang Q N Xu X Y 2013 Chin. Phys. Lett. 30 040501
[36] Mitchell K A Steck D A 2007 Phys. Rev. 76 031403
[37] Friedman N Kaplan A Carasso D Davidson N 2001 Phys. Rev. Lett. 86 1518
[38] Mitchell K A Delos J B 2006 Physica D: Nonlinear Phenomena 221 170