High-Harmonic Generation Spectroscopy of Gas-Phase Bromoform

High-Harmonic Generation (HHG) spectra of randomly aligned bromoform (CHBr3) molecules have been experimentally measured and theoretically simulated at various laser pulse intensities. From the experiments, we obtained a significant number of harmonics that goes beyond the cutoff limit predicted by the three-step model (3SM) with ionization from HOMO. To interpret the experiment, we resorted to real-time time-dependent configuration interaction with single excitations. We found that electronic bound states provide an appreciable contribution to the harmonics. More in detail, we analyzed the electron dynamics by decomposing the HHG signal in terms of single molecular-orbital contributions, to explain the appearance of harmonics around 20-30 eV beyond the expected cutoff due to HOMO. HHG spectra can be therefore explained by considering the contribution at high energy of HOMO-6 and HOMO-9, thus indicating a complex multiple-orbital strong-field dynamics. However, even though the presence of the bromoform cation should be not enough to produce such a signal, we could not exclude a priori that the origin of harmonics in the H29-H45 to be due to the cation, which has more energetic ionization channels.

Implementation of the Time-Dependent Complete-Active-Space Self-Consistent-Field Method for Diatomic Molecules

We present a computational approach that implements the time-dependent complete-active-space self-consistent-field method, as introduced in [Phys. Rev. A 88, 023402 (2013)]. Our implementation addresses the challenge of diatomic molecules subjected to an intense laser pulse by considering the full dimensionality of the problem using prolate spheroidal coordinates. The method incorporates the gauge-invariant frozen-core approximation, boosts the evaluation of the electron-electron interaction term using finite-element discrete-variable representation with Neumann expansion, and utilizes an exponential time differencing scheme tailored for the stable propagation of the stiff nonlinear orbital functions. We have successfully applied this methodology to study high-harmonic generation in diatomic molecules such as H2, LiH, and N2, shedding light on the impact of electron correlations in these systems.

Quantitative performance analysis and comparison of optimal-continuum Gaussian basis sets for high-harmonic generation spectra

Quantum-chemistry methods in the time domain with Gaussian basis sets are increasingly used to compute high-harmonic generation (HHG) spectra of atomic and molecular systems. The quality of these approaches is limited by the accuracy of Gaussian basis sets to describe continuum energy states. In the literature, optimal-continuum Gaussian basis sets have been proposed: Kaufmann et al. [J. Phys. B: At., Mol. Opt. Phys. 22, 2223 (1989)], Woźniak et al. [J. Chem. Phys. 154, 094111 (2021)], Nestmann and Peyerimhoff [J. Phys. B: At., Mol. Opt. Phys. 23, L773 (1990)], Faure et al. [Comput. Phys. Commun. 144, 224 (2002)], and Krause et al. [J. Chem. Phys. 140, 174113 (2014)]. In this work, we have compared the performances of these basis sets to simulate HHG spectra of H atom at different laser intensities. We have also investigated different strategies to balance basis sets with these continuum functions, together with the role of angular momentum. To quantify the performance of the different basis sets, we introduce local and global HHG descriptors. Comparisons with the grid and exact calculations are also provided.

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.

Quantum-Path-Resolved Attosecond High-Harmonic Spectroscopy

Strong-field ionization of molecules releases electrons which can be accelerated and driven back to recombine with their parent ion, emitting high-order harmonics. This ionization also initiates attosecond electronic and vibrational dynamics in the ion, evolving during the electron travel in the continuum. Revealing this subcycle dynamics from the emitted radiation usually requires advanced theoretical modeling. We show that this can be avoided by resolving the emission from two families of electronic quantum paths in the generation process. The corresponding electrons have the same kinetic energy, and thus the same structural sensitivity, but differ by the travel time between ionization and recombination-the pump-probe delay in this attosecond self-probing scheme. We measure the harmonic amplitude and phase in aligned CO_{2} and N_{2} molecules and observe a strong influence of laser-induced dynamics on two characteristic spectroscopic features: a shape resonance and multichannel interference. This quantum-path-resolved spectroscopy thus opens wide prospects for the investigation of ultrafast ionic dynamics, such as charge migration.

High harmonic spectra computed using time-dependent Kohn-Sham theory with Gaussian orbitals and a complex absorbing potential

High harmonic spectra for H₂ and H₂ ⁺ are simulated by solving the time-dependent Kohn-Sham equation in the presence of a strong laser field using an atom-centered Gaussian representation of the density and a complex absorbing potential. The latter serves to mitigate artifacts associated with the finite extent of the basis functions, including spurious reflection of the outgoing electronic wave packet. Interference between the outgoing and reflected waves manifests as peak broadening in the spectrum as well as the appearance of spurious high-energy peaks after the harmonic progression has terminated. We demonstrate that well-resolved spectra can be obtained through the use of an atom-centered absorbing potential. As compared to grid-based algorithms, the present approach is more readily extensible to larger molecules.

Effects of electronic correlation on the high harmonic generation in helium: a time-dependent configuration interaction singles vs time-dependent full configuration interaction study

In this paper, we investigate the effects of full electronic correlation on the high harmonic generation in the helium atom subjected to laser pulses of extremely high intensity. To do this, we perform real-time propagations of the helium atom wavefunction using quantum chemistry methods coupled to Gaussian basis sets. The calculations are done within the real-time time-dependent configuration interaction framework, at two levels of theory: time-dependent configuration interation with single excitations (TD-CIS, uncorrelated method) and time-dependent full configuration interaction (TD-FCI, fully correlated method). The electronic wavefunction is expanded in Dunning basis sets supplemented with functions adapted to describing highly excited and continuum states. We also compare the TD-CI results with grid-based propagations of the helium atom within the single-active-electron approximation. Our results show that when including the dynamical electron correlation, a noticeable improvement to the description of HHG can be achieved, in terms of e.g. a more constant intensity in the lower energy part of the harmonic plateau. However, such effects can be captured only if the basis set used suffices to reproduce the most basic features, such as the HHG cutoff position, at the uncorrelated level of theory.

Continuum Electronic States: The Tiresia Code

A multicenter (LCAO) B-spline basis is described in detail, and its capabilities concerning affording convergent solutions for electronic continuum states and wavepacket propagation are presented. It forms the core of the Tiresia code, which implements static-DFT and TDDFT hamiltonians, as well as single channel Dyson-DFT and Dyson-TDDFT descriptions to include correlation in the bound states. Together they afford accurate and computationally efficient descriptions of photoionization properties of complex systems, both in the single photon and strong field environments. A number of examples are provided.

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.

Dialogue on analytical and ab initio methods in attoscience

The perceived dichotomy between analytical and ab initio approaches to theory in attosecond science is often seen as a source of tension and misconceptions. This Topical Review compiles the discussions held during a round-table panel at the 'Quantum Battles in Attoscience' cecam virtual workshop, to explore the sources of tension and attempt to dispel them. We survey the main theoretical tools of attoscience-covering both analytical and numerical methods-and we examine common misconceptions, including the relationship between ab initio approaches and the broader numerical methods, as well as the role of numerical methods in 'analytical' techniques. We also evaluate the relative advantages and disadvantages of analytical as well as numerical and ab initio methods, together with their role in scientific discovery, told through the case studies of two representative attosecond processes: non-sequential double ionisation and resonant high-harmonic generation. We present the discussion in the form of a dialogue between two hypothetical theoreticians, a numericist and an analytician, who introduce and challenge the broader opinions expressed in the attoscience community.

Quantum electronic coherences by attosecond transient absorption spectroscopy: ab initio B-spline RCS-ADC study

Here I present a fully ab initio time-resolved study of X-ray attosecond transient absorption spectroscopy (ATAS) in a prototypical polyatomic molecule, pyrazine, and demonstrate the possibility of retrieving the many-electron quantum ionic coherences arising in attosecond molecular photoionisation and pre-determining the subsequent charge-directed photochemical reactivity. Advanced first-principles many-electron simulations are performed, within a hybrid XUV pump/X-ray probe setup, to describe the interaction of pyrazine with both XUV pump and X-ray probe pulses, and study the triggered correlated many-electron dynamics. The calculations are carried out by means of the recently-developed ab initio method for many-electron dynamics in polyatomic molecules, the time-dependent (TD) B-spline Restricted Correlation Space-Algebraic Diagrammatic Construction (RCS-ADC). RCS-ADC simulates molecular ionisation from first principles, combining the accurate description of electron correlation of quantum chemistry with the full account of the continuum dynamics of the photoelectron. Complete theoretical characterisation of the atto-ionised many-electron state and photo-induced attosecond charge dynamics is achieved by calculating the reduced ionic density matrix (R-IDM) of the bipartite ion-photoelectron system, with full inclusion of the correlated shakeup states. Deviations from the sudden approximation picture of photoionisation, in the low-photon-energy limit, are presented. The effect of the multi-channel interaction between the parent-ion and the emitted photoelectron on the onset of the quantum electronic coherences is analysed. Moreover, I show how the Schmidt decomposition of the R-IDM unravels the many-electron dynamics triggered by the pump, allowing for the identification of the key channels involved. Finally, I calculate the X-ray attosecond transient absorption spectra of XUV-ionised pyrazine. The results unveil the mapping of the ATAS measurement onto the quantum electronic coherences, and related non-zero R-IDM matrix elements, produced upon ionisation by the XUV pump laser pulse.

A systematic construction of Gaussian basis sets for the description of laser field ionization and high-harmonic generation

A precise understanding of mechanisms governing the dynamics of electrons in atoms and molecules subjected to intense laser fields has a key importance for the description of attosecond processes such as the high-harmonic generation and ionization. From the theoretical point of view, this is still a challenging task, as new approaches to solve the time-dependent Schrödinger equation with both good accuracy and efficiency are still emerging. Until recently, the purely numerical methods of real-time propagation of the wavefunction using finite grids have been frequently and successfully used to capture the electron dynamics in small one- or two-electron systems. However, as the main focus of attoscience shifts toward many-electron systems, such techniques are no longer effective and need to be replaced by more approximate but computationally efficient ones. In this paper, we explore the increasingly popular method of expanding the wavefunction of the examined system into a linear combination of atomic orbitals and present a novel systematic scheme for constructing an optimal Gaussian basis set suitable for the description of excited and continuum atomic or molecular states. We analyze the performance of the proposed basis sets by carrying out a series of time-dependent configuration interaction calculations for the hydrogen atom in fields of intensity varying from 5 × 10¹³ W/cm² to 5 × 10¹⁴ W/cm². We also compare the results with the data obtained using Gaussian basis sets proposed previously by other authors.

Probing the molecular frame of uracil and thymine with high-harmonic generation spectroscopy

In this work, we present the computed high-harmonic generation (HHG) spectra of uracil and thymine molecules, by means of the real-time time-dependent formulation of Gaussian-based configuration interaction with single excitations (RT-TD-CIS). According to the experimental work [Hutchison et al., Comparison of high-order harmonic generation in uracil and thymine ablation plumes, Phys. Chem. Chem. Phys., 2013, 15, 12308], a pulse wavelength of 780 nm has been used, together with an intensity of 1014 W cm-2 and a pulse duration of 23 optical cycles. In order to examine the effect of pulse polarisation, rotationally averaged (to mimic the gas-phase sample) and single-polarisations have been computed for both molecules. Our results show that the HHG signal for both molecules possibly originates from different ionisation channels, involving HOMO, HOMO-1, HOMO-2 and HOMO-3 orbitals, which lie within 4 eV. We characterize the HHG spectrum of thymine, supporting the idea that the absence of the thymine signal in the original work does not depend on the single-molecule behaviour. Present results for uracil are consistent with the experimental data. Moreover, we have observed that states below and above the chosen ionisation threshold provide different contributions to the HHG spectrum in averaged and single-polarisation calculations.

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.

Time-resolved recombination by attosecond-controlled high harmonic generation

We theoretically investigate the coherent control of strong-field high-harmonic generation in the presence of an isolated attosecond pulse. It is found that the rapid modulation of the controlled signal exhibits interference fringe structures in the delay-dependent spectra. By comparing the classical trajectory model with quantum mechanical calculation, it is demonstrated that the fringes are resulted from the interference between the photon- and the tunnelling-initiated recombination pathways. The relative recombination times for the two paths are reconstructed from the interference fringes, which provides a novel scheme for optical observation of the interplay of the photionization and tunneling ionization electron dynamics in attosecond resolution.

Ultrafast Quantum Interference in the Charge Migration of Tryptophan

Extreme-ultraviolet-induced charge migration in biorelevant molecules is a fundamental step in the complex path leading to photodamage. In this work we propose a simple interpretation of the charge migration recently observed in an attosecond pump-probe experiment on the amino acid tryptophan. We find that the decay of the prominent low-frequency spectral structure with increasing pump-probe delay is due to a quantum beating between two geometrically distinct, almost degenerate charge oscillations. Quantum beating is ubiquitous in these systems, and at least on the few-to-tens of femtosecond time scales, it may dominate over decoherence the line intensities of time-resolved spectra. We also address the experimentally observed phase shift in the charge oscillations of two different amino acids, tryptophan and phenylalanine. Our results indicate that a beyond mean-field treatment of the electron dynamics is necessary to reproduce the correct behavior.

Onset of ionic coherence and ultrafast charge dynamics in attosecond molecular ionisation

Here is presented a fully ab initio theoretical framework for simulating the correlated many-electron dynamics occurring during and emerging from molecular ionisation by attosecond laser pulses. This is based on the time-dependent (TD) version of the B-spline restricted correlation space (RCS)-algebraic diagrammatic construction (ADC) method, with the full description of the photoelectron and inclusion of electron correlation effects, such as shakeup processes and inter-channel couplings. The nature of the ultrafast charge dynamics in the molecular ion is elucidated by quantitatively predicting the degree of electronic coherence and eigenstate content of the prepared molecular cationic state, beyond the commonly used sudden approximation. The results presented here for the acetylene and ethylene molecules show that even in the high photon energy regime the simulated hole dynamics is quantitatively different from the prediction of the sudden approximation. Moreover, for high-bandwidth ionising pulse, the residual interaction between the cation, in highly-excited shake-up states, and the emitted slow photoelectron gives rise to a loss of coherence in the ionic system which can persist for the first few femtoseconds after ionisation.

Restricted Correlation Space B-Spline ADC Approach to Molecular Ionization: Theory and Applications to Total Photoionization Cross-Sections

Herein is presented a new approach to the ab initio algebraic diagrammatic construction (ADC) schemes for the polarization propagator, which is explicitly designed to accurately and efficiently describe molecular ionization. The restricted correlation space (RCS) version of the ADC methods up to second order of perturbation theory is derived via the intermediate state representation (ISR) and implemented in the multicenter B-spline basis set for the electronic continuum. Remarkably a general close-coupling structure of the RCS-ADC many-electron wave function, connecting the N-particle to the ( N - 1)-particle ADC intermediate states, emerges naturally as a nontrivial result of the RCS ansatz. Moreover, the introduced RCS-ADC schemes prove to be particularly manageable from a computational point of view, overcoming the practical limitations of the conventional ADC approaches. The quality of the new RCS-ADC( n) approaches is verified by performing a series of total photoionization cross-section calculations on a test set of molecules. The excellent agreement of the results with existing accurate benchmarks demonstrates that the RCS versions of the ADC schemes are optimal and quantitatively accurate methods for studying multichannel molecular photoionization.

CHEERS: a tool for correlated hole-electron evolution from real-time simulations

We put forward a practical nonequilibrium Green's function (NEGF) scheme to perform real-time evolutions of many-body interacting systems driven out of equilibrium by external fields. CHEERS is a computational tool to solve the NEGF equation of motion in the so called generalized Kadanoff-Baym ansatz and it can be used for model systems as well as first-principles Hamiltonians. Dynamical correlation (or memory) effects are added to the Hartree-Fock dynamics through a many-body self-energy. Applications to time-dependent quantum transport, time-resolved photoabsorption and other ultrafast phenomena are discussed.

Zombie states for description of structure and dynamics of multi-electron systems

Canonical Coherent States (CSs) of Harmonic Oscillator have been extensively used as a basis in a number of computational methods of quantum dynamics. However, generalising such techniques for fermionic systems is difficult because Fermionic Coherent States (FCSs) require complicated algebra of Grassmann numbers not well suited for numerical calculations. This paper introduces a coherent antisymmetrised superposition of "dead" and "alive" electronic states called here Zombie State (ZS), which can be used in a manner of FCSs but without Grassmann algebra. Instead, for Zombie States, a very simple sign-changing rule is used in the definition of creation and annihilation operators. Then, calculation of electronic structure Hamiltonian matrix elements between two ZSs becomes very simple and a straightforward technique for time propagation of fermionic wave functions can be developed. By analogy with the existing methods based on Canonical Coherent States of Harmonic Oscillator, fermionic wave functions can be propagated using a set of randomly selected Zombie States as a basis. As a proof of principles, the proposed Coupled Zombie States approach is tested on a simple example showing that the technique is exact.