Dynamic two-center interference in high-order harmonic generation from molecules with attosecond nuclear motion

Ultrafast Picometer-Resolved Molecular Structure Imaging by Laser-Induced High-Order Harmonics

Real-time visualization of molecular transformations is a captivating yet challenging frontier of ultrafast optical science and physical chemistry. While ultrafast x-ray and electron diffraction methods can achieve the needed subangstrom spatial resolution, their temporal resolution is still limited to hundreds of femtoseconds, much longer than the few femtoseconds required to probe real-time molecular dynamics. Here, we show that high-order harmonics generated by intense femtosecond lasers can be used to image molecules with few-ten-attosecond temporal resolution and few-picometer spatial resolution. This is achieved by exploiting the sensitive dependence of molecular recombination dipole moment to the geometry of the molecule at the time of harmonic emission. In a proof-of-principle experiment, we have applied this high-harmonic structure imaging (HHSI) method to monitor the structural rearrangement in NH_{3}, ND_{3}, and N_{2} from one to a few femtoseconds after the molecule is ionized by an intense laser. Our findings establish HHSI as an effective approach to resolve molecular dynamics with unprecedented spatiotemporal resolution, which can be extended to trace photochemical reactions in the future.

Coulomb-induced emission time shifts in high-order harmonic generation from H2

Accurate emission times of high-order harmonic generation (HHG) are vital for high-precision ultrafast detection in attosecond science, but a quantitative analysis of Coulomb effects on this time is absent in the molecular HHG. Here, we investigate the Coulomb-induced emission-time shift in HHG of H2+ with two different internuclear distances R, where the times obtained via the Gabor transform of numerical data from solving the time-dependent Schrödinger equation are used as simulation experiment results. Based on the molecular strong-field approximation, we develop a trajectory-resolved classical model that takes into account the molecular two-center structure. By selecting appropriate electron trajectories and including Coulomb interactions, the classical trajectory method can reproduce Gabor emission times well. This consistence reveals that Coulomb tails cause an emission-time shift of ∼35 as at the R = 2.0 a.u. case and of ∼40-60 as at the R = 2.6 a.u. case under the present laser parameters when compared to the Coulomb-free quantum-orbit model. Our results are of significance to probe the attosecond dynamics via two-center interference.

Vibronic Correlations in Molecular Strong-Field Dynamics

We investigate the ultrafast vibronic dynamics triggered by intense femtosecond infrared pulses in small molecules. Our study is based on numerical simulations performed with 2D model molecules and analyzed in the perspective of the renowned Lochfrass and bond-softening models. We give a new interpretation of the observed nuclear wave packet dynamics with a focus on the phase of the bond oscillations. Our simulations also reveal intricate features in the field-induced nuclear motion that are not accounted for by existing models. Our analyses assign these features to strong dynamical correlations between the active electron and the nuclei, which significantly depend on the carrier envelope phase of the pulse, even for relatively "long" pulses, which should make them experimentally observable.

Real-Time Extension of TAO-DFT

Thermally assisted occupation density functional theory (TAO-DFT) has been an efficient electronic structure method for studying the ground-state properties of large electronic systems with multi-reference character over the past few years. To explore the time-dependent (TD) properties of electronic systems (e.g., subject to an intense laser pulse), in this work, we propose a real-time (RT) extension of TAO-DFT, denoted as RT-TAO-DFT. Moreover, we employ RT-TAO-DFT to study the high-order harmonic generation (HHG) spectra and related TD properties of molecular hydrogen H2 at the equilibrium and stretched geometries, aligned along the polarization of an intense linearly polarized laser pulse. The TD properties obtained with RT-TAO-DFT are compared with those obtained with the widely used time-dependent Kohn-Sham (TDKS) method. In addition, issues related to the possible spin-symmetry breaking effects in the TD properties are discussed.

Minimum structure of high-order harmonic spectrum from molecular multi-orbital effects involving inner-shell orbitals

The spectral features of high-order harmonic spectra can provide rich information for probing the structure and dynamics of molecules in intense laser fields. We theoretically study the high harmonic spectrum with the laser polarization direction perpendicular to the N2O molecule and find a minimum structure in the plateau region of the harmonic spectrum. Through analyzing the time-dependent survival probability of different electronic orbitals and the time-dependent wave packet evolution, it is found that this minimum position is caused by the harmonic interference of HOMO a, HOMO-1, and HOMO-3 a orbitals. Moreover, this interference minimum is discovered over a wide frequency range of 0.087 a.u. to 0.093 a.u., as well as a range of driving laser intensities with peak amplitudes between 0.056 a.u. and 0.059 a.u.. This study sheds light on the multi-electron effects and ultrafast dynamics of inner-shell electrons in intense laser pulses, which are crucial for understanding and controlling chemical reactions in molecules.

Two-Center Interference in the Photoionization Delays of Kr_{2}

We present the experimental observation of two-center interference in the ionization time delays of Kr_{2}. Using attosecond electron-ion-coincidence spectroscopy, we simultaneously measure the photoionization delays of krypton monomer and dimer. The relative time delay is found to oscillate as a function of the electron kinetic energy, an effect that is traced back to constructive and destructive interference of the photoelectron wave packets that are emitted or scattered from the two atomic centers. Our interpretation of the experimental results is supported by solving the time-independent Schrödinger equation of a 1D double-well potential, as well as coupled-channel multiconfigurational quantum-scattering calculations of Kr_{2}. This work opens the door to the study of a broad class of quantum-interference effects in photoionization delays and demonstrates the potential of attosecond coincidence spectroscopy for studying weakly bound systems.

Channel Coupling Dynamics of Deep-Lying Orbitals in Molecular High-Harmonic Generation

Investigation on structures in the high-harmonic spectrum has provided profuse information of molecular structure and dynamics in intense laser fields, based on which techniques of molecular ultrafast dynamics imaging have been developed. Combining ab initio calculations and experimental measurements on the high-harmonic spectrum of the CO_{2} molecule, we find a novel dip structure in the low-energy region of the harmonic spectrum which is identified as fingerprints of participation of deeper-lying molecular orbitals in the process and decodes the underlying attosecond multichannel coupling dynamics. Our work sheds new light on the ultrafast dynamics of molecules in intense laser fields.

Optical solitons in media with focusing and defocusing saturable nonlinearity and a parity-time-symmetric external potential

We report results for solitons in models of waveguides with focusing or defocusing saturable nonlinearity and a parity-time ([Formula: see text])-symmetric complex-valued external potential of the Scarf-II type. The model applies to the nonlinear wave propagation in graded-index optical waveguides with balanced gain and loss. We find both fundamental and multipole solitons for both focusing and defocusing signs of the saturable nonlinearity in such [Formula: see text]-symmetric waveguides. The dependence of the propagation constant on the soliton's power is presented for different strengths of the nonlinearity saturation, S The stability of fundamental, dipole, tripole and quadrupole solitons is investigated by means of the linear-stability analysis and direct numerical simulations of the corresponding (1+1)-dimensional nonlinear Schrödinger-type equation. The results show that the instability of the stationary solutions can be mitigated or completely suppressed, increasing the value of SThis article is part of the theme issue 'Dissipative structures in matter out of equilibrium: from chemistry, photonics and biology (part 1)'.

Monitoring ultrafast vibrational dynamics of isotopic molecules with frequency modulation of high-order harmonics

Molecules constituted by different isotopes are different in vibrational modes, making it possible to elucidate the mechanism of a chemical reaction via the kinetic isotope effect. However, the real-time observation of the vibrational motion of isotopic nuclei in molecules is still challenging due to its ultrashort time scale. Here we demonstrate a method to monitor the nuclear vibration of isotopic molecules with the frequency modulation of high-order harmonic generation (HHG) during the laser-molecule interaction. In the proof-of-principle experiment, we report a red shift in HHG from H₂ and D₂. The red shift is ascribed to dominant HHG from the stretched isotopic molecules at the trailing edge of the laser pulse. By utilizing the observed frequency shift, the laser-driven nuclear vibrations of H₂ and D₂ are retrieved. These findings pave an accessible route toward monitoring the ultrafast nuclear dynamics and even tracing a chemical reaction in real time.

Previous studies on high harmonic generation from molecules have been used to identify the spectral properties and orbital contributions. Here the authors measure the isotopic effects in the energy shift of the HHG spectra caused by the nuclear motion of the molecules.

An apparatus for quantitative high-harmonic generation spectroscopy in molecular vapours

We present an apparatus for performing gas phase high-harmonic generation spectroscopy of molecules primarily found in the liquid phase. Liquid molecular samples are heated in a temperature controlled bath and their vapour is used to back a continuous flow gas jet, with vapour pressures of over 1 bar possible. In order to demonstrate the system, we perform high harmonic spectroscopy experiments in benzene with a 1.8 μm driving field. Using the unique capabilities of the system, we obtain spectra that are nearly free from the effects of longitudinal phase-matching, amenable to comparison with advanced numerical modelling.

Resonance enhanced high-order harmonic generation in H2+ by two sequential laser pulses

We investigate high-order harmonic generation in H2+ by using two sequential laser pulses, which consist of a 800-nm pump pulse and a time-delayed 1600-nm probe pulse. Based on the solution of the time-dependent Schrödinger equation, we demonstrate that the harmonic cutoff in our two-pulse scheme is significantly extended compared to that in the 1600-nm probe pulse alone. Meanwhile, the harmonic efficiency is enhanced by 2-3 orders of magnitude due to charge-resonance-enhanced ionization steered by the 800-nm pump pulse. By using a probe pulse with longer wavelength, our scheme can be used for efficient high harmonic generation in the water window region. In addition, the influence of the intensity of the pump pulse and the relative time delay of the two laser pulses on the harmonic generation are also investigated.

Attosecond Probing of Nuclear Dynamics with Trajectory-Resolved High-Harmonic Spectroscopy

We report attosecond-scale probing of the laser-induced dynamics in molecules. We apply the method of high-harmonic spectroscopy, where laser-driven recolliding electrons on various trajectories record the motion of their parent ion. Based on the transient phase-matching mechanism of high-order harmonic generation, short and long trajectories contributing to the same harmonic order are distinguishable in both the spatial and frequency domains, giving rise to a one-to-one map between time and photon energy for each trajectory. The short and long trajectories in H_{2} and D_{2} are used simultaneously to retrieve the nuclear dynamics on the attosecond and ångström scale. Compared to using only short trajectories, this extends the temporal range of the measurement to one optical cycle. The experiment is also applied to methane and ammonia molecules.

Role of tunnel ionization in high harmonic generation from substituted benzenes

We theoretically study high-harmonic generation in toluene, ortho-xylene and fluorobenzene driven by a 1.8 μm ultrashort pulse. We find that the chemical substitutions have a strong influence on the amplitude and phase of the emission from the highest occupied molecular orbital, despite having a small influence on the orbital itself. We show that this influence is due to the tunnel ionization step, which depends critically on the sign and amplitude of the asymptotic part of the wave function. We discuss how these effects would manifest in phase-sensitive high-harmonic generation spectroscopy experiments.

Assessment of the LFAs-PBE exchange–correlation potential for high-order harmonic generation of aligned H2+ molecules

We examine the performance of the LFAs-PBE exchange-correlation potential for the HHG spectra and related properties of H₂⁺ molecules aligned parallel and perpendicular to the polarization of an intense linearly polarized laser pulse.

Carrier-envelope phase dependence of molecular harmonic spectral minima induced by mid-infrared laser pulses

The spectral minima in harmonic spectra of H2+ induced by mid-infrared laser pulses are numerically investigated based on two models of Born-Oppenheimer (BO) and non-Born-Oppenheimer (NBO) approximations. The simulation results show that, with the variation of the mid-infrared laser's carrier-envelope phase (CEP), the spectral minima positions (SMPs) are fixed for the BO model, while oscillate periodically for the NBO model. This can be understood by the two-center-destructive-interference theory via the detailed investigation to several physical quantities for each CEP case, such as SMPs, effective potential, internuclear separation and the electron's de Broglie wavelength at the time for interference occurring. The fittings to these quantities' CEP-dependent curves demonstrate that they follow a variation law in the form of a sine function.

Theoretical study of the inversion motion of the ammonia cation with subfemtosecond resolution for high-harmonic spectroscopy

In a recent PACER (Probing Attosecond dynamics with Chirp-Encoded Recollisions) experiment on ammonia that comprises a comparison of the high-harmonic spectra of the isotopes NH3 and ND3, the nuclear dynamics of the created ammonia cation is traced with a time resolution of about 100 attoseconds. For modelling the experiment the autocorrelation functions between the neutral initial state and the ionic wave packet are extracted from experimental photoelectron spectra incorporating a correction for the geometry-dependent strong-field ionisation probability. Good agreement is found between model and experiment, but in addition an unexpected maximum in the autocorrelation ratio is predicted by the model, however occurring at 5 fs and thus outside the experimentally covered time interval. In this work the autocorrelation functions are calculated explicitly using a one-dimensional model for describing the inversion motion of ammonia and its cation, adopting a position-dependent mass for considering the coupling with the stretching mode of the hydrogen atoms in neutral ammonia. This results in a clear physical picture explaining the occurrence of the previously predicted maximum in the ratio of the autocorrelation functions. Furthermore, different initial states and two different ways of incorporating strong-field corrections to the Franck-Condon approximation are briefly discussed.

Attosecond nuclear dynamics in the ammonia cation: relation between high-harmonic and photoelectron spectroscopies

We report measurements of the umbrella motion in the ammonia cation on the attosecond time scale. The motion is prepared by strong-field ionization and probed by photorecombination through the process of high-harmonic generation. Performing such measurements at multiple wavelengths (0.8, 1.44, 1.8 μm) enables us to follow the nuclear dynamics over a broad temporal range (0.8-3.8 fs). The intensity of the driving field is found to have a significant impact on the observed dynamics through the vibrational-state dependence of the strong-field ionization rates. We derive a general model that includes these effects and establishes a new link between high-harmonic spectroscopy and classical photoelectron spectroscopy. Our model reproduces the observed dynamics and their dependence on the intensity of the driving field. Moreover, the model predicts much richer nuclear dynamics on the few-fs timescale than most previous theories. The newly predicted features are shown to reflect the quantized vibronic level structure of the molecular cation.

Application of discrete variable representation to planar H2+ in strong xuv laser fields

We present an efficient and accurate grid method to study the strong field dynamics of planar H(2)(+) under Born-Oppenheimer approximation. After introducing the elliptical coordinates to the planar H(2)(+), we show that the Coulomb singularities at the nuclei can be successfully overcome so that both bound and continuum states can be accurately calculated by the method of separation of variables. The time-dependent Schrödinger equation (TDSE) can be accurately solved by a two-dimensional discrete variable representation (DVR) method, where the radial coordinate is discretized with the finite-element discrete variable representation for easy parallel computation and the angular coordinate with the trigonometric DVR which can describe the periodicity in this direction. The bound states energies can be accurately calculated by the imaginary time propagation of TDSE, which agree very well with those computed by the separation of variables. We apply the TDSE to study the ionization dynamics of the planar H(2)(+) by short extreme ultra-violet (xuv) pulses, in which case the differential momentum distributions from both the length and the velocity gauge agree very well with those calculated by the lowest order perturbation theory.

Strong field ionization to multiple electronic states in water

High harmonic spectra show that laser-induced strong field ionization of water has a significant contribution from an inner-valence orbital. Our experiment uses the ratio of H(2)O and D(2)O high harmonic yields to isolate the characteristic nuclear motion of the molecular ionic states. The nuclear motion initiated via ionization of the highest occupied molecular orbital (HOMO) is small and is expected to lead to similar harmonic yields for the two isotopes. In contrast, ionization of the second least bound orbital (HOMO-1) exhibits itself via a strong bending motion which creates a significant isotope effect. We elaborate on this interpretation by solving the time-dependent Schrödinger equation to simulate strong field ionization and high harmonic generation from the water isotopes. We expect that this isotope marking scheme for probing excited ionic states in strong field processes can be generalized to other molecules.

Two-center interference in high-order harmonic generation from heteronuclear diatomic molecules

Two-center interference for heteronuclear diatomic molecules (HeDM) is investigated. The minimum in the high-order harmonic spectrum, as a consequence of the destructive interference, is shifted to lower harmonic orders compared with that in a homonuclear case. This phenomenon is explained by performing phase analysis. It is found that, for an HeDM, the high harmonic spectrum contains information not only on the internuclear separation but also on the properties of the two separate centers, which implies the potential application of estimating the asymmetry of molecules and judging the linear combination of atomic orbitals (LCAO) for the highest occupied molecular orbital (HOMO). Moreover, the possibility to monitor the evolution of HOMO itself in molecular dynamics is also promised.

High-harmonic homodyne detection of the ultrafast dissociation of Br2 molecules

We report the time-resolved observation of the photodissociation of Br2 using high-harmonic generation (HHG) as a probe. The simultaneous measurement of the high-harmonic and ion yields shows that high harmonics generated by the electronically excited state interfere with harmonics generated by the ground state. The resulting homodyne effect provides a high sensitivity to the excited state dynamics. We present a simple theoretical model that accounts for the main observations. Our experiment paves the way towards the dynamic imaging of molecules using HHG.

Mapping molecular orbital symmetry on high-order harmonic generation spectrum using two-color laser fields

We have measured high-order harmonic generation spectra of D2, N2, and CO2 by mixing orthogonally polarized 800 and 400 nm laser fields. The intensity of the high-harmonic spectrum is modulated as we change the relative phase of the two pulses. For randomly orientated molecules, the phase of the intensity modulation depends on the symmetry of the molecular orbitals from which the high harmonics are emitted. This allows us to identify the symmetry of any orbital that contributes to high-harmonic generation, even without aligning the molecule. Our approach can be a route to imaging dynamical changes in three-dimensional molecular orbitals on a time scale as short as a few hundred attoseconds.

Controlling the interference of multiple molecular orbitals in high-harmonic generation

We demonstrate a new method to investigate the origin of spectral structures in high-harmonic generation. We report detailed measurements of high-harmonic spectra in aligned nitrogen and carbon dioxide molecules. Varying the wavelength and intensity of the generating laser field, we show that the minimum in aligned N2 molecules is nearly unaffected, whereas the minimum in aligned CO2 molecules shifts over more than 15 eV. Our quantitative analysis shows that both the interference of multiple orbitals and their structural characteristics affect the position of the minimum. Our method provides a simple approach to the investigation of the high-harmonic generation process in more complex molecules.

High harmonic spectroscopy of multichannel dynamics in strong-field ionization

We perform high harmonic generation spectroscopy of aligned nitrogen molecules to characterize the attosecond dynamics of multielectron rearrangement during strong-field ionization. We use the spectrum and ellipticity of the harmonic light to reconstruct the relative phase between different ionization continua participating in the ionization, and thus determine the shape and location of the hole left in the molecule by strong-field ionization. Our interferometric technique uses transitions between the ionic states, induced by the laser field on the subcycle time scale.

Extension of high harmonic spectroscopy in molecules by a 1300 nm laser field

The emerging techniques of molecular spectroscopy by high order harmonic generation have hitherto been conducted only with Ti:Sapphire lasers which are restricted to molecules with high ionization potentials. In order to gain information on the molecular structure, a broad enough range of harmonics is required. This implies using high laser intensities which would saturate the ionization of most molecular systems of interest, e.g. organic molecules. Using a laser at 1300 nm, we are able to extend the technique to molecules with relatively low ionization potentials (approximately 11 eV), observing wide harmonic spectra reaching up to 60 eV. This energy range improves spatial resolution of the high harmonic spectroscopy to the point where interference minima in harmonic spectra of N(2)O and C(2)H(2) can be observed.

High harmonic interferometry of multi-electron dynamics in molecules

High harmonic emission occurs when an electron, liberated from a molecule by an incident intense laser field, gains energy from the field and recombines with the parent molecular ion. The emission provides a snapshot of the structure and dynamics of the recombining system, encoded in the amplitudes, phases and polarization of the harmonic light. Here we show with CO(2) molecules that high harmonic interferometry can retrieve this structural and dynamic information: by measuring the phases and amplitudes of the harmonic emission, we reveal 'fingerprints' of multiple molecular orbitals participating in the process and decode the underlying attosecond multi-electron dynamics, including the dynamics of electron rearrangement upon ionization. These findings establish high harmonic interferometry as an effective approach to resolving multi-electron dynamics with sub-Angström spatial resolution arising from the de Broglie wavelength of the recombining electron, and attosecond temporal resolution arising from the timescale of the recombination event.

Attosecond circular dichroism spectroscopy of polyatomic molecules

We describe the roles of multiple electronic continua in high-harmonic generation from aligned molecules. First, we show how the circularity of emitted harmonics tracks the interplay of different electronic continua participating in the nonlinear response. Second, we show that the interplay of different continua can lead to large variations of harmonic phases. Finally, we show how multiple electronic continua allow one to shape the polarization of high harmonics and attosecond pulses.