First-order nonadiabatic couplings from time-dependent hybrid density functional response theory: Consistent formalism, implementation, and performance

Time-Resolved Circular Dichroism in Molecules: Experimental and Theoretical Advances

Following changes in chirality can give access to relevant information on the function or reactivity of molecular systems. Time-resolved circular dichroism (TRCD) spectroscopy proves to be a valid tool to achieve this goal. Depending on the class of molecules, different temporal ranges, spanning from seconds to femtoseconds, need to be investigated to observe such chiroptical changes. Therefore, over the years, several approaches have been adopted to cover the timescale of interest, especially based on pump-probe schemes. Moreover, various theoretical approaches have been proposed to simulate and explain TRCD spectra, including linear and non-linear response methods as well as non-adiabatic molecular dynamics. In this review, an overview on both experimental and theoretical advances in the TRCD field is provided, together with selected applications. A discussion on future theoretical developments for TRCD is also given.

Ultrafast photochemistry and electron diffraction for cyclobutanone in the S2 state: Surface hopping with time-dependent density functional theory

We simulate the photodynamics of gas-phase cyclobutanone excited to the S2 state using fewest switches surface hopping (FSSH) dynamics powered by time-dependent density functional theory (TDDFT). We predict a total photoproduct yield of 8%, with a C3:C2 product ratio of 0 trajectories to 8 trajectories. One primary S2 → S1 conical intersection is identified involving the compression of an α-carbon-carbon-hydrogen bond angle. Excited state lifetimes computed with respect to electronic state populations were found to be 3.96 ps (S2 → S1) and 498 fs (S1 → S0). We also generate time-resolved difference pair distribution functions (ΔPDFs) from our TDDFT-FSSH dynamics results in order to generate direct comparisons with ultrafast electron diffraction experiment observables. Global and target analysis of time-resolved ΔPDFs produced a distinct set of lifetimes: (i) a 0.548 ps decay and (ii) a 1.69 ps decay, both resembling the S2 minimum, as well as (iii) a long decay that resembles the S1 minimum geometry and the fully separated C2 products. Finally, we contextualize our results by considering the impact of the most likely sources of significant errors.

On the Topological Phase around Conical Intersections with Tamm-Dancoff Linear-Response Time-Dependent Density Functional Theory

Regions of nuclear-configuration space away from the Franck-Condon geometry can prove problematic for some electronic structure methods, given the propensity of such regions to possess conical intersections, i.e., (highly connected) points of degeneracy between potential energy surfaces. With the likelihood (perhaps even inevitability) for nonadiabatic dynamics simulations to explore molecular geometries in close proximity to conical intersections, it is vital that the performance of electronic structure methods is routinely examined in this context. In a recent paper [Taylor, J. T. J. Chem. Phys. 2023, 159, 214115.], the ability of linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) to provide a proper description of conical intersections, in terms of their topology and topography, was investigated, with particular attention paid to conical intersections between two excited electronic states. For the same prototypical molecules, protonated formaldimine and pyrazine, we herein consider whether AA LR-TDDFT can correctly reproduce the topological phase accumulated by the adiabatic electronic wave function upon traversing a closed path around an excited-to-excited state conical intersection despite not using the appropriate quadratic-response nonadiabatic coupling vectors. Equally, we probe the ability of the ground-to-excited state intersection ring exhibited by AA LR-TDDFT in protonated formaldimine to give rise to a similar topological phase in spite of its incorrect dimensionality.

The Excited-State Lifetime of Poly(NDI2OD-T2) Is Intrinsically Short

Conjugated polymers composed of alternating electron donor and acceptor segments have come to dominate the materials being considered for organic photoelectrodes and solar cells, in large part because of their favorable near-infrared absorption. The prototypical electron-transporting push-pull polymer poly(NDI2OD-T2) (N2200) is one such material. While reasonably efficient organic solar cells can be fabricated with N2200 as the acceptor, it generally fails to contribute as much photocurrent from its absorption bands as the donor with which it is paired. Moreover, transient absorption studies have shown N2200 to have a consistently short excited-state lifetime (∼100 ps) that is dominated by a ground-state recovery. In this paper, we investigate whether these characteristics are intrinsic to the backbone structure of this polymer or if these are extrinsic effects from ubiquitous solution-phase and thin-film aggregates. We compare the solution-phase photophysics of N2200 with those of a pair of model compounds composed of alternating bithiophene (T2) donor and naphthalene diimide (NDI) acceptor units, NDI-T2-NDI and T2-NDI-T2, in a dilute solution. We find that the model compounds have even faster ground-state recovery dynamics (τ = 45, 27 ps) than the polymer (τ = 133 ps), despite remaining molecularly isolated in solution. In these molecules, as in the case of the N2200 polymer, the lowest excited state has a T2 to NDI charge-transfer (CT) character. Electronic-structure calculations indicate that the short lifetime of this state is due to fast nonradiative decay to the ground state (GS) promoted by strong CT-GS electronic coupling and strong electron-vibrational coupling with high-frequency (quantum) normal modes.

Practical phase-space electronic Hamiltonians for ab initio dynamics

Modern electronic structure theory is built around the Born-Oppenheimer approximation and the construction of an electronic Hamiltonian Ĥel(X) that depends on the nuclear position X (and not the nuclear momentum P). In this article, using the well-known theory of electron translation (Γ') and rotational (Γ″) factors to couple electronic transitions to nuclear motion, we construct a practical phase-space electronic Hamiltonian that depends on both nuclear position and momentum, ĤPS(X,P). While classical Born-Oppenheimer dynamics that run along the eigensurfaces of the operator Ĥel(X) can recover many nuclear properties correctly, we present some evidence that motion along the eigensurfaces of ĤPS(X,P) can better capture both nuclear and electronic properties (including the elusive electronic momentum studied by Nafie). Moreover, only the latter (as opposed to the former) conserves the total linear and angular momentum in general.

Highly efficient implementation of analytic nonadiabatic derivative couplings within the pseudospectral method

A pseudospectral implementation of nonadiabatic derivative couplings in the Tamm-Dancoff approximation is reported, and the accuracy and efficiency of the pseudospectral nonadiabatic derivative couplings are studied. Our results demonstrate that the pseudospectral method provides mean absolute errors of 0.2%-1.9%, while providing a significant speedup. Benchmark calculations on fullerenes (Cn, n up to 100) using B3LYP achieved 10- to 15-fold, 8- to 17-fold, and 43- to 75-fold speedups for 6-31G**, 6-31++G**, and cc-pVTZ basis sets, respectively, when compared to the conventional spectral method.

Estimating Nonradiative Excited-State Lifetimes in Photoactive Semiconducting Nanostructures

The time evolution of the exciton generated by light adsorption in a photocatalyst is an important feature that can be approached from full nonadiabatic molecular dynamics simulations. Here, a crucial parameter is the nonradiative recombination rate between the hole and the electron that form the exciton. In the present work, we explore the performance of a Fermi's golden rule-based approach on predicting the recombination rate in a set of photoactive titania nanostructures, relying solely on the coupling of the ground and first excited state. In this scheme the analysis of the first excited state is carried out by invoking Kasha's rule thus avoiding computationally expensive nonadiabatic molecular dynamics simulations and resulting in an affordable estimate of the recombination rate. Our results show that, compared to previous ones from nonadiabatic molecular dynamics simulations, semiquantitative recombination rates can be predicted for the smaller titania nanostructures, and qualitative values are obtained from the larger ones. The present scheme is expected to be useful in the field of computational heterogeneous photocatalysis whenever a complex and computationally expensive full nonadiabatic molecular dynamics cannot be carried out.

Which Electronic Structure Method to Choose in Trajectory Surface Hopping Dynamics Simulations? Azomethane as a Case Study

Nonadiabatic dynamics simulations have become a standard approach to explore photochemical reactions. Such simulations require underlying potential energy surfaces and couplings between them, calculated at a chosen level of theory, yet this aspect is rarely assessed. Here, in combination with the popular trajectory surface hopping dynamics method, we use a high-accuracy XMS-CASPT2 electronic structure level as a benchmark for assessing the performances of various post-Hartree-Fock methods (namely, CIS, ADC(2), CC2, and CASSCF) and exchange-correlation functionals (PBE, PBE0, and CAM-B3LYP) in a TD-DFT/TDA context, using the isomerization around a double bond as test case. Different relaxation pathways are identified, and the ability of the different methods to reproduce their relative importance and time scale is discussed. The results show that multireference electronic structure methods should be preferred, when studying nonadiabatic decay between excited and ground states. If not affordable, TD-DFT with TDA and hybrid functionals and ADC(2) are efficient alternatives but overestimate the nonradiative decay yield and thus may miss deexcitation pathways.

Modeling the Fluorescence Quantum Yields of Aromatic Compounds: Benchmarking the Machinery to Compute Intersystem Crossing Rates

The from-first-principles calculation of fluorescence quantum yields (FQYs) and lifetimes of organic dyes remains very challenging. In this article, we extensively test the machinery to calculate FQYs. Specifically, we perform an extensive analysis on the parameters influencing the intersystem crossing (ISC), internal conversion (IC), and fluorescence rate constants calculations. The impact of (i) the electronic structure (chosen exchange-correlation functional and spin-orbit Hamiltonian), (ii) the vibronic parameters (coordinate system, broadening function, and dipole expansion), and (iii) the excited-state kinetic models are systematically assessed for a series of seven rigid aromatic molecules. Our studies provide more insights into the choice of parameters and the expected accuracy for the computational protocols aiming to deliver FQY values. Some challenges are highlighted, such as, on the one hand, the difficulty to benchmark against the experimental nonradiative rate constants, for which the separation between the IC and ISC contributions is often not provided in the literature and, on the other hand, the need to go beyond the harmonic approximation for the calculation of the IC rates.

On the description of conical intersections between excited electronic states with LR-TDDFT and ADC(2)

Conical intersections constitute the conceptual bedrock of our working understanding of ultrafast, nonadiabatic processes within photochemistry (and photophysics). Accurate calculation of potential energy surfaces within the vicinity of conical intersections, however, still poses a serious challenge to many popular electronic structure methods. Multiple works have reported on the deficiency of methods like linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) or algebraic diagrammatic construction to second-order [ADC(2)]-approaches often used in excited-state molecular dynamics simulations-to describe conical intersections between the ground and excited electronic states. In the present study, we focus our attention on conical intersections between excited electronic states and probe the ability of AA LR-TDDFT and ADC(2) to describe their topology and topography, using protonated formaldimine and pyrazine as two exemplar molecules. We also take the opportunity to revisit the performance of these methods in describing conical intersections involving the ground electronic state in protonated formaldimine-highlighting in particular how the intersection ring exhibited by AA LR-TDDFT can be perceived either as a (near-to-linear) seam of intersection or two interpenetrating cones, depending on the magnitude of molecular distortions within the branching space.

Nonadiabatic Dynamics Simulations for Photoinduced Processes in Molecules and Semiconductors: Methodologies and Applications

Nonadiabatic dynamics (NAMD) simulations have become powerful tools for elucidating complicated photoinduced processes in various systems from molecules to semiconductor materials. In this review, we present an overview of our recent research on photophysics of molecular systems and periodic semiconductor materials with the aid of ab initio NAMD simulation methods implemented in the generalized trajectory surface-hopping (GTSH) package. Both theoretical backgrounds and applications of the developed NAMD methods are presented in detail. For molecular systems, the linear-response time-dependent density functional theory (LR-TDDFT) method is primarily used to model electronic structures in NAMD simulations owing to its balanced efficiency and accuracy. Moreover, the efficient algorithms for calculating nonadiabatic coupling terms (NACTs) and spin-orbit couplings (SOCs) have been coded into the package to increase the simulation efficiency. In combination with various analysis techniques, we can explore the mechanistic details of the photoinduced dynamics of a range of molecular systems, including charge separation and energy transfer processes in organic donor-acceptor structures, ultrafast intersystem crossing (ISC) processes in transition metal complexes (TMCs), and exciton dynamics in molecular aggregates. For semiconductor materials, we developed the NAMD methods for simulating the photoinduced carrier dynamics within the framework of the Kohn-Sham density functional theory (KS-DFT), in which SOC effects are explicitly accounted for using the two-component, noncollinear DFT method. Using this method, we have investigated the photoinduced carrier dynamics at the interface of a variety of van der Waals (vdW) heterojunctions, such as two-dimensional transition metal dichalcogenides (TMDs), carbon nanotubes (CNTs), and perovskites-related systems. Recently, we extended the LR-TDDFT-based NAMD method for semiconductor materials, allowing us to study the excitonic effects in the photoinduced energy transfer process. These results demonstrate that the NAMD simulations are powerful tools for exploring the photodynamics of molecular systems and semiconductor materials. In future studies, the NAMD simulation methods can be employed to elucidate experimental phenomena and reveal microscopic details as well as rationally design novel photofunctional materials with desired properties.

The carbonyl-lock mechanism underlying non-aromatic fluorescence in biological matter

Challenging the basis of our chemical intuition, recent experimental evidence reveals the presence of a new type of intrinsic fluorescence in biomolecules that exists even in the absence of aromatic or electronically conjugated chemical compounds. The origin of this phenomenon has remained elusive so far. In the present study, we identify a mechanism underlying this new type of fluorescence in different biological aggregates. By employing non-adiabatic ab initio molecular dynamics simulations combined with a data-driven approach, we characterize the typical ultrafast non-radiative relaxation pathways active in non-fluorescent peptides. We show that the key vibrational mode for the non-radiative decay towards the ground state is the carbonyl elongation. Non-aromatic fluorescence appears to emerge from blocking this mode with strong local interactions such as hydrogen bonds. While we cannot rule out the existence of alternative non-aromatic fluorescence mechanisms in other systems, we demonstrate that this carbonyl-lock mechanism for trapping the excited state leads to the fluorescence yield increase observed experimentally, and set the stage for design principles to realize novel non-invasive biocompatible probes with applications in bioimaging, sensing, and biophotonics.

Exploring Thioxanthone Derivatives as Singlet Oxygen Photosensitizers for Photodynamic Therapy at the Near-IR Region

In the lowest excited triplet state, the excited photosensitizer reacts with tissue oxygen and forms reactive oxygen species (ROS), which kills tissue cells in photodynamic therapy (PDT). Metal-free thio-based pure organic molecules and analogous nucleobases can be used as photosensitizers for PDT applications. Using quantum chemical methods, we studied one- and two-photon optical absorptions, fluorescence, and other excited-state properties of substituted thioxanthone derivatives for their potential as photosensitizers for PDT. Our calculated values were compared with the available experimental data. The calculation of the intersystem crossing rate constant for these photosensitizers explains the high quantum yield of the formation of ROS, as reported experimentally. The excited triplet-state population of the photosensitizer occurs through the 1π-π* → 3n-π* channel of intersystem crossing and increases in the presence of halogen substitution.

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Surface hopping, electron translation factors, electron rotation factors, momentum conservation, and size consistency

For a system without spin-orbit coupling, the (i) nuclear plus electronic linear momentum and (ii) nuclear plus orbital electronic angular momentum are good quantum numbers. Thus, when a molecular system undergoes a nonadiabatic transition, there should be no change in the total linear or angular momentum. Now, the standard surface hopping algorithm ignores the electronic momentum and indirectly equates the momentum of the nuclear degrees of freedom to the total momentum. However, even with this simplification, the algorithm still does not conserve either the nuclear linear or the nuclear angular momenta. Here, we show that one way to address these failures is to dress the derivative couplings (i.e., the hopping directions) in two ways: (i) we disallow changes in the nuclear linear momentum by working in a translating basis (which is well known and leads to electron translation factors) and (ii) we disallow changes in the nuclear angular momentum by working in a basis that rotates around the center of mass [which is not well-known and leads to a novel, rotationally removable component of the derivative coupling that we will call electron rotation factors below, cf. Eq. (96)]. The present findings should be helpful in the short term as far as interpreting surface hopping calculations for singlet systems (without spin) and then developing the new surface hopping algorithm in the long term for systems where one cannot ignore the electronic orbital and/or spin angular momentum.

Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics

When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.

Ab Initio Simulation of the Ultrafast Circular Dichroism Spectrum of Provitamin D Ring-Opening

We present a method to simulate ultrafast pump-probe time-resolved circular dichroism (TRCD) spectra based on time-dependent density functional theory trajectory surface hopping. The method is applied to simulate the TRCD spectrum along the photoinduced ring-opening of provitamin D. Simulations reveal that the initial decay of the signal is due to excited state relaxation, forming the rotationally flexible previtamin D. We further show that oscillations in the experimental TRCD spectrum arise from isomerizations between previtamin D rotamers with different chirality, which are associated with the helical conformation of the triene unit. We give a detailed description of the formation dynamics of different rotamers, playing a key role in the natural regulation of vitamin D photosynthesis. Going beyond the sole extraction of decay rates, simulations greatly increase the amount of information that can be retrieved from ultrafast TRCD, making it a sensitive tool to unravel details in the subpicosecond dynamics of photoinduced chirality changes.

Calculation of Excited State Internal Conversion Rate Constant Using the One-Effective Mode Marcus-Jortner-Levich Theory

In this article, the one-effective mode Marcus-Jortner-Levich (MJL) theory and the classical Marcus theory for electron transfer were applied to estimate the internal conversion rate constant, k_IC, of organic molecules and a Ru-based complex, all belonging to the Marcus inverted region. For this, the reorganization energy was calculated using the minimum energy conical intersection point to account for more vibrational levels, correcting the density of states. The results showed good agreement with experimental and theoretically determined k_IC, with a small overestimation by the Marcus theory. Also, molecules less dependent on the solvent effects, like benzophenone, presented better results than molecules with an expressive dependence, like 1-aminonaphthalene. Moreover, the results suggest that each molecule possesses unique normal modes leading to the excited state deactivation that does not necessarily match the X-H bond stretching, as previously suggested.

Internal conversion rate constant calculations considering Duschinsky, anharmonic and Herzberg-Teller effects

A novel method for calculating rate constants for internal conversion (k_IC) that simultaneously accounts for Duschinsky, anharmonic and Herzberg-Teller effects has been developed and implemented. This method has been applied to robust planar molecules like tetraoxa[8]circulene (4B), free-base porphyrin (H2P) and pyrometene (PM567) with small Duschinsky rotation (i.e. with almost identical normal coordinates in the ground and excited states) and to poly[n]fluorenes (P[n]F) (n = 2-14) with a substantial Duschinsky rotation. The obtained results show that the Duschinsky effect is large in the harmonic approximation, whereas it is in general much smaller in the anharmonic approximation. The Duschinsky effect is found to be large for high frequency vibrational modes with energies of ∼3300 cm^-1 such as the X-H (X = C, N and O) stretching modes that mix in the S₁ → S₀ electronic transition. However, even in this case, the increase in k_IC due to the Duschinsky effect does not exceed one order of magnitude. The calculations show that anharmonic contributions to k_IC are larger than Herzberg-Teller contributions which in turn are larger than contributions from the Duschinsky effect ANH > HT > Du. We also show that an approximation, where only X-H bonds are considered in the k_IC calculation, is accurate even for P[n]F (n = 2-14).

Exact non-adiabatic coupling vectors for the time-dependent density functional based tight-binding method

We report on non-adiabatic coupling vectors between electronic excited states for the time-dependent-density functional theory based tight-binding (TD-DFTB) method. The implementation includes orbital relaxation effects that have been previously neglected and covers also the case of range-separated exchange-correlation functionals. Benchmark calculations with respect to first principles TD-DFT highlight the large dependence of non-adiabatic couplings on the functional. Closer investigations of the topology around a conical intersection between excited states show that TD-DFTB delivers near-exact values of the Berry phase, which paves the way for consistent non-adiabatic molecular dynamics simulations for large systems.

Analytical derivative couplings within the framework of time-dependent density functional theory coupled with conductor-like polarizable continuum model: Formalism, implementation, and applications

The nonadiabatic phenomena, which are characterized by a strong coupling between electronic and nuclear motions, are ubiquitous. The nonadiabatic effect of the studied system can be significantly affected by the surrounding environment, such as solvents, in which such nonadiabatic process takes place. It is essential to develop the theoretical models to simulate these processes while accurately modeling the solvent environment. The time-dependent density functional theory (TDDFT) is currently the most efficient approach to describe the electronic structures and dynamics of complex systems, while the polarizable continuum model (PCM) represents one of the most successful examples among continuum solvation models. Here, we formulate the first-order derivative couplings (DCs) between the ground and excited states as well as between two excited states by utilizing time-independent equation of motion formalism within the framework of both linear response and spin flip formulations of TDDFT/CPCM (the conductor-like PCM), and implement the analytical DCs into the Q-CHEM electronic structure software package. The analytic implementation is validated by the comparison of the analytical and finite-difference results, and reproducing geometric phase effect in the protonated formaldimine test case. Taking 4-(N,N-dimethylamino)benzonitrile and uracil in the gas phase and solution as an example, we demonstrate that the solvent effect is essential not only for the excitation energies of the low-lying excited-states but also for the DCs between these states. Finally, we calculate the internal conversion rate of benzophenone in a solvent with DC being used. The current implementation of analytical DCs together with the existing analytical gradient and Hessian of TDDFT/PCM excited states allows one to study the nonadiabatic effects of relatively large systems in solutions with low computational cost.

Exciton dispersion and exciton-phonon interaction in solids by time-dependent density functional theory

Understanding, predicting, and ultimately controlling exciton band structure and exciton dynamics are central to diverse chemical and materials problems. Here, we have developed a first-principles method to determine exciton dispersion and exciton-phonon interaction in semiconducting and insulating solids based on time-dependent density functional theory. The first-principles method is formulated in planewave bases and pseudopotentials and can be used to compute exciton band structures, exciton charge density, ionic forces, the non-adiabatic coupling matrix between excitonic states, and the exciton-phonon coupling matrix. Based on the spinor formulation, the method enables self-consistent noncollinear calculations to capture spin-orbital coupling. Hybrid exchange-correlation functionals are incorporated to deal with long-range electron-hole interactions in solids. A sub-Hilbert space approximation is introduced to reduce the computational cost without loss of accuracy. For validations, we have applied the method to compute the exciton band structure and exciton-phonon coupling strength in transition metal dichalcogenide monolayers; both agree very well with the previous GW-Bethe-Salpeter equation and experimental results. This development paves the way for accurate determinations of exciton dynamics in a wide range of solid-state materials.

Interexcited State Photophysics I: Benchmarking Density Functionals for Computing Nonadiabatic Couplings and Internal Conversion Rate Constants

We present the first benchmarking study of nonadiabatic matrix coupling elements (NACMEs) calculated using different density functionals. Using the S₁ → S₀ transition in perylene solvated in toluene as a case study, we calculate the photophysical properties and corresponding rate constants for a variety of density functionals from each rung of Jacob's ladder. The singlet photoluminescence quantum yield (sPLQY) is taken as a measure of accuracy, measured experimentally here as 0.955. Important quantum chemical parameters such as geometries, absorption, emission, and adiabatic energies, NACMEs, Hessians, and transition dipole moments were calculated for each density functional basis set combination (data set) using density functional theory based multireference configuration interaction (DFT/MRCI) and compared to experiment where possible. We were able to derive simple relations between the TDDFT and DFT/MRCI photophysical properties; with semiempirical damping factors of ∼0.843 ± 0.017 and ∼0.954 ± 0.064 for TDDFT transition dipole moments and energies to DFT/MRCI level approximations, respectively. NACMEs were dominated by out-of-plane derivative components belonging to the center-most ring atoms with weaker contributions from perturbations along the transverse and longitudinal axes. Calculated theoretical spectra compared well to both experiment and literature, with fluorescence lifetimes between 7.1 and 12.5 ns, agreeing within a factor of 2 with experiment. Internal conversion (IC) rates were then calculated and were found to vary wildly between 10⁶-10¹⁶ s^-1 compared with an experimental rate of the order 10⁷ s^-1. Following further testing by mixing data sets, we found a strong dependence on the method used to obtain the Hessian. The 5 characterized data sets ranked in order of most promising are PBE0/def2-TZVP, ωB97XD/def2-TZVP, HCTH407/TZVP, PBE/TZVP, and PBE/def2-TZVP.

Investigation of plasmon relaxation mechanisms using nonadiabatic molecular dynamics

Hot carriers generated from the decay of plasmon excitation can be harvested to drive a wide range of physical or chemical processes. However, their generation efficiency is limited by the concomitant phonon-induced relaxation processes by which the energy in excited carriers is transformed into heat. However, simulations of dynamics of nanoscale clusters are challenging due to the computational complexity involved. Here, we adopt our newly developed Trajectory Surface Hopping (TSH) nonadiabatic molecular dynamics algorithm to simulate plasmon relaxation in Au20 clusters, taking the atomistic details into account. The electronic properties are treated within the Linear Response Time-Dependent Tight-binding Density Functional Theory (LR-TDDFTB) framework. The relaxation of plasmon due to coupling to phonon modes in Au20 beyond the Born–Oppenheimer approximation is described by the TSH algorithm. The numerically efficient LR-TDDFTB method allows us to address a dense manifold of excited states to ensure the inclusion of plasmon excitation. Starting from the photoexcited plasmon states in Au20 cluster, we find that the time constant for relaxation from plasmon excited states to the lowest excited states is about 2.7 ps, mainly resulting from a stepwise decay process caused by low-frequency phonons of the Au20 cluster. Furthermore, our simulations show that the lifetime of the phonon-induced plasmon dephasing process is ∼10.4 fs and that such a swift process can be attributed to the strong nonadiabatic effect in small clusters. Our simulations demonstrate a detailed description of the dynamic processes in nanoclusters, including plasmon excitation, hot carrier generation from plasmon excitation dephasing, and the subsequent phonon-induced relaxation process.

Modeling Photolytic Decomposition of Energetically Functionalized Dodecanes

The photolytic stability of explosives and energetic functional groups is of importance for those who regularly handle or are exposed to explosives in typical environmental conditions. This study models the photolytic degradation of dodecane substituted with various energetic functional groups: azide, nitro, nitrate ester, and nitramine. For the studied molecules, it was found that excitons localize on the energetic functional group, no matter where they were initially formed, and thus, the predominant degradation pathway involves the degradation of the energetic functional group. The relative trends for both 4 and 8 eV excitation energies followed with what is expected from the relative stability of the energetic functional groups to thermal and sub-shock degradation. The one notable exception was the azide functional group; more work should be done to further understand the photolytic effects on the azide functional group.

Accurate Non-Adiabatic Couplings from Optimally-Tuned Range-Separated Hybrid Functionals

Precise theoretical calculations of non-adiabatic couplings, which describe the interaction between two Born-Oppenheimer surfaces, are important for the modeling of radiationless decay mechanisms in photochemical processes. Here, we demonstrate that accurate non-adiabatic couplings can be calculated in the framework of linear-response time-dependent density functional theory by using non-empirical, optimally-tuned range-separated hybrid (OT-RSH) functionals. We focus on molecular radicals, in which ultrafast non-radiative decay plays a crucial role, to find that the OT-RSH functional compares well to wave-function based reference data and competes with the accuracy of semi-empirical CAM-B3LYP calculations. Our findings show that the OT-RSH approach provides very accurate non-adiabatic couplings and, therefore, provides a computationally efficient alternative to wave-function based techniques.

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. The nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

Evaluation of molecular photophysical and photochemical properties using linear response time-dependent density functional theory with classical embedding: Successes and challenges

Time-dependent density functional theory (TDDFT) based approaches have been developed in recent years to model the excited-state properties and transition processes of the molecules in the gas-phase and in a condensed medium, such as in a solution and protein microenvironment or near semiconductor and metal surfaces. In the latter case, usually, classical embedding models have been adopted to account for the molecular environmental effects, leading to the multi-scale approaches of TDDFT/polarizable continuum model (PCM) and TDDFT/molecular mechanics (MM), where a molecular system of interest is designated as the quantum mechanical region and treated with TDDFT, while the environment is usually described using either a PCM or (non-polarizable or polarizable) MM force fields. In this Perspective, we briefly review these TDDFT-related multi-scale models with a specific emphasis on the implementation of analytical energy derivatives, such as the energy gradient and Hessian, the nonadiabatic coupling, the spin-orbit coupling, and the transition dipole moment as well as their nuclear derivatives for various radiative and radiativeless transition processes among electronic states. Three variations of the TDDFT method, the Tamm-Dancoff approximation to TDDFT, spin-flip DFT, and spin-adiabatic TDDFT, are discussed. Moreover, using a model system (pyridine-Ag20 complex), we emphasize that caution is needed to properly account for system-environment interactions within the TDDFT/MM models. Specifically, one should appropriately damp the electrostatic embedding potential from MM atoms and carefully tune the van der Waals interaction potential between the system and the environment. We also highlight the lack of proper treatment of charge transfer between the quantum mechanics and MM regions as well as the need for accelerated TDDFT modelings and interpretability, which calls for new method developments.

Ultrafast coherent photoexcited dynamics in a trimeric dendrimer probed by X-ray stimulated-Raman signals

The photoinduced ultrafast coherent inter-chromophore energy redistribution in a triarylamine trimer is explored using nonadiabatic excited state molecular dynamics followed by simulations of X-ray Raman signals. The nitrogencentered system ensures strong interchromophore interactions and, thus, the presence of coherences. Nevertheless, the multitude of non-deterministic photoinduced pathways during the ultrafast inter-branch migration of the excitation results in random confinement on some branches and, therefore, spatial exciton scrambling and loss of phase information at long times. We show that the vibronic coherence dynamics evolving into the incoherent scrambling mechanism on ultrafast 50 fs timescale, is accurately probed by the TRUECARS X-ray stimulated Raman signal. In combination with previous results, where the technique has revealed long-lived coherences in a rigid heterodimer, the signal is most valuable for detecting ultrafast molecular coherences or their absence. We demonstrate that X-ray Raman spectroscopy is a useful tool in the chemical design of functional molecular building blocks.

The photoinduced ultrafast coherent inter-chromophore energy redistribution in a triarylamine trimer is explored using nonadiabatic excited state molecular dynamics followed by simulations of X-ray Raman signals.

Fewest switches surface hopping with Baeck-An couplings

In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.

Fewest switches surface hopping with Baeck-An couplings

In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.

Exciton Dynamics of a Diketo-Pyrrolopyrrole Core for All Low-Lying Electronic Excited States Using Density Functional Theory-Based Methods

Ab initio treatments of interexcited state internal conversion (IC) are more often than not missing from exciton dynamic descriptions, because of their inherent complexity. Here, we define "interexcited state IC" as a same-spin nonradiative transition between states i and j, where i ≠ j ≠ 0. Competing directly with multiexciton processes such as singlet fission or triplet photoupconversion, inclusion of this mechanism in the narrative of molecular photophysics would allow for strategic synthesis of chromophores for more efficient photon-harvesting applications. Herein, we present a robust formalism which can model these rates using density functional theory (DFT)-based methods within the Franck-Condon and Herzberg-Teller regime. Using an unsubstituted diketo-pyrrolopyrrole (DPP) core as a case study, we illustrate the exciton dynamics along the first four excited states for both singlet and triplet manifolds, showing ultrafast same-spin transfer mechanisms due to all excited states, excluding the first triplet level, being in close energetic proximity (within 0.8 eV of each other). The resulting electron same-spin rates outcompete the electron spin-flipping intersystem crossing (ISC) rates, with excitons firmly obeying Kasha's rule as they cascade down from the high-lying excited states toward the lower states. Furthermore, we calculated that only the first singlet excited state displayed a reasonable probability of triplet exciton generation, of ∼40%, with a near-zero chance of the exciton reverting to the singlet manifold once the electron-hole pair are of parallel spin.

A theoretical study on the donor ability adjustment of tris(2,4,6-trichlorophenyl)methyl-triarylamine (TTM-TPA) radicals aiming to develop better organic luminescent materials

Spin-unrestricted DFT and spin-unrestricted TDDFT calculations were performed to systematically investigate the correlation between the electron donating ability of donors and photophysical properties in D–A luminescent radicals.

The combination of skeleton-engineering and periphery-engineering: a design strategy for organic doublet emitters

A series of radicals based on tris(2,4,6-trichlorophenyl)methyl (TTM) were theoretically designed and evaluated by combining skeleton-engineering and periphery-engineering strategies.

Theoretical Investigation of the Effect of Alkylation and Bromination on Intersystem Crossing in BODIPY-Based Photosensitizers

Halogenated and alkylated BODIPY derivatives are reported as suitable candidates for their use as photosensitizers in photodynamic therapy due to their efficient intersystem crossing (ISC) between states of different spin multiplicities. Spin-orbit couplings (SOCs) are evaluated using an effective one-electron spin-orbit Hamiltonian for brominated and alkylated BODIPY derivatives to investigate the quantitative effect of alkyl and bromine substituents on ISC. BODIPY derivatives containing bromine atoms have been found to have significantly stronger SOCs than alkylated BODIPY derivatives outside the Frank-Condon region while they are nearly the same at local minima. Based on calculated time-dependent density functional theory (TD-DFT) vertical excitation energies and SOCs, excited-state dynamics of three BODIPY derivatives were further explored with TD-DFT surface hopping molecular dynamics employing a simple accelerated approach. Derivatives containing bromine atoms have been found to have very similar lifetimes, which are much shorter than those of the derivatives possessing just the alkyl moieties. However, both bromine atoms and alkyl moieties reduce the HOMO/LUMO gap, thus assisting the derivatives to behave as efficient photosensitizers.

On the inclusion of one double within CIS and TDDFT

We present an improved approach for generating a set of optimized frontier orbitals (HOMO and LUMO) that minimizes the energy of one double configuration. We further benchmark the effect of including such a double within a rigorous configuration interaction singles or a parameterized semi-empirical time-dependent density functional theory Hamiltonian for a set of test cases. Although we cannot quite achieve quantitative accuracy, the algorithm is quite robust and routinely delivers an enormous qualitative improvement to standard single-reference electronic structure calculations.

Nonadiabatic dynamics with spin-flip vs linear-response time-dependent density functional theory: A case study for the protonated Schiff base C5H6NH2. J Chem Phys 2021;155:124111. [PMID: 34598550 DOI: 10.1063/5.0062757]

Nonadiabatic trajectory surface hopping simulations are reported for trans-C₅H₆NH₂ ⁺, a model of the rhodopsin chromophore, using the augmented fewest-switches algorithm. Electronic structure calculations were performed using time-dependent density functional theory (TDDFT) in both its conventional linear-response (LR) and its spin-flip (SF) formulations. In the SF-TDDFT case, spin contamination in the low-lying singlet states is removed by projecting out the lowest triplet component during iterative solution of the TDDFT eigenvalue problem. The results show that SF-TDDFT qualitatively describes the photoisomerization of trans-C₅H₆NH₂ ⁺, with favorable comparison to previous studies using multireference electronic structure methods. In contrast, conventional LR-TDDFT affords qualitatively different photodynamics due to an incorrect excited-state potential surface near the Franck-Condon region. In addition, the photochemistry (involving pre-twisting of the central double bond) appears to be different for SF- and LR-TDDFT, which may be a consequence of different conical intersection topographies afforded by these two methods. The present results contrast with previous surface-hopping studies suggesting that the LR-TDDFT method's incorrect topology around S₁/S₀ conical intersections is immaterial to the photodynamics.

Reduced-dimensional surface hopping with offline-online computations

Molecular dynamics simulations often classically evolve the nuclear geometry on adiabatic potential energy surfaces (PESs), punctuated by random hops between energy levels in regions of strong coupling, in an algorithm known as surface hopping. However, the computational expense of integrating the geometry on a full-dimensional PES and computing the required couplings can quickly become prohibitive as the number of atoms increases. In this work, we describe a method for surface hopping that uses only important reaction coordinates, performs all expensive evaluations of the true PESs and couplings only once before simulating dynamics (offline), and then queries the stored values during the surface hopping simulation (online). Our Python codes are freely available on GitHub. Using photodissociation of azomethane as a test case, this method is able to reproduce experimental results that have thus far eluded ab initio surface hopping studies.

NAC-TDDFT: Time-Dependent Density Functional Theory for Nonadiabatic Couplings

First-order nonadiabatic coupling (NAC) matrix elements (fo-NACMEs) are the basic quantities in theoretical descriptions of electronically nonadiabatic processes that are ubiquitous in molecular physics and chemistry. Given the large size of systems of chemical interests, time-dependent density functional theory (TDDFT) is usually the first choice of methods. However, the lack of many-electron wave functions in TDDFT renders the formulation of NAC-TDDFT for fo-NACMEs conceptually difficult. Because of this, various variants of NAC-TDDFT have been proposed in the literature from different standing points, including the Hellmann-Feynman-like expression and auxiliary/pseudo-wave function (AWF)-, equation-of-motion (EOM)-, and time-dependent perturbation theory (TDPT)-based formulations. Based on critical analyses, the following conclusions are made here: (1) The Hellmann-Feynman-like expression, which is rooted in exact wave function theory, is hardly useful due to huge demand on basis sets. (2) Although most popular, the AWF variants of NAC-TDDFT are not theoretically founded and become ambiguous particularly for the fo-NACMEs between two excited states, although they do agree with the EOM and TDPT variants under the Tamm-Dancoff approximation. (3) The TDPT variant of NAC-TDDFT is theoretically most rigorous but suffers from numerical instabilities on the one hand and does not differ to a significant extent from the EOM variant on the other hand. (4) As such, the EOM variant of NAC-TDDFT for the fo-NACMEs between the ground and excited states and between two excited states is solely the right choice in practice. These formal analyses are fully supported by numerical experimentations, taking azulene as a showcase. The proper implementation of the EOM variant of NAC-TDDFT is also highlighted, showing that the fo-NACMEs between the ground and excited states and between two excited states are computationally very much the same as the analytic energy gradients of DFT and TDDFT, respectively. Possible future developments of the EOM variant of NAC-TDDFT are also highlighted. Its extensions to spin-adapted open-shell TDDFT and proper treatment of spin-orbit couplings (which are another source of force for electronically nonadiabatic processes) are particularly warranted in the near future.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections

Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.

Highly efficient implementation of the analytical gradients of pseudospectral time-dependent density functional theory

The accuracy and efficiency of time-dependent density functional theory (TDDFT) excited state gradient calculations using the pseudospectral method are presented. TDDFT excited state geometry optimizations of the G2 test set molecules, the organic fluorophores with large Stokes shifts, and the Pt-complexes show that the pseudospectral method gives average errors of 0.01-0.1 kcal/mol for the TDDFT excited state energy, 0.02-0.06 pm for the bond length, and 0.02-0.12° for the bond angle when compared to the results from conventional TDDFT. TDDFT gradient calculations of fullerenes (C_n, n up to 540) with the B3LYP functional and 6-31G** basis set show that the pseudospectral method provides 8- to 14-fold speedups in the total wall clock time over the conventional methods. The pseudospectral TDDFT gradient calculations with a diffuse basis set give higher speedups than the calculations for the same basis set without diffuse functions included.

Internal conversion of singlet and triplet states employing numerical DFT/MRCI derivative couplings: Implementation, tests, and application to xanthone

We present an efficient implementation of nonadiabatic coupling matrix elements (NACMEs) for density functional theory/multireference configuration interaction (DFT/MRCI) wave functions of singlet and triplet multiplicity and an extension of the Vibes program that allows us to determine rate constants for internal conversion (IC) in addition to intersystem crossing (ISC) nonradiative transitions. Following the suggestion of Plasser et al. [J. Chem. Theory Comput. 12, 1207 (2016)], the derivative couplings are computed as finite differences of wave function overlaps. Several measures have been taken to speed up the calculation of the NACMEs. Schur's determinant complement is employed to build up the determinant of the full matrix of spin-blocked orbital overlaps from precomputed spin factors with fixed orbital occupation. Test calculations on formaldehyde, pyrazine, and xanthone show that the mutual excitation level of the configurations at the reference and displaced geometries can be restricted to 1. In combination with a cutoff parameter of t_norm = 10^-8 for the DFT/MRCI wave function expansion, this approximation leads to substantial savings of cpu time without essential loss of precision. With regard to applications, the photoexcitation decay kinetics of xanthone in apolar media and in aqueous solution is in the focus of the present work. The results of our computational study substantiate the conjecture that S₁ T₂ reverse ISC outcompetes the T₂ ↝ T₁ IC in aqueous solution, thus explaining the occurrence of delayed fluorescence in addition to prompt fluorescence.

Interfacial photoinduced carrier dynamics tuned by polymerization of coronene molecules encapsulated in carbon nanotubes: bridging type-I and type-II heterojunctions

Carbon nanomaterials play important roles in modern scientific research. Integrating different carbon-based building blocks into nano-hybrid architectures not only takes full advantage of each component, but also brings in novel interfacial properties. Herein, we have employed density functional theory (DFT) calculations to investigate the effects of polymerization degree of coronene molecules encapsulated in single-walled carbon nanotubes (SWNTs) (19,0) on their interfacial properties. The present results reveal that the interfacial properties of the formed heterojunctions are remarkably regulated by the polymerization degree. For example, monomer- and dimer-encapsulated SWNTs are type-I heterojunctions in which interfacial excitation energy transfer is preferred, whereas interfacial charge carrier transfer is favorable in trimer- and polymer-encapsulated SWNTs because they are type-II heterojunctions. On the other hand, we have employed the time-domain nonadiabatic dynamics simulation approach to explore the interfacial carrier dynamics in type-II polymer-encapsulated SWNT heterojunctions. It is found that the electron and hole transfer processes are asymmetric and occur in opposite directions and at different rates. The former takes place from polymers to SWNTs in an ultrafast way (ca. 370 fs), whereas the latter occurs slowly from SWNTs to polymers (ca. 24 ps). A closer analysis uncovers the fact that the different carrier transfer rates mainly originate from the different densities of the acceptor states, energy differences and inter-state couplings between the donor and acceptor states. Finally, the present work demonstrates that the polymerization degree could act as a new regulating strategy to tune the interfacial properties of molecule-encapsulated SWNT heterojunctions.

Fast and Accurate Computation of Nonadiabatic Coupling Matrix Elements Using the Truncated Leibniz Formula and Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory

We present a fast and accurate numerical algorithm for computing the first-order nonadiabatic coupling matrix element (NACME). The algorithm employs the truncated Leibniz formula (TLF) approximation within the finite-difference method, which makes it easily applicable in connection with any wave function-based methodology. In this work, we used the algorithm in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT, MRSF for brevity). The accuracy is assessed for NACME between the singlet electronic states of a dissociating hydrogen molecule. It is demonstrated that an intermediate approximation, TLF(1), affords a negligible numeric error on the order of ∼10^-10 a.u. while enabling a fast computation of NACME. As the MRSF method yields the correct description of the dissociation curves of H₂ for all the electronic states involved, the numeric TLF(1)/MRSF NACME values are in excellent agreement with the reference analytical values obtained by the full configuration interaction. For polyatomic molecules, the MRSF NAC vectors agree very closely with the MRCISD NAC vectors. Hence, the proposed protocol is a promising tool for the evaluation of NACMEs.

Exploring the Franck-Condon region of a photoexcited charge transfer complex in solution to interpret femtosecond stimulated Raman spectroscopy: excited state electronic structure methods to unveil non-radiative pathways

We present electronic structure methods to unveil the non-radiative pathways of photoinduced charge transfer (CT) reactions that play a main role in photophysics and light harvesting technologies. A prototypical π-stacked molecular complex consisting of an electron donor (1-chloronaphthalene, 1ClN) and an electron acceptor (tetracyanoethylene, TCNE) was investigated in dichloromethane solution for this purpose. The characterization of TCNE:π:1ClN in both its equilibrium ground and photoinduced low-lying CT electronic states was performed by using a reliable and accurate theoretical-computational methodology exploiting ab initio molecular dynamics simulations. The structural and vibrational time evolution of key vibrational modes is found to be in excellent agreement with femtosecond stimulated Raman spectroscopy experiments [R. A. Mathies et al., J. Phys. Chem. A, 2018, 122, 14, 3594], unveiling a correlation between vibrational fingerprints and electronic properties. The evaluation of nonadiabatic coupling matrix elements along generalized normal modes has made possible the interpretation on the molecular scale of the activation of nonradiative relaxation pathways towards the ground electronic state. In particular, two low frequency vibrational modes such as the out of plane bending and dimer breathing and the TCNE central C[double bond, length as m-dash]C stretching play a prominent role in relaxation phenomena from the electronic CT state to the ground state one.

Efficient enumeration of bosonic configurations with applications to the calculation of non-radiative rates

This work presents algorithms for the efficient enumeration of configuration spaces following Boltzmann-like statistics, with example applications to the calculation of non-radiative rates, and an open-source implementation. Configuration spaces are found in several areas of physics, particularly wherever there are energy levels that possess variable occupations. In bosonic systems, where there are no upper limits on the occupation of each level, enumeration of all possible configurations is an exceptionally hard problem. We look at the case where the levels need to be filled to satisfy an energy criterion, for example, a target excitation energy, which is a type of knapsack problem as found in combinatorics. We present analyses of the density of configuration spaces in arbitrary dimensions and how particular forms of kernel can be used to envelope the important regions. In this way, we arrive at three new algorithms for enumeration of such spaces that are several orders of magnitude more efficient than the naive brute force approach. Finally, we show how these can be applied to the particular case of internal conversion rates in a selection of molecules and discuss how a stochastic approach can, in principle, reduce the computational complexity to polynomial time.