Spectral Shifts of Nonadiabatic High-Order Harmonic Generation

High-order harmonic generation (HHG) is a nonlinear nonperturbative process in ultrashort intense laser-matter interaction. It is the main source of coherent attosecond (1 as = 10−18 s) laser pulses to investigate ultrafast electron dynamics. HHG has become an important table-top source covering a spectral range from infrared to extreme ultraviolet (XUV). One way to extend the cutoff energy of HHG is to increase the intensity of the laser pulses. A consequence of HHG in such intense short laser fields is the characteristic nonadiabatic red and blue shifts of the spectrum, which are reviewed in the present work. An example of this nonperturbative light-matter interaction is presented for the one-electron nonsymmetric molecular ion HeH2+, as molecular systems allow for the study of the laser-molecule orientation dependence of such new effects including a four-step model of MHOHG (Molecular High-order Harmonic Generation).


Introduction
Atomic high-order harmonic generation (HHG) was discovered in 1980s [1,2].It is a nonlinear nonperturbative process, which converts low-frequency fundamental femtosecond (1 fs = 10 −15 s) laser fields to coherent high-energy attosecond (1 as = 10 −18 s) pulse trains or an isolated attosecond laser pulse [3].The typical laser intensity I used for HHG is around 10 14 W/cm 2 .The general feature of atomic HHG which also appears in molecular high-order harmonic generation (MHOHG) is a fast OPEN ACCESS decay of the lower order harmonics, then a long plateau, and a short cutoff [4].The mechanism of HHG is interpreted by a three-step model [5][6][7] neglecting the role of excited states due to their short lifetime in strong laser field.First, the bound state electron is freed by tunnel ionization through the static barrier created by the instantaneous laser field and Coulomb potential.Then, the electron oscillates in the laser field and acquires kinetic energy.It may finally return to the parent ion and recombine to the ground state to emit HHG photons.The maximal kinetic energy k E of electrons returning to the origin is 3.17 p U , where 3.17 In MHOHG, the ionized electron can recombine with neighboring ions in the same molecule leading to a much larger cutoff energy depending on the internuclear distance R [8][9][10].This larger cutoff energy comes from longer acceleration distance in the above process.
From Equation ( 1), one has to increase the ionization potential p I , laser intensity I and wavelength  to extend the cutoff energy cutoff


of the HHG spectrum.In short laser fields with very high intensity, the depletion of the ground state and rapid change of intensity ( ) I t will lead to nonadiabatic effects in HHG, such as red [11] and blue [12] shifts of HHG.We review the recent progress in these nonadiabatic spectral shifts of HHG and in particular MHOHG in this article.

Blue Shift in HHG
Earlier work has studied the blue shift of HHG.It is the result of mainly two mechanisms.One is the propagation effect of a fundamental laser pulse in an ionized medium [13,14].The other is the nonadiabatic response of the electronic dipole to the rapid change of laser fields [15][16][17].

Blue Shift by Propagating Effect
Free electrons created during the laser pulse induce a change of the refractive index [13

Nonadiabatic Blue Shift in HHG
Rapid change of a laser pulse shape will lead to a chirp of HHG.The electronic dipole response to the external laser field is nonadiabatic.This effect is intuitive in the above three-step model.As shown in Figure 1, the laser intensity is rising in the first half of the pulse.The effective amplitude of the successive laser cycle is higher than the previous one.The electron gains more energy from the latter laser field compared to the previous cycle.In other words, the rising laser envelope (RLE) with t < 0 induces a positive-chirp role in HHG, thus leading to a blue shift in HHG [16].In contrast, a negative chirp of HHG occurs on the falling laser envelope (FLE) with t > 0. This causes a red shift of HHG.However, no net shift of HHG is observed with low laser intensity since the blue shift generated on the RLE and the red shift produced on the FLE will result in broadening of the harmonics, but no net change in the central frequency.When the laser intensity is above the ionization saturation threshold, the depletion of the ground state has to be included.Especially for long intense laser pulses, the sample is highly ionized on the RLE, little HHG is emitted on the FLE.This imbalance of HHG will lead to a net blue shift in HHG.This has been experimentally confirmed in [12,19].5)).The total pulse duration τ is 20 optical cycles.

Red Shift of HHG in Resonant Systems
From the discussion in the above section, the depletion of the ground state on the RLE will prevent the HHG emission on the FLE.It is therefore hard to observe a noticeable red shift of HHG directly in such a case.However, enhanced excitation (EE) in resonant systems will change the above three-step model to a four-step model.We take the resonant nonsymmetric system HeH 2+ as an example to explore the nonadiabatic red shift in MHOHG [11,20,21].This one-electron model also emphasizes the importance of laser-induced electron transfer, LIET [20], in molecular spectra.

Four-Step Model of MHOHG in Resonant Systems
We study numerically the following time-dependent Schrödinger equation (TDSE): where 0 H is the field-free Hamiltionian [22], χ is the angle between the laser polarization and the molecular axis z in the xz plane.The electric field ( ) t E of the laser pulse is given by: 0 where the pulse shape The initial state for time evolution is the ground state 1s .The laser polarization is along the z molecular axis.The power spectra of MHOHG is calculated by a Fourier transformation of the dipole momentum in the acceleration form ( ) A d t , as it is the most reliable numerical method for strong field interactions avoiding after pulse transient effects such as in the dipole form [24].To further probe the temporal structures of MHOHG, we perform a time profile analysis of the harmonic spectra.The time profile of harmonic q  is obtained by a wavelet analysis as in signal processing [25][26][27]: with the wavelet kernel ' ' , ( ) ( ( )) The mother wavelet we use is a Morlet wavelet: The MHOHG spectrum of HeH 2+ in a laser field at wavelength 400 nm and intensity 15 3.5 10 I   W/cm 2 with a duration of 15 optical cycles is shown in Figure 2a.The corresponding time profile of the spectrum is presented in Figure 2b.Due to the break of the symmetry of the system, both odd and even order harmonics are generated as illustrated in Figure 2a.A typical difference from the usual HHG spectra of neutral atoms and molecules is a strong resonance around harmonic order 15 in Figure 2a which agrees well with the dressed energy difference between the first excited state and the ground state [20].One can also observe a strong resonance range around harmonic order 15 in the time profile analysis in Figure 2b.The conversion efficiency of the resonant harmonics is about 3 orders higher than the harmonics in the long plateau, which can be used to generate bright attosecond pulses (trains).The strong resonance in MHOHG (Figure 2a) cannot be explained consistently by the three-step model [5].
We take the resonant EE of long lifetime excited states into account explicitly and generalize the three-step model to a four-step model as presented in Figure 3.The electron is preexcited to a localized long lifetime resonant bound state, tunnel ionized from the resonant excited state, driven by the laser field, followed by a recombination to the ground state to emit MHOHG.This four-step model is different from the four-step mechanism proposed by Strelkov [28].In the latter model, the resonant state is an autoionizing state embedded in the continuum.The sequence of the above two models is also different.In the latter, the electron is ionized from the ground state first, then oscillation in the laser field, transition the resonant state, and recombination to the ground state finally.
Resonant MHOHG described by the above four-step model have novel features [20,21], which do not occur in the three-step model of atomic HHG.For example, resonance enhanced harmonic generation, multiple frequency harmonic series, double plateau, and a larger cutoff energy.Another important feature presented in the inset of Figure 2 is the redshift of MHOHG above the ionization threshold.No such shift appears in the lower-order MHOHG in the perturbative regime with photon energy less than the ionization potential p I since it is of pure molecular bound state character [29].
It clearly shows that the red shift is a nonadiabatic effect.
To understand the red shift in detail, we explore the role of EE and enhanced ionization (EI) [30,31] in MHOHG.The latter is a nonlinear characteristic effect of molecules in intense laser fields.We present the population of the resonant 2 p state and the total ionization rate ( ) t  as a function of time in Figure 4, from which one concludes that EE and EI are not synchronous.For laser intensity below the saturation threshold, the excitation to the localized resonant state 2 p on the RLE will strongly enhance the ionization rate of the system on the FLE.In this way, the red shift generated on the FLE will be larger than the blue shift on the RLE.An overall noticeable red shift can be observed in the MHOHG process.The time delay between EE and subsequent EI plays an important role in producing a noticeable "redshift" in MHOHG for very short pulses.Step (1): the laser field pumps the electron from the ground 1s state to the localized long lifetime 2 p excited state; part of the electron is ionized from the excited state in step (2), then accelerated by the laser field in step (3), and then transits back to the ground state to emit photons in step (4) when the laser field changes its direction.Part of the electron in the excited state will directly transit to the ground state in step (2*) when the laser field changes its phase.This leads to a strong resonance in the harmonic spectra.The Stark shift of the energy levels 1 E and 2 E is included [11]. 1 E  and 2 E  refer to the dressed energies with the field along +z direction, while 1 E  and 2 E  refer to the dressed energies with the field along −z direction.

Pulse Duration Dependence of Red Shift in MHOHG
The red shift in MHOHG is a nonadiabatic process as it depends on the pulse shape ( ) f t in Figure 2b.An important parameter in this process is the pulse duration in addition to the laser intensity.This nonadiabatic effect is strongly related to ( ) / dI t dt .Shorter pulses enhance this effect, while this nonadiabatic effect will be less noticeable for longer pulses.To emphasize this effect, MHOHG spectra for harmonic order 40 to 42 in the plateau produced by different duration of laser pulses are illustrated in Figure 5a.The laser intensity and wavelength are kept the same as those in Figure 2. The corresponding time profiles of harmonic with order 40 are presented in Figure 5b.The nonadiabatic spectral red shift of MHOHG gradually decreases as the laser-pulse duration increases from 10 to 20 optical cycles.From the time profile analysis, the harmonic generation process is gradually shifted from the falling part of the pulse to the rising part in agreement with the above prediction.

Orientation Dependence of Red Shift in MHOHG
From the above discussion, the resonance of 2 p state with long life time plays a key role in the red shift of MHOHG.The radiative transition from the ground 1s state to the resonant 2 p state is however restricted to the parallel component of the laser field.In this subsection, we explore the orientation dependence of the red shift in MHOHG [22].The numerical results of MHOHG spectra with different orientation angles χ and the corresponding time profiles are presented in Figure 6a,b, respectively.
An obvious feature in Figure 6a is that the overall intensity of MHOHG is strongly dependent on the resonance of the excited state 2 p .The intensity of MHOHG in the plateau with orientation angle =0   (parallel) is more than two orders higher than that with orientation angle =90   (perpendicular).
The reason is that the direct ionization rate from the ground state 1s is very low due to the large ionization potential.However, the EE to the resonant excited state 2 p will lead to EI, which greatly increases the ionization rate of the system.For the parallel orientation, the amplitude of radiative transition to the 2 p state is maximal.This corresponds to an efficient laser-induced electron transfer (LIET) in a nonsymmetric system [20].One can observe a strong resonance in the MHOHG spectra.
When we gradually increase the orientation angle χ between the molecular axis and laser polarization, the transition amplitude along the molecular axis direction diminishes, and the intensity of the resonance and the plateau in MHOHG decreases.It can be observed in the corresponding time profiles in Figure 6b.For perpendicular orientation with =90   , the transition from the 1s state to 2 p state is forbidden by selection rules.The electron may transit from the 1s state to 2 p state (−0.7388 a.u.(20.1 eV)).However, the 2 p state is not a resonant state with long life time.Its role in the MHOHG is therefore negligible.As a consequence, the above four-step model gradually changes to a three-step model with increasing angle χ, and the time profiles of higher-order harmonics illustrated in Figure 6b becomes more symmetric with respect to the time when the laser pulse reaches its peak intensity, i.e., 0 t  .The red shift of MHOHG decreases with the increase of the orientation angle χ.No noticeable red shift occurs in the MHOHG spectra when =90   as illustrated in the calculations presented in Figure 6a for the intensities and in Figure 6b for the time profiles.One may expect only odd-order harmonic generation due to the symmetry of the system with =90   .However, a little population on the 2 p state will break the symmetry of the system, which will lead to even-order harmonic generation [22].
energy of the electrons oscillating in the laser field, E and  are the amplitude and frequency of the laser field, respectively (atomic units e e m  =  = 1 are used throughout the paper unless otherwise specified).As a consequence, the cutoff energy of HHG, i.e., maximal kinetic energy k E plus ionization potential p I , is:

Figure 1 .
Figure 1.Illustration of the positive (solid line) and negative (dotted line) chirp of High-order harmonic generation (HHG) by the rapid change of the electric field E(t) (Equation (5)).The total pulse duration τ is 20 optical cycles. 0

Figure 2 .
Figure 2. (a) Molecular high-order harmonic generation (MHOHG) spectra of HeH 2+ in linearly polarized laser fields; (b) wavelet analysis of the harmonic generation.The insets in (a) are the enlarged areas of the MHOHG spectrum.The left inset shows the harmonics in the purtabative regime without shift.The right inset presents the harmonics near the cutoff with noticable red shift.The internuclear distance R = 4 a.u. is fixed.The initial state is the ground 1s state.The laser polarization is along the z axis, i.e., the orientation angle 0    .The peak laser intensity is

Figure 3 .
Figure 3. Four-step model for MHOHG in the resonant molecular ion system HeH 2+ .Step (1): the laser field pumps the electron from the ground 1s state to the localized long lifetime 2 p excited state; part of the electron is ionized from the excited state in step (2),

- 1 HHGFigure 4 .
Figure 4. Population of the first excited state 2 p (solid black line) and the total ionization rate (dotted red line) as a function of time.The laser parameters are the same as those in Figure 2.

Figure 5 .
Figure 5. (a) Red shifted harmonics in the plateau with order N = 40-42 at different laser pulse durations τ.With increasing the pulse duration, the nonadiabatic effect is reducing and the redshift of the harmonics becomes smaller.(b) The corresponding time profiles of harmonic order N = 40.When the laser pulse duartion is increased, more harmonics generated with 0 t  , and the asymmetry between 0 t  and 0 t  is less, which reduces the red shift of the harmonics.The curves have been shifted vertically for clarity.

Figure 6 .
Figure 6.(a) Orientation dependence of MHOHG in linearly polarized laser fields.(b) The corresponding time profile analysis.With increasing orientation angle, the intensity of the resonance in MHOHG around order 18 becomes weaker and weaker.The redshift in the inset of (a) is becoming less noticeable.For higher order MHOHG (larger than 40), the time-profile analysis in (b) shows increasing symmetry with increasing orientation angle.
[21]the life time is around 4 ns[23].)resonantwiththe ground state.The MHOHG spectra of HeH 2+ are obtained by numerically solving the TDSE.For the computational details of this one-electron system, we refer to reference[21].The internuclear distance R is fixed at 4 a.u.(near the excited state minimum R =3.89 a.u.).The energy of the first excited state 2 p is −1.03 a.u.(28 eV).