Three-step model for high-harmonic generation in many-electron systems

Time-dependent ab initio molecular-orbital decomposition for high-harmonic generation spectroscopy

We propose a real-time time-dependent ab initio approach within a configuration-interaction-singles ansatz to decompose the high-harmonic generation (HHG) signal of molecules in terms of individual molecular-orbital (MO) contributions. Calculations have been performed by propagating the time-dependent Schrödinger equation with complex energies, in order to account for ionization of the system, and by using tailored Gaussian basis sets for high-energy and continuum states. We have studied the strong-field electron dynamics and the HHG spectra in aligned CO2 and H2O molecules. Contribution from MOs in the strong-field dynamics depends on the interplay between the MO ionization energy and the coupling between the MO and the laser-pulse symmetries. Such contributions characterize different portions of the HHG spectrum, indicating that the orbital decomposition encodes nontrivial information on the modulation of the strong-field dynamics. Our results correctly reproduce the MO contributions to HHG for CO2 as described in the literature experimental and theoretical data and lead to an original analysis of the role of the highest occupied molecular orbitals HOMO, HOMO-1, and HOMO-2 of H2O according to the polarization direction of the laser pulse.

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.

Role of Inner Molecular Orbitals in High-Harmonic Generation Spectra of Aligned Uracil

In this work, we decompose the high-harmonic generation (HHG) signal of aligned gas-phase uracil into single molecular-orbital (MO) contributions. We compute HHG spectra for a pulse linearly polarized perpendicular to the molecular plane, with an intensity of 0.6 and 0.85 × 1014 W/cm2 and a wavelength of 800 nm. We use the real-time time-dependent Configuration Interaction with singles method, coupled to a Gaussian-based representation of the time-dependent wavefunction. The strong-field dynamics is affected by the energy of the ionization/recombination channels and by the coupling between the orbital symmetry and laser polarization. In the configuration studied here, we expect that π-type MOs favorably couple with the incoming pulse and play a substantial role in generating the HHG spectrum. Indeed, we show that HOMO, HOMO - 1, and HOMO - 4, which all are π-like, determine the intensity of harmonic peaks at different energies, while HOMO - 2 and HOMO - 3 provide a smaller contribution. It is worth mentioning that HOMO - 4 produces a stronger signal than that from HOMO - 1, even though the corresponding ionization energy, in an one-electron picture, is around 2.5 eV larger and more than 4 eV larger than the HOMO one.

Extreme ultraviolet pellicle wrinkles influence on mask 3D effects: experimental demonstration

Extreme ultraviolet (EUV) lithography uses reflective optics and a thick mask absorber, leading to mask 3D (M3D) effects. These M3D effects cause disparities in the amplitudes and phases of EUV mask diffractions, impacting mask imaging performance and reducing process yields. Our findings demonstrate that wrinkles in the EUV pellicle can exacerbate M3D effects. This imbalance results in critical dimension variation, image contrast loss, and pattern shift in mask images. Therefore, the use of a pellicle material with thermodynamic characteristics that minimize wrinkles when exposed to EUV rays is imperative.

Efficient low-density grating setup for monochromatization of XUV ultrafast light sources

Ultrafast light sources have become an indispensable tool to access and understand transient phenomenon in material science. However, a simple and easy-to-implement method for harmonic selection, with high transmission efficiency and pulse duration conservation, is still a challenge. Here we showcase and compare two approaches for selecting the desired harmonic from a high harmonic generation source while achieving the above goals. The first approach is the combination of extreme ultraviolet spherical mirrors with transmission filters and the second approach uses a normal-incidence spherical grating. Both solutions target time- and angle-resolved photoemission spectroscopy with photon energies in the 10-20 eV range but are relevant for other experimental techniques as well. The two approaches for harmonic selection are characterized in terms of focusing quality, photon flux, and temporal broadening. It is demonstrated that a focusing grating is able to provide much higher transmission as compared to the mirror+filter approach (3.3 times higher for 10.8 eV and 12.9 times higher for 18.1 eV), with only a slight temporal broadening (6.8% increase) and a somewhat larger spot size (∼30% increase). Overall, our study establishes an experimental perspective on the trade-off between a single grating normal incidence monochromator design and the use of filters. As such, it provides a basis for selecting the most appropriate approach in various fields where an easy-to-implement harmonic selection from high harmonic generation is needed.

Investigation of the Spatio-Temporal Characteristics of High-Order Harmonic Generation Using a Bohmian Trajectory Scheme

High-order harmonic generation of atoms irradiated by an ultrashort laser pulse was calculated by numerically solving the time-dependent Schrödinger equation and the Bohmian trajectory scheme. The harmonic spectra with the two schemes are quantitatively consistent. Using the wavelet behavior of the Bohmian trajectory, the spatio-temporal features of harmonic emission from different energy regions are analyzed. It is found that the spatio-temporal distribution of the harmonic well revealed the physical mechanism of harmonic generation. This method will contribute to the understanding of harmonic emission mechanisms in complex systems, which include many atoms.

Quantum-Compute Algorithm for Exact Laser-Driven Electron Dynamics in Molecules

In this work, we investigate the capability of known quantum computing algorithms for fault-tolerant quantum computing to simulate the laser-driven electron dynamics of excitation and ionization processes in small molecules such as lithium hydride, which can be benchmarked against the most accurate time-dependent full configuration interaction (TD-FCI) calculations. The conventional TD-FCI wave packet propagation is reproduced using the Jordan-Wigner transformation for wave function and operators and the Trotter product formula for expressing the propagator. In addition, the time-dependent dipole moment, as an example of a time-dependent expectation value, is calculated using the Hadamard test. To include non-Hermitian operators in the ionization dynamics, a similar approach to the quantum imaginary time evolution (QITE) algorithm is employed to translate the propagator, including a complex absorption potential, into quantum gates. The computations are executed on a quantum computer simulator. By construction, all quantum computer algorithms, except for the QITE algorithm used only for ionization but not for excitation dynamics, would scale polynomially on a quantum computer with fully entangled qubits. In contrast, TD-FCI scales exponentially. Hence, quantum computation holds promises for substantial progress in the understanding of electron dynamics of excitation processes in increasingly large molecular systems, as has already been witnessed in electronic structure theory.

Effect of wrinkles on extreme ultraviolet pellicle reflectivity and local critical dimension

Extreme ultraviolet (EUV) pellicles must have an EUV reflectance (EUVR) below 0.04% to prevent the reduction of critical dimension (CD). However, pellicle wrinkles cause localized CD variation by locally amplifying the EUVR. This study demonstrates that wrinkles can increase the pellicle's EUVR by approximately four times, and the CD drop depends on the relative position of the reflected light from the wrinkle to the 0th- or 1st-order diffracted light. The CD decreases by 6 nm. Therefore, even if the pellicle satisfies the requirement for the EUVR, we need to tightly control the generation of wrinkles to suppress CD variation during the entire exposure process.

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.

Time-dependentab initioapproaches for high-harmonic generation spectroscopy

High-harmonic generation (HHG) is a nonlinear physical process used for the production of ultrashort pulses in XUV region, which are then used for investigating ultrafast phenomena in time-resolved spectroscopies. Moreover, HHG signal itself encodes information on electronic structure and dynamics of the target, possibly coupled to the nuclear degrees of freedom. Investigating HHG signal leads to HHG spectroscopy, which is applied to atoms, molecules, solids and recently also to liquids. Analysing the number of generated harmonics, their intensity and shape gives a detailed insight of, e.g., ionisation and recombination channels occurring in the strong-field dynamics. A number of valuable theoretical models has been developed over the years to explain and interpret HHG features, with the three-step model being the most known one. Originally, these models neglect the complexity of the propagating electronic wavefunction, by only using an approximated formulation of ground and continuum states. Many effects unravelled by HHG spectroscopy are instead due to electron correlation effects, quantum interference, and Rydberg-state contributions, which are all properly captured by anab initioelectronic-structure approach. In this review we have collected recent advances in modelling HHG by means ofab initiotime-dependent approaches relying on the propagation of the time-dependent Schrödinger equation (or derived equations) in presence of a very intense electromagnetic field. We limit ourselves to gas-phase atomic and molecular targets, and to solids. We focus on the various levels of theory employed for describing the electronic structure of the target, coupled with strong-field dynamics and ionisation approaches, and on the basis used to represent electronic states. Selected applications and perspectives for future developments are also given.

Role of exchange and correlation in high-harmonic generation spectra of H2, N2, and CO2: Real-time time-dependent electronic-structure approaches

This study arises from the attempt to answer the following question: how different descriptions of electronic exchange and correlation affect the high-harmonic generation (HHG) spectroscopy of H₂, N₂, and CO₂ molecules? We compare HHG spectra for H₂, N₂, and CO₂ with different ab initio electronic structure methods: real-time time-dependent configuration interaction and real-time time-dependent density functional theory (RT-TDDFT) using truncated basis sets composed of correlated wave functions expanded on Gaussian basis sets. In the framework of RT-TDDFT, we employ Perdew-Burke-Ernzerhof (PBE) and long-range corrected Perdew-Burke-Ernzerhof (LC-ωPBE) functionals. We study HHG spectroscopy by disentangling the effect of electronic exchange and correlation. We first analyze the electronic exchange alone, and in the case of RT-TDDFT with LC-ωPBE, we use ω = 0.3 and ω = 0.4 to tune the percentage of long-range Hartree-Fock exchange and short-range exchange PBE. Then, we added the correlation as described by the PBE functional. All the methods give very similar HHG spectra, and they seem not to be particularly sensitive to the different description of exchange and correlation or to the correct asymptotic behavior of the Coulomb potential. Despite this general trend, some differences are found in the region connecting the cutoff and the background. Here, the harmonics can be resolved with different accuracy depending on the theoretical schemes used. We believe that the investigation of the molecular continuum and its coupling with strong fields merits further theoretical investigations in the near future.

Photoinduced valence tautomerism of a cobalt-dioxolene complex revealed with femtosecond M-edge XANES

Cobalt complexes that undergo charge-transfer induced spin-transitions or valence tautomerism from low spin Co^III to high spin (HS) Co^II are potential candidates for magneto-optical switches. We use M_2,3-edge X-ray absorption near-edge structure (XANES) spectroscopy with 40 fs time resolution to measure the excited-state dynamics of Co^III(Cat-N-SQ)(Cat-N-BQ), where Cat-N-BQ and Cat-N-SQ are the singly and doubly reduced forms of the 2-(2-hydroxy-3,5-di-tert-butylphenyl-imino)-4,6-di-tert-butylcyclohexa-3,5-dienone ligand. The extreme ultraviolet probe pulses, produced using a tabletop high-harmonic generation light source, measure 3p → 3d transitions and are sensitive to the spin and oxidation state of the Co center. Photoexcitation at 525 nm produces a low-spin Co^II ligand-to-metal charge transfer state which undergoes intersystem crossing to high-spin Co^II in 67 fs. Vibrational cooling from this hot HS Co^II state competes on the hundreds-of-fs time scale with back-intersystem crossing to the ground state, with 60% of the population trapped in a cold HS Co^II state for 24 ps. Ligand field multiplet simulations accurately reproduce the ground-state spectra and support the excited-state assignments. This work demonstrates the ability of M_2,3-edge XANES to measure ultrafast photophysics of molecular Co complexes.

Charge migration and charge transfer in molecular systems

The transfer of charge at the molecular level plays a fundamental role in many areas of chemistry, physics, biology and materials science. Today, more than 60 years after the seminal work of R. A. Marcus, charge transfer is still a very active field of research. An important recent impetus comes from the ability to resolve ever faster temporal events, down to the attosecond time scale. Such a high temporal resolution now offers the possibility to unravel the most elementary quantum dynamics of both electrons and nuclei that participate in the complex process of charge transfer. This review covers recent research that addresses the following questions. Can we reconstruct the migration of charge across a molecule on the atomic length and electronic time scales? Can we use strong laser fields to control charge migration? Can we temporally resolve and understand intramolecular charge transfer in dissociative ionization of small molecules, in transition-metal complexes and in conjugated polymers? Can we tailor molecular systems towards specific charge-transfer processes? What are the time scales of the elementary steps of charge transfer in liquids and nanoparticles? Important new insights into each of these topics, obtained from state-of-the-art ultrafast spectroscopy and/or theoretical methods, are summarized in this review.

Attosecond physics at the nanoscale

Recently two emerging areas of research, attosecond and nanoscale physics, have started to come together. Attosecond physics deals with phenomena occurring when ultrashort laser pulses, with duration on the femto- and sub-femtosecond time scales, interact with atoms, molecules or solids. The laser-induced electron dynamics occurs natively on a timescale down to a few hundred or even tens of attoseconds (1 attosecond  =  1 as  =  10^-18 s), which is comparable with the optical field. For comparison, the revolution of an electron on a 1s orbital of a hydrogen atom is  ∼152 as. On the other hand, the second branch involves the manipulation and engineering of mesoscopic systems, such as solids, metals and dielectrics, with nanometric precision. Although nano-engineering is a vast and well-established research field on its own, the merger with intense laser physics is relatively recent. In this report on progress we present a comprehensive experimental and theoretical overview of physics that takes place when short and intense laser pulses interact with nanosystems, such as metallic and dielectric nanostructures. In particular we elucidate how the spatially inhomogeneous laser induced fields at a nanometer scale modify the laser-driven electron dynamics. Consequently, this has important impact on pivotal processes such as above-threshold ionization and high-order harmonic generation. The deep understanding of the coupled dynamics between these spatially inhomogeneous fields and matter configures a promising way to new avenues of research and applications. Thanks to the maturity that attosecond physics has reached, together with the tremendous advance in material engineering and manipulation techniques, the age of atto-nanophysics has begun, but it is in the initial stage. We present thus some of the open questions, challenges and prospects for experimental confirmation of theoretical predictions, as well as experiments aimed at characterizing the induced fields and the unique electron dynamics initiated by them with high temporal and spatial resolution.

Chirp-free isolated attosecond pulse generation from an atom irradiated by a fundamental terahertz pulse synchronizing an infrared laser pulse

We theoretically study high-order harmonic generation (HHG) and attosecond pulses from an atom irradiated synchronically by a terahertz (THz) pulse and an infrared laser pulse. For the HHG spectrum from the THz pulse alone and the combined pulse, an apparent peak-valley structure appears the platform region. Specially, for the periodic structure generated by an atom under the action of the combined pulse is originated from the interference between the electrons ionized at different instants in the laser field, which undergo many recollision and return to the core at the same time. Therefore, continuum harmonics with few chirps from the interference enhancement region can be achieved, which result in a chirp-free isolated attosecond pulse.

Influence of many-body interactions during the ionization of gases by short intense optical pulses

The excitation of atomic gases by short high-intensity optical pulses leads to significant electron ionization. In dilute systems, the generated distribution of ionized electrons is highly anisotropic, reflecting the quantum mechanical properties of the atomic states involved in the many photon transitions. For higher atomic densities, the Coulomb interaction in the electron-ion system leads to the development of an isotropic electron plasma. To study the ionization process in the presence of the many-body interaction, a fully microscopic model is developed that combines a generalized version of the optical Bloch equations describing the optical excitation with a microscopic description of the many-body interactions. The numerical evaluation shows that the Coulomb interaction significantly modifies the distribution anisotropy already during the excitation process. Whereas a reduced anisotropy is still present after the pulse for low ionization degrees and pressures, it is completely absent for elevated gas densities. An ionization degree is predicted that is significantly enhanced by the many-body interactions.

High-harmonic spectroscopy of oriented OCS molecules: emission of even and odd harmonics

We study the emission of even and odd high-harmonic orders from oriented OCS molecules. We use an intense, nonresonant femtosecond laser pulse superimposed with its phase-controlled second harmonic field to impulsively align and orient a dense sample of molecules from which we subsequently generate high-order harmonics. The even harmonics appear around the full revivals of the rotational dynamics. We demonstrate perfect coherent control over their intensity through the subcycle delay of the two-color fields. The odd harmonics are insensitive to the degree of orientation, but modulate with the degree of axis alignment, in agreement with calculated photorecombination dipole moments. We further compare the shape of the even and odd harmonic spectra with our calculations and determine the degree of orientation.

Imaging orbitals with attosecond and Ångström resolutions: toward attochemistry?

The recently developed attosecond light sources make the investigation of ultrafast processes in matter possible with unprecedented time resolution. It has been proposed that the very mechanism underlying the attosecond emission allows the imaging of valence orbitals with Ångström space resolution. This controversial idea together with the possibility of combining attosecond and Ångström resolutions in the same measurements has become a hot topic in strong-field science. Indeed, this could provide a new way to image the evolution of the molecular electron cloud during, e.g. a chemical reaction in 'real time'. Here we review both experimental and theoretical challenges raised by the implementation of these prospects. In particular, we show how the valence orbital structure is encoded in the spectral phase of the recombination dipole moment calculated for Coulomb scattering states, which allows a tomographic reconstruction of the orbital using first-order corrections to the plane-wave approach. The possibility of disentangling multi-channel contributions to the attosecond emission is discussed as well as the necessary compromise between the temporal and spatial resolutions.

Analytic scaling analysis of high harmonic generation conversion efficiency

Closed form expressions for the high harmonic generation (HHG) conversion efficiency are obtained for the plateau and cutoff regions. The presented formulas eliminate most of the computational complexity related to HHG simulations, and enable a detailed scaling analysis of HHG efficiency as a function of drive laser parameters and material properties. Moreover, in the total absence of any fitting procedure, the results show excellent agreement with experimental data reported in the literature. Thus, this paper opens new pathways for the global optimization problem of extreme ultraviolet (EUV) sources based on HHG.

Multielectron correlation in high-harmonic generation: a 2D model analysis

We analyze the role of multielectron dynamics in high-harmonic generation spectroscopy, using an example of a two-electron system. We identify and systematically quantify the importance of correlation and exchange effects. One of the main sources for correlation is identified to be the polarization of the ion by the recombining continuum electron. This effect, which plays an important qualitative and quantitative role, seriously undermines the validity of the standard approaches to high-harmonic generation, which ignore the contribution of excited ionic states to the radiative recombination of the continuum electron.

Observation of electronic structure minima in high-harmonic generation

We report detailed measurements of the high-harmonic spectra generated from argon atoms. The spectra exhibit a deep minimum that is shown to be independent of the laser intensity, and is thus a clear measure of the electronic structure of the atom. We show that exact field-free continuum wave functions reproduce the minimum, but plane wave and Coulomb wave functions do not. This remarkable observation suggests that electronic structure can be accurately determined in high-harmonic experiments despite the presence of the strong laser field. Our results clarify the relation between high-harmonic generation and photoelectron spectroscopy. The use of exact continuum functions also resolves the ambiguity associated with the choice of the dispersion relation.

High harmonic generation and molecular orbital tomography in multielectron systems

High harmonic radiation is produced when atoms or molecules are ionized by an intense femtosecond laser pulse. The radiated spectrum has been shown experimentally to contain information on the electronic structure of the molecule, which can be interpreted as an image of a single molecular orbital. Previous theory for high harmonic generation has been limited to the single-active-electron approximation. Utilizing semisudden approximation, the authors develop a theory of the recombination step in high harmonic generation and tomographic reconstruction in multielectron systems, taking into account electron spin statistics and electron-electron correlations within the parent molecule and the ion. They show that the resulting corrections significantly modify the theoretical predictions, and bring them in a better agreement with experiment. They further show that exchange contributions to harmonic radiation can be used to extract additional information on the electronic wave function.

Nonsequential double recombination in intense laser fields

A second plateau in the harmonic spectra of laser-driven two-electron atoms is observed both in the numerical solution of a low-dimensional model helium atom and using an extended strong field approximation. It is shown that the harmonics well beyond the usual cutoff are due to the simultaneous recombination of the two electrons, which were emitted during different, previous half-cycles. The new cutoff is explained in terms of classical trajectories. Classical predictions and the time-frequency analysis of the ab initio quantum results are in excellent agreement. The mechanism corresponds to the inverse single photon double ionization process in the presence of a (low frequency) laser field.

Role of many-electron dynamics in high harmonic generation

High harmonic generation (HHG) in many-electron atoms is studied theoretically. The breakdown of the frozen-core single active electron approximation is demonstrated, as it predicts roughly the same radiation amplitude in all noble gases. This is in contradiction with experiments, where heavier noble gases are known to emit much stronger HHG radiation than lighter ones. This experimental behavior of the noble gases can be qualitatively reproduced when many-electron dynamics, within a simple approximation, is taken into account.