Many-electron dynamics in laser-driven molecules: wavefunction theory vs. density functional theory

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.

Visualizing and characterizing excited states from time-dependent density functional theory

Time-dependent density functional theory (TD-DFT) is the most widely-used electronic structure method for excited states, due to a favorable combination of low cost and semi-quantitative accuracy in many contexts, even if there are well recognized limitations. This Perspective describes various ways in which excited states from TD-DFT calculations can be visualized and analyzed, both qualitatively and quantitatively. This includes not just orbitals and densities but also well-defined statistical measures of electron-hole separation and of Frenkel-type exciton delocalization. Emphasis is placed on mathematical connections between methods that have often been discussed separately. Particular attention is paid to charge-transfer diagnostics, which provide indicators of when TD-DFT may not be trustworthy due to its categorical failure to describe long-range electron transfer. Measures of exciton size and charge separation that are directly connected to the underlying transition density are recommended over more ad hoc metrics for quantifying charge-transfer character.

Tailoring light-induced charge transfer and intersystem crossing in FeCO using time-dependent spin-orbit configuration interaction

Real-time (RT) electronic structure methods provide a natural framework for describing light-matter interactions in arbitrary time-dependent electromagnetic fields (EMF). Optically induced excited state transitions are of particular interest, which require tuned EMF to drive population transfer to and from the specific state(s) of interest. Intersystem crossing, or spin-flip, may be driven through shaped EMF or laser pulses. These transitions can result in long-lived "spin-trapped" excited states, which are especially useful for materials requiring charge separation or protracted excited state lifetimes. Time-dependent configuration interaction (TDCI) is unique among RT methods in that it may be implemented in a basis of eigenstates, allowing for rapid propagation of the time-dependent Schrödinger equation. The recent spin-orbit TDCI (TD-SOCI) enables a real-time description of spin-flip dynamics in an arbitrary EMF and, therefore, provides an ideal framework for rational pulse design. The present study explores the mechanism of multiple spin-flip pathways for a model transition metal complex, FeCO, using shaped pulses designed to drive controlled intersystem crossing and charge transfer. These results show that extremely tunable excited state dynamics can be achieved by considering the dipole transition matrix elements between the states of interest.

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.

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.

Electronic circular dichroism from real-time propagation in state space

In this paper, we propose to compute the electronic circular dichroism (ECD) spectra of chiral molecules using a real-time propagation of the time-dependent Schrödinger equation (TDSE) in the space of electronic field-free eigenstates, by coupling TDSE with a given treatment of the electronic structure of the target. The time-dependent induced magnetic moment is used to compute the ECD spectrum from an explicit electric perturbation. The full matrix representing the transition magnetic moment in the space of electronic states is generated from that among pairs of molecular orbitals. In the present work, we show the ECD spectra of methyloxirane, of several conformers of L-alanine, and of the Λ-Co(acac)₃ complex, computed from a singly excited ansatz of time-dependent density functional theory eigenstates. The time-domain ECD spectra properly reproduce the frequency-domain ones obtained in the linear-response regime and quantitatively agree with the available experimental data. Moreover, the time-domain approach to ECD allows us to naturally go beyond the ground-state rotationally averaged ECD spectrum, which is the standard outcome of the linear-response theory, e.g., by computing the ECD spectra from electronic excited states.

Following the density evolution using real time density functional theory and density based indexes: Application to model push-pull molecules

Considering as test case a family of organic rod like push-pull molecules, we derived and applied density based index enabling the description and diagnostic of the electronic density evolution in real time-time dependent density functional theory (RT-TDDFT) simulations. In particular, both the charge transfer (CT) distance and a diagnostic index, the D_CT and M_AC ^RT respectively, were computed on the fly from the density distribution obtained at a given time and the reference ground state density and their mean values were compared with what obtained at Linear Response-TDDFT level. Besides giving a way of analyzing the density redistribution occurring in time, these tools allowed to show how RT-TDDFT, which is definitely a powerful method to model the evolution of the density in CT or charge separation processes, can be affected by the same artifacts known for LR-TDDFT approaches and, particularly, to those related to the use of approximate exchange correlation functionals. The analysis here performed allowed to identify and discard on fly the electronic configurations corresponding to spurious situations.

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.

Resonantly Enhanced Difference-Frequency Generation in the Core X-ray Absorption of Molecules

We use real-time time-dependent density functional theory simulations to numerically demonstrate that resonantly enhanced difference-frequency generation (re-DFG) involving intense ultrashort coherent X-ray pulses can selectively excite core states of atoms in molecules. As a model case, we evaluate the spectral selectivity of re-DFG excitation of the oxygen K-edge by illumination of a single gas-phase water molecule with two-color X-ray pulses of different photon energies and durations. The re-DFG excitation is further probed by a small delayed pulse with central photon energy resonant with the oxygen K-edge peak absorption line. Based on these results, we anticipate that highly selective excitation by re-DFG X-ray nonlinear processes might be achieved in more complex molecular systems and bulk materials by using highly penetrating two-color hard X-ray pulses, with extensive applications.

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.