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

Modeling of High-Harmonic Generation in the C60 Fullerene Using Ab Initio, DFT-Based, and Semiempirical Methods

We report calculations of the high-harmonic generation spectra of the C60 fullerene molecule carried out by employing a diverse set of real-time time-dependent quantum chemical methods. All methodologies involve expanding the propagated electronic wave function in bases consisting of the ground and singly excited time-independent eigenstates obtained through the solution of the corresponding linear-response equations. We identify the correlation and exchange effect in the spectra by comparing the results from methods relying on the Hartree-Fock reference determinant with those obtained using approaches based on the density functional theory with different exchange-correlation functionals. The effect of the full random-phase approximation treatment of the excited electronic states is also analyzed and compared with the configuration interaction singles and the Tamm-Dancoff approximation. We also showcase the fact that the real-time extension of the semiempirical method INDO/S can be effectively applied for an approximate description of laser-driven dynamics in large systems.

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.

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.

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.

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.

A complex Gaussian approach to molecular photoionization

We develop and implement a Gaussian approach to calculate partial cross-sections and asymmetry parameters for molecular photoionization. Optimal sets of complex Gaussian-type orbitals (cGTOs) are first obtained by nonlinear optimization, to best fit sets of Coulomb or distorted continuum wave functions for relevant orbital quantum numbers. This allows us to represent the radial wavefunction for the outgoing electron with accurate cGTO expansions. Within a time-independent partial wave approach, we show that all the necessary transition integrals become analytical, in both length and velocity gauges, thus facilitating the numerical evaluation of photoionization observables. Illustrative results, presented for NH₃ and H₂ O within a one-active-electron monocentric model, validate numerically the proposed strategy based on a complex Gaussian representation of continuum states.

Gaussian-Type Orbital Calculations for High Harmonic Generation in Vibrating Molecules: Benchmarks for H2. J Chem Theory Comput 2021;17:7353-7365. [PMID: 34747605 DOI: 10.1021/acs.jctc.1c00837]

The response of the hydrogen molecular ion, H₂⁺, to few-cycle laser pulses of different intensities is simulated. To treat the coupled electron-nuclear motion, we use adiabatic potentials computed with Gaussian-type basis sets together with a heuristic ionization model for the electron and a grid representation for the nuclei. Using this mixed-basis approach, the time-dependent Schrödinger equation is solved, either within the Born-Oppenheimer approximation or with nonadiabatic couplings included. The dipole response spectra are compared to all-grid-based solutions for the three-body problem, which we take as a reference to benchmark the Gaussian-type basis set approaches. Also, calculations employing the fixed-nuclei approximation are performed, to quantify effects due to nuclear motion. For low intensities and small ionization probabilities, we get excellent agreement of the dynamics using Gaussian-type basis sets with the all-grid solutions. Our investigations suggest that high harmonic generation (HHG) and high-frequency response, in general, can be reliably modeled using Gaussian-type basis sets for the electrons for not too high harmonics. Further, nuclear motion destroys electronic coherences in the response spectra even on the time scale of about 30 fs and affects HHG intensities, which reflect the electron dynamics occurring on the attosecond time scale. For the present system, non-Born-Oppenheimer effects are small. The Gaussian-based, nonadiabatically coupled, time-dependent multisurface approach to treat quantum electron-nuclear motion beyond the non-Born-Oppenheimer approximation can be easily extended to approximate wavefunction methods, such as time-dependent configuration interaction singles (TD-CIS), for systems where no benchmarks are available.

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.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.