⇦Chinese Physics Letters, 2017, Vol. 34, No. 1, Article code 014206 Exploration of High-Harmonic Generation from the CS$_2$ Molecule by the Lewenstein Method in Two-Color Circularly Polarized Laser Field * Hong-Dan Zhang(张宏丹)1, Jing Guo(郭静)1,2**, Yan Shi(史岩)1,3, Hui Du(杜慧)1, Hai-Feng Liu(刘海凤)1, Xu-Ri Huang(黄旭日)2, Xue-Shen Liu(刘学深)1**, Jun Jing(景俊)1 Affiliations 1Institute of Atomic and Molecular Physics, Jilin University, Changchun 130012 2Institute of Theoretical Chemistry, Jilin University, Changchun 130012 3National Center of Quality Supervision and Inspection of Automobile Spare Parts, Changchun 130012 Received 6 September 2016 *Supported by the National Natural Science Foundation of China under Grant Nos 11574117, 11271158, 61575077 and 11575071.
**Corresponding author. Email: gjing@jlu.edu.cn; liuxs@jlu.edu.cn Citation Text: Zhang H D, Guo J, Shi Y, Du H and Liu H F et al 2017 Chin. Phys. Lett. 34 014206 Abstract The high harmonic generation (HHG) from the CS$_2$ molecule in intense laser fields is investigated using the extended Lewenstein method. The initial state is the highest-occupied molecular orbital of the CS$_2$ molecule, which can be well described by Gaussian wave packet using GAMESS-UK package. Compared with the case of the elliptical laser, the HHG can be extended in two-color circularly polarized laser field. The time-frequency analysis and classical electron trajectory as well as the ionization yield curve are also presented to further explain the underlying mechanism. After adding a static electric field on the $z$-direction, the single quantum path control is realized and the supercontinuum spectra are obtained. Moreover, an isolated 110 as pulse can be obtained by superposing the harmonics from 130th to 180th order. DOI:10.1088/0256-307X/34/1/014206 PACS:42.65.Ky, 42.65.Re, 32.80.Rm © 2017 Chinese Physics Society Article Text In strong-field physics, high harmonic generation (HHG) from molecules such as N$_2$, O$_2$ and CO$_2$ has attracted much attention in recent years.[1-3] It can be used as tools for measuring and controlling the dynamic process of electrons in molecules, especially in ultrafast electron dynamics. Explicitly, such a phenomenon can be described as follows:[4] a molecule emits its electron to the continuous state through tunneling or multi-photon ionization, then the electron is propagated in the continuum and accelerated in the laser field; finally, it comes back to the core and recombines with the parent ion. Recently, the HHG of atoms has been almost well understood theoretically and experimentally by scientists. The interaction between laser and molecules shows more extensive phenomena, such as the multicenter structure,[5-7] molecular orbital tomographic imaging,[8] and the quantum interference between high-order harmonics originated from different cores. Furthermore, much effort has been given to broadening the bandwidth of the supercontinuum of HHG and shortening the pulse duration for the potential applications with much higher time resolutions.[9-11] It has been shown that HHG depends on the ellipticity of the pulse sensitively, thus the recombination electrons can be gated into one half cycle by using the few-cycle polarization gating (PG) technique, and then a broadband supercontinuum in the plateau region can be produced.[11-14] With the development of the laser technology, the mid-infrared wavelength laser can be generated.[13-15] Under this condition, the tunneling process is predominant, thus the dynamics of molecules and atoms shows some new phenomena, which attract a great deal of interest and have been widely investigated recently.[16-18] In 2005, Itatani et al.[8] realized control of the high-order harmonic generation process by changing the alignment of molecules. Recently, the influence of nuclear motion in a molecule has been widely studied.[19-23] In 2014, Guo et al.[24] investigated the influence of different vibrational states on HHG and isolated attosecond pulse from a molecule. Recently, the CS$_2$ molecule has attracted some attention in theory and experiment.[25,26] The strong-field Lewenstein model for molecules has been widely used.[28-31] In this work, we choose the system of the CS$_2$ molecule and consider its highest occupied molecular orbital (HOMO) as the initial state, which is constructed by the GAMESS-UK package. The transition amplitude is $$\begin{align} r(t)=\,&i\int^{\infty}_{0}d\tau[\pi/(\varepsilon+i\tau/2)]^{3/2}d^{*}_{\rm rec}[{\boldsymbol p}_{\rm st}(t,\tau)\\ &-{\boldsymbol A}(t)]\times d_{\rm ion}[{\boldsymbol p}_{\rm st}(t,\tau)-{\boldsymbol A}(t-\tau)]\\ & \times \exp[-i{\boldsymbol S}(p_{\rm st},t,\tau)]{\boldsymbol E}(t-\tau)+{\rm c.c.},~~ \tag {1} \end{align} $$ where $d_{\rm ion}[p_{\rm st}(t,\tau)-{\boldsymbol A}(t-\tau)]$ is the transition amplitude of ionization, $d^{*}_{\rm rec}[{\boldsymbol p}_{\rm st}(t,\tau)-{\boldsymbol A}(t)]$ is the transition amplitude of recombination, $\varepsilon$, $\tau$, ${\boldsymbol E}(t)$ and ${\boldsymbol A}(t)$ give the positive regularization constant, the traveling time of free electron in the laser, the electric field and the vector potential of the electric field, respectively, ${\boldsymbol p}_{\rm st}(t,\tau)$ is the stationary momentum, and ${\boldsymbol S}({\boldsymbol p}_{\rm st},t,\tau)$ is the quasiclassical action of the electron, which can be calculated as follows: $$\begin{alignat}{1} {\boldsymbol S}({\boldsymbol p}_{\rm st},t,\tau)=\,&\tau I_{\rm p}-\frac{1}{2}[{\boldsymbol p}^{2}_{{\rm st}y}(t,\tau)+{\boldsymbol p}^{2}_{{\rm st}z}(t,\tau)]\\ &+\frac{1}{2}\int^{t}_{t-\tau}[{\boldsymbol A}^{2}_{y}(t'')+{\boldsymbol A}^{2}_{z}(t'')]dt,~~ \tag {2} \end{alignat} $$ where $I_{\rm p}$ is the ionization potential of the CS$_{2}$ molecule, ${\boldsymbol p}_{{\rm st}y}(t,\tau)$, ${\boldsymbol A}_{y}(t'')$ and ${\boldsymbol p}_{{\rm st}z}(t,\tau)$, ${\boldsymbol A}_{z}(t'')$ are the stationary momentum and vector potential in the $y$ and $z$ directions, respectively. In our calculation, the lasers are set on the $y$ and $z$ directions. Thus the transition amplitude would be written as $$\begin{align} r_{y}(t)=\,&i\int^{\infty}_{0}d\tau[\pi(\varepsilon+i\tau/2)]^{3/2}d^{*}_{{\rm rec}y}[{\boldsymbol p}_{\rm st}(t,\tau)-{\boldsymbol A}(t)]\\ &\times \{d_{{\rm ion}y}[{\boldsymbol p}_{\rm st}(t,\tau)-{\boldsymbol A}(t,\tau)]{\boldsymbol E}_{y}(t-\tau)\\ &+d_{{\rm ion}z}[{\boldsymbol p}_{\rm st}(t,\tau)-{\boldsymbol A}(t,\tau)]{\boldsymbol E}_{z}(t-\tau)\}\\ &\times \exp[-i{\boldsymbol S}({\boldsymbol p}_{\rm st},t,\tau)]+{\rm c.c.},~~ \tag {3} \end{align} $$ $$\begin{align} r_{z}(t)=\,&i\int^{\infty}_{0}d\tau[\pi(\varepsilon+i\tau/2)]^{3/2}d^{*}_{recz}[{\boldsymbol p}_{\rm st}(t,\tau)-{\boldsymbol A}(t)]\\ &\times \{d_{{\rm ion}y}[{\boldsymbol p}_{\rm st}(t,\tau)-{\boldsymbol A}(t,\tau)]{\boldsymbol E}_{y}(t-\tau)\\ &+d_{{\rm ion}z}[{\boldsymbol p}_{\rm st}(t,\tau)-{\boldsymbol A}(t,\tau)]{\boldsymbol E}_{z}(t-\tau)\}\\ &\times \exp[-i{\boldsymbol S}({\boldsymbol p}_{\rm st},t,\tau)]+{\rm c.c.}~~ \tag {4} \end{align} $$ The corresponding HHG power spectrum is proportional to the time-dependent dipole moment (in acceleration form) $$ P_{\rm a}(\omega)\sim|\int[r_{y}(t){{\boldsymbol e}}_{y}+r_{z}(t){{\boldsymbol e}}_{z}]\exp(-i\omega t)dt|^{2}.~~ \tag {5} $$ The attosecond pulse can be calculated by superposing several harmonics as follows: $$ I(t)=|\sum\limits_{q}a_{q}\exp(iq\omega t)|^{2},~~ \tag {6} $$ where $a_{q}=\int r(t)\exp(-iq\omega t)dt$.
cpl-34-1-014206-fig1.png
Fig. 1. The initial density function of the HOMO orbital of CS$_2$ with symmetric stretch mode calculated by GAMESS-UK.
We first present the initial density function of the CS$_2$ molecule with symmetric stretch mode, which is calculated by GAMESS-UK in Fig. 1. As can be seen, there are four parts distributed like a butterfly, and the initial density function is symmetric with respect to the axes $y$ and $z$. In our simulation the expression of the elliptical field is $$ {\boldsymbol E}(t)=E_{0}f(t)[\varepsilon\sin(\omega_1 t){{\boldsymbol e}_{y}}+\cos(\omega_1 t){{\boldsymbol e}_{z}}],~~ \tag {7} $$ and the two-color circularly polarized field can be represented as $$\begin{align} E_{y}(t)=\,&E_{0}f(t)[(\sin(\omega_1 t)-\sin(\omega_2 t))],~~ \tag {8} \end{align} $$ $$\begin{align} E_{z}(t)=\,&E_{0}f(t)[(\cos(\omega_1 t)+\cos(\omega_2 t))],~~ \tag {9} \end{align} $$ where $f(t)=\exp[-2\ln2(t-T/2)^2/\tau^2]$ is the pulse envelope, $E_{0}=1\times10^{14}$ W/cm$^{2}$ is the amplitude of laser field, $\omega_1=0.038$ a.u. and $\omega_2=0.076$ a.u. are the frequencies corresponding to the laser field of 1200 nm and 600 nm, respectively, and $\varepsilon=0.1$ is the ellipticity of the elliptical field. The durations of $t$ for the time evolution are $T=12 T_0$ and $\tau_{\rm d}=5$ fs, where $T_0$ is the optical cycle (o.c.) of a 1200-nm pulse. Figure 2(a) is the electric field of elliptical laser field (red solid line) which shows the clockwise spiral along the axis of time, and the black dashed line is the projection to the $y$–$z$ plane. Figure 2(b) shows the two-color circularly polarized laser field (red solid line). The projection to the $y$–$z$ plane represents a three-fold spatiotemporal symmetry like a clover leaf. The laser has three peaks per cycle, which means that the electron can recombine three times per laser cycle. Figure 3 shows the HHG spectra from the CS$_2$ molecule driven by the elliptical laser and counter-rotating laser as well as the H atom driven by the counter-rotating laser, respectively. The black solid line is the spectra driven by the elliptical laser with the cutoff of about 50th order. The red dotted line is the one driven by a two-color circularly polarized laser with extended cutoff of about 80th order. Moreover, the intensity has been increased by three orders of magnitude. The blue dashed line is the spectra driven by the two-color laser of the H atom, which shows the similar cutoff order with the decreasing emission efficiency of HHG above the 40th order.
cpl-34-1-014206-fig2.png
Fig. 2. The electric field of elliptical laser (red solid line) (a) and two-color circularly polarized laser (red solid line) (b). The black dashed lines in (a) and (b) are the projections to the $y$–$z$ plane of these two fields.
cpl-34-1-014206-fig3.png
Fig. 3. Harmonic spectra of the CS$_2$ molecule induced by elliptical laser (black solid line) and two-color circularly polarized laser (red dotted line), and the harmonic spectra of the H atom induced by two-color circularly polarized laser (blue dashed line).
To interpret the mechanism of the extension in the HHG plateau, we present the time-frequency analysis along the $y$ and $z$ directions for both the lasers (see Figs. 4(a)–4(d)). Figures 4(a) and 4(b) show the distribution in the $y$ and $z$ directions of the elliptical laser, respectively. As can be seen in Fig. 4(a), there are three peaks in the $y$-direction with a relatively weak peak at about $-$0.5 o.c. Thus the dominant contribution to the emission of HHG on the $y$-direction comes from the other two peaks. However in Fig. 4(b), there are two peaks contributed to HHG, which shows the same tendency as that in Fig. 4(a). Moreover, the intensity of HHG in the $y$-direction is higher than that in the $z$-direction. In Figs. 4(c) and 4(d), there are five similar peaks and the intensity in the $y$-direction is higher than that in the $z$-direction. Since the molecular axis is along the $z$-direction, we will take the results on the $z$-direction as an example. As shown in Fig. 4(b), there are two peaks contributed to the emission of HHG, and the maximum order of emitted HHG located at about the 55th order, which is consistent with that in Fig. 3. However in Fig. 4(d), there are about five peaks, and the peak intensity about $-$0.5 o.c. is relatively weak, which can be ignored. Thus there are four peaks with comparable intensity contributed to the photon emission. That is the reason why the HHG has so many modulations. Furthermore, the maximum order of HHG has been extended to about the 80th order, which is in accordance with that in Fig. 3.
cpl-34-1-014206-fig4.png
Fig. 4. The time-frequency analysis of the harmonic spectrum along (a) the $y$-direction and (b) the $z$-direction in the elliptical laser field. The time-frequency analysis of harmonic spectrum along (c) the $y$-direction and (d) the $z$-direction in the two-color circularly polarized laser field.
To further understand this phenomenon, we draw the trajectory of the electron moving in the laser fields. The trajectories in Figs. 5(a) and 5(b) are under the elliptical laser and the two-color circularly polarized laser and the combined field with static electric field, respectively. Figure 5(a) shows the paths with different initial transverse velocities $v_{0}$ and emission times $t_{\rm r}$, and they appear to be the elliptical shape. The trajectory with emission time $t_{\rm r}= 0.41$ o.c. moves further than the trajectory with emission time $t_{\rm r}= -0.05$ o.c., which is in accordance with the fact that the harmonic order located at time $t_{\rm r}= 0.41$ o.c. is higher than that located at time $t_{\rm r}= -0.05$ o.c. in Fig. 4(b). However in Fig. 5(b), the trajectories are also with different proper initial transverse velocities and emission times, which appear to be the triangular way. It can be visibly seen that electrons move further in the latter case, which can well explain the extension of the plateau of harmonic spectra.
cpl-34-1-014206-fig5.png
Fig. 5. The classical electron trajectories in (a) the elliptical laser and (b) the two-color circularly polarized laser and the two-color circularly polarized laser combined with static electric field, respectively.
cpl-34-1-014206-fig6.png
Fig. 6. Ionization yield of the CS$_2$ molecule in the elliptical laser (black solid line) and the two-color circularly polarized laser (blue dotted line).
cpl-34-1-014206-fig7.png
Fig. 7. Harmonic spectra generated by the two-color circularly polarized laser combined with the static electric field of $\alpha=0.1$, $\alpha=0.3$ and $\alpha=0.35$, respectively.
Figure 6 shows the ionization yield of the CS$_2$ molecule under these two lasers. The ionization yield of the CS$_2$ molecule in the case of the two-color circularly polarized field (blue dotted line) is higher than that of the elliptical laser (black solid line), which is consistently with the higher intensity of HHG of the two-color circularly polarized laser case (red dotted line) in Fig. 3.
cpl-34-1-014206-fig8.png
Fig. 8. The time-frequency analysis of harmonic spectra generated by the two-color circularly polarized laser combined with the static electric field of $\alpha=0.35$.
cpl-34-1-014206-fig9.png
Fig. 9. The temporal profiles of the attosecond pulses generated by superposing some orders of harmonic spectra in the two-color circularly polarized laser (red solid line) and the combined laser with $\alpha=0.35$ (blue dashed line).
To optimize the harmonic spectra and to extend the plateau of HHG, we add the static electric field to the $z$-direction with the form $E_{\rm static}=\alpha E_{0}$, where $\alpha$ is the static parameter. In our calculation, we choose $\alpha=0.1$, $\alpha=0.3$ and $\alpha=0.35$, respectively. Figure 7 shows the harmonic spectra of these three cases. As can be seen, the plateau of HHG can be extended with the increasing parameter $\alpha$. When $\alpha=0.35$, the harmonic spectra becomes supercontinuum with a bandwidth from 120th to 180th order. We represent the time-frequency analysis of the $\alpha=0.35$ case in Fig. 8. There are five peaks contributing to the emission of HHG under the 80th order. That is the reason why the harmonic spectra less than the 80th order has so many modulations. However, there is only one peak contributing to the emission of harmonic spectra above the 80th order. It can be clearly seen that the long trajectory is suppressed, which realizes the single quantum path control. After adding the static electric field, the electron will have a transverse velocity (along the $z$-direction) and will obtain energy from the static electric field, which leads to the enhancement of the short trajectory and the suppression of the long trajectory. The red dashed-dotted line in Fig. 5(b) shows the trajectory of an electron with emission time $t_{\rm r}=0.29$ o.c. within the long path domain. The trajectory can be nearly ignored relative to others which is in accordance with the above long path suppression statement. Figure 9 shows the attosecond pulse generated from the two-color circularly polarized laser and the two-color circularly polarized laser with static electric field of $\alpha=0.35$, respectively. By superposing some harmonic spectra from the 50th to the 70th order, we gain the attosecond pulse with duration of 190 as, as shown in Fig. 9. However, with the existence of static electric field, the isolated attosecond pulse is realized and the duration has reduced to 110 as by superposing the orders from 130th to 150th. In summary, we have investigated the HHG and attosecond pulse from the CS$_{2}$ molecule with the elliptical laser and the two-color circularly polarized laser by the Lewenstein model. The results show that the plateau of HHG generated by the two-color polarized laser is extended and the intensity is about three orders higher than that of the elliptical laser. Both the time-frequency analysis and the classical trajectory as well as the ionization yield can explain the underlying mechanism well. By adding the static electric field of $\alpha=0.35$ to the $z$ direction, a supercontinuum spectra can be gained and the single quantum path control is realized. Finally, by superposing proper orders of harmonic spectra from the CS$_{2}$ molecule with the combined laser, an isolated attosecond pulse with duration of 110 as can be obtained. References Controlling the Interference of Multiple Molecular Orbitals in High-Harmonic GenerationHigh-order harmonic generation of N 2 molecule in two-color circularly polarized laser fieldsAccurate prediction of interference minima in linear molecular harmonic spectra by a modified two-center modelPlasma perspective on strong field multiphoton ionizationAlignment dependence of high-order harmonic generation from N 2 and O 2 molecules in intense laser fieldsTwo-center interference in high-order harmonic generation from heteronuclear diatomic moleculesNonsequential double ionization of triatomic molecules in strong laser fieldsTomographic imaging of molecular orbitalsAttosecond physicsIsolated Single-Cycle Attosecond PulsesMid-infrared modulated polarization gating for ultra-broadband supercontinuum generationChirp of the single attosecond pulse generated by a polarization gatingSingle attosecond pulse and xuv supercontinuum in the high-order harmonic plateauBroadband supercontinuum generation method combining mid-infrared chirped-pulse modulation and generalized polarization gatingMacroscopic effects for quantum control of broadband isolated attosecond pulse generation with a two-color fieldDouble Optical Gating of High-Order Harmonic Generation with Carrier-Envelope Phase Stabilized LasersInterference effects in high-order harmonic generation by homonuclear diatomic moleculesAlignment-dependent nonsequential double ionization of N 2 in intense laser fields: The role of different valence orbitalsResonance effects in two-photon double ionization of H 2 by femtosecond XUV laser pulsesSignatures of nuclear motion in molecular high-order harmonics and in the generation of attosecond pulse trains by ultrashort intense laser pulsesCorrelated Electron and Nuclear Dynamics in Strong Field Photoionization of H 2 + Probing H 2 + vibrational motions with high-order harmonic generationQuantum-path control and isolated-attosecond-pulse generation using H 2 + molecules with moving nuclei in few-cycle laser pulsesInfluence of vibrational states on high-order-harmonic generation and an isolated attosecond pulse from a N 2 moleculeExperimental and theoretical study on nonsequential double ionization of carbon disulfide in strong near-IR laser fieldsFilamentation of light in carbon disulfideThe GAMESS-UK electronic structure package: algorithms, developments and applicationsHIGH-ORDER HARMONIC GENERATION IN DIATOMIC MOLECULES: QUANTUM INTERFERENCE, NODAL STRUCTURES AND MULTIPLE ORBITALSLaser-induced nonsequential double ionization in diatomic molecules: one- and two-centre rescattering scenariosTheory of high-harmonic generation by low-frequency laser fieldsElliptical High-Order Harmonic Generation from H + 2 in Linearly Polarized Laser Fields
[1] Worner H J et al 2010 Phys. Rev. Lett. 104 233904
[2] Du H et al 2016 Chin. Phys. B 25 043202
[3] Cui X et al 2016 Chin. Phys. B 25 033205
[4] Corkum P B 1993 Phys. Rev. Lett. 71 1994
[5] Zhou X X et al 2005 Phys. Rev. A 72 033412
[6] Zhu X et al 2011 Opt. Express 19 436
[7] Pei L and Guo C 2010 Phys. Rev. A 82 021401(R)
[8] Itatani J et al 2004 Nature 432 867
[9] krausz F and Ivanov M 2009 Rev. Mod. Phys. 81 163
[10] Sansone G et al 2006 Science 314 443
[11] Hong W et al 2010 Opt. Express 18 11308
[12] Chang Z 2005 Phys. Rev. A 71 023813
[13] Chang Z 2004 Phys. Rev. A 70 043802
[14] Du H and Hu B 2010 Opt. Express 18 25958
[15] Lan P et al 2009 Phys. Rev. A 79 043413
[16] Mashiko H et al 2008 Phys. Rev. Lett. 100 103906
[17] Odžak S and Miloševi D B 2009 Phys. Rev. A 79 023414
[18] Jia X Y et al 2009 Phys. Rev. A 80 053405
[19] Guan X et al 2013 Phys. Rev. A 88 043402
[20] Bandrauk A D et al 2009 J. Phys. B 42 075602
[21] Silva R E F et al 2013 Phys. Rev. Lett. 110 113001
[22] Zhao J and Zhao Z 2008 Phys. Rev. A 78 053414
[23] Ge X L et al 2014 Phys. Rev. A 89 023424
[24] Guo J et al 2014 Phys. Rev. A 90 053410
[25] Zuo W L et al 2016 Phys. Rev. A 93 053402
[26] Besse V et al 2014 Phys. Rev. A 89 043840
[27]GAMESS-UK is a package of ab initio programs (see http://www.cfs.dl.ac.uk/gamess-uk/index.shtml)
Guest M F et al 2005 Mol. Phys. 103 719
[28] Augstein B B and Figueira de M F C 2012 Mod. Phys. Lett. B 26 1130002
[29] Figueira de M F C 2009 J. Phys. B 42 105602
[30] Lenwenstein M et al 1994 Phys. Rev. A 49 2117
[31] Zhang B and Zhao Z X 2013 Chin. Phys. Lett. 30 023202