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

With recent experimental advances in laser-driven electron dynamics in polyatomic molecules, the need arises for their reliable theoretical modelling. Among efficient, yet fairly accurate methods for many-electron dynamics are Time-Dependent Configuration Interaction Singles (TD-CIS) (a Wave Function Theory (WFT) method), and Real-Time Time-Dependent Density Functional Theory (RT-TD-DFT), respectively. Here we compare TD-CIS combined with extended Atomic Orbital (AO) bases, TD-CIS/AO, with RT-TD-DFT in a grid representation of the Kohn-Sham orbitals, RT-TD-DFT/Grid. Possible ionization losses are treated by complex absorbing potentials in energy space (for TD-CIS/AO) or real space (for RT-TD-DFT), respectively. The comparison is made for two test cases: (i) state-to-state transitions using resonant lasers (π-pulses), i.e., bound electron motion, and (ii) large-amplitude electron motion leading to High Harmonic Generation (HHG). Test systems are a H2 molecule and cis- and trans-1,2-dichlorethene, C2H2Cl2, (DCE). From time-dependent electronic energies, dipole moments and from HHG spectra, the following observations are made: first, for bound state-to-state transitions enforced by π-pulses, TD-CIS nicely accounts for the expected population inversion in contrast to RT-TD-DFT, in agreement with earlier findings. Secondly, when using laser pulses under non-resonant conditions, dipole moments and lower harmonics in HHG spectra are obtained by TD-CIS/AO which are in good agreement with those obtained with RT-TD-DFT/Grid. Deviations become larger for higher harmonics and at low laser intensities, i.e., for low-intensity HHG signals. We also carefully test effects of basis sets for TD-CIS/AO and grid size for RT-TD-DFT/Grid, different exchange-correlation functionals in RT-TD-DFT, and absorbing boundaries. Finally, for the present examples, TD-CIS/AO is observed to be at least an order of magnitude more computationally efficient than RT-TD-DFT/Grid.

SIMULATIONS OF PUMP-PROBE EXCITATIONS OF ELECTRONIC WAVE PACKETS FOR A LARGE QUASI-RIGID MOLECULAR SYSTEM BY MEANS OF AN EXTENSION TO THE TIME-DEPENDENT CONFIGURATION INTERACTION SINGLES METHOD

In this paper, we report simulations of laser-driven many-electron dynamics by means of the time-dependent configuration interaction singles (TD-CIS) approach. Photoionization is included by a heuristic model within calculations employing standard Gaussian basis sets. Benzo[g]-N-methyl-quinolinium-7-hydroxylate (BMQ7H) serves as a test system to generate predefined wave packets, i.e. a superposition between the ground and fifth excited state, in a large molecule. For this molecule, these two states have a very similar geometry, which enables us to use the fixed nuclei approximation. Furthermore, this geometric stability would also prevent a dephasing of the electron wave packet due to nuclear dynamics in an experimental realization of our simulations. We also simulate the possible detection of such a wave packet by ultra short probe laser pulses, i.e. pump-probe spectra.

Molecular response properties from explicitly time-dependent configuration interaction methods

In this paper we report the calculation of molecular electric response properties with the help of explicitly time-dependent configuration interaction (TD-CI) methods. These methods have the advantage of being applicable (within the limitations of the time-dependent Schrodinger equation) to time-dependent perturbations of arbitrary shape and strength. Three variants are used to solve the time-dependent electronic Schrodinger equation, namely, the TD-CIS (inclusion of single excitations only), TD-CISD (inclusion of single and double excitations), and TD-CIS(D) (single excitations and perturbative treatment of double excitations) methods and applied for illustration to small molecules, H(2) and H(2)O. In the calculation, slowly varying off-resonant electric fields are applied to the molecules and linear (polarizabilities) and nonlinear (hyperpolarizabilities, harmonic generation) response properties are determined from the time-dependent dipole moments.