Home Nonsequential double ionization channels control of CO2 molecules with counter-rotating two-color circularly polarized laser field by laser wavelength
Article Open Access

Nonsequential double ionization channels control of CO2 molecules with counter-rotating two-color circularly polarized laser field by laser wavelength

  • Zhenjie Ge , Lihua Bai EMAIL logo , Xu Su and Keying Liu
Published/Copyright: October 7, 2023

Abstract

Using the classical ensemble model, we investigate the effect of laser wavelength on the electron dynamics process of nonsequential double ionization (NSDI) for linear triatomic molecules driven by a counter-rotating two-color circularly polarized laser field. Based on the delay time between recollision and final double ionization, two particular ionization channels are separated: recollision-impact ionization (RII) and recollision-induced excitation with subsequent ionization (RESI). Numerical results show that with the increase of the laser wavelength, the triangle structure of the ion momentum distribution becomes more obvious, which indicates that the electron–electron correlation of NSDI is enhanced. In addition, we find that the ratio of the RESI channel gradually decreases with the increase of the laser wavelength, while the ratio of the RII channel is opposite. However, the dominant channel is still RESI. It means that the two ionization channels can be controlled effectively by changing the laser wavelength.

1 Introduction

With the development of laser technology, the interaction of atoms, molecules, and strong laser fields has generated many interesting nonlinear phenomena, such as above-threshold ionization [1,2,3], nonsequential double ionization (NSDI) [4,5,6], and high-order harmonic generation [7]. Theoretical and experimental work has been devoted to observing these physical phenomena. Among them, NSDI in strong field physics has attracted much attention [4,8]. Because the electron–electron correlation plays a significant role in the NSDI process, a large number of studies about NSDI are related to the electron–electron correlation [9], such as angular correlation [10,11,12] and recollision [13,14,15].

The NSDI process can be explained by a three-step model [16,17]. First, when an atom or a molecule is exposed to a strong laser field, an electron can be freed by tunneling ionization [18]. Then, the released electron is accelerated by the intense laser field, and it can be driven back to the parent ion when the laser field is reversed [18]. Therefore, the returning electron recollides with the parent ion in-elastically, resulting in NSDI [5,6,19]. Due to the recollision of the two electrons during the NSDI process, the two electrons are highly correlated. A large number of theoretical and experimental studies have been performed to explore the correlated dynamics of the two electrons in NSDI. Two recollision mechanisms, recollision-impact ionization (RII) [20,21] and recollision-induced excitation with subsequent ionization (RESI) [22], in the NSDI process have been found. These two recollision mechanisms can be distinguished by the delay time. When the energy of the returning electron is greater than the ionization potential energy of the bound electron, the bound electron will be ionized directly [23,24]. The process for the shorter delay time is defined as RII [18,25]. The bound electron cannot be ionized directly if the energy of the returning electron is less than the ionization potential energy of the bound electron [26,27] and the process for the longer delay time is defined as RESI [18,25].

The development of laser technology has also enabled people to achieve a lot of different laser pulses, such as elliptically polarized (EP) laser pulses, circularly polarized (CP) laser pulses, and two-color circularly polarized (TCCP) laser pulses [28]. Recently, it has been reported experimentally that the NSDI is studied in counter-rotating two-color circularly polarized (CRTC) laser fields [29,30,31,32]. The CRTC laser field is widely used in the research of NSDI due to its special electric field structure. In recent years, Zhang et al. have found that the CRTC laser field is beneficial in increasing the double ionization (DI) probability for O2 molecules [33]. Few-cycle laser pulses can simplify the NSDI process by achieving a single recollision event [21]. It has been found that the relative phase and the polarization of the TCCP laser field can affect the DI probability and the correlation effect [18,34,35]. Currently, we have found that NSDI in the CRTC laser fields for different laser wavelengths has not received much attention.

In this article, we investigate the effect of the laser wavelength of NSDI for CO2 molecules in the CRTC laser field with a classical ensemble method [36,37]. It is shown that the laser wavelength influences not only the yield of DI but also two particular channels (RII and RESI) in NSDI. In addition, we study the electron–electron correlation in NSDI at different laser wavelengths. Control between the two channels can be achieved by varying the laser wavelength.

2 Theoretical model

In this work, we apply the two-dimensional (2D) classical ensemble model to investigate the effect of the laser wavelength in NSDI by the CRTC laser field. The classical ensemble model has been widely used in strong laser field electron dynamics. The underlying recollision processes and mechanisms in NSDI can be presented intuitively by tracing the classical trajectories [38,39].

In our calculation, the CO2 molecule is arranged on the x-axis. The two oxygen atoms are placed at coordinates (−R,0) and (R,0). The Hamiltonian of the CO2 molecules without an electric field can be written as follows (in atomic units):

(1) H e = T ( p i ) + V ( r i ) ,

where the kinetic energy T and potential energy V of the CO2 molecules can be represented as follows:

(2) T ( p i ) = i = 1 , 2 p i 2 2 ,

(3) V ( r i ) = i =1,2 2 / 3 x i 2 + y i 2 + a c 2 i = 1 , 2 2 / 3 ( x i R ) 2 + y i 2 + a o 2 i = 1 , 2 2 / 3 ( x i + R ) 2 + y i 2 + a o 2 + 1 ( x 1 x 2 ) 2 + ( y 1 y 2 ) 2 + b 2 ,

the Hamiltonian of the CO2 molecules is the total energy without an external electric field (the sum of the first ionization potential energy and the second ionization potential energy). Here, p i and r i represent the momentum and the position of the two electrons [34]. The internuclear R = 2.19 a .u . . In this article, a c = 1 , a o = 0.8 , and b = 0.05 ( a c and a o are the softcore coefficients used to adjust the bond potential of the atom, shielding the attraction of the atomic nucleus to electrons, and b is the softcore coefficient used to adjust the Coulomb potential between electrons, that is, the Coulomb repulsion of electrons to electrons). Then, we switch on the external laser field. The Hamiltonian of the two electrons in the laser field is expressed as follows:

(4) H = H e + ( r 1 + r 2 ) E ( t ) ,

where E ( t ) = E r ( t ) + E b ( t ) is the CRTC laser field, where E r ( t ) is the fundamental (red) laser pulse, and E b ( t ) is the second harmonic (blue) laser pulse [28]. It is expressed as follows:

(5) E r ( t ) = E 0 1 + γ E f ( t ) [ cos ( ω r t ) x + sin ( ω r t ) y ] ,

(6) E b ( t ) = γ E E 0 1 + γ E f ( t ) [ cos ( ω b t ) x sin ( ω b t ) y ] .

With the maximum combined electric field amplitude E 0 , γ E is equal to 1, which means that the laser intensities of two CP laser fields are the same. ω r is the fundamental frequency (corresponding to the laser wavelength of the fundamental laser field), and ω b is the second harmonic frequency (corresponding to the laser wavelength of the second harmonic laser field). f ( t ) = sin 2 ( π t / N T ) is the pulse envelope function. T is an optical cycle (o.c.), and the number of optical cycles is four optical cycles ( N = 4 ). The motion of two electrons in the laser field needs to solve Hamilton’s equations, which is written as follows:

(7) d r i d t = H p i , d p i d t = H r i .

The aforementioned equations can be solved by the standard fourth to fifth Runge–Kutta algorithm. In order to ensure that the total energy is positive, the initial positions and momentum distributions of the two electrons satisfy the Gaussian distribution. Then, let the electrons evolve freely for a long enough time (100 a.u.) to obtain the initial ensemble when the electron position and momentum distribution are stable after getting a stable initial ensemble through screening. Then, the laser field is added, and all electrons begin to evolve under the action of the Coulomb potential and electric field until the pulse finishes. If the total energy of the two electrons is greater than zero at the end of the pulse, a DI event is recorded [29].

3 Results and discussion

In order to obtain the maximum recollision, we investigate the DI probability of CO2 molecules as a function of the laser intensity in the CRTC laser field with different laser wavelengths, as shown in Figure 1. We can find the “knee” structure, which means a significant electron–electron correlation effect in NSDI. The DI probability is changed with different laser wavelengths under a certain laser intensity, which indicates that the laser wavelength can affect the DI probability of CO2 molecules in the few-cycle CRTC laser field. However, it is a little bit different for linearly polarized and EP laser fields [21,40] because their curves of DI probability versus laser intensity twist together. Meanwhile, we choose the laser field intensity of I 0 = 2 × 10 14 W / cm 2 .

Figure 1 
               The DI probability of CO2 molecules as a function of the laser intensity in the CRTC laser field. The curves of red circle, green triangle, and blue square represent laser fields at laser wavelengths of 800 nm/400 nm, 1,200 nm/600 nm, and 1,600 nm/800 nm, respectively.
Figure 1

The DI probability of CO2 molecules as a function of the laser intensity in the CRTC laser field. The curves of red circle, green triangle, and blue square represent laser fields at laser wavelengths of 800 nm/400 nm, 1,200 nm/600 nm, and 1,600 nm/800 nm, respectively.

One of the important significances of NSDI is to explore the electron–electron correlation. Since there is no simple and clear polarized direction for the CRTC laser field, we explore the correlation between two electrons through recoil momentum distributions of the ion [18]. Because the net momentum of the electroneutral ion-electron system is zero, the relationship P CO 2 2+ = ( P e 1 + P e 2 ) can be obtained [29]. Therefore, Figure 2 shows the ion momentum distribution of CO2 molecules under different laser wavelengths in CRTC laser fields at the laser intensity 2 × 10 14 W / cm 2 . We can find that the shape of the negative vector potential of the electric field is a triangle structure. The shape of the ion momentum distribution is consistent with the shape of the negative vector potential of the electric field. Similar structures have been extensively studied both experimentally and theoretically [19,22]. In these studies, the results and analysis show that the RII is an important reason for the structure. Therefore, the RII forms ion momentum distributions far from the origin due to the short delay time and strong electron–electron correlation. The RESI forms ion momentum distributions near the origin due to the long delay time and weak electron–electron correlation. Figure 2 also shows that the ion momentum distribution gradually spreads out from the origin as the laser wavelength increases, which allows us to infer that the ionization channel has been changed. We will continue to analyze the details and separate the ionization channels by the delay time.

Figure 2 
               Ion momentum distributions in the CRTC laser field. The laser intensity is 
                     
                        
                        
                           2
                           ×
                           
                              
                                 10
                              
                              
                                 14
                              
                           
                           
                           
                              W
                              /
                              
                                 
                                    
                                       cm
                                    
                                    2
                                 
                              
                           
                        
                        2\times {10}^{14}\hspace{.25em}\text{W}/{\text{cm}}^{\text{2}}
                     
                  , and the laser wavelengths are chosen as (a) 800 nm/400 nm, (b) 1,200 nm/600 nm, and (c) 1,600 nm/800 nm, respectively. The solid red line represents the negative vector potential 
                     
                        
                        
                           ‒
                           A
                           (
                           t
                           )
                        
                        ‒{\boldsymbol{A}}(t)
                     
                   of the CRTC laser field.
Figure 2

Ion momentum distributions in the CRTC laser field. The laser intensity is 2 × 10 14 W / cm 2 , and the laser wavelengths are chosen as (a) 800 nm/400 nm, (b) 1,200 nm/600 nm, and (c) 1,600 nm/800 nm, respectively. The solid red line represents the negative vector potential A ( t ) of the CRTC laser field.

We investigate the DI time (t DI) and the recollision time (t RC) in the CRTC laser field with different laser wavelengths. When the recollision of the two electrons occurs, the Coulomb repulsion energy increases instantaneously [21]. We define t RC as the moment when the Coulomb potential energy between two electrons reaches its maximum [18]. Since the two electrons are indistinguishable, we define the final ionized electron as the second ionized electron after the recollision. Therefore, the corresponding t DI is defined as the second ionization moment after the recollision. The delay time (t DIt RC) is defined as the time interval between the recollision and the final DI [21,40].

Figure 3 shows that t DI versus t RC for different laser wavelengths in CRTC laser fields. The white diagonals of the pictures represent the moments when t DI and t RC are equal, and regions close to the white diagonal are marked by red dashed boxes. Therefore, the part far from the white diagonal means that the delay time is longer, corresponding to the RESI channel, whereas the part near the white diagonal means that the delay time is shorter, corresponding to the RII channel. As shown in Figure 3(a), when the laser wavelength of the combined laser field is 800 and 400 nm, the whole cluster is far away from the diagonal, and the proportion of the time distribution in the red dashed box is about 17.15%. As shown in Figure 3(b), when the laser wavelength of the combined laser field is 1,200 and 600 nm, the whole cluster is close to the diagonal, and the proportion of the time distribution in the red dashed box rises to 24.22%, which means that the delay time and the ionization channel have been changed. In Figure 3(c), we can also see that more of the whole cluster is close to the diagonal when the laser wavelength of the combined laser field is 1,600 and 800 nm, and the proportion of the time distribution in the red dashed box rises to 33.46%. From Figure 3, we can conclude that the laser wavelength will affect t DI and t RC. The underlying mechanism of t DI versus t RC means more RII channels and fewer RESI channels.

Figure 3 
               The DI time (t
                  DI) versus the recollision time (t
                  RC) for different laser wavelengths in the CRTC laser field. The laser intensity is 
                     
                        
                        
                           2
                           ×
                           
                              
                                 10
                              
                              
                                 14
                              
                           
                           
                           
                              W
                              /
                              
                                 
                                    
                                       cm
                                    
                                    2
                                 
                              
                           
                        
                        2\times {10}^{14}\hspace{.25em}\text{W}/{\text{cm}}^{\text{2}}
                     
                  , and the laser wavelengths are chosen as (a) 800 nm/400 nm, (b) 1,200 nm/600 nm, and (c) 1,600 nm/800 nm, respectively.
Figure 3

The DI time (t DI) versus the recollision time (t RC) for different laser wavelengths in the CRTC laser field. The laser intensity is 2 × 10 14 W / cm 2 , and the laser wavelengths are chosen as (a) 800 nm/400 nm, (b) 1,200 nm/600 nm, and (c) 1,600 nm/800 nm, respectively.

Figure 4 shows the counts of the delay-time distribution with different laser wavelengths in CRTC laser fields. We can find that the counts of delay time reach their first minimum around 0.12o.c–0.2o.c. with the increase of the laser wavelength, which clearly separates shorter and longer delay time [21,40]. Therefore, we define the RESI channel with a delay time between recollision and final DI longer than 0.16o.c. The others with a shorter delay time are classified as the RII channel.

Figure 4 
               Counts of the delay time for different laser wavelengths in the CRTC laser field. The laser intensity is 
                     
                        
                        
                           2
                           ×
                           
                              
                                 10
                              
                              
                                 14
                              
                           
                           
                              W
                              /
                              
                                 
                                    
                                       cm
                                    
                                    2
                                 
                              
                           
                        
                        2\times {10}^{14}\text{W}/{\text{cm}}^{\text{2}}
                     
                  , and the laser wavelengths are chosen as (a) 800 nm/400 nm, (b) 1,200 nm/600 nm, and (c) 1,600 nm/800 nm, respectively. The RII and RESI channels can be distinguished by the first minimum of the statistic results.
Figure 4

Counts of the delay time for different laser wavelengths in the CRTC laser field. The laser intensity is 2 × 10 14 W / cm 2 , and the laser wavelengths are chosen as (a) 800 nm/400 nm, (b) 1,200 nm/600 nm, and (c) 1,600 nm/800 nm, respectively. The RII and RESI channels can be distinguished by the first minimum of the statistic results.

The ion momentum distributions shown in Figure 2 can be separated into two ion momentum distributions (RESI and RII). Figure 5 shows the ion momentum distribution of RESI (RII) with the increase of the laser wavelength in the CRTC laser field at the laser intensity of 2 × 10 14 W / cm 2 . In Figure 2, the ion momentum distributions consist of two parts: the ion momentum distributions of the RESI channel and the ion momentum distributions of the RII channel. We simulate the above two parts separately and then obtain the ion momentum distributions ionized by the RESI channel, as shown in Figure 5(a)–(c), and the ion momentum distributions ionized by the RII channel, as shown in Figure 5(d)–(f). As shown in Figure 5(a)–(c), although the ion momentum distributions of RESI gradually diffuse to the surroundings with the increase of the laser wavelength, they are still mainly distributed in the central area of the triangular structure. This is because when the first electron is ionized, the second electron is not ionized immediately. Therefore, it means a weak electron–electron correlation in NSDI. In Figure 5(d)–(f), we can see that the ion momentum distributions of RII are mainly distributed in the edge region of the triangular structure with the increase of the laser wavelength. In this case, two electrons are ionized from the core at almost the same time. And it means a strong electron–electron correlation in NSDI.

Figure 5 
               Ion momentum distributions for different laser wavelengths in the CRTC laser field. The laser intensity is 
                     
                        
                        
                           2
                           ×
                           
                              
                                 10
                              
                              
                                 14
                              
                           
                           
                              W
                              /
                              
                                 
                                    
                                       cm
                                    
                                    2
                                 
                              
                           
                        
                        2\times {10}^{14}\text{W}/{\text{cm}}^{\text{2}}
                     
                  , and the laser wavelengths are chosen as (a) and (d) 800 nm/400 nm, (b) and (e) 1,200 nm/600 nm, and (c) and (f) 1,600 nm/800 nm, respectively. The upper row shows the RESI channel, and the lower row shows the RII channel.
Figure 5

Ion momentum distributions for different laser wavelengths in the CRTC laser field. The laser intensity is 2 × 10 14 W / cm 2 , and the laser wavelengths are chosen as (a) and (d) 800 nm/400 nm, (b) and (e) 1,200 nm/600 nm, and (c) and (f) 1,600 nm/800 nm, respectively. The upper row shows the RESI channel, and the lower row shows the RII channel.

To distinguish the RESI channel and RII channel clearly, Figure 6 shows the typical energy trajectories in the combined laser field of 800 and 400 nm. Among them, the black solid line is the repulsion energy of two electrons, and the red dashed line and the blue dotted line are the energies of the two electrons, respectively. In Figure 6, the recollision occurs when the repulsion energy is maximum. From Figure 6(a), we can find that the energy of the returning electron is greater than zero, and the bound electron oscillates at the excited bound state after recollision. Then, the DI occurs at the maximum of the field (the energy of both electrons is greater than zero). There is a longer delay time in this process, corresponding to the RESI. From Figure 6(b), we can find that two electrons are ionized directly after recollision (the energy of both electrons is greater than zero), which means that the returning electron collides strongly with the bound electron and transfers enough energy to make the bound electron ionized. The delay time is shorter in this process, corresponding to the RII.

Figure 6 
               Two typical energy trajectories of (a) RESI channel and (b) RII channel. The laser wavelength is 800 nm/400 nm, and the laser intensity is 
                     
                        
                        
                           2
                           ×
                           
                              
                                 10
                              
                              
                                 14
                              
                           
                           
                              W
                              /
                              
                                 
                                    
                                       cm
                                    
                                    2
                                 
                              
                           
                        
                        2\times {10}^{14}\text{W}/{\text{cm}}^{\text{2}}
                     
                  . The black solid line marks the repulsion energy between the two electrons. The red dashed line and blue dotted line mark the energies of the two electrons, respectively.
Figure 6

Two typical energy trajectories of (a) RESI channel and (b) RII channel. The laser wavelength is 800 nm/400 nm, and the laser intensity is 2 × 10 14 W / cm 2 . The black solid line marks the repulsion energy between the two electrons. The red dashed line and blue dotted line mark the energies of the two electrons, respectively.

In addition, the delay time is different in the two ionization channels, so there must be different trajectories. Figure 7 shows two typical classical trajectories in the combined laser field of 800 and 400 nm. Among them, the red solid line represents the trajectory of the returning electron, and the blue solid line represents the trajectory of the bound electron. In Figure 7(a), we can see that the returning electron has a shorter trajectory in the CRTC laser field, and then it collides with the bound electron. At this time, the bound electron oscillates at the excited bound state. The bound electrons ionize in the opposite direction to that of the returning electrons after a period of time. It corresponds to the RESI in Figure 6(a). However, in Figure 7(b), we can see that the returning electron has a longer trajectory in the CRTC laser field, then the returning electron collides with the bound electron. Due to the shorter delay time, both electrons are ionized in the same direction. It corresponds to the RII in Figure 6(b).

Figure 7 
               The trajectory of the two electrons in the 2D plane for (a) the RESI channel and (b) the RII channel. The laser wavelength is 800 nm/400 nm, and the laser intensity is 
                     
                        
                        
                           2
                           ×
                           
                              
                                 10
                              
                              
                                 14
                              
                           
                           
                              W
                              /
                              
                                 
                                    
                                       cm
                                    
                                    2
                                 
                              
                           
                        
                        2\times {10}^{14}\text{W}/{\text{cm}}^{\text{2}}
                     
                  . The red curve represents the trajectory of the returning electron, and the blue curve represents the trajectory of the bound electron.
Figure 7

The trajectory of the two electrons in the 2D plane for (a) the RESI channel and (b) the RII channel. The laser wavelength is 800 nm/400 nm, and the laser intensity is 2 × 10 14 W / cm 2 . The red curve represents the trajectory of the returning electron, and the blue curve represents the trajectory of the bound electron.

Finally, we discuss the ratio of the two ionization channels in the total ionization channels as the laser wavelength increases. We define the ratio of the RESI as N RESI / ( N RESI + N RII ) and the ratio of the RII as N RII / ( N RESI + N RII ) . N RESI and N RII are the counts of RESI and RII. In Figure 8, we can see that when the laser wavelength of the combined laser field is 400 and 200 nm, the ratio of the RESI channel is much higher than the ratio of the RII channel. We can also see that the two ionization channels show different trends with the change in the laser wavelength. As the laser wavelength increases, the ratio of the RESI channel decreases gradually while the ratio of the RII channel increases gradually. It means that the two channels can be controlled effectively by changing the laser wavelength. In addition, we find that although the RII channel increases with the increase of the laser wavelength, the dominant channel is still RESI.

Figure 8 
               Ratio between different ionization channels and total ionization versus laser wavelengths. The laser intensity is 
                     
                        
                        
                           2
                           ×
                           
                              
                                 10
                              
                              
                                 14
                              
                           
                           
                              W
                              /
                              
                                 
                                    
                                       cm
                                    
                                    2
                                 
                              
                           
                        
                        2\times {10}^{14}\text{W}/{\text{cm}}^{\text{2}}
                     
                  , and the dominance of the RESI is obvious.
Figure 8

Ratio between different ionization channels and total ionization versus laser wavelengths. The laser intensity is 2 × 10 14 W / cm 2 , and the dominance of the RESI is obvious.

4 Conclusions

In conclusion, the two ionization channels in NSDI of CO2 molecules in the few-cycle CRTC laser field are investigated with the 2D classical ensemble model. The numerical results show that the shape of the ion momentum distribution is consistent with the shape of the negative vector potential of the CRTC laser field, and the ion momentum distribution gradually spreads out from the origin as the laser wavelength increases. In addition, the delay time decreases gradually with the increase of the laser wavelength, but the electron–electron correlation is stronger. To further investigate the effect of the laser wavelength, we separate the ion momentum distribution of the RII channel and the RESI channel by the delay time. Finally, the ratio of the two ionization channels can be influenced by the laser wavelength, while the dominance of the RESI channel is obvious. This article further complemented the research of the NSDI electron dynamics process of linear triatomic molecules in CRTC laser fields, and the results also provide references to explore the electron–electron correlation in experiments.

  1. Funding information: The authors state no funding involved.

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

References

[1] Paulus GG, Nicklich W, Xu H, Lambropoulos P, Walther H. Plateau in above threshold ionization spectra. Phys Rev Lett. 1994;72(18):2851–4.10.1103/PhysRevLett.72.2851Search in Google Scholar PubMed

[2] Agostini P, Fabre F, Mainfray G, Petite G, Rahman NK. Free-free transitions following six-photon ionization of xenon atoms. Phys Rev Lett. 1979;42(17):1127–30.10.1103/PhysRevLett.42.1127Search in Google Scholar

[3] Kruit P, Kimman J, Muller HG, van der Wiel MJ. Electron spectra from multiphoton ionization of xenon at 1064, 532, and 355 nm. Phys Rev A. 1983;28(1):248–55.10.1103/PhysRevA.28.248Search in Google Scholar

[4] Song KL, Yu WW, Ben S, Xu TT, Zhang HD, Guo PY, et al. Theoretical study on non-sequential double ionization of carbon disulfide with different bond lengths in linearly polarized laser fields. Chin Phys B. 2017;26(2):023204.10.1088/1674-1056/26/2/023204Search in Google Scholar

[5] Walker B, Sheehy B, DiMauro LF, Agostini P, Schafer KJ, Klander KC. Precision measurement of strong field double ionization of helium. Phys Rev Lett. 1994;73(9):1227–30.10.1103/PhysRevLett.73.1227Search in Google Scholar PubMed

[6] Fittinghoff DN, Bolton PR, Chang B, Kulander KC. Observation of nonsequential double ionization of helium with optical tunneling. Phys Rev Lett. 1992;69(18):2642–5.10.1103/PhysRevLett.69.2642Search in Google Scholar PubMed

[7] Krausz F, Ivanov M. Attosecond physics. Rev Mod Phys. 2009;81(1):163–234.10.1103/RevModPhys.81.163Search in Google Scholar

[8] Li HY, Wang B, Chen J, Jiang HB, Li XF, Liu J, et al. Effects of a static electric field on nonsequential double ionization. Phys Rev A. 2007;76(3):033405.10.1103/PhysRevA.76.033405Search in Google Scholar

[9] Becker W, Liu X, Ho PJ, Eberly JH. Theories of photoelectron correlation in laser-driven multiple atomic ionization. Rev Mod Phys. 2012;84(3):1011–43.10.1103/RevModPhys.84.1011Search in Google Scholar

[10] Fleischer A, Wörner HJ, Arissian L, Liu LR, Meckel M, Pippert A, et al. Probing angular correlations in sequential double ionization. Phys Rev Lett. 2011;107(11):113003.10.1103/PhysRevLett.107.113003Search in Google Scholar PubMed

[11] Wang X, Tian J, Eberly JH. Angular correlation in strong-field double ionization under circular polarization. Phys Rev Lett. 2013;110(7):073001.10.1103/PhysRevLett.110.073001Search in Google Scholar PubMed

[12] Guo J, Wang T, Liu XS, Sun JZ. Non-sequential double ionization of Mg atoms in elliptically polarized laser fields. Laser Physics. 2013;23(5):055303.10.1088/1054-660X/23/5/055303Search in Google Scholar

[13] Mauger F, Chandre C, Uzer T. Recollisions and correlated double ionization with circularly polarized light. Phys Rev Lett. 2010;105(8):083002.10.1103/PhysRevLett.105.083002Search in Google Scholar PubMed

[14] Fu LB, Xin GG, Ye DF, Liu J. Recollision dynamics and phase diagram for nonsequential double ionization with circularly polarized laser fields. Phys Rev Lett. 2012;108(10):103601.10.1103/PhysRevLett.108.103601Search in Google Scholar PubMed

[15] Hao XL, Wang GQ, Jia XY, Li WD. Nonsequential double ionization of Ne in an elliptically polarized intense laser field. Phys Rev A. 2009;80(2):023408.10.1103/PhysRevA.80.023408Search in Google Scholar

[16] Corkum PB. Plasma perspective on strong field multiphoton ionization. Phys Rev Lett. 1993;71(13):1994–7.10.1103/PhysRevLett.71.1994Search in Google Scholar PubMed

[17] Kulander KC, Cooper J, Schafer KJ. Laser-assisted inelastic rescattering during above-threshold ionization. Phys Rev A. 1995;51(1):561–8.10.1103/PhysRevA.51.561Search in Google Scholar PubMed

[18] Chen JH, Xu TT, Han T, Sun Y, Xu QY, Liu XS. Relative phase effect of nonsequential double ionization in Ar by two-color elliptically polarized laser field. Chin Phys B. 2020;29(1):013203.10.1088/1674-1056/ab5d03Search in Google Scholar

[19] Huang C, Zhou Y, Zhang Q, Lu P. Contribution of recollision ionization to the cross-shaped structure in nonsequential double ionization. Opt Express. 2013;21(9):11382.10.1364/OE.21.011382Search in Google Scholar PubMed

[20] Ye DF, Liu J. Strong-field double ionization at the transition to below the recollision threshold. Phys Rev A. 2010;81(4):043402.10.1103/PhysRevA.81.043402Search in Google Scholar

[21] Ben S, Wang T, Xu TT, Guo J, Liu XS. Nonsequential double ionization channels control of Ar with few-cycle elliptically polarized laser pulse by carrier-envelope-phase. Opt Express. 2016;24(7):7525.10.1364/OE.24.007525Search in Google Scholar PubMed

[22] Bergues B, Kübel M, Johnson NG, Fischer B, Camus N, Betsch KJ, et al. Attosecond tracing of correlated electron-emission in non-sequential double ionization. Nat Commun. 2012;3(1):813.10.1038/ncomms1807Search in Google Scholar PubMed

[23] Maxwell AS, Faria CFdeM. Controlling below-threshold nonsequential double ionization via quantum interference. Phys Rev Lett. 2016;116(14):143001.10.1103/PhysRevLett.116.143001Search in Google Scholar PubMed

[24] Ye D, Li M, Fu L, Liu J, Gong QH, Liu YQ, et al. Scaling laws of the two-electron sum-energy spectrum in strong-field double ionization. Phys Rev Lett. 2015;115(12):123001.10.1103/PhysRevLett.115.123001Search in Google Scholar PubMed

[25] Hao JX, Hao XL, Li WD, Hu SL, Chen J. Controlling three-dimensional electron–electron correlation via elliptically polarized intense laser field. Chin Phys Lett. 2017;34(4):043201.10.1088/0256-307X/34/4/043201Search in Google Scholar

[26] Huang C, Guo W, Zhou Y, Wu Z. Role of coulomb repulsion in correlated-electron emission from a doubly excited state in nonsequential double ionization of molecules. Phys Rev A. 2016;93(1):013416.10.1103/PhysRevA.93.013416Search in Google Scholar

[27] Liu YQ, Fu LB, Ye DF, Liu J, Li M, Wu CY, et al. Strong-field double ionization through sequential release from double excitation with subsequent coulomb scattering. Phys Rev Lett. 2014;112(1):013003.10.1103/PhysRevLett.112.013003Search in Google Scholar PubMed

[28] Xu TT, Chen JH, Pan XF, Zhang HD, Ben S, Liu XS. Effect of elliptical polarizations on nonsequential double ionization in two-color elliptically polarized laser fields. Chin Phys B. 2018;27(9):093201.10.1088/1674-1056/27/9/093201Search in Google Scholar

[29] Ben S, Guo PY, Pan XF, Xu TT, Song KL, Liu XS. Recollision induced excitation-ionization with counter-rotating two-color circularly polarized laser field. Chem Phys Lett. 2017;679:38–44.10.1016/j.cplett.2017.04.076Search in Google Scholar

[30] Eckart S, Richter M, Kunitski M, Hartung A, Rist J, Henrichs K, et al. Nonsequential double ionization by counterrotating circularly polarized two-color laser fields. Phys Rev Lett. 2016;117(13):133202.10.1103/PhysRevLett.117.133202Search in Google Scholar PubMed

[31] Mancuso CA, Dorney KM, Dorney DD, Chaloupka JL, Ellis JL, Dollar FJ, et al. Controlling nonsequential double ionization in two-color circularly polarized femtosecond laser fields. Phys Rev Lett. 2016;117(13):133201.10.1103/PhysRevLett.117.133201Search in Google Scholar PubMed

[32] Chaloupka JL, Hickstein DD. Dynamics of strong-field double ionization in two-color counterrotating fields. Phys Rev Lett. 2016;116(14):143005.10.1103/PhysRevLett.116.143005Search in Google Scholar PubMed

[33] Li BQ, Yang X, Ren XH, Zhang JT. Enhanced double ionization rate from O2 molecules driven by counter-rotating circularly polarized two-color laser fields. Opt Express. 2019;27(22):32700–08.10.1364/OE.27.032700Search in Google Scholar PubMed

[34] Peng M, Bai LH, Guo Z. Influence of relative phase on nonsequential double ionization process of CO2 molecules by counter-rotating two-color circularly polarized laser fields. CommunTheorPhys. 2021;73(7):075501.10.1088/1572-9494/abfb60Search in Google Scholar

[35] Su J, Liu ZC, Liao JY, Huang XF, Li YB, Huang C. Electron correlation and recollision dynamics in nonsequential double ionization by counter-rotating two-color elliptically polarized laser fields. Opt Express. 2022;30(14):24898.10.1364/OE.462022Search in Google Scholar PubMed

[36] Panfili R, Eberly JH, Haan SL. Comparing classical and quantum dynamics of strong-field double ionization. Opt Express. 2001;8(7):431.10.1364/OE.8.000431Search in Google Scholar

[37] Haan SL, Breen L, Karim A, Eberly JH. Variable time lag and backward ejection in full-dimensional analysis of strong-field double ionization. Phys Rev Lett. 2006;97(10):103008.10.1103/PhysRevLett.97.103008Search in Google Scholar PubMed

[38] Xu TT, Ben S, Wang T, Zhang J, Guo J, Liu XS. Exploration of the nonsequential double-ionization process of a Mg atom with different delay time in few-cycle circularly polarized laser fields. Phys Rev A. 2015;92(3):033405.10.1103/PhysRevA.92.033405Search in Google Scholar

[39] Huang C, Zhong M, Wu Z. Intensity-dependent two-electron emission dynamics in nonsequential double ionization by counter-rotating two-color circularly polarized laser fields. Opt Express. 2018;26(20):26045.10.1364/OE.26.026045Search in Google Scholar PubMed

[40] Zhu QY, Xu TT, Ben S, Chen JH, Song KL, Liu XS The carrier-envelope phase effect in non-sequential double ionization by few-cycle linearly polarized laser field. Opt Commun 2018;426:602–6.10.1016/j.optcom.2018.06.022Search in Google Scholar

Received: 2023-07-16
Revised: 2023-08-29
Accepted: 2023-09-04
Published Online: 2023-10-07

© 2023 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Dynamic properties of the attachment oscillator arising in the nanophysics
  3. Parametric simulation of stagnation point flow of motile microorganism hybrid nanofluid across a circular cylinder with sinusoidal radius
  4. Fractal-fractional advection–diffusion–reaction equations by Ritz approximation approach
  5. Behaviour and onset of low-dimensional chaos with a periodically varying loss in single-mode homogeneously broadened laser
  6. Ammonia gas-sensing behavior of uniform nanostructured PPy film prepared by simple-straightforward in situ chemical vapor oxidation
  7. Analysis of the working mechanism and detection sensitivity of a flash detector
  8. Flat and bent branes with inner structure in two-field mimetic gravity
  9. Heat transfer analysis of the MHD stagnation-point flow of third-grade fluid over a porous sheet with thermal radiation effect: An algorithmic approach
  10. Weighted survival functional entropy and its properties
  11. Bioconvection effect in the Carreau nanofluid with Cattaneo–Christov heat flux using stagnation point flow in the entropy generation: Micromachines level study
  12. Study on the impulse mechanism of optical films formed by laser plasma shock waves
  13. Analysis of sweeping jet and film composite cooling using the decoupled model
  14. Research on the influence of trapezoidal magnetization of bonded magnetic ring on cogging torque
  15. Tripartite entanglement and entanglement transfer in a hybrid cavity magnomechanical system
  16. Compounded Bell-G class of statistical models with applications to COVID-19 and actuarial data
  17. Degradation of Vibrio cholerae from drinking water by the underwater capillary discharge
  18. Multiple Lie symmetry solutions for effects of viscous on magnetohydrodynamic flow and heat transfer in non-Newtonian thin film
  19. Thermal characterization of heat source (sink) on hybridized (Cu–Ag/EG) nanofluid flow via solid stretchable sheet
  20. Optimizing condition monitoring of ball bearings: An integrated approach using decision tree and extreme learning machine for effective decision-making
  21. Study on the inter-porosity transfer rate and producing degree of matrix in fractured-porous gas reservoirs
  22. Interstellar radiation as a Maxwell field: Improved numerical scheme and application to the spectral energy density
  23. Numerical study of hybridized Williamson nanofluid flow with TC4 and Nichrome over an extending surface
  24. Controlling the physical field using the shape function technique
  25. Significance of heat and mass transport in peristaltic flow of Jeffrey material subject to chemical reaction and radiation phenomenon through a tapered channel
  26. Complex dynamics of a sub-quadratic Lorenz-like system
  27. Stability control in a helicoidal spin–orbit-coupled open Bose–Bose mixture
  28. Research on WPD and DBSCAN-L-ISOMAP for circuit fault feature extraction
  29. Simulation for formation process of atomic orbitals by the finite difference time domain method based on the eight-element Dirac equation
  30. A modified power-law model: Properties, estimation, and applications
  31. Bayesian and non-Bayesian estimation of dynamic cumulative residual Tsallis entropy for moment exponential distribution under progressive censored type II
  32. Computational analysis and biomechanical study of Oldroyd-B fluid with homogeneous and heterogeneous reactions through a vertical non-uniform channel
  33. Predictability of machine learning framework in cross-section data
  34. Chaotic characteristics and mixing performance of pseudoplastic fluids in a stirred tank
  35. Isomorphic shut form valuation for quantum field theory and biological population models
  36. Vibration sensitivity minimization of an ultra-stable optical reference cavity based on orthogonal experimental design
  37. Effect of dysprosium on the radiation-shielding features of SiO2–PbO–B2O3 glasses
  38. Asymptotic formulations of anti-plane problems in pre-stressed compressible elastic laminates
  39. A study on soliton, lump solutions to a generalized (3+1)-dimensional Hirota--Satsuma--Ito equation
  40. Tangential electrostatic field at metal surfaces
  41. Bioconvective gyrotactic microorganisms in third-grade nanofluid flow over a Riga surface with stratification: An approach to entropy minimization
  42. Infrared spectroscopy for ageing assessment of insulating oils via dielectric loss factor and interfacial tension
  43. Influence of cationic surfactants on the growth of gypsum crystals
  44. Study on instability mechanism of KCl/PHPA drilling waste fluid
  45. Analytical solutions of the extended Kadomtsev–Petviashvili equation in nonlinear media
  46. A novel compact highly sensitive non-invasive microwave antenna sensor for blood glucose monitoring
  47. Inspection of Couette and pressure-driven Poiseuille entropy-optimized dissipated flow in a suction/injection horizontal channel: Analytical solutions
  48. Conserved vectors and solutions of the two-dimensional potential KP equation
  49. The reciprocal linear effect, a new optical effect of the Sagnac type
  50. Optimal interatomic potentials using modified method of least squares: Optimal form of interatomic potentials
  51. The soliton solutions for stochastic Calogero–Bogoyavlenskii Schiff equation in plasma physics/fluid mechanics
  52. Research on absolute ranging technology of resampling phase comparison method based on FMCW
  53. Analysis of Cu and Zn contents in aluminum alloys by femtosecond laser-ablation spark-induced breakdown spectroscopy
  54. Nonsequential double ionization channels control of CO2 molecules with counter-rotating two-color circularly polarized laser field by laser wavelength
  55. Fractional-order modeling: Analysis of foam drainage and Fisher's equations
  56. Thermo-solutal Marangoni convective Darcy-Forchheimer bio-hybrid nanofluid flow over a permeable disk with activation energy: Analysis of interfacial nanolayer thickness
  57. Investigation on topology-optimized compressor piston by metal additive manufacturing technique: Analytical and numeric computational modeling using finite element analysis in ANSYS
  58. Breast cancer segmentation using a hybrid AttendSeg architecture combined with a gravitational clustering optimization algorithm using mathematical modelling
  59. On the localized and periodic solutions to the time-fractional Klein-Gordan equations: Optimal additive function method and new iterative method
  60. 3D thin-film nanofluid flow with heat transfer on an inclined disc by using HWCM
  61. Numerical study of static pressure on the sonochemistry characteristics of the gas bubble under acoustic excitation
  62. Optimal auxiliary function method for analyzing nonlinear system of coupled Schrödinger–KdV equation with Caputo operator
  63. Analysis of magnetized micropolar fluid subjected to generalized heat-mass transfer theories
  64. Does the Mott problem extend to Geiger counters?
  65. Stability analysis, phase plane analysis, and isolated soliton solution to the LGH equation in mathematical physics
  66. Effects of Joule heating and reaction mechanisms on couple stress fluid flow with peristalsis in the presence of a porous material through an inclined channel
  67. Bayesian and E-Bayesian estimation based on constant-stress partially accelerated life testing for inverted Topp–Leone distribution
  68. Dynamical and physical characteristics of soliton solutions to the (2+1)-dimensional Konopelchenko–Dubrovsky system
  69. Study of fractional variable order COVID-19 environmental transformation model
  70. Sisko nanofluid flow through exponential stretching sheet with swimming of motile gyrotactic microorganisms: An application to nanoengineering
  71. Influence of the regularization scheme in the QCD phase diagram in the PNJL model
  72. Fixed-point theory and numerical analysis of an epidemic model with fractional calculus: Exploring dynamical behavior
  73. Computational analysis of reconstructing current and sag of three-phase overhead line based on the TMR sensor array
  74. Investigation of tripled sine-Gordon equation: Localized modes in multi-stacked long Josephson junctions
  75. High-sensitivity on-chip temperature sensor based on cascaded microring resonators
  76. Pathological study on uncertain numbers and proposed solutions for discrete fuzzy fractional order calculus
  77. Bifurcation, chaotic behavior, and traveling wave solution of stochastic coupled Konno–Oono equation with multiplicative noise in the Stratonovich sense
  78. Thermal radiation and heat generation on three-dimensional Casson fluid motion via porous stretching surface with variable thermal conductivity
  79. Numerical simulation and analysis of Airy's-type equation
  80. A homotopy perturbation method with Elzaki transformation for solving the fractional Biswas–Milovic model
  81. Heat transfer performance of magnetohydrodynamic multiphase nanofluid flow of Cu–Al2O3/H2O over a stretching cylinder
  82. ΛCDM and the principle of equivalence
  83. Axisymmetric stagnation-point flow of non-Newtonian nanomaterial and heat transport over a lubricated surface: Hybrid homotopy analysis method simulations
  84. HAM simulation for bioconvective magnetohydrodynamic flow of Walters-B fluid containing nanoparticles and microorganisms past a stretching sheet with velocity slip and convective conditions
  85. Coupled heat and mass transfer mathematical study for lubricated non-Newtonian nanomaterial conveying oblique stagnation point flow: A comparison of viscous and viscoelastic nanofluid model
  86. Power Topp–Leone exponential negative family of distributions with numerical illustrations to engineering and biological data
  87. Extracting solitary solutions of the nonlinear Kaup–Kupershmidt (KK) equation by analytical method
  88. A case study on the environmental and economic impact of photovoltaic systems in wastewater treatment plants
  89. Application of IoT network for marine wildlife surveillance
  90. Non-similar modeling and numerical simulations of microploar hybrid nanofluid adjacent to isothermal sphere
  91. Joint optimization of two-dimensional warranty period and maintenance strategy considering availability and cost constraints
  92. Numerical investigation of the flow characteristics involving dissipation and slip effects in a convectively nanofluid within a porous medium
  93. Spectral uncertainty analysis of grassland and its camouflage materials based on land-based hyperspectral images
  94. Application of low-altitude wind shear recognition algorithm and laser wind radar in aviation meteorological services
  95. Investigation of different structures of screw extruders on the flow in direct ink writing SiC slurry based on LBM
  96. Harmonic current suppression method of virtual DC motor based on fuzzy sliding mode
  97. Micropolar flow and heat transfer within a permeable channel using the successive linearization method
  98. Different lump k-soliton solutions to (2+1)-dimensional KdV system using Hirota binary Bell polynomials
  99. Investigation of nanomaterials in flow of non-Newtonian liquid toward a stretchable surface
  100. Weak beat frequency extraction method for photon Doppler signal with low signal-to-noise ratio
  101. Electrokinetic energy conversion of nanofluids in porous microtubes with Green’s function
  102. Examining the role of activation energy and convective boundary conditions in nanofluid behavior of Couette-Poiseuille flow
  103. Review Article
  104. Effects of stretching on phase transformation of PVDF and its copolymers: A review
  105. Special Issue on Transport phenomena and thermal analysis in micro/nano-scale structure surfaces - Part IV
  106. Prediction and monitoring model for farmland environmental system using soil sensor and neural network algorithm
  107. Special Issue on Advanced Topics on the Modelling and Assessment of Complicated Physical Phenomena - Part III
  108. Some standard and nonstandard finite difference schemes for a reaction–diffusion–chemotaxis model
  109. Special Issue on Advanced Energy Materials - Part II
  110. Rapid productivity prediction method for frac hits affected wells based on gas reservoir numerical simulation and probability method
  111. Special Issue on Novel Numerical and Analytical Techniques for Fractional Nonlinear Schrodinger Type - Part III
  112. Adomian decomposition method for solution of fourteenth order boundary value problems
  113. New soliton solutions of modified (3+1)-D Wazwaz–Benjamin–Bona–Mahony and (2+1)-D cubic Klein–Gordon equations using first integral method
  114. On traveling wave solutions to Manakov model with variable coefficients
  115. Rational approximation for solving Fredholm integro-differential equations by new algorithm
  116. Special Issue on Predicting pattern alterations in nature - Part I
  117. Modeling the monkeypox infection using the Mittag–Leffler kernel
  118. Spectral analysis of variable-order multi-terms fractional differential equations
  119. Special Issue on Nanomaterial utilization and structural optimization - Part I
  120. Heat treatment and tensile test of 3D-printed parts manufactured at different build orientations
Downloaded on 14.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2023-0114/html
Scroll to top button