# Multi-channel electronic and vibrational dynamics in polyatomic resonant high-order harmonic generation

A. Ferré, A. E. Boguslavskiy, M. Dagan, Valérie Blanchet, B. D. Bruner, F. Burgy, A. Camper, D. Descamps, B. Fabre, N. Fedorov, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

A. Ferré, A. E. Boguslavskiy, M. Dagan, Valérie Blanchet, B. D. Bruner, et al. Multi-channel electronic and vibrational dynamics in polyatomic resonant high-order harmonic generation. Nature Communications, 2015, 6, pp.5952. 10.1038/ncomms6952. cea-01244391

HAL Id: cea-01244391

https://cea.hal.science/cea-01244391

Submitted on 7 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

Received 28 May 2014 | Accepted 25 Nov 2014 | Published 22 Jan 2015

DOI: 10.1038/ncomms6952

OPEN

# Multi-channel electronic and vibrational dynamics in polyatomic resonant high-order harmonic generation

A. Ferré<fnref n="1" />, A.E. Boguslavskiy²,³, M. Dagan⁴, V. Blanchet<fnref n="1" />, B.D. Bruner⁴, F. Burgy<fnref n="1" />, A. Camper⁵, D. Descamps<fnref n="1" />, B. Fabre<fnref n="1" />, N. Fedorov<fnref n="1" />, J. Gaudin<fnref n="1" />, G. Geoffroy<fnref n="1" />, J. Mikosch²,⁶, S. Patchkovskii²,⁶, S. Petit<fnref n="1" />, T. Ruchon⁵, H. Soifer⁴, D. Staedter⁷, I. Wilkinson², A. Stolow²,³,⁸, N. Dudovich⁴ & Y. Mairesse<fnref n="1" />

High-order harmonic generation in polyatomic molecules generally involves multiple channels of ionization. Their relative contribution can be strongly influenced by the presence of resonances, whose assignment remains a major challenge for high-harmonic spectroscopy. Here we present a multi-modal approach for the investigation of unaligned polyatomic molecules, using SF₆ as an example. We combine methods from extreme-ultraviolet spectroscopy, above-threshold ionization and attosecond metrology. Fragment-resolved above-threshold ionization measurements reveal that strong-field ionization opens at least three channels. A shape resonance in one of them is found to dominate the signal in the 20-26 eV range. This resonance induces a phase jump in the harmonic emission, a switch in the polarization state and different dynamical responses to molecular vibrations. This study demonstrates a method for extending high-harmonic spectroscopy to polyatomic molecules, where complex attosecond dynamics are expected.

<footnote n="1">¹Université de Bordeaux—CNRS—CEA, CELIA, UMR5107, F33405 Talence, France. ²National Research Council Canada, 100 Sussex Drive, Ottawa, K1A 0R6 Ontario, Canada. ³Department of Physics, University of Ottawa, 150 Louis Pasteur, Ottawa, K1N 6N5 Ontario, Canada. ⁴Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot 76100, Israel. ⁵CEA, IRAMIS, Lasers, Interactions and Dynamics Laboratory—LIDyL, CEA-SACLAY, F-91191 Gif-sur-Yvette, France. ⁶Max-Born-Institute, Max-Born-Strasse 2A, 12489 Berlin, Germany. ⁷Université de Toulouse—CNRS, LCAR-IRSAMC, 31062 Toulouse, France. ⁸Department of Chemistry, University of Ottawa, 10 Marie Curie, Ottawa K1N 6N5 Ontario, Canada K1N 6N5. Correspondence and requests for materials should be addressed to Y.M. (email: mairesse@celia.u-bordeaux1.fr).</footnote>
In the simplest view of molecular strong-field ionization, the most probable valence electron to be removed by a strong laser field is the most weakly bound one, namely that from the highest occupied molecular orbital. This assumption is based on quasi-static, single active electron pictures and the exponential decay of the tunnel ionization probability through a field-lowered Coulomb barrier¹. However, even in a quasi-static picture, the tunnel ionization probability also depends on the spatial structure of the molecular orbital²,³. Furthermore, more than one electron could get driven coherently by the laser field. The timescales of electron motions in molecules may not be fast compared with the laser period and the adiabatic approximation inherent to the quasi-static picture may fail: non-adiabatic multi-electron dynamics can ensue⁴. Consequently, strong-field ionization is significantly more complex for molecules than for atoms, with possible contributions from multiple channels, leaving the ion in a coherent superposition of electronic states. Such ionic excitation can result in multi-electron processes and ultra-fast charge migrations on an attosecond timescale⁵,⁶.

The signatures of multi-channel ionization in above-threshold ionization (ATI)⁷ and high-harmonic generation (HHG)⁸⁻¹⁰ are fundamentally different. The ATI process for each ionization channel can be understood as a coherent superposition of continuum photoelectron wavepackets, which are created during every half cycle of the laser field. However, the ionization pathways leading to different ionization continua are distinguishable and will hence simply add up incoherently. ATI is a single-channel interferometer, different ionization channels produce independent ATI progressions in the photoelectron kinetic energy spectrum. In contrast, in HHG, the different ionization channels connect the same initial and final states (that is, neutral ground electronic state) and therefore their contributions are coherently summed in the harmonic emission. The HHG mechanism can thus be described as a multi-channel interferometer on the timescale of half a laser cycle. The weights of the interfering ionization channels are mainly determined by the ionization yields and the recombination dipole moments¹¹,¹².

In many molecular systems, the recombination process can involve resonant transitions due to the presence of autoionizing states or shape resonances. In that case, the recolliding electron can get trapped in the vicinity of the parent ion for typically several hundreds of attosecond before recombining to the ground state¹³,¹⁴. On the contrary to autoionization, shape resonances have an energy bandwidth—typically a few eV—that depends strongly on the non-spherical and non-local characters of the ionic potential. Consequently, their energy and lifetime can vary as function of the molecular geometry. These scattering patterns related to the molecular field enhance the photoabsorption cross section and are imprinted in the photoelectron asymmetry parameter. In HHG, the basic properties of the harmonic—their amplitude, phases and polarization state—may be dramatically influenced by such peculiar resonant electron scattering. Their presence may significantly increase the contribution of inner ionization channels, commonly hidden in HHG experiments. Despite their importance in molecular HHG, very few experiments have reported signatures of resonances. A maximum in the harmonic spectrum from xenon atoms has been attributed to a giant resonance¹⁵, but other experiments have shown that the shape of the harmonic emission in xenon critically depended on focusing conditions, such that maxima could appear or disappear in the spectrum¹⁶. Clearly, the identification of the resonances' contribution requires the measurement of additional experimental parameters. Following this approach, Haessler et al.¹⁷ reported the measurement of a phase jump across a narrow resonance in ablation plasma plumes. An alternative is to add control parameters to the interaction and

measure the response of the system to these parameters. For instance, molecular alignment is broadly used in HHG experiments for channel identification in diatomic and linear triatomic molecules⁵,⁶,¹⁸⁻²⁰. It has allowed to interpret a maximum in the intensity from aligned $N_2$ as a signature of a shape resonance in the lowest ionization channel²¹. However, for polyatomic molecules, field-free three-dimensional alignment remains challenging in general²². Therefore, we need a more direct control of the interaction.

In this letter, we combine complementary strong-field measurements to develop a multi-modal high-harmonic spectroscopy. To illustrate this, we chose a polyatomic system in which many channels can contribute to the harmonic emission: $SF_6$, where the four lowest states of the cation (each multiply degenerate) lie within less than 3 eV as illustrated in Fig. 1. Each cation state of $SF_6$ can be described in a first approximation by a single-hole configuration of the neutral molecule (Dyson orbital). Therefore, each subshell ionization channel corresponds to one cation state only and we will label the channels via their associated cation state. Due to its symmetry ($O_h$ point group), this system cannot be spatially laser aligned in the laboratory frame. The identification of the channels contributing to HHG and the measurement of multi-channel ionization dynamics is thus particularly challenging. To discern the ionization mechanism, which is the first step involved in the harmonic emission, we performed fragment-resolved ATI measurements²³ on $SF_6$ for laser wavelengths ranging from 800 to 400 nm, revealing that multiple channels contribute to the strong-field ionization. Using high-harmonic spectroscopy, we performed a complete characterization of the emission—spectrum, spectral phase, polarization state—and used the laser ellipticity and vibrational excitation as control parameters. This approach reveals a clear transition between resonant and non-resonant HHG around 26 eV. Vibrational dynamics in the ground electronic state of $SF_6$ (via stimulated Raman pumping) is a molecular control parameters that allows to explore the dependency of the resonance on the molecular geometry. The results of these combined experiments, corroborated by state-of-the-art ab initio calculations, allowed us to identify the inner ionization channels

<figure><img src="image_3.png" /><figcaption>**Figure 1 |** Schematic view of multi-channel HHG. Multiple open channels in strong-field ionization are identified with fragment-resolved ATI measurements. The field-free, vertical ionization potential of each channel is indicated next to the schematic view of the associated valence subshell of ionization. The radiative recombination process is modified by the presence of resonances, increasing some of the channels' contribution to the harmonic emission.</figcaption></figure>
as the dominant contributors to HHG in SF₆ and the prominent role of shape resonances in this mechanism.

## Results

Figure 1 shows the different ionization channels connecting the initial and final states, as well as their energies and symmetries. This figure also summarizes the different experiments described in the following.

**Strong-field ionization channels in SF₆.** ATI and HHG are closely related subcycle ionization phenomena. In ATI, strong-field ionization (SFI) releases electrons into the continuum around every maximum of the laser electric field. The emission is repeated every optical cycle throughout the pulse. As a result, the photoelectron spectrum shows discrete peaks separated by the photon energy ħω of the driving laser. The absolute position of the peaks of this ATI comb is specific to each ionization channel:

$$KE_e^j = nħω - I_p^{(s)} - U_p,$$

as it is determined by the Stark shifted (s) ionization potential I_p(s) of ionization channel j and the channel-independent ponderomotive shift U_p (n is an integer). At a constant intensity, the ATI combs of different channels give a direct measurement of their relative ionization potentials modulo the photon energy. In polyatomic molecules, one can exploit field-free fragmentation to better identify the different SFI channels. Using covariance or coincidence detection, the ATI spectra associated with different fragments are measured. How well different channels can be separated depends on the details of the unimolecular fragmentation process. In cases where specific fragments or a parent ion single out individual ionization continua, fragment-resolved ATI measurements become channel-resolved (the CRATI technique)²³,²⁴.

The position of an ATI comb in the photoelectron spectrum reflects the channel I_p(s) modulo the photon energy. In general, if the energy difference between two ionization channels matches the photon energy of the driving laser, the ATI combs will overlap. For SFI of SF₆, the situation is complicated by the fact that the states involved are significantly Stark shifted in the laser field. Moreover, these Stark shifts depend sensitively on the wavelength. We have estimated the ionization channel-specific dynamic Stark shifts computationally (see Supplementary Fig. 1). Importantly, if there is a strong resonance for a specific ionization channel, the Stark shifts might also diverge. As seen from Supplementary Fig. 1, there are channel-specific divergences for photon energies around 1.2, 1.4 and 2.7 eV, which corresponds to the energetic separation of the field-free X and A/B, A/B and C and X and C states of the SF₆ cation. For photon energies in the vicinity of such a Stark shift divergence, the respective comb might be washed out entirely due to averaging over the laser bandwidth.

To resolve different SFI channels, we varied the laser wavelength between 400 and 800 nm and performed covariance measurements, producing photoelectron spectra associated with the various fragments of SF₆ (see Methods). In Fig. 2, we show spectra for the most abundant fragment SF₅⁺ (note that the ground state of SF₆⁺ is unstable). In the 800-nm photoelectron spectrum displayed in Fig. 2a, a single ATI comb is observed. However, if the photon energy is increased, the single comb starts to split (Fig. 2b-d), and at 400 nm three progressions of peaks, each spaced by the photon energy, are clearly visible (Fig. 2e). This demonstrates that at least three ionization continua contribute to the SFI of SF₆ and that their yields must be comparable. For the 800-nm driving field, the multi-channel nature of the strong-field ionization process is obscured by the

<figure><img src="image_4.png" /><figcaption>**Figure 2 | SF₅⁺ fragment photoelectron spectra from SFI of SF₆, at different laser wavelengths and intensities (note that the ground state of SF₆⁺ is unstable).** While a single ATI comb is observed for 800 nm (a), it splits up into two (b,d) and ultimately three (e) progressions for shorter wavelengths. The different ATI combs are marked by the sets of vertical lines, which are spaced by the respective photon energy. It can hence be concluded that multiple ionization continua with similar yields must be present in the SFI of SF₆. For 800 nm, the multi-channel nature is obscured by the cation electronic states having energy level spacings, which happen to be multiples of the photon energy. (f) Calculated ionization yield for different channels at 800 nm. The X, A, B and C states are the ground and first three electronically (doublet) excited states of the cation.</figcaption></figure>

Stark shifted cation electronic states having ionization potentials congruent modulo the photon energy. The variation of the driving field frequency was required to resolve the multiple SFI channels involved.

An assignment of the ionization channels in Fig. 2 would require more reliable Stark shift calculations and a precise intensity calibration, since Stark shifts are in the eV range. This is
beyond the scope of the present work. To quantify the relative yields of the active SFI channels, we would require knowledge of the branching ratios for ionic fragmentation from each of the populated cation states, which are unknown, along with a detailed analysis of the photoelectron spectra associated with all ionic fragments. Instead we present computationally determined channel-specific SFI yields for 800 nm by calculating the tunnel ionization probability $P_{\text{ion}}^j$ for each channel $j$ by using a previously described approach$^{25,26}$ and averaging over all molecular orientations. Ionization calculations used a $\sin^2$ pulse matching the peak of a 1/2-cycle 800-nm pulse at $10^{14} \text{ W cm}^{-2}$ peak intensity. The resulting yields are shown in Fig. 2f. Three main factors determine the channel SFI rates: (i) the ionization potential of the channel—the ionization rate decreases exponentially with $I_p$; (ii) the topography of the ionizing molecular orbital—nodes in the plane perpendicular to the laser polarization tend to suppress the ionization probability; (iii) interchannel coupling, due to both the laser-induced dipole transitions and electron correlation—subcycle electronic transitions within the ion core can redistribute the final ionization amplitude between coupled channels. The first effect was long considered to be dominant: SFI mainly occurs from the highest occupied molecular orbital, even though deeper channels can be observed under specific conditions, such as in aligned molecules. However, in a CRATI study of *n*-butane$^{23}$, it was shown that excited-state channels dominate over the ground-state channel, despite their higher $I_p$s. The results from Fig. 2f are alike due to geometrical effects and interchannel coupling in SFI of SF$_6$ which overcome the exponential dependence on $I_p$, even for unaligned targets.

**High-harmonic spectrum.** SFI involving multiple channels generates multiple interfering pathways in HHG. The harmonic signal is determined by the coherent sum over all possible channels (see Fig. 1). Figure 3a shows an experimental harmonic spectrum generated in SF$_6$ using 800 nm pulses at $I \approx 1.3 \times 10^{14} \text{ W cm}^{-2}$. To distinguish the simple, single-orbital non-resonant picture from the multi-channel case, we compared our measurement with argon, an atomic system with a similar ionization

<figure><img src="image_5.png" /><figcaption>**Figure 3** | HHG from SF<sub>6</sub>. (a) Harmonic spectra generated in Ar (grey) and SF<sub>6</sub> (black) at *I* ≈ 1.3 × 10<sup>14</sup> W cm<sup>-2</sup>. The argon spectrum is vertically shifted to improve visibility. (b) Harmonic spectra obtained after normalizing out the recolliding electron wavepacket. Experiment (crosses) and contribution from different channels.</figcaption></figure>

potential to the first $I_p$ of SF$_6$. The SF$_6$ spectrum extends to higher orders and shows a pronounced minimum at harmonic 17, as was previously observed at a lower laser intensity$^{27}$. Varying the laser intensity in our experiment does not modify the position of this minimum.

The intensity of the harmonic emission at frequency $\omega$ for channel $j$ is determined by the SFI probability $P_{\text{ion}}^j$, the energy spectrum of the recolliding electron wavepacket in the continuum $S_{\text{elec}}^j(\omega)$ and the recombination cross section $d_{\text{rec}}^j(\omega)^{11,12}$: $I_j(\omega) = S_{\text{elec}}^j(\omega) |d_{\text{rec}}^j(\omega)|^2$. We assume that the recolliding electron wavepacket spectrum has the same energy structure for all channels as in argon since they have similar ionization potentials: $S_{\text{elec}}^j(\omega) = S^{\text{Ar}}(\omega) \propto I^{\text{Ar}}(\omega) / (|d_{\text{rec}}^{\text{Ar}}(\omega)|^2)$. Using the measured harmonic spectrum and the experimental photoionization cross section of Ar from ref. 28, we obtain the wavepacket spectrum. Normalizing the harmonic spectrum of SF$_6$ by $S_{\text{elec}}(\omega)$ provides a quantity which is only related to ionization probabilities and recombination cross-sections: $\tilde{I}(\omega) = P_{\text{ion}} |d_{\text{rec}}(\omega)|^2$. This normalized spectrum shows a maximum at harmonic 15 and a minimum at harmonic 17 (Fig. 3b).

Calculating the high-order harmonic spectrum from SF$_6$ is quite challenging due to the number of associated channels and the presence of resonances in some of them. As a first approximation, we assume that there are no subcycle non-adiabatic multi-electron transitions, which would coherently couple cation states: an electron leaving a given cation electronic state will recombine to the same state. We will not try to simulate the harmonic emission resulting from the coherent superposition of different channels, but simply to estimate the contribution of each channel.

We make an approximation that the recombination dipole moment $d_{\text{rec}}^j$ is equal to the field-free XUV photoionization dipole moment. This was studied experimentally using synchrotron radiation, but the A and B channels are not distinguished$^{29}$. We thus prefer to rely on theoretical calculations, which provide differential cross-sections for all channels. Modelling the photoionization cross-sections of SF$_6$ is difficult and requires taking into account accurately intra- and interchannel couplings$^{30,31}$ (Supplementary Fig. 2). The recent calculations from Jose et al. show good agreement with experiments, once downshifted by 2 eV (Supplementary Fig. 3). The normalized harmonic intensity $\tilde{I}^j$ calculated for each channel using these calculated cross-sections $d_{\text{rec}}^j$ and the calculated $P_{\text{ion}}^j$ from Fig. 2b is shown in Fig. 3b. The picture, which is very complex in the case of XUV photoionization where many channels have similar cross-sections (Supplementary Fig. 3), appears to be simplified by SFI. The tunnelling probability from the ground-state channel X is significantly reduced by the presence of nodes in the corresponding orbital, such that the deeper A channel is favoured.

Below 26 eV, there are broad maxima in the cross-sections from A, B and C channels, which significantly increases their contributions. These maxima are related to the following shape resonances as calculated by Jose et al.: $5t_{1u} \rightarrow \epsilon t_{2g}$ (A channel-40 Mbarn) at 24 eV, $1t_{2u} \rightarrow \epsilon t_{2g}$ (B channel-90 Mbarn) as well at 24 eV and finally the $3e_g \rightarrow \epsilon t_{1u}$ (C channel-15 Mbarn) at 22 eV (ref. 31, see Supplementary Discussion for details). Due to the coupling between X and (A, B) channels, a slight enhancement in the cross-section at the energy of the shape resonance is as well observed in the X channel (37 Mbarn)$^{31}$. Last, another broad peak is predicted in the A channel around 38 eV, associated to the $5t_{1u} \rightarrow \epsilon e_g$ shape resonance (40 Mbarn).

Comparison with the measured harmonic signal (crosses in Fig. 3) indicates that the maximum at H15 may be the signature of the $5t_{1u} \rightarrow \epsilon t_{2g}$ shape resonance in the A channel. By contrast, there is no clear signature around H25 of the resonance predicted
at 38 eV ($5t_{1u} \rightarrow \epsilon e_g$) in the harmonic spectrum. Noteworthy, experimental photoionization cross-sections hardly show any feature related to this resonance either$^{29}$. However, the presence or absence of a maximum in the harmonic spectrum is insufficient to draw any definitive conclusion on the involvement of resonances in the process$^{16}$.

**Revealing resonant HHG using elliptical laser light.** While the qualitative study of harmonic spectra has indicated the possible role of the $5t_{1u} \rightarrow \epsilon t_{2g}$ shape resonance in the emission, a definitive assignment of the structures observed in the experimental spectrum would require a full calculation of the emission, taking into account the complex photorecombination cross-sections and including the effects of laser-induced couplings. How can we identify the presence of resonances without needing a full theoretical model of the harmonic emission? Such identification requires a differential measurement, a measurement that reveals the response of the molecular system to a variation of the laser parameters. This can be achieved by measuring the variation of the harmonic signal with laser ellipticity, $\epsilon$.

High-order harmonic emission in an elliptical laser field can be described in a quasi-static semi-classical picture, as follows. First, electrons tunnel out of the molecule with a finite transverse momentum distribution. This distribution is dictated by the molecular orbital structure, filtered through the tunnelling process$^{32}$. If the electron wavepacket at the exit of the tunnel is spatially localized, then the momentum distribution is broad. During propagation in the continuum, the electron wavepacket spreads transversely. The spread of the returning electron wavepacket depends on the time spent in the continuum (and thus on the harmonic order, the lowest harmonics being emitted before the highest ones$^{33}$), but also on the initial transverse momentum distribution: the broader the initial transverse momentum distribution, the broader the returning wavepacket. With an elliptically polarized driving field, the orthogonal component of the field shifts the electron in the lateral direction, leading to a suppression of the recollision probability. Such suppression depends on both the trajectory length and the initial channel-specific transverse momentum.

The measured harmonic signal as a function of laser ellipticity in argon and in SF$_6$ is shown in Fig. 4a,b. The harmonic intensity decays exponentially with $\epsilon^2$ (ref. 34) and can be fitted by a Gaussian $I_q(\epsilon) = I_0 e^{-\beta_q \epsilon^2}$, where $q$ is the harmonic order and $\beta_q$ the variation with ellipticity. The ellipticity-dependent responses of argon and SF$_6$ are remarkably similar for harmonic 21, whereas that of argon is higher for harmonic 15. The variation of the decay rate $\beta$ with harmonic order is shown in Fig. 4 for Ar and SF$_6$(c). A clear transition is visible about harmonic 17: the harmonic signal decays slower in SF$_6$ in the lower energy range.

The trajectory lengths are similar in the Ar and SF$_6$ measurements (same intensity, close ionization potentials). The similarity of the decay rates with ellipticity in these two species above harmonic 17 shows that the transverse momentum distributions of the involved channels are very similar. Calculations of the momentum distributions indeed give identical results for Ar and the A channel (Supplementary Fig. 4). This measurement thus tends to indicate the predominance of the A channel in the emission above H17. The sudden transition to a regime of lower decay rate with ellipticity below H17 shows a switch in the generating process. This could either be due to a transition to another ionizing state, or more likely to the influence of the resonance in recombination.

The high-order harmonic emission from an elliptical laser field is a vectorial quantity: the harmonic electric field can be decomposed in two projections, parallel and perpendicular to

<figure><img src="image_6.png" /><figcaption>**Figure 4** | HHG in elliptical laser fields. Intensity of harmonic 15 (a) and 21 (b) generated in Ar (grey) and SF$_6$ (black) at $I \approx 1.3 \times 10^{14} \text{ W cm}^{-2}$ as a function of ellipticity. (c) Decay rate $\beta$ of the harmonic signal with ellipticity as a function of harmonic order.</figcaption></figure>

the main axis of the fundamental ellipse. Experimentally, these two projections are resolved by introducing an XUV polarizer between generation and detection. Rotating the polarizer (or equivalently the fundamental ellipse) enables measuring modulations of the harmonic signal, which are characteristic of the Malus law$^{35}$. Assuming that the radiation contains no unpolarized component, we can extract the polarization state of the harmonics as a function of the laser ellipticity, and therefore the evolution of the two projections of the harmonic signal. The results are shown in Fig. 5a,b. The signal of each harmonic is normalized by the sum of the parallel and perpendicular intensities over all ellipticities. The lower decay rate of the harmonic signal with ellipticity is accompanied by the creation of a strong orthogonal component. Figure 5c shows the polarization ellipses of the harmonics obtained with a fundamental ellipticity, $\epsilon_0 = 20\%$. Harmonic 13 to 17 show high degrees of ellipticity and the direction of the polarization ellipse changes sign between H15 and H17. These features were found to be robust against variation of the laser intensity.

We have recently performed a theoretical and experimental study of HHG in Ar atoms, which revealed that resonances could induce strong orthogonal components of the harmonic emission and high degrees of ellipticity$^{36}$. This study also confirmed that the radiation emitted by SF$_6$ in the 20–25 eV range was not unpolarized but highly elliptical. The polarization-resolved measurements presented in Fig. 5 thus constitute a clear demonstration of the contribution of a resonance in the HHG process for harmonic 13 to 15. Comparison with the calculations of Fig. 3b enables us to assign this resonance as the $5t_{1u} \rightarrow \epsilon t_{2g}$ shape resonance in the A channel. We also note that there is no sign of the other shape resonance ($5t_{1u} \rightarrow \epsilon e_g$) predicted around H25 in Fig. 3b.

**Measurement of the recombination phase.** Elliptically driven HHG in SF$_6$ revealed a clear transition in the emission process: harmonics 13 to 17 are produced by resonant emission, probably mostly through the A channel, whereas harmonics above 17 are produced by non-resonant emission. This resonance is expected
<figure><img src="image_7.png" /><figcaption>**Figure 5** | **High-harmonic polarimetry.** Parallel (a) and perpendicular (b) components of the harmonic emission as a function of laser ellipticity. (c) Polarization ellipses of the harmonics generated by a 20% elliptical laser field.</figcaption></figure>

to induce a phase shift in the emission but to the best of our knowledge no data is available on this matter in the literature. Indeed, measuring phase shifts in photoionization matrix elements is a complicated task. Attosecond pulse trains can be used to determine such shifts, the main difficulty being the deciphering of photoelectron spectra containing many components³⁷. By contrast, measuring the phase shift in the harmonic emission is quite straightforward.

The RABBIT technique³⁸ is an attosecond pulse characterization method that measures the relative phase between consecutive harmonics. RABBIT is based on photoionization of a target gas by high-order harmonics, in the presence of a weak fundamental laser field. The photoelectron spectrum shows main lines, resulting from the absorption of harmonics, and sidebands, corresponding to absorption of one harmonic and absorption or emission of an infrared photon (Fig. 6a). Two quantum paths involving two consecutive harmonics lead to the same sideband, which results in an interference process. The sideband intensity thus oscillates as a function of delay between harmonics and the IR. The phase of this oscillation $\Phi_q^{SB}$ encodes the relative phase between the neighbouring harmonics $q+1$ and $q-1$: $\Phi_q^{SB} = \Delta\varphi_q + \Delta\phi_q^{at}$, where $\Delta\varphi_q = \varphi_{q+1} - \varphi_{q-1}$ is the phase difference between harmonics $q+1$ and $q-1$ and $\Delta\phi_q^{at}$ is a scattering phase intrinsic to the target (the atomic phase, that is,

<figure><img src="image_8.png" /><figcaption>**Figure 6** | **Measurement of the harmonic phase.** (a) Principle of the RABBIT method: high harmonics photoionize argon atoms in the presence of an infrared (800 nm) field. Absorption of harmonic 15 and one IR photon or absorption of harmonic 17 and stimulated emission of one IR photon both lead to the production of the same sideband (SB16) in the photoelectron spectrum. (b) Phase of the sideband oscillations as a function of the delay between IR and harmonics. The dotted line is the linear contribution of the group delay dispersion. (c) Recombination phase as a function of harmonic order (light), and recombination phase of the perpendicular (dark) component of the harmonic emission in a 20% elliptical laser field.</figcaption></figure>

the phase accumulated by the photoionized electron in the atomic scattering potential³⁸).

In Fig. 6b, we show the measured harmonic phase difference obtained using SF₆ as a generating gas and argon as a detecting gas. In a simple single-channel non-resonant emission case, the harmonic phase difference is expected to increase linearly with harmonic order, reflecting the group delay dispersion of the attosecond pulse emission. This linear behaviour is visible on the three highest sidebands. A clear deviation is observed at sideband 16, reflecting a specific phase difference between harmonic 15 and 17. A similar behaviour was reported in ref. 39. We decompose the harmonic phase as the sum of a continuum phase, reflecting the group delay and a recombination phase: $\varphi_q = \varphi_q^{continuum} + \varphi_q^{recomb}$. Integrating the phase difference after subtracting the linear contribution of the continuum phase, we obtain the recombination phase plotted in Fig. 6c (for the parallel component). There is a 1.2-rad phase jump through the resonance, with the major part (1 rad) between H15 and H17.

The RABBIT measurement was performed in linear laser polarization and thus measured a phase jump in the parallel recombination dipole moment. This phase jump occurs at the same position as the sign change in the polarization ellipse direction (Fig. 5c), which indicates a phase jump also in the orthogonal dipole moment. Assuming that the phase of the parallel dipole moment does not vary as the laser ellipticity increases from 0 to 20%, one can combine the polarimetry and RABBIT measurements to extract the recombination phase of the perpendicular component (by adding the relative phase between the two vectorial components obtained in polarimetry to the phase of the parallel component obtained in RABBIT). The
results, shown in Fig. 6c, show a larger phase jump (1.7 rad) across the resonance in the perpendicular component. This determination of the recombination phase along the two components of the harmonic emission demonstrates the strength of the combination of RABBIT and polarimetry measurements.

**Response to molecular vibrations.** We have identified a transition between two regimes in HHG from SF₆—resonant (below H17) and non-resonant (above H17)—by manipulating the electron trajectories in the generating mechanism. The emission through a shape resonance involves a peculiar process of electron scattering within the molecular potential, which is expected to be sensitive to changes of the molecular geometry. In the following, we vary the nuclear degrees of freedom via vibrational excitation. As discussed below, this further highlights the specific nature of the resonant HHG process.

Small-molecular structural changes in the generating medium can induce significant modifications in the high-harmonic emission. For instance, Raman pump—strong-field probe experiments in SF₆ showed that the intensity of harmonic 25 to 45 varies by ~10% on molecular vibration⁴⁰,⁴¹, even if the bond length variations remain small (~1%). First, changes in molecular geometry are expected to alter the vertical ionization potentials. Such modification influences the relative amplitudes and phases of the channels. For example, a 65-meV change in Iₚ results in a 100-mrad phase change for a plateau harmonic. Second, molecular vibrations can lift the electronic degeneracy of intermediate cation states, increasing the number of distinct channels. Finally, any distortions of the electronic wavefunctions associated with these vibrations may affect both the ionization and recombination dipole moments, in particular in the vicinity of resonances.

We have performed pump-probe experiments using an intense (~7 × 10¹³ W cm⁻²) 800 nm pulse to Raman-pump vibrational modes in the SF₆ molecules and produce high harmonics in the excited medium. In addition to conventional experiments in which the harmonic intensity is measured (Fig. 7a,c), we have used two-source interferometry to determine the evolution of the harmonic phase with pump-probe delay (Fig. 7b,d)⁵,⁴² (see Methods). The harmonic emission is periodically modulated as a function of pump-probe delay. The typical modulation depth of the harmonic signals are 5% for the intensity and 100 mrad for the harmonic phase. Three of the six normal modes of SF₆ are Raman active vibrations: v₁—fully symmetric and strongly active mode A₁g with a quantum of 774 cm⁻¹ (vibrational period of ≈43 fs), v₂—the doubly degenerate mode E_g with a quantum of 643 cm⁻¹ (≈52 fs) and finally the triply degenerate v₅ – T₂g mode with a quantum of 524 cm⁻¹ (≈63 fs)⁴³. A Fourier analysis is performed to determine which vibrational modes modulate the signal (Fig. 7a,b).

The weight of each vibrational mode in the intensity and phase modulations is extracted by integrating the corresponding peak and normalizing to the total signal (Fig. 7c,d). The results show a clear transition around harmonic 17, consistent with the switching from resonant to non-resonant emission measured in static molecules. Below harmonic 17, in the region of the shape resonance, the harmonic intensity is exclusively modulated by the v₁ mode. Clearly the asymmetric distorting mode v₅ does not impact the amplitude of the harmonic emission through the shape resonance. Above harmonic 17, the v₅ mode becomes slightly predominant. There is no significant contribution of the v₂ mode except for harmonic 13. Interestingly, the behaviour of the harmonic phase is very different but also shows a transition. The phase of the resonant emission varies similarly for v₁ and v₅ nuclear dynamics and shows a weak contribution from v₂. The phase of the emission above H17 seems insensitive to v₂ and

dominated by v₁. These results represent to our knowledge the first dynamical study of a shape resonance in HHG. Shape resonances are known to depend on bond length⁴⁴ and our results reflect this dependence as the molecule vibrates.

## Discussion

Multi-channel HHG in polyatomic molecules is a rich and complex process that initiates a broad range of multi-electron phenomena. The identification and isolation of the different channels is a major challenge, which requires the development and integration of new experimental tools. In this work we introduced and combined several experimental schemes in an advanced spectroscopic study, enabling us to elucidate clear features related to resonant and non-resonant HHG in SF₆. First, we identified the appearance of different SFI channels by performing wavelength-dependent fragment-resolved ATI measurements. Varying the wavelength addressed the issue of overlapping ATI combs that can occur at specific wavelengths, notably at the 800-nm driving frequency implemented in the majority of this work. We introduced a new scheme to HHG spectroscopy by resolving the spectral dependence of the normalized ellipticity response. In non-resonant cases, this measure enables isolation of the transverse shape of the electronic wavefunction at the moment of tunnel ionization, removing the complexity imposed by the additional steps of the interaction. In resonant cases, the laser ellipticity induces the creation of a strong orthogonal dipole moment and a large ellipticity of the harmonic emission, which were detected by polarimetry. We made use of the coherent properties of the process via RABBIT measurements, which enabled the measurement of the recombination phase through the resonance. Finally, we have made use of the nuclear degrees of freedom—via vibrational dynamics—as an important control knob that manipulates both the phase and amplitude of the emission, advancing our ability to identify the specific nature of the resonant generation process.

All the experimental observations lead to the assignment of the 5t₁u → ct₂g shape resonance in the A channel as the main contributor to the generation of harmonics 13 to 17. This result is particularly striking compared with photoionization experiments in which it is the resonance in the B channel that dominates the process. In HHG, this channel is suppressed by the inefficiency of strong-field ionization due to a large number of nodal planes in the orbital profile. No feature related to the 5t₁u → ce_g shape resonance in the A channel was observed in any of the measurements. This absence, which is consistent with the very weak contribution of this resonance to experimental photoionization cross-sections, calls for additional theoretical studies. In particular, Jose et al.³¹ suggested that vibrational excitation could drastically affect the shape resonance. We believe that high-harmonic spectroscopy is the ideal tool to study such effects. On one hand, the vibration of the neutral can be controlled using stimulated impulsive Raman scattering as we have shown. On the other hand, the dynamics of the ion between ionization and recombination are encoded in the harmonic spectrum⁴⁵, and can be probed by measuring the wavelength dependence of the harmonic signal. Last but no least, the coherence of the HHG process enables performing phase-resolved measurements and accessing the temporal properties of the emission, observing the decay of autoionizing states in the time domain.

The identification of the basic mechanism associated with multiple channels is the first fundamental step in characterizing time-resolved multi-electron processes, the long term goal being to resolve charge migrations or field-induced coupling between the channels. We believe that it will be possible to isolate such dynamics by increasing the dimensionality of the experiment by
<figure><img src="image_9.png" /></figure>

**Figure 7** | The modulation of HHG by specific vibrations in ground-state SF₆ molecules. Fourier transform of the intensity (a) and phase (b) modulations with pump-probe delay for harmonic 13 (dots) and 19 (line) showing the contribution of different Raman modes: $v_5 - T_{2g}$ (524 cm⁻¹, yellow), $v_2 - E_g$ (643 cm⁻¹, green) and $v_1 - A_{1g}$ (774 cm⁻¹, blue). The contribution of each mode, obtained by integration of the peak in the FFT, is shown as a function of harmonic order in (c) intensity and (d) phase.

scanning the fundamental wavelength of the strong laser field in HHG. Such scanning will scale the interferometer, presented in Fig. 1, enabling the decoupling of the different degrees of freedom involved in the process. In summary, we believe that the combination of methods presented here constitute an important building block for strong-field spectroscopy and will lead to the discoveries of a broad range of multi-electron strong-field phenomena in complex targets such as polyatomic molecules.

**Methods**

**Fragment-resolved ATI.** Experiments were performed using a Legend Elite Duo laser system (Coherent Inc.) that delivered 35 fs, 3.1 mJ, 800 nm pulses at 1 kHz (Molecular Photonics Laboratory, NRC Ottawa). Laser fundamental pulses were implemented directly in the 800-nm experiments, whereas ~40-fs 400-nm pulses were generated through second harmonic generation in thin β-BaBO₄ crystals (150-μm thick). Approximately, 40 fs 480, 532 and 590 nm pulses were generated through sum frequency generation of the signal pulses from an optical parametric amplifier (Light Conversion, TOPAS-C) and the fundamental beam in thin β-BaBO₄ crystals (250 μm thick). Laser pulses were focused with a $f=50$ cm spherical mirror onto a continuous, neat expansion of SF₆ molecules (expanded through a 50-μm pinhole ~60 cm from the laser-molecular beam interaction region in a differentially pumped chamber). The molecular beam PEPICO (PhotoElectron PhotoIon COincidence) spectrometer used to collect the resulting photoelectrons and photoioins has been described elsewhere²³. In brief, the kinetic energies of the photoelectrons were resolved by time of flight using a wide-bore permanent magnet, magnetic bottle spectrometer. Photoioins were collected using a coaxial pulsed-field Wiley-McLaren time-of-flight mass spectrometer. The photoelectron and photoion time-of-flight signals were post processed using a covariance analysis⁴⁶,⁴⁷ to correlate specific electron signals with the dominant SF₆ ions. The typical exponential ATI decay in energy has been subtracted for clarity. The energy and intensity calibration is a multi-step process, which uses the ATI progressions themselves as a ruler, and also involves investigating the intensity-dependent ATI spectra of O₂ and H₂O impurities, where Stark shifts can be neglected.

**Calculations of SFI rates.** Bound-state neutral and cation wavefunctions were calculated using aug-cc-pVTZ basis set⁴⁸, large-core effective core potentials⁴⁹ and RASSCF wavefunctions with 35 (cation) or 36 (neutral) electrons in 25 orbitals, allowing single and double excitations. For the cation, the active space was optimized for an equally weighted average of the six lowest states of the cation (X through D). The RASSCF vertical IPs were shifted uniformly by -0.512 eV. At equilibrium, $O_h$ geometry (R(S-F) = 1.575 Å), the resulting cation state energies are within 0.3 eV of the correlated MCQDPT2(27/28,21)/aug-cc-pVTZ values.

The magnitude of the vibrational distortions induced by an impulsive Raman pump was determined⁵⁰ from normal vibrational modes and dipole polarizabilities calculated at the CCSD/aug-cc-pVTZ level⁵¹. Calculated displacements assumed pump energy flux of ≈8.4 J cm⁻². Pump and probe polarizations were assumed to be collinear and were integrated on an order-9 Lebedev grid⁵². RASSCF and MCQDPT2 (ref. 53 calculations were performed with GAMESS-US⁵⁴. CCSD calculations used CFOUR⁵⁵. Strong-field ionization calculations²⁵,²⁶ used Cartesian product grid with 0.18 Bohr spacing, extending to ±18 Bohr from the origin, with an absorbing boundary⁵⁶ of 9.4 Bohr at the box edges. The simulation used 14 spatial channels, corresponding to the X, A, B, C and D final states of the cation. Both laser and correlation coupling terms were included. Leap-frog time evolution used a time step of 0.003 AU.

**High-order harmonic spectroscopy.** The high-harmonic spectroscopy experiments were performed using the Aurore laser system at CELIA, which delivers 28 fs, 7 mJ, 800 nm pulses at 1 kHz. High harmonics were generated by focusing the pulses with a $f=50$ cm lens in a continuous jet produced by a 250-μm nozzle backed by a few hundred mbars of gas. The emitted radiation was analysed by an extreme-ultraviolet spectrometer (1,200 mm⁻¹ grating, dual Micro Channel Plates and a CCD camera).

To manipulate the laser ellipticity, two broadband zero-order waveplates were inserted in the laser beam: a motorized half waveplate in front of a fixed quarter waveplate. The harmonic signal was recorded as a function of the angle $\alpha$ between the half and quarter waveplate axis that determines the laser ellipticity $\epsilon = \tan(2\alpha)$. The harmonic polarization state was characterized using an XUV polarizer consisting of three bare gold mirrors under 75–60–75° incidence, which leads to a typical contrast of 10–20 between S and P polarizations in the studied spectral
range. A second half-wave plate was used to rotate the laser polarization, which is equivalent to rotating the XUV polarizer, to record the Malus' law.

For RABBIT measurements, we used the experimental set-up described in ref. 57. In brief, the harmonic and infrared beams were delayed using a beam-splitting dichroic delay-line. They were focused by a toroidal mirror in a magnetic bottle time-of-flight spectrometer that recorded the electron spectrum produced by photoionizing argon atoms, as a function of delay between XUV and infrared beams.

Phase-resolved pump-probe experiments were performed using a 0–π phase mask⁴² inserted on the harmonic generating beam, before focusing³,⁵,⁸,⁹. This mask created an interference pattern at focus, which resulted in the creation of two harmonic sources spatially separated by ~100 µm. In the far field, these sources produced an interference pattern, similar to Young's slits. An 800-nm pump beam was superimposed on one of the two sources, exciting the molecules by impulsive stimulated Raman scattering. The evolution of the phase of the harmonic emission as a function of pump-probe delay was measured by determining the shift of the interference fringes, using a spatial Fourier transform.

## References

1. Keldysh, L. Ionization in the field of a strong electromagnetic wave. Sov. JETP 20, 1307 (1965).

2. Muth-Böhm, J., Becker, A. & Faisal, F. H. M. Suppressed molecular ionization for a class of diatomics in intense femtosecond laser fields. Phys. Rev. Lett. 85, 2280–2283 (2000).

3. Tong, X. M., Zhao, Z. X. & Lin, C. D. Theory of molecular tunneling ionization. Phys. Rev. A 66, 033402 (2002).

4. Lezius, M. et al. Nonadiabatic multielectron dynamics in strong field molecular ionization. Phys. Rev. Lett. 86, 51–54 (2001).

5. Smirnova, O., Patchkovskii, S., Mairesse, Y., Dudovich, N. & Ivanov, M. Y. Strong-field control and spectroscopy of attosecond electron-hole dynamics in molecules. Proc. Natl Acad. Sci. USA 106, 16556–16561 (2009).

6. Haessler, S. et al. Attosecond imaging of molecular electronic wavepackets. Nat. Phys. 6, 200–206 (2010).

7. Agostini, P., Fabre, F., Mainfray, G., Petite, G. & Rahman, N. K. Free-free transitions following six-photon ionization of xenon atoms. Phys. Rev. Lett. 42, 1127–1130 (1979).

8. McPherson, A. et al. Studies of multiphoton production of vacuum-ultraviolet radiation in the rare gases. J. Opt. Soc. Am. B 4, 595 (1987).

9. Corkum, P. B. Plasma perspective on strong field multiphoton ionization. Phys. Rev. Lett. 71, 1994–1997 (1993).

10. Krause, J. L., Schafer, K. J. & Kulander, K. C. High-order harmonic generation from atoms and ions in the high intensity regime. Phys. Rev. Lett. 68, 3535–3538 (1992).

11. Lewenstein, M., Balcou, P., Ivanov, M. Y., L'Huillier, A. & Corkum, P. B. Theory of high-harmonic generation by low-frequency laser fields. Phys. Rev. A 49, 2117–2132 (1994).

12. Le, A.-T., Lucchese, R., Tonzani, S., Morishita, T. & Lin, C. Quantitative rescattering theory for high-order harmonic generation from molecules. Phys. Rev. A 80, 013401 (2009).

13. Strelkov, V. Role of autoionizing state in resonant high-order harmonic generation and attosecond pulse production. Phys. Rev. Lett. 104, 123901 (2010).

14. Tudorovskaya, M. & Lein, M. High-order harmonic generation in the presence of a resonance. Phys. Rev. A 84, 013430 (2011).

15. Shiner, A. D. et al. Probing collective multi-electron dynamics in xenon with high-harmonic spectroscopy. Nat. Phys. 7, 464467 (2011).

16. Vozzi, C. et al. Phase-matching effects in the generation of high-energy photons by mid-infrared few-cycle laser pulses. New J. Phys. 13, 073003 (2011).

17. Haessler, S. et al. Phase distortions of attosecond pulses produced by resonance-enhanced high harmonic generation. New J. Phys. 15, 013051 (2013).

18. McFarland, B., Farrell, J., Bucksbaum, P. & Gühr, M. High-order harmonic phase in molecular nitrogen. Phys. Rev. A 80, 033412 (2009).

19. Mairesse, Y. et al. High harmonic spectroscopy of multichannel dynamics in strong-field ionization. Phys. Rev. Lett. 104, 213601 (2010).

20. Shafir, D. et al. Resolving the time when an electron exits a tunneling barrier. Nature 485, 343–346 (2012).

21. Jin *et al.* Intensity dependence of multiple orbital contributions and shape resonance in high-order harmonic generation of aligned N₂ molecules. Phys. Rev. A 85, 013405 (2012).

22. Lee, K. F., Villeneuve, D. M., Corkum, P. B., Stolow, A. & Underwood, J. G. Field-free three-dimensional alignment of polyatomic molecules. Phys. Rev. Lett. 97, 173001 (2006).

23. Boguslavskiy, A. E. et al. The multielectron ionization dynamics underlying attosecond strong-field spectroscopies. Science 335, 1336–1340 (2012).

24. Mikosch, J. *et al.* Channel- and angle-resolved above threshold ionization in the molecular frame. Phys. Rev. Lett. 110, 023004 (2013).

25. Spanner, M. & Patchkovskii, S. One-electron ionization of multielectron systems in strong nonresonant laser fields. Phys. Rev. A 80, 063411 (2009).

26. Spanner, M. & Patchkovskii, S. Molecular strong field ionization and high harmonic generation: a selection of computational illustrations. Chem. Phys. 414, 10–19 (2013).

27. Lyngå, C., L'Huillier, A. & Wahlström, C.-G. High-order harmonic generation in molecular gases. J. Phys. B 29, 3293–3302 (1996).

28. Samson, J. A. R. & Stolte, W. C. Precision measurements of the total photoionization cross-sections of He, Ne, Ar, Kr, and Xe. J. Electron Spectrosc. Relat. Phenom. 123, 265 (2002).

29. Holland, D. M. P. et al. An experimental and theoretical study of the valence shell photoelectron spectrum of sulphur hexafluoride. Chem. Phys. 192, 333 (1995).

30. Stener, M., Toffoli, D., Fronzoni, G. & Decleva, P. Time dependent density functional study of the photoionization dynamics of SF6. J. Chem. Phys. 124, 114306–114306-13 (2006).

31. Jose, J., Lucchese, R. R. & Rescigno, T. N. Interchannel coupling effects in the valence photoionization of SF6. J. Chem. Phys. 140, 204305 (2014).

32. Murray, R., Spanner, M., Patchkovskii, S. & Ivanov, M. Tunnel ionization of molecules and orbital imaging. Phys. Rev. Lett. 106, 173001 (2011).

33. Mairesse, Y. et al. Attosecond synchronization of high-harmonic soft x-rays. Science 302, 1540 (2003).

34. Ivanov, M. Y., Brabec, T. & Burnett, N. Coulomb corrections and polarization effects in high-intensity high-harmonic emission. Phys. Rev. A 54, 742 (1996).

35. Antoine, P., Carré, B., L'Huillier, A. & Lewenstein, M. Polarization of high-order harmonics. Phys. Rev. A 55, 13141324 (1997).

36. Ferré, A. et al. A table-top ultrashort light source in the extreme-ultraviolet for time-resolved circular dichroism experiments. Nat. Phot., doi:10.1038/NPHOTON.2014.314 (2014).

37. Haessler, S. et al. Phase-resolved attosecond near-threshold photoionization of molecular nitrogen. Phys. Rev. A 80, 011404(R) (2009).

38. Paul, P. M. et al. Observation of a train of attosecond pulses from high harmonic generation. Science 292, 1689 (2001).

39. Rothhardt, J. et al. in Laser Science. OSA Technical Digest (online), LW1H.6 (Optical Society of America, 2012).

40. Wagner, N. L. et al. Monitoring molecular dynamics using coherent electrons from high harmonic generation. Proc. Natl Acad. Sci. USA 103, 13279–13285 (2006).

41. Walters, Z. B., Tonzani, S. & Greene, C. H. Vibrational interference of raman and high harmonic generation pathways. Chem. Phys. 366, 103–114 (2009).

42. Camper, A. et al. High-harmonic phase spectroscopy using a binary diffractive optical element. Phys. Rev. A 89, 043943 (2014).

43. Boudon, V., Domènech, J. L., Bermejo, D. & Willner, H. High-resolution raman spectroscopy of the v1 region and Raman-Raman double resonance spectroscopy of the 2v1-v1 band of 32SF6 and 34SF6. determination of the equilibrium bond length of sulfur hexafluoride. J. Mol. Spectrosc. 228, 392–400 (2004).

44. Piancastelli, M. The neverending story of shape resonances. J. Electron Spectrosc. Relat. Phenom. 100, 167190 (1999).

45. Baker, S. et al. Probing proton dynamics in molecules on an attosecond time scale. Science 312, 424427 (2006).

46. Frasinski, L. et al. Molecular H₂ in intense laser fields probed by electron-electron, electron-ion, and ion-ion covariance techniques. Phys. Rev. A 46, R6789 (1992).

47. Mikosch, J. & Patchkovskii, S. Coincidence and covariance data acquisition in photoelectron and -ion spectroscopy. i. formal theory. J. Mod. Optic. 60, 1426–1438 (2013).

48. Woon, D. E. & Dunning, T. H. Gaussian basis sets for use in correlated molecular calculations. III. The atoms aluminum through argon. J. Chem. Phys. 98, 1358 (1993).

49. Bergner, A. et al. Ab initio energy-adjusted pseudopotentials for elements of groups 1317. Mol. Phys. 6, 1431 (1993).

50. De Silvestri, S. et al. Femtosecond time-resolved measurements of optic phonon dephasing by impulsive stimulated raman scattering in ¹Î -perylene crystal from 20 to 300 K. Chem. Phys. Lett. 116, 145 (1985).

51. Kallay, M. & Gauss, J. Calculation of frequency-dependent polarizabilities using general coupled-cluster models. J. Mol. Struc. Theochem. 768, 71 (2006).

52. Lebedev, V. I. Values of the nodes and weights of ninth to seventeenth order gauss-markov quadrature formulae invariant under the octahedron group with inversion. USSR Comp. Math. Math. Phys. 15, 48 (1975).

53. Nakano, H. Quasidegenerate perturbation theory with multiconfigurational self-consistent-field reference functions. J. Chem. Phys. 99, 7983 (1993).

54. Schmidt, M. W. et al. General atomic and molecular electronic structure system. J. Comput. Chem. 14, 1347 (1993).

55. Stanton, J. F. et al. CFOUR, coupled-cluster techniques for computational chemistry, a quantum-chemical program package. Available at http://www.cfour.de.

56. Manolopoulos, D. E. Derivation and reflection properties of a transmission-free absorbing potential. J. Chem. Phys. 117, 9552 (2002).
57. Loch, R. A. et al. Phase characterization of the reflection on an extreme UV multilayer: comparison between attosecond metrology and standing wave measurements. *Opt. Lett.* **36**, 3386–3388 (2011).

58. Kanai, T., Takahashi, E. J., Nabekawa, Y. & Midorikawa, K. Destructive interference during high harmonic generation in mixed gases. *Phys. Rev. Lett.* **98**, 153904 (2007).

59. Yang, L. et al. Energy-dependent valence photoelectron spectra of SF6. ab initio calculations and measurements. *J. Electron Spectrosc. Relat. Phenom.* **94**, 163–179 (1998).

## Acknowledgements

We thank R. Lausten (NRC), R. Bouillaud, C. Medina and L. Merzeau for technical assistance and S. Haessler for discussions. S. Patchkovskii acknowledges helpful discussions with and suggestions by M. Spanner (NRC). A.S. acknowledges NSERC (Canada) for financial support via the Discovery Grant program. N.D. acknowledges support from the Israeli Centers of Research Excellence program, the Crown photonics center, the Israeli Science Foundation, the Minerva Foundation and the European Research Council starting investigator grant MIDAS. We acknowledge financial support of the Conseil Regional d'Aquitaine (20091304003 ATTOMOL and COLA 2 No. 2.1.3-09010502), l'Agence Nationale pour la Recherche (ANR-14-CE32-0014 MISFITS), the European Union (Laserlab-Europe II NÅ°228334), l'Ambassade de France au Canada and NSERC (Canada).

## Author contributions

A.E.B., J.M. and I.W. carried out the strong-field ionization measurements, S. Patchkovskii the strong-field ionization calculations, V.B., A.F. and Y.M. the

high-harmonic polarimetry measurements, N.F., J.G., G.G. andf Y.M. the RABBIT measurements, V.B., B.D.B, A.C., M.D., A.F., Y.M., T.R., H.S. and D.S the molecular vibration measurements. F.B., D.D. and S. Petit operated the laser system at CELIA. B.F., A.F. and Y.M. analysed the high-harmonic data. A.E.B., J.M., S. Patchkovskii, A.S. and I.W. analysed the strong-field ionization data. N.D. and Y.M. designed the manuscript. All authors contributed to the interpretation of the data and writing of the manuscript.

## Additional information

**Supplementary Information** accompanies this paper at http://www.nature.com/naturecommunications

**Competing financial interests:** The authors declare no competing financial interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**How to cite this article:** Ferré, A. et al. Multi-channel electronic and vibrational dynamics in polyatomic resonant high-order harmonic generation. *Nat. Commun.* 6:5952 doi: 10.1038/ncomms6952 (2015).

<figure><img src="image_10.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/