## Abstract

We experimentally investigate the angular distribution and the laser intensity dependence of the two-center interference in high-order harmonic generation (HHG) from O_{2} molecules by comparing with CO_{2} molecules. Through the measurement of both the temporal evolution and the angular distribution of HHG, the characteristic enhancement and suppression are observed, which can be well explained by the modified interference model. Furthermore, we demonstrate that the spectral region of the constructive enhancement in aligned O_{2} molecules can be shifted by tuning the driving laser intensity.

©2010 Optical Society of America

## 1. Introduction

In the semiclassical three-step model of high-order harmonic generation (HHG), active electrons first tunnel through the potential barrier, are then accelerated by laser fields, and finally recombine with parent ions to emit high-energy photons [1–2]. This process can provide the approach to probe electronic and nuclear dynamics of atoms and molecules with an unprecedented fast temporal resolution [3–6]. For example, recently, Wörner *et al*. [7] followed a chemical reaction of less than an attosecond between the stretched and compressed geometry of weakly vibrationally excited Br_{2} using high-harmonic interferometry, and Goulielmakis *et al*. [8] probed the real-time valence electron motion through the spectrally resolved absorption of an attosecond extreme-ultraviolet pulse.

The harmonic emission from impulsively aligned CO_{2} shows a special feature that the harmonic yield in some spectral region exhibits an inverted modulation versus the molecular alignment parameter, and the modulation inversion is attributed to the interference of recombining electron from the two oxygen atoms [9–16]. Boutu *et al.* [17] even gave the phase measurement of harmonics generated in aligned CO_{2} molecules, which provided proof for an interpretation in terms of destructive quantum interference. However, the harmonic orders that are inversely modulated differed in different experiments from CO_{2} molecules, and the effect of the molecular potential on the recolliding electron wave packet was treated differently in the interpretation of experimental data [9–13]. Therefore, we attributed such phenomena to the laser field distortion of the molecular HOMO, which should play a key role in such quantum system [16]. Moreover, the role of the multiple HOMOs effect was also proposed for the laser intensity dependence in aligned molecules recently [18–23]. Probably the field distortion of the molecular HOMO and the multiple HOMOs effect both play roles in such system. We would like to point out that, among the explanations for the laser intensity dependence on the modulation inversion of the harmonic yield, our model focuses on the field distortion effect on the molecular HOMO and has the merit of simplicity. In this letter, we further represent the new experimental results on O_{2} molecules which can be well explained by our model.

The previous results [9–10] [12–23] about the interference (or the laser intensity dependence) mainly focus on CO_{2} and N_{2} molecules, and there is a comparative lack of the investigation on O_{2} molecule. Infrequently for O_{2} molecule, Itatani *et al.* [3] have showed that the high harmonics are enhanced when aligned near 45° with respect to the laser polarization due to the π* _{g}* symmetry of its HOMO, and Vozzi

*et al.*[11] reported the order dependence of the harmonic yield due to the two-center interference. Here we further measure the angular distribution and the laser intensity dependence of the harmonic yield from O

_{2}molecules.

In this work, we investigate the two-center interference in aligned O_{2} molecules and perform a more comprehensive study of the interference characteristics by comparing with CO_{2} molecules. Through the measurement of both the temporal evolution and the angular distribution of harmonic intensity, the detailed characteristics of intensity enhancement and suppression are observed, and the effects of the laser intensity and the two-centre interference are summarized. Moreover, the shift of the constructive spectral region in aligned O_{2} molecules is demonstrated by tuning the driving laser intensity.

## 2. Theoretical model and experimental setup

#### 2.1 Theoretical model

For the high-order harmonic emission from impulsively aligned O_{2} or CO_{2} molecules irradiated by an ultra short laser pulse, the quantum interference is attributed to the interference of recombining electron from the two oxygen atoms [10–11]. For such antisymmetric atomic orbitals of the linear molecules, the interference factor is defined as$\left[1-\mathrm{cos}(2\pi R\mathrm{cos}\theta /{\lambda}_{B})\right]$, and the harmonic ratio *S/S _{0}* between the signal from aligned molecules and that from an isotropic ensemble, at the considered time delay

*t*

_{0}, can be simply approximated as follows:

where *θ* is the angle between the laser polarization and the internuclear axis, *P _{Θ}*(

*θ, t*) [10, 11, 15] is the angle distribution of the molecular ensemble at time

_{0}*t*, probed at an angle

_{0}*Θ*from the aligning field direction;

*n*is the harmonic order,

*R*is the internuclear separation,

*C*is a fitting constant;

*λ*

_{B}(

*n*) is the de Broglie wavelength of an electron responsible for the emission of the

*n*th harmonic, given here by ${\lambda}_{B}(n)=h/\sqrt{2{m}_{e}{E}_{k}}$, where

*m*

_{e}is the electron mass,

*E*is the electron kinetic energy. Considering the field distortion of the molecular potential, we treat the kinetic energy of the returning electron as ${E}_{k}=nh{v}_{0}-\delta {I}_{p}$ (δ≥0, variable) [16], where

_{k}*nhv*

_{0}is the emitted photon energy (

*hv*

_{0}=1.55eV), and $\delta {I}_{p}$is the effective molecular potential which is closely related to the laser field. The effective ionization potentials could also be expanded to the new understanding of the continuum electron wave functions, and the field-free scattering states are now known to give good agreement with experiments [23] [25].

In this model, we directly and simply convolve the molecular distribution with the interference intensity$\left[1-\mathrm{cos}(2\pi R\mathrm{cos}\theta /{\lambda}_{B})\right]$ for its simplification, and this is a rather crude approximation and a strongly simplified model to give a qualitative explanation. The laser field distortion of the molecular HOMO should play a key role through the effective molecular potential. The order dependence is shown as a function of the harmonic order *n*; the angle dependence is shown as a function of the pump-probe polarization angle *Θ*; and the laser intensity dependence is shown as the effective molecular potential$\delta {I}_{p}$ through by${\lambda}_{B}(n)=h/\sqrt{2{m}_{e}(nh{\nu}_{0}-\delta {I}_{p})}$.

#### 2.2 Experimental setup

The experiments are performed by using a Ti:sapphire-based chirped pulse amplification laser system (Spectra-Physics, TSA-25), which produces 50 fs laser pulses at 800 nm center wavelength with a repetition rate of 10 Hz. The laser pulse is split into two beams, one used as the pump pulse (for aligning molecules) and the other as the probe pulse (for driving HHG from the molecules). The two beams are collinearly focused with variable pump-probe time delay on a pulsed supersonic molecular beam located in a high-vacuum interaction chamber. The aligning laser energy is kept constant and the intensity is estimated to be 4 × 10^{13} W/cm^{2} at the interaction region within the gas jet. The probe laser energy is adjustable by using a half-wave plate and a high-extinction film polarizer. The stagnation pressure of CO_{2} gas is around 2 bars, leading to a rotational temperature of tens of kelvins. The generated high-order harmonics are detected by a home-made flat-field grating spectrometer equipped with a soft X-ray charge-coupled device camera (CCD, Princeton Instruments, SX 400).

The experimentally measured 23rd harmonic ratio S/S_{0}, with (S) and without (S_{0}) the aligning laser pulses, as functions of pump-probe delay are represented in Fig. 1
(the black solid lines and left axis, the probe laser intensity is 2.8 × 10^{14} W/cm^{2}, the polarizations of the aligning and probe pulse are parallel). The calculated time-dependent alignment parameters $<{\mathrm{cos}}^{2}\theta >(t)$ [15–16] are also shown in Fig. 1 (the red dashed curves and right axis, the initial rotational temperature is taken to be 80 K). It is evident that the modulation of harmonic signal for CO_{2} molecules is reversely matched with that of the molecular alignment parameter$<{\mathrm{cos}}^{2}\theta >$, while the modulation of harmonic signal for O_{2} molecules is matched with the molecular alignment parameter $<{\mathrm{cos}}^{2}\theta >$ exactly. The phenomena are regarded as the evidence of interference of the recombining electron originated from two oxygen atoms [10–11].

## 3. Results and discussion

#### 3.1 Experimental results

For the detailed observation of the interference characteristics from different harmonic orders, we present the harmonic ratio S/S_{0} at the half and quarter revival from three different harmonic orders, which are shown in Fig. 2
together with the calculated time-dependent alignment parameters$<{\mathrm{cos}}^{2}\theta >$. For O_{2} molecules shown in Figs. 2 (a) - (c), the results are measured at the driving laser intensity of 2.8 × 10^{14} W/cm^{2} for 21st, 29th, and 37th harmonic orders, and they display that the constructive interference gradually weakens and even inverses from 21st to 37th harmonic order. For CO_{2} molecules shown in Figs. 2 (d) - (f), the results are measured at the driving laser intensity of 2.0 × 10^{14} W/cm^{2} for 21st, 25th, and 29th harmonic orders, and they display that the destructive interference gradually weakens and even inverses from 21st to 29th harmonic order. Both the two gases display the order dependence of the two-center interference.

We also measure the angular distributions of HHG as functions of the pump-probe polarization angle (the polarization of the aligning pulse is fixed) at the top alignment of the half revival, and the results are shown in Fig. 3
. Due to the symmetry, we only measure the angle in 0° – 90° and make the polar plot. The angular distributions of the 19th (red star) and 27th (black dot) harmonic orders from aligned O_{2} molecules are shown in Fig. 3(a), and the same distributions from aligned CO_{2} molecules are shown in Fig. 3(b). The driving laser intensities for both the two gases are estimated to be 2.8 × 10^{14} W/cm^{2}. These phenomena of the two gases at the certain laser intensity are similar to the observation of Mairesse *et al*. [24], showing that the harmonic emission of O_{2} molecules is maximized at 0° while the emission of CO_{2} molecules is maximized at 90°. The simulated results fitted with Eq. (1) are also shown in Fig. 3(c) for O_{2} molecules and in Fig. 3(d) for CO_{2} molecules. The simulated results in Fig. 3(c) for O_{2} molecules are using a slightly longer *R* = 0.133nm and $\delta {I}_{p}$ = 0eV, and the simulated results in Fig. 3 (d) for CO_{2} molecules are using a slightly longer *R* = 0.245nm and $\delta {I}_{p}$ = 13.77eV. These figures show that the simulations can well reproduce the experimental results, and the angular distributions are dependent on the harmonic orders, which result from the two-center interference.

The order dependence of the characteristic enhancement and suppression are obvious in both gases through the temporal evolution (shown in Fig. 2) and the angular distribution (shown in Fig. 3). Due to that the internuclear separation of CO_{2} molecule is much longer than that of O_{2} molecule, the destructive interference first appears at the 21st harmonic order of CO_{2} molecules, and then appears at the 37th harmonic order of O_{2} molecules in our observation window.

We then further investigate the role of the laser intensity, which can shift the spectral region of the two-center interference and can inverse the modulation of the harmonic emission. For O_{2} molecules shown in Fig. 4(a)
, we measure the experimental spectral distributions of the harmonic ratio *S/S _{0}* with three different driving laser intensities at

*Θ*= 0°, together with the calculations by numerically solving Eq. (1) using the internuclear separation

*R*= 0.121nm (the solid curves), 0.133nm (a slightly longer

*R,*the dashed curves) and the effective molecular potential $\delta {I}_{p}$ = 0eV (the thick curves), 12.1eV (the thin curves). We can observe that the spectral region of the constructive enhancement (

*S/S*> 1) in aligned O

_{0}_{2}molecules can be shifted to the higher harmonic order by enhancing the driving laser intensity: the constructive enhancement appears for the 19th - 23rd harmonics at the driving laser intensity of 1.8 × 10

^{14}W/cm

^{2}(circles), for the 19th - 29th harmonics at 2.3 × 10

^{14}W/cm

^{2}(squares), and for the 19th – 31st harmonics at 2.8 × 10

^{14}W/cm

^{2}(triangles). However, the destructive suppression (

*S/S*< 1) in aligned O

_{0}_{2}molecules should occur at the higher orders of harmonics, which need the much higher intensity of the driving laser field. At the maximal driving laser intensity (2.8 × 10

^{14}W/cm

^{2}) under our experimental condition, the destructive suppression (

*S/S*< 1) in aligned O

_{0}_{2}molecules appears for the 33rd-37th harmonic orders, and it is different from the observation of Vozzi

*et al.*[11] that the suppression appears for the 41st–47th harmonic orders at the driving laser intensity of 3.5 × 10

^{14}W/cm

^{2}. Furthermore, our experimental results of O

_{2}molecules are best fitted by the electron kinetic energy as ${E}_{k}=nh{v}_{0}$ (the thick blue curve), while the results of Vozzi

*et al.*are more fitted by the electron kinetic energy as ${E}_{k}=nh{v}_{0}-{I}_{p}$ (the red thin curve) at the comparative high laser intensity. These results of O

_{2}molecules show that the higher driving laser intensity corresponds to the larger effective molecular potential. As a comparison, we also plot the spectral distributions in aligned CO

_{2}molecules at the different driving laser intensity, shown in Fig. 4 (b) [16]. They also show that the spectral region of the destructive suppression (

*S/S*< 1) can be shifted to the higher harmonic order by enhancing the driving laser intensity. Therefore, the results of both two gases indicate that the two-center interference does not follow a relationship corresponding to a fixed molecular potential when the driving laser intensity is varied, which can be explained by our simple model.

_{0}#### 3.2 Discussion

We use the slightly longer *R* in simulations for two reasons: (1) the field distortion of the molecular HOMO cannot only change the effective molecular potential$\delta {I}_{p}$, but also stretch the effective distance *R* between the centers of electron density; (2) If $\delta {I}_{p}$is restricted to being 0, the simulation using the slightly longer *R* more fits the experimental data, although the same effect can be reached by further decreasing $\delta {I}_{p}$<0.

There is almost ‘no destructive effect’ visible in aligned O_{2} molecules due to that: (1) the destructive effect visible should appear at the much higher harmonic order, therefore one can’t see the effect visible at the low harmonic orders with decreasing intensity for O_{2} molecules; (2) learning from Fig. 2, the Signal Noise Ratio from O_{2} molecules is much smaller than that from CO_{2} molecules, due to that it’s harder to get an impulsive alignment for O_{2} molecules. Therefore, if decreasing the driving laser intensity, the Signal Noise Ratio is too small to identify the Signal under our experimental condition.

The dispersion relation could also be written as${E}_{k}=nh{v}_{0}-{I}_{p}+{}^{*}I{}_{p}$, where ${I}_{p}$ is the original molecular potential, and ${}^{*}I{}_{p}$is the new active molecular potential considering the field distortion. If the driving laser intensity is very low and the field distortion of HOMO can be neglected (i.e. ${I}_{p}={}^{*}I{}_{p}$), the dispersion relation will be close to${E}_{k}=nh{v}_{0}$ (i.e. $\delta {I}_{p}={I}_{p}-{}^{*}I{}_{p}=0$). This is consistent with the experimental observation that the lower driving laser intensity is close to *δI _{p} =* 0. Therefore, $\delta {I}_{p}={I}_{p}-{}^{*}I{}_{p}$ indicates the depth of the field distortion of the molecular potential, which should have a close relation to the driving laser intensity (i.e. they have the same trend —— the higher driving laser intensity, the larger depth$\delta {I}_{p}$). Whether $\delta {I}_{p}$>

*I*or $\delta {I}_{p}$<0, we think they are possible at the extreme laser condition.

_{p}The experimental and calculated results of these two kinds of the π_{g}-type molecules indicate that the general feature of modulation inversion is governed by the interference effect, but the spectral region of modulation inversion can be shifted by tuning the driving laser intensity. Therefore, the two-center interference in the recombination process can be manipulated and even inverted from enhancement to suppression by changing the laser intensity through control over the effective molecular potential. The actual physical effect is that the higher driving laser intensity, the deeper field distortion of the molecular HOMO (The magnitude of the coefficient *δ* shows the depth of the field distortion of the molecular HOMO). These observations are helpful for the future experiment on such molecular system and important for the theoretical study on this topic.

Recently, through varying the wavelength and intensity of the generating laser field, Wörner *et al.* [23] have showed that the minimum in aligned N_{2} molecules is nearly unaffected, whereas the minimum in aligned CO_{2} molecules shifts over more than 15eV. The minimum in aligned N_{2} molecules is a result of the orbital structure, whereas the minimum in aligned CO_{2} molecules is dynamical and results from the interference between the multiple orbitals. These results directly present that the laser intensity dependence in aligned N_{2} molecules is different from the laser intensity dependence in aligned π_{g}-type molecules.

## 4. Conclusions

In summary, we investigate the processes of the two-center interference in aligned O_{2} molecules by comparing with CO_{2} molecules, and observe the detailed characteristic enhancement and suppression in the measurement of the temporal evolution and the angular distribution of the harmonic intensity. For the both two gases, the angular distribution and the laser intensity dependence can be well explained by the modified interference model. Moreover, we demonstrate that the spectral region of the constructive enhancement in aligned O_{2} molecules can be shifted by tuning the driving laser intensity, which is also supported by the modified interference model. Finally, the effects of quantum interference and laser intensity on HHG from the π_{g}-type molecules are summarized.

## Acknowledgements

We acknowledge the support from National Natural Science Foundation (Grants Nos. 10734080, 60578049 and 10523003), National Basic Research Program of China (Grant No. 2006CB806000), Chinese Academy of Sciences, and Shanghai Commission of Science and Technology (Grants Nos. 06DZ22015, 0652nm005 and 07pj14091).

## References and links

**1. **P. B. Corkum, “Plasma perspective on strong field multiphoton ionization,” Phys. Rev. Lett. **71**(13), 1994–1997 (1993). [CrossRef] [PubMed]

**2. **K. C. Kulander, J. Cooper, and K. J. Schafer, “Laser-assisted inelastic rescattering during above-threshold ionization,” Phys. Rev. A **51**(1), 561–568 (1995). [CrossRef] [PubMed]

**3. **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] [PubMed]

**4. **K. Miyazaki, M. Kaku, G. Miyaji, A. Abdurrouf, and F. H. M. Faisal, “Field-free alignment of molecules observed with high-order harmonic generation,” Phys. Rev. Lett. **95**(24), 243903 (2005). [CrossRef] [PubMed]

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

**6. **S. Baker, J. S. Robinson, C. A. Haworth, H. Teng, R. A. Smith, C. C. Chirilã, M. Lein, J. W. Tisch, and J. P. Marangos, “Probing proton dynamics in molecules on an attosecond time scale,” Science **312**(5772), 424–427 (2006). [CrossRef] [PubMed]

**7. **H. J. Wörner, J. B. Bertrand, D. V. Kartashov, P. B. Corkum, and D. M. Villeneuve, “Following a chemical reaction using high-harmonic interferometry,” Nature **466**(7306), 604–607 (2010). [CrossRef] [PubMed]

**8. **E. Goulielmakis, Z. H. Loh, A. Wirth, R. Santra, N. Rohringer, V. S. Yakovlev, S. Zherebtsov, T. Pfeifer, A. M. Azzeer, M. F. Kling, S. R. Leone, and F. Krausz, “Real-time observation of valence electron motion,” Nature **466**(7307), 739–743 (2010). [CrossRef] [PubMed]

**9. **T. Kanai, S. Minemoto, and H. Sakai, “Quantum interference during high-order harmonic generation from aligned molecules,” Nature **435**(7041), 470–474 (2005). [CrossRef] [PubMed]

**10. **C. Vozzi, F. Calegari, E. Benedetti, J. P. Caumes, G. Sansone, S. Stagira, M. Nisoli, R. Torres, E. Heesel, N. Kajumba, J. P. Marangos, C. Altucci, and R. Velotta, “Controlling two-center interference in molecular high harmonic generation,” Phys. Rev. Lett. **95**(15), 153902 (2005). [CrossRef] [PubMed]

**11. **C. Vozzi, F. Calegari, E. Benedetti, R. Berlasso, G. Sansone, S. Stagira, M. Nisoli, C. Altucci, R. Velotta, R. Torres, E. Heesel, N. Kajumba, and J. Marangos, “Probing two-centre interference in molecular high harmonic generation,” J. Phys. At. Mol. Opt. Phys. **39**(13), 457 (2006). [CrossRef]

**12. **X. Zhou, R. Lock, W. Li, N. Wagner, M. M. Murnane, and H. C. Kapteyn, “Molecular recollision interferometry in high harmonic generation,” Phys. Rev. Lett. **100**(7), 073902 (2008). [CrossRef] [PubMed]

**13. **T. Kanai, E. Takahashi, Y. Nabekawa, and K. Midorikawa, “Observing molecular structures by using high-order harmonic generation in mixed gases,” Phys. Rev. A **77**(4), 041402 (2008). [CrossRef]

**14. **A. Le, X. Tong, and C. Lin, “Evidence of two-center interference in high-order harmonic generation from CO_{2},” Phys. Rev. A **73**(4), 041402 (2006). [CrossRef]

**15. **P. Liu, P. Yu, Z. Zeng, H. Xiong, X. Ge, R. Li, and Z. Xu, “Laser intensity dependence of high-order harmonic generation from aligned CO_{2} molecules,” Phys. Rev. A **78**(1), 015802 (2008). [CrossRef]

**16. **P. Wei, P. Liu, J. Chen, Z. Zeng, X. Guo, X. Ge, R. Li, and Z. Xu, “Laser-field-related recombination interference in high-order harmonic generation from CO_{2} molecules,” Phys. Rev. A **79**(5), 053814 (2009). [CrossRef]

**17. **W. Boutu, S. Haessler, H. Merdji, P. Breger, G. Waters, M. Stankiewicz, L. J. Frasinski, R. Taieb, J. Caillat, A. Maquet, P. Monchicourt, B. Carre, and P. Salieres, “Coherent control of attosecond emission from aligned molecules,” Nat. Phys. **4**(7), 545–549 (2008). [CrossRef]

**18. **B. K. McFarland, J. P. Farrell, P. H. Bucksbaum, and M. Gühr, “High harmonic generation from multiple orbitals in N_{2.},” Science **322**(5905), 1232–1235 (2008). [CrossRef] [PubMed]

**19. **O. Smirnova, Y. Mairesse, S. Patchkovskii, N. Dudovich, D. M. Villeneuve, P. B. Corkum, and M. Y. Ivanov, “High harmonic interferometry of multi-electron dynamics in molecules,” Nature **460**(7258), 972–977 (2009). [CrossRef] [PubMed]

**20. **G. Lee, I. Kim, S. Park, T. Kim, Y. Lee, and C. Nam, “Alignment dependence of high harmonics contributed from HOMO and HOMO-1 orbitals of N_{2} molecules,” J. Phys. B **43**(20), 205602 (2010). [CrossRef]

**21. **S. Haessler, J. Caillat, W. Boutu, C. Giovanetti-Teixeira, T. Ruchon, T. Auguste, Z. Diveki, P. Breger, A. Maquet, B. Carré, R. Taïeb, and P. Salières, “Attosecond imaging of molecular electronic wavepackets,” Nat. Phys. **6**(3), 200–206 (2010). [CrossRef]

**22. **C. Figueira de Morisson Faria and B. B. Augstein, “Molecular high-order harmonic generation with more than one active orbital: Quantum interference effects,” Phys. Rev. A **81**(4), 043409 (2010). [CrossRef]

**23. **H. J. Wörner, J. B. Bertrand, P. Hockett, P. B. Corkum, and D. M. Villeneuve, “Controlling the interference of multiple molecular orbitals in high-harmonic generation,” Phys. Rev. Lett. **104**(23), 233904 (2010). [CrossRef] [PubMed]

**24. **Y. Mairesse, J. Levesque, N. Dudovich, P. B. Corkum, and D. M. Villeneuve, “High harmonic generation from aligned molecules-amplitude and polarization,” J. Mod. Opt. **55**(16), 2591–2602 (2008). [CrossRef]

**25. **T. Morishita, A.-T. Le, Z. Chen, and C. D. Lin, “Accurate retrieval of structural information from laser-induced photoelectron and high-order harmonic spectra by few-cycle laser pulses,” Phys. Rev. Lett. **100**(1), 013903 (2008). [CrossRef] [PubMed]