Modulations to molecular high order harmonic generation by electron de Broglie wave
^{1}
Abstract
We present a new theory that the molecular high order harmonic generation in an intense laser field is determined by molecular internal symmetry and momentum distribution of the tunnelingionized electron. The molecular internal symmetry determines the quantum interference form of the returning electron inside the molecule. The electron momentum distribution determines the relative interference strength of each individual electron de Broglie wave. All individual electron de Broglie wave interferences add together to collectively modulate the molecular high harmonic generation. We specifically discuss the suppression of the generation on adjacent harmonic orders and the dependence of molecular high harmonic generation on laser intensities and molecular axis alignment. Our theoretical results are in good consistency with the experimental observations.
pacs:
42.65.Ky, 32.80.RmLasermolecule interaction has become an active research area in high field physics. The molecular multicenter structure strongly affects its interaction with the external laser field, resulting in abundant new physics(1); (2); (3); (4); (5); (6); (7); (10); (9); (11); (8); (12); (13). Recent studies show that the molecular orbital symmetry strongly affects the single ionization(1); (2); (3); the nuclei motion in molecular dissociation and Coulomb explosion can be monitored by studying the electron rescattering process(4); (5); the molecular orbital can be imaged by studying the angular distribution of coincident ion pairs in the molecular double ionization(6); by studying the spectral profile of molecular high order harmonic emission, one is able to reconstruct the tomography of the molecular structure(9); (11); (12); (13). Of particular interest is the molecular high order harmonic generation (HHG). The molecular multicenter structure provides the possibility to actively modulate the molecular HHG to generate a single attosecond pulse.
The ionization process of a molecule in an intense laser field is described using the wellknown simpleman’s model(14). First one molecular electron tunnels through the potential barrier assisted by the laser field, then the electron in the continuum state may return to the molecular ion as the laser field reverses its phase. Different from the rescattering process in the laseratom interaction, the multiple atomic sites inside the molecular ion provide multiple options for the electron recombination. The pathdifferences of the returning electron to the multiple atomic sites result in the interference of electron de Broglie wave () which significantly affects the molecular ionization and HHG. This was first noticed by Lein and his colleagues in studying the interaction of laser with molecules(15). Lein et al. showed that when recaptured into the bonding orbit of a molecule, the electron waves from the two atomic site emitters (with positions and ) interfere as,
(1) 
which results in constructive or destructive interference when () is a multiple or a multiple and a half of the laser wavelength. As a result, the order molecular harmonic emission (, ) is enhanced (or suppressed).
Recently two experiments were carried out independently to study the HHG of in an intense laser field with similar experimental configurations only differing in laser intensities. The suppression of HHG was observed for a range of adjacent harmonic orders in both experiments(16); (17); while the suppressed harmonic spectral ranges were displaced by about 8 harmonic orders with respect to each other. This cannot be explained by Eq. (1).
In this letter, we present a new theory based on the modified molecular Lewenstein model to describe the lasermolecular interaction. We show that similar to the Young’s twoslit interference, the returning electron (tunnelingionized in the external laser field) wave interferes on the multiple atomic sites inside the molecule ion. The returning electron generally has a momentum () distribution. Each single electron momentum state (the electron de Broglie wave) selfinterferes on the multiple atomic sites. All interferences of different electron momentum states () add together to affect the molecular ionization and HHG. The interference structure is uniquely determined by the molecular structure. The electron momentum distribution is a function of laser intensity and molecular axis alignment angle in the external laser field. Thus the molecular HHG changes as laser intensity or molecular axis alignment varies. We also show that molecular HHG at adjacent harmonic orders has similar dependence on electron momentum distribution, thus the suppression or enhancement on molecular HHG occurs over a broad range of adjacent harmonic orders. Our theoretical results agree well with experimental observations.
To directly compare with experimental observations, we study the HHG of which has antibonding orbital as the highestoccupied state and can be approximately considered as a twocenter system when it interacts with the external laser field(16); (17). Using the molecular Lewenstein model, the dipole acceleration of to yield the HHG radiation in an external laser field (polarized along ) is written as(18); (19); (20); (21),
(2) 
where and are vector potential and electric field of the laser pulse. Beside the action , two additional action terms and are introduced accounting for the effects induced by the two Oatomic sites. The three phase terms carrying these actions interfere with each other and modulate the molecular HHG(19); (20); (21); (22). is the Fourier transformation of the atomic wavefunction under linear combination of atomic orbitalsmolecular orbitals (LCAOMO) approximation in which the molecular wavefunction is expressed as .
After adopting pole approximation(20) and performing timeintegral in Eq. (2), the Fourier transformation of Eq. (2) gives the HHG spectral amplitude of (21); (22)
(3) 
where is the order Bessel function, and are numbers of photons that the electron absorbs or emits at ionization or recombination, with . The electron’s momentum is now given as , where is the ionization potential of the molecule. Considering the acceleration effect by bound potential when the electron is in the vicinity of the core, the modified electron momentum (parallel to ) is given as (22).
Eq. (3) shows that the molecular HHG is completely determined by the structural factor( ) which corresponds to the molecular internal symmetry (here it describes the antibonding orbital of ) and the momentum distribution [the rest part of Eq. (3)] of the rescattering electron which is a function of the laser intensity and molecular axis alignment angle (with respect to the direction of the external laser field). Eq. (3) also shows that the electron momentum distribution and consequently the molecular HHG spectral profile vary for different laser intensity and molecular alignment angle. In the following, we numerically study the HHG of under various experimental conditions and compare the results with available experimental data.
We first study the contribution to molecular HHG from electron momentum states with momentum parallel to , because they were considered to contribute the most to the HHG. As shown in Fig. 1, for both harmonic orders [Fig. 1(a)] and [Fig. 1(b)], contributions from individual electron momentum states (with parallel to ) are significantly suppressed for smaller but remain unchanged for big , when considering the interference effect that is due to the antibonding symmetry of . An integral of was then conducted to consider contributions from all electron momentum states. Both [Fig. 1(c)] and [Fig. 1(d)] are constant for big , indicating that the contribution from big electron momentum states is negligible. When considering the interference effect, is relatively suppressed while is relatively enhanced.
Next we numerically study the variation of HHG spectral profile of as a function of laser intensity and molecular alignment angle . For fixed laser intensity, the HHG spectral profiles vary for different values of as shown in Figs. 2(a) and (b) in which the plots are generated with and , respectively. Eq. (3) shows that harmonic emissions at adjacent orders have similar dependence on electron momentum contribution. Consequently similar modulations occur on adjacent HHG orders. As shown in Fig. 2(c), with and , HHG is suppressed at a range of adjacent harmonic orders centering at . The center shifts to for [Fig. 2(d)], qualitatively agreeing with recent experimental observations(16); (17). (Intensity stabilization and calibration of intense ultrashort laser pulse remain to be a challenge in the experiment.)
To clearly exhibit the dependence of molecular HHG on the laser intensity, we numerically study the HHG at two harmonic orders. For , our numerical results show that the harmonic generation is suppressed at [Fig.3(a)] and enhanced at [Fig.3(b)]. For , we show the suppression of harmonic generation at [Fig.3(c)] and enhancement at [Fig.3(d)].
Lastly, it is interesting to find out under what kind of condition the single electron momentum picture given by Lein et al. may be used in predicting the suppression position in the HHG spectrum. Examining Eq. (2) shows that if (describing the molecular structure) is much smaller than the action (corresponding to the electron trajectory), then can be extracted out as a common factor.
The Fouriertransformation of Eq. (2) can now be simplified into Eq. (1). Assuming , we have , where is the period of the laser field and is the average kinetic energy of the electron with and . Then we reach such a relation
(4) 
where is the electron’s quiver amplitude in the laser field. Eq. (4) shows that if the electron’s quiver amplitude is much larger than the molecular internuclear distance, or in other words, the time difference for the rescattering electron to recombine with the different atomic sites are negligible comparing to the time for electron being in the continuum state, then Eq. (1) can be used to approximately describe the molecular HHG process. Comparing the internuclear distances for and for with the electron’s quiver amplitude of in a laser field with intensity of , it is easy to understand that Eq. (1) gives a better estimation to the suppression of HHG of (15) than to the suppression of HHG of .
In conclusion, we develop a new theory about the molecular HHG in an intense laser field. We illustrate that the molecular HHG at each order is contributed by a group of interference effects on the Young’s slit (multipleatomic sites) with each produced by an independent electron de Broglie wave (eigen electron momentum state) and weighed by the electron momentum state amplitude. We show that the single electron momentum picture illustrated by Lein et al. is only an approximation of our theory at high laser intensity where the electron’s quiver motion amplitude is much larger than the molecular internuclear distance. Our theoretical analysis shows that the difference in the observed suppressed HHG spectral ranges in recent experiments is attributed to the different laser intensities applied in the measurement. We further suggest that the molecular HHG in an intense laser field may be generally described in the form of Eq. (3) by simply replacing the interference factor with the appropriate structure factor of the molecule in use.
This work was supported by the National Natural Science Foundation of China under Grant No. 10574019, 973 research program No. 2006CB806000, ICF research fund under Grant No. 2004AA84ts08, CAEP Foundation No. 2006z0202 and the Multidisciplinary University Research Initiative Center for Photonic Quantum Information Systems (Army Research Office/DTO program DAAD190310199).
Footnotes
 footnotetext: chenjing@iapcm.ac.cn; jingyun.fan@nist.gov
References
 Talebpour A, Chien C Y and Chin S L 1996 J. Phys. B 29 L677
 Guo C, Li M, Nibarger J P and Gibson G N 1998 Phys. Rev. A 58 R4271
 MuthBöhm J, Becker A and Faisal F H M 2000 Phys. Rev. Lett. 85 2280
 Niikura H et al 2003 Nature 421 826
 Alnaser A S et al 2003 Phys. Rev. Lett. 91 163002
 Alnaser A S et al 2004 Phys. Rev. Lett. 93 113003
 Litvinyuk I V et al 2003 Phys. Rev. Lett. 90 233003
 Zeidler D et al 2005 Phys. Rev. Lett. 95 203003
 Itatani J et al 2004 Nature 432 867
 de Nalda R et al 2004 Phys. Rev. A 69 031804(R)
 Itatani J et al 2005 Phys. Rev. Lett. 94 123902
 Levesque Jet al 2007 Phys. Rev. Lett. 98 183903
 Le VanHoang et al 2007 Phy. Rev. A 76 013414
 Corkum P B 1993 Phys. Rev. Lett. 71 1994

Lein M et al 2002 Phys. Rev. Lett. 88
183903
Lein M et al 2002 Phys. Rev. A 66 023805
Lein M et al 2003 Phys. Rev. A 67 023819  Kanai T, Minemoto S and Sakai H 2005 Nature 435 470473

Vozzi C et al 2005 Phys. Rev. Lett. 95
153902
Vozzi C et al 2006 J. Phys. B 39 S457  Le AnhThu, Tong X M and Lin C D 2006 Phy. Rev. A 73 041402(R)
 Chirila C C and Lein M 2006 Phys. Rev. A 73 023410
 Chen J, Chu ShihI and Liu J 2006 J. Phys. B. 39 47474758
 Chen Y J, Chen J and Liu J 2006 Phys. Rev. A. 74 063405
 Chen Y J et al Model analysis of twocenter interference on high harmonic generation, unpublished.