Modeling Fast Electron Dynamics with Real-Time Time-Dependent Density Functional Theory: Application to Small Molecules and Chromophores

Assessment of Real-Time Time-Dependent Density Functional Theory (RT-TDDFT) in Radiation Chemistry: Ionized Water Dimer

Ionization in the condensed phase and molecular clusters leads to a complicated chain of processes with coupled electron-nuclear dynamics. It is difficult to describe such dynamics with conventional nonadiabatic molecular dynamics schemes since the number of states swiftly increases as the molecular system grows. It is therefore attractive to use a direct electron and nuclear propagation such as the real-time time-dependent density functional theory (RT-TDDFT). Here we report a RT-TDDFT benchmark study on simulations of singly and doubly ionized states of a water monomer and dimer as a prototype for more complex processes in a condensed phase. We employed the RT-TDDFT based Ehrenfest molecular dynamics with a generalized gradient approximate (GGA) functional and compared it with wave-function-based surface hopping (SH) simulations. We found that the initial dynamics of a singly HOMO ionized water dimer is similar for both the RT-TDDFT/GGA and the SH simulations but leads to completely different reaction channels on a longer time scale. This failure is attributed to the self-interaction error in the GGA functionals and it can be avoided by using hybrid functionals with large fraction of exact exchange (represented here by the BHandHLYP functional). The simulations of doubly ionized states are reasonably described already at the GGA level. This suggests that the RT-TDDFT/GGA method could describe processes following the autoionization processes such as Auger emission, while its applicability to more complex processes such as intermolecular Coulombic decay remains limited.

Chemical Reactivity and Spectroscopy Explored From QM/MM Molecular Dynamics Simulations Using the LIO Code

In this work we present the current advances in the development and the applications of LIO, a lab-made code designed for density functional theory calculations in graphical processing units (GPU), that can be coupled with different classical molecular dynamics engines. This code has been thoroughly optimized to perform efficient molecular dynamics simulations at the QM/MM DFT level, allowing for an exhaustive sampling of the configurational space. Selected examples are presented for the description of chemical reactivity in terms of free energy profiles, and also for the computation of optical properties, such as vibrational and electronic spectra in solvent and protein environments.

Phase Cycling RT-TDDFT Simulation Protocol for Nonlinear XUV and X-ray Molecular Spectroscopy

Real-time time-dependent density functional theory (RT-TDDFT) provides a practical algorithm for propagating a many-electron system driven by external laser fields. The fields are included nonperturbatively in the propagation, and the molecular reduced single-electron density operator and various spectroscopic and diffraction signals can be computed directly, avoiding the expensive calculation of many-body states. Nonlinear optical signals contain contributions of multiple pathways. A phase cycling protocol is implemented in order to separate these pathways. Simulations of XUV four-wave mixing signals in the CO molecule are compared with ab initio sum-over-states calculations.

Excited-state absorption in tetrapyridyl porphyrins: comparing real-time and quadratic-response time-dependent density functional theory

Three meso-substituted tetrapyridyl porphyrins (free base, Ni(ii), and Cu(ii)) were investigated for their optical limiting (OL) capabilities using real-time (RT-), linear-response (LR-), and quadratic-response (QR-) time-dependent density functional theory (TDDFT) methods. These species are experimentally known to display a prominent reverse saturable absorption feature between the Q and B bands of the ground-state absorption (GSA), which has been attributed to increased excited-state absorption (ESA) relative to GSA. A recently developed RT-TDDFT based method for calculating ESA from a LR-TDDFT density was utilized with eight exchange-correlation functionals (BLYP, PBE, B3LYP, CAM-B3LYP, PBE0, M06, BHLYP, and BHandH) and contrasted with calculations of ESA using QR-TDDFT with five exchange-correlation functionals (BLYP, B3LYP, CAM-B3LYP, BHLYP, and BHandH). This allowed for comparison between functionals with varying amounts of exact exchange as well as between the ability of RT-TDDFT and QR-TDDFT to reproduce OL behavior in porphyrin systems. The absorption peak positions and intensities for GSA and ESA are significantly impacted by the choice of DFT functional, with the most critical factor identified as the amount of exact exchange in the functional form. Calculating ESA with QR-TDDFT is found to be significantly more sensitive to the amount of exact exchange than GSA and ESA with RT-TDDFT, as well as GSA with LR-TDDFT. An analogous behavior is also demonstrated for the polycyclic aromatic hydrocarbon coronene. This is problematic when using the same approximate functional for calculation of both GSA and ESA, as the LR- and QR-TDDFT excitation energies will not have similar errors. Overall, the RT-TDDFT method with hybrid functionals reproduces the OL features for the porphyrin systems studied here and is a viable computational approach for efficient screening of molecular complexes for OL properties.

Candidate eco-friendly gas mixtures for MPGDs

Modern gas detectors for detection of particles require F-based gases for optimal performance. Recent regulations demand the use of environmentally unfriendly F-based gases to be limited or banned. This review studies properties of potential eco-friendly gas candidate replacements.

Charge transfer in time-dependent density functional theory

Charge transfer plays a crucial role in many processes of interest in physics, chemistry, and bio-chemistry. In many applications the size of the systems involved calls for time-dependent density functional theory (TDDFT) to be used in their computational modeling, due to its unprecedented balance between accuracy and efficiency. However, although exact in principle, in practise approximations must be made for the exchange-correlation functional in this theory, and the standard functional approximations perform poorly for excitations which have a long-range charge-transfer component. Intense progress has been made in developing more sophisticated functionals for this problem, which we review. We point out an essential difference between the properties of the exchange-correlation kernel needed for an accurate description of charge-transfer between open-shell fragments and between closed-shell fragments. We then turn to charge-transfer dynamics, which, in contrast to the excitation problem, is a highly non-equilibrium, non-perturbative, process involving a transfer of one full electron in space. This turns out to be a much more challenging problem for TDDFT functionals. We describe dynamical step and peak features in the exact functional evolving over time, that are missing in the functionals currently used. The latter underestimate the amount of charge transferred and manifest a spurious shift in the charge transfer resonance position. We discuss some explicit examples.

Attosecond Charge Migration with TDDFT: Accurate Dynamics from a Well-Defined Initial State

We investigate the ability of time-dependent density functional theory (TDDFT) to capture attosecond valence electron dynamics resulting from sudden X-ray ionization of a core electron. In this special case the initial state can be constructed unambiguously, allowing for a simple test of the accuracy of the dynamics. The response following nitrogen K-edge ionization in nitrosobenzene shows excellent agreement with fourth-order algebraic diagrammatic construction (ADC(4)) results, suggesting that a properly chosen initial state allows TDDFT to adequately capture attosecond charge migration. Visualizing hole motion using an electron localization picture (ELF), we provide an intuitive chemical interpretation of the charge migration as a superposition of Lewis dot resonance structures.

X-ray circular dichroism signals: a unique probe of local molecular chirality

Core-resonant circular dichroism (CD) signals are induced by molecular chirality and vanish for achiral molecules and racemic mixtures. The highly localized nature of core excitations makes them ideal probes of local chirality within molecules. Simulations of the circular dichroism spectra of several molecular families illustrate how these signals vary with the electronic coupling to substitution groups, the distance between the X-ray chromophore and the chiral center, geometry, and chemical structure. Clear insight into the molecular structure is obtained through analysis of the X-ray CD spectra.

Probing Zeolite Crystal Architecture and Structural Imperfections using Differently Sized Fluorescent Organic Probe Molecules

A micro‐spectroscopic method has been developed to probe the accessibility of zeolite crystals using a series of fluorescent 4‐(4‐diethylaminostyryl)‐1‐methylpyridinium iodide (DAMPI) probes of increasing molecular size. Staining large zeolite crystals with MFI (ZSM‐5) topology and subsequent mapping of the resulting fluorescence using confocal fluorescence microscopy reveal differences in structural integrity: the 90° intergrowth sections of MFI crystals are prone to develop structural imperfections, which act as entrance routes for the probes into the zeolite crystal. Polarization‐dependent measurements provide evidence for the probe molecule's alignment within the MFI zeolite pore system. The developed method was extended to BEA (Beta) crystals, showing that the previously observed hourglass pattern is a general feature of BEA crystals with this morphology. Furthermore, the probes can accurately identify at which crystal faces of BEA straight or sinusoidal pores open to the surface. The results show this method can spatially resolve the architecture‐dependent internal pore structure of microporous materials, which is difficult to assess using other characterization techniques such as X‐ray diffraction.

Ultra-fast charge migration competes with proton transfer in the early chemistry of H2O˙+

Oxidation by the ultra-short lived radical cation of water, H₂O˙⁺, can potentially take place at the interface of water and numerous heterogeneous systems involved in radiation therapy, energy and environmental industries.

Can Quantized Vibrational Effects Be Obtained from Ehrenfest Mixed Quantum-Classical Dynamics?

We explore the question of whether mean-field or "Ehrenfest" mixed quantum-classical dynamics is capable of capturing the quantized vibrational features in photoabsorption spectra that result from infrared and Raman-active vibrational transitions. We show that vibrational and electronic absorption spectra can indeed be obtained together within a single Ehrenfest simulation. Furthermore, the electronic transitions show new sidebands that are absent in electronic dynamics simulations with fixed nuclei. Inspection of the electronic sidebands reveals that the spacing corresponds to vibrational frequencies of totally symmetric vibrational modes of the ground electronic state. A simple derivation of the time-evolving dipole in the presence of external fields and vibrational motion shows the origin of these features, demonstrating that mixed quantum-classical Ehrenfest dynamics is capable of producing infrared, Raman, and electronic absorption spectra from a single simulation.

Wire and extended ladder model predict THz oscillations in DNA monomers, dimers and trimers

We call monomer a B-DNA base pair and study, analytically and numerically, electron or hole oscillations in monomers, dimers and trimers. We employ two tight binding (TB) approaches: (I) at the base-pair level, using the on-site energies of the base pairs and the hopping parameters between successive base pairs i.e. a wire model, and (II) at the single-base level, using the on-site energies of the bases and the hopping parameters between neighbouring bases, specifically between (a) two successive bases in the same strand, (b) complementary bases that define a base pair, and (c) diagonally located bases of successive base pairs, i.e. an extended ladder model since it also includes the diagonal hoppings (c). For monomers, with TB II, we predict periodic carrier oscillations with frequency [Formula: see text]-550 THz. For dimers, with TB I, we predict periodic carrier oscillations with [Formula: see text]-100 THz. For trimers made of identical monomers, with TB I, we predict periodic carrier oscillations with [Formula: see text]-33 THz. In other cases, either with TB I or TB II, the oscillations may be not strictly periodic, but Fourier analysis shows similar frequency content. For dimers and trimers, TB I and TB II are successfully compared giving complementary aspects of the oscillations.

Employing Range Separation on the meta-GGA Rung: New Functional Suitable for Both Covalent and Noncovalent Interactions

We devise a scheme for converting an existing exchange functional into its range-separated hybrid variant. The underlying exchange hole of the Becke-Roussel type has the exact second-order expansion in the interelectron distance. The short-range part of the resulting range-separated exchange energy depends on the kinetic energy density and the Laplacian even if the base functional lacks the dependence on these variables. The most successful practical realization of the scheme, named LC-PBETPSS, combines the range-separated Perdew-Burke-Ernzerhof (PBE) exchange lifted to the hybrid meta-generalized gradient approximation rung and the Tao-Perdew-Staroverov-Scuseria (TPSS) correlation. The value of the range-separation parameter is estimated theoretically and confirmed by empirical optimization. The D3 dispersion correction is recommended for all energy computations employing the presented functional. Numerical tests show remarkably robust performance of the method for noncovalent interaction energies, barrier heights, main-group thermochemistry, and excitation energies.

Excited-State Absorption from Real-Time Time-Dependent Density Functional Theory: Optical Limiting in Zinc Phthalocyanine

Optical-limiting materials are capable of attenuating light to protect delicate equipment from high-intensity light sources. Phthalocyanines have attracted a lot of attention for optical-limiting applications due to their versatility and large nonlinear absorption. With excited-state absorption (ESA) being the primary mechanism for optical limiting behavior in phthalocyanines, the ability to tune the optical absorption of ground and excited states in phthalocyanines would allow for the development of advanced optical limiters. We recently developed a method for the calculation of ESA based on real-time time-dependent density functional theory propagation of an excited-state density. In this work, we apply the approach to zinc phthalocyanine, demonstrating the ability of our method to efficiently identify the optical limiting potential of a molecular complex.

Efficient Algorithms for Estimating the Absorption Spectrum within Linear Response TDDFT

We present a special symmetric Lanczos algorithm and a kernel polynomial method (KPM) for approximating the absorption spectrum of molecules within the linear response time-dependent density functional theory (TDDFT) framework in the product form. In contrast to existing algorithms, the new algorithms are based on reformulating the original non-Hermitian eigenvalue problem as a product eigenvalue problem and the observation that the product eigenvalue problem is self-adjoint with respect to an appropriately chosen inner product. This allows a simple symmetric Lanczos algorithm to be used to compute the desired absorption spectrum. The use of a symmetric Lanczos algorithm only requires half of the memory compared with the nonsymmetric variant of the Lanczos algorithm. The symmetric Lanczos algorithm is also numerically more stable than the nonsymmetric version. The KPM algorithm is also presented as a low-memory alternative to the Lanczos approach, but the algorithm may require more matrix-vector multiplications in practice. We discuss the pros and cons of these methods in terms of their accuracy as well as their computational and storage cost. Applications to a set of small and medium-sized molecules are also presented.

Interplay between plasmon and single-particle excitations in a metal nanocluster

Plasmon-generated hot carriers are used in photovoltaic or photochemical applications. However, the interplays between the plasmon and single-particle excitations in nanosystems have not been theoretically addressed using ab initio methods. Here we show such interplays in a Ag₅₅ nanocluster using real-time time-dependent density functional theory simulations. We find that the disappearance of the zero-frequency peak in the Fourier transform of the band-to-band transition coefficient is a hallmark of the plasmon. We show the importance of the d-states for hot-carrier generations. If the single-particle d-to-s excitations are resonant to the plasmon frequency, the majority of the plasmon energy will be converted into hot carriers, and the overall hot-carrier generation is enhanced by the plasmon; if such resonance does not exist, we observe an intriguing Rabi oscillation between the plasmon and hot carriers. Phonons play a minor role in plasmonic dynamics in such small systems. This study provides guidance on improving plasmonic applications.

Plasmons can enhance hot-carrier generation for efficient photochemical reactions, but the interplay between plasmons and single-particle excitations are difficult to capture in models. Here, the authors use real-time time-dependent density functional theory to study these interactions in silver nanocrystals.

Unbiased charge oscillations in B-DNA: monomer polymers and dimer polymers

We call monomer a B-DNA base pair and examine, analytically and numerically, electron or hole oscillations in monomer and dimer polymers, i.e., periodic sequences with repetition unit made of one or two monomers. We employ a tight-binding (TB) approach at the base-pair level to readily determine the spatiotemporal evolution of a single extra carrier along a N base-pair B-DNA segment. We study highest occupied molecular orbital and lowest unoccupied molecular orbital eigenspectra as well as the mean over time probabilities to find the carrier at a particular monomer. We use the pure mean transfer rate k to evaluate the easiness of charge transfer. The inverse decay length β for exponential fits k(d), where d is the charge transfer distance, and the exponent η for power-law fits k(N) are computed; generally power-law fits are better. We illustrate that increasing the number of different parameters involved in the TB description, the fall of k(d) or k(N) becomes steeper and show the range covered by β and η. Finally, for both the time-independent and the time-dependent problems, we analyze the palindromicity and the degree of eigenspectrum dependence of the probabilities to find the carrier at a particular monomer.

Development of efficient time-evolution method based on three-term recurrence relation

The advantage of the real-time (RT) propagation method is a direct solution of the time-dependent Schrödinger equation which describes frequency properties as well as all dynamics of a molecular system composed of electrons and nuclei in quantum physics and chemistry. Its applications have been limited by computational feasibility, as the evaluation of the time-evolution operator is computationally demanding. In this article, a new efficient time-evolution method based on the three-term recurrence relation (3TRR) was proposed to reduce the time-consuming numerical procedure. The basic formula of this approach was derived by introducing a transformation of the operator using the arcsine function. Since this operator transformation causes transformation of time, we derived the relation between original and transformed time. The formula was adapted to assess the performance of the RT time-dependent Hartree-Fock (RT-TDHF) method and the time-dependent density functional theory. Compared to the commonly used fourth-order Runge-Kutta method, our new approach decreased computational time of the RT-TDHF calculation by about factor of four, showing the 3TRR formula to be an efficient time-evolution method for reducing computational cost.