Communication: Exciton analysis in time-dependent density functional theory: How functionals shape excited-state characters

Optimal Transport Distances to Characterize Electronic Excitations

Understanding the character of electronic excitations is important in computational and reaction mechanistic studies, but their classification from simulations remains an open problem. Distances based on optimal transport have proven very useful in a plethora of classification problems and, therefore, seem a natural tool to try to tackle this challenge. We propose and investigate a new diagnostic Θ based on the Sinkhorn divergence from optimal transport. We evaluate a k-NN classification algorithm on Θ, the popular Λ diagnostic, and their combination, and assess their performance in labeling excitations, finding that (i) the combination only slightly improves the classification, (ii) Rydberg excitations are not separated well in any setting, and (iii) Θ breaks down for charge transfer in small molecules. We then define a length-scale-normalized version of Θ and show that the result correlates closely with Λ for results obtained with Gaussian basis functions. Finally, we discuss the orbital dependence of our approach and explore an orbital-independent version. Using an optimized combination of the optimal transport and overlap diagnostics together with a different metric is in our opinion the most promising for future classification studies.

PyrAtes: Modular Organic Salts with Large Stokes Shifts for Fluo-rescence Microscopy

The deployment of small-molecule fluorescent agents plays an ever-growing role in medicine and drug development. Herein, we complement the portfolio of powerful fluorophores, reporting the serendipitous discovery and development of a novel class with an imidazo[1,2-a]pyridinium triflate core, which we term PyrAtes. These fluorophores are synthesized in a single step from readily available materials (>60 examples) and display Stokes shifts as large as 240 nm, while also reaching NIR-I emissions at λmax as long as 720 nm. Computational studies allow the development of a platform for the prediction of λmax and λEm. Furthermore, we demonstrate the compatibility of these novel fluorophores with live cell imaging in HEK293 cells, suggesting PyrAtes as potent intracellular markers.

Earth Mover's Charge Transfer Distance: A General and Robust Approach for Describing Excited State Locality

A novel approach for assessing the extent of electron displacement in optical transitions is proposed by implementing the Earth Mover's Distance (EMD) method, which quantifies the spatial dissimilarity between ground and excited state electron density distributions. In contrast to previous descriptors, this index provides a representative and intuitively understandable distance under a robust and computationally efficient scheme for all possible forms of locality, even in the most difficult to dissect topological cases. The theoretical differences among the existing indices and our method are first illustrated with the help of a simplified model system, followed by a benchmarking of several partial atomic charge models using experimentally relevant push-pull compounds with diverse symmetries. These same molecules are finally employed to further demonstrate the principal advantages of the EMD index and its capabilities in rationalizing charge transfer phenomena.

Visualizing and characterizing excited states from time-dependent density functional theory

Time-dependent density functional theory (TD-DFT) is the most widely-used electronic structure method for excited states, due to a favorable combination of low cost and semi-quantitative accuracy in many contexts, even if there are well recognized limitations. This Perspective describes various ways in which excited states from TD-DFT calculations can be visualized and analyzed, both qualitatively and quantitatively. This includes not just orbitals and densities but also well-defined statistical measures of electron-hole separation and of Frenkel-type exciton delocalization. Emphasis is placed on mathematical connections between methods that have often been discussed separately. Particular attention is paid to charge-transfer diagnostics, which provide indicators of when TD-DFT may not be trustworthy due to its categorical failure to describe long-range electron transfer. Measures of exciton size and charge separation that are directly connected to the underlying transition density are recommended over more ad hoc metrics for quantifying charge-transfer character.

Tailoring D-π-A architectures with hybridized local and charge transfer fluorophores exhibiting high electroluminescence exciton utilization and low threshold amplified spontaneous emission

Novel amorphous compounds which could simultaneously use 25% singlet excitons and 75% triplet excitons as the energy source for light amplification enable the reduction of the threshold current density for electrically pumped organic semiconductor laser diodes (OSLDs); however, there is always a trade-off between the high radiative decay rate of the local excited (LE) state that is required for amplified spontaneous emission (ASE) and high exciton utilization benefiting from the charge-transfer (CT) state during electroluminescence (EL). Herein, we have explored a delicate balance to achieve both low ASE threshold and high EL exciton utilization by adopting a carefully tailored hybridized local and charge-transfer (HLCT) molecular design. A series of donor-π-acceptor (D-π-A) molecules (SBz-1, SBz-2 and SBz-3) are synthesized, and the structural change mainly refers to the spatial distance between D and A which could regulate the excited-state character via adjusting the CT length. Notably, the ASE phenomenon with a low threshold (2.97 μJ cm-2) and a high exciton utilization of 57.6% are achieved at the same time for SBz-2 with an appropriate CT length. The results provide guidance for molecular design toward harvesting triplet excitons in organic laser materials.

Photophysics of Anionic Bis(4H-imidazolato)CuI Complexes

In this paper, the photophysical behavior of four panchromatically absorbing, homoleptic bis(4H-imidazolato)CuI complexes, with a systematic variation in the electron-withdrawing properties of the imidazolate ligand, were studied by wavelength-dependent time-resolved femtosecond transient absorption spectroscopy. Excitation at 400, 480, and 630 nm populates metal-to-ligand charge transfer, intraligand charge transfer, and mixed-character singlet states. The pump wavelength-dependent transient absorption data were analyzed by a recently established 2D correlation approach. Data analysis revealed that all excitation conditions yield similar excited-state dynamics. Key to the excited-state relaxation is fast, sub-picosecond pseudo-Jahn-Teller distortion, which is accompanied by the relocalization of electron density onto a single ligand from the initially delocalized state at Franck-Condon geometry. Subsequent intersystem crossing to the triplet manifold is followed by a sub-100 ps decay to the ground state. The fast, nonradiative decay is rationalized by the low triplet-state energy as found by DFT calculations, which suggest perspective treatment at the strong coupling limit of the energy gap law.

Absorption behavior of doxorubicin hydrochloride in visible region in different environments: a combined experimental and computational study

The experimental absorption measurements in the interval 350-600 nm (Vis), molecular dynamics simulations, quantum-mechanics calculations and an advanced molecular treatment of simulation data are here combined to provide a complete picture of the absorption behavior in the visible portion of the electromagnetic spectrum of the doxorubicin hydrochloride (DX) molecule in different environments. By such an approach, we have shown that it is possible to characterize the effect of the environment on the DX absorption behavior - including the vibronic contributions - as well as to interpret such differences in terms of molecular electronic excited states, which are found to be strongly influenced by the environment.

Robust ΔSCF calculations with direct energy functional minimization methods and STEP for molecules and materials

The direct energy functional minimization method using the orbital transformation (OT) scheme in the program package CP2K has been employed for Δ self-consistent field (ΔSCF) calculations. The OT method for non-uniform molecular orbitals occupations allows us to apply the ΔSCF method for various kinds of molecules and periodic systems. Vertical excitation energies of heteroaromatic molecules and condensed phase systems, such as solvated ethylene and solvated uracil obeying periodic boundary conditions, are reported using the ΔSCF method. In addition, a Re-phosphate molecule attached to the surface of anatase (TiO₂) has been investigated. Additionally, we have implemented a recently proposed state-targeted energy projection ΔSCF algorithm [K. Carter-Fenk and J. M. Herbert, J. Chem. Theory Comput. 16(8), 5067-5082 (2020)] for diagonalization based SCF in CP2K. It is found that the OT scheme provides a smooth and robust SCF convergence for all investigated excitation energies and (non-)periodic systems.

The role of excited-state character, structural relaxation, and symmetry breaking in enabling delayed fluorescence activity in push-pull chromophores

Thermally activated delayed fluorescence (TADF) is a current promising route for generating highly efficient light-emitting devices. However, the design process of new chromophores is hampered by the complicated underlying photophysics. In this work, four closely related donor-π-acceptor-π-donor systems are investigated, two of which were synthesised previously, with the aim of elucidating their varying effectiveness for TADF. We outline that the frontier orbitals are insufficient for discriminating between the molecules. Subsequently, a detailed analysis of the excited states at a correlated ab initio level highlights the presence of a number of closely spaced singlet and triplet states of varying character. Results from five density functionals are compared against this reference revealing dramatic changes in, both, excited state energies and wavefunctions following variations in the amount of Hartree-Fock exchange included. Excited-state minima are optimised in solution showing the crucial role of structural variations and symmetry breaking for producing a strongly emissive S₁ state. The adiabatic singlet-triplet gaps thus obtained depend strongly on the range separation parameter used in the hybrid density functional calculations. More generally, this work highlights intricate differences present between singlet and triplet excited state wavefunctions and the challenges in describing them accurately.

Delocalized exciton formation in C60 linear molecular aggregates

Organic semiconducting materials containing C60 molecules are efficient acceptors for planar perovskite solar cells. In this work, we theoretically investigate the optical and excitonic properties of C60 linear molecular aggregates (composed of 1 to 7 C60 molecules) via the real-time-propagation rt-TDDFT technique. In the case of a single C60 molecule, the photoabsorption peaks are dominated by localized molecular excitons. We furthermore demonstrate that, in the case of linear molecular aggregates, the photoabsorption peaks are contributed by localized molecular excitons, charge transfer excitons, and Wannier-like delocalized excitons. This result is different to the accepted theory that only localized molecular excitons or charge transfer excitons can be produced in organic semiconducting materials. This work provides additional insights into the exciton formation in C60 molecular aggregates and may help in the rational design of efficient solar cells.

Q-Force: Quantum Mechanically Augmented Molecular Force Fields

The quality of molecular dynamics simulations strongly depends on the accuracy of the underlying force fields (FFs) that determine all intra- and intermolecular interactions of the system. Commonly, transferable FF parameters are determined based on a representative set of small molecules. However, such an approach sacrifices accuracy in favor of generality. In this work, an open-source and automated toolkit named Q-Force is presented, which augments these transferable FFs with molecule-specific bonded parameters and atomic charges that are derived from quantum mechanical (QM) calculations. The molecular fragmentation procedure allows treatment of large molecules (>200 atoms) with a low computational cost. The generated Q-Force FFs can be used at the same computational cost as transferable FFs, but with improved accuracy: We demonstrate this for the vibrational properties on a set of small molecules and for the potential energy surface on a complex molecule (186 atoms) with photovoltaic applications. Overall, the accuracy, user-friendliness, and minimal computational overhead of the Q-Force protocol make it widely applicable for atomistic molecular dynamics simulations.

Embelin's Versatile Photochemistry Makes It a Potent Photosensitizer for Photodynamic Therapy

Embelin, a natural product isolated from Embelia ribes, is a promising small-molecular drug for photodynamic anticancer therapy. Using modern quantum chemical methodology, embelin is shown to possess a versatile photochemistry comprising the capability of singlet oxygen generation, excited-state proton transfer, and oxidation. In particular, the detailed molecular mechanisms of singlet oxygen generation and proton transfer upon excitation are studied in great detail. While excited-state proton transfer can damage the protein itself, it also mediates intersystem crossing along its reaction pathway, thus facilitating singlet oxygen generation. When embelin is bound to proteins, all these processes can lead to protein damage and the desired phototoxicity.

Modulating Charge Separation and Intersystem Crossing in Donor-Switch-Acceptor Systems: A Computational Study

Charge separation and intersystem crossing play critical roles in various applications of organic long persistent luminescence materials, including light-emitting diodes, chemical sensors, theranostics, and many biomedical and information applications. Using first-principles calculations, we demonstrate that an azobenzene acting as a photoswitch can be used for altering the configuration of a donor-switch-acceptor (D-S-A) molecular system to ensure charge separation and promote intersystem crossing upon photoexcitation. The trans to cis photoisomerization of an azobenzene switch creates an electron trap that stabilizes the charge-separated state. The cis conformation further facilitates the singlet to triplet intersystem crossing in the excited state. Our theoretical study of the D-S-A system may help the design of long persistent luminescent organic devices.

Influence of the Environment on Shaping the Absorption of Monomeric Infrared Fluorescent Proteins

Infrared fluorescent proteins (iRFPs) are potential candidates for deep-tissue in vivo imaging. Here, we provide molecular-level insights into the role of the protein environment in the structural stability of the chromophore within the protein binding pocket through the flexible hydrogen-bonding network using molecular dynamics simulation. Furthermore, we present systematic excited-state analysis to characterize the nature of the first two excited states and the role of the environment in shaping the nature of the chromophore's excited states within the hybrid quantum mechanics/molecular mechanics framework. Our results reveal that the environment red-shifts the absorption of the chromophore by about 0.32 eV compared to the isolated counterpart, and besides the structural stability, the protein environment does not alter the nature of the excited state of the chromophore significantly. Our study contributes to the fundamental understanding of the excited-state processes of iRFPs in a complex environment and provides a design principle for developing iRFPs with desired spectral properties.

Evaluation of charge-transfer rates in fullerene-based donor-acceptor dyads with different density functional approximations

The shift towards renewable energy is one of the main challenges of this generation. Dye-sensitized solar cells (DSSCs), based on donor-acceptor architectures, can help in this transition as they present excellent photovoltaic efficiencies yet cheap and simple manufacturing. For molecular heterojunction DSSCs, donor-acceptor pairs are linked in a covalent manner, which facilitates their tailoring and rational design. Nevertheless, reliable computational characterization of charge transfer rate constants (kCT) is needed to speed this development process up. In this context, the performance of time-dependent density functional theory for the calculation of kCT values in donor-acceptor fullerene-based dyads has not been benchmarked yet. Herein, we present a detailed analysis on the performance of seven well-known density functional approximations (DFAs) for this type of system, focusing on several parameters such as the reorganization energies (λ), electronic couplings (VDA), and Gibbs energies (ΔG0CT), as well as the final rate constants. The amount of exact exchange at short range (SR) and long range (LR) electron-electron distances (and the transition from the SR to LR) turned out to be key for the success of the prediction. The tuning of these parameters improves significantly the performance of current DFAs.

Influence of the crystal packing in singlet fission: one step beyond the gas phase approximation

Singlet fission (SF), a multiexciton generation process, has been proposed as an alternative to enhance the performance of solar cells. The gas phase dimer model has shown its utility to study this process, but it does not always cover all the physics and the effect of the surrounding atoms has to be included in such cases. In this contribution, we explore the influence of crystal packing on the electronic couplings, and on the so-called exciton descriptors and electron–hole correlation plots. We have studied three tetracene dimers extracted from the crystal structure, as well as several dimers and trimers of the α and β polymorphs of 1,3-diphenylisobenzofuran (DPBF). These polymorphs show different SF yields. Our results highlight that the character of the excited states of tetracene depends on both the mutual disposition of molecules and inclusion of the environment. The latter does however not change significantly the interpretation of the SF mechanism in the studied systems. For DPBF, we establish how the excited state analysis is able to pinpoint differences between the polymorphs. We observe strongly bound correlated excitons in the β polymorph which might hinder the formation of the ¹TT state and, consequently, explain its low SF yield.

Singlet fission (SF), a multiexciton generation process, has been proposed as an alternative to enhance the performance of solar cells.

Information theoretical and thermodynamic view of the excited-state density functional theory of Coulomb systems

Excited states of Coulomb systems are studied within density functional theory with information theoretical quantities. The Ghosh-Berkowitz-Parr thermodynamic transcription is extended to excited states, and the concept of the local temperature is introduced. It is shown that extremization of information entropy or Fisher information results in a constant temperature. For Coulomb systems, there is a simple relation between the total energy and phase-space Fisher information. The phase-space fidelity between excited states is proportional to the position-space fidelity, with a factor of proportionality depending on total energies. The phase-space relative entropy is equal to the position-space relative entropy plus a term depending only on the total energies. The relationship between the phase-space fidelity susceptibility and Fisher information is also presented.

Spectral shift, electronic coupling and exciton delocalization in nanocrystal dimers: insights from all-atom electronic structure computations

Delocalization of excitons promoted by electronic coupling between clusters or quantum dots (QD) changes the dynamical processes in nanostructured aggregates enhancing energy transport. A spectroscopic shift of the absorption spectrum upon QD aggregation is commonly observed and ascribed to quantum mechanical coupling between neighbouring dots but also to exciton delocalization over the sulphur-based ligand shell or to other mechanisms as a change in the dielectric constant of the surrounding medium. We address the question of electronic coupling and exciton delocalization in nanocrystal aggregates by performing all-atom electronic structure calculations in models of colloidal QD dimers. The relation between spectral shift, interdot coupling and exciton delocalization is investigated in atomistic detail in models of dimers formed by CdSe clusters kept together by bridging organic ligands. Our results support the possibility of obtaining exciton delocalization over the dimer and point out the crucial role of the bridging ligand in enhancing interdot electronic coupling.

TheoDORE: A toolbox for a detailed and automated analysis of electronic excited state computations

The advent of ever more powerful excited-state electronic structure methods has led to a tremendous increase in the predictive power of computation, but it has also rendered the analysis of these computations much more challenging and time-consuming. TheoDORE tackles this problem through providing tools for post-processing excited-state computations, which automate repetitive tasks and provide rigorous and reproducible descriptors. Interfaces are available for ten different quantum chemistry codes and a range of excited-state methods implemented therein. This article provides an overview of three popular functionalities within TheoDORE, a fragment-based analysis for assigning state character, the computation of exciton sizes for measuring charge transfer, and the natural transition orbitals used not only for visualization but also for quantifying multiconfigurational character. Using the examples of an organic push-pull chromophore and a transition metal complex, it is shown how these tools can be used for a rigorous and automated assignment of excited-state character. In the case of a conjugated polymer, we venture beyond the limits of the traditional molecular orbital picture to uncover spatial correlation effects using electron-hole correlation plots and conditional densities.

Distortion-Controlled Redshift of Organic Dye Molecules

It is shown, quantum chemically, how structural distortion of an aromatic dye molecule can be leveraged to rationally tune its optoelectronic properties. By using a quantitative Kohn-Sham molecular orbital (KS-MO) approach, in combination with time-dependent DFT (TD-DFT), the influence of various structural and electronic tuning parameters on the HOMO-LUMO gap of a benzenoid model dye have been investigated. These parameters include 1) out-of-plane bending of the aromatic core, 2) bending of the bridge with respect to the core, 3) the nature of the bridge itself, and 4) π-π stacking. The study reveals the coupling of multiple structural distortions as a function of bridge length and number of bridges in benzene to be chiefly responsible for a decreased HOMO-LUMO gap, and consequently, red-shifting of the absorption wavelength associated with the lowest singlet excitation (λ≈560 nm) in the model cyclophane systems. These physical insights together with a rational approach for tuning the oscillator strength were leveraged for the proof-of-concept design of an intense near-infrared (NIR) absorbing cyclophane dye at λ=785 nm. This design may contribute to a new class of distortion-controlled NIR absorbing organic dye molecules.

Toward an understanding of electronic excitation energies beyond the molecular orbital picture

Can we gain an intuitive understanding of excitation energies beyond the molecular picture?

Synthesis and Cyclization-Induced Charge Transfer of Rectangular Bisterthiophenesiloxanes

Cyclization-modified terthiophene displays the change of emission behavior from locally excited (LE) to the intramolecular charge transfer (ICT) state. The rectangular bisterthiophenesiloxanes (DSiTh) was successfully prepared by π-π-stacking-aided hydrogen-bonding interactions. Cyclization-induced ICT in DSiTh could be observed, which was confirmed by absorption spectra, fluorescence spectra, and quantum chemistry analysis. The cyclization produces a strong intramolecular electron redistribution of a highly packed π-conjugated terthiophene. Thus, a distinctive variation of the dipole moment and a through-space ICT happen.

Electron-Hole Correlation as Unambiguous and Universal Classification for the Nature of Low-Lying ππ* States of Nitrogen Heterocycles

The ¹L_a and ¹L_b classification of electronically excited states of cata-condensed hydrocarbons proposed by Platt in 1949 ( Platt , J. R. J. Chem. Phys. 1949 , 17 , 484 ) is challenged by investigating a series of N-heteronaphthalenes and comparison of their low-lying ππ* excited states to those of naphthalene. The breakdown of Platt's classification scheme for N-heterocycles is highlighted, and a reliable and versatile alternative using exciton analyses is presented. The strength of electron-hole correlation turns out to be the most reliable distinguishing feature, and thus, an alternative nomenclature of ¹L_w (weakly correlated) and ¹L_s (strongly correlated) is proposed. Furthermore, fundamental guidelines for their property modulation through N-atom substitution patterns are discussed.

Structure of Electronically Reduced N-Donor Bidentate Ligands and Their Heteroleptic Four-Coordinate Zinc Complexes: A Survey of Density Functional Theory Results

The role of Hartree-Fock exchange in describing the structural changes occurring upon reduction of bipyridine-based ligands and their complexes is investigated within the framework of density functional theory (DFT) calculations. A set of four free ligands in their neutral and radical anionic forms, and two of their zinc complexes in their dicationic and monocationic radical forms, is used to compare a large panel of pure, conventional, and long-range corrected hybrid DFT functionals; coupled cluster single and double calculations are used alongside experimental results as benchmarks. Particular attention has been devoted to the magnitude of the change, upon reduction, of the Δ-parameter, which measures the difference between the C_py-C_py and the C-N bond lengths in bipyridine ligand and is known to experimentally correlate with the charge of the ligands. Our results indicate that the structural changes significantly depend on the amount of exact exchange included in the functional. A progressive evolution is observed for the free ligands, whereas two distinct sets of results are obtained for the complexes. Functionals with a small degree of HF exchange, e.g., B3LYP, do not adequately describe geometric changes for the considered species, and, quite surprisingly, the same holds for the CC2 method. The best agreement to experimental and CCSD values is obtained with functionals that include a significant but not excessive part of exact exchange, e.g., CAM-B3LYP, M06-2X, and ωB97X-D. The calculated localization of the added electron after reduction, which depends on the self-interaction error, is used to rationalize these outcomes. Static correlation is also shown to play a role in the accurate description of the electronic structure.

Density-based descriptors and exciton analyses for visualizing and understanding the electronic structure of excited states

Analysis and interpretation of the electronic structure of excited electronic states are prerequisites for developing a fundamental understanding of photochemistry and optical properties of molecular systems and an everyday task for a computational photochemist. Hence, wavefunction-based and density-based analysis tools have been devised over the last decades, and most recently also a family of quantitative exciton-wavefunction based descriptors has been developed. While the latter represent the main focus of this perspective, they are also discussed in the context of other existing analysis methods. Exciton analysis bridges the gap between the physically intuitive exciton picture and complex quantum-chemical wavefunctions by yielding insightful quantitative descriptors like exciton size, hole and electron size, electron-hole distance and exciton correlation. Thereby, not only a comprehensive characterization of the electronic structure is provided, but moreover, the formalism is automatizable and thus also optimally suited for benchmarking excited-state electronic structure methods.

Rational design of near-infrared absorbing organic dyes: Controlling the HOMO-LUMO gap using quantitative molecular orbital theory

Principles are presented for the design of functional near-infrared (NIR) organic dye molecules composed of simple donor (D), spacer (π), and acceptor (A) building blocks in a D-π-A fashion. Quantitative Kohn-Sham molecular orbital analysis enables accurate fine-tuning of the electronic properties of the π-conjugated aromatic cores by effecting their size, including silaaromatics, adding donor and acceptor substituents, and manipulating the D-π-A torsional angle. The trends in HOMO-LUMO gaps of the model dyes correlate with the excitation energies computed with time-dependent density functional theory at CAMY-B3LYP. Design principles could be developed from these analyses, which led to a proof-of-concept linear D-π-A with a strong excited-state intramolecular charge transfer and a NIR absorption at 879 nm. © 2018 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Dipole-Induced Transition Orbitals: A Novel Tool for Investigating Optical Transitions in Extended Systems

Optical absorption spectra for nanostructures and solids can be obtained from the macroscopic dielectric function within the random phase approximation. While experimental spectra can be reproduced with good accuracy, important properties, such as the charge-transfer character associated with a particular transition, are not retrievable. This contribution presents a computationally inexpensive method for the analysis of optical and excitonic properties for extended systems based on solely their electronic ground-state structure. We formulate a perturbative orbital transformation theory based on dipole-induced transition moments between orbitals, which yields correlated pairs of particle and hole functions. To demonstrate the potency of this new transformation formalism, we investigate the nature of excitations in inorganic molecular complexes and in extended systems. With our method, it is possible to extract mechanistic insights from the transitions observed in the optical spectrum, without requiring explicit calculation of the many-electron excited states.

Can TDDFT Describe Excited Electronic States of Naphthol Photoacids? A Closer Look with EOM-CCSD

The ¹L_b and ¹L_a excited states of naphthols are characterized by using time-dependent density functional theory (TDDFT), configuration interaction with singles (CIS), and equation-of-motion coupled cluster singles and doubles (EOM-CCSD) methods. TDDFT fails dramatically at predicting the energy and ordering of the ¹L_a and ¹L_b excited states as observed experimentally, while EOM-CCSD accurately predicts the excited states as characterized by natural transition orbital analysis. The limitations of TDDFT are attributed to the absence of correlation from doubly excited configurations as well as the inconsistent description of excited electronic states of naphthol photoacids revealed by excitation analysis based on the one-electron transition density matrix.

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.

Hot and Cold Charge-Transfer Mechanisms in Organic Photovoltaics: Insights into the Excited States of Donor/Acceptor Interfaces

The evolution of the excited-state manifold in organic D/A aggregates (e.g., the prototypical P3HT/PCBM) is investigated through a bottom-up approach via first-principles calculations. We show how the excited-state energies, the charge transfer (CT) states, and the electron-hole density distributions are strongly influenced by the size, the orientation, and the position (i.e., on-top versus on-edge phases) of P3HT/PCBM domains. We discuss how the structural order influences the excited-state electronic structure, providing an atomistic interpretation of the photophysics of organic blends. We show how the simultaneous presence of on-top and on-edge phases does not alter the optical absorption spectrum of the blend but does affect the photophysics. Photovoltaic processes such as (i) the simultaneous charge generation obtained from hot and cold excitations, (ii) the instantaneous and delayed charge separation, and (iii) the pump-push-probe charge generation can be interpreted based on our study.