Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Widely tunable XUV harmonics using double IR pulses

Open Access Open Access

Abstract

Tunable attosecond pulses are necessary for various attosecond resolved spectroscopic applications, which can potentially be obtained through the tuning of high harmonic generation. Here we show theoretically, using the time-dependent Schrödinger equation and strong field approximation, a continuously tunable spectral shift of high-order harmonics by exploiting the interaction of two delayed identical infrared (IR) pulses within the single-atom response. The tuning spans more than twice the driving frequency (∼2ω) range, for several near-cutoff harmonics, with respect to only one control parameter: the change in delay between the two IR pulses. We show that two distinct mechanisms contribute to the spectral shift of the harmonic spectra. The dominant part of the spectral shift of the harmonics is due to the modulation of the central frequency of the composite IR-IR pulse with respect to delay. The second contribution comes from the non-adiabatic phase-shift of the recolliding electron wavepacket due to the change in amplitude of the subcycle electric field within the double pulse envelope. For optical few-cycle pulses this scheme can produce tunable attosecond pulse trains (APT), and in the single-cycle regime the same can be used for tuning isolated attosecond pulses (IAP). We quantify the dependence of tuning range and tuning rate on the laser pulse duration. We envision that the proposed scheme can be easily implemented with compact in-line setups for generating frequency tunable APT/IAP.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Finely tunable ultrashort XUV pulses are necessary for various applications, such as in ultrafast photoelectron spectroscopy [1,2], time- and angle-resolved photoemission spectroscopy [3,4], high resolution transient absorption spectroscopy [5] etc. High intensity tunable XUV and Xray pulses of attosecond temporal width can be obtained readily from large scale synchrotron sources [6], however for a table-top experimental set up one may rely on high harmonic generation (HHG) based XUV sources [7].

HHG based sources are tuned by changing the fundamental frequency of the femtosecond laser system, or by selecting the desired portion of the spectrum by using suitable filter [8]. These techniques may suffer from the constraint of limited tunability or loss of coherence. One efficient way of obtaining a tunable table-top XUV and soft-Xray source could be through varying the spectral shift of harmonics in the HHG process.

Previous studies have shown spectral shifts in HHG due to different mechanisms. For molecular systems spectral red-shift has been demonstrated which originated from the reduction in ionization energy between recombination time and ionization time [9], a non-adiabatic red-shift mechanism has been theoretically demonstrated in [10] where enhanced excitation of localized long lifetime excited states shifts the harmonic generation in the falling edge of the short laser pulses. For solids red-shift in the higher order plateau has been demonstrated theoretically which comes from the indirect step-by-step excitation of higher conduction bands [11].

For atomic systems spectral tuning of high-order harmonics has been achieved experimentally and theoretically through various mechanisms such as, the control of laser energy and chirp [12], using two-color laser pulses [1317], nonadiabatic blueshift of harmonics due to the generation in the rising edge of the pulse [18], due to ionization of the medium [19], and using the carrier-envelope phase (CEP) [20]. Attosecond twin-pulse control and generation of fractional higher order harmonics with energy tuning have also been demonstrated by generalized kinetic heterodyne mixing [2123]. Most of these methods are complex to implement and offer limited tuning range and only discrete tunability.

As HHG in atoms is found to be the most promising table-top candidate in terms of intensity and coherence, it is desirable to achieve the maximum possible tuning of harmonics through the temporal shaping of the driving pulse. Schemes of spectral tuning of harmonics by using two delayed IR pulses have been proposed and experimentally verified recently [2430]. Tuning of cutoff order occurs due to delay-controlled variation in the amplitude of the composite IR-IR pulses [2426,28]. Furthermore, tuning of spectral shift and width of harmonics have been shown to arise from a similar tuning of the composite pulses [27,29,30]. However, the full potential of IR-IR pulse for generating widely tunable harmonics and its capability to generate tunable isolated attosecond pulses have not yet been explored.

Here we demonstrate the tuning of harmonic spectra over $2\omega$ in the cutoff harmonics theoretically, using two delayed replica IR pulses interacting with an atomic system. This can be utilized to readily produce a tunable XUV source in a very simple manner as demonstrated in [28]. Employing the time-dependent Schrödinger equation and the strong field approximation we compute the HHG spectra with respect to the delay between the two IR pulses which depicts the prominent spectral shift. Two distinct mechanisms contribute to the spectral shift of the harmonics: the frequency modulation of the driving IR-IR pulse with respect to delay and the non-adiabatic response of the electron to the shape of the delayed IR-IR laser pulse envelope. For optical few-cycle pulses this scheme is found to give tuning of attosecond pulse trains (APT), and in the single-cycle regime the same can be used for tuning isolated attosecond pulses (IAP). We also quantify the dependence of tuning range and tuning rate on the duration of the IR pulse. The paper in organized as follows: after briefly introducing our model in section 2 we present the results and the discussion in section 3 which is followed by a concluding section.

2. Theoretical model

We solve the time-dependent Schrödinger equation (TDSE) in reduced dimensionality (1d)

$$i\frac{\partial \psi(x,t)}{\partial t} ={-}\frac{1}{2} \frac{\partial^2}{\partial x^2} \psi(x,t) + V(x) \psi(x,t) + H_I \psi(x,t),$$
where $\psi (x,t)$ is time-dependent wavefunction at position $x$ at time $t$, $V(x)$ is the atomic binding potential, $H_I=-xE(t)$ describes the atom-field interaction Hamiltonian within dipole approximation. We consider two different model potentials: (i) a Gaussian potential of the form $V(x) = -1.1 exp(-x^2/\sigma ^{2})$ with $\sigma ^{2}=1.21$ [31] and (ii) a Coulomb-like soft-core potential $V(x) = -1/\sqrt {x^2+a_0^{2}}$, with $a_0^{2}=1.4$ [32,33]. The ground state in both potentials is $-0.58$ au corresponding to the ionization potential of argon. The Gaussian potential has only a single bound-state while the soft-core Coulomb potential allows for a Rydberg series.

As the combined delayed IR-IR pulse is linearly polarized the 1d TDSE is expected to reveal the mechanism in the studied parameter regime. In addition, we have also verified our results with the calculation employing the strong field approximation due to Lewenstein for the harmonic spectra and the spectra-map. To observe the harmonic production with respect to time we have used the short time Fourier transformation, with a Gabor window. Propagation effects (plasma generation, self focusing, absorption and dispersion) are not considered here. The presented single atom response calculations could be compared with experiments with thin gas jet and loose focusing geometry [34,35].

The electric field of the combined delayed IR-IR pulse is defined as

$$E(t)=E_{01}f(t)\cos(\omega t+\phi_{ce})+E_{02}f(t-\tau)\cos(\omega (t-\tau)+\phi_{ce}),$$
where, the temporal envelope of each IR pulse is $f(t)=\cos ^2\left (\frac {t}{T_{pulse}}\right )$ and their amplitudes are $E_{01}$ and $E_{02}$. The CEP of both IR pulses is $\phi _{ce}=\pi /2$ (sin pulse) throughout this work.

3. Results and discussions

Figure 1 shows the spectral intensity of the harmonics vs the delay ($\tau$) between the two IR replica pulses (hereafter this will be called as spectra-map) for 5 optical cycles pulse (the number of complete optical oscillation is 5 in each replica IR pulse) with $E_{01}=E_{02}=0.06$ au. Figure 1(a) shows the spectra-map near the cutoff regime computed using the TDSE with the Gaussian potential and with the soft-core Coulomb potential as shown in Fig. 1(e). To independently verify the observed systematic features in the spectra-map we have performed a calculation using strong-field approximation (SFA) [36,37] as presented in Fig. 1(b). The overall agreement between the two methods is very good in this regime of harmonic order and other part of the spectra-map (not shown here) except for the fact that the SFA slightly overestimates the cut-off and hence shows a larger width of the “harmonic band”. Note that, for further analysis we take recourse of TDSE simulations for both soft-core and Gaussian potentials. In Fig. 2 we have plotted the harmonic spectra for different $\tau$ which falls in the range of delay as shown in the spectra-map (Fig. 1(a)) in the harmonic order region: 51-71 of the fundamental angular frequency ($\omega =0.057$ au).

 figure: Fig. 1.

Fig. 1. The harmonic spectra-map for a delay range $\tau = 0.8 - 1.2 T_{oc}$. (a) HHG for the Gaussian potential, and (b) SFA calculation. (c) Peak positions of the harmonic bands (solid line) vs delay, obtained from (a). The dashed lines represent delay-dependent spectral position of $N^{th}$ harmonic with respect to the central frequency of the driver. (d) The difference between the calculated maxima positions and that of the $N^{th}$ multiple of the central frequency of IR-IR driver. (e) and (f) are same as (a) and (c) for soft-core Coulomb potential.

Download Full Size | PDF

 figure: Fig. 2.

Fig. 2. The harmonic spectra (a) for different delays $\tau =$ 0.9 $T_{oc}$, 1.0 $T_{oc}$ and 1.1 $T_{oc}$ for a range N=51–71 orders, (b) $\tau =$ 0.0, 0.1 $T_{oc}$, 0.2 $T_{oc}$ and 0.3 $T_{oc}$ for a range N=61–81 harmonic order, (c) $\tau =$ 1.7 $T_{oc}$, 1.8 $T_{oc}$, 1.9 $T_{oc}$, 2.0 $T_{oc}$ and 2.1 $T_{oc}$ over N=31–45 harmonic order.

Download Full Size | PDF

For each harmonic peak shown in Fig. 1(a) we have obtained the position of the maxima ($\sim$ centroid) and plotted it with respect to the delay as shown in Fig. 1(c) (in solid lines). One can see from Fig. 1(a)-(c): (1) clear harmonic comb structures separated by two harmonic orders, and (2) a widely tunable spectral shift of more than two harmonic orders for a change in delay of $0.2~T_{oc}$ as we vary the delay $\tau =0.9~T_{oc}$ to $\tau = 1.1~T_{oc}$. Thus the harmonic comb is spanning the complete spectral range ($2\omega$) in this cut-off region. Both the central frequency and bandwidths of the IR-IR driver pulse manifest in a rather complicated manner in the HHG spectra. In this article we will mostly focus on understanding the origin of the wide tunability of the spectral shift. The shift of the harmonics with respect to the delay can be understood in terms of two important contributing factors: (a) adiabatic contribution, arising from a variation in the central frequency of the composite IR-IR pulse, and (b) non-adiabatic contribution, related to the phase shift of the recolliding electron wavepacket due to the subcycle change in amplitude of the electric field. In the following we will discuss these two contributions in detail.

3.1 Adiabatic contribution

To quantify the variation in the central frequency of the composite driving pulse, the spectrum of the IR-IR pulse for a few sample delays are shown in Fig. 3(c). In Fig. 3(a) we show the spectrum of IR-IR pulses for a delay range of $\pm 0.2~T_{oc}$ around the delay position $\tau =1T_{oc}$ in color map. The peak frequency [frequency at which the intensity is maximum] and central frequency [defined as $\omega _{IR-IR}(\tau )=\int \omega |E_{IR-IR}^{\tau }(\omega )|^2d\omega /\int |E_{IR-IR}^{\tau }(\omega )|^2d\omega$] (see cross and the dashed lines in Fig. 3(a),(d),(e) on top of the color maps) are plotted on the color map. Note that, the central and peak frequency are nearly identical in the chosen parameter range, we have used the peak frequency for further analysis in this paper. For $\tau =1T_{oc}$ the central frequency is the same as that of the single IR pulse ($\omega =0.057~$au or $\lambda =800~$nm). The continuous shift of the central frequency of the IR-IR pulse versus the the IR-IR delay near $\tau =1T_{oc}$ can be clearly seen here. For $\tau < T_{oc}$ the central frequency shifts towards the blue side (higher than $\omega =0.057~au$) and for $\tau > T_{oc}$ the central frequency shifts towards red side (lower than $\omega =0.057~au$). Intuitively, the time delay between the two IR pulses corresponds to a phase-shift between their spectral components in the Fourier-domain. The interference effect in the spectral domain modifies the resulting spectrum of the composite pulse leading to a shift in the central frequency [27]. This frequency shift of the IR-IR composite driving pulse with respect to the delay clearly explains the main contribution to the shift of the harmonic comb near-cutoff. In Fig. 1(c) we have plotted the peak of the harmonic bands as marked in Fig. 1(a) along with the harmonic frequency position of the driver IR-IR pulse with dashed lines. Quantitatively, the shift in frequency of the driver from $\tau =0.9T_{oc}$ to $\tau =1.1T_{oc}$ is $\delta \omega =\sim 0.0025~au$ which translates to a much higher shift in the harmonic orders. For example, for $n=58$ order of the fundamental IR ($\omega =0.057 ~au$) the shift is $n*\delta \omega \sim 2.5$, as observed in Fig. 1(c). For the soft-core Coulomb potential the maxima positions of the harmonic bands along with the $n$ times the central frequency of the driver IR-IR pulse are shown in Fig. 1(f) in solid and dashed lines, respectively. The results for the Gaussian and soft-Coulomb potentials are qualitatively same.

 figure: Fig. 3.

Fig. 3. (a) Frequency spectra vs delay of the driving IR-IR pulse near the first maxima (i.e., around $\tau =1 T_{oc}$). (b) Electric-field profile of the IR-IR pulse for different delays $0.8T_{oc}, 0.9T_{oc}, 1.0T_{oc}, 1.1T_{oc}$ and $1.2T_{oc}$, and (c) corresponding spectra. (d),(e) Same as (a) near $\tau =0$ and around $\tau =2 T_{oc}$. Central (peak) frequency of the delayed pulses are marked as dashed line (cross) on the color maps.

Download Full Size | PDF

A similar mechanism can be observed for the near zero delay position i.e., for $\tau =0$ and for $\tau =2T_{oc}$ as shown in Fig. 2(b),(c), respectively. In Fig. 3(d),(e) we have presented the spectrum of IR-IR pulses for a delay range of $0.2~oc$ near the central delay position $\tau =0$ and for $\tau =2T_{oc}$, respectively, along with the peak and central frequencies. Clearly for $\tau =0$ the delay dependent shift of the driving IR-IR pulse is smaller, so as the shift of the harmonics (see Fig. 2(b)). On the other hand the central frequency shift of the driving IR-IR pulse is higher with respect to the delay near $\tau =2T_{oc}$ as shown in Fig. 3(e) which is reflected as pronounced shifts of the generated harmonics (see Fig. 2(c)).

Thus we see that the modulation of the central frequency of the driving IR-IR pulse partially explains the total frequency tuning of the generated high harmonics. In Fig. 1(d) we have plotted the difference of the peak/central frequencies of the generated harmonic comb with respect to the harmonics of the driver along the harmonic bands (for example, as marked by 1, 2 along the representative dashed lines in Fig. 1(a),(b),(e)). It is clearly evident from this graph that there is a second important contributing factor to the frequency tuning of the generated harmonics. Specifically (see in Fig. 4(d)-(f)) (i) there is a red-shift of the harmonics throughout the generating range near $1T_{oc}$, (ii) the red shift increases with increasing the offset ($\Delta \tau =|\tau -T_{oc}|$) delay from $\tau =1T_{oc}$.

 figure: Fig. 4.

Fig. 4. (a)–(c) The peak positions of the harmonic bands in the adiabatic driving limit (see text for detail). The corresponding plots for the original driving are shown in (d)–(f).

Download Full Size | PDF

3.2 Non-adiabatic contribution

To understand the origin of the non-adiabatic shift of the harmonic spectra, we have computed the harmonic spectra with respect to different delay in the adiabatic limit. To this end we have designed the driving pulse with a trapezoidal envelope of 12oc length with 1oc rising time from zero to maximum amplitude, 10oc constant amplitude and 1oc falling time back to zero. We take the frequency of the driver as $\omega _{driver}(\tau )=\omega _{IR-IR}^{center}(\tau )$ and the intensity of the driver as $I_{driver}(\tau )=I_{IR-IR}^{center}(\tau )$. Therefore, now we drive the atomic system with the modified frequency, i.e., the central frequency of the composite IR-IR driver at given $\tau$ value, along with the reduced contribution of the non-adiabatic effect due to the change in intensity over the pulse length. We have extracted the central (peak) position of the harmonic bands (as explained earlier) over the range of delays and plotted along with the $n^{th}$ order of the driving central frequency ($\omega _{IR-IR}^{center}(\tau )$) in Fig. 4. For a comparison we have plotted the same for the harmonics generated with the original IR-IR driving. Clearly the shift of the harmonic band peaks (centers) other than that due to the shift of the driving pulse’s frequency is reduced near zero for $0$, $1T_{oc}$ and $2T_{oc}$ ((a)-(c) of Fig. 4) in the adiabatic limit. This explains that the red (blue in some cases) shifts are due to the non-adiabatic dynamics of HHG in the short pulse limit.

Further, the temporal origin of the red-shift can be identified from the time-frequency analysis of the harmonic generation. The time-frequency response of of the near-cutoff harmonics generated at different delays around $1T_{oc}$ are presented in Fig. 5. The pulse centers vary from $t= 2.9-3.1 ~T_{oc}$ in this delay range from $0.8-1.2 ~T_{oc}$ (see the vertical dashed line at $3.0T_{oc}$) From the time-frequency response (TFR) it is clear that the harmonics are dominantly generated in the “falling edge” of the IR-IR pulse. In this scenario, the major part of the contributing return electron trajectories experience a drop in the electric field within its return (with respect to the conventional flat pulse) time, $\sim T/2$. This leads to an additional phase shift of the return electron amplitude which in turn leads to the shift of the harmonic peaks towards the red side [10,18].

 figure: Fig. 5.

Fig. 5. Time-frequency responses of the harmonic generation for different delays are shown in color map (log scale). Color scales are given on the right-side of each map. Corresponding driving electric fields are shown in dashed lines in each plot (electric field scales are within $\pm 0.15$ au in y-axis).

Download Full Size | PDF

3.3 Tuning of IAP with a near single-cycle driver

We further show that the proposed scheme of double IR pulses can produce tunable isolated attosecond pulses in the regime of single-cycle carriers. To this end we have presented the spectra-map in Fig. 6(a), where each IR pulse is of 1.5 optical cycles long and the carrier-envelope phase is kept same as before (i.e., $\pi /2$).

 figure: Fig. 6.

Fig. 6. (a) Spectra-map for IR-IR driving obtained through TDSE calculation with Gaussian potential. Each IR pulse is 1.5 optical cycles long with other parameters kept same as in Eq. (1). In the inset of (a) IR-IR driving pulses for different delays are presented. (b) The XUV pulse profiles for the delay range $0-0.35$ optical cycles are presented as color map (log scale). Color scales are given on the right side of each map.

Download Full Size | PDF

The spectral tunability of the cutoff regime of the harmonic spectra-map is clearly visible (denoted as “abcdef” for guiding the eye) which ranges from a continuum of $\sim 15$ orders near zero-delay $\tau =0$ to $\sim 5$ orders near $\tau =0.35T_{oc}$ spanning over 20 harmonic orders. We have computed the corresponding temporal profiles of the pulses generated in this spectral range. In Fig. 6(b) the color map of these temporal profiles in the delay range is shown. The temporal width of these IAP changes from $\sim 0.005T_{oc} (13~as)$ to $\sim 0.02T_{oc} (53~as)$ in this delay range.

3.4 Variation of tunability with respect to pulse length

A quantification of the tuning range with respect to pulse length is very important. To this end we have computed the variation of the tuning range vs delay for several pulse durations from 40 cycle to 3 optical cycle. We have shown the spectral peak positions of the harmonic band corresponding to $39^{th}$ harmonic order near $1T_{oc}$ delay between the two IR pulses using both soft-core Coulomb and Gaussian potentials in Fig. 7(a). In order to make a comparison of the tuning range of near-cutoff harmonics, we have kept the cut-off harmonic nearly identical, for which the pulse amplitude for the 3 cycle pulse was increased appropriately. To quantify, we have computed the tuning range (defined as the change in harmonic peak position for a delay range, $\Delta \tau =0.3T_{oc}$ about delay position $1T_{oc}$, marked as $\Delta \Omega$ in Fig. 7(a) for a five-cycle pulse) and the tuning rate (defined as $\frac {1}{\Omega _{mean}}\Delta \Omega /\Delta \tau$, where $\Omega _{mean}$ is the mean position of the maxima of the harmonic band in the range as shown in Fig. 7(a)) as shown in Fig. 7(b),(c), respectively. Comparison of tuning range and tuning rate of a representative 39 order harmonics shows that tuning range increases as one decreases the pulse length. It is also clear from the plots that the tuning range and tuning rate are nearly the same for both Gaussian and soft-core Coulomb potentials, although the detailed structure of HHG may vary for the two potentials [38]. The range of tuning varies from $\sim 0.04$ harmonic order for long pulses ($\sim 40$ optical cycles long) to $\sim 7$ harmonic order near 3-cycle pulse with the corresponding tuning rate from $\sim 0.003-\sim 0.6$ per optical cycles, respectively.

4. Conclusions

We have theoretically shown that the time-delayed identical IR pulses leads to a wide spectral tuning of atomic HHG in excess of $\sim 2\omega$ frequency range, for several near-cutoff harmonics. We found that the observed spectral shift results from two contributions: (a) from the shift in the central frequency of the driving IR-IR pulse with varying delay, and (b) from to the non-adiabatic response of the recolliding electron wavepacket to the change in the laser amplitude within an optical cycle. Both, the adiabatic and nonadiabatic component, depend on the delay. For single-cycle pulses, this double-pulse scheme could be used as a source for tunable IAP whereas for the few-cycle regime this results in tunable APT. Finally we have quantified the tuning range and tuning rate with respect to pulse duration. Such tunability of harmonics would be highly useful in application of attosecond pulses and pulse trains to probe atomic and molecular dynamics. Of course, one has to remember that these results have been obtained within a single atom picture. The macroscopic propagation and phase matching can alter the experimental result. It would be interesting to explore how the IR-IR pulse-shaping scheme controls the generation of lower-order harmonics for various driving frequencies. We envision that the work will find applications in various time-resolved spectroscopic studies including tunable seeding of XFEL sources, tunable RABBITT and will motivate experiments to generate widely tunable HHG.

 figure: Fig. 7.

Fig. 7. (a) Peak position of the 39 harmonic order band vs IR-IR time delay near $1T_{oc}$ delay. Filled circles and corresponding fits in dashed line (soft-core Coulomb) and open squares and fit (Gaussian potential). Pulse durations from 40 to 3 optical cycles are marked on the plot. (b), (c) Tuning range and tuning rate for several pulse durations for Gaussian and soft-core Coulomb potentials, as extracted from the data in (a). Shaded regions show sub-5 optical cycle regime.

Download Full Size | PDF

Funding

Science and Engineering Research Board.

Acknowledgments

AM would like to acknowledge the IISER Mohali post doctoral fellowship and SERB post doctoral fellowship. KPS acknowledges Max-Planck Society for funding through Partners group and DST India for SERB/CRG grant.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper may be obtained from the authors upon reasonable request.

References

1. P. Wernet, J. Gaudin, K. Godehusen, O. Schwarzkopf, and W. Eberhardt, “Femtosecond time-resolved photoelectron spectroscopy with a vacuum-ultraviolet photon source based on laser high-order harmonic generation,” Rev. Sci. Instrum. 82(6), 063114 (2011). [CrossRef]  

2. C. Corder, P. Zhao, X. Li, M. D. Kershis, M. G. White, and T. K. Allison, “Development of a tunable high repetition rate XUV source for time-resolved photoemission studies of ultrafast dynamics at surfaces,” In Laser Applications in Microelectronic and Optoelectronic Manufacturing (LAMOM) XXIII, International Society for Optics and Photonics (SPIE, 2018), pp. 52–68.

3. G. L. Dakovski, Y. Li, T. Durakiewicz, and G. Rodriguez, “Tunable ultrafast extreme ultraviolet source for time-and angle-resolved photoemission spectroscopy,” Rev. Sci. Instrum. 81(7), 073108 (2010). [CrossRef]  

4. E. J. Sie, T. Rohwer, C. Lee, and N. Gedik, “Time-resolved XUV ARPES with tunable 24–33 eV laser pulses at 30 meV resolution,” Nat. Commun. 10(1), 3535 (2019). [CrossRef]  

5. W. Ubachs, “High resolution molecular physics studies using a laser-based tunable XUV source at 107 resolving power,” J. Electron Spectrosc. Relat. Phenom. 144-147, 97–101 (2005). [CrossRef]  

6. J. Duris, S. Li, T. Driver, et al., “Tunable isolated attosecond X-ray pulses with gigawatt peak power from a free-electron laser,” Nat. Photonics 14(1), 30–36 (2020). [CrossRef]  

7. C. Winterfeldt, C. Spielmann, and G. Gerber, “Colloquium: Optimal control of high-harmonic generation,” Rev. Mod. Phys. 80(1), 117–140 (2008). [CrossRef]  

8. M. Chini, K. Zhao, and Z. Chang, “The generation, characterization and applications of broadband isolated attosecond pulses,” Nat. Photonics 8(3), 178–186 (2014). [CrossRef]  

9. S. Yue, S. Fu, J. Li, X. Zhang, Y. Feng, B. Hu, and H. Du, “A redshift mechanism of high-order harmonics: Change of ionization energy,” J. Chem. Phys. 148(23), 234304 (2018). [CrossRef]  

10. X.-B. Bian and A. D. Bandrauk, “Nonadiabatic molecular high-order harmonic generation from polar molecules: Spectral redshift,” Phys. Rev. A 83(4), 041403 (2011). [CrossRef]  

11. G.-R. Jia, X.-H. Huang, and X.-B. Bian, “Nonadiabatic redshifts in high-order harmonic generation from solids,” Opt. Express 25(20), 23654–23662 (2017). [CrossRef]  

12. H. T. Kim, D. G. Lee, K.-H. Hong, J.-H. Kim, I. W. Choi, and C. H. Nam, “Continuously tunable high-order harmonics from atoms in an intense femtosecond laser field,” Phys. Rev. A 67(5), 051801 (2003). [CrossRef]  

13. F. Calegari, C. Vozzi, M. Negro, G. Sansone, F. Frassetto, L. Poletto, P. Villoresi, M. Nisoli, S. De Silvestri, and S. Stagira, “Efficient continuum generation exceeding 200 eV by intense ultrashort two-color driver,” Opt. Lett. 34(20), 3125–3127 (2009). [CrossRef]  

14. T. Siegel, R. Torres, D. J. Hoffmann, L. Brugnera, I. Procino, A. Zaïr, J. G. Underwood, E. Springate, I. C. E. Turcu, L. E. Chipperfield, and J. P. Marangos, “High harmonic emission from a superposition of multiple unrelated frequency fields,” Opt. Express 18(7), 6853–6862 (2010). [CrossRef]  

15. B. Schütte, P. Weber, K. Kovács, E. Balogh, B. Major, V. Tosa, S. Han, M. J. J. Vrakking, K. Varjú, and A. Rouzée, “Bright attosecond soft X-ray pulse trains by transient phase-matching in two-color high-order harmonic generation,” Opt. Express 23(26), 33947–33955 (2015). [CrossRef]  

16. B. Major, E. Balogh, K. Kovács, S. Han, B. Schátte, P. Weber, M. J. J. Vrakking, V. Tosa, A. Rouzée, and K. Varjú, “Spectral shifts and asymmetries in mid-infrared assisted high-order harmonic generation,” J. Opt. Soc. Am. B 35(4), A32–A38 (2018). [CrossRef]  

17. G. Li, Y. Zheng, X. Ge, Z. Zeng, and R. Li, “Frequency modulation of high-order harmonic generation in an orthogonally polarized two-color laser field,” Opt. Express 24(16), 18685–18694 (2016). [CrossRef]  

18. H. J. Shin, D. G. Lee, Y. H. Cha, K. H. Hong, and C. H. Nam, “Generation of nonadiabatic blueshift of high harmonics in an intense femtosecond laser field,” Phys. Rev. Lett. 83(13), 2544–2547 (1999). [CrossRef]  

19. C.-G. Wahlström, J. Larsson, A. Persson, T. Starczewski, S. Svanberg, P. Salieres, P. Balcou, and A. L’Huillier, “High-order harmonic generation in rare gases with an intense short-pulse laser,” Phys. Rev. A 48(6), 4709–4720 (1993). [CrossRef]  

20. M. V. Frolov, N. L. Manakov, A. M. Popov, O. V. Tikhonova, E. A. Volkova, A. A. Silaev, N. V. Vvedenskii, and A. F. Starace, “Analytic theory of high-order-harmonic generation by an intense few-cycle laser pulse,” Phys. Rev. A 85(3), 033416 (2012). [CrossRef]  

21. P. Raith, C. Ott, and T. Pfeifer, “Attosecond twin-pulse control by generalized kinetic heterodyne mixing,” Opt. Lett. 36(2), 283–285 (2011). [CrossRef]  

22. P. Raith, C. Ott, C. P. Anderson, A. Kaldun, K. Meyer, M. Laux, Y. Zhang, and T. Pfeifer, “Fractional high-order harmonic combs and energy tuning by attosecond-precision split-spectrum pulse control,” Appl. Phys. Lett. 100(12), 121104 (2012). [CrossRef]  

23. P. Raith, C. Ott, K. Meyer, A. Kaldun, M. Laux, M. Ceci, C. P. Anderson, and T. Pfeifer, “Carrier-envelope phase-and spectral control of fractional high-harmonic combs,” J. Appl. Phys. 114(17), 173102 (2013). [CrossRef]  

24. J. A. Pérez-Hernández, D. J. Hoffmann, A. Zaïr, L. E. Chipperfield, L. Plaja, C. Ruiz, J. P. Marangos, and L. Roso, “Extension of the cut-off in high-harmonic generation using two delayed pulses of the same colour,” J. Phys. B 42(13), 134004 (2009). [CrossRef]  

25. A. Mandal and K. P. Singh, “IR-IR control of High Harmonic Generation,” presented at Ultra-Fast Sciences 2019, IIT Bombay, India, 7-9 Nov. 2019

26. A. Mandal, S. Dahiya, A. Tyagi, M. S. Sidhu, and K. P. Singh, “Control of High Harmonic Generation in two delayed IR pulses,” presented at CAMNP 2019, Department of Applied Physics, Delhi Technological University, Delhi, India, 18-20 Dec. 2019.

27. L. G. Oldal, T. Csizmadia, P. Ye, N. G. Harshitha, A. Zaïr, S. Kahaly, K. Varjú, M. Fále, and B. Major, “Generation of high-order harmonics with tunable photon energy and spectral width using double pulses,” Phys. Rev. A 102(1), 013504 (2020). [CrossRef]  

28. S. Dahiya, M. S. Sidhu, A. Tyagi, A. Mandal, B. Nandy, J. M. Rost, T. Pfeifer, and K. P. Singh, “In-line ultra-thin attosecond delay line with direct absolute-zero delay reference and high stability,” Opt. Lett. 45(18), 5266–5269 (2020). [CrossRef]  

29. V. Schuster, V. Hilbert, R. Klas, C. Liu, M. Tschernajew, B. Bernhardt, J. Rothhardt, and J. Limpert, “Agile spectral tuning of high-order harmonics by interference of two driving pulses,” Opt. Express 29(14), 22117–22126 (2021). [CrossRef]  

30. L. G. Oldal, P. Ye, Z. Filus, T. Csizmadia, T. Grósz, M. de Marco, Z. Bengery, I. Seres, B. Gilicze, P. Jójárt, K. Varjú, S. Kahaly, and B. Major, “All-optical experimental control of high harmonic photon energy,” Phys. Rev. Appl. 16(1), L011001 (2021). [CrossRef]  

31. C. F. de M. Faria and J. M. Rost, “High-order harmonic generation from a confined atom,” Phys. Rev. A 62(5), 051402 (2000). [CrossRef]  

32. Q. Su and J. H. Eberly, “Model atom for multiphoton physics,” Phys. Rev. A 44(9), 5997–6008 (1991). [CrossRef]  

33. A. A. Silaev, M. Yu. Ryabikin, and N. V. Vvedenskii, “Strong-field phenomena caused by ultrashort laser pulses: Effective one- and two-dimensional quantum-mechanical descriptions,” Phys. Rev. A 82(3), 033416 (2010). [CrossRef]  

34. J. Itatani, J. Levesque, D. Zeidler, H. Pépin, N. Hiromichi, J. C. Kieffer, P. B. Corkum, and D. M. Villeneuve, “Tomographic imaging of molecular orbitals,” Nature (London) 432(7019), 867–871 (2004). [CrossRef]  

35. J. Itatani, D. Zeidler, J. Levesque, M. Spanner, D. M. Villeneuve, and P. B. Corkum, “Controlling High Harmonic Generation with Molecular Wave Packets,” Phys. Rev. Lett. 94(12), 123902 (2005). [CrossRef]  

36. M. Lewenstein, P. Balcou, M. Y. Ivanov, A. L’Huillier, and P. B. Corkum, “Theory of high-harmonic generation by low-frequency laser fields,” Phys. Rev. A 49(3), 2117–2132 (1994). [CrossRef]  

37. W. C. Chu, J. C. Travers, and P. St J Russell, “Near-ionization-threshold emission in atomic gases driven by intense sub-cycle pulses,” New J. Phys. 18(2), 023018 (2016). [CrossRef]  

38. S. V. Popruzhenko, “Coulomb phase in high harmonic generation,” J. Phys. B: At., Mol. Opt. Phys. 51(14), 144006 (2018). [CrossRef]  

Data availability

Data underlying the results presented in this paper may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. The harmonic spectra-map for a delay range $\tau = 0.8 - 1.2 T_{oc}$. (a) HHG for the Gaussian potential, and (b) SFA calculation. (c) Peak positions of the harmonic bands (solid line) vs delay, obtained from (a). The dashed lines represent delay-dependent spectral position of $N^{th}$ harmonic with respect to the central frequency of the driver. (d) The difference between the calculated maxima positions and that of the $N^{th}$ multiple of the central frequency of IR-IR driver. (e) and (f) are same as (a) and (c) for soft-core Coulomb potential.
Fig. 2.
Fig. 2. The harmonic spectra (a) for different delays $\tau =$ 0.9 $T_{oc}$, 1.0 $T_{oc}$ and 1.1 $T_{oc}$ for a range N=51–71 orders, (b) $\tau =$ 0.0, 0.1 $T_{oc}$, 0.2 $T_{oc}$ and 0.3 $T_{oc}$ for a range N=61–81 harmonic order, (c) $\tau =$ 1.7 $T_{oc}$, 1.8 $T_{oc}$, 1.9 $T_{oc}$, 2.0 $T_{oc}$ and 2.1 $T_{oc}$ over N=31–45 harmonic order.
Fig. 3.
Fig. 3. (a) Frequency spectra vs delay of the driving IR-IR pulse near the first maxima (i.e., around $\tau =1 T_{oc}$). (b) Electric-field profile of the IR-IR pulse for different delays $0.8T_{oc}, 0.9T_{oc}, 1.0T_{oc}, 1.1T_{oc}$ and $1.2T_{oc}$, and (c) corresponding spectra. (d),(e) Same as (a) near $\tau =0$ and around $\tau =2 T_{oc}$. Central (peak) frequency of the delayed pulses are marked as dashed line (cross) on the color maps.
Fig. 4.
Fig. 4. (a)–(c) The peak positions of the harmonic bands in the adiabatic driving limit (see text for detail). The corresponding plots for the original driving are shown in (d)–(f).
Fig. 5.
Fig. 5. Time-frequency responses of the harmonic generation for different delays are shown in color map (log scale). Color scales are given on the right-side of each map. Corresponding driving electric fields are shown in dashed lines in each plot (electric field scales are within $\pm 0.15$ au in y-axis).
Fig. 6.
Fig. 6. (a) Spectra-map for IR-IR driving obtained through TDSE calculation with Gaussian potential. Each IR pulse is 1.5 optical cycles long with other parameters kept same as in Eq. (1). In the inset of (a) IR-IR driving pulses for different delays are presented. (b) The XUV pulse profiles for the delay range $0-0.35$ optical cycles are presented as color map (log scale). Color scales are given on the right side of each map.
Fig. 7.
Fig. 7. (a) Peak position of the 39 harmonic order band vs IR-IR time delay near $1T_{oc}$ delay. Filled circles and corresponding fits in dashed line (soft-core Coulomb) and open squares and fit (Gaussian potential). Pulse durations from 40 to 3 optical cycles are marked on the plot. (b), (c) Tuning range and tuning rate for several pulse durations for Gaussian and soft-core Coulomb potentials, as extracted from the data in (a). Shaded regions show sub-5 optical cycle regime.

Equations (2)

Equations on this page are rendered with MathJax. Learn more.

i ψ ( x , t ) t = 1 2 2 x 2 ψ ( x , t ) + V ( x ) ψ ( x , t ) + H I ψ ( x , t ) ,
E ( t ) = E 01 f ( t ) cos ( ω t + ϕ c e ) + E 02 f ( t τ ) cos ( ω ( t τ ) + ϕ c e ) ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.