Asymmetric structure of atomic above-threshold ionization spectrum in two-color elliptically polarized laser fields∗

Xu-Cong Zhou(周旭聪),Shang Shi(石尚),Fei Li(李飞),Yu-Jun Yang(杨玉军),Jing Chen(陈京),Qing-Tian Meng(孟庆田),and Bing-Bing Wang(王兵兵),†

1School of Physics and Electronics,Shandong Normal University,Jinan 250358,China

2Laboratory of Optical Physics,Beijing National Laboratory for Condensed Matter Physics,Institute of Physics,Chinese Academy of Sciences,Beijing 100190,China

3Insititute of Atomic and Molecular Physics,Jilin University,Changchun 130012,China

4HEDPS,Center for Applied Physics and Technology,Peking University,Beijing 100871,China

5Institute of Applied Physics and Computational Mathematics,Beijing 100088,China

6University of Chinese Academy of Sciences,Beijing 100049,China

Keywords:above-threshold ionization,elliptically polarized,two color laser fields,asymmetric structure

1.Introduction

In the last decades,the rapid development of ultrashort laser pulse technology has promoted the study of light-matter interactions. Especially when an atom or molecule is exposed to an intense laser field,a lot of novel phenomena,such as above-threshold ionization(ATI),[1,2]high-order abovethreshold ionization(HATI),[3,5]high harmonic generation(HHG),[6,7]and non-sequential double ionization(NSDI),[8,9]may occur. These processes can be described by three-step model.[10]In this model,the bound electron in an atom or molecule can firstly absorb photons and be ionized by the strong laser.In the next step,the ionized electron is accelerated to escape from the parent ion.Finally,when the direction of laser field is reversed,part of ionized electrons can go back to the parent ion and collide with it.Some electrons can combine with the parent ion and then emit high energy photons,which process is called HHG;some may have an elastic collision with the parent ion,which may lead to the HATI;while part of electrons have inelastic collisions and cause another electron to be ionized,which is called NSDI.By using elliptical polarized laser fields,it was found that the probability and the cutoff of HATI spectrum decrease with the degree of the laser ellipticity increasing,which demonstrates the correction to the three-step model.[11]

Recently,the ionization of an atom in two-color laser fields is becoming a hotspot.Kazansky et al.[12]studied the angle-resolved spectra of photoelectrons by infrared and extreme ultraviolet(IR+XUV)two-color laser pulses.Radcliffe et al.[13]found a plateau-shape envelope of the ATI spectrum with increasing the intensity of the IR laser field. What is more,Zhang et al.[14]observed a step-like structure of ATI spectrum in IR+XUV two-color laser fields. Most of these studies are based on linearly polarized(LP)laser fields,because EP laser field is difficult to deal with experimentally and theoretically.Therefore,the ionization by an elliptically polarized(EP)laser pulse has attracted the attention of many scientists.[15,18]Busuladžić et al.[15]investigated the symmetry properties of molecular HATI spectra by changing the molecular orientation in an EP laser field.Li et al.[16]investigated the strong-field ionization dynamics of atoms in circularly polarized laser fields by a three-dimensional(3D)electron ensemble method.Furthermore,Yu et al.[17]explored the ionization dynamics of high-energy photoelectrons in EP laser fields.In this paper,we investigate the asymmetric structure for ATI of an atom in two-color EP laser fields.

2.Theoretical method

Here the ionization of the atom in two-color laser fields is taken as a research object.The laser field includes an EP laser field and an LP laser field.This atom–laser system can be regarded as an isolated system,where the total energy is conserved.Based on the frequency-domain theory,the Hamiltonian of this system is

where H0=(−i∇)2/2me+ω1Na1+ω2Na2represents the electron–photon energy operator,with Na1and Na2being the photon number operators of the two laser fields with different frequencies ω1and ω2,respectively,andwith i=1 and 2,where a and a†are the annihilation and creation operators of the photon mode respectively;U represents the atomic binding potential,and V is the electron–laser interaction potential and expressed as

We now consider the ATI of an atom in a two-color laser field as a genuine scattering process.Using the formal scattering method,[20]the transition matrix element can be written as

Here,the first term and the second term on the right-hand side of Eq.(3)correspond to the direct ATI and the rescattering ATI(i.e.,HATI),respectively.So equation(3)can also be written as Tfi=Td+Tr,where Tdis the process of the electrons being ionized by directly absorbing photons,while Tris the process of the electrons absorbing photons then going back to collide with the parent ion before it has been ionized.

The initial quantum state is expressed aswhich corresponds to the eigenstate of the Hamiltonian H0+U with the associated energy Ei=−Ip+(l1+1/2)ω1+(l2+1/2)ω2,in which Φi(r)is the wavefunction describing the ground state of an atom,is the initial state of the laser field for(j=1,2),and Iprepresents the threshold of the atom.The final state of the system can be denoted aswhich is the Volkov state of the electron in the two laser fields,which can be expressed as

where Veis the normalized volume,is the final state momentum of ionized electron,up1=Up1/ω1(up2=Up2/ω2),with Up1(Up2)being the ponderomotive potential in the first(second)laser field,andis the photon momentum of the first(second)laser field.The total energy of final state is

where m1and m2are photon numbers of the two laser fields,respectively.

The generalized Bessel functioncan be written as

On a whole,the direct ATI can be expressed as

where j1=l1−m1,j2=l2−m2,with j1and j2being the photon numbers of the electron absorbing from the two laser fields.And the transition matrix element of HATI is given by

In order to obtain Eq.(8),we utilize the completeness relation of the Volkov statesTherefore,the rescattering process can be divided into two parts,i.e.,the ATI process and the laser-assisted LAC process.[21]

with s1=l1−n1,s2=l2−n2.

Finally,we can obtain[22]

where q1=n1−m1,q2=n2−m2represent the number of the electrons absorbed from different fields in the LAC process. We assume that the initial phases of the two laser fields are both 0,which is φ1=φ2=0 for simplicity. The s wave function in the ground state of the momentum space iswith α=2meEB. In this work,we employ the short potential U=−exp(−r)/r to represent the interaction between the ionized electron and its parent ion.Atomic units are used throughout this work unless stated otherwise.

3.Numerical results

Firstly,we study the total angle-resolved ATI spectrum of hydrogen atom in two-color laser fields,which includes a beam of EP IR laser and a beam of LP XUV laser. The major axis of the EP laser field is along the z axis,and the minor axis is along the x axis. The LP laser field is along the z axis.The frequency of IR laser field is ω1=1.165 eV,and the field intensity is I1=1×1013W/cm2,while the frequency of the XUV laser is ω2=50ω1and the field intensity is I2=I1=1×1013W/cm2.

Fig.1.Total angle-resolved ATI spectrum of H atom with ε1=0(a),0.27(b),0.58(c),and 1(d),and other parameters being IP=13.6 eV,ω2=50ω1,ω1=1.165 eV,and I1=I2=1×1013 W/cm2.

Figure 1 shows the total angle-resolved ATI spectrum of hydrogen atom for different values of right-handed EP degree ε1of IR laser field,where θfis the angle between the momentum direction of the ionized electron and z axis.From Fig.1,we can clearly see that the spectrum shows two plateaus,in which the probability of the first plateau is much larger than that of the second plateau.Especially,the total ATI spectrum is symmetric about the angle of π/2 in Fig.1(a),when the IR laser is linearly polarized.However,the spectrum shows asymmetric for the IR laser becomes an EP field.Meanwhile,we can find that the spectrum shifts toward the right with the increase of the EP degree of the IR laser from Figs.1(b)–1(d).This phenomenon can be easily understood as follows. For the EP IR laser,the ionization can be regarded as a tunneling ionization process,where the largest ionization probability is along the direction of the IR laser polarization,i.e.,θf=ξ1/2;while for the XUV laser,the ionization can be regarded as a multiphoton ionization process,hence the ionization probability is not sensitive to the momentum direction of the ionized electron.As a result,the total angle-resolved ATI spectrum is mainly shifted by the EP laser field.Therefore,the greater the polarization degree of the EP laser field,the more the spectrum moves to the right.

By setting the spectral minimum point of the first plateau to be the reference point(marked by the arrows in Fig.1),we may illustrate the shift of the spectrum versus the degree of the elliptical polarization of IR laser in Fig.2,where the reference points are denoted by black square dots.Since the degree of the elliptical polarization is defined as ε1=tan(ξ1/2),we can obtain the shift angle of the spectrum(denoted by red circle dots)from the formula ∆ξ=arctanε1. Figure 2 shows that the reference point curve coincide well with the curve of ∆ξ versus ε1.Moreover,the spectrum shift is strongly dependent on the polarization degree of EP laser field,and the spectrum shifts about π/4 as the polarization degree of the IR laser field changes from ε1=0 to 1.These results indicate that we may measure the degree of the laser elliptical polarization by angleresolved ATI spectrum of atoms.

Moreover,figure 1 also shows the interference fringes on the two plateaus of the ATI spectrum.To explain these interference fringes on the spectrum,we analyze the total spectrum by two processes in the following:direct ATI and rescattering ATI.

Fig.2.Spectrum shifting with polarization degree of EP laser field increasing.

Figure 3 depicts the direct ATI spectrum with different EP degrees of IR laser.The shift of the spectrum with EP degree is consistent with what we discussed above. We now focus on finding the cause of the interference fringes on the spectra.Under our present laser conditions,the generalized Bessel function in Eq.(7)can be reduced into

Fig.3.Angle-resolved direct ATI spectrum of H atom with ε1=0(a),0.27(b),0.58(c),and 1(d).

According to the property of the Bessel function,the Bessel functioncan be zero when ξ2=0 for the orderwhere ξ2=0 corresponds to the case that the direction of emitted electrons is perpendicular to the XUV laser polarization(i.e.,θf=π/2).Therefore,the direct ATI spectra present the interference fringes for q2=1 and 2 at the angleMoreover,since the Bessel function J−q6(ζ6)provides a small modulation in Eq.(11),the position of the interference fringes spreads a width around θf=π/2.Figure 4 showsfor the two plateaus of the direct ATI spectra by Eq.(12)with the sum for q6=−2,−1,and 0.One may find that the spectra in Fig.4 are consistent well with those of Fig.3.

Fig.4.Values of by Eq.(12)as a function of θf with sum for q6=−2,−1,and 0.

On the other hand,the width of each plateau on the direct ATI spectrum is found to be determined by the Bessel function J−q1(ζ1,ζ3)under the present laser conditions.This Bessel function can be expressed as

with T=2π/ω1.Here,the IR laser filed can be regarded as a classical field,whereis the laser’s vector potential with E1being the amplitude of IR laser field andthe direction of the laser polarization.So the classical action of an electron in IR laser field is

where q2represents the number of XUV photons absorbed by the electron.Furthermore,we can obtain the final energy of ionized electron as follows:

Fig.5.Angle-resolved direct ATI spectrum of H atom with ε1=0.27,beginning value(dashed line),and cutoff value(solid line of classical orbits at tc=cos(ω1t0)(c=1,2).

By the saddle-point approximation,the beginning and the cutoff curves can be simulated by Eq.(17),and the results are shown in Fig.5,indicating that they are also consistent well with the angle-resolved direct ATI spectrum.

In the following part,we investigate the angle-resolved rescattering ATI spectrum with different EP degrees of IR laser and LP XUV laser as shown in Fig.6.

Obviously,the angle-resolved HATI spectrum of H atom also shows clearly dip structures.In order to study the cause of the dip structures on the spectrum in detail,we divide the spectrum of HATI for ε1=0.27 into two parts depending on the different XUV photons absorbed by the electron as shown in Fig.7.

Figure 7(a)shows the first plateau where the electron absorbs one XUV photon,while figure 7(b)presents the second plateau where the electron absorbs two XUV photons.Obviously,the dip structure exists mainly on the second plateau.

According to the frequency-domain theory,we can regard HATI as two processes:ATI and LAC shown by Eqs.(9)and(10),respectively. In the following,we investigate the second plateau especially.We now define the channel(s2,q2),in which s2and q2refer to the electron-absorbed number of XUV photons in the ATI process and the LAC process,respectively.According to the above definition,the second plateau consists of three channels,(0,2),(1,1),and(2,0),where channels(1,1)and(2,0)are shown in Fig.8.Channel(0,2)has little contribution to HATI,hence it may be ignored.

Fig.6.Angle-resolved HATI spectrum of H atom with ε1=0(a),0.27(b),0.28(c),and 1(d),other parameters are IP=13.6 eV,ω2=50ω1,ω1=1.165 eV,and I1=I2=1×1013 W/cm2.

Fig.7.Angle-resolved HATI spectrum for(a)the first plateau and(b)the second plateau.

From Fig.8,we find that the dip structure of the second plateau mainly comes from the contribution of channel(1,1).In order to explain these dip structures,we investigate further channel(1,1)for different values of angle θnbetween the direction of electron momentum Pnand the polarization direction of laser field,with θn=20◦,40◦,60◦,and 160◦as shown in Figs.9(a)–9(d).It is evident that all of the spectra show dip structures for four different angles,and the dip structure shifts rightwards with the increase of the angle.

Subsequently,we define sub-channel(s2|s1,q2),where s1is the number of IR photons absorbed by the electron in the ATI process.We show the spectra of different sub-channels in Fig.9 with sub-channel(1|−10,1)(Figs.9(a1)–9(d1)),(1|2,1)(Figs.9(a2)–9(d2)),and(1|10,1)(Figs.9(a3)–9(d3)).By analyzing these sub-channels in detail,we find that the kinetic energy of the ionized electron increases with the increase of s1value.It is easy to understand that the more the IR photons that are absorbed by the electron in the ATI process,the greater the kinetic energy that can be obtained by the electron.Besides,we can see that the location of dip structure in channel(1,1)corresponds to the waist structure of the sub-channels.Therefore,we find that the waist structure on the spectrum of channel(1,1)comes from the interference of all sub-channels,while the dip structure on it is attributed to the sum of the waist structures of all sub-channels.

Fig.8.Angle-resolved HATI spectrum for(a)channel(1,1)and(b)channel(2,0).

In order to further understand the interference fringes of the sub-channel in Fig.9,we use the saddle-point approximation and analyze the Bessel functionin the LAC process in Eq.(10).[23]The Bessel function can be reduced into For the Bessel function Jq1(ζn1−ζf1),it can be expressed in an integral form

where T=2π/ω1.Like Eq.(15),now equation(19)becomes

with g(t)=Scl(t,Pn)−Scl(t,Pf)+q2ω2t. By using the saddle-point approximation,equation(20)can be rewritten as

where G(t0)=(ζn1−ζf1)sin(ω1t0)−q1ω1t0,and the saddle point t0satisfiesand

Hence we obtain the energy conservation in the LAC process to be

where q2represents the absorbed number of XUV photons in the LAC process. Finally,we derive the electron energy as follows:

Fig.10.Comparison between interference fringes of sub-channel(1|2,1)and classical orbits with different values of collision moment tc=cos(ω1t0)and θn=20◦.

Additionally,we may find that the destructive interference appears in the angular distribution if cos[F(t0)−π/4]=0,which indicates that the interference pattern is attributed to the interference of different orbits at collision moments t0and 2π/ω1−t0in the LAC process.

The above analysis shows that the classical orbits are consistent well with the interference fringes,and these energy orbits come across the position of waist structure,where q1=0 and the LAC process requires the electron momentum to be unchanged. However,electron absorbs one XUV photon in the LAC process in the sub-channel(1|2,1).Therefore,it obviously violates the momentum conservation law and this is the reason why the destructive interference fringes occur at the waist structure.Meanwhile,equation(23)infers the waist structure shift rightwards with the increase of the laser ellipticity at the same θn.

4.Conclusions

In this work,we investigate the angle-resolved ATI spectra of H atom in different cases of EP two-color laser fields by using the frequency-domain theory.The spectrum presents asymmetric multi-plateau structure with interference fringes.It is found that the spectrum shifts with the increase of the EP degree of the IR laser filed.Furthermore,we study in detail the interference fringes on direct ATI spectrum and rescattering ATI spectrum.For the direct ATI,the fringes on the spectrum are mainly attributed to the fact that the ionization probability is very small when the direction of emitted electrons is perpendicular to the XUV laser polarization.While the interference fringes in the rescattering ATI spectrum mainly comes from the superposition of the waist structures of all sub-channels.This is because the electron momentum cannot keep conserved before and after the collision at some directions that the electron emits in the LAC process of the sub-channels. On the whole,the direct ATI process makes a dominate contribution to the whole ATI process.

Acknowledgment

We thank all the members of strong-field atomic and molecular physics(SFAMP)club for their helpful discussion and suggestion.