Quantum path control and isolated attosecond pulse generation in the combination of near-infrared and terahertz pulses
Zhong Hui-Ying, Guo Jing†, Zhang Hong-Dan, Du Hui, Liu Xue-Shen‡
Institute of Atomic and Molecular Physics, Jilin University, Changchun 130012, China

Corresponding author. E-mail: gjing@jlu.edu.cn

Corresponding author. E-mail: liuxs@jlu.edu.cn

*Project supported by the National Natural Science Foundation of China (Grant Nos. 11174108, 11104108, and 11271158).

Abstract

We present an efficient and realizable scheme for the generation of an ultrashort single attosecond (as) pulse from H atom with a 800-nm fundamental laser field combined with a terahertz (THz) field. The high-order harmonic generation (HHG) can be obtained by solving the time-dependent Schrödinger equation accurately and efficiently with time-dependent generalized pseudo-spectral (TDGPS) method. The result shows that the plateau of high-order harmonics is extended and the broadband spectra can be produced by the combined laser pulse, which can be explained by the corresponding ionization probability. The time–frequency analysis and semi-classical three-step model are also presented to further investigate this mechanism. Besides, by the superposition of the harmonics near the cutoff region, an isolated 133-as pulse can be obtained.

PACS: 32.80.Rm; 42.65.Ky; 42.65.Re
Keyword: high-order harmonic generation; generalized pseudospectral method; attosecond pulse
1. Introduction

Attosecond pulses can trace the ultrafast electronic dynamics in atoms and molecules with high accuracy, which has been a forefront subject that attracts much attention in science and technology.[15] Since high harmonic spectrum presents plateau area which is distributed regularly with equal frequencies, which makes it become first choice to produce attosecond light.[68] A typical description of HHG spectrum presents a sharp decrease in the initial harmonics, then the spectrum shows a plateau where all the harmonics have a similar amplitude, and finally a cutoff at the energy Ip + 3.17Up. Here, Ip is the atomic ionization potential, and Up is the ponderomotive energy which only depends on the laser parameters. The essence of HHG process can be described by the semiclassical three-step model.[911] First, the free electrons tunnel from the ground state to continuous states by tunneling or multi-photon ionization. Second, the electrons are then accelerated by a strong oscillating laser field. Finally, the electrons recombine with the parent ion and emit photons. There are many techniques for the generation of attosecond pulses, such as scheme of a single-cycle pulse, [12] a few-cycle laser field, [1315] double optical gating, [16] and multicycle driving laser field directly from an amplifier.[17] By applying these techniques, isolated ultrashort XUV pulses have been obtained by several laboratories worldwidely.[18, 19] Recently, Luo et al.[20] have shown dramatic cutoff extension and broadband supercontinuum generation in multi-cycle two-color pulses, which can directly generate an isolated sub-100-as pulse. Hu et al.[21] proposed effect of electron correlation on high-order harmonic generation of He model atom. Li et al.[22] theoretically achieved an isolated 26-as pulse of He+ by means of an intense few-cycle chirped laser field. Recently, a scheme have been used to generate an isolated attosecond pulse by using a Terahertz (THz) field, which attracts much attention due to its potential applications.[2326]

In this paper, we will produce an isolated ultrashort attosecond laser pulse by controlling the short and long trajectories in a 800-nm fundamental laser field combined with a THz field. The high-order harmonic generation (HHG) can be obtained by solving the time-dependent Schrö dinger equation (TDSE) accurately with time-dependent generalized pseudospectral (TDGPS) method.[27, 28] The result shows that the plateau of high-order harmonics can be extended largely in the combined laser field. With such a spectrum, the generation of isolated attosecond pulses would be obtained straightforward, the shortest of which is 133  as.

2. Model

The TDSE for H atom can be written as (in atomic units)

Here, Ĥ 0 is the unperturbed Hamiltonian of H atom, and is the potential energy of H atom in an intense laser field, which is given by

The TDSE can be solved accurately and efficiently by the TDGPS method.[27, 28] This method has two key steps. (i) The spatial coordinates are optimally discretized by pseudospectral discretization, which allows us to obtain denser grids near the nuclear origin, leading to more accurate initial wavefunction by using smaller number of grid points. (ii) A second-order split-operator technique in the energy representation is used for the efficient and accurate three-step time propagation of the wave function, which is given by

Once the time-dependent wave function is determined, we can compute the expectation value of the induced dipole moment in the acceleration form:

The corresponding HHG power spectrum can be obtained by Fourier transformation of the time-dependent dipole moment

The attosecond pulse can be calculated by superposing several harmonics as follows:

where aq = d (t) e− iqω t dt. To study the detailed spectral and temporal structures of HHG, we perform time-frequency analysis by means of the wavelet transform,

where wt0, ω (t) is the wavelet transform kernel. It can be shown as follows:

where is the Morlet wavelet.

3. Numerical results

In our calculation, a 800-nm fundamental laser field only and a 800-nm fundamental laser field combined with a THz field can be expressed as follows:

Here, E1 = 0.10676  a.u. and ETHz = 0.01947  a.u. are the electric-field amplitudes, f1 (t) and fTHz (t) are the Gaussian envelopes, ω 1 = 0.0569  a.u. (λ = 800  nm) and ω THz = 0.00949  a.u. (62.5  THz) are the frequencies, φ 1 = 0.25π is the carrier-envelope phase (CEP), the durations of the Gaussian pulses are 5  fs and 16.022  fs, respectively. In our calculation, we use N = 1200 radial grid points and partial wave L = 120, the total time duration T = 16T0, where T0 = 2π /ω 1. Figure  1(a) shows E0(t) and Ec(t) as functions of time, respectively. Figure  1(b) presents the acceleration-form power spectra of H atom in fundamental field and combined field, respectively. From Fig.  1(b) we can see the HHG spectrum in the 800-nm fundamental laser field has a cutoff of 60th order. When a THz field is added, it can be apparently seen that the plateau in the combined fields is extended to about 95th order.

Fig.  1. (a) A 800-nm fundamental laser field only (black dotted curve) and a 800-nm fundamental laser field combined with a THz field (red solid curve). (b) The corresponding acceleration form power spectra in a 800-nm fundamental laser field only and a 800-nm fundamental laser field combined with a THz field, respectively.

Since HHG is very sensitive to ionization probability, the corresponding ionization probability in a 800-nm fundamental laser field only and a 800-nm fundamental laser field combined with a THz field are also calculated in Fig.  2. From Fig.  2 we can see that the ionization probability in both cases presents the same tendency: the curve increases rapidly at first, then almost keeps the same. While after adding the THz field, ionization occurs earlier, and ionization probability is also higher in combined laser field case, especially at T = − 0.67  o.c.

Fig.  2. Ionization probability in a 800-nm fundamental laser field only and a 800-nm fundamental laser field combined with a THz field, respectively.

Figure  3 shows the time– frequency distribution of the HHG ((a1), (a2)) and corresponding three-step models ((b1), (b2)) in a 800-nm fundamental laser field only and a 800-nm fundamental laser field combined with a THz field, respectively. In Fig.  3(a1), there are four main peaks contributing to the HHG, which are marked as P1, P2, P3, and P4. For harmonic orders above 35th, the peak of P4 almost makes no contribution to the HHG. Thus, only peaks P1, P2, and P3 are left for the HHG. The maximum harmonic orders of P1, P2, and P3 are about 62th, 53th, and 50th, respectively. The three peaks with both long and short trajectories (a trajectory with earlier ionization time but later emission time is called the long trajectory, and a trajectory with later ionization time but earlier emission time is called the short trajectory) all contribute to the HHG, and the intensity of long and short trajectories is comparable, leading to modulations in the harmonic spectrum, which is disadvantageous to synthesis single attosecond pulse. From Fig.  3(a2) we can see that there are also four peaks marked as and in the combined laser pulse case, the three peaks and are suppressed, and only the peak is left to contribute to the HHG for the harmonics from 60th to 95th order. For peak , two dominant quantum paths (short trajectory and long trajectory) with different emission times contribute to the same harmonic in half optical cycle of the fundamental pulse. It is apparently seen that the long trajectory is suppressed and the short one is relatively enhanced. Then for the harmonics in the plateau, the time of two paths are different and the harmonics are not in phase. Apparently, the contribution of short trajectory is higher than the long one, thus a single quantum path control is realized, which is advantageous to produce an isolated attosecond pulse.

Fig.  3. Time– frequency distribution of the HHG ((a1), (a2)) and corresponding three-step models ((b1), (b2)) in a 800-nm fundamental laser field only and a 800-nm fundamental laser field combined with a THz field, respectively.

Figure  3(b1) is the three-step model corresponding to a 800-nm fundamental laser field only case, and there are four major peaks of emission time marked as A, B, C, and D. The peaks A, B, and C all contribute to HHG. The maximum harmonic frequency of 60ω 1 is consistent with the cutoff position shown in Fig.  1(b) with the black dotted curve. In Fig.  3(b2), after adding a THz field in the 800-nm fundamental laser field, only the peak A′ is left to contribute to HHG for the harmonics over order of 60th, with the maximum harmonic frequency equal to 95ω 1. Thus, the cutoff of harmonics can be extended, which is consistent with the result shown by the red solid curve in Fig.  1(b). For peak A′ , electrons with the short trajectory are mainly ionized from T = − 0.67  o.c. to T = − 0.48  o.c., and electrons with the long trajectory are mainly ionized from T = − 0.80  o.c. to T = − 0.67  o.c. However, the ionization rate of electrons corresponding to the short trajectory is much higher than that corresponding to the long trajectory (see the red solid curve in Fig.  2). Consequently, the harmonic intensities for the short trajectory are much higher than that for the long trajectory, which is consistent with the time– frequency distribution (see Fig.  3(a2)).

Figure  4 shows the case of a 800-nm fundamental laser field, by superposing the HHG from the 40th to the 60th order, and a 252-as pulse and a relative weak 201-as pulse can be achieved. For the case of the combination of a 800-nm fundamental laser field combined with a THz field, the HHG spectrum of plateau is smoother and with much fewer modulation. Thus, by superposing the harmonics in the range from 75th to 95th order, an isolated attosecond pulse with duration of about 133 as can be obtained.

Fig.  4. Temporal profiles of the attosecond pulses generated from H atom: the black dotted curve describes the temporal profiles by superposing the harmonics from 40th to 60th order in a 800-nm fundamental laser field (5  fs); the red solid curve describes the temporal profiles by superposing the harmonics from 75th to 95th order in a 800-nm fundamental laser field combined with a THz field (5  fs/800 nm, 16  fs/4800 nm).

4. Conclusions

In conclusion, we theoretically investigate the high-order harmonic generation of H atom in a 800-nm fundamental laser field combined with a THz field by solving the time-dependent Schrö dinger equation accurately with TDGPS method. The result shows that the plateau of high-order harmonic is extended in the presence of THz field, which can be explained by the corresponding ionization probability. With the combined laser field, the long trajectory can be suppressed and the short trajectory can be enhanced. As a result, the optimizing combined field allows electrons to acquire more kinetic energy from the laser field, leading to the extension of the HHG plateau, and an isolated 133-as pulse can be generated.

Acknowledgments

The authors sincerely thank Prof. P. C. Li for helpful discussion.

Reference
1 Goulielmakis E, Loh Z H and Wirth Aet al. 2010 Nature 466 739 DOI:10.1038/nature09212 [Cited within:1]
2 Itatani J, Levesque J, Zeidler D, Hiromichi N, Pépin H, Kieffer J C, Corkum P B and Villeneuve D M 2004 Nature 432 867 DOI:10.1038/nature03183 [Cited within:1]
3 Paul P M, Toma E S, Breger P, Mullot G, Augé F, Balcou Ph, Muller H G and Agostini P 2001 Science 292 1689 DOI:10.1126/science.1059413 [Cited within:1]
4 Schafer K J, Gaarde M B, Heinrich A, Biegert J and Keller U 2004 Phys. Rev. Lett. 92 023003 DOI:10.1103/PhysRevLett.92.023003 [Cited within:1]
5 Johnsson P, López-Martens R and Kazamias Set al. 2005 Phys. Rev. Lett. 95 013001 DOI:10.1103/PhysRevLett.95.013001 [Cited within:1]
6 Mauritsson J, Johnsson P, Gustafsson E, L’Huillier A, Schafer K J and Gaarde M B 2006 Phys. Rev. Lett. 97 013001 DOI:10.1103/PhysRevLett.97.013001 [Cited within:1]
7 Chang Z 2004 Phys. Rev. A 70 043802 DOI:10.1103/PhysRevA.70.043802 [Cited within:1]
8 Carrera J J, Tong X M and Chu S I 2006 Phys. Rev. A 74 023404 DOI:10.1103/PhysRevA.74.023404 [Cited within:1]
9 Corkum P B 1993 Phys. Rev. Lett. 71 1994 DOI:10.1103/PhysRevLett.71.1994 [Cited within:1]
10 Lewenstein M, Balcou P, Ivanov M Y, L’Huillier A and Corkum P B 1994 Phys. Rev. A 49 2117 DOI:10.1103/PhysRevA.49.2117 [Cited within:1]
11 Krause J L, Schafer K J and Kuland er K C 1992 Phys. Rev. Lett. 68 3535 DOI:10.1103/PhysRevLett.68.3535 [Cited within:1]
12 Christov I P, Murnane M M and Kapteyn H C 1997 Phys. Rev. Lett. 78 1251 DOI:10.1103/PhysRevLett.78.1251 [Cited within:1]
13 Solanpää J, Budagosky J A, Shvetsov-Shilovski N I, Castro A, Rubio A and Räsänen E 2014 Phys. Rev. A 90 053402 DOI:10.1103/PhysRevA.90.053402 [Cited within:1]
14 Ge X L, Wang T, Guo J and Liu X S 2014 Phys. Rev. A 89 023424 DOI:10.1103/PhysRevA.89.023424 [Cited within:1]
15 Wu J, Zhang G T, Xia C L and Liu X S 2010 Phys. Rev. A 82 013411 DOI:10.1103/PhysRevA.82.013411 [Cited within:1]
16 Mashiko H, Gilbertson S, Li C, Khan S D, Shakya M M, Moon E and Chang Z H 2008 Phys. Rev. Lett. 100 103906 DOI:10.1103/PhysRevLett.100.103906 [Cited within:1]
17 Tang S S and Chen X F 2012 Phys. Rev. A 85 063816 DOI:10.1103/PhysRevA.85.063816 [Cited within:1]
18 Kaindl R A, Carnahan M A, Orenstein J, Chemla D S, Christen H M, Zhai H Y, Paranthaman M and Lowndes D 2001 Phys. Rev. Lett. 88 027003 DOI:10.1103/PhysRevLett.88.027003 [Cited within:1]
19 Zeng Z N, Cheng Y, Song X H, Li R X and Xu Z Z 2007 Phys. Rev. Lett. 98 203901 DOI:10.1103/PhysRevLett.98.203901 [Cited within:1]
20 Luo J H, Hong W Y, Zhang Q B, Liu K L and Lu P X 2012 Opt. Express 20 9801 DOI:10.1364/OE.20.009801 [Cited within:1]
21 Hu S L and Shi T Y 2013 Chin. Phys. B 22 013101 DOI:10.1088/1674-1056/22/1/013101 [Cited within:1]
22 Li P C, Zhou X X, Wang G L and Zhao Z X 2009 Phys. Rev. A 80 053825 DOI:10.1103/PhysRevA.80.053825 [Cited within:1]
23 Löffler T, Bauer T, Siebert K and Roskos H G 2001 Opt. Express 9 616 DOI:10.1364/OE.9.000616 [Cited within:1]
24 Ferguson B, Wang S, Gray D, Abbot D and Zhang X 2002 Opt. Lett. 27 13122 [Cited within:1]
25 Hong W Y, Lu P X, Lan P F, Zhang Q B and Wang X B 2009 Opt. Express 17 5139 DOI:10.1364/OE.17.005139 [Cited within:1]
26 Yuan K J and Band rauk A D 2013 Phys. Rev. Lett. 110 023003 DOI:10.1103/PhysRevLett.110.023003 [Cited within:1]
27 Tong X M and Chu S I 1997 Chem. Phys. 217 119 DOI:10.1016/S0301-0104(97)00063-3 [Cited within:2]
28 Li P C Laughlin C Chu S I 2014 Phys. Rev. A 89 023431 DOI:10.1103/PhysRevA.89.023431 [Cited within:2]