Odd–even harmonic emission from asymmetric molecules: Identifying the mechanism?

Jianguo Chen(陳建國(guó)),Shujuan Yu(于術(shù)娟),Yanpeng Li(李雁鵬), Shang Wang(王賞),and Yanjun Chen(陳彥軍)

College of Physics and Information Technology,Shaan’xi Normal University,Xi’an 710119,China

Odd–even harmonic emission from asymmetric molecules: Identifying the mechanism?

Jianguo Chen(陳建國(guó)),Shujuan Yu(于術(shù)娟),Yanpeng Li(李雁鵬), Shang Wang(王賞),and Yanjun Chen(陳彥軍)?

College of Physics and Information Technology,Shaan’xi Normal University,Xi’an 710119,China

We study odd–even high-order harmonic generation(HHG)from oriented asymmetric molecules HeH2+numerically and analytically.The variational method is used to improve the analytical description of the ground-state wave function for the asymmetric system,with which the ground-state-continuum-state transition dipole is evaluated.The comparison between the odd–even HHG spectra and the improved dipoles allows us to identify and clarify the complex generation mechanism of odd–even harmonics from asymmetric molecules,providing deep insights into the relation between the odd–even HHG and the asymmetric molecular orbital.

odd–even harmonics,asymmetric molecules,dipoles,variational method

1.Introduction

In recent years,high-order harmonic generation (HHG)[1,2]has been a subject of great research interest in strong laser–matter interactions.The HHG has important applications in attosecond science.[4]It can be well described by a three-step model,[5,6]in which harmonics are emitted through the processes of tunneling ionization,propagation,and recombination of the electron with the nuclei.For atoms and symmetric molecules,the emission of harmonics includes only the odd component due to the central symmetry of the atomic and molecular potentials.[7–13]For asymmetric molecules,[14–18]due to the absence of the inversion symmetry,both odd and even harmonics are emitted.[19–24]

Present studies show that odd and even harmonics from asymmetric molecules possess different spectral properties[25–27]and carry different information of the target.[28–33]They therefore need to be studied separately.A simple model[28]has been proposed to describe the generation mechanism of the odd and even harmonics,where the emissions of odd(even)harmonics from asymmetric HeH2+molecules are considered to be closely associated with the socalled odd(even)transition dipole between the even(odd) component of the asymmetric orbital and the continuum state.This model has been applied to tomographic reconstruction of the asymmetric orbital of the CO molecule with HHG,[33]the analysis of two-center-interference-induced minima in the odd–even HHG spectra,[34]the polarization properties of the odd–even harmonics,[35]the time-resolved study of the odd–even HHG process,[36]and attosecond probing of the vibrational dynamics of asymmetric molecules with odd– even HHG,[37]etc.This model works better for asymmetric molecules with small internuclear distances.For large ones, the resonance between the ground state and the first excited state of the asymmetric system also plays an important role in the HHG process.[38]In this situation,to describe the odd–even HHG,a more complex model which considers this resonance effect is needed.

On the other hand,the HHG spectra from symmetric molecules show a striking minimum,[39]which has been identified as arising from the effect of two-center interference and has attracted broad interest in recent years.[40–44]It has been shown that this minimum in the HHG spectrum corresponds to the minimum in the bound-continuum transition dipole and therefore has important applications in molecular high-order harmonic spectroscopy(HHS).For example,this minimum in the HHG spectrum can be used to read the bond length of the molecule[40,41]and judge the phase of the dipole which is important in the molecular orbital tomography procedure.[45]

For asymmetric molecules,however,in many cases,the striking minimum disappears in the odd or even HHG spectrum and the reason has been attributed to the interplay of different recombination routes originating from the asymmetry of the molecular orbital.[34]Fortunately,with the simple odd–even-HHG model,[28]it is shown that when the minimum cannot be read from the odd–even HHG spectra directly,it can be probed[35]through the polarization measurement[46–53]of the odd–even HHG.Specifically,the position of the minimum in the odd or even dipole corresponds to the harmonic order for the maximal ellipticity of the odd or even harmonics.This polarization measurement of odd–even HHG alreadyshows promise for use in probing the vibrational motion of the asymmetric molecule.[37]However,the prediction of the dipole minimum by the polarization measurement does not always agree well with the theoretical evaluation.In some cases,a remarkable difference is also observed.[35]It is then natural to ask whether this remarkable difference arises from other mechanisms that go beyond the description of the simple model,or if it arises from only the inaccurate evaluation of the dipole in relevant theoretical treatments where some rough approximations are used?

In this paper,we answer the question with improving the evaluation of the dipole.We choose the simplest asymmetric diatomic molecule HeH2+as the target molecule,with varying molecular parameters such as the effective charge and the internuclear distance.In the frame of linear combination of atomic-orbitals-molecular-orbitals(LCAO-MO)approximation,we use the variational method to obtain the analytical expression for the ground-state wave function of model HeH2+.Using the developed ground-state wave function,we calculate the transition dipole between the ground state and the continuum state.Then we compare the calculated odd–even dipoles with the spectra and ellipticity of odd–even harmonics from the three-dimensional(3D)model HeH2+obtained through numerical solution of the time-dependent Schr?dinger equation(TDSE).

Our simulations show that the use of the variational method to the ground-state wave function remarkably improves the evaluation of the odd–even dipoles,in comparison with the exact ones obtained through diagonalizing the field free Hamiltonian in one-dimensional(1D)cases.In particular,the minimum in the improved odd–even dipoles agrees well with that predicted by the polarization measurements of odd–even harmonics for different molecular parameters.This agreement also holds as we compare the improved dipoles with the spectra,which are obtained with considering the transition of the continuum electrons back to only the ground state where a striking minimum also emerges.Our results verify the one-to-one matching between the odd–even HHG spectra and the odd–even dipoles,and give important suggestions on the generation mechanism of the odd–even harmonics.

2.Theoretical descriptions

2.1.TDSE simulations of odd–even HHG

The Hamiltonian of the model HeH2+studied here is H(t)=P2/2+V(r)+r·E(t)(in atomic units ofˉh=e= me=1).We assume that the molecular axis is located in the xoy plane and the laser field is linearly polarized along a direction parallel to the x axis.The potential used here has the formHere(r?R2)2=(x+R2cosθ)2+(y+R2sinθ)2+z2.Z1and Z2are the effective charges.R1and R2are the positions of the He and H nuclei to the origin,respectively,with R1= Z2R/(Z1+Z2)and R2=Z1R/(Z1+Z2).R is the internuclear separation.ξ=0.5 is the smoothing parameter,and θ denotes the angle between the molecular axis and the laser polarization.For different R,the effective charges Z1and Z2are adjusted in such a manner that the ionization potential of model HeH2+reproduced in our 3D simulations is 1.1 a.u. E(t)=exE(t),where E(t)=f(t)E sinω0t is the external electric field,and exis the unit vector along the x axis.f(t)is the envelope function.E and ω0are the amplitude and the frequency of the external electric field,respectively.In our calculations,we use trapezoidally shaped laser pulses with a total duration of 10 optical cycles and linear ramps of three optical cycles.Numerically,the above Schr?dinger equation is solved by the spectral method.[54]We work with a grid of sizes Lx×Ly×Lz=409.6 a.u.×51.2 a.u.×51.2 a.u.for the x,y,and z axes,respectively.The laser wavelength used here is λ=800 nm and the laser intensity is I=8×1014W/cm2. In each time step,the TDSE wave function ψ(t)of H(t)is multiplied by a mask function to absorb the continuum wave packet at the boundary.The mask function along the x axis has the form F(x)=cos1/8[π(|x|?x0)/(Lx?2x0)]for|x|≥x0and F(x)=1 for|x|＜x0.Here,x0is the boundary of the absorbing procedure along the x axis.For the accurate simulations,we have used x0=Lx/8.The situation is similar for other dimensions of y and z.With the present laser parameters, the maximal classical displacement[5]of the electron along the laser polarization is xmax=2E/w20≈93 a.u.,which can easily be represented in our numerical grids.The absorbing procedure of x0=Lx/8 is also sufficient to retain the contributions of long–short trajectories and multiple returns[6]to HHG.

Alternatively,we can set the boundary of the absorbing procedure along the x direction(i.e.,the direction of the laser polarization)as x0=1.2E/w20with y0=Ly/8 and z0=Lz/8 in the y and z directions unchanged.This treatment removes the contributions of the long trajectory and multiple returns,and the short-trajectory contributions are not influenced.We therefore can analyze the ellipticity of harmonics only arising from the short trajectory.The value of x0=1.2E/w20corresponds to the maximal displacement of the electron as it travels in the laser field following the short trajectory.

According to the three-step model,it is well known that in each laser cycle,the long and the short electron trajectories contribute significantly to the HHG,which can be modulated through the propagation effect[55]or using phase-stabilized driving pulses.[56]It has been shown that the interference of the long and short electron trajectories[57]greatly influences the ellipticity of the harmonics.[58]As a result,the ellipticity of odd or even harmonics from asymmetric molecules shows strong oscillation and it is difficult to identify the angle dependence of the ellipticity.[35]For these reasons,in this paper,we consider only the short-trajectory TDSE results,for which the elipticity of the odd–even HHG spectra shows a clear dependence on the orientation angle.

Once the TDSE wave function ψ(t)is obtained,the coherent part of the HHG spectrum,parallel or perpendicular to the laser polarization,can be evaluated using

where ω is the emitted-proton frequency.To study the influence of the excited states on the HHG,the following expression is also used in the evaluation of the HHG spectrum:

which denotes the transition of the continuum electron back to the ground state|0〉with the amplitude a0(t)=〈0|ψ(t)〉. Here,the contributions of the excited states are excluded.

The ellipticity of HHG is determined by the amplitude ratio and the phase difference of the parallel and perpendicular harmonics

whereμ=S⊥/S‖and δ=φ⊥?φ‖.The intensity and phase of the harmonic components are given by S‖(⊥)=|F‖(⊥)(ω)|2and φ‖(⊥)(ω)=arg[F‖(⊥)(ω)].The range of the ellipticity is 0≤ε≤1.The linear,elliptical,and circular polarizations correspond to ε=0,0＜ε＜1,and ε=1,respectively.

2.2.Analytical description of odd–even HHG

According to the simple model,[28]the dipole Dodd(ω,θ) that is mainly responsible for the emission of odd harmonics from HeH2+along the laser polarization e‖can be written as

with Godd(ω,θ)=a1cos(pkR1cosθ)+a2cos(pkR2cosθ). Similarly,the dipole Deven≡Deven(ω,θ)that is related to the emission of even harmonics along the laser polarization can be written as

with Geven(ω,θ)=a1sin(pkR1cosθ)?a2sin(pkR2cosθ). Here,R1=Z2R/(Z1+Z2)and R2=Z1R/(Z1+Z2).R is the internuclear separation.Z1and Z2are the effective changes of the He and H nuclei.Ipis the ground-state ionization potential of HeH2+.pkis the effective momentum of the continuum state,|p〉∝|eipk·r〉,with pk=|pk|=[2(Ip+Ep)]1/2that considers the Coulomb acceleration.[45]Epis the energy of the continuum state|p〉,agreeing with the energy conservation relation Ep=ω?Ip.In the expressions of Eqs.(4)and (5),the ground-state wave function φ1σ(r)of HeH2+in the mass-center coordinate is roughly approximated by φ1σ(r)= Nf[a1e?κra+a2e?κrb].Here,Nfis the normalization factor, ra=|r?R1|,and rb=|r?R2|.

In this paper,we improve the description of φ1σ(r)with the variational method.According to the LCAO-MO approximation,we assume that the ground-state wave function φ1σ(r)of HeH2+in the mass-center coordinate has the form φ1σ(r)=α?a+β?b.Here,α and β are the coefficients.?a=are the wave functions of He+ion and H atom,respectively,with κ1=Z1and κ2=Z2.We approximate the continuum state|p〉by the plane wave as in Eqs.(4)and(5).Then the dipole Dodd(even)(ω,θ) can be written as

2.3.Application of variational method to HeH2+

The Hamiltonian of the asymmetric molecule HeH2+is

Here,V(r)=?Z1/ra?Z2/rbis the coulomb potential with raand rbbeing the positions of the electron to the two nuclei as defined following Eq.(5).Z1and Z2are the effective charges,and R is the internuclear separation.In the frame of LCAO-MO approximation,the bound wave function of HeH2+in the mass-center coordinate can be written as

The energy of the asymmetric system can be evaluated using

Substituting Eqs.(8)and(9)into Eq.(10),we have

with

where Ea=?Z12/2 and Eb=?Z22/2 are the ground-state energies of He+ion and H atom,respectively.According to the normalization condition for ?aand ?b,we have Saa=Sbb=1 and Sab=Sba.Taking the derivative of α or β in Eq.(11),we obtain the following equations:

The secular equation for Eq.(12)is

To solve the secular equation,we have AE2+BE+C=0, withandThe solutions of the above equation are the energies of the 1σ and 2σ states

According to the normalization condition of Eq.(9),we have

Then we arrive at

Substituting Eqs.(12)and(13)(Eq.(14))into Eq.(15),we have

3.Results and discussion

In the following,we apply the variational method to model HeH2+with different molecular parameters such as the effective charges Z1and Z2and the internuclear distance R. We make comparisons between the previous dipoles and the improved dipoles.Then we use these dipoles as a benchmark to analyze the HHG mechanism of odd–even harmonics from asymmetric molecules.

3.1.1σ-state and 2σ-state energy curves and improved ground-state wave function

In Fig.1,we compare the R-dependent ground-state energy of HeH2+obtained by the variational method(blacktriangle)with that obtained by the exact numerical simulations (red-circle).The results of the two different methods show good agreement with each other.

In Table 1,we show the coefficients of α and β obtained with the variational method for HeH2+at different R. As shown in Table 1,the parameter α1is much larger than β1,implying that for the 1σ state of HeH2+,the probability density distribution around the helium nucleus is much higher than that around the hydrogen nucleus.In contrast,for the 2σ state,the electron is mainly located around the H nucleus(the parameter α2is smaller than β2).

Fig.1.(color online)R-dependent 1σ-state(a)and 2σ-state(b)energy of model HeH2+with Z1=2 and Z2=1 obtained by two different methods:the variational method with Eqs.(13)and(14)(black-triangle)and the exact numerical simulation(red-circle).

Table 1.Coefficients of the wave functions of Eq.(18)for HeH2+obtained using the variational method with Z1=2 and Z2=1 at different internuclear separations R.α1 and β1 obtained using Eq.(16)are the coefficients for the 1σ state.α2 and β2 obtained using Eq.(17)are the coefficients for the 2σ state.

These coefficients also show that the 1σ wave function φ1σ(r)=α?a+β?bobtained with the variational method differs remarkably from the previous one φ1σ(r)=Nf[a1e?κra+ a2e?κrb].One can expect that this difference will play an important role in the calculation of the odd–even dipoles,as shown below.

The accurate calculation of the bound-continuum transition dipole〈0|r|p〉in 3D cases is difficult,as it needs the knowledge of both the bound state|0〉and the continuum state|p〉.Here,we evaluate the exact dipole in 1D cases where all of the bound and continuum states of the system can be obtained by diagonalizing the field-free Hamiltonianat the fixed-nuclei approximation.The 1D asymmetric Coulomb potential used here has the form V(x)=with ξ=0.5 and Z2=Z1/2.R1=Z2R/(Z1+Z2)and R2=Z1R/(Z1+Z2). For different internuclear distances R,we adjust the parameter Z1in such a manner that the ionization potential of the 1D asymmetric system reproduced in our numerical simulations equals 2.25 a.u.(the ionization potential of).When the calculated continuum state has an odd-like(even-like)parity, through the expression〈0|r|p〉,we obtain the corresponding odd(even)exact dipole,[28]as shown in Figs.2(a)and 2(b). The parameters Z1and Z2used in the numerical simulations are also used in evaluating the analytical dipoles with Eqs.(4)–(7)in Fig.2.

Fig.2.(color online)Comparisons of odd(solid-black)and even (dashed-red)dipoles for HeH2+with I p=2.25 a.u.,Z1/Z2=2 at R=2 a.u.(the left column)and R=4 a.u.(right),obtained by different methods.In panels(a)and(b),we show the 1D exact dipoles |〈0|x|p〉|2/ω4 for the continuum state|p〉having the odd-like(corresponding to odd dipoles)or even-like(even dipoles)parity.The continuum state|p〉has the energy E p=ω?I p.Here,ω=nω0 is the assumed harmonic energy with ω0=0.057 a.u.(λ=800 nm).In panels(c)and(d),we show the dipoles|D odd(even)(θ,ω)|2/ω4 with θ=0° obtained using Eqs.(4)and(5)(the previous dipoles).In panels(e) and(f),we show those obtained using Eqs.(6)and(7)(the improved dipoles).The log10 scale is used here.

One can observe from the intersections of the odd(solid black)versus even(dashed-red)dipoles in Fig.2 that the improved dipoles agree better with the exact ones.We use the vertical-dashed arrows to indicate several of these intersections.For the positions of the minima in the odd or even dipoles,the improved ones are also closer to the exact ones, especially for the minima in the odd dipoles,as indicated by the vertical-solid arrows in the right column of Fig.2.These results verify the applicability of the ground-state wave function obtained with the variational method in calculating the odd–even dipoles for HeH2+.

It should be stressed that although the use of the variational method improves the agreement between the exact and the analytical odd–even dipoles.There are still some differences between them,both in the positions of the minima in the dipoles and the intersections of the odd versus even dipoles. These differences can arise from the plane-wave approximation with the effective momentum pk[59]for the continuum electron in calculating the dipoles,where the Coulomb effect is not well described.One can expect that this Coulomb effect is stronger for molecules with larger ionization potentials.

3.2.Comparison between dipoles and spectra

It has been shown[35]that for model HeH2+with small internuclear distances R,the previous dipole of Eq.(4)associated with odd harmonics shows a striking minimum.However, the minimum is absent in the odd TDSE spectrum.By contrast,the ellipticity curve of odd harmonics shows a maximum at one harmonic order at which the minimum appears in the odd dipole.The ellipticity of harmonics therefore can be used as a tool to probe the position of the minimum in the dipole. This dipole minimum is important,as it encodes the information of the molecular structure and has potential applications in asymmetric molecular orbital imaging.[33]However,in some cases,the position of the minimum in the dipole does not agree well with the harmonic order for maximal ellipticity.To use the ellipticity measurement of harmonics as a tool to judge the dipole minimum,it is necessary to clarify this disagreement.

In Fig.3,we compare the spectra,ellipticity,and dipoles of odd harmonics for a model HeH2+molecule with Ip= 1.1 a.u.,Z1/Z2=1.5,and R=1.5 a.u.at different orientation angles θ.The spectra are obtained through 3D TDSE simulations.

Fig.3.(color online)Comparisons of(a)–(c)spectra,(d)–(f)ellipticity,and(g)–(i)relevant dipoles of odd harmonics for 3D model HeH2+ with I p=1.1 a.u.,Z1/Z2=1.5,and R=1.5 a.u.at θ=10°(the left column),θ=30°(middle),and θ=50°(right).The odd spectra are calculated using Eq.(1)(accurate spectra,solid-black)and Eq.(2)(spectra associated with the transition of the electron back to only the ground state,dashed-red)with short-trajectory 3D TDSE simulations.The odd dipoles are calculated using Eq.(4)(previous dipoles,solid-black)andEq.(6)(improved dipoles,dashed-red).The log10 scale is used to show the spectra and dipoles.

First,one can observe that the accurate TDSE spectra of Eq.(1)(solid-black)in the first row of Fig.3 do not show a striking minimum.By contrast,the spectra approximated with Eq.(2)(dashed-red),where the transition of the continuum electron back to only the ground state is considered,show a striking minimum.The position of the spectral minimum shifts towards higher harmonic orders as the orientation angle increases.The comparisons imply that the transition of the continuum electron back to the excited states plays an important role in the HHG of asymmetric molecules,resulting in the disappearance of the minimum in the accurate TDSE spectrum of Eq.(1).

Secondly,corresponding to the spectral minima in the first row of Fig.3,the ellipticity curves of harmonics in the second row of Fig.3 show a striking maximum.In particular, the maximal ellipticity appears at the harmonic order at which the improved dipoles of Eq.(6)(dashed-red)in the third row of Fig.3 show a striking minimum,as indicated by the vertical arrows.In contrast,the positions of the minima in the previous dipoles of Eq.(4)(solid-black)in the third row of Fig.3 differ remarkably from the improved ones and this difference is more remarkable for larger orientation angles θ.These comparisons reveal that there is a one-to-one matching between the dipole minimum and the ellipticity maximum.The previous dipoles of Eq.(4)with a rough approximation for the ground-state wave function underestimate the positions of the minima for about 10 to 30 harmonic orders here.

To check our results,we also perform simulations at other molecular parameters,as shown in Figs.4 and 5,where wefix the value of Z1/Z2and change the value of the internuclear distance R.On the whole,in all cases,the position of the minimum in the improved dipole agrees with the harmonic order at which the maximal ellipticity emerges.By contrast,the position of the minimum in the pervious dipoles differs remarkably from the harmonic order of the maximal ellipticity in some cases.It should be mentioned that for the case of the large angle of θ=50°in Fig.4(i),the improved dipole has no minimum.The corresponding HHG spectra in Fig.4(c)also do not show a minimum.In this case,the correspondence between the improved dipole and the ellipticity is not very obvious.However,as we increase the laser wavelength(such as λ=900 nm)with extending HHG cutoff,both minima appear in the HHG spectra and the improved dipole, and the correspondence between the improved dipole and the ellipticity becomes striking.In addition,we also fix the value of R and change the value of Z1/Z2,the close relation between the dipoles,spectra,and ellipticity is also observed in our extended simulations.For simplicity,in the above discussion, we have chosen the HeH2+with an active electron as the target molecule.This close relation is also expected to appear for other asymmetric molecules such as CO with more electrons,for which the odd–even dipoles also give a good prediction of the relative yields of odd–even harmonics.[33]When applying the variational method to CO to obtain the analytical expression of the 1σ ground-state wave function,the multiple-electron effect needs to be considered.

Fig.4.(color online)Same as Fig.3,but for modelwith R=1.7 a.u.and Z1/Z2=2.

Fig.5.(color online)Same as Fig.3,but for modelwith R=2 a.u.and Z1/Z2=2.

4.Conclusion

We have studied the HHG from asymmetric molecules, focusing on the generation mechanism of odd–even harmonics.In the frame of LCAO-MO approximation,we applied the variational method to obtain the analytical expression for the ground-state wave function of model HeH2+.The obtained ground-state wave function remarkably improves the agreement between odd–even dipoles and ellipticity of odd–even harmonics.Specifically,the position of the minimum in the dipoles agrees with the harmonic order for the maximal ellipticity.As the minimum is usually absent in the HHG spectra of asymmetric molecules,our work supports the conclusion that the ellipticity measurement of harmonics can be used as a tool to probe the position of the minimum.The one-to-one matching between the dipole minimum and the ellipticity maximum for odd or even harmonics also sheds light on the complex generation mechanism of odd–even harmonics from asymmetric molecules.

[1]McPherson A,Gibson G,Jara H,Johann U,Luk T S,McIntyre I A, Boyer K and Rhodes C K 1987 J.Opt.Soc.Am.B 4 595

[2]L’Huillier A,Schafer K J and Kulander K C 1991 J.Phys.B 24 3315

[3]Antoine P,Anne L and Lewenstein M 1996 Phys.Rev.Lett.77 1234

[4]Krausz F and Ivanov M 2009 Rev.Mod.Phys.81 163

[5]Corkum P B 1993 Phys.Rev.Lett.71 1994

[6]Lewenstein M,Balcou P,Ivanov M Y,L’Huillier A and Corkum P B 1994 Phys.Rev.A 49 2117

[7]Lein M,Hay N,Velotta R,Marangos J P and Knight P L 2002 Phys. Rev.A 66 023805

[8]Lein M,Corso P P,Marangos J P and Knight P L 2003 Phys.Rev.A 67 023819

[9]Lagmago K G and Bandrauk A D 2005 Phys.Rev.A 71 053407

[10]Zhou X X,Tong X M,Zhao Z X and Lin C D 2005 Phys.Rev.A 72 033412

[11]Wang B,Cheng T,Li X,Fu P,Chen S and Liu J 2005 Phys.Rev.A 72 063412

[12]Chen J,Chu S I and Liu J 2006 J.Phys.B 39 4747

[13]Song Y,Li S Y,Liu X S,Guo F M and Yang Y J 2013 Phys.Rev.A 88 053419

[14]Kamta G L and Bandrauk A D 2005 Phys.Rev.Lett.94 203003

[15]Wu J,Zeng H and Guo Cl2006 Phys.Rev.A 74 031404

[16]Rupenyan A,Kraus P M,Schneider J and W?rner H J 2013 Phys.Rev. A 87 031401

[17]Akagi H,Otobe T,Staudte A,Shiner A,Turner F,D?rner R,Villeneuve D M and Corkum P B 2009 Science 325 1364

[18]Shi Y Z,Zhang B,Li W Y,Yu S J and Chen Y J 2017 Phys.Rev.A 95 033406

[19]Etches A and Madsen L B 2010 J.Phys.B 43 155602

[20]Augstein B B and Faria C F D M 2011 J.Mod.Opt.58 1173

[21]Pan Y,Zhao S F and Zhou X X 2013 Phys.Rev.A 87 035805

[22]Zhu XS,Zhang QB,Hong WY,Lan P F and Lu P X 2011 Opt.Express 19 436

[23]Du H C,Luo L Y,Wang X S and Hu B T 2012 Phys.Rev.A 86 013846

[24]Miao X Y and Du H N 2013 Phys.Rev.A 87 053403

[25]Bian X B and Bandrauk A D 2010 Phys.Rev.Lett.105 093903

[26]Etches A,Gaarde M B and Madsen L B 2011 Phys.Rev.A 84 023418

[27]Heslar J,Telnov D and Chu S I 2011 Phys.Rev.A 83 043414

[28]Chen Y J and Zhang B 2011 Phys.Rev.A 84 053402

[29]Frumker E,Hebeisen C T,Kajumba N,Bertrand JB,W?rner HJ,Spanner M,Villeneuve D M,Naumov A and Corkum P B 2012 Phys.Rev. Lett.109 113901

[30]Kraus P M,Rupenyan A and W?rner H J 2012 Phys.Rev.Lett.109 233903

[31]Frumker E,Kajumba N,Bertrand J B,W?rner H J,Hebeisen C T, Hockett P,Spanner M,Patchkovskii S,Paulus G G,Villeneuve D M, Naumov A and Corkum P B 2012 Phys.Rev.Lett.109 233904

[32]Kraus P M,Baykusheva D and W?rner H J 2014 Phys.Rev.Lett.113 023001

[33]Chen Y J,Fu L B and Liu J 2013 Phys.Rev.Lett.111 073902

[34]Zhang B,Chen Y J,Jiang X Q and Sun X D 2013 Phys.Rev.A.88 053428

[35]Yu S J,Zhang B,Li Y,Yang S and Chen Y J 2014 Phys.Rev.A 90 053844

[36]Zhang B,Yu S J,Chen Y J,Jiang X Q and Sun X D 2015 Phys.Rev.A. 92 053833

[37]Li W Y,Yu S J,Wang S,and Chen Y J 2016 Phys.Rev.A.94 053407

[38]Chen Y J and Zhang B 2012 Phys.Rev.A.86 023415

[39]Lein M,Hay N,Velotta R,Marangos J P and Knight P L 2002 Phys. Rev.Lett.88 183903

[40]Kanai T,Minemoto S and Sakai H 2005 Nature 435 470

[41]Vozzi C,Calegari F,Benedetti E,Caumes J P,Sansone G,Stagira S, Nisoli M,Torres R,Heesel E,Kajumba N,Marangos J P,Altucci C and Velotta R 2005 Phys.Rev.Lett.95 153902

[42]Le A T,Tong X M and Lin C D 2006 Phys.Rev.A 73 041402

[43]Chen Y J,Liu J and Hu B 2009 J.Chem.Phys.130 044311

[44]Wu Y,Zhang J,Ye H and Xu Z 2011 Phys.Rev.A 83 023417

[45]Itatani J,Levesque J,Zeidler D,Niikura H,Pepin H,Kieffer J C, Corkum P B and Villeneuve D M 2004 Nature 432 867

[46]Son S K,Telnov D A and Chu S I 2010 Phys.Rev.A 82 043829

[47]Zhang X F,Zhu X S,Liu X,Wang D,Zhang Q B,Lan P F and Lu P X 2017 Opt.Lett.42 1027

[48]Zhou X,Lock R,Wagner N,Li W,Kapteyn H C and Murnane M M 2009 Phys.Rev.Lett.102 073902

[49]Li L,Wang Z,Li F and Long H 2017 Opt.Quant.Electron.49 73

[50]Levesque J,Mairesse Y,Dudovich N,Pépin H,Kieffer J C,Corkum P B and Villeneuve D M 2007 Phys.Rev.Lett.99 243001

[51]Zhai C Y,Zhu X S,Lan P F,Wang F,He L X,Shi W J,Li Y,Li M, Zhang Q B and Lu P X 2017 Phys.Rev.A 95 033420

[52]Ramakrishna S,Sherratt P A J,Dutoi A D and Seideman T 2010 Phys. Rev.A 81 021802

[53]Sherratt P A J,Ramakrishna S and Seideman T 2011 Phys.Rev.A 83 053425

[54]Feit M D,Fleck J J A and Steiger A 1982 J.Comput.Phys.47 412

[55]Antoine P,L’Huillier A and Lewenstein M 1996 Phys.Rev.Lett.77 1234

[56]Sansone G,Benedetti E,Caumes J P,Stagira S,Vozzi C,Silvestri S D and Nisoli M 2006 Phys.Rev.A 73 053408

[57]Za?r A,Holler M,Guandalini A,Schapper F,Biegert J,Gallmann L, Keller U,Wyatt A S,Monmayrant A,Walmsley I A,Cormier E,Auguste T,Caumes J P and Salières P 2008 Phys.Rev.Lett.100 143902

[58]Strelkov V V,Gonoskov A A,Gonoskov I A and Ryabikin M Y 2011 Phys.Rev.Lett.107 043902

[59]Chen Y J and Hu B 2009 Phys.Rev.A 80 033408

22 March 2017;revised manuscript

25 April 2017;published online 18 July 2017)

10.1088/1674-1056/26/9/094209

?Project supported by the National Natural Science Foundation of China(Grant No.11274090)and the Fundamental Research Funds for the Central Universities, China(Grant No.SNNU.GK201403002).

?Corresponding author.E-mail:chenyanjun@snnu.edu.cn

?2017 Chinese Physical Society and IOP Publishing Ltd http://iopscience.iop.org/cpb http://cpb.iphy.ac.cn