Laser-induced blurring of molecular structure information in high harmonic spectroscopy

High harmonic spectroscopy gives access to molecular structure with Angström resolution. Such information is encoded in the destructive interferences occurring between the harmonic emissions from the different parts of the molecule. By solving the time-dependent Schrödinger equation, either numerically or with the molecular strong-field approximation, we show that the electron dynamics in the emission process generally results in a strong spectral smoothing of the interferences, blurring the structural information. However we identify specific generation conditions where they are unaffected. These findings have important consequences for molecular imaging and orbital tomography using high harmonic spectroscopy.

High harmonic spectroscopy gives access to molecular structure with Angström resolution. Such information is encoded in the destructive interferences occurring between the harmonic emissions from the different parts of the molecule. By solving the time-dependent Schrödinger equation, either numerically or with the molecular strong-field approximation, we show that the electron dynamics in the emission process generally results in a strong spectral smoothing of the interferences, blurring the structural information. However we identify specific generation conditions where they are unaffected. These findings have important consequences for molecular imaging and orbital tomography using high harmonic spectroscopy.
High Harmonic Spectroscopy (HHS) is a powerful technique in which the process of High Harmonic Generation (HHG) is used to probe the structure and dynamics of the generating medium with Angström and attosecond resolution. In HHG, an Electron Wave-Packet (EWP) is prepared by tunnel ionization of the atomic/molecular target gas, it is then accelerated by the strong laser field and finally driven back to the core 1,2 . In the recombination process, a burst of extreme ultraviolet photons is emitted, encoding a wealth of information on the target. First, the recombination dipole moment imprints the structure of the molecular orbital involved in the emission, under the form of structural interferences 3 equivalent to Cohen-Fano interferences in the photo-ionization dipole 4 . The accurate spectral position and shape of the signature of these interferences in the HHG spectrum (amplitude and phase) is the required information for performing molecular imaging 3 and orbital quantum tomography [5][6][7] . Second, any dynamics occurring in the core during the EWP continuum excursion strongly affects the recombination. For instance, tunnel ionization from different valence orbitals results in interfering channels in the recombination [8][9][10] .
However, it is crucial to disentangle the different effects in order to access this rich information [11][12][13] . This was generally performed by recording the dependence on the laser parameters (intensity, wavelength), since it was widely accepted that structural interferences do not depend on them 14 . This appears clearly in the Quantitative ReScattering theory (QRS) 15 where the harmonic dipole is expressed as the product of the returning EWP, containing all the laser parameters, and the target-specific field-free recombination dipole, only dependent on the EWP energy and recollision direction. Such a factorization has been very successful in predicting and explaining resonant features in the harmonic intensity, e.g., in 16 . It was further confirmed by analytic calculations in the case of a model potential 17 .
Here, we study the signatures of the electron dynamics on the structural interference occurring in the harmonic emission from diatomic molecules. Based on Time-dependent Schrödinger Equation (TDSE) calculations resolving the two shortest EWP trajectories for each harmonic, we show that the value of the laser field at recombination plays a crucial role on the shape of the two-center destructive interference. In most cases, the intensity minimum and the associated phase-jump are spread over a broad spectral range resulting in a very smooth behavior, and the phase-jumps for the short and long EWP trajectories are in opposite directions. Furthermore, we identify very peculiar situations for which this phase-jump is very sharp, providing ideal conditions for performing molecular imaging and quantum tomography. Using the Strong-Field Approximation (SFA) approach 18 modified for molecules 19,20 , we explain these features by the modification of the saddle-point trajectories induced by the fast phase variation of the recombination dipole. Our findings thus question the strict separation between continuum dynamics and recombination advocated, e.g., in QRS, with strong implications for HHS.
We chose to perform our study on a one dimensional (1D) system that contains the essential physics of two-center interference in HHG. This has two advantages: (i) we can easily achieve numerous computations both within the TDSE and the SFA framework, and implement analytical developments; (ii) we avoid additional effects like Coulomb refocusing or EWP spreading that was previously invoked to explain phase smoothing effects 21,22 . Note that we carried out preliminary 2D computations that are in overall agreement with our 1D results, eventhough they present extra features, that we plan to address in a forthcoming paper. Our model system represents a diatomic molecule with fixed nuclei, under the single-active electron approximation.
First, we solved the TDSE to simulate the "exact" electron dynamics under the influence of a low-frequency strong laser field. The electron-nuclei interaction is represented by: where R is the internuclear distance. This is valid for molecules with "heavy" nuclei that exhibit a strong two-center character, such as CO 2 7,23,24 . In order to study the R-dependence of the structural interferences, the regularization parameter a is adapted to maintain the same ground state energy (−I p = −0.567 a.u = −15.43 eV) at each considered value of R. The electric field E(t), of amplitude E L , has a sine-square envelop lasting two optical cycles with frequency ω L corresponding to a 800-nm Ti:sapphire laser. The carrier -envelope phase is set equal to zero in order to obtain the contributions of only one set of short and long trajectories in the HHG spectrum. The t → ω Fourier transform of the acceleration 25 provides the harmonic field. We calibrate the latter with the one calculated for a reference "atom" (R = 0 in Eq. 1) with same ionization potential, i.e., we divide the harmonic field calculated for the molecule by the one simulated in the reference atom. This standard procedure in HHS allows removing both the plateau-cutoff shape in the spectrum and the group delay dispersion associated with the attochirp 26 , in order to evidence the two-center signatures we are interested in. We numerically discriminate the short and long trajectory contributions with an absorber at a distance E L /ω 2 , using the fact that the short trajectories never go above this limit while the long ones always cross it 27 . Figure 1 reports the calibrated harmonic intensity and phase for several internuclear distances between 1.4 and 1.7 a.u., and a laser peak intensity of 2.55 × 10 14 W.cm −2 . As expected 3 , the possibility for the electron to recombine on either center results in a destructive two-center interference that appears as a minimum in the spectrum accompanied by a ≈ π phase-jump. Its position moves to high harmonic orders with decreasing R, since it requires a shorter EWP de Broglie wavelength λ dB (R ≈ λ dB /2). However, we observe in Fig. 1 many remarkable unexpected features: (i) for most R values, the intensity minimum and associated phase-jump are spread over a large spectral range; (ii) the phase-jumps for the short and long trajectories are in opposite directions for high values of R; (iii) the shape of the interference strongly depends on R. In particular for the long trajectory, at a critical internuclear distance R c ≈ 1.55 a.u., the interference becomes extremely destructive at harmonic 41 (H41) with a very steep π phase-jump. For smaller R values, the interference becomes smooth again, but with a change of sign of the phase-jump that now resembles that of the short trajectory. Indeed, the interference is now positioned in the cutoff region (>H41) where the two trajectories coalesce.
Within the QRS, the calibrated emission should be directly proportional to the molecular recombination dipole matrix element. We computed it with the numerically exact scattering waves of the potential V(x) and found an excellent agreement only at the critical distance R c and for the long trajectory.
In order to shed light on the above features, we performed a time-frequency analysis of the harmonic dipole using Gabor transforms as in 28 . We retrieved the emission times of the harmonics for which the destructive interferences occur and found that, at the particular internuclear distance R c , it corresponds to an almost zero instantaneous electric field for the long trajectory. Furthermore, larger (lower) values of R lead to recollision with a negative (positive) electric field for the long trajectories while it is always positive for the short trajectories. Thus, the sign and shape of the phase-jump seem to be strongly correlated to the instantaneous value of the electric field when the harmonics are emitted.
To further examine the influence of the electric field, we varied the peak laser intensity I L from 2 to 4 × 10 14 W.cm −2 at fixed R = 1.4 a.u. and show the results in Fig. 2(a,b). While the position of the destructive interference is relatively independent of I L (i.e., around H49), the phase-jump behaves similarly when varying R at fixed I L . We found a critical intensity I c = 3.24 × 10 14 W.cm −2 for which we observe the inversion of direction with a very sharp πphase-jump. Here again, this value directly corresponds to an instantaneous electric field ≈ 0 when H49 is emitted.
We carried out the same study within the molecular SFA framework, derived following 20 , the results of which are reported in Fig. 2(c,d). The same conclusions concerning the phase-jump can be drawn, with minor differences on the position of the destructive interference and the laser intensity range. These differences are due to the way the continuum is described in SFA [plane-wave (PW) approximation], and also to the position of the cutoff beyond which harmonics cannot be computed within SFA. These SFA computations reproduce remarkably well the smoothing observed in the TDSE simulations. This evidences that the main origin of the smoothing is not related to Coulomb scattering effects 29 . A major advantage of the SFA lies in its semi-analytical form. It allows for a semi-classical analysis by searching the saddle-points of the action that determine the main quantum orbits/  trajectories contributing to HHG 18,30 . Thus, we performed a Taylor expansion of the stationary solutions for the molecular SFA around the ones for the reference atom to get physical insight into the electron dynamics. Following 19,20 , we approximate the symmetric electronic ground state of a homonuclear diatomic molecule by the sum of two symmetric atomic orbitals centered on each nucleus. The PW-recombination and -ionization dipole matrix elements become proportional to cos(pR/2), with = Ω − p I 2( ) p the continuum electron momentum associated with each harmonic frequency Ω. By expressing cos(pR/2) in terms of exp[±ipR/2], the action reads: , 2} label the atomic centers and A is the vector potential associated with the electric field. The first two terms are identical to the atomic SFA while the last two are induced by the fast phase variations of the molecular ionization and recombination dipoles. The saddle-point equations associated with this modified action 19,20 allow to determine the stationary solutions for the momentum p stat , the ionization and recombination times ′ t stat and t stat respectively. Compared to their atomic equivalents, we found that p stat is unchanged to first order in the molecular case, while the times ′ t stat and t stat differ by: respectively. The physical meaning of these quantities is rather straightforward. They correspond to an increase (decrease) of either the barrier to tunnel through or the path in the continuum, for the electron, whether the nucleus is located at −R/2 (+R/2) for j or k equal 1 or 2, respectively. As the former is classically forbidden, the change in ionization time is purely imaginary, while the latter, classically allowed, is purely real. The saddle-point harmonic dipole can then be factorized in the spirit of the QRS but with a modified recombination dipole, which reads to first order: rec s tat where ζ is a constant and  p ( ) is proportional to the overlap of the atomic orbitals with the PWs representing the continuum 27 . One clearly sees that in addition to the "bare" recombination dipole (first term in Eq. (4), here purely imaginary),  d rec contains a real valued part varying like sin(pR/2) and depending on the electric field at recombination time t. We have checked numerically that ζ, which originates from the rigorous handling of the prefactors arising from the saddle-point approximations within SFA, is almost constant and small compared to the other term in the bracket 27 . It explains why the discontinuous π-jump is observed at an almost but not exactly zero electric field at recombination time.
Therefore, (i)  d rec is not restricted to the imaginary axis across the HHG spectra and travels in the complex plane, (ii) it does not vanish when the two-center interference term, proportional to cos(pR/2), is zero, except when the bracket in Eq. (4) is also zero, (iii) it has a real part which is almost proportional to E(t) (ζ is significantly small with respect to the other term in the brackets 27 ) and which value and sign drives the direction and the smoothness of the phase-jump. Thus, =  d 0 rec happens only when E(t) ≈ 0 at emission time of H49 of the long trajectory in the conditions of Fig. 2, confirming our findings from both the TDSE and the full molecular SFA simulations.
One should note that, if we neglect ζ in Eq. (4), we recover an expression for  d rec close to the one found by including the ad hoc dressing of the molecular ground state by the electric field, in the light of ref. 31 . Indeed, in the presence of a laser field, we may approximate the dressed fundamental wave-function as a linear combination: 0 g u of a symmetric state σ g (field-free molecular ground-state), sum of two symmetric atomic orbitals centered at ±R/2, and an antisymmetric state σ u (field-free first excited state), difference of the same two atomic orbitals 32 . This decomposition is a very good approximation even at short distance R 27 . Within the time-dependent perturbation theory, c(t) depends linearly on the instantaneous electric field E(t), i.e., c(t) = χE(t), and we checked, by computing the exact dressed ground-state 33 , that this is still valid for fields strength ≈10 14 W.cm −2 . Thus, the modified recombination dipole associated with Φ 0 within the PW approximation, becomes: rec which resembles the expression given in Eq. (4). Thus, the signature of the electron dynamics in  d rec was attributed to the manifestation of the dressing of the electronic ground-state.
In practice, the influence of the electron dynamics on the structural interference could be observed by following the phase difference between the short/long trajectory contributions when varying the laser intensity, e.g., by recording the Quantum Path Interferences in the total harmonic dipole 34,35 .
In summary, we have shown that the electron dynamics in the HHG process may blur the two-center interferences encoding the molecular orbital structure. The value and sign of the driving laser field at recombination become critical in spreading the intensity minimum and associated phase-jump over a large spectral range and inducing opposite phase-jumps for the short and long trajectories. Except in specific situations, our study demonstrates that the harmonic emission is not entirely proportional to the field-free recombination dipole moment. This raises questions since such factorization is central in HHS, e.g., in the QRS theory. This should be taken into account when comparing the recombination dipole extracted from HHG to the photoionization dipole, or using it for molecular imaging and quantum tomography. Our conclusions could be generalized to other types of structures occurring in the recombination dipole, like shape or autoionization resonances. If narrow enough, the corresponding fast phase variations could modify the electron dynamics and result in a smoother behavior in the harmonic emission. Further studies are in progress to uncover all the implications of the reported findings.