Modulation of multiphoton resonant high-order harmonic spectra driven by two frequency-comb fields
Zhao Yuanyuan, Zhao Di, Jiang Chen-Wei, Li Fu-li
School of Science, Xi’an Jiaotong University, Xi’an 710049, China

 

† Corresponding author. E-mail: d.zhao@mail.xjtu.edu.cn

Abstract

The frequency-comb structure in the extreme ultraviolet (XUV) and vacuum ultraviolet (VUV) regions can be realized by the high-order harmonic generation (HHG) process driven by frequency-comb fields, providing an alternative approach for the measurement of an unknown frequency in XUV or VUV. We consider the case of two driving frequency-comb fields with the same repetition frequency and the carrier frequencies of fundamental- and third-harmonics, respectively. The many-mode Floquet theorem (MMFT) is employed to provide a nonperturbative and exact treatment of the interaction between a quantum system and the frequency-comb laser fields. Multiphoton transition paths involving both fundamental- and third-harmonic photons are opened due to the coupling of the third-harmonic frequency-comb field. The multiphoton transition paths are superpositioned when the carrier-envelope-phase shifts (CEPs) fulfill the matching condition. And the interference of the multiphoton transition paths can be controlled by tuning the relative envelope delay between the fields. We find that the quasienergy structure, as well as the multiphoton resonant high-order harmonic generation (HHG) spectra, driven by the two frequency-comb fields can be coherently controlled via the interference of multiphoton transition paths. It is also found that the spectral intensities of the generated harmonics can be modulated, and the modulation behavior is harmonic-sensitive.

1. Introduction

In the past two decades, high-order harmonic generation (HHG) from the interaction between matter and intense laser fields has been proposed as an alternative method to obtain coherent light sources in the extreme ultraviolet (XUV) and vacuum ultraviolet (VUV) regions.[14] Besides the fundamental physics of the laser–matter interaction, the manipulate of the features of the generated spectra has also drown researchers’ attention due to potential applications.[510] When the matter–laser interaction is driven by a train of identical laser pulses with a fixed time interval, which corresponds to a frequency-comb structure in the frequency-domain, the frequency-comb structure is preserved in the harmonics in the HHG spectrum.[1114] It provides an approach to bring the frequency-comb structure into XUV and VUV regions, which will allow the measurement of an unknown frequency in XUV or VUV owing to the precise and direct link between known and unknown frequencies. The generation of XUV and VUV frequency comb spectroscopy has led to the advancements in development of optical atomic clock,[1520] high-precision optical frequency spectroscopy,[2128] and test of quantum electrodynamics.[14,2931]

It has been theoretically proved that HHG driven by a frequency-comb field can be dramatically enhanced via multiphoton resonance by tuning the carrier-envelope phase (CEP) shift of the driving field.[32,33] The employment of different color laser pulses has also been proved to be one of the schemes to enhance the intensities of the HHG spectra.[3448] And the spectra can be modulated by tuning the relative envelope delay τd and/or the relative phase between pulses φ. For the case of two multiple-cycle pulses, the modulation features of the spectra by tuning τd are similar to those by tuning φ. When two trains of pulses are employed, the modulation features of the spectra by tuning τd may be different from those by tuning φ. Recently, a two-level system driven by intense fundamental- and weak second-harmonic frequency-comb fields was investigated in Ref. [49]. It was shown that the frequency-comb HHG spectra can be enhanced via constructively interference between the multiphoton transition paths.

In this paper, we theoretically investigate the coherent control of the quasienergy structure and multiphoton resonant high-order harmonic generation (HHG) driven by two frequency-comb fields. Additional multiphoton transition paths involving photons from both fields are opened, and the coherent control is achieved via interference of the multiphoton transition paths by tuning the relative envelope delay between the fields. We consider the case of two frequency-comb fields with the same repetition frequency and the carrier frequencies of fundamental- and third-harmonics, respectively. When the CEP shifts of the two fields are matched, different multiphoton transition paths can interfere with each other. Our calculation indicates that the quasienergy structure of the driven system can be modulated by tuning the relative envelope delay. Meanwhile, the spectral intensities of the generated harmonics can be modulated, and the modulation behavior is harmonic-dependent, and related to the multiphoton transition paths that contribute to the spectrum.

The paper is organized as follows. We present the many-mode Floquet theorem (MMFT) for the treatment of the interaction between quantum systems and two frequency-comb laser fields with different carrier frequencies and the same repetition frequency in Section 2. The coherent control of the quasienergy structure and HHG power spectra is investigated in Section 3. And the conclusion is presented in Section 4.

2. Theoretical method

For a quantum system driven by laser fields, the time-dependent Schrödinger equation can be written as

with the Hamiltonian
where is the unperturbed Hamiltonian of the quantum system, and is the electric dipole moment operator. For the case of two driving frequency-comb fields, the laser field , where[32]
with the relative carrier phase φ and the relative envelop delay τd between the two fields. The comb frequencies are given as
with the main angular frequencies
where [ ] is the round function, ωc and are the carrier frequencies, ωr is the repetition frequency, and are the offset frequencies, and and are the CEP shifts of the two fields, respectively. The field amplitudes can be expressed as
with the peak amplitudes f0 and , and the standard deviation of a Gaussian function σ.

The time-dependent Schrödinger equation with Hamiltonian (2) is polychromatic, and it can be solved by employing the MMFT.[5054] For the case of fields with the same repetition frequency, there are three independent frequencies ω0, Ω0, and ωr. And only two independent frequencies (ω0 and ωr) are left when the specific matching condition for the two CEP shifts is fulfilled.[49] Then the time-dependent Schrödinger equation can be converted into an equivalent two-mode Floquet matrix eigenvalue problem. And the two-mode Floquet matrix eigenvalue equation can be constructed as[49]

where the basis vectors are employed. The two-mode Floquet matrix can be constructed by
with
where , , , and . More details of the structure of can be found in Refs. [32] and [49].

Solving the matrix eigenvalue problem (9) with the Floquet matrix (10), we can obtain a set of quasienergies and the corresponding eigenvectors which satisfy the orthonormality condition. And the harmonic generation spectra can be expressed as

in which each harmonic element .

3. Results and discussion

In this section, we present a case study of the modulation of quasienergy structures and HHG power spectra generated from a two-level system driven by fundamental- and third-harmonic frequency-comb laser fields by varying the envelope delay. The parameters of the fundamental frequency-comb field are peak intensity , carrier frequency 374.7 THz (corresponding to a.u. and wavelength 800 nm), and repetition frequency 10 THz (corresponding to a.u.) generated from a train of Gaussian pulses with 20 fs full width at half maximum (FWHM). The third-harmonic frequency-comb field has peak intensity , which is 4% of the fundamental field, and its carrier frequency is 1124.1 THz (corresponding to a.u.), i.e., the third harmonic of the fundamental field. The third-harmonic field is generated from a train of Gaussian pulses with 20 fs FWHM, which has repetition frequency 10 THz. The CEP shifts are set as , such that the matching condition is fulfilled.[49] For the two-level system, the energy separation between the levels is a.u., corresponding to the five-photon dominant resonance regime of the fundamental field ( ), and the transition dipole moment a.u. is used.

Under the matching condition , we calculate the quasienergies and time-averaged transition probabilities as a function of when the delay is set as , which are presented in Figs. 1(a) and 1(b), respectively. The resonance can be achieved five times by varying the CEP shift, and these five peaks are separated exactly by . The resonance peaks of different multiphoton transition paths are coincident due to the superposition of the paths.[49] In Fig. 1(c), we also present the time-averaged transition probabilities as a function of for cases with different delays , , ( ). The resonance peak position changes as τd is varied, i.e., the quasienergy structure is modified. The modification of the quasienergy structure is due to the change of the overlap between pulses from the two laser fields, which mathematically is embodied in the phase of the off-diagonal elements with and in the Floquet matrix (10). Physically, when the pulses from the fundamental field overlap with those from the third-harmonic field, the total amplitude of the driving fields is modulated by varying τd. Thus, the Stark energy shifts of the levels, i.e., the quasienergy structure, are altered. It is reasonable to expect that the modification of the quasienergy structure is large when varying τd in the vicinity of 0, i.e., when the pulses from the two fields mostly overlap with each other. And the quasienergy structure is nearly unchanged with large values of τd where the pulses are hardly overlap with each other. We plot the CEP shifts of the resonance peaks as a function of τd in Fig. 2. It is shown that the resonance peak position oscillates with a period about . And over multiple oscillation periods, the oscillation amplitude decreases with increasing τd.

Fig. 1. (color online) (a) Plots of quasienergies and (b) time-averaged transition probabilities as a function of the CEP shift for the two-level system driven by fundamental and third-harmonic frequency-comb fields with their relative envelope delay . It has five resonance positions separated by due to the five-photon dominant resonance. (c) Enlarged plots of the time-averaged transition probabilities as a function of the CEP shift, with envelope delay (black solid line), (red dash line), and (blue dotted line). The fundamental field has peak intensity , carrier frequency 374.7 THz, and repetition frequency 10 THz, and is generated from a train of Gaussian pulses with 20 fs FWHM. The third-harmonic field has peak intensity , carrier frequency 1124.1 THz, and repetition frequency 10 THz, and is generated from another train of Gaussian pulses with 20 fs FWHM. is set to the fulfill the matching condition. The relative carrier phase is set as φ = 0.
Fig. 2. (color online) The resonant CEP shifts as a function of the envelope delay τd from 0 to . at are represented with red circles and those at are represented with black squares. The resonant CEP shifts as a function of τd from 0 to are shown in the inset. The other parameters are the same as those in Fig. 1.

Besides the quasienergy structure, the harmonic spectra can also be modulated by varying the envelope delay. Owing to the coupling of the third-harmonic driving field, additional multiphoton transition paths involving both fundamental- and third-harmonic photons are opened, and they contribute to the HHG power spectra, as shown in Fig. 3(a). Figure 3(a) and 3(b) present the HHG power spectra driven by the frequency-comb laser fields with varying peak intensity of the third-harmonic driving field at and , respectively. Note that, for different values of and τd, the CEP shifts are all tuned to achieve the multiphoton resonances. One can see that the spectral intensities of the harmonics are enhanced with increasing . Generally, it is due to the strengthen of the multiphoton transition paths involving both fundamental- and third-harmonic photons. However, the enhancing feature of the case with is quite different from that with . For the case of , the enhancement factors of specific harmonics 11th and 17th are larger than those of the others, and multiple plateaus are formed in the spectrum. It suggests that the overlap in the time domain between pulses from the two driving fields impacts the harmonic spectra in the frequency domain, and this should be related to the opened multiphoton-transition paths.

Fig. 3. (color online) HHG power spectra with (a) and (b) by varying the peak intensity of the third-harmonic frequency-comb field. For the achievement of multiphoton resonance, (a) the CEP shift is tuned to for (black solid line), 0.101534 for (red dash line), and 0.113857 for (blue dotted line), and (b) the CEP shift is tuned to for (black solid line), 0.101703 for (red dash line), and 0.114485 for (blue dotted line). For comparison, the spectrum with is presented wit the black dash-dotted line in (a). The other parameters are the same as those in Fig. 1.

For further investigating the dependence of the spectra on the envelope delay, the HHG power spectra by varying the envelope delays τd with and are calculated and presented in Figs. 4(a) and 4(b), respectively. For each case, the CEP shift is tuned to achieve multiphoton resonances. One can see that, as τd is varied from to , the spectral intensities of the 11th and 17th harmonics are nearly unchanged, while those of the 9th, 13th, 15th, and 19th harmonics are suppressed. Considering the multiphoton transition paths under the multiphoton resonant condition, the electron is firstly promoted to the excited level via multiphoton resonant transition. After that, it absorbs additional multiple photons and transits back to the ground level, emitting a harmonic photon. The paths involving both fundamental- and third-harmonic photons are dependent on the envelope delay, since the multiphoton transition probabilities are determined by the field intensities in the overlap parts.

Fig. 4. (color online) HHG power spectra with (a) and (b) by varying the envelope delay τd. In each plot, the spectra with , , and are presented by black dash-dotted line, red dash line, and blue solid line, respectively. In panel (a), the CEP shift is tuned to for , 0.101645 for , and 0.101703 for for the achievement of multiphoton resonance. In panel (b), the CEP shift is tuned to for , 0.114277 for , and 0.114485 for for the achievement of multiphoton resonance. The other parameters are the same as those in Fig. 1.

For instance, the 11th harmonic is mainly contributed by three multiphoton-transition paths taking the excited level as the intermediate state, the first one via the absorption of 6 fundamental-harmonic photons, the second one via the absorption of 2 third-harmonic photons, and the third one via the absorption of 3 fundamental-harmonic photons and 1 third-harmonic photon. The first two paths are the main contribution and they are independent of the envelope delay, while the third path is dependent on the envelope delay. The 13th harmonic is also mainly contributed by three paths, the first one via absorption of 8 fundamental-harmonic photons, the second one via the absorption of 5 fundamental-harmonic photons and 1 third-harmonic photon, and the third one via the absorption of 2 fundamental-harmonic photons and 2 third-harmonic photons. However, only the first path is independent of the envelope delay. As a result, the spectral intensities of the 11th and 13th harmonics show different modulation behaviors as varying the envelope delay.

4. Conclusion

We present a theoretical investigation of the coherent control of the quasienergy structure and HHG power spectra driven by two frequency-comb fields via interference of multiphoton transition paths by tuning the relative envelope delay between the fields. We consider the case of two frequency-comb fields with the same repetition frequency and the carrier frequencies of fundamental- and third-harmonic frequencies. For accurate treatment of the interaction between a quantum system and the frequency-comb fields, the many-mode Floquet theory is employed. Multiphoton transition paths involving both fundamental- and third-harmonic photons are opened, and they can interfere with each other when the CEP shifts satisfy the matching condition. It is found that, not only the HHG power spectra, but also the quasienergy structure of the driven system can be coherently controlled. Our calculation results show that the quasienergies of the driven system can be tuned by varying the relative envelope delay when the pulses of the two fields overlap with each other. Meanwhile, owing to the dependence of the multiphoton transition paths on the relative envelope delay, the spectral intensities of the generated harmonics show different modulation behaviors when varying the relative envelope delay. Note that the effects of ionization and decay of levels, which are not included in this study of closed two-level system, also impact the quasienergy structure. And the study of realistic atomic systems including these effects is in process.

Reference
[1] L’Huillier A Balcou P 1993 Phys. Rev. Lett. 70 774
[2] Ravasio A Gauthier D Maia F R N C et al. 2009 Phys. Rev. Lett. 103 028104
[3] Winterfeldt C Spielmann C Gerber G 2008 Rev. Mod. Phys. 80 117
[4] Kapteyn H C Murnane M M Christov I P 2005 Phys. Today 58 39
[5] Hentschel M 2001 Nature 414 509
[6] Tzallas P 2003 Nature 426 267
[7] Kienberger R 2004 Nature 427 817
[8] Nabekawa Y Hasegawa H Takahashi E J Midorikawa K 2005 Phys. Rev. Lett. 94 043001
[9] Sansone G 2006 Science 314 443
[10] Goulielmakis E 2008 Science 320 1614
[11] Carrera J J Son S K Chu S I 2008 Phys. Rev. 77 031401
[12] Carrera J J Chu S I 2009 Phys. Rev. 79 063410
[13] Yost D C Schibli T R Ye J Tate J L Hostetter J Gaarde M B Schafer K J 2009 Nat. Phys. 5 815
[14] Kandula D Z Gohle C Pinkert T J Ubachs W Eikema K S E 2010 Phys. Rev. Lett. 105 063001
[15] Peik E Tamm Chr 2003 Europhys. Lett. 61 181
[16] Takamoto M Hong F L Higashi R Katori H 2005 Nature 435 321
[17] Gerginov V Tanner C E Diddams S A Bartels A Hollberg L 2005 Opt. Lett. 30 1734
[18] Rellergert W G DeMille D Greco R R Hehlen M P Torgerson J R Hudson E R 2010 Phys. Rev. Lett. 104 200802
[19] Campbell C J Radnaev A G Kuzmich A 2011 Phys. Rev. Lett. 106 223001
[20] Diddams S A Udem Th Bergquist J C Curtis E A Drullinger R E Hollberg L Itano W M Lee W D Oates C W Vogel K R Wineland D J 2001 Science 293 825
[21] Stowe M C Cruz F C Marian A Ye J 2006 Phys. Rev. Lett. 96 153001
[22] Witte S Zinkstok R Th Ubachs W Hogervorst W Eikema K S E 2005 Science 307 400
[23] Cingöz A Yost D C Allison Th K Ruehl A Fermann M E Hartl I Ye J 2012 Nature 482 68
[24] Stalnaker J E Coq Y Le Fortier T M Diddams S A Oates C W Hollberg L 2007 Phys. Rev. 75 040502
[25] Holzwarth R Udem Th Hänsch T W Knight J C Wadsworth W J Russell P St J 2000 Phys. Rev. Lett. 85 2264
[26] Niering M Holzwarth R Reichert J Pokasov P Udem Th Weitz M Hänsch T W Lemonde P Santarelli G Abgrall M Laurent P Salomon C Clairon A 2000 Phys. Rev. Lett. 84 5496
[27] Marian A Stowe M C Felinto D Ye J 2005 Phys. Rev. Lett. 95 023001
[28] Margolis H S Barwood G P Huang G Klein H Lea S N Szymaniec K Gill P 2004 Science 306 1355
[29] Herrmann M 2009 Phys. Rev. 79 052505
[30] Mabuchi H Ye J Kimble H J 1999 Appl. Phys. 68 1095
[31] Eyler E E 2008 Euro Phys. J. 48 43
[32] Son S K Chu S I 2008 Phys. Rev. 77 063406
[33] Zhao D Li F L Chu S I 2013 Phys. Rev. 87 043407
[34] Pfeifer T Gallmann L Abel M J Nagel P M Neumark D M Leone S R 2006 Phys. Rev. Lett. 97 163901
[35] Hong W Lu P Lan P Yang Z Li Y Liao Q 2008 Phys. Rev. 77 033410
[36] Zhang G T Wu J Xia C L Liu X S 2009 Phys. Rev. 80 055404
[37] Song X Zeng Z Fu Y Cai B Li R Cheng Y Xu Z 2007 Phys. Rev. 76 043830
[38] Wu J Zhang G T Xia C L Liu X S 2010 Phys. Rev. 82 013411
[39] Zeng Z Cheng Y Song X Li R Xu Z 2007 Phys. Rev. Lett. 98 203901
[40] Zhai Z Liu X S 2008 J. Phys. B: At. Mol. Opt. Phys. 41 125602
[41] Lan P Lu P Cao W Li Y Wang X 2007 Phys. Rev. 76 011402
[42] Faria C F Dörr M Becker W Sandner W 1999 Phys. Rev. 60 1377
[43] Li X X Xu Z Z Mao L Q 1997 Acta Phys. Sin. 6 197 in Chinese
[44] Ye X L Zhou X X Zhao S F Li P C 2009 Acta Phys. Sin. 58 1579 in Chinese
[45] Xiang Y Niu Y Gong S 2009 Phys. Rev. 79 053419
[46] Carrera J J Chu S I 2007 Phys. Rev. 75 033807
[47] Xu J Zeng B Yu Y 2010 Phys. Rev. 82 053822
[48] Xiao J Sun Z R Wang Y F Deng L Zhang X Y Chen G L Zhang W P Wang Z G Xu Z Z Li R X 2005 Chin. Phys. Lett. 22 873
[49] Zhao D Jiang C W Li F L 2015 Phys. Rev. 92 043413
[50] Chu S I Telnov D A 2004 Phys. Rep. 390 1
[51] Ho T S Chu S I Tietz J V 1983 Chem. Phys. Lett. 96 464
[52] Ho T S Chu S I 1984 J. Phys. 17 2101
[53] Ho T S Chu S I 1985 Phys. Rev. 31 659
[54] Ho T S Chu S I 1985 Phys. Rev. 32 377