Time-dependent Configuration Interaction Research Articles

Light-induced electron dynamics of molecules in cavities: comparison of model Hamiltonians

The rapidly expanding field of polaritonic chemistry requires accurate theoretical simulations to understand new phenomena at the atomic scale. Computing the optoelectronic properties of molecules using established electronic structure methods is a careful balance of accuracy and computational expense, and expanding these methods to quantum electrodynamics to describe coupled cavity-molecule systems is an active topic of development. Key to these methods are the Hamiltonian operators representing the photon cavity modes. The recently introduced quantum electrodynamics time-dependent configuration interaction (QED-TDCI) method allows for the combination of electron dynamics simulations with quantum electrodynamics, enabling the simulation of time-dependent optoelectronic properties of cavity-molecule systems. Using this method, a comparison of two many-state QED Hamiltonians—the Pauli-Fierz and quantum Rabi model Hamiltonians—is presented, with a particular focus on time-dependent properties in applied electric fields.

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.

Cavity-modified molecular dipole switching dynamics.

Polaritonic states, which are formed by resonances between a molecular excitation and the photonic mode of a cavity, have a number of useful properties that offer new routes to control molecular photochemistry using electric fields. To provide a theoretical description of how polaritonic states affect the real-time electron dynamics in molecules, a new method is described where the effects of strong light-molecule coupling are implemented using real-time electronic structure theory. The coupling between the molecular electronic states and the cavity is described by the Pauli-Fierz Hamiltonian, and transitions between polaritonic states are induced via an external time-dependent electric field using time-dependent configuration interaction (TDCI) theory, producing quantum electrodynamics TDCI (QED-TDCI). This method is used to study laser-induced ultrafast charge transfer and dipole-switching dynamics of the LiCN molecule inside a cavity. The increase in cavity coupling strength is found to have a significant impact on the energies and transition dipole moments of the molecule-cavity system. The convergence of the polaritonic state energies as a function of the number of included electronic and photonic basis states is discussed.

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.

Ultrafast internal conversion and photochromism in gas-phase salicylideneaniline.

Salicylideneaniline (SA) is an archetypal system for excited-state intramolecular proton transfer (ESIPT) in non-planar systems. Multiple channels for relaxation involving both the keto and enol forms have been proposed after excitation to S1 with near-UV light. Here, we present transient absorption measurements of hot gas-phase SA, jet-cooled SA, and SA in Ar clusters using cavity-enhanced transient absorption spectroscopy (CE-TAS). Assignment of the spectra is aided by simulated TAS spectra, computed by applying time-dependent complete active space configuration interaction (TD-CASCI) to structures drawn from nonadiabatic molecular dynamics simulations. We find prompt ESIPT in all conditions followed by the rapid generation of the trans keto metastable photochrome state and fluorescent keto state in parallel. Increasing the internal energy increases the photochrome yield and decreases the fluorescent yield and fluorescent state lifetime observed in TAS. In Ar clusters, internal conversion of SA is severely hindered, but the photochrome yield is unchanged. Taken together, these results are consistent with the photochrome being produced via the vibrationally excited keto population after ESIPT.

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.

Charge Migration in HCCI Cations Probed by Strong Field Ionization: Time-Dependent Configuration Interaction and Vibrational Wavepacket Simulations.

Strong field ionization of neutral iodoacetylene (HCCI) can produce a coherent superposition of the X and A cations and results in charge migration between the CC π orbital and the iodine π-type lone pair. This charge migration causes oscillations in the rate of strong field ionization of the cation to the dication that can be monitored using intense few-cycle probe pulses. The dynamics and strong field ionization of the coherent superposition the X and A states of HCCI+ have been modeled by time-dependent configuration interaction (TDCI) simulations. When the nuclei are allowed to move, the electronic wavefunctions need to be multiplied by vibrational wavefunctions. Nuclear motion has been modeled by vibrational packets moving on quadratic approximations to the potential energy surfaces for the X and A states of the cation. The overlap of the vibrational wavepackets decays in about 10-15 fs. Consequently, the oscillations in the strong field ionization decay on the same time scale. A revival of the vibrational overlap and in the oscillations of the strong field ionization is seen at 60-110 fs. TDCI simulations show that the decay and revival of the charge migration can be monitored by strong field ionization with intense 2- and 4-cycle linearly polarized 800 nm pulses. The revival is also seen with 7-cycle pulses.

Approximation Schemes to Include Nuclear Motion in Laser-Driven Ab Initio Electron Dynamics: Application to High Harmonic Generation.

The quantum mechanical description of many-electron dynamics in molecules driven by short laser pulses is at the heart of theoretical attochemistry. In addition to the formidable time-dependent electronic structure problem, the field faces the challenge that nuclear motion, ideally also treated quantum mechanically, may not be negligible, but scales enormously in effort. As a consequence, most first-principles calculations on ultrafast electron dynamics in molecules are done within the fixed-nuclear approximation. For laser-pulse excitation in H2+, where an "exact" treatment of the coupled nuclear-electron dynamics is possible, it has been shown that nuclear motion can have a nonnegligible impact on high harmonic generation (HHG) spectra (Witzorky et al., J. Chem. Theor. Comput. 2021, 17, 7353-7365). It is not so clear, however, how to include (quantum) nuclear motion also for more complicated molecules, with more electrons and/or nuclei, in particular when the electronic structure is described by correlated, multistate wavefunction methods such as the time-dependent configuration interaction (TD-CI). In this work, we suggest a scheme in which the Born-Oppenheimer potential energy surfaces of a molecule are approximated by model potentials (harmonic and asymptotic, as an expansion in 1/R), obtained from only a few ab initio calculations, with the prospect to treat complex molecular systems. The method is tested successfully for HHG by few-cycle laser pulses for the "exact" H2+ reference. It is then applied for diatomic molecules with more electrons and for a two-dimensional model of the water molecule using TD-CIS (S = single) for the electronic structure part.

Energy natural orbital characterization of nonadiabatic electron wavepackets in the densely quasi-degenerate electronic state manifold.

Dynamics and energetic structure of largely fluctuating nonadiabatic electron wavepackets are studied in terms of Energy Natural Orbitals (ENOs) [K. Takatsuka and Y. Arasaki, J. Chem. Phys. 154, 094103 (2021)]. Such huge fluctuating states are sampled from the highly excited states of clusters of 12 boron atoms (B12), which have densely quasi-degenerate electronic excited-state manifold, each adiabatic state of which gets promptly mixed with other states through the frequent and enduring nonadiabatic interactions within the manifold. Yet, the wavepacket states are expected to be of very long lifetimes. This excited-state electronic wavepacket dynamics is extremely interesting but very hard to analyze since they are usually represented in large time-dependent configuration interaction wavefunctions and/or in some other complicated forms. We have found that ENO gives an invariant energy orbital picture to characterize not only the static highly correlated electronic wavefunctions but also those time-dependent electronic wavefunctions. Hence, we first demonstrate how the ENO representation works for some general cases, choosing proton transfer in water dimer and electron-deficient multicenter chemical bonding in diborane in the ground state. We then penetrate with ENO deep into the analysis of the essential nature of nonadiabatic electron wavepacket dynamics in the excited states and show the mechanism of the coexistence of huge electronic fluctuation and rather strong chemical bonds under very random electron flows within the molecule. To quantify the intra-molecular energy flow associated with the huge electronic-state fluctuation, we define and numerically demonstrate what we call the electronic energy flux.

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.

Strong-field molecular ionization beyond the single active electron approximation.

This work explores quantitative limits to the single-active electron approximation, often used to deal with strong-field ionization and subsequent attosecond dynamics. Using a time-dependent, multiconfiguration approach, specifically the time-dependent configuration interaction method, we solve the time-dependent Schrödinger equation for the two-electron dihydrogen molecule with the possibility of tuning at will the electron-electron interaction by an adiabatic switch-on/switch-off function. We focus on signals of the single ionization of H2 under a strong near-infrared, four-cycle, linearly polarized laser pulse of varying intensity and within a vibrationally frozen molecular model. The observables we address are post-pulse total ionization probability profiles as a function of the laser peak intensity. Three values of the internuclear distance R taken as a parameter are considered, R = Req = 1.4 a.u. for the equilibrium geometry of the molecule, R = 5.0 a.u. for an elongated molecule, and R = 10.2 a.u. for a dissociating molecule. The most striking observation is the non-monotonous behavior of the ionization probability profiles at intermediate elongation distances with an instance of enhanced ionization and one of partial ionization quenching. We give an interpretation of this in terms of a resonance-enhanced-multiphoton ionization mechanism with interfering overlapping resonances resulting from excited electronic states.

A TD-CIS study of high-harmonic generation of uracil cation fragments

High-harmonic generation (HHG) spectroscopy was investigated from laser ablation plasmas of uracil targets in the experimental work of Hutchinson et al., [Phys. Chem. Chem. Phys. Comparison of high-order harmonic generation in uracil and thymine ablation plumes, 2013, 15, 12308] and of Lopez-Quintas et al., [Appl. Surf. Sci., Characterization of laser-induced plasmas of nucleobases: Uracil and thymine, 2014, 302, 299]. Inspired by these works we computed HHG spectra of uracil cation fragments by means of real-time time-dependent configuration interaction with single excitations (RT-TD-CIS) and continuum-optimal Gaussian basis sets. We calculated the rotationally-averaged HHG spectra of different cation fragments analysing the relation between HHG, geometry of the fragments (linear/nonlinear) and electron dynamics (bound and continuum energy states). Moreover, to simulate the complexity of a sample containing fragments upon ablation, we calculated weighted-rotationally-averaged HHG spectra where the contribution of each fragment is included through weights assigned accordingly to the experimental mass spectra.

Electronic Dynamics of a Molecular System Coupled to a Plasmonic Nanoparticle Combining the Polarizable Continuum Model and Many-Body Perturbation Theory.

The efficiency of plasmonic metallic nanoparticles in harvesting and concentrating light energy in their proximity triggers a wealth of important and intriguing phenomena. For example, spectroscopies are able to reach single-molecule and intramolecule sensitivities, and important chemical reactions can be effectively photocatalyzed. For the real-time description of the coupled dynamics of a molecule’s electronic system and of a plasmonic nanoparticle, a methodology has been recently proposed (J. Phys. Chem. C. 120, 2016, 28774−28781) which combines the classical description of the nanoparticle as a polarizable continuum medium with a quantum-mechanical description of the molecule treated at the time-dependent configuration interaction (TDCI) level. In this work, we extend this methodology by describing the molecule using many-body perturbation theory: the molecule’s excitation energies, transition dipoles, and potentials computed at the GW/Bethe–Salpeter equation (BSE) level. This allows us to overcome current limitations of TDCI in terms of achievable accuracy without compromising on the accessible molecular sizes. We illustrate the developed scheme by characterizing the coupled nanoparticle/molecule dynamics of two prototype molecules, LiCN and p-nitroaniline.

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 high harmonic generation in the helium atom subjected to laser pulses of extremely high intensity. To do this, we perform real-time propagations of helium atom wavefunction using quantum chemistry methods coupled to Gaussian basis sets. Calculations are performed within the real-time time-dependent configuration interaction framework at two levels of theory: time-dependent configuration interaction with single excitations (uncorrelated method) and time-dependent full configuration interaction (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 time-dependent configuration interaction 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 high harmonic generation (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.

Making optical excitations visible – An exciton wavefunction extension to the time-dependent configuration interaction method

We report the animation and analysis of laser-driven many-electron dynamics by means of time-dependent local densities and quantities derived from the one-particle transition density matrix, such as time-dependent natural transition orbital densities, particle and hole positions and exciton size. The time-dependent configuration interaction method was used to revisit studies on hydrogen molecule and lithium cyanide by Saalfrank and coworkers and our study shines new light on the optical transitions in these small molecules to benchmark the implemented time-dependent exciton properties. One focus of our simulations is the comparison of the local densities to the quantities of the two-body exciton wavefunction for resonant excitations leading either to state-to-state transitions or to the creation of a wave packet.

Ionization of HCCI Neutral and Cations by Strong Laser Fields Simulated With Time Dependent Configuration Interaction.

Strong field ionization of neutral iodoacetylene (HCCI) can produce a coherent superposition of the X and A cations. This superposition results in charge migration between the CC π orbital and the iodine π-type lone pair which can be monitored by strong field ionization with short, intense probe pulses. Strong field ionization of the X and A states of HCCI cation was simulated with time-dependent configuration interaction using singly ionized configurations and singly excited, singly ionized configurations (TD-CISD-IP) and an absorbing boundary. Studies with static fields were used to obtain the 3-dimensional angular dependence of instantaneous ionization rates by strong fields and the orbitals involved in producing the cations and dications. The frequency of charge oscillation is determined by the energy separation of the X and A states; this separation can change depending on the direction and strength of the field. Furthermore, fields along the molecular axis can cause extensive mixing between the field-free X and A configurations. For coherent superpositions of the X and A states, the charge oscillations are characterized by two frequencies–the driving frequency of the laser field of the probe pulse and the intrinsic frequency due to the energy separation between the X and A states. For linear and circularly polarized pulses, the ionization rates show marked differences that depend on the polarization direction of the pulse, the carrier envelope phase and initial phase of the superposition. Varying the initial phase of the superposition at the beginning of the probe pulse is analogous to changing the delay between the pump and probe pulses. The charge oscillation in the coherent superposition of the X and A states results in maxima and minima in the ionization yield as a function of the superposition phase.

Floquet Time-Dependent Configuration Interaction for Modeling Ultrafast Electron Dynamics

Time-dependent electronic structure methods are a valuable tool for simulating spectroscopic experiments. Recent advances in time-dependent configuration interaction (TDCI) algorithms have made them an attractive means of modeling many-electron dynamics, particularly for cases where multireference effects are essential. Here we present an extension to TDCI, Floquet TDCI, where the electronic wave function is expanded in a basis of light-dressed determinants. Our approach is based on our high-performance graphics processing unit (GPU) accelerated implementation of complete active space configuration interaction (CASCI). Simulations of two-photon absorption demonstrate that Floquet TDCI is well-suited for modeling dynamics in intense, ultrashort laser pulses. Accurate results are obtained for pulse energies up to ∼4 × 10-4 J/cm2 per pulse in the most difficult case explored here. By simulation of a set of molecules under continuous wave coupling, we demonstrate the ability of Floquet to describe the entanglement of light and multiple molecules in a cavity (i.e., a cavity polariton). Excellent computational performance is observed: a 320 fs propagation of a large dye (C30N2H22) with a 2 as timestep and a large active space (10 electrons in 11 orbitals), including a monochromatic pulse with three photon states, was performed in 3 h 6 min on a single Tesla V100 GPU. Our Floquet TDCI algorithm scales linearly with the number of photon states and exponentially with the number of photon colors included in the calculation. We argue that its energy-conserving nature makes Floquet TDCI well-suited to drive nonadiabatic molecular dynamics simulations.

Relativistic time-dependent configuration-interaction singles method

In this work, a derivation and implementation of the relativistic time-dependent configuration interaction singles (RTDCIS) method is presented. Various observables for krypton and xenon atoms obtained by RTDCIS are compared with experimental data and alternative relativistic calculations. This includes energies of occupied orbitals in the Dirac-Fock ground state, Rydberg state energies, Fano resonances and photoionization cross sections. Diagrammatic many-body perturbation theory, based on the relativistic random phase approximation, is used as a benchmark with excellent agreement between RTDCIS reported at the Tamm-Dancoff level. Results from RTDCIS are computed in the length gauge, here the negative energy states can be omitted with acceptable loss of accuracy. A complex absorbing potential, that is used to remove photoelectrons far from the ion, is implemented as a scalar potential and validated for RTDCIS. The RTDCIS methodology presented here opens for future studies of strong-field processes, such as attosecond transient absorption and high-order harmonic generation, with electron and hole spin dynamics and other relativistic effects described by first principle via the Dirac equation.

RhoDyn: A ρ-TD-RASCI Framework to Study Ultrafast Electron Dynamics in Molecules.

This article presents the program module RhoDyn as part of the OpenMOLCAS project intended to study ultrafast electron dynamics within the density-matrix-based time-dependent restricted active space configuration interaction framework (ρ-TD-RASCI). The formalism allows for the treatment of spin-orbit coupling effects, accounts for nuclear vibrations in the form of a vibrational heat bath, and naturally incorporates (auto)ionization effects. Apart from describing the theory behind and the program workflow, the paper also contains examples of its application to the simulations of the linear L2,3 absorption spectra of a titanium complex, high harmonic generation in the hydrogen molecule, ultrafast charge migration in benzene and iodoacetylene, and spin-flip dynamics in the core excited states of iron complexes.

Gaussian-Type Orbital Calculations for High Harmonic Generation in Vibrating Molecules: Benchmarks for H2.

The response of the hydrogen molecular ion, H2+, 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.