Understanding real-time time-dependent density-functional theory simulations of ultrafast laser-induced dynamics in organic molecules

Real-Time Simulation of Ultrafast Electronic Dynamics of Nanoscale Systems Involving an Organic Molecule and a Nanoparticle Dimer

Understanding and predicting the behavior of nanomaterials composed of plasmons interacting with quantum emitters at ultrafast timescales is crucial for the better manipulation of light at the nanoscale and advancing technologies like ultrafast communication and computing. Here we perform a simulation of the "real-time" electronic dynamics of a coupled molecule-metal nanoparticle dimer interacting with an ultrashort resonant laser pulse by combining the real-time time-dependent density functional theory (RT-TDDFT) approach with the time-domain frequency-dependent fluctuating charge (TD-ωFQ) model, an atomistic electromagnetic (AEM) model for the dynamic plasmonic response of nanoparticles. It is shown that the induced dipoles evolve from an exponential decay pattern to a beat pattern with an increase in coupling strength, which is altered by changing the molecular orientation relative to the dimer axis. It is further shown that in the strong coupling regime, both the excited molecule and the plasmon relax rapidly due to the molecule-plasmon interaction, and the efficient coherent energy exchange between the interacting molecule and plasmon modes occurs on a femtosecond (fs) timescale. This work provides guidance on manipulating light-matter interaction and studying molecular plasmonics at extremely fast timescales.

On-the-fly simulation of time-resolved fluorescence spectra and anisotropy

We combine on-the-fly trajectory surface hopping simulations and the doorway-window representation of nonlinear optical response functions to create an efficient protocol for the evaluation of time- and frequency-resolved fluorescence (TFRF) spectra and anisotropies of the realistic polyatomic systems. This approach gives the effective description of the proper (e.g., experimental) pulse envelopes, laser field polarizations, and the proper orientational averaging of TFRF signals directly from the well-established on-the-fly nonadiabatic dynamic simulations without extra computational cost. To discuss the implementation details of the developed protocol, we chose cis-azobenzene as a prototype to simulate the time evolution of the TFRF spectra governed by its nonadiabatic dynamics. The results show that the TFRF is determined by the interplay of several key factors, i.e., decays of excited-state populations, evolution of the transition dipole moments along with the dynamic propagation, and scaling factor of the TFRF signals associated with the cube of emission frequency. This work not only provides an efficient and effective approach to simulate the TFRF and anisotropies of realistic polyatomic systems but also discusses the important relationship between the TFRF signals and the underlining nonadiabatic dynamics.

Pulse-Induced Dynamics of a Charge-Transfer Complex from First Principles

The ultrafast dynamics of charge carriers in organic donor-acceptor interfaces are of primary importance to understanding the fundamental properties of these systems. In this work, we focus on a charge-transfer complex formed by quaterthiophene p-doped by tetrafluoro-tetracyanoquinodimethane and investigate electron dynamics and vibronic interactions also at finite temperatures by applying a femtosecond pulse in resonance with the two lowest energy excitations of the system with perpendicular and parallel polarization with respect to the interface. The adopted ab initio formalism based on real-time time-dependent density-functional theory coupled to Ehrenfest dynamics enables monitoring the dynamical charge transfer across the interface and assessing the role played by the nuclear motion. Our results show that the strong intermolecular interactions binding the complex already in the ground state influence the dynamics, too. The analysis of the nuclear motion involved in these processes reveals the participation of different vibrational modes depending on the electronic states stimulated by the resonant pulse. Coupled donor-acceptor modes mostly influence the excited state polarized across the interface, while intramolecular vibrations in the donor molecule dominate the excitation in the orthogonal direction. The results obtained at finite temperatures are overall consistent with this picture, although thermal disorder contributes to slightly decreasing interfacial charge transfer.

Size-dependent errors in real-time electron density propagation

Real-time (RT) electron density propagation with time-dependent density functional theory (TDDFT) or Hartree-Fock (TDHF) is one of the most popular methods to model the charge transfer in molecules and materials. However, both RT-TDHF and RT-TDDFT within the adiabatic approximation are known to produce inaccurate evolution of the electron density away from the ground state in model systems, leading to large errors in charge transfer and erroneous shifting of peaks in absorption spectra. Given the poor performance of these methods with small model systems and the widespread use of the methods with larger molecular and material systems, here we bridge the gap in our understanding of these methods and examine the size-dependence of errors in RT density propagation. We analyze the performance of RT density propagation for systems of increasing size during the application of a continuous resonant field to induce Rabi-like oscillations, during charge-transfer dynamics, and for peak shifting in simulated absorption spectra. We find that the errors in the electron dynamics are indeed size dependent for these phenomena, with the largest system producing the results most aligned with those expected from linear response theory. The results suggest that although the RT-TDHF and RT-TDDFT methods may produce severe errors for model systems, the errors in charge transfer and resonantly driven electron dynamics may be much less significant for more realistic, large-scale molecules and materials.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Analyzing Excitation-Energy Transfer Based on the Time-Dependent Density Functional Theory in Real Time

Excitation-energy transfer is a key step in processes such as photosynthesis that convert light into other forms of energy. Time-dependent density functional theory (DFT) in real time is ideal for the first-principles simulation of such processes due to its computational efficiency. We here demonstrate how real-time DFT can be used for analyzing excitation-energy transfer from first-principles. We discuss several measures of energy transfer that are based solely on the time-dependent density, are well founded in the DFT framework, allow for intuitive understanding and visualization, and reproduce important limiting cases of an analytical model. We demonstrate their usefulness in calculations for model systems, both with static nuclei and in the context of DFT-based Ehrenfest dynamics.

Ultrafast Internal Conversion Dynamics through the on-the-Fly Simulation of Transient Absorption Pump-Probe Spectra with Different Electronic Structure Methods

An on-the-fly surface-hopping simulation protocol is developed for the evaluation of transient absorption (TA) pump-probe (PP) signals of molecular systems exhibiting internal conversion to the electronic ground state. We study the nonadiabatic dynamics of azomethane and the associating TA PP spectra at three levels of the electronic-structure theory, OM2/MRCI, SA-CASSCF, and XMS-CASPT2. The impact of these methods on the population dynamics and time-resolved TA PP signals is substantially different. This difference is attributed to the strong non-Condon effects that must be taken into account for the proper understanding and interpretation of time-resolved TA PP signals of nonadiabatic polyatomic systems. This shows that the combination of the dynamical and spectral simulations definitely provides more accurate and detailed information on the microscopic mechanisms of photophysical and photochemical processes. Hence the simulation of time-resolved spectroscopic signals provides another important dimension to examine the accuracy of quantum chemistry methods.

Understanding the mechanism of plasmon-driven water splitting: hot electron injection and a near field enhancement effect

Utilizing plasmon-generated hot carriers to drive chemical reactions has currently become an active area of research in solar photocatalysis at the nanoscale. However, the mechanism underlying exact transfer and the generation dynamics of hot carriers, and the strategies used to further improve the quantum efficiency of the photocatalytic reaction still deserve further investigation. In this work, we perform a nonadiabatic excited-state dynamics study to depict the correlation between the reaction rate of plasmon-driven water splitting (PDWS) and the sizes of gold particles, the incident light frequency and intensity, and the near-field spatial distribution. Four model systems, H₂O and Au₂₀@H₂O separately interacting with the laser field and the near field generated by the Au nanoparticle (NP) with a few nanometers in size, have been investigated. Our simulated results clearly unveil the mechanism of PDWS and hot-electron injection in a Schottky-free junction: the electrons populated on the antibonding orbitals of H₂O are mandatory to drive the OH bond breaking and the strong orbital hybridization between Au₂₀ and H₂O creates the conditions for direct electron injection. We further find that the linear dependence of the reaction rate and the field amplitude only holds at a relatively weak field and it breaks down when the second OH bond begins to dissociate and field-induced water fragmentation occurs at a very intensive field, and that with the guarantee of electron injection, the water splitting rate increases with an increase in the NP size. This study will be helpful for further improving the efficiency of photochemical reactions involving plasmon-generated hot carriers and expanding the applications of hot carriers in a variety of chemical reactions.

Laser-Induced Electronic and Vibronic Dynamics in the Pyrene Molecule and Its Cation

Among polycyclic aromatic hydrocarbons, pyrene is widely used as an optical probe thanks to its peculiar ultraviolet absorption and infrared emission features. Interestingly, this molecule is also an abundant component of the interstellar medium, where it is detected via its unique spectral fingerprints. In this work, we present a comprehensive first-principles study on the electronic and vibrational response of pyrene and its cation to ultrafast, coherent pulses in resonance with their optically active excitations in the ultraviolet region. The analysis of molecular symmetries, electronic structure, and linear optical spectra is used to interpret transient absorption spectra and kinetic energy spectral densities computed for the systems excited by ultrashort laser fields. By disentangling the effects of the electronic and vibrational dynamics via ad hoc simulations with stationary and moving ions, and, in specific cases, with the aid of auxiliary model systems, we rationalize that the nuclear motion is mainly harmonic in the neutral species, while strong anharmonic oscillations emerge in the cation, driven by electronic coherence. Our results provide additional insights into the ultrafast vibronic dynamics of pyrene and related compounds and set the stage for future investigations on more complex carbon-conjugated molecules.

Spectral Fingerprint of Excited-State Energy Transfer in Dendrimers through Polarization-Sensitive Transient-Absorption Pump-Probe Signals: On-the-Fly Nonadiabatic Dynamics Simulations

The time-resolved polarization-sensitive transient-absorption (TA) pump-probe (PP) spectra are simulated using on-the-fly surface-hopping nonadiabatic dynamics and the doorway-window representation of nonlinear spectroscopy. A dendrimer model system composed of two linear phenylene ethynylene units (2-ring and 3-ring) is taken as an example. The ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA) contributions as well as the total TA PP signals are obtained and carefully analyzed. It is shown that intramolecular excited-state energy transfer from the 2-ring unit to the 3-ring unit can be conveniently identified by employing pump and probe pulses with different polarizations. Our results demonstrate that time-resolved polarization-sensitive TA PP signals provide a powerful tool for the elucidation of excited-state energy-transfer pathways, notably in molecular systems possessing several optically bright nonadiabatically coupled electronic states with different orientations of transition dipole moments.

LayerPCM: An implicit scheme for dielectric screening from layered substrates

We present LayerPCM, an extension of the polarizable-continuum model coupled to real-time time-dependent density-functional theory, for an efficient and accurate description of the electrostatic interactions between molecules and multilayered dielectric substrates on which they are physisorbed. The former are modeled quantum-mechanically, while the latter are treated as polarizable continua characterized by their dielectric constants. The proposed approach is purposely designed to simulate complex hybrid heterostructures with nano-engineered substrates including a stack of anisotropic layers. LayerPCM is suitable for describing the polarization-induced renormalization of frontier energy levels of the adsorbates in the static regime. Moreover, it can be reliably applied to simulating laser-induced ultrafast dynamics of molecules through the inclusion of electric fields generated by Fresnel-reflection at the substrate. Depending on the complexity of the underlying layer structure, such reflected fields can assume non-trivial shapes and profoundly affect the dynamics of the photo-excited charge carriers in the molecule. In particular, the interaction with the substrate can give rise to strong delayed fields, which lead to interference effects resembling those of multi-pulse-based spectroscopy. The robustness of the implementation and the above-mentioned features are demonstrated with a number of examples, ranging from intuitive models to realistic systems.

Nonlinear light absorption in many-electron systems excited by an instantaneous electric field: a non-perturbative approach

Applications of low-cost non-perturbative approaches in real time, such as time-dependent density functional theory, for the study of nonlinear optical properties of large and complex systems are gaining increasing popularity. However, their assessment still requires the analysis and understanding of elementary dynamical processes in simple model systems. Motivated by the aim of simulating optical nonlinearities in molecules, here exemplified by the case of the quaterthiophene oligomer, we investigate light absorption in many-electron interacting systems beyond the linear regime by using a single broadband impulse of an electric field; i.e. an electrical impulse in the instantaneous limit. We determine non-pertubatively the absorption cross section from the Fourier transform of the time-dependent induced dipole moment, which can be obtained from the time evolution of the wavefunction. We discuss the dependence of the resulting cross section on the magnitude of the impulse and we highlight the advantages of this method in comparison with perturbation theory by working on a one-dimensional model system for which numerically exact solutions are accessible. Thus, we demonstrate that the considered non-pertubative approach provides us with an effective tool for investigating fluence-dependent nonlinear optical excitations.

Simulating Vibronic Spectra without Born-Oppenheimer Surfaces

We show how linear vibronic spectra in molecular systems can be simulated efficiently using first-principles approaches without relying on the explicit use of multiple Born-Oppenheimer potential energy surfaces. We demonstrate and analyze the performance of mean-field and beyond-mean-field dynamics techniques for the H₂ molecule in one dimension, in the later case capturing the vibronic structure quite accurately, including quantum Franck-Condon effects. In a practical application of this methodology we simulate the absorption spectrum of benzene in full dimensionality using time-dependent density functional theory at the multitrajectory Ehrenfest level, finding good qualitative agreement with experiment and significant spectral reweighting compared to commonly used single-trajectory Ehrenfest dynamics. These results form the foundation for nonlinear spectral calculations and show promise for future application in capturing phenomena associated with vibronic coupling in more complex molecular and potentially condensed phase systems.

Simulation of photo-electron spectrum and electron scattering by dual time propagation

A dual time propagation approach is introduced to describe electron scattering and ionization. The space is divided into two regions, a central region with a full time-dependent Hamiltonian and an outer region where the kinetic operator and the laser field dominate. The two regions are connected by a source term. Time-dependent density functional theory calculations of wave packet scattering on molecules and photoelectron spectrum due to circularly polarized laser are presented to illustrate the efficiency and applicability of the approach.

Exploring organic semiconductors in solution: the effects of solvation, alkylization, and doping

The first-principles simulation of the electronic structure of organic semiconductors in solution poses a number of challenges that are not trivial to address simultaneously. In this work, we investigate the effects and the mutual interplay of solvation, alkylization, and doping on the structural, electronic, and optical properties of sexithiophene, a representative organic semiconductor molecule. To this end, we employ (time-dependent) density functional theory in conjunction with the polarizable-continuum model. We find that the torsion between adjacent monomer units plays a key role, as it strongly influences the electronic structure of the molecule, including energy gap, ionization potential, and band widths. Alkylization promotes delocalization of the molecular orbitals up to the first methyl unit, regardless of the chain length, leading to an overall shift of the energy levels. The alterations in the electronic structure are reflected in the optical absorption, which is additionally affected by dynamical solute-solvent interactions. Taking all these effects into account, solvents decrease the optical gap by an amount that depends on its polarity, and concomitantly increase the oscillator strength of the first excitation. The interaction with a dopant molecule promotes planarization. In such scenario, solvation and alkylization enhance charge transfer both in the ground state and in the excited state.