Attosecond pulse generation from two-electron harmonic emission spectrum*
Feng Li-Qianga),b)†, Liu Hangc)‡
College of Science, Liaoning University of Technology, Jinzhou 121000, China
State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
School of Chemical and Environmental Engineering, Liaoning University of Technology, Jinzhou 121000, China

Corresponding author. E-mail: lqfeng_lngy@126.com

Corresponding author. E-mail: liuhang20@126.com

*Project supported by the Scientific Research Fund of Liaoning Provincial Education Department, China (No. L2014242) and the Scientific Research Fund of Liaoning University of Technology, China (Grant Nos. X201319 and X201312).

Abstract

In this paper, we theoretically investigate the high-order harmonic generation and attosecond pulse generation when a two-electron He atom is exposed to the intense laser pulse. It shows that due to the two-electron double recombination mechanism, an extended plateau beyond the classical single-electron harmonic has been obtained on the two-electron harmonic spectrum. Further by using this two-electron harmonic extension scheme combined with the two-color field, two supercontinuum bandwidths with 200 eV have been obtained. As a result, a series of sub-60 as extreme ultraviolet (XUV) pulses have been directly generated.

Keyword: 42.65.Ky; 42.65.Re; 32.80.Fb; high-order harmonic generation; attosecond extreme ultraviolet pulse; two-electron model
1. Introduction

The attosecond pulse, as a powerful tool to trace and control the dynamics processes deep inside atoms and molecules, has opened a new time scale of the ultrafast electronic dynamics with unprecedented resolution.[13] Currently, one of the most successful ways towards achieving the attosecond pulse is through high-order harmonic generation (HHG) from the rare gases[46] that can be well understood by means of the classical “ three-step” model consisting of ionization, acceleration and recombination.[7] Further, an attosecond pulse can be produced through superposing the harmonic spectrum in the supercontinuum region.[8] Moreover, to generate the shorter attosecond pulse, many works have been done, such as the few-cycle driving pulse scheme, [9] the polarization gate scheme, [10, 11] the chirp pulse scheme, [12, 13] the two-color or multi-color pulses scheme, [1418] etc. Among them the shortest pulse duration achieved by these means is about 78 as[11] and sub-10 as[12] from the experimental and the theoretical sides, respectively. However, most of the above harmonics and attosecond pulses[918] are generated from the single-electron-approximation (SEA) model.

Recently, Koval et al.[19] and Feng et al.[20] respectively investigated the two-electron harmonic emission process. It is found that due to the two-electron sequential and nonsequential double recombination, there is a new harmonic plateau on the two-electron harmonic spectrum beyond the classical single-electron harmonic. As we know that the broader harmonic plateau we get, the shorter attosecond pulse will be obtained. Therefore, we see that the two-electron harmonic extension scheme is a promising way to generate the ultrashort extreme ultraviolet (XUV) pulse.

Thus, in this paper, by solving the two-electron time-dependent Schrö dinger equation (TDSE), we further investigate the two-electron harmonic emission and attosecond pulse generation when a two-electron He atom is exposed to the intense laser pulse. The results show that due to the two-electron double recombination, many new harmonic cutoffs beyond the classical SEA model have been found. Further, by using the two-color field combined with this harmonic extension scheme, two bandwidths with 200  eV have been obtained, which will support a series of sub-60 as XUV pulses.

2. Theoretical methods

We consider a two-electron model where each electron is allowed to move along the laser polarization direction. Then, the TDSE with the dipole approximation in the length gauge can be written as (atomic units are used throughout this paper unless stated otherwise)

where

and are the electron– nucleus and the electron– electron interactions, respectively. x1 and x2 are the electronic coordinates of the two electrons, which are defined by − 250  a.u. < x1, x2 < 250  a.u. with Δ x1 = Δ x2 = 0.5  a.u. The absorbing regions on x1 and x2 extend over the last 100 grid points. The time space is chosen to be Δ t = 0.1  a.u. The electronic wave function φ (x1, x2, t) can be carried out using the standard second-order split-operator method.[2124] The laser field can be expressed as

for the single-color field and

for the two-color field. Here, Ei, ω i, and τ i, (i = 1, 2) are the amplitudes, the frequencies, and the pulse durations (FWHM) of the 5  fs/800  nm and the 10  fs/1600  nm pulses. φ is the relative phase between the two pulses.

According to the Ehrenfest theorem, [25] the time-dependent dipole acceleration can be calculated by

The HHG spectra can be obtained by Fourier transforming the time-dependent dipole acceleration a(t)

Finally, the attosecond pulse can be obtained by harmonic superposing as follows:

where .

3. Results and discussion

Figure  1 shows the harmonic spectra of the two-electron He model (solid black line) and the SEA He model (dash red line), respectively. The laser field is the single-color 5  fs/800  nm pulse with I1 = 5.0 × 1014  W/cm2 (I = E2). Clearly, for the SEA model, the harmonic cutoff (Ecutoff = 76ω 1) is in good agreement with the classical value Ip1 + 3.17Up, where Ip1 = 24.6  eV is the first ionization potential of the He atom and Up = I/4ω 2 is the ponderomotive energy of the free electron in the laser field. For the two-electron He model, it shows that except for the classical SEA harmonic cutoff, there is a weak extended plateau on the harmonic spectrum.

Fig.  1. HHG spectra from the two-electron He atom (solid black line) and the SEA He model (dash red line). The laser field is 5  fs/800  nm, I1 = 5.0 × 1014  W/cm2.

According to the former investigations, [19, 20] we know that the harmonic extension is caused by the two-electron double recombination. Here, we illuminate it in brief. As we know, the harmonic follows ionization, acceleration and recombination in every half of the laser cycle, and the emitted energies are estimated by EHC = v2/2, v = − ∫ E(t)dt. For instance, under the present single 5-fs laser pulse, there are three main half cycles, as shown in Fig.  2(a), marked as A– B– C, B– C– D, and C– D– E processes, respectively. Thus, recombination of the single electron will occur three times to produce three different cutoffs EHC1 = 44ω 1, EHC2 = 60ω 1, and EHC3 = 38ω 1, as shown in Fig.  2(b). Now, if a two-electron model is considered, each of the two electrons will occur three times recombination under the present laser pulse. If the two electrons are simultaneously recombined (Fig.  2(c)), then the possible return energies will be Ip1(16ω 1) + Ip2(35ω 1) + 2EHC1 = 139ω 1, Ip1 + Ip2 + 2 × EHC2 = 171ω 1, and Ip1 + Ip2 + 2EHC3 = 127ω 1, which is called the sequential double recombination.[20] Secondly, if the double recombination of the two-electron has a time delay, i.e., the double recombination occurs in a half cycle delay (Fig.  2(d)) or one cycle delay (Fig.  2(e)), then the possible return energies will be Ip1 + Ip2 + EHC2 + EHC1 = 155ω 1, Ip1 + Ip2 + EHC3 + EHC2 = 149ω 1 (for half cycle delay), and Ip1 + Ip2 + EHC3 + EHC1 = 133ω 1 (for one cycle delay), which is called the nonsequential double recombination.[19, 20]

Fig.  2. (a)  Laser profile of the 5  fs/800  nm, I = 5.0 × 1014  W/cm2 pulse. (b)  The SEA recombination energies under the above laser pulse. Time is in units of optical cycle (o.c.) of an 800-nm pulse in all the following figures unless stated otherwise. (c)– (e) Two-electron double recombination schemes with recombination delay time equaling to zero, half laser cycle, and one laser cycle, respectively. (c)  Ecutoff = Ip1 + Ip2 + 2EHC1, Ip1 + Ip2 + 2EHC2, Ip1 + Ip2 + 2EHC3. (d)  Ecutoff = Ip1 + Ip2 + EHC2 + EHC1, Ip1 + Ip2 + EHC3 + EHC2. (e)  Ecutoff = Ip1 + Ip2 + EHC3 + EHC1.

From the above investigations, we can confirm that the two-electron harmonic spectra can remarkably enhance the harmonic cutoff. Thus, for further extension and optimization of the harmonic spectra, the two-color field scheme has been employed to produce the two-electron harmonic spectra.

Figure  3(a) shows the two-electron (solid black line) and SEA (dash red line) harmonic spectra of the He atom driven by the two-color laser field, consisting of the 5  fs/800  nm pulse with I1 = 5.0 × 1014  W/cm2 and the 10  fs/1600  nm pulse with I2 = 1.0 × 1014  W/cm2, φ = 0.0π . Clearly, the harmonic cutoff of the two-color two-electron model is remarkably enhanced in comparison with the single-color field case (see Fig.  1) and the two-color SEA mode. As a result, two supercontinuum bandwidths with 200  eV and 208  eV on the first and the second plateau have been obtained, which will support the generation of the favorable attosecond pulses. Figures  3(b) and 3(c) show the controlling wavelength and pulse intensity effects on the two-electron harmonic emission. The other parameters are the same as those in Fig.  3(a). As we know that the maximum harmonic cutoffs (Emax = Ip + nUp) are decided by either the amplitude intensities (relating to the pulse intensities) or the widths (relating to the frequencies) of the half-cycle laser profiles.[26] Thus, with the increasing of the controlling wavelength and the enhancing of the controlling pulse intensity, the harmonic cutoffs have been further extended due to the instantaneous reduction of the central frequency as well as the enhancement of the amplitude of the laser field. However, the large modulations on the harmonic plateaus are unbeneficial to the attosecond pulse generation. Thus, through our calculations, λ 2 = 1600  nm combined with I2 = 1.0 × 1014  W/cm2 are the suitable controlling wavelength and intensity for harmonic extension. Figure  3(d) shows the relative phase effect on the harmonic spectra. Clearly, the maximum cutoff energy is achieved when the relative phase φ = 0.0π . Increasing of the relative phase φ , only leads to harmonic cutoffs move towards the low harmonic orders.

Fig.  3. (a) Two-electron (solid black line) and SEA (dash red line) harmonic emission spectra of He atom driven by the two-color field, consisting of the 5  fs/800  nm, I1 = 5.0 × 1014  W/cm2 pulse and the 10  fs/1600  nm, I2 = 1.0 × 1014  W/cm2, φ = 0.0π pulse. (b)– (d) Controlling wavelength, controlling pulse intensity, and relative phase effects on the two-electron harmonic spectra, respectively.

Figure  4 shows the ionization probabilities of the total ionization, single ionization and the double ionization for the present laser pulse. We see that the single ionization is almost dominated by the ionization process at the beginning of the laser. However, after about 1.0 o.c. (o.c. means the optical cycle of the 800-nm pulse), the double ionization begins to enhance, while the single ionization is decreased but still larger than the double ionization at the end of the laser, which is responsible for different yields between the SEA plateau and the extended new plateau on the harmonic spectra.

Fig.  4. The single, the double, and the total ionization probabilities for the present laser pulse.

To better understand the two-electron harmonic emission characteristic in Fig.  5, we present the optimal two-color two-electron time-frequency harmonic distribution (Fig.  3(a) solid black line), obtained by using the wavelet transformation of the dipole acceleration a(t)[27, 28]

where W(ω 1(t′ − t)) is the Morlet wavelet, which can be expressed as

Here, α = 30 in our calculations. Clearly, the harmonic energies exhibit two parts. One is the SEA region under the 197ω 1 harmonic order, which agrees well with the first plateau on the two-electron spectrum shown in Fig.  3(a). Moreover, each energy receives two similar contributions, named the short and the long quantum paths, [29] which are responsible for the interference on the harmonic spectrum. While, the other one is the extended weak energies higher than 197ω 1, which is caused by the second electron recombination. As we know that the higher harmonic cutoff we get, the higher photon energy of the generated pulse will be obtained. Thus, we see that the extended harmonic plateau caused by the two-electron harmonic emission scheme is much better for the generation of the attosecond pulse with higher photon energy.

Fig.  5. Time-frequency distribution of the optimal two-color two-electron harmonic spectrum.

Figure  6 shows the temporal profiles of the generated attosecond pulses from the optimal two-color two-electron harmonic emission spectrum (Fig.  3(a) solid black line). Particularly, on the first plateau, by properly superposing the harmonics from the 70th to the 110th order or from the 110th to the 150th order, two isolated XUV attosecond pulses with durations of 55  as and 56  as have been obtained as shown in Fig.  6(a). On the second plateau, by properly superposing the harmonics from the 215th to the 255th order, a 52  as XUV pulse train can be obtained as shown in Fig.  6(b). However, due to the fact that the secondary highest attosecond pulse is too weak, this generated attosecond train can be viewed as an isolated attosecond pulse. But if we further superpose the higher harmonics, the secondary highest attosecond pulse will be enhanced, which is unbeneficial to generate the isolated attosecond pulse, as shown in Fig.  6(c) for the example case of superposing the harmonics from the 255th to the 295th order.

Fig.  6. Temporal profiles of the XUV pulses by superposing the harmonics (a) from the 70th to the 110th order or from the 110th to the 150th order, (b) from the 215th to the 255th order, and (c) from the 255th to the 295th order.

4. Conclusion

In conclusion, we theoretically investigate the two-electron harmonic emission and the attosecond pulse generation when the He atom is exposed to the intense laser pulse. The results show an extended harmonic plateau caused by the two-electron double recombination scheme. Further, by using and optimizing the two-color field combined with this two-electron harmonic extension scheme, a series of isolated sub-60as XUV pulses can be obtained.

Reference
1 Brabec T and Krausz F 2000 Rev. Mod. Phys. 72 545 DOI:10.1103/RevModPhys.72.545 [Cited within:1] [JCR: 44.982]
2 Osborne I and Yeston J 2007 Science 317 765 DOI:10.1126/science.317.5839.765 [Cited within:1]
3 Krausz F and Ivanov M 2009 Rev. Mod. Phys. 81 163 DOI:10.1103/RevModPhys.81.163 [Cited within:1] [JCR: 44.982]
4 Zhang G T, Wu J, Xia C L and Liu X S 2009 Phys. Rev. A 80 055404 DOI:10.1103/PhysRevA.80.055404 [Cited within:1] [JCR: 3.042]
5 Zhou Z Y and Yan J M 2008 Chin. Phys. B 17 4523 DOI:10.1088/1674-1056/17/12/032 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429]
6 Feng L Q and Chu T S 2013 Commun. Comput. Chem. 1 52 DOI:10.4208/cicc.2013.v1.n1.6 [Cited within:1]
7 Corkum P B 1993 Phys. Rev. Lett. 71 1994 DOI:10.1103/PhysRevLett.71.1994 [Cited within:1] [JCR: 7.943]
8 Miao X Y and Du H N 2013 Phys. Rev. A 87 053403 DOI:10.1103/PhysRevA.87.053403 [Cited within:1] [JCR: 3.042]
9 Goulielmakis E, Schultze M, Hofstetter M, Yakovlev V S, Gagnon J, Uiberacker M, Aquila A L, Gullikson E M, Attwood D T, Kienberger R, Krausz F and Kleineberg U 2008 Science 320 1614 DOI:10.1126/science.1157846 [Cited within:2]
10 Xia C L, Zhang G T, Wu J and Liu X S 2010 Phys. Rev. A 81 043420 DOI:10.1103/PhysRevA.81.043420 [Cited within:1] [JCR: 3.042]
11 Zhao K, Zhang Q, Chini M, Wu Y, Wang X W and Chang Z H 2012 Opt. Lett. 37 3891 DOI:10.1364/OL.37.003891 [Cited within:2] [JCR: 3.385]
12 Feng L Q and Chu T S 2011 Phys. Rev. A 84 053853 DOI:10.1103/PhysRevA.84.053853 [Cited within:2] [JCR: 3.042]
13 Xu J J, Zeng B and Yu Y L 2010 Phys. Rev. A 82 053822 DOI:10.1103/PhysRevA.82.053822 [Cited within:1] [JCR: 3.042]
14 Zeng Z, Cheng Y, Song X, Li R and Xu Z 2007 Phys. Rev. Lett. 98 203901 DOI:10.1103/PhysRevLett.98.203901 [Cited within:1] [JCR: 7.943]
15 Feng L Q and Chu T S 2012 Chin. Phys. B 21 124204 DOI:10.1088/1674-1056/21/12/124204 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429]
16 Jin C, Wang G L, Wei H, Le A T and Lin C D 2014 Nat. Commun. 5 4003 DOI:10.1038/ncomms5003 [Cited within:1] [JCR: 10.015]
17 Chipperfield L E, Robinson J S, Tisch J W G and Marangos J P 2009 Phys. Rev. Lett. 102 063003 DOI:10.1103/PhysRevLett.102.063003 [Cited within:1] [JCR: 7.943]
18 Band ulet H C, Comtois D, Bisson E, Fleischer A, Pépin H, Kieffer J C, Corkum P B and Villeneuve D M 2010 Phys. Rev. A 81 013803 DOI:10.1103/PhysRevA.81.013803 [Cited within:2] [JCR: 3.042]
19 Koval P, Wilken F, Bauer D and Keitel C H 2007 Phys. Rev. Lett. 98 043904 DOI:10.1103/PhysRevLett.98.043904 [Cited within:3] [JCR: 7.943]
20 Feng L Q, Li W L, Yuan M H, Duan Y B and Chu T S 2014 Phys. Lett. A DOI:10.1016/j.physleta.2014.05.035 [Cited within:4] [JCR: 1.11]
21 Li W L and Han K L 2013 J. Math. Chem. 51 1293 DOI:10.1007/s10910-013-0145-8 [Cited within:1] [JCR: 1.226]
22 Hu J, Han K L and He G Z 2005 Phys. Rev. Lett. 95 123001 DOI:10.1103/PhysRevLett.95.123001 [Cited within:1] [JCR: 7.943]
23 Chu T S, Zhang Y and Han K L 2006 Int. Rev. Phys. Chem. 25 201 DOI:10.1080/01442350600677929 [Cited within:1] [JCR: 4.92]
24 Lu R F, Zhang P Y and Han K L 2008 Phys. Rev. E 77 066701 DOI:10.1103/PhysRevE.77.066701 [Cited within:1] [JCR: 2.313]
25 Burnett K, Reed V C, Cooper J and Knight P L 1992 Phys. Rev. A 45 3347 DOI:10.1103/PhysRevA.45.3347 [Cited within:1] [JCR: 3.042]
26 Feng L Q, Yuan M H and Chu T S 2013 Phys. Plasmas 20 122307 DOI:10.1063/1.4848757 [Cited within:1] [JCR: 2.376]
27 Antoine P, Piraux B and Maquet A 1995 Phys. Rev. A 51R1750 DOI:10.1103/PhysRevA.51.R1750 [Cited within:1] [JCR: 3.042]
28 Feng L Q and Liu H 2014 Commun. Comput. Chem. 2 47 DOI:10.4208/cicc.2014.v2.n2.1 [Cited within:1]
29 Mairesse Y, Bohan A D, Frasinski L J, Merdji H, Dinu L C, Monchicourt P, Breger P, Kovačev M, Taïeb R, Carré B, Muller H G, Agostini P and Salières P 2003 Science 302 1540 DOI:10.1126/science.1090277 [Cited within:1]